Aleksandr Pankov

DNA Methylation and Somatic Mutations Converge on the Cell Cycle and Define Similar Evolutionary Histories in Brain Tumors

The evolutionary history of tumor cell populations can be reconstructed from patterns of genetic alterations. In contrast to stable genetic events, epigenetic states are reversible and sensitive to the microenvironment, prompting the question whether epigenetic information can similarly be used to discover tumor phylogeny. We examined the spatial and temporal dynamics of DNA methylation in a cohort of low-grade gliomas and their patient-matched recurrences. Genes transcriptionally upregulated through promoter hypomethylation during malignant progression to high-grade glioblastoma were enriched in cell cycle function, evolving in parallel with genetic alterations that deregulate the G1/S cell cycle checkpoint. Moreover, phyloepigenetic relationships robustly recapitulated phylogenetic patterns inferred from somatic mutations. These findings highlight widespread co-dependency of genetic and epigenetic events throughout brain tumor evolution.

Intratumoral Heterogeneity of the Epigenome

Investigation into intratumoral heterogeneity (ITH) of the epigenome is in a formative stage. The patterns of tumor evolution inferred from epigenetic ITH and genetic ITH are remarkably similar, suggesting widespread co-dependency of these disparate mechanisms. The biological and clinical relevance of epigenetic ITH are becoming more apparent. Rare tumor cells with unique and reversible epigenetic states may drive drug resistance, and the degree of epigenetic ITH at diagnosis may predict patient outcome. This perspective presents these current concepts and clinical implications of epigenetic ITH, and the experimental and computational techniques at the forefront of ITH exploration.

Clonal expansion and epigenetic reprogramming following deletion or amplification of mutant IDH1

Significance Identifying the drivers of tumorigenesis provides insight into mechanisms of transformation and can suggest novel therapeutic targets. IDH1 mutations in gliomas are one such promising target. Drivers of tumor initiation may be distinct from those at tumor recurrence, however. Here, we demonstrate that in a subset of initially IDH1 mutant gliomas IDH1 is deleted or amplified at recurrence, yielding a higher grade tumor with a reprogrammed epigenome. We also report systematic selection for cells with IDH1 CNA in vitro and in vivo. Thus, while IDH1 mutation likely initiates gliomagenesis, neither mutant IDH1 nor the oncometabolite 2HG that it produces are required at recurrence. These findings have important implications for emerging therapeutic strategies targeting mutant IDH1. IDH1 mutation is the earliest genetic alteration in low-grade gliomas (LGGs), but its role in tumor recurrence is unclear. Mutant IDH1 drives overproduction of the oncometabolite d-2-hydroxyglutarate (2HG) and a CpG island (CGI) hypermethylation phenotype (G-CIMP). To investigate the role of mutant IDH1 at recurrence, we performed a longitudinal analysis of 50 IDH1 mutant LGGs. We discovered six cases with copy number alterations (CNAs) at the IDH1 locus at recurrence. Deletion or amplification of IDH1 was followed by clonal expansion and recurrence at a higher grade. Successful cultures derived from IDH1 mutant, but not IDH1 wild type, gliomas systematically deleted IDH1 in vitro and in vivo, further suggestive of selection against the heterozygous mutant state as tumors progress. Tumors and cultures with IDH1 CNA had decreased 2HG, maintenance of G-CIMP, and DNA methylation reprogramming outside CGI. Thus, while IDH1 mutation initiates gliomagenesis, in some patients mutant IDH1 and 2HG are not required for later clonal expansions.

Data-driven regions of interest for longitudinal change in frontotemporal lobar degeneration

Current research is investigating the potential utility of longitudinal measurement of brain structure as a marker of drug effect in clinical trials for neurodegenerative disease. Recent studies in Alzheimers disease (AD) have shown that measurement of change in empirically derived regions of interest (ROIs) allows more reliable measurement of change over time compared with regions chosen a-priori based on known effects of AD on brain anatomy. Frontotemporal lobar degeneration (FTLD) is a devastating neurodegenerative disorder for which there are no approved treatments. The goal of this study was to identify an empirical ROI that maximizes the effect size for the annual rate of brain atrophy in FTLD compared with healthy age matched controls, and to estimate the effect size and associated power estimates for a theoretical study that would use change within this ROI as an outcome measure. Eighty six patients with FTLD were studied, including 43 who were imaged twice at 1.5 T and 43 at 3 T, along with 105 controls (37 imaged at 1.5 T and 67 at 3 T). Empirically-derived maps of change were generated separately for each field strength and included the bilateral insula, dorsolateral, medial and orbital frontal, basal ganglia and lateral and inferior temporal regions. The extent of regions included in the 3 T map was larger than that in the 1.5 T map. At both field strengths, the effect sizes for imaging were larger than for any clinical measures. At 3 T, the effect size for longitudinal change measured within the empirically derived ROI was larger than the effect sizes derived from frontal lobe, temporal lobe or whole brain ROIs. The effect size derived from the data-driven 1.5 T map was smaller than at 3 T, and was not larger than the effect size derived from a-priori ROIs. It was estimated that measurement of longitudinal change using 1.5 T MR systems requires approximately a 3-fold increase in sample size to obtain effect sizes equivalent to those seen at 3 T. While the results should be confirmed in additional datasets, these results indicate that empirically derived ROIs can reduce the number of subjects needed for a longitudinal study of drug effects in FTLD compared with a-priori ROIs. Field strength may have a significant impact on the utility of imaging for measuring longitudinal change.

NSeq: a multithreaded Java application for finding positioned nucleosomes from sequencing data

We introduce NSeq, a fast and efficient Java application for finding positioned nucleosomes from the high-throughput sequencing of MNase-digested mononucleosomal DNA. NSeq includes a user-friendly graphical interface, computes false discovery rates (FDRs) for candidate nucleosomes from Monte Carlo simulations, plots nucleosome coverage and centers, and exploits the availability of multiple processor cores by parallelizing its computations. Java binaries and source code are freely available at https://github.com/songlab/NSeq. The software is supported on all major platforms equipped with Java Runtime Environment 6 or later.

Mutant IDH1 expression is associated with down-regulation of monocarboxylate transporters

Mutations in isocitrate dehydrogenase 1 (IDH1) are characteristic of low-grade gliomas. We recently showed that mutant IDH1 cells reprogram cellular metabolism by down-regulating pyruvate dehydrogenase (PDH) activity. Reduced pyruvate metabolism via PDH could lead to increased pyruvate conversion to lactate. The goal of this study was therefore to investigate the impact of the IDH1 mutation on the pyruvate-to-lactate flux. We used 13C magnetic resonance spectroscopy and compared the conversion of hyperpolarized [1-13C]-pyruvate to [1-13C]-lactate in immortalized normal human astrocytes expressing mutant or wild-type IDH1 (NHAIDHmut and NHAIDHwt). Our results indicate that hyperpolarized lactate production is reduced in NHAIDHmut cells compared to NHAIDHwt. This reduction was associated with lower expression of the monocarboxylate transporters MCT1 and MCT4 in NHAIDHmut cells. Furthermore, hyperpolarized lactate production was comparable in lysates of NHAIDHmut and NHAIDHwt cells, wherein MCTs do not impact hyperpolarized pyruvate delivery and lactate production. Collectively, our findings indicated that lower MCT expression was a key contributor to lower hyperpolarized lactate production in NHAIDHmut cells. The SLC16A3 (MCT4) promoter but not SLC16A1 (MCT1) promoter was hypermethylated in NHAIDHmut cells, pointing to possibly different mechanisms mediating reduced MCT expression. Finally analysis of low-grade glioma patient biopsy data from The Cancer Genome Atlas revealed that MCT1 and MCT4 expression was significantly reduced in mutant IDH1 tumors compared to wild-type. Taken together, our study shows that reduced MCT expression is part of the metabolic reprogramming of mutant IDH1 gliomas. This finding could impact treatment and has important implications for metabolic imaging of mutant IDH1 gliomas.

An Atlas of the Human Kinome Reveals the Mutational Landscape Underlying Dysregulated Phosphorylation Cascades in Cancer

Kinase inhibitors are used widely to treat various cancers, but adaptive reprogramming of kinase cascades and activation of feedback loop mechanisms often contribute to therapeutic resistance. Determining comprehensive, accurate maps of kinase circuits may therefore help elucidate mechanisms of response and resistance to kinase inhibitor therapies. In this study, we identified and validated phosphorylatable target sites across human cell and tissue types to generate PhosphoAtlas, a map of 1,733 functionally interconnected proteins comprising the human phospho-reactome. A systematic curation approach was used to distill protein phosphorylation data cross-referenced from 38 public resources. We demonstrated how a catalog of 2,617 stringently verified heptameric peptide regions at the catalytic interface of kinases and substrates could expose mutations that recurrently perturb specific phospho-hubs. In silico mapping of 2,896 nonsynonymous tumor variants identified from thousands of tumor tissues also revealed that normal and aberrant catalytic interactions co-occur frequently, showing how tumors systematically hijack, as well as spare, particular subnetworks. Overall, our work provides an important new resource for interrogating the human tumor kinome to strategically identify therapeutically actionable kinase networks that drive tumorigenesis. Cancer Res; 76(7); 1733-45. ©2016 AACR.

NF45 and NF90/NF110 coordinately regulate ESC pluripotency and differentiation

While years of investigation have elucidated many aspects of embryonic stem cell (ESC) regulation, the contributions of post-transcriptional and translational mechanisms to the pluripotency network remain largely unexplored. In particular, little is known in ESCs about the function of RNA binding proteins (RBPs), the protein agents of post-transcriptional regulation. We performed an unbiased RNAi screen of RBPs in an ESC differentiation assay and identified two related genes, NF45 (Ilf2) and NF90/NF110 (Ilf3), whose knockdown promoted differentiation to an epiblast-like state. Characterization of NF45 KO, NF90 + NF110 KO, and NF110 KO ESCs showed that loss of NF45 or NF90 + NF110 impaired ESC proliferation and led to dysregulated differentiation down embryonic lineages. Additionally, we found that NF45 and NF90/NF110 physically interact and influence the expression of each other at different levels of regulation. Globally across the transcriptome, NF45 KO ESCs and NF90 + NF110 KO ESCs show similar expression changes. Moreover, NF90 + NF110 RNA immunoprecipitation (RIP)-seq in ESCs suggested that NF90/NF110 directly regulate proliferation, differentiation, and RNA-processing genes. Our data support a model in which NF45, NF90, and NF110 operate in feedback loops that enable them, through both overlapping and independent targets, to help balance the push and pull of pluripotency and differentiation cues.

Data-driven regions of interest for longitudinal change in three variants of frontotemporal lobar degeneration

Longitudinal imaging of neurodegenerative disorders is a potentially powerful biomarker for use in clinical trials. In Alzheimers disease, studies have demonstrated that empirically derived regions of interest (ROIs) can provide more reliable measurement of disease progression compared with anatomically defined ROIs.

Abstract IA08: An epigenome perspective of human tumor evolution

Cancers develop through a process of clonal evolution in which ongoing genetic and epigenetic diversification allows for repeated cycles of sub-clonal selection and expansion (Greaves and Maley, 2012; Nowell, 1976). As a result, human tumors can display substantial intratumoral heterogeneity, including discordant genetic alterations between the initial tumor and local recurrence or distant metastases (Gerlinger et al., 2012; Okosun et al., 2014; Wu et al., 2012; Yachida et al., 2010). As mutations accumulate over time, genomic profiling of spatially or temporally separated tumor samples can be used to reconstruct the evolutionary history and underlying clonal architectures of individual tumors (Gerlinger et al., 2014). In contrast to stable genetic events, epigenetic states are reversible and may change in response to the tumor microenvironment, prompting the question whether epigenetic information can similarly be used to discover tumor phylogeny. In low-grade glioma, the course of tumor evolution is particularly clinically significant. Low-grade gliomas are diffuse, infiltrative tumors that recur locally and may unpredictably undergo malignant progression to a higher grade with a worse prognosis (Sanai et al., 2011). Recurrences that progress to highly malignant WHO grade IV glioblastoma (GBM) acquire genetic alterations in the RB and AKT-mTOR pathways (Johnson et al., 2014; Louis, 2006). In fact, adjuvant treatment with alkylating chemotherapeutics such as temozolomide (TMZ) can induce hypermutation that emerges in recurrent tumors (Bodell et al, 2003; Hunter et al., 2006), and we recently linked treatment-associated driver mutations in these two pathways to malignant progression of grade II glioma to GBM (Johnson et al., 2014). The treatment associated malignant progression follows selection of tumor cells with epigenetic silencing of the DNA repair protein MGMT (van Thiujl et al, 2015). The critical role that epigenetic alterations play in the development and therapeutic response of gliomas is increasingly being appreciated (Fouse and Costello, 2009). Somatic mutations in IDH1 or IDH2 may be the first genetic driver in the development of many low-grade gliomas (Johnson et al., 2014; Lai et al., 2011; Watanabe et al., 2009). Genetic mutations in IDH genes induce a pattern of early epigenetic alterations known as the glioma CpG island methylator phenotype (G-CIMP) characterized by extensive remodeling of the DNA methylome (Hill et al., 2014; Noushmehr et al., 2010; Toyota et al., 1999; Turcan et al., 2012). The inactivation of other genes mutated in low-grade gliomas, such as ATRX (Jiao et al., 2012) and SMARCA4 (Johnson et al., 2014), is known to induce specific DNA methylation changes as well (Banine et al., 2005; Gibbons et al., 2000). Although there has been extensive characterization of tumor methylomes using a single sampling per tumor, little is known about intratumoral heterogeneity at the epigenetic level or of temporal evolution of the low-grade glioma methylome and its relationship to the genome. In this presentation, I will show that phylogenic analysis of spatial and temporal patterns of either reversible DNA methylation or irreversible somatic mutations independently yield remarkably similar evolutionary histories. I will also incorporate this inter-dependent evolution into a detailed model of brain tumorigenesis that extends from the first mutation and epimutations through tumor recurrence. Citation Format: Tali Mazor, Aleksandr Pankov, Jun Song, Joseph Costello. An epigenome perspective of human tumor evolution. [abstract]. In: Proceedings of the AACR Special Conference on Chromatin and Epigenetics in Cancer; Sep 24-27, 2015; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2016;76(2 Suppl):Abstract nr IA08.