Thais S. Sabedot

Molecular Profiling Reveals Biologically Discrete Subsets and Pathways of Progression in Diffuse Glioma

Therapy development for adult diffuse glioma is hindered by incomplete knowledge of somatic glioma driving alterations and suboptimal disease classification. We defined the complete set of genes associated with 1,122 diffuse grade II-III-IV gliomas from The Cancer Genome Atlas and used molecular profiles to improve disease classification, identify molecular correlations, and provide insights into the progression from low- to high-grade disease. Whole-genome sequencing data analysis determined that ATRX but not TERT promoter mutations are associated with increased telomere length. Recent advances in glioma classification based on IDH mutation and 1p/19q co-deletion status were recapitulated through analysis of DNA methylation profiles, which identified clinically relevant molecular subsets. A subtype of IDH mutant glioma was associated with DNA demethylation and poor outcome; a group of IDH-wild-type diffuse glioma showed molecular similarity to pilocytic astrocytoma and relatively favorable survival. Understanding of cohesive disease groups may aid improved clinical outcomes.

TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data

The Cancer Genome Atlas (TCGA) research network has made public a large collection of clinical and molecular phenotypes of more than 10 000 tumor patients across 33 different tumor types. Using this cohort, TCGA has published over 20 marker papers detailing the genomic and epigenomic alterations associated with these tumor types. Although many important discoveries have been made by TCGAs research network, opportunities still exist to implement novel methods, thereby elucidating new biological pathways and diagnostic markers. However, mining the TCGA data presents several bioinformatics challenges, such as data retrieval and integration with clinical data and other molecular data types (e.g. RNA and DNA methylation). We developed an R/Bioconductor package called TCGAbiolinks to address these challenges and offer bioinformatics solutions by using a guided workflow to allow users to query, download and perform integrative analyses of TCGA data. We combined methods from computer science and statistics into the pipeline and incorporated methodologies developed in previous TCGA marker studies and in our own group. Using four different TCGA tumor types (Kidney, Brain, Breast and Colon) as examples, we provide case studies to illustrate examples of reproducibility, integrative analysis and utilization of different Bioconductor packages to advance and accelerate novel discoveries.

RGBM: regularized gradient boosting machines for identification of the transcriptional regulators of discrete glioma subtypes

Abstract We propose a generic framework for gene regulatory network (GRN) inference approached as a feature selection problem. GRNs obtained using Machine Learning techniques are often dense, whereas real GRNs are rather sparse. We use a Tikonov regularization inspired optimal L-curve criterion that utilizes the edge weight distribution for a given target gene to determine the optimal set of TFs associated with it. Our proposed framework allows to incorporate a mechanistic active biding network based on cis-regulatory motif analysis. We evaluate our regularization framework in conjunction with two non-linear ML techniques, namely gradient boosting machines (GBM) and random-forests (GENIE), resulting in a regularized feature selection based method specifically called RGBM and RGENIE respectively. RGBM has been used to identify the main transcription factors that are causally involved as master regulators of the gene expression signature activated in the FGFR3-TACC3-positive glioblastoma. Here, we illustrate that RGBM identifies the main regulators of the molecular subtypes of brain tumors. Our analysis reveals the identity and corresponding biological activities of the master regulators characterizing the difference between G-CIMP-high and G-CIMP-low subtypes and between PA-like and LGm6-GBM, thus providing a clue to the yet undetermined nature of the transcriptional events among these subtypes.

SIRT1 regulates Mxd1 during malignant melanoma progression

In a murine melanoma model, malignant transformation promoted by a sustained stress condition was causally related to increased levels of reactive oxygen species resulting in DNA damage and massive epigenetic alterations. Since the chromatin modifier Sirtuin-1 (SIRT1) is a protein attracted to double-stranded DNA break (DSB) sites and can recruit other components of the epigenetic machinery, we aimed to define the role of SIRT1 in melanomagenesis through our melanoma model. The DNA damage marker, γH2AX was found increased in melanocytes after 24 hours of deadhesion, accompanied by increased SIRT1 expression and decreased levels of its target, H4K16ac. Moreover, SIRT1 started to be associated to DNMT3B during the stress condition, and this complex was maintained along malignant progression. Mxd1 was identified by ChIP-seq among the DNA sequences differentially associated with SIRT1 during deadhesion and was shown to be a common target of both, SIRT1 and DNMT3B. In addition, Mxd1 was found downregulated from pre-malignant melanocytes to metastatic melanoma cells. Treatment with DNMT inhibitor 5AzaCdR reversed the Mxd1 expression. Sirt1 stable silencing increased Mxd1 mRNA expression and led to down-regulation of MYC targets, such as Cdkn1a, Bcl2 and Psen2, whose upregulation is associated with human melanoma aggressiveness and poor prognosis. We demonstrated a novel role of the stress responsive protein SIRT1 in malignant transformation of melanocytes associated with deadhesion. Mxd1 was identified as a new SIRT1 target gene. SIRT1 promoted Mxd1 silencing, which led to increased activity of MYC oncogene contributing to melanoma progression.

A non-functional galanin receptor-2 in a multiple sclerosis patient

Multiple Sclerosis (MS) is an inflammatory neurodegenerative disease that affects approximately 2.5 million people globally. Even though the etiology of MS remains unknown, it is accepted that it involves a combination of genetic alterations and environmental factors. Here, after performing whole exome sequencing, we found a MS patient harboring a rare and homozygous single nucleotide variant (SNV; rs61745847) of the G-protein coupled receptor (GPCR) galanin-receptor 2 (GALR2) that alters an important amino acid in the TM6 molecular toggle switch region (W249L). Nuclear magnetic resonance imaging showed that the hypothalamus (an area rich in GALR2) of this patient exhibited an important volumetric reduction leading to an enlarged third ventricle. Ex vivo experiments with patient-derived blood cells (AKT phosphorylation), as well as studies in recombinant cell lines expressing the human GALR2 (calcium mobilization and NFAT mediated gene transcription), showed that galanin (GAL) was unable to stimulate cell signaling in cells expressing the variant GALR2 allele. Live cell confocal microscopy showed that the GALR2 mutant receptor was primarily localized to intracellular endosomes. We conclude that the W249L SNV is likely to abrogate GAL-mediated signaling through GALR2 due to the spontaneous internalization of this receptor in this patient. Although this homozygous SNV was rare in our MS cohort (1:262 cases), our findings raise the potential importance of impaired neuroregenerative pathways in the pathogenesis of MS, warrant future studies into the relevance of the GAL/GALR2 axis in MS and further suggest the activation of GALR2 as a potential therapeutic route for this disease.

Abstract 780: Multi-omic profiling of gliomas reveals distinct DNA methylation changes at tumor recurrence

Varying possibilities of tumor relapse for lower grade glioma (LGG, WHO grade II and III) and glioblastoma (GBM, WHO grade IV) have led to heterogenous clinical outcomes for patients. Our current study aims to establish a molecular profile of recurrence of matched primary and recurrent LGG (n = 28) and recurrent GBM (n = 24) tumor samples. The Cancer Genome Atlas (TCGA) has comprehensively profiled these matched primary/recurrent tumor sets; whole genomes, coding exomes, methylomes, and transcriptomes have been sequenced, and the samples have undergone targeted profiling of the proteome. While unsupervised analysis techniques have not led to a clear recurrence signature, supervised analysis methods have revealed interesting patterns. Protein profiling has shown that recurrent gliomas retain the overall molecular signature of their primary counterpart, but the DNA damage response, apoptosis and RTK pathways are downregulated in the recurrent gliomas, in contrast to RAS/MAPK, PI3K/AKT, and EMT pathways, which are upregulated. Whole genome sequencing and rearrangement analysis have revealed increased genomic complexity among most recurrent gliomas as well as new fusions of interest in recurrent LGG samples (PTPRZ1-MET and involving ATRX). Using genome-wide Illumina HumanMethylation 450K data we observed that 78.6% of LGGs showed depletion of DNA methylation at recurrence and 50% of GBM tumors showed an enrichment of DNA methylation at recurrence. Patient centric enrichment analysis allowed us to discover a candidate biological subgroup characterized by a subset of LGG recurrences (50%) exhibiting an aberrant CpG methylation loss at inintergenic opensea regions when compared with canonical CpG islands and shores (fold > 1.3 and confidence intervals of 99%). Importantly, inspection of CpG sites significantly hypomethylated at openseas showed that this pronounced epigenetic signature maps to candidate TSS distal and hypomethylated enhancers. The gene-targets of these hypomethylated CpG sites show a corresponding up-regulation of expression. Pathway analysis has demonstrated that these upregulated genes are involved in cellular growth and proliferation, cellular function and maintenance, and cell cycle regulation. Our results provide evidence that DNA methylation may represent a stable signature of glioma recurrence and that the crosstalk between DNA hypomethylation at openseas and chromosomal instability may be involved in glioma recurrence. We plan to further integrate our findings between data types and correlate with treatment and patient clinical outcome. Citation Format: Lindsay C. Stetson, Camila Ferreira de Souza, Tathiane Maistro Malta, Thais Sarraf Sabedot, Quinn Ostrom, Peter Liao, Daniela Pretti da Cunha Tirapelli, Luciano Neder, Carlos Gilberto Carlotti, Rehan Akbani, Sofie Salama, Laila Poisson, Daniel Brat, The Cancer Genome Atlas Network, Houtan Noushmehr, Jill Barnholtz-Sloan. Multi-omic profiling of gliomas reveals distinct DNA methylation changes at tumor recurrence. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 780.

Abstract 4767: DNA sequences differentially associated with Sirt1 and Dnmt3b during melanoma progression

Proceedings: AACR 106th Annual Meeting 2015; April 18-22, 2015; Philadelphia, PA Melanoma is one of the most worrying tumors due to its high capacity to cause metastasis. Sirt1, a class III histone deacetylase, was shown to be attracted to sites of double strand breaks in DNA, favoring, by recruiting components of the epigenetic machinery, aberrant epigenetic silencing in that region (O’HAGAN et al., 2008). Beyond that, subsequent studies have attributed important roles of Sirt1 and other components of the epigenetic machinery, such as DNA methyltransferases (such as Dnmt3b), in tumor initiation and progression. In this sense, using a murine model of malignant transformation of melanocytes, associated with sustained stress conditions, our goal was to identify genomic sequences differentially associated to Sirt1 and/or Dnmt3B before and after stress condition, represented in our model by anchorage blockade, and analyze if these sequences are targeted by aberrant methylation in cell lines obtained by this process that represent different steps in melanoma progression. We have showed increased DNA damage after 24 hours of anchorage impediment followed by increased Sirt1 expression. In addition, we have showed Sirt1/Dnmt3b association just after this stress condition that leads to malignant transformation. It was identified several DNA sequences differentially associated to Sirt1 and also to Dnmt3b after stress condition which are also associated to genes related to different metabolic processes. Some of these sequences remain associated to Sirt1 and Dnmt3b during melanoma progression, indicating that Dnmt3b remains associated to Sirt1 during malignant transformation. Among the sequences, Tiam1 and Mxd1 are potential targets of Sirt1 and Dnmt3b that have important roles in tumor progression. We have also identified motifs where Sirt1 and Dnmt3b are probably bound before and after stress condition. The main motif related to Sirt1 was coincident with Boris motif. Subsequent analysis are been developed to elucidate the importance of these interesting discoveries. Supported by FAPESP and DAAD. Citation Format: Fabiana Marcelino Meliso, Camila Tainah da Silva, Simon Coetzee, Thais Sarraf Sabedot, Houtan Noushmehr, Regine Schneider-Stock, Miriam Galvonas Jasiulionis. DNA sequences differentially associated with Sirt1 and Dnmt3b during melanoma progression. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4767. doi:10.1158/1538-7445.AM2015-4767