Search

Activating mutations in the Telomerase Reverse Transcriptase (TERT) promoter are the single most common non-coding mutation in cancer and enable limitless cell division characteristic of immortal cells1–12. Solving the immortality mechanism represents a major step towards selectively reversing it in cancer cells. TERT promoter mutations create a de novo E26 transformation specific (ETS) transcription factor binding motif, however most of the 28 ETS factors and many other transcriptional regulators have been implicated13–22. Cancer type and mutation specific mechanisms have also been proposed. Here, we uniformly and robustly analyzed fifty-three cell lines representing sixteen cancer types and six recurrent mutations and found that a tetramer of the GA-binding protein (GABP) is specifically responsible for mutant TERT promoter activation in all cases, with no such role in TERT promoter wild type cells. Strikingly, TERT expression is maintained in tetramer depleted tumor cells. We show how and why the tetramer is serially replaced, not by other transcription factors, but by GABP dimers and then weakly by a paralogous tetramer complex. Elimination of the tetramer and dimer reinstates epigenetic repression of TERT, activates checkpoint programs and prevents cancer cell division. We conclude that unique features of the GABP tetramer must therefore determine positive selection of nearly all TERT promoter mutations in human cancer. Furthermore, domains shared among the three GABP complexes present pan-cancer vulnerabilities.

Purpose: To address the association between sequence variants within the MGMT (O6-methylguanine-DNA methyltransferase) promoter–enhancer region and methylation of MGMT in premalignant lesions from smokers and lung adenocarcinomas, their biological effects on gene regulation, and targeting MGMT for therapy.Experimental Design: Single nucleotide polymorphisms (SNP) identified through sequencing a 1.9 kb fragment 5′ of MGMT were examined in relation to MGMT methylation in 169 lung adenocarcinomas and 1,731 sputum samples from smokers. The effect of promoter haplotypes on MGMT expression was tested using a luciferase reporter assay and cDNA expression analysis along with allele-specific sequencing for methylation. The response of MGMT methylated lung cancer cell lines to the alkylating agent temozolomide (TMZ) was assessed.Results: The A allele of rs16906252 and the haplotype containing this SNP were strongly associated with increased risk for MGMT methylation in adenocarcinomas (ORs ≥ 94). This association was observed to a lesser extent in sputum samples in both smoker cohorts. The A allele was selectively methylated in primary lung tumors and cell lines heterozygous for rs16906252. With the most common haplotype as the reference, a 20 to 41% reduction in promoter activity was seen for the haplotype carrying the A allele that correlated with lower MGMT expression. The sensitivity of lung cancer cell lines to TMZ was strongly correlated with levels of MGMT methylation and expression.Conclusions: These studies provide strong evidence that the A allele of a MGMT promoter–enhancer SNP is a key determinant for MGMT methylation in lung carcinogenesis. Moreover, TMZ treatment may benefit a subset of lung cancer patients methylated for MGMT. Clin Cancer Res; 17(7); 2014–23. ©2011 AACR.

Tumors arise in part from the deleterious effects of genetic and epigenetic mechanisms on gene expression. In several mouse models of human tumors, the tumorigenic phenotype is reversible, suggesting that epigenetic mechanisms also contribute significantly to tumorigenesis in mice. It is not known whether these are the same epigenetic mechanisms in human and mouse tumors or whether they affect homologous genes. Using an integrated approach for genome-wide methylation and copy number analyses, we identified SLC5A8 on chromosome 12q23.1 that was affected frequently by aberrant methylation in human astrocytomas and oligodendrogliomas. SLC5A8 encodes a sodium monocarboxylate cotransporter that was highly expressed in normal brain but was significant down-regulated in primary gliomas. Bisulfite sequencing analysis showed that the CpG island was unmethylated in normal brain but frequently localized methylated in brain tumors, consistent with the tumor-specific loss of gene expression. In glioma cell lines, SLC5A8 expression was also suppressed but could be reactivated with a methylation inhibitor. Expression of exogenous SLC5A8 in LN229 and LN443 glioma cells inhibited colony formation, suggesting that it may function as a growth suppressor in normal brain cells. Remarkably, 9 of 10 murine oligodendroglial tumors (from p53+/− or ink4a/arf+/− animals transgenic for S100β-v-erbB) showed a similar tumor-specific down-regulation of mSLC5A8, the highly conserved mouse homologue. Taken together, these data suggest that SLC5A8 functions as a growth suppressor gene in vitro and that it is silenced frequently by epigenetic mechanisms in primary gliomas. The shared epigenetic inactivation of mSLC5A8 in mouse gliomas indicates an additional degree of commonality in the origin and/or pathway to tumorigenesis between primary human tumors and these mouse models of gliomas.

Supplementary Material from Shared Epigenetic Mechanisms in Human and Mouse Gliomas Inactivate Expression of the Growth Suppressor SLC5A8

Supplementary Figures S1 and S2 from Shared Epigenetic Mechanisms in Human and Mouse Gliomas Inactivate Expression of the Growth Suppressor SLC5A8

Background TERT promoter mutations are observed in 80% of wild-type IDH glioblastoma (GBM). Moreover, the upstream TERT transcription factor GABPB1 was recently identified as a cancer-specific therapeutic target for tumors harboring a TERT promoter mutation. In that context, noninvasive imaging biomarkers are needed for the detection of TERT modulation. Methods Multiple GBM models were investigated as cells and in vivo tumors and the impact of TERT silencing, either directly or by targeting GABPB1, was determined using 1H and hyperpolarized 13C magnetic resonance spectroscopy (MRS). Changes in associated metabolic enzymes were also investigated. Results 1H-MRS revealed that lactate and glutathione (GSH) were the most significantly altered metabolites when either TERT or GABPB1 was silenced, and lactate and GSH levels were correlated with cellular TERT expression. Consistent with the drop in lactate, 13C-MRS showed that hyperpolarized [1-13C]lactate production from [1-13C]pyruvate was also reduced when TERT was silenced. Mechanistically, the reduction in GSH was associated with a reduction in pentose phosphate pathway flux, reduced activity of glucose-6-phosphate dehydrogenase, and reduced NADPH. The drop in lactate and hyperpolarized lactate were associated with reductions in glycolytic flux, NADH, and expression/activity of GLUT1, monocarboxylate transporters, and lactate dehydrogenase A. Conclusions Our study indicates that MRS-detectable GSH, lactate, and lactate production could serve as metabolic biomarkers of response to emerging TERT-targeted therapies for GBM with activating TERT promoter mutations. Importantly these biomarkers are readily translatable to the clinic, and thus could ultimately improve GBM patient management.

BACKGROUND: The efficacy of immunotherapy in gliomas is limited by tumor heterogeneity and low mutational burden. Recent studies revealed increased alternative splicing (AS) in various cancer types that potentially translate into targetable neoantigens. We developed a novel comprehensive in silico pipeline for detecting tumor-specific AS events (neojunctions), and successfully identified tumor-wide, public, alternatively spliced neoantigens (ASNs) that elicited CD8+ T-cell-mediated immune responses. METHODS: Our computational pipeline first identified recurring splice junctions in high purity TCGA LGG/GBM bulk RNA-seq samples (n=429, PSR > 10%) and not in GTEx normal tissue bulk RNA-seq data (n=9166, PSR < 1%). We subsequently employed two independent algorithms to predict the peptide processing likelihood and the HLA-binding affinity of ASN candidates before validating neojunction and ASN expression in RNA-seq data from patient-derived glioma cell lines (n=68) and spatially-mapped glioma samples (n=535), and in the Clinical Proteomic Tumor Analysis Consortium (CPTAC) GBM mass spectrometry data (n=99). In vitro sensitization of healthy donor-derived CD8+ T-cells against our top ASN candidates followed by 10x VDJ scRNA-seq was performed to identify ASN-specific TCR sequences. We transduced these TCRs into Jurkat76s and co-cultured them with ASN-pulsed T2 cells and ASN and HLA-expressing COS7 cells to determine whether the TCR recognized the presented HLA:ASN. RESULTS: Our analysis identified 249 putative neojunctions that translated into 222 cancer-specific peptides sequences. Iterating across these peptides yielded 17,562 glioma-specific n-mer sequences between 8 to 11 amino acids in length. Both prediction algorithms concurrently identified 636 n-mers likely to be presented by demographically common HLA haplotypes. Through transcriptomic and proteomic validation with patient-derived samples and cell lines, we identified 8 ASNs that were tumor-wide, 4 of which were predicted to be HLA*A0201-presented. In vitro sensitization of healthy-donor derived CD8+ T-cells against these 4 ASNs resulted in an expansion of reactive CD8+ T-cells against ASNs derived from neojunctions within S100A6, RPL22, and GNAS. Subsequent 10x VDJ scRNA-seq on the expanded CD8+ population identified mutRPL22-reactive TCR clonotypes with neoantigen-mediated increases in IFNG and GZMB signatures. When transduced and expressed in Jurkat76 cells, one particular TCR clonotype demonstrated recognition and immunogenic activation against mutRPL22-pulsed T2 cells and COS7 cells endogenously expressing both HLA*A0201 and mutRPL22. CONCLUSION: Our unique integrative pipeline detected novel tumor-wide splice-derived neoantigen candidates, and ASN-reactive TCRs identified through our pipeline offer a new avenue for TCR-based therapies across cancer types. Citation Format: Darwin W. Kwok, Takahide Nejo, Nicholas Stevers, Lee H. Chen, Kaori Okada, Chibo Hong, Gary K. Chan, Akane Yamamichi, Aidan Du, Maggie Colton, James Woo, Joseph Costello, Hideho Okada. Targetable public neoantigens are generated by tumor-wide cancer-specific splicing events in gliomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3152.

Background Chemotherapy improves overall survival after surgery and radiotherapy for newly diagnosed high-risk IDH-mutant low-grade gliomas (LGGs), but a proportion of patients treated with temozolomide (TMZ) will develop recurrent tumors with TMZ-induced hypermutation. We aimed to determine the prevalence of TMZ-induced hypermutation at recurrence and prognostic implications. Methods We sequenced recurrent tumors from 82 patients with initially low-grade IDH-mutant gliomas who underwent reoperation and correlated hypermutation status with grade at recurrence and subsequent clinical outcomes. Results Hypermutation was associated with high-grade disease at the time of reoperation (OR 12.0 95% CI 2.5-115.5, P = .002) and was identified at transformation in 57% of recurrent LGGs previously exposed to TMZ. After anaplastic (grade III) transformation, hypermutation was associated with shorter survival on univariate and multivariate analysis (HR 3.4, 95% CI 1.2-9.9, P = .024), controlling for tumor grade, subtype, age, and prior radiotherapy. The effect of hypermutation on survival after transformation was validated in an independent, published dataset. Hypermutated (HM) tumors were more likely to develop discontiguous foci of disease in the brain and spine (P = .003). To estimate the overall incidence of high-grade transformation among low-grade IDH-mutant tumors, data from a phase II trial of TMZ for LGG were analyzed. Eight-year transformation-free survival was 53.8% (95% CI 42.8-69.2), and 61% of analyzed transformed cases were HM. Conclusions TMZ-induced hypermutation is a common event in transformed LGG previously treated with TMZ and is associated with worse prognosis and development of discontiguous disease after recurrence. These findings impact tumor classification at recurrence, prognostication, and clinical trial design.

The factors driving therapy resistance in diffuse glioma remain poorly understood. To identify treatment-associated cellular and genetic changes, we analyzed RNA and/or DNA sequencing data from the temporally separated tumor pairs of 304 adult patients with isocitrate dehydrogenase (IDH)-wild-type and IDH-mutant glioma. Tumors recurred in distinct manners that were dependent on IDH mutation status and attributable to changes in histological feature composition, somatic alterations, and microenvironment interactions. Hypermutation and acquired CDKN2A deletions were associated with an increase in proliferating neoplastic cells at recurrence in both glioma subtypes, reflecting active tumor growth. IDH-wild-type tumors were more invasive at recurrence, and their neoplastic cells exhibited increased expression of neuronal signaling programs that reflected a possible role for neuronal interactions in promoting glioma progression. Mesenchymal transition was associated with the presence of a myeloid cell state defined by specific ligand-receptor interactions with neoplastic cells. Collectively, these recurrence-associated phenotypes represent potential targets to alter disease progression.

Tumor cell immortality is a fundamental hallmark of human cancers. Normally silenced during somatic cell differentiation, 90% of human tumors reactivate Telomerase Reverse Transcriptase (TERT) expression to achieve cellular immortality. TERT, the catalytic subunit of telomerase, complexes with the RNA template molecule TERC to maintain telomeres. Mutations in the TERT promoter (TERTp) are the most common non-coding mutation across all cancer types and the most frequent mutation within many cancers, such as IDH wildtype glioblastoma (GBM), Melanoma, and Bladder Cancer. TERTp mutations generate de novo E26 Transformation Specific (ETS) binding motifs that are spaced full helical turns from TERTp native ETS sites. Together the de novo and native ETS motifs specifically recruit the GABP tetrameric complex but not the GABP dimer. CRISPR-cas9 mediated insertion/deletion mutagenesis of the unique exon of GABP tetramer forming subunit, GABPB1L (B1L), reduces TERT transcriptional activity in a TERT promoter-mutation dependent manner. Here we show that GABPB1S (B1S), the GABP dimer restricted alternative isoform of GABPB1, is consistently and significantly increased following B1L reduction, a process we have determined to be driven by a conserved homeostatic mechanism whereby the GABP tetramer suppresses expression of one of its own components, GABPB1. In contrast to the native setting, in the absence of B1L the elevated B1S expression leads to dimer binding to the mutant TERTp and maintenance of TERT expression. Indeed, co-targeting B1L and B1S together, but not B1L alone, via CRISPR-cas9 knockout resulted in a near complete elimination of GABP recruitment to the TERTp and TERT expression, and lead to tumor cell death and eventual senescence in a telomere length dependent manner. Together, this data suggests a new model of the TERT-GABP axis involving the tetramer and dimer and highlights a new and potentially more potent therapeutic strategy to eliminate TERT expression and reverse tumor cell immortality.

The prospect of using immunotherapy for IDH-mutant LGGs has been daunting given the immune-poor microenvironment and low mutational burden. We hypothesized that LGG-targeting T cells might still be present at low frequency and with limited regional infiltration into the tumor. To improve sensitivity, we combined high-density multi-region tumor sampling with high-throughput neoantigen-T cell screening for a patient with WHO Grade II diffuse astrocytoma who eventually progressed, at second recurrence, to anaplastic astrocytoma. We performed maximal-anatomic sampling from 10 distinct regions of the tumor at the initial resection, as well as single sampling at first recurrence, for exome-based prediction of clonal and subclonal expressed neoantigens, RNAseq-based estimation of regional immune cell composition, and T cell receptor (TCR) beta deep sequencing. Based on our predictions, we then generated a barcoded library of patient-specific peptide-HLA multimers loaded with predicted neopeptides. Using this library, neoantigen-specific CD8 T cells were captured and isolated from patient peripheral blood and subjected to single cell TCR sequencing. We screened patient-derived peripheral blood drawn two years after initial resection and identified five T cell clones recognizing three LGG neoepitopes. Two neoepitopes were derived from truncal, tumor-wide mutations, including a truncating splice-site mutation in TP53 and a missense mutation in MRPL46. Each of these neoepitopes were recognized by two distinct TCRs, consistent with TCR convergence. A third neoepitope was derived from a subclonal MRPL46 mutation. Using the TCR beta sequence as a molecular barcode, TCRs specific to the 2 clonal neoepitopes, but not the subclonal neoepitope, were detectable in the glioma and blood at initial resection, as well as in the recurrent glioma. In summary, we demonstrate the existence and persistence of neoantigen-targeting T cells within the blood and tumor of an IDH-mutant LGG patient. These findings suggest a feasible methodology to develop personalized T cell-based immunotherapies for IDH-mutant LGGs.

Mutations in the TERT promoter are the single most common non-coding mutation in cancer and represent the genetic underpinnings of tumor cell immortality. Beyond the two most common point mutations, G228A and G250A, which selectively recruit the ETS factor GABP to activate TERT, the significance of other variants in the TERT promoter are unknown. We identified duplications of wildtype sequence within the core promoter region of TERT in 7 different cancer types that have strikingly similar features including size, insertion position, and inclusion of one of the native ETS motifs. Each duplication activates the TERT promoter to a similar level as G228A and G250A and is critically dependent on the insertion site. The GABP tetramer binds to the TERT duplicated promoter sequence by virtue of the native ETS motif and its duplicated version with precise spacing, and it is necessary for the transcriptional activation by all duplications tested. Spatiotemporal analysis in a multifocal glioblastoma shows the duplication is clonal and its activation of TERT is readily detectable at the single cell level and in bulk tumor tissue. We conclude that recurrent TERT promoter duplications of the native ETS sequence are functionally and mechanistically equivalent to the hotspot mutations that confer tumor cell immortality. The shared mechanism of these divergent somatic genetic alterations suggests a strong selective pressure for recruitment of the GABP tetramer to activate TERT.

Oligodendroglioma is a subtype of diffuse glioma defined by a mutation in the isocitrate dehydrogenase (IDH) genes and a co-deletion of chromosome arms 1p and 19q. These tumors primarily occur in adult patients in their third and fourth decade of life and are universally fatal due to an inevitable recurrence that follows a treatment regimen of surgical resection and an optional combination of alkylating chemotherapy and/or radiation therapy. While initially slow growing, recurrent tumors exhibit increasingly aggressive phenotypes that become progressively more difficult to treat with conventional therapy. Currently, the molecular mechanisms and cellular phenotypes that drive this recurrence remain unknown. To understand these factors, we assembled a cohort of matched initial and recurrent oligodendroglioma samples from over 100 patients and performed whole-genome sequencing and whole-exome sequencing on each of them. To link these molecular profiles to cell state changes, we additionally performed bulk and single-nucleus RNA-sequencing on a subset of these tumor pairs. In nearly 40% of alkylating chemotherapy-treated patients, recurrent tumors presented with hypermutation that corresponded with an increase in neoplastic cell proliferation. Additionally, while individual somatic alterations specific to recurrence were relatively rare, we observed a subset of tumors that acquired deletions in the cell cycle regulator CDKN2A following treatment with radiotherapy. Acquisition of either of these features associated with shorter patient survival and higher grade at recurrence, implicating cell cycle dysregulation as a mechanism of treatment resistance and increased tumor severity. Together, these results indicate that oligodendrogliomas evolve in a treatment-specific manner following chemo- and radiation therapy and highlight key pathways that can be targeted to delay the onset of recurrence.

TERT promoter mutation (TPM) is found in over 80% of IDH-wildtype glioblastomas (GBMs) and IDH-mutant oligodendrogliomas (ODs). As TPM-mediated reactivation of telomerase is tumor specific, TPM- tumor cells could be selectively targeted to reverse tumor cell immortalization, especially if this mutation is present throughout the tumor. Previous studies on the clonality of TPM present conflicting results and have relied on minimal tumor sampling. In this study, we use a 3D maximal tumor sampling approach to address this critical issue. We investigated TPM clonality in 253 primary and recurrent tumor samples, from 19 IDH-wildtype GBMs and 9 IDH-mutant ODs. An average of nine regions per tumor were intraoperatively biopsied in a manner that maximally represents tumor geography. The 3D location of each biopsy was then recorded and mapped back to the patient’s pre- and post-operative MRI, allowing 3D characterization of each tumor. For initial screening of TPM, PCR and Sanger sequencing were performed on all 253 samples. Deep amplicon sequencing (Amp-seq) was performed on 179 of these samples, including all those in which TPM was not detectable by Sanger. Each Amp-seq reaction included a TPM amplicon and an amplicon with IDH1 mutation for ODs or an amplicon spanning a chromosome 10 SNP to evaluate chromosome 10 loss in GBMs, which are early clonal events for these tumor types. Tumor purity, estimated by the FACETS algorithm applied to whole exome sequencing, was available for all 253 tumor samples and their patient-matched blood cell DNA. TERT expression was measured by RNA-seq (144 samples) and RNA Scope (49 samples). Correlations between TPM and clonal alterations and TPM and TERT expression were analyzed for significance. Of 253 tumor samples, tumor purity could be estimated for 210 (83.0%), and of these, 100% had TPM detected by Sanger and Amp-seq. For samples whose tumor purity could not be estimated (assigned a value of NA by FACETS), TPM was still detected in 36 of 43 (83.7%) samples. Variant allele frequencies (VAFs) of TPM showed high positive correlation with those of clonal alterations in GBMs (R = 0.85, p < 0.0001) and ODs (R = 0.86, p < 0.0001) as well as with tumor purity in both GBMs (R = 0.91, p < 0.0001) and ODs (R = 0.89, p < 0.0001). TPM VAF showed moderate positive correlation with TERT expression in GBMs (R = 0.52, p < 0.0001) and ODs (R = 0.62, p < 0.0001). TPM also showed a nonlinear relationship with TERT expression only in ODs. RNA Scope detected TERT expression in a subset of cells in both GBMs and ODs. To conclude, TPM is a tumor-wide, clonal mutation in primary and recurrent GBMs and ODs, making TPM tumor cells a possible therapeutic target. TPM VAF is moderately correlated with TERT expression, with a nonlinear relationship seen in ODs, suggesting other variables affect TERT expression. RNA Scope detected TERT expression at the single cell level in situ in GBMs and ODs and therefore could serve as a tumor cell-specific marker, although follow up studies are needed. Citation Format: Christina Leann Appin, Abigail K. Suwala, Stephanie Hilz, Radhika Mathur, Ivan V. Smirnov, Chibo Hong, Nicholas O. Stevers, Anny Shai, Albert Wang, Mitchel S. Berger, Susan M. Chang, Joanna J. Phillips, Joseph F. Costello. 3D whole tumor analysis of the TERT promoter mutation in gliomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3797.

Ninety percent of human tumors reactivate telomerase reverse transcriptase (TERT) to achieve cellular immortality. The route to immortality found in over 50 cancer types is mutation of the TERT promoter (TERTp) mutations. Studies in glioblastoma have shown the two most common TERTp mutations, G228A and G250A, reactivate TERT expression through de novo E26 transformation specific (ETS) site generation that, in tandem with a native TERTp ETS site, recruits the tetrameric GA-binding protein (GABP) complex formed by two GABP alpha (GABPA) subunits and two tetramer specific GABPB1 (B1) subunits, GABPB1L (B1L). While ten additional TERTp mutations have been described across cancer types, functional studies have been limited to select mutations and cancer types. We found that each de novo ETS motif generating TERTp mutation increases TERTp transcriptional activity in a GABP-dependent manner. Furthermore, our pan-cancer analysis demonstrated selective enrichment of GABPA at the mutant TERTp across 14 cancer types and further determined that this recruitment is necessary for TERT expression. Prior investigations found that CRISPR-cas9 mediated mutagenesis of the B1 ninth exon, specific to B1L, reduces TERT transcriptional activity in a TERTp-mutation dependent manner. We explored B1L targeting in the reversal of tumor cell immortality in several cancers. We found that the B1 dimer restricted isoform, GABPB1S (B1S), is consistently and significantly increased following B1L reduction, which we determined to be the consequence of a disengagement of a GABP tetramer mediated negative feedback loop acting on the B1 promoter. Interestingly, this upregulated B1S expression enables GABP dimer binding and activation of the mutant TERTp thereby compensating for B1L loss. Indeed, knockout of B1L and B1S together resulted in a near complete elimination of GABPA TERTp recruitment and a striking reduction in TERT expression, leading to telomere shortening and ultimately resulting in tumor cell death and senescence. This data identifies GABP as the master regulator of the mutant TERTp across cancer type and de novo ETS mutation type and suggests a new model of the GABP-TERT axis involving both the GABP tetramer, and the GABP dimer. Citation Format: Nicholas O. Stevers, Carter Barger, Olivia Lenzo, Chibo Hong, Katarzyna Soczek, Samuel H. Wu, Andrew M. McKinney, Abigail Suwala, Jennifer A. Doudna, Joseph F. Costello. The pan-cancer regulator of the mutant TERT promoter and a new model of the GABP-TERT axis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2302.

Treatment failure in glioblastoma, the most common and lethal primary brain tumor in adults, is often attributed to intratumoral heterogeneity as it fosters tumor evolution and selection of resistant clones. Sources of intratumoral heterogeneity may include genomic alterations, such as single-nucleotide and structural variants, and epigenomic alterations, such as changes in chromatin structure and transcriptional regulation. The relative extent and functional significance of these contributors to intratumoral heterogeneity in glioblastoma have yet to be elucidated. In collaboration with neurosurgeons and neuro-imaging experts, we have established a novel approach towards characterizing intratumoral heterogeneity in three-dimensional (3D) space. For patients undergoing tumor resection, we utilize 3D surgical neuronavigation to safely acquire ~10 samples representing maximal anatomical diversity. Samples are mapped by 3D spatial coordinates and integrated with patient MRI scans for 360º visualization of sample location in context of the brain. We have now conducted whole-exome sequencing (Exome-Seq), assay for transposase-accessible chromatin (ATAC-Seq), and RNA-sequencing (RNA-Seq) for 83 spatially mapped samples obtained from 8 patients with primary IDH-WT glioblastoma. Integrative data analysis provides unprecedented insight into sources of intratumoral heterogeneity in glioblastoma and their 3D spatial patterning. We find that tumor cells show aberrant transcription factor activity, activation of fetal brain programs, and variable chromatin accessibility at CTCF-binding loop anchors indicating plasticity in higher-order chromatin structure. Chromosome conformation capture analysis by Hi-C extends these findings and reveals intratumoral differences in long-range chromatin interactions due to structural variants and enhancer hijacking. Further, we use deconvolution to identify microenvironmental contributors to epigenomic intratumoral heterogeneity including neural, glial, and immune populations. We define chromatin signatures associated with microenvironmental cell types and states, revealing their 3D spatial patterning, and validate these findings by single-cell ATAC-Seq. Our work thus establishes both tumor and microenvironmental sources of intratumoral heterogeneity in glioblastoma, revealing their chromatin programs and 3D spatial patterning. As a resource for further investigation, we have developed an interactive data sharing platform that enables visualization of brain tumor intratumoral heterogeneity in 360 degrees. Citation Format: Radhika Mathur, Qixuan Wang, Patrick Schupp, Ana Nikolic, Takafumi Yamaguchi, Stephanie Hilz, Chibo Hong, Ivan Smirnov, Marisa LaFontaine, Joanna Phillips, Susan Chang, Yan Li, Janine Lupo, Paul Boutros, Marco Gallo, Michael Oldham, Feng Yue, Joseph Costello. 3D spatial sampling and integrated omics reveal sources and patterning of intratumoral heterogeneity in glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3621.

Treatment failure in glioblastoma is often attributed to intratumoral heterogeneity (ITH), which fosters tumor evolution and generation of therapy-resistant clones. While ITH in glioblastoma has been well-characterized at the genomic and transcriptomic levels, the extent of ITH at the epigenomic level and its biological and clinical significance are not well understood. In collaboration with neurosurgeons, neuropathologists, and biomedical imaging experts, we have established a novel topographical approach towards characterizing epigenomic ITH in three-dimensional (3-D) space. We utilize pre-operative MRI scans to define tumor volume and then utilize 3-D surgical neuro-navigation to intra-operatively acquire 10+ samples representing maximal anatomical diversity. The precise spatial location of each sample is mapped by 3-D coordinates, enabling tumors to be visualized in 360-degrees and providing unprecedented insight into their spatial organization and patterning. For each sample, we conduct assay for transposase-accessible chromatin using sequencing (ATAC-Seq), which provides information on the genomic locations of open chromatin, DNA-binding proteins, and individual nucleosomes at nucleotide resolution. We additionally conduct whole-exome sequencing and RNA sequencing for each spatially mapped sample. Integrative analysis of these datasets reveals distinct patterns of chromatin accessibility within glioblastoma tumors, as well as their associations with genetically defined clonal expansions. Our analysis further reveals how differences in chromatin accessibility within tumors reflect underlying transcription factor activity at gene regulatory elements, including both promoters and enhancers, and drive expression of particular gene expression sets, including neuronal and immune programs. Collectively, this work provides the most comprehensive characterization of epigenomic ITH to date, establishing its importance for driving tumor evolution and therapy resistance in glioblastoma. As a resource for further investigation, we have provided our datasets on an interactive data sharing platform – The 3D Glioma Atlas – that enables 360-degree visualization of both genomic and epigenomic ITH.

Aim: To provide a comprehensive understanding of gene regulatory networks in the developing human brain and a foundation for interpreting pathogenic deregulation. Materials & methods: We generated reference epigenomes and transcriptomes of dissected brain regions and primary neural progenitor cells (NPCs) derived from cortical and ganglionic eminence tissues of four normal human fetuses. Results: Integration of these data across developmental stages revealed a directional increase in active regulatory states, transcription factor activities and gene transcription with developmental stage. Consistent with differences in their biology, NPCs derived from cortical and ganglionic eminence regions contained common, region specific, and gestational week specific regulatory states. Conclusion: We provide a high-resolution regulatory network for NPCs from different brain regions as a comprehensive reference for future studies.

Background Emerging data suggest that a subset of patients with diffuse isocitrate dehydrogenase (IDH)-mutant low-grade glioma (LGG) who receive adjuvant temozolomide (TMZ) recur with hypermutation in association with malignant progression to higher-grade tumors. It is currently unclear why some TMZ-treated LGG patients recur with hypermutation while others do not. MGMT encodes O6-methylguanine-DNA methyltransferase, a DNA repair protein that removes cytotoxic and potentially mutagenic lesions induced by TMZ. Here, we hypothesize that epigenetic silencing of MGMT by promoter methylation facilitates TMZ-induced mutagenesis in LGG patients and contributes to development of hypermutation at recurrence. Methods We utilize a quantitative deep sequencing assay to characterize MGMT promoter methylation in 109 surgical tissue specimens from initial tumors and post-treatment recurrences of 37 TMZ-treated LGG patients. We utilize methylation arrays to validate our sequencing assay, RNA sequencing to assess the relationship between methylation and gene expression, and exome sequencing to determine hypermutation status. Results Methylation level at the MGMT promoter is significantly higher in initial tumors of patients that develop hypermutation at recurrence relative to initial tumors of patients that do not (45.7% vs 34.8%, P = 0.027). Methylation level in initial tumors can predict hypermutation at recurrence in univariate models and multivariate models that incorporate patient age and molecular subtype. Conclusions These findings reveal a mechanistic basis for observed differences in patient susceptibility to TMZ-driven hypermutation. Furthermore, they establish MGMT promoter methylation level as a potential biomarker to inform clinical management of LGG patients, including monitoring and treatment decisions, by predicting risk of hypermutation at recurrence.

The design and evaluation of immunotherapies in IDH-mutant lower grade gliomas (LGG) is hindered by a poor understanding of the LGG T-cell repertoire. We present data on the temporal evolution, intratumoral spatial distribution, and prognostic value of the T-cell repertoire in IDH-mutant LGGs. We performed immunogenomic profiling using T-cell receptor beta-chain sequencing of 163 glioma and peripheral blood samples from 33 immunotherapy-naive glioma patients (22 astrocytomas, 11 oligodendrogliomas). T-cell repertoire evolution was analyzed in a subset of 26 patients (69 samples) with matched primary (WHO grade II) and recurrent (WHO grade II-IV) glioma samples. T-cell repertoire diversity was defined as the number of unique T-cell clonotypes by V-gene, J-gene, and CDR3 nucleotide sequences. Malignant transformed (Grade III or IV) recurrent gliomas demonstrated increased T-cell repertoire diversity compared to their patient-matched primary tumors (p=0.0023), but grade II recurrences did not show the same increased diversity (p=0.26). This increase in T-cell repertoire diversity was greater in patients who underwent transformation in the context of TMZ-associated hypermutation compared to spontaneously transformed counterparts (p=0.035). In grade II primary astrocytomas (n=17), T-cell repertoire diversity above the median (186 unique T-cell clonotypes per sample) was associated with worse transformation-free (HR=4.2, p=0.045) and overall survival (HR=6.4, p=0.025). Next, we evaluated intratumoral immune heterogeneity in 7 patients by sampling from up to 10 distinct and maximally-separated intratumoral sites per LGG (64 samples). Eighty-two to 96% of unique clonotypes within a given tumor were present only within a single sampled site. Despite this heterogeneity, six LGG patients harbored T-cell clonotypes present tumor-wide across all sampled sites within a given tumor. Ten of 24 (42%) tumor-wide T-cell clonotypes were enriched in the glioma compared to matched peripheral blood, suggesting glioma-specificity. Taken together, T-cell receptor profiling in LGGs may have utility both as a prognostic biomarker and to identify glioma-specific T-cells.

TERT promoter mutation (TPM), found in over 80% of IDH-wildtype glioblastomas (GBMs) and oligodendrogliomas, leads to reactivation of telomerase and consequently tumor cell immortalization, which is potentially reversible. TERT could therefore serve as an effective target in treating tumors with TPM, if TPM is present throughout the tumor. Previous studies using a single sample or minimal sampling per tumor have shown potentially conflicting results, suggesting TPM is clonal in some cases and subclonal in others. Here we use spatially mapped tumor samples representing maximal tumor sampling to address this critical issue. Sanger sequencing was performed on 311 newly diagnosed and recurrent tumor samples from 19 IDH-wildtype GBMs and 10 oligodendrogliomas. To validate Sanger sequencing and resolve potentially ambiguous samples, deep amplicon sequencing was performed on 164 samples. To determine tumor purity and TERT expression levels, whole exome sequencing (164 samples) and RNA-Seq (129 samples) data sets were analyzed computationally. Sanger and amplicon sequencing showed that TPM was present in 305 of 311 samples (98.1%). TPM was not detected in 6 samples which had tumor purity estimates too low to be accurately determined by FACETS and lacked evidence of any driver mutation. Variant allele frequencies (VAFs) of TPM showed high positive correlation with those of clonal alterations in GBMs (r(90) = .93, p < .0001) and oligodendrogliomas (r(48) = .96, p < .0001). TPM VAFs also showed high positive correlation with tumor purity in both GBMs (r(112) = .92, p < .0001) and oligodendrogliomas (r(48) = .89, p < .0001). TPM VAF showed a moderate positive correlation with TERT expression in GBMs (r(78) = .40, p < .001) and oligodendrogliomas (r(47) = .49, p < .001). Therefore, TPM is a tumor-wide, clonal mutation in both newly diagnosed and recurrent GBMs and oligodendrogliomas. TPM VAF is moderately correlated with TERT expression.

Treatment failure in glioblastoma is often attributed to intratumoral heterogeneity (ITH), which fosters tumor evolution and selection of therapy-resistant clones. While genomic alterations are known contributors to ITH, emerging studies highlight functional roles for epigenomic ITH which integrates differentiation status, stochastic events, and microenvironmental inputs. Here, we have established a novel platform for integrative characterization of genomic and epigenomic ITH of glioblastoma in three-dimensional (3-D) space. In collaboration with neurosurgeons and biomedical imaging experts, we utilize 3-D surgical neuro-navigation to safely acquire ~10 tumor samples per patient representing maximal anatomical diversity. We conduct whole-exome sequencing, RNA sequencing, and assay for transposase-accessible chromatin using sequencing (ATAC-Seq) on each sample. The spatial location of each sample is mapped by its 3-D coordinates, allowing 360-degree visualization of genomic and epigenomic ITH for each patient. We demonstrate this approach on 8 patients with primary IDH-WT glioblastoma (83 spatially mapped samples), providing unprecedented insight into their spatial organization at the genomic and epigenomic levels. We link genetically defined tumor subclones to patterns of open chromatin and gene regulation, revealing underlying transcription factor binding at active promoters and enhancers. We also identify ITH in whole-genome doubling and focal oncogene amplification events in multiple patients, which we then link with epigenomic ITH. Further, to study microenvironmental inputs and their contribution to epigenomic ITH, we conduct deconvolution of RNA sequencing and ATAC-Seq data by analyzing feature co-variation. We resolve the 3-D spatial organization of immune, neural, and other nontumor cell types present in glioblastoma, characterizing their functional states and interactions with tumor cells. This work provides the most comprehensive spatial characterization of genomic and epigenomic ITH to date in glioblastoma. As a resource for further investigation, we have developed an interactive data sharing platform – The 3D Glioma Atlas – that enables 360-degree visualization of both genomic and epigenomic ITH.

The human placenta and its specialized cytotrophoblasts rapidly develop, have a compressed lifespan, govern pregnancy outcomes, and program the offspring's health. Understanding the molecular underpinnings of these behaviors informs development and disease. Profiling the extraembryonic epigenome and transcriptome during the 2nd and 3rd trimesters revealed H3K9 trimethylation overlapping deeply DNA hypomethylated domains with reduced gene expression and compartment-specific patterns that illuminated their functions. Cytotrophoblast DNA methylation increased, and several key histone modifications decreased across the genome as pregnancy advanced. Cytotrophoblasts from severe preeclampsia had substantially increased H3K27 acetylation globally and at genes that are normally downregulated at term but upregulated in this syndrome. In addition, some cases had an immature pattern of H3K27ac peaks, and others showed evidence of accelerated aging, suggesting subtype-specific alterations in severe preeclampsia. Thus, the cytotrophoblast epigenome dramatically reprograms during pregnancy, placental disease is associated with failures in this process, and H3K27 hyperacetylation is a feature of severe preeclampsia.

BACKGROUND Previous studies of solid tumors have been restricted in their ability to map how heterogeneous cell populations evolved within the tumor in three-dimensional (3D) space due to insufficient sampling, typically one sample per tumor, and a lack of knowledge of where within the tumor the sample was obtained. Knowledge of the extensivity of heterogeneity and how it is spatially distributed is crucial for assessing whether a therapeutic target is truly tumor-wide, and for exploring how mutations relate to heterogeneity in the local microenvironment. METHODS We developed a novel platform to integrate and visualize in 3D multi-omics data generated from each of 8–10 spatially mapped samples per tumor. Together, the 171 samples collected using this approach from 16 adult diffuse glioma at diagnosis and recurrence form a novel resource – the 3D Glioma Atlas. RESULTS By maximally sampling the tumor geography without excluding samples based on low cancer cell fraction (CCF), we identify a subpopulation of glioblastoma with pervasively lower CCF likely excluded by other atlases, such as the TCGA, that used stringent CCF cutoffs. Exome sequencing of 3D-mapped samples from lower-grade tumors revealed that clonal expansions are typically spatially segregated, implying minimal tumor-wide intermixing of genetically heterogenous cells. Heterogeneity is less spatially segregated for faster-growing high-grade tumors, suggesting that cell populations expand in these tumors differently. Recurrent low-grade tumors have greater intratumoral mutational heterogeneity than newly diagnosed tumors, though this did not translate into greater dissimilarity in gene expression profiles for recurrent tumors, suggesting minimal functional impact of this additional mutational diversity on gene expression. CONCLUSIONS The delineation of spatial patterns of heterogeneity that our work provides enables more informed interpretation of biopsies and greater insight into the factors shaping intratumoral variation of gene expression patterns. Ongoing work is exploring the spatial patterning of amplification events and the tumor microenvironment.

// Pavithra Viswanath 1 , Chloe Najac 1 , Jose L. Izquierdo-Garcia 1 , Aleksandr Pankov 2 , Chibo Hong 2 , Pia Eriksson 1 , Joseph F. Costello 2 , Russell O. Pieper 2 , Sabrina M. Ronen 1 1 Department of Radiology and Biomedical Imaging, University of California San Francisco, San Francisco, CA 94143, USA 2 Department of Neurological Surgery, Helen Diller Research Center, University of California San Francisco, San Francisco, CA 94143, USA Correspondence to: Sabrina Ronen, e-mail: sabrina.ronen@ucsf.edu Keywords: MCT1, MCT4, mutant IDH1, metabolic reprogramming, magnetic resonance spectroscopy Received: February 05, 2016 Accepted: April 10, 2016 Published: April 26, 2016 ABSTRACT 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 13 C magnetic resonance spectroscopy and compared the conversion of hyperpolarized [1- 13 C]-pyruvate to [1- 13 C]-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.

2506 Background: Temozolomide, a commonly used alkylating agent to treat gliomas, can induce somatic hypermutation. The prevalence and clinical implications of this phenomenon are not well characterized. Methods: We used targeted and whole exome sequencing from a cohort of 82 patients with recurrent IDH-mutant low grade gliomas undergoing re-operation to evaluate the prevalence as well as the clinical implications of hypermutation. Results: Hypermutation was identified at transformation in 57% of recurrent gliomas exposed to Temozolomide, 94% of which were transformed to higher WHO grades. All patients who developed hypermutation were exposed to Temozolomide. Hypermutation was associated with transformation to higher WHO grade (OR 12.0 95% CI 2.5 – 115.5, p = 0.002) and shorter survival after transformation (HR 2.1, 95% CI 1.1-4.0, p = 0.018) compared with non-hypermutated transformed tumors, controlling for grade, molecular subtype, age, and prior radiotherapy. Patients with transformation to glioblastoma had poor survival regardless of hypermutation (p = 0.78). Hypermutated tumors were associated with development of discontiguous disease at a significantly higher frequency (p = 0.003), including four cases with spinal dissemination. Conclusions: TMZ-induced hypermutation is associated with high grade transformation, unique patterns of dissemination and shortened survival after transformation. Next generation sequencing should be considered in this patient population. These data have important implications for the management of newly diagnosed and recurrent IDH-mutant low grade gliomas.

Low-grade gliomas (LGGs), which include grade II astrocytoma and grade II oligodendroglioma, inevitably recur despite aggressive treatment with surgery, and sometimes, with radiation and the chemotherapeutic agent temozolomide (TMZ). The clinical benefit of TMZ in LGG is unclear, and a subset of TMZ-treated LGGs recur with hypermutation in association with malignant progression to high-grade tumors. It is currently unknown why some TMZ-treated LGGs recur with hypermutation while others do not. O6-methylguanine-DNA methyltransferase (MGMT) is a DNA repair protein that reverses mutagenic lesions induced by TMZ. The amount of MGMT protein in a cell is regulated at the epigenetic level by promoter methylation. Here, we hypothesize that epigenetic silencing of MGMT by promoter methylation facilitates TMZ-induced mutagenesis and contributes to the development of hypermutation. We demonstrate in a cohort of 37 TMZ-treated patients with an initial diagnosis of IDH-mutant LGG that methylation level of the MGMT promoter in initial untreated tumors is significantly associated with hypermutation at recurrence. We also confirm our previous finding that methylation level of the MGMT promoter in recurrent hypermutated tumors is higher than in recurrent tumors that are not hypermutated. These results provide a plausible mechanistic basis for observed differences in propensity of TMZ-treated LGG patients to develop hypermutation at recurrence. Furthermore, they establish the potential of MGMT promoter methylation level to inform treatment decisions in the clinic for patients with newly diagnosed LGG.

Most cancers result from multiple clonal expansion events, resulting in tumors comprised of many subclones. The coexistence of these genetically distinct populations within a single tumor is thought to underlie the failure of most current treatments. Ascertaining how diverse populations of cells evolve within tumors, and the functional implications of this diversity, are the next steps to developing more informed treatments that overcome the barrier currently posed by intratumoral heterogeneity (ITH). Previous studies of adult diffuse glioma have been restricted in their ability to provide a detailed understanding of the spatial evolution of solid tumors due to insufficient sampling, typically just one biopsy per patient, and a lack of knowledge of where within the heterogeneous tumor the sample was obtained. We have developed a unique sample collection protocol, in which we record the precise, three-dimensional (3D) coordinate for each of 10 spatially-distinct samples obtained from each patient’s tumor. Whole exome data from IDH1 mutant grade II astrocytomas analyzed in this manner have revealed that expansion events are typically spatially segregated, suggesting minimal tumor-wide intermixing of subclonal cell populations in these low grade gliomas. Expression profiling of one particular case has furthermore revealed ITH of the immune microenvironment, uncovering a region of increased PD-L1 expression and IFN-gamma signaling. Ongoing analysis is focused on integrating genetic and expression profiles from these and additional spatially-mapped cases. This integrative analysis aims to dissect the relative contributions of subclonal variation and location within the tumor on gene expression patterns and ITH of the immune microenvironment. The delineation of the spatial patterns of ITH our work provides will enable more informed selection of biopsy targets in the future. It also lends greater insight into the factors shaping intratumoral variation of gene expression patterns and the immune microenvironment in adult diffuse glioma.

Background Rare multicentric lower-grade gliomas (LGGs) represent a unique opportunity to study the heterogeneity among distinct tumor foci in a single patient and to infer their origins and parallel patterns of evolution. Methods In this study, we integrate clinical features, histology, and immunohistochemistry for 4 patients with multicentric LGG, arising both synchronously and metachronously. For 3 patients we analyze the phylogeny of the lesions using exome sequencing, including one case with a total of 8 samples from the 2 lesions. Results One patient was diagnosed with multicentric isocitrate dehydrogenase 1 (IDH1) mutated diffuse astrocytomas harboring distinct IDH1 mutations, R132H and R132C; the latter mutation has been associated with Li-Fraumeni syndrome, which was subsequently confirmed in the patient's germline DNA and shown in additional cases with The Cancer Genome Atlas data. In another patient, phylogenetic analysis of synchronously arising grade II and grade III diffuse astrocytomas demonstrated a single shared mutation, IDH1 R132H, and revealed convergent evolution via non-overlapping mutations in ATRX and TP53. In 2 cases, there was divergent evolution of IDH1-mutated and 1p/19q-codeleted oligodendroglioma and IDH1-mutated and 1p/19q-intact diffuse astrocytoma, occurring synchronously in one case and metachronously in a second. Conclusions Each tumor in multicentric LGG cases may arise independently or may diverge very early in their development, presenting as genetically and histologically distinct tumors. Comprehensive sampling of these lesions can therefore significantly alter diagnosis and management. Additionally, somatic IDH1 R132C mutation in either multicentric or solitary LGG identifies unsuspected germline TP53 mutation, validating the limited number of published cases.

TERT promoter mutations reactivate telomerase, allowing for indefinite telomere maintenance and enabling cellular immortalization. These mutations specifically recruit the multimeric ETS factor GABP, which can form two functionally independent transcription factor species: a dimer or a tetramer. We show that genetic disruption of GABPβ1L (β1L), a tetramer-forming isoform of GABP that is dispensable for normal development, results in TERT silencing in a TERT promoter mutation-dependent manner. Reducing TERT expression by disrupting β1L culminates in telomere loss and cell death exclusively in TERT promoter mutant cells. Orthotopic xenografting of β1L-reduced, TERT promoter mutant glioblastoma cells rendered lower tumor burden and longer overall survival in mice. These results highlight the critical role of GABPβ1L in enabling immortality in TERT promoter mutant glioblastoma., This work was supported by a generous gift from the Dabbiere family (J.F.C.), the Hana Jabsheh Research Initiative (J.F.C.), NIH grant NCI P50CA097257 (J.F.C. and J.A.D.), NCI P01CA118816-06 (J.F.C.), T32 GM008568 and T32 CA151022 (A.M.), and NCI R01CA163336 (J.S.S.), and the Sontag Foundation Distinguished Scientist Award (J.S.S.). C.F. is supported by a US NIH K99/R00 Pathway to Independence Award (K99GM118909) from the National Institute of General Medical Sciences. Additional support was provided by Fundação para a Ciência e Tecnologia SFRH/BD/88220/2012 (A.X.-M.) and IF/00601/2012 (B.M.C.). J.A.D. is an investigator of the Howard Hughes Medical Institute., info:eu-repo/semantics/publishedVersion

Contains fulltext : 155226.pdf (Publisher’s version ) (Closed access) Temozolomide (TMZ) increases the overall survival of patients with glioblastoma (GBM), but its role in the clinical management of diffuse low-grade gliomas (LGG) is still being defined. DNA hypermethylation of the O (6) -methylguanine-DNA methyltransferase (MGMT) promoter is associated with an improved response to TMZ treatment, while inactivation of the DNA mismatch repair (MMR) pathway is associated with therapeutic resistance and TMZ-induced mutagenesis. We previously demonstrated that TMZ treatment of LGG induces driver mutations in the RB and AKT-mTOR pathways, which may drive malignant progression to secondary GBM. To better understand the mechanisms underlying TMZ-induced mutagenesis and malignant progression, we explored the evolution of MGMT methylation and genetic alterations affecting MMR genes in a cohort of 34 treatment-naive LGGs and their recurrences. Recurrences with TMZ-associated hypermutation had increased MGMT methylation compared to their untreated initial tumors and higher overall MGMT methylation compared to TMZ-treated non-hypermutated recurrences. A TMZ-associated mutation in one or more MMR genes was observed in five out of six TMZ-treated hypermutated recurrences. In two cases, pre-existing heterozygous deletions encompassing MGMT, or an MMR gene, were followed by TMZ-associated mutations in one of the genes of interest. These results suggest that tumor cells with methylated MGMT may undergo positive selection during TMZ treatment in the context of MMR deficiency.

SummaryThe 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.

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

Transposable element (TE)-derived sequences comprise half of the human genome and DNA methylome and are presumed to be densely methylated and inactive. Examination of genome-wide DNA methylation status within 928 TE subfamilies in human embryonic and adult tissues identified unexpected tissue-specific and subfamily-specific hypomethylation signatures. Genes proximal to tissue-specific hypomethylated TE sequences were enriched for functions important for the relevant tissue type, and their expression correlated strongly with hypomethylation within the TEs. When hypomethylated, these TE sequences gained tissue-specific enhancer marks, including monomethylation of histone H3 at lysine 4 (H3K4me1) and occupancy by p300, and a majority exhibited enhancer activity in reporter gene assays. Many such TEs also harbored binding sites for transcription factors that are important for tissue-specific functions and showed evidence of evolutionary selection. These data suggest that sequences derived from TEs may be responsible for wiring tissue type-specific regulatory networks and may have acquired tissue-specific epigenetic regulation.

Human organic cation transporter 3 (OCT3 and SLC22A3) mediates the uptake of many important endogenous amines and basic drugs in a variety of tissues. OCT3 is identified as one of the important risk loci for prostate cancer, and is markedly underexpressed in aggressive prostate cancers. The goal of this study was to identify genetic and epigenetic factors in the promoter region that influence the expression level of OCT3. Haplotypes that contained the common variants, g.-81G>delGA (rs60515630) (minor allele frequency 11.5% in African American) and g.-2G>A (rs555754) (minor allele frequency>30% in all ethnic groups) showed significant increases in luciferase reporter activities and exhibited stronger transcription factor-binding affinity than the haplotypes that contained the major alleles. Consistent with the reporter assays, OCT3 messenger RNA expression levels were significantly higher in Asian (P

While patients with primary glioblastoma (GBM) have benefited from increased overall survival due to treatment with the alkylating chemotherapy temozolomide (TMZ), the benefit of TMZ for patients with diffuse low-grade gliomas (LGG) remains unknown. We previously demonstrated that among ten TMZ-treated patients with initial LGG (IDH1 mutant, 1p19q intact), treatment induced hypermutation in the recurrent tumors of six patients, all six of whom underwent malignant progression to secondary GBM (sGBM) (Johnson & Mazor, Science, 2014). To further explore the relationships among TMZ treatment, hypermutation and malignant progression, we studied tumor evolution in an expanded cohort of 35 TMZ-treated patients by exome-sequencing of paired IDH1/2-mutant LGGs and their post-TMZ recurrences. This expanded cohort allowed us to contrast patterns of hypermutation between 1p19q intact versus co-deleted subgroups. Hypermutated recurrences were common in both subgroups and consistently underwent malignant progression to the highest grade while acquiring somatic mutations in mismatch repair (MMR) genes. In contrast to the 1p19q intact subgroup in which TP53 mutations were present at diagnosis, the hypermutated recurrences of the 1p19q co-deleted subgroup frequently acquired TMZ-associated mutations in TP53. This suggests that TP53 mutations may cooperate with MMR mutations to convert TMZ-induced cytotoxicity to mutagenicity. Once hypermutated clones arose, we found that they dominated subsequent recurrences. In one case, a hypermutated sGBM was resected, but soon recurred with leptomeningeal spread, at which point hypermutated clones were found in the spinal cord. We examined the clonality of hypermutated tumors through exome-sequencing of multiple spatially distinct samples, from which we estimate that hypermutated recurrent tumors can originate from a single initiating cell. Together, these findings further our understanding of the molecular features and clinical behavior of hypermutated clones, raising concerns about the potent mutagenic activity of TMZ and other alkylating agents in LGG patients.

The human telomerase is an enzyme critical for maintaining telomere length and chromosomal stability in stem cells. The transcriptional regulation of the telomerase reverse transcriptase (TERT) gene, encoding the catalytic subunit of telomerase, is a rate-limiting step in modulating telomerase activity. Although normally silenced in somatic cells, TERT is aberrantly expressed in 90% of aggressive cancers, highlighting this event as a hallmark of tumorigenesis. Reactivating telomerase helps cells with finite lifespan to achieve limitless proliferative potential and bypass cellular senescence induced by DNA replication-associated telomere shortening. Recently discovered non-coding mutations in the TERT promoter are among the most common genetic alterations observed across multiple cancer types. Specifically, one of two positions, G228A or G250A, is mutated in 20% of medulloblastomas, 44% of hepatocellular carcinomas (HCC), 66% of urothelial carcinomas of the bladder, 71% of melanomas, 78% of oligodendrogliomas, and 83% of primary glioblastomas (GBM). The common mutation sites may create de-novo ETS family transcription factor binding sites, but the precise mechanism of how these mutations confer increased TERT expression has remained elusive. Here, a focused siRNA screen of the many ETS transcription factors expressed in GBM identifies GABPA as the single ETS factor to selectively regulate the mutant but not the wild type TERT promoter. Single molecule binding assays and ChIP-qPCR analysis reveal that GABPA is exclusively recruited to the mutant allele in vitro and in vivo respectively. Furthermore, this allelic recruitment is consistent across four tested cancer types, highlighting a shared mechanism underlying mutant TERT promoter activation. Tandem flanking native ETS motifs critically cooperate with these mutations to activate TERT, likely by facilitating GABP heterotetramer binding. GABP thus directly links TERT promoter mutations to aberrant expression across multiple cancers.

Although it is known that the methylation of DNA in 5' promoters suppresses gene expression, the role of DNA methylation in gene bodies is unclear. In mammals, tissue- and cell type-specific methylation is present in a small percentage of 5' CpG island (CGI) promoters, whereas a far greater proportion occurs across gene bodies, coinciding with highly conserved sequences. Tissue-specific intragenic methylation might reduce, or, paradoxically, enhance transcription elongation efficiency. Capped analysis of gene expression (CAGE) experiments also indicate that transcription commonly initiates within and between genes. To investigate the role of intragenic methylation, we generated a map of DNA methylation from the human brain encompassing 24.7 million of the 28 million CpG sites. From the dense, high-resolution coverage of CpG islands, the majority of methylated CpG islands were shown to be in intragenic and intergenic regions, whereas less than 3% of CpG islands in 5' promoters were methylated. The CpG islands in all three locations overlapped with RNA markers of transcription initiation, and unmethylated CpG islands also overlapped significantly with trimethylation of H3K4, a histone modification enriched at promoters. The general and CpG-island-specific patterns of methylation are conserved in mouse tissues. An in-depth investigation of the human SHANK3 locus and its mouse homologue demonstrated that this tissue-specific DNA methylation regulates intragenic promoter activity in vitro and in vivo. These methylation-regulated, alternative transcripts are expressed in a tissue- and cell type-specific manner, and are expressed differentially within a single cell type from distinct brain regions. These results support a major role for intragenic methylation in regulating cell context-specific alternative promoters in gene bodies.

A mutant promoter's partner in crime Telomerase is an enzyme that maintains the ends of chromosomes. TERT , the gene coding for the enzyme's catalytic subunit, is not expressed in healthy somatic cells, but its expression is reactivated in the majority of human cancers. The resultant high levels of telomerase help cancer cells survive and multiply. Recurrent mutations in the promoter region of TERT are associated with high telomerase levels in multiple cancer types. Bell et al. show that a specific transcription factor called GABP is selectively recruited to the mutant form of the TERT promoter, which activates TERT gene expression Science , this issue p. 1036

While significant effort has been dedicated to the characterization of epigenetic changes associated with prenatal differentiation, relatively little is known about the epigenetic changes that accompany post-natal differentiation where fully functional differentiated cell types with limited lifespans arise. Here we sought to address this gap by generating epigenomic and transcriptional profiles from primary human breast cell types isolated from disease-free human subjects. From these data we define a comprehensive human breast transcriptional network, including a set of myoepithelial- and luminal epithelial-specific intronic retention events. Intersection of epigenetic states with RNA expression from distinct breast epithelium lineages demonstrates that mCpG provides a stable record of exonic and intronic usage, whereas H3K36me3 is dynamic. We find a striking asymmetry in epigenomic reprogramming between luminal and myoepithelial cell types, with the genomes of luminal cells harbouring more than twice the number of hypomethylated enhancer elements compared with myoepithelial cells., Epigenetic changes associated with post-natal differentiation have been characterized. Here the authors generate epigenomic and transcriptional profiles from primary human breast cells, providing insights into the transcriptional and epigenetic events that define post-natal cell differentiation in vivo.

The role of intermediate methylation states in DNA is unclear. Here, to comprehensively identify regions of intermediate methylation and their quantitative relationship with gene activity, we apply integrative and comparative epigenomics to 25 human primary cell and tissue samples. We report 18,452 intermediate methylation regions located near 36% of genes and enriched at enhancers, exons and DNase I hypersensitivity sites. Intermediate methylation regions average 57% methylation, are predominantly allele-independent and are conserved across individuals and between mouse and human, suggesting a conserved function. These regions have an intermediate level of active chromatin marks and their associated genes have intermediate transcriptional activity. Exonic intermediate methylation correlates with exon inclusion at a level between that of fully methylated and unmethylated exons, highlighting gene context-dependent functions. We conclude that intermediate DNA methylation is a conserved signature of gene regulation and exon usage., Many loci in the mammalian genome are intermediately methylated. Here, by comprehensively identifying these loci and quantifying their relationship with gene activity, the authors show that intermediate methylation is an evolutionarily conserved epigenomic signature of gene regulation.

The clonal evolution of tumor cell populations can be reconstructed from patterns of genetic alterations. In contrast, tumor epigenetic states, including DNA methylation, are reversible and sensitive to the tumor microenvironment, presumably precluding the use of epigenetics to discover tumor phylogeny. Here we examined the spatial and temporal dynamics of DNA methylation in a clinically and genetically characterized cohort of IDH1-mutant low-grade gliomas and their patient-matched recurrences. WHO grade II gliomas are diffuse, infiltrative tumors that frequently recur and may undergo malignant progression to a higher grade with a worse prognosis. The extent to which epigenetic alterations contribute to the evolution of low-grade gliomas, including malignant progression, is unknown. While all gliomas in the cohort exhibited the hypermethylation signature associated with IDH1 mutation, low-grade gliomas that underwent malignant progression to high-grade glioblastoma (GBM) had a unique signature of DNA hypomethylation enriched for active enhancers, as well as sites of age-related hypermethylation in the brain. Genes with promoter hypomethylation and concordant transcriptional upregulation during evolution to GBM were enriched in cell cycle function, evolving in concert with genetic alterations that deregulate the G1/S cell cycle checkpoint. Despite the plasticity of tumor epigenetic states, phyloepigenetic trees robustly recapitulated phylogenetic trees derived from somatic mutations in the same patients. These findings highlight widespread co-dependency of genetic and epigenetic events throughout the clonal evolution of initial and recurrent glioma.

CpG islands are present in one-half of all human and mouse genes and typically overlap with promoters or exons. We developed a method for high-resolution analysis of the methylation status of CpG islands genome-wide, using arrays of BAC clones and the methylation-sensitive restriction enzyme NotI. Here we demonstrate the accuracy and specificity of the method. By computationally mapping all NotI sites, methylation events can be defined with single-nucleotide precision throughout the genome. We also demonstrate the unique expandability of the array method using a different methylation-sensitive restriction enzyme, BssHII. We identified and validated new CpG island loci that are methylated in a tissue-specific manner in normal human tissues. The methylation status of the CpG islands is associated with gene expression for several genes, including SHANK3, which encodes a structural protein in neuronal postsynaptic densities. Defects in SHANK3 seem to underlie human 22q13 deletion syndrome. Furthermore, these patterns for SHANK3 are conserved in mice and rats.

Studies have linked the consumption of broccoli and other cruciferous vegetables to a reduced risk of breast cancer. The phytochemical indole-3-carbinol (I3C), present in cruciferous vegetables, and its major acid-catalyzed reaction product 3,3'-diindolylmethane (DIM) have bioactivities relevant to the inhibition of carcinogenesis. In this study, the effect of DIM on angiogenesis and tumorigenesis in a rodent model was investigated. We found that DIM produced a concentration-dependent decrease in proliferation, migration, invasion and capillary tube formation of cultured human umbilical vein endothelial cells (HUVECs). Consistent with its antiproliferative effect, which was significant at only 5 microM DIM, this indole caused a G1 cell cycle arrest in actively proliferating HUVECs. Furthermore, DIM downregulated the expression of cyclin-dependent kinases 2 and 6 (CDK2, CDK6), and upregulated the expression of CDK inhibitor, p27(Kip1), in HUVECs. We observed further in a complementary in vivo Matrigel plug angiogenesis assay that, compared with vehicle control, neovascularization was inhibited up to 76% following the administration of 5 mg/kg DIM to female C57BL/6 mice. Finally, this dose of DIM also inhibited the growth of human MCF-7 cell tumor xenografts by up to 64% in female athymic (nu/nu) mice, compared with the vehicle control. This is the first study to show that DIM can strongly inhibit the development of human breast tumor in a xenograft model and to provide evidence for the antiangiogenic properties of this dietary indole.