Search

BackgroundDespite recent advances in the biology of IDH-wildtype glioblastoma, it remains a devastating disease with median survival of less than 2 years. However, the molecular underpinnings of the heterogeneous response to the current standard-of-care treatment regimen consisting of maximal safe resection, adjuvant radiation, and chemotherapy with temozolomide remain unknown.MethodsComprehensive histopathologic, genomic, and epigenomic evaluation of paired initial and recurrent glioblastoma specimens from 106 patients was performed to investigate the molecular evolution and cellular phenotypes underlying differential treatment responses.ResultsWhile TERT promoter mutation and CDKN2A homozygous deletion were early events during gliomagenesis shared by initial and recurrent tumors, most other recurrent genetic alterations (eg, EGFR, PTEN, and NF1) were commonly private to initial or recurrent tumors indicating acquisition later during clonal evolution. Furthermore, glioblastomas exhibited heterogeneous epigenomic evolution with subsets becoming more globally hypermethylated, hypomethylated, or remaining stable. Glioblastoma that underwent sarcomatous transformation had shorter interval to recurrence and were significantly enriched in NF1, TP53, and RB1 alterations and the mesenchymal epigenetic class. Patients who developed somatic hypermutation following temozolomide treatment had significantly longer interval to disease recurrence and prolonged overall survival, and increased methylation at 4 specific CpG sites in the promoter region of MGMT was significantly associated with this development of hypermutation. Finally, an epigenomic evolution signature incorporating change in DNA methylation levels across 347 critical CpG sites was developed that significantly correlated with clinical outcomes.ConclusionsGlioblastoma undergoes heterogeneous genetic, epigenetic, and cellular evolution that underlies prognostically different treatment responses.

Glioblastoma is a clinically and molecularly heterogeneous disease, and new predictive biomarkers are needed to identify those patients most likely to respond to specific treatments. Through prospective genomic profiling of 459 consecutive primary treatment-naïve IDH-wildtype glioblastomas in adults, we identified a unique subgroup (2%, 9/459) defined by somatic hypermutation and DNA replication repair deficiency due to biallelic inactivation of a canonical mismatch repair gene. The deleterious mutations in mismatch repair genes were often present in the germline in the heterozygous state with somatic inactivation of the remaining allele, consistent with glioblastomas arising due to underlying Lynch syndrome. A subset of tumors had accompanying proofreading domain mutations in the DNA polymerase POLE and resultant "ultrahypermutation". The median age at diagnosis was 50 years (range 27-78), compared with 63 years for the other 450 patients with conventional glioblastoma (p

BackgroundThe TERT promoter mutation (TPM) is acquired in most IDH-wildtype glioblastomas (GBM) and IDH-mutant oligodendrogliomas (OD) enabling tumor cell immortality. Previous studies on TPM clonality show conflicting results. This study was performed to determine whether TPM is clonal on a tumor-wide scale.MethodsWe investigated TPM clonality in relation to presumed early events in 19 IDH-wildtype GBM and 10 IDH-mutant OD using 3-dimensional comprehensive tumor sampling. We performed Sanger sequencing on 264 tumor samples and deep amplicon sequencing on 187 tumor samples. We obtained tumor purity and copy number estimates from whole exome sequencing. TERT expression was assessed by RNA-seq and RNAscope.ResultsWe detected TPM in 100% of tumor samples with quantifiable tumor purity (219 samples). Variant allele frequencies (VAF) of TPM correlate positively with chromosome 10 loss in GBM (R = 0.85), IDH1 mutation in OD (R = 0.87), and with tumor purity (R = 0.91 for GBM; R = 0.90 for OD). In comparison, oncogene amplification was tumor-wide for MDM4- and most EGFR-amplified cases but heterogeneous for MYCN and PDGFRA, and strikingly high in low-purity samples. TPM VAF was moderately correlated with TERT expression (R = 0.52 for GBM; R = 0.65 for OD). TERT expression was detected in a subset of cells, solely in TPM-positive samples, including samples equivocal for tumor.ConclusionsOn a tumor-wide scale, TPM is among the earliest events in glioma evolution. Intercellular heterogeneity of TERT expression, however, suggests dynamic regulation during tumor growth. TERT expression may be a tumor cell-specific biomarker.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Tumors may contain billions of cells, including distinct malignant clones and nonmalignant cell types. Clarifying the evolutionary histories, prevalence, and defining molecular features of these cells is essential for improving clinical outcomes, since intratumoral heterogeneity provides fuel for acquired resistance to targeted therapies. Here we present a statistically motivated strategy for deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial tumor sections (MOMA). By combining deep sampling of IDH-mutant astrocytomas with integrative analysis of single-nucleotide variants, copy-number variants, and gene expression, we reconstruct and validate the phylogenies, spatial distributions, and transcriptional profiles of distinct malignant clones. By genotyping nuclei analyzed by single-nucleus RNA-seq for truncal mutations, we further show that commonly used algorithms for identifying cancer cells from single-cell transcriptomes may be inaccurate. We also demonstrate that correlating gene expression with tumor purity in bulk samples can reveal optimal markers of malignant cells and use this approach to identify a core set of genes that are consistently expressed by astrocytoma truncal clones, including AKR1C3, whose expression is associated with poor outcomes in several types of cancer. In summary, MOMA provides a robust and flexible strategy for precisely deconstructing intratumoral heterogeneity and clarifying the core molecular properties of distinct cellular populations in solid tumors.

Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter- and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of antigens. In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface antigens that could be targeted by antibodies and chimeric antigen receptor-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas) and 9166 normal tissue samples (from the Genotype-Tissue Expression project), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN, which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative antigens. Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

BackgroundSchwannomas are common peripheral nerve sheath tumors that can cause severe morbidity given their stereotypic intracranial and paraspinal locations. Similar to many solid tumors, schwannomas and other nerve sheath tumors are primarily thought to arise due to aberrant hyperactivation of the RAS growth factor signaling pathway. Here, we sought to further define the molecular pathogenesis of schwannomas.MethodsWe performed comprehensive genomic profiling on a cohort of 96 human schwannomas, as well as DNA methylation profiling on a subset. Functional studies including RNA sequencing, chromatin immunoprecipitation-DNA sequencing, electrophoretic mobility shift assay, and luciferase reporter assays were performed in a fetal glial cell model following transduction with wildtype and tumor-derived mutant isoforms of SOX10.ResultsWe identified that nearly one-third of sporadic schwannomas lack alterations in known nerve sheath tumor genes and instead harbor novel recurrent in-frame insertion/deletion mutations in SOX10, which encodes a transcription factor responsible for controlling Schwann cell differentiation and myelination. SOX10 indel mutations were highly enriched in schwannomas arising from nonvestibular cranial nerves (eg facial, trigeminal, vagus) and were absent from vestibular nerve schwannomas driven by NF2 mutation. Functional studies revealed these SOX10 indel mutations have retained DNA binding capacity but impaired transactivation of glial differentiation and myelination gene programs.ConclusionsWe thus speculate that SOX10 indel mutations drive a unique subtype of schwannomas by impeding proper differentiation of immature Schwann cells.

Cancer-associated fibroblasts (CAFs) were presumed absent in glioblastoma given the lack of brain fibroblasts. Serial trypsinization of glioblastoma specimens yielded cells with CAF morphology and single-cell transcriptomic profiles based on their lack of copy number variations (CNVs) and elevated individual cell CAF probability scores derived from the expression of 9 CAF markers and absence of 5 markers from non-CAF stromal cells sharing features with CAFs. Cells without CNVs and with high CAF probability scores were identified in single-cell RNA-Seq of 12 patient glioblastomas. Pseudotime reconstruction revealed that immature CAFs evolved into subtypes, with mature CAFs expressing actin alpha 2, smooth muscle (ACTA2). Spatial transcriptomics from 16 patient glioblastomas confirmed CAF proximity to mesenchymal glioblastoma stem cells (GSCs), endothelial cells, and M2 macrophages. CAFs were chemotactically attracted to GSCs, and CAFs enriched GSCs. We created a resource of inferred crosstalk by mapping expression of receptors to their cognate ligands, identifying PDGF and TGF-β as mediators of GSC effects on CAFs and osteopontin and HGF as mediators of CAF-induced GSC enrichment. CAFs induced M2 macrophage polarization by producing the extra domain A (EDA) fibronectin variant that binds macrophage TLR4. Supplementing GSC-derived xenografts with CAFs enhanced in vivo tumor growth. These findings are among the first to identify glioblastoma CAFs and their GSC interactions, making them an intriguing target.

BackgroundIDH-mutant diffuse gliomas are heterogeneous, and improved methods for optimal patient therapeutic stratification are needed. PI3K/AKT/mTOR signaling activity can drive disease progression and potential therapeutic inhibitors of the pathway are available. Yet, the prevalence of PI3K/AKT/mTOR signaling pathway activity in IDH-mutant glioma is unclear and few robust strategies to assess activity in clinical samples exist.MethodsPI3K/AKT/mTOR signaling pathway activity was evaluated in a retrospective cohort of 132 IDH-mutant diffuse glioma (91 astrocytoma and 41 oligodendroglioma, 1p/19q-codeleted) through quantitative multiplex immunoprofiling using phospho-specific antibodies for PI3K/AKT/mTOR pathway members, PRAS40, RPS6, and 4EBP1, and tumor-specific anti-IDH1 R132H. Expression levels were correlated with genomic evaluation of pathway intrinsic genes and univariate and multivariate Cox proportional hazard regression models were used to evaluate the relationship with outcome.ResultsTumor-specific expression of p-PRAS40, p-RPS6, and p-4EBP1 was common in IDH-mutant diffuse glioma and increased with CNS WHO grade from 2 to 3. Genomic analysis predicted pathway activity in 21.7% (13/60) while protein evaluation identified active PI3K/AKT/mTOR signaling in 56.6% (34/60). Comparison of expression in male versus female patients suggested sexual dimorphism. Of particular interest, when adjusting for clinical prognostic factors, the level of phosphorylation of RPS6 was strongly associated with PFS (P

Meningiomas are the most common primary intracranial tumors. There are no effective medical therapies for meningioma patients, and new treatments have been encumbered by limited understanding of meningioma biology. Here, we use DNA methylation profiling on 565 meningiomas integrated with genetic, transcriptomic, biochemical, proteomic and single-cell approaches to show meningiomas are composed of three DNA methylation groups with distinct clinical outcomes, biological drivers and therapeutic vulnerabilities. Merlin-intact meningiomas (34%) have the best outcomes and are distinguished by NF2/Merlin regulation of susceptibility to cytotoxic therapy. Immune-enriched meningiomas (38%) have intermediate outcomes and are distinguished by immune infiltration, HLA expression and lymphatic vessels. Hypermitotic meningiomas (28%) have the worst outcomes and are distinguished by convergent genetic and epigenetic mechanisms driving the cell cycle and resistance to cytotoxic therapy. To translate these findings into clinical practice, we show cytostatic cell cycle inhibitors attenuate meningioma growth in cell culture, organoids, xenografts and patients.

BackgroundIntratumoral heterogeneity is a hallmark of diffuse gliomas. DNA methylation profiling is an emerging approach in the clinical classification of brain tumors. The goal of this study is to investigate the effects of intratumoral heterogeneity on classification confidence.MethodsWe used neuronavigation to acquire 133 image-guided and spatially separated stereotactic biopsy samples from 16 adult patients with a diffuse glioma (7 IDH-wildtype and 2 IDH-mutant glioblastoma, 6 diffuse astrocytoma, IDH-mutant and 1 oligodendroglioma, IDH-mutant and 1p19q codeleted), which we characterized using DNA methylation arrays. Samples were obtained from regions with and without abnormalities on contrast-enhanced T1-weighted and fluid-attenuated inversion recovery MRI. Methylation profiles were analyzed to devise a 3-dimensional reconstruction of (epi)genetic heterogeneity. Tumor purity was assessed from clonal methylation sites.ResultsMolecular aberrations indicated that tumor was found outside imaging abnormalities, underlining the infiltrative nature of this tumor and the limitations of current routine imaging modalities. We demonstrate that tumor purity is highly variable between samples and explains a substantial part of apparent epigenetic spatial heterogeneity. We observed that DNA methylation subtypes are often, but not always, conserved in space taking tumor purity and prediction accuracy into account.ConclusionOur results underscore the infiltrative nature of diffuse gliomas and suggest that DNA methylation subtypes are relatively concordant in this tumor type, although some heterogeneity exists.

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.

BackgroundLarge-scale genome-wide association studies (GWAS) have implicated thousands of germline genetic variants in modulating individuals' risk to various diseases, including cancer. At least 25 risk loci have been identified for low-grade gliomas (LGGs), but their molecular functions remain largely unknown.MethodsWe hypothesized that GWAS loci contain causal single nucleotide polymorphisms (SNPs) that reside in accessible open chromatin regions and modulate the expression of target genes by perturbing the binding affinity of transcription factors (TFs). We performed an integrative analysis of genomic and epigenomic data from The Cancer Genome Atlas and other public repositories to identify candidate causal SNPs within linkage disequilibrium blocks of LGG GWAS loci. We assessed their potential regulatory role via in silico TF binding sequence perturbations, convolutional neural network trained on TF binding data, and simulated annealing-based interpretation methods.ResultsWe built an interactive website (http://education.knoweng.org/alg3/) summarizing the functional footprinting of 280 variants in 25 LGG GWAS regions, providing rich information for further computational and experimental scrutiny. We identified as case studies PHLDB1 and SLC25A26 as candidate target genes of rs12803321 and rs11706832, respectively, and predicted the GWAS variant rs648044 to be the causal SNP modulating ZBTB16, a known tumor suppressor in multiple cancers. We showed that rs648044 likely perturbed the binding affinity of the TF MAFF, as supported by RNA interference and in vitro MAFF binding experiments.ConclusionsThe identified candidate (causal SNP, target gene, TF) triplets and the accompanying resource will help accelerate our understanding of the molecular mechanisms underlying genetic risk factors for gliomas.

Among the large, diverse set of mammalian long noncoding RNAs (lncRNAs), long noncoding primary microRNAs (lnc-pri-miRNAs) are those that host miRNAs. Whether lnc-pri-miRNA loci have important biological function independent of their cognate miRNAs is poorly understood. From a genome-scale lncRNA screen, lnc-pri-miRNA loci were enriched for function in cell proliferation, and in glioblastoma (i.e., GBM) cells with DGCR8 or DROSHA knockdown, lnc-pri-miRNA screen hits still regulated cell growth. To molecularly dissect the function of a lnc-pri-miRNA locus, we studied LOC646329 (also known as MIR29HG), which hosts the miR-29a/b1 cluster. In GBM cells, LOC646329 knockdown reduced miR-29a/b1 levels, and these cells exhibited decreased growth. However, genetic deletion of the miR-29a/b1 cluster (LOC646329-miR29Δ) did not decrease cell growth, while knockdown of LOC646329-miR29Δ transcripts reduced cell proliferation. The miR-29a/b1-independent activity of LOC646329 corresponded to enhancer-like activation of a neighboring oncogene (MKLN1), regulating cell propagation. The LOC646329 locus interacts with the MKLN1 promoter, and antisense oligonucleotide knockdown of the lncRNA disrupts these interactions and reduces the enhancer-like activity. More broadly, analysis of genome-wide data from multiple human cell types showed that lnc-pri-miRNA loci are significantly enriched for DNA looping interactions with gene promoters as well as genomic and epigenetic characteristics of transcriptional enhancers. Functional studies of additional lnc-pri-miRNA loci demonstrated cognate miRNA-independent enhancer-like activity. Together, these data demonstrate that lnc-pri-miRNA loci can regulate cell biology via both miRNA-dependent and miRNA-independent mechanisms.

BackgroundEmerging 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.MethodsWe 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.ResultsMethylation 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.ConclusionsThese 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.

Meningiomas are the most common primary intracranial tumors, but the molecular drivers of meningioma tumorigenesis are poorly understood. We hypothesized that investigating intratumor heterogeneity in meningiomas would elucidate biologic drivers and reveal new targets for molecular therapy. To test this hypothesis, here we perform multiplatform molecular profiling of 86 spatially-distinct samples from 13 human meningiomas. Our data reveal that regional alterations in chromosome structure underlie clonal transcriptomic, epigenomic, and histopathologic signatures in meningioma. Stereotactic co-registration of sample coordinates to preoperative magnetic resonance images further suggest that high apparent diffusion coefficient (ADC) distinguishes meningioma regions with proliferating cells enriched for developmental gene expression programs. To understand the function of these genes in meningioma, we develop a human cerebral organoid model of meningioma and validate the high ADC marker genes CDH2 and PTPRZ1 as potential targets for meningioma therapy using live imaging, single cell RNA sequencing, CRISPR interference, and pharmacology.

Although DNA methylation is a key regulator of gene expression, the comprehensive methylation landscape of metastatic cancer has never been defined. Through whole-genome bisulfite sequencing paired with deep whole-genome and transcriptome sequencing of 100 castration-resistant prostate metastases, we discovered alterations affecting driver genes that were detectable only with integrated whole-genome approaches. Notably, we observed that 22% of tumors exhibited a novel epigenomic subtype associated with hypermethylation and somatic mutations in TET2, DNMT3B, IDH1 and BRAF. We also identified intergenic regions where methylation is associated with RNA expression of the oncogenic driver genes AR, MYC and ERG. Finally, we showed that differential methylation during progression preferentially occurs at somatic mutational hotspots and putative regulatory regions. This study is a large integrated study of whole-genome, whole-methylome and whole-transcriptome sequencing in metastatic cancer that provides a comprehensive overview of the important regulatory role of methylation in metastatic castration-resistant prostate cancer.

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.

Ependymomas exist within distinct genetic subgroups, but the molecular diversity within individual ependymomas is unknown. We perform multiplatform molecular profiling of 6 spatially distinct samples from an ependymoma with C11orf95-RELA fusion. DNA methylation and RNA sequencing distinguish clusters of samples according to neuronal development gene expression programs that could also be delineated by differences in magnetic resonance blood perfusion. Exome sequencing and phylogenetic analysis reveal epigenomic intratumor heterogeneity and suggest that chromosomal structural alterations may precede accumulation of single-nucleotide variants during ependymoma tumorigenesis. In sum, these findings shed light on the oncogenesis and intratumor heterogeneity of ependymoma.

BackgroundIDH-mutant lower-grade gliomas (LGGs) evolve under the selective pressure of therapy, but well-characterized patient-derived cells (PDCs) modeling evolutionary stages are lacking. IDH-mutant LGGs may develop therapeutic resistance associated with chemotherapy-driven hypermutation and malignant progression. The aim of this study was to establish and characterize PDCs, single-cell-derived PDCs (scPDCs), and xenografts (PDX) of IDH1-mutant recurrences representing distinct stages of tumor evolution.MethodsWe derived and validated cell cultures from IDH1-mutant recurrences of astrocytoma and oligodendroglioma. We used exome sequencing and phylogenetic reconstruction to examine the evolutionary stage represented by PDCs, scPDCs, and PDX relative to corresponding spatiotemporal tumor tissue and germline DNA. PDCs were also characterized for growth and tumor immortality phenotypes, and PDX were examined histologically.ResultsThe integrated astrocytoma phylogeny revealed 2 independent founder clonal expansions of hypermutated (HM) cells in tumor tissue that are faithfully represented by independent PDCs. The oligodendroglioma phylogeny showed more than 4000 temozolomide-associated mutations shared among tumor samples, PDCs, scPDCs, and PDX, suggesting a shared monoclonal origin. The PDCs from both subtypes exhibited hallmarks of tumorigenesis, retention of subtype-defining genomic features, production of 2-hydroxyglutarate, and subtype-specific telomere maintenance mechanisms that confer tumor cell immortality. The oligodendroglioma PDCs formed infiltrative intracranial tumors with characteristic histology.ConclusionsThese PDCs, scPDCs, and PDX are unique and versatile community resources that model the heterogeneous clonal origins and functions of recurrent IDH1-mutant LGGs. The integrated phylogenies advance our knowledge of the complex evolution and immense mutational load of IDH1-mutant HM glioma.

The majority of common variants associated with common diseases, as well as an unknown proportion of causal mutations for rare diseases, fall in noncoding regions of the genome. Although catalogs of noncoding regulatory elements are steadily improving, we have a limited understanding of the functional effects of mutations within them. Here, we perform saturation mutagenesis in conjunction with massively parallel reporter assays on 20 disease-associated gene promoters and enhancers, generating functional measurements for over 30,000 single nucleotide substitutions and deletions. We find that the density of putative transcription factor binding sites varies widely between regulatory elements, as does the extent to which evolutionary conservation or integrative scores predict functional effects. These data provide a powerful resource for interpreting the pathogenicity of clinically observed mutations in these disease-associated regulatory elements, and comprise a rich dataset for the further development of algorithms that aim to predict the regulatory effects of noncoding mutations.

Adult diffuse gliomas are a diverse group of brain neoplasms that inflict a high emotional toll on patients and their families. The Cancer Genome Atlas and similar projects have provided a comprehensive understanding of the somatic alterations and molecular subtypes of glioma at diagnosis. However, gliomas undergo significant cellular and molecular evolution during disease progression. We review the current knowledge on the genomic and epigenetic abnormalities in primary tumors and after disease recurrence, highlight the gaps in the literature, and elaborate on the need for a new multi-institutional effort to bridge these knowledge gaps and how the Glioma Longitudinal Analysis Consortium (GLASS) aims to systemically catalog the longitudinal changes in gliomas. The GLASS initiative will provide essential insights into the evolution of glioma toward a lethal phenotype, with the potential to reveal targetable vulnerabilities and, ultimately, improved outcomes for a patient population in need.

BackgroundRare 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.MethodsIn 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.ResultsOne 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.ConclusionsEach 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.

The lncRNA HOTAIR has been implicated in several human cancers. Here, we evaluated the molecular alterations and upstream regulatory mechanisms of HOTAIR in glioma, the most common primary brain tumors, and its clinical relevance. HOTAIR gene expression, methylation, copy-number and prognostic value were investigated in human gliomas integrating data from online datasets and our cohorts. High levels of HOTAIR were associated with higher grades of glioma, particularly IDH wild-type cases. Mechanistically, HOTAIR was overexpressed in a gene dosage-independent manner, while DNA methylation levels of particular CpGs in HOTAIR locus were associated with HOTAIR expression levels in GBM clinical specimens and cell lines. Concordantly, the demethylating agent 5-Aza-2'-deoxycytidine affected HOTAIR transcriptional levels in a cell line-dependent manner. Importantly, HOTAIR was frequently co-expressed with HOXA9 in high-grade gliomas from TCGA, Oncomine, and our Portuguese and French datasets. Integrated in silico analyses, chromatin immunoprecipitation, and qPCR data showed that HOXA9 binds directly to the promoter of HOTAIR. Clinically, GBM patients with high HOTAIR expression had a significantly reduced overall survival, independently of other prognostic variables. In summary, this work reveals HOXA9 as a novel direct regulator of HOTAIR, and establishes HOTAIR as an independent prognostic marker, providing new therapeutic opportunities to treat this highly aggressive cancer.

Approximately 8% of the human genome is comprised of endogenous retroviral insertions (ERVs) originating from historic retroviral integration into germ cells. The function of ERVs as regulators of gene expression is well established. Less well studied are insertional polymorphisms of ERVs and their contribution to the heritability of complex phenotypes. The most recent integration of ERV, HERV-K, is expressed in a range of complex human conditions from cancer to neurologic diseases. Using an in-house computational pipeline and whole-genome sequencing data from the diverse 1,000 Genomes Phase 3 population (n = 2,504), we identified 46 polymorphic HERV-K insertions that are tagged by adjacent single nucleotide polymorphisms (SNPs). To test the potential role of polymorphic HERV-K in the heritability of complex diseases, existing databases were queried for enrichment of established relationships between the HERV-K insertion-associated SNPs (hiSNPs), and tissue specific gene expression and disease phenotypes. Overall, hiSNPs for the 46 polymorphic HERV-K sites were statistically enriched (p < 1.0E-16) for eQTLs across 44 human tissues. Fifteen of the 46 HERV-K insertions had hiSNPs annotated in the EMBL-EBI GWAS Catalog and cumulatively associated with >100 phenotypes. Experimental factor ontology enrichment analysis suggests that polymorphic HERV-K specifically contribute to neurologic and immunologic disease phenotypes, including traits related to intra cranial volume (FDR 2.00E-09), Parkinson's disease (FDR 1.80E-09), and autoimmune diseases (FDR 1.80E-09). These results provide strong candidates for context-specific study of polymorphic HERV-K insertions in disease-related traits, serving as a roadmap for future studies of the heritability of complex disease.

We present an extensive assessment of mutation burden through sequencing analysis of >81,000 tumors from pediatric and adult patients, including tumors with hypermutation caused by chemotherapy, carcinogens, or germline alterations. Hypermutation was detected in tumor types not previously associated with high mutation burden. Replication repair deficiency was a major contributing factor. We uncovered new driver mutations in the replication-repair-associated DNA polymerases and a distinct impact of microsatellite instability and replication repair deficiency on the scale of mutation load. Unbiased clustering, based on mutational context, revealed clinically relevant subgroups regardless of the tumors' tissue of origin, highlighting similarities in evolutionary dynamics leading to hypermutation. Mutagens, such as UV light, were implicated in unexpected cancers, including sarcomas and lung tumors. The order of mutational signatures identified previous treatment and germline replication repair deficiency, which improved management of patients and families. These data will inform tumor classification, genetic testing, and clinical trial design.

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.

Despite the great progress in our understanding of the molecular basis of human cancer, the heterogeneity of individual tumors and the evolutionary pressures imposed by therapy have hampered our ability to effectively eradicate and control this disease. How, therefore, do cancers evolve under the selective pressures of cancer therapy? Recent studies have linked both primary (or de novo) and acquired treatment resistance to intratumor heterogeneity and clonal evolution. Resistance to targeted therapies often includes mutation of the drug target itself and aberrations of pathways upstream of, downstream from, or parallel to the drug target. For systemic chemotherapies, discrete and recurrent resistance-conferring genetic aberrations have eluded the community, due in part to their wide-ranging mutagenic effects. In this review, we discuss different patterns of clonal evolution during treatment-specific selective pressures and focus on the genetic mechanisms of treatment resistance that have emerged to both targeted therapies and chemotherapies.

Mutations in the isocitrate dehydrogenase genes IDH1 and IDH2 are among the first genetic alterations observed during the development of lower-grade glioma (LGG). LGG-associated IDH mutations confer gain-of-function activity by converting α-ketoglutarate to the oncometabolite R-2-hydroxyglutarate (2HG). Clinical samples and gene expression data from The Cancer Genome Atlas (TCGA) demonstrate reduced expression of cytotoxic T lymphocyte-associated genes and IFN-γ-inducible chemokines, including CXCL10, in IDH-mutated (IDH-MUT) tumors compared with IDH-WT tumors. Given these findings, we have investigated the impact of IDH mutations on the immunological milieu in LGG. In immortalized normal human astrocytes (NHAs) and syngeneic mouse glioma models, the introduction of mutant IDH1 or treatment with 2HG reduced levels of CXCL10, which was associated with decreased production of STAT1, a regulator of CXCL10. Expression of mutant IDH1 also suppressed the accumulation of T cells in tumor sites. Reductions in CXCL10 and T cell accumulation were reversed by IDH-C35, a specific inhibitor of mutant IDH1. Furthermore, IDH-C35 enhanced the efficacy of vaccine immunotherapy in mice bearing IDH-MUT gliomas. Our findings demonstrate a mechanism of immune evasion in IDH-MUT gliomas and suggest that specific inhibitors of mutant IDH may improve the efficacy of immunotherapy in patients with IDH-MUT gliomas.

Purpose: Second malignant neoplasms (SMNs) are severe late complications that occur in pediatric cancer survivors exposed to radiotherapy and other genotoxic treatments. To characterize the mutational landscape of treatment-induced sarcomas and to identify candidate SMN-predisposing variants, we analyzed germline and SMN samples from pediatric cancer survivors.Experimental Design: We performed whole-exome sequencing (WES) and RNA sequencing on radiation-induced sarcomas arising from two pediatric cancer survivors. To assess the frequency of germline TP53 variants in SMNs, Sanger sequencing was performed to analyze germline TP53 in 37 pediatric cancer survivors from the Childhood Cancer Survivor Study (CCSS) without any history of a familial cancer predisposition syndrome but known to have developed SMNs.Results: WES revealed TP53 mutations involving p53's DNA-binding domain in both index cases, one of which was also present in the germline. The germline and somatic TP53-mutant variants were enriched in the transcriptomes for both sarcomas. Analysis of TP53-coding exons in germline specimens from the CCSS survivor cohort identified a G215C variant encoding an R72P amino acid substitution in 6 patients and a synonymous SNP A639G in 4 others, resulting in 10 of 37 evaluable patients (27%) harboring a germline TP53 variant.Conclusions: Currently, germline TP53 is not routinely assessed in patients with pediatric cancer. These data support the concept that identifying germline TP53 variants at the time a primary cancer is diagnosed may identify patients at high risk for SMN development, who could benefit from modified therapeutic strategies and/or intensive posttreatment monitoring. Clin Cancer Res; 23(7); 1852-61. ©2016 AACR.

Juvenile myelomonocytic leukemia (JMML) is a myeloproliferative disorder of childhood caused by mutations in the Ras pathway. Outcomes in JMML vary markedly from spontaneous resolution to rapid relapse after hematopoietic stem cell transplantation. Here, we hypothesized that DNA methylation patterns would help predict disease outcome and therefore performed genome-wide DNA methylation profiling in a cohort of 39 patients. Unsupervised hierarchical clustering identifies three clusters of patients. Importantly, these clusters differ significantly in terms of 4-year event-free survival, with the lowest methylation cluster having the highest rates of survival. These findings were validated in an independent cohort of 40 patients. Notably, all but one of 14 patients experiencing spontaneous resolution cluster together and closer to 22 healthy controls than to other JMML cases. Thus, we show that DNA methylation patterns in JMML are predictive of outcome and can identify the patients most likely to experience spontaneous resolution.

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.

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.

UnlabelledRegulation of mRNA splicing, a critical and tightly regulated cellular function, underlies the majority of proteomic diversity and is frequently disrupted in disease. Using an integrative genomics approach, we combined both genomic data and exon-level transcriptome data in two somatic tissues (cerebella and peripheral ganglia) from a transgenic mouse model of neuroblastoma, a tumor that arises from the peripheral neural crest. Here, we describe splicing quantitative trait loci associated with differential splicing across the genome that we use to identify genes with previously unknown functions within the splicing pathway and to define de novo intronic splicing motifs that influence splicing from hundreds of bases away. Our results show that these splicing motifs represent sites for functional recurrent mutations and highlight novel candidate genes in human cancers, including childhood neuroblastoma.SignificanceSomatic mutations with predictable downstream effects are largely relegated to coding regions, which comprise less than 2% of the human genome. Using an unbiased in vivo analysis of a mouse model of neuroblastoma, we have identified intronic splicing motifs that translate into sites for recurrent somatic mutations in human cancers.

Somatic mutations in cancer are more frequent in heterochromatic and late-replicating regions of the genome. We report that regional disparities in mutation density are virtually abolished within transcriptionally silent genomic regions of cutaneous squamous cell carcinomas (cSCCs) arising in an XPC(-/-) background. XPC(-/-) cells lack global genome nucleotide excision repair (GG-NER), thus establishing differential access of DNA repair machinery within chromatin-rich regions of the genome as the primary cause for the regional disparity. Strikingly, we find that increasing levels of transcription reduce mutation prevalence on both strands of gene bodies embedded within H3K9me3-dense regions, and only to those levels observed in H3K9me3-sparse regions, also in an XPC-dependent manner. Therefore, transcription appears to reduce mutation prevalence specifically by relieving the constraints imposed by chromatin structure on DNA repair. We model this relationship among transcription, chromatin state, and DNA repair, revealing a new, personalized determinant of cancer risk.

SummaryWe present methylC track, an efficient mechanism for visualizing single-base resolution DNA methylation data on a genome browser. The methylC track dynamically integrates the level of methylation, the position and context of the methylated cytosine (i.e. CG, CHG and CHH), strand and confidence level (e.g. read coverage depth in the case of whole-genome bisulfite sequencing data). Investigators can access and integrate these information visually at specific locus or at the genome-wide level on the WashU EpiGenome Browser in the context of other rich epigenomic datasets.Availability and implementationThe methylC track is part of the WashU EpiGenome Browser, which is open source and freely available at http://epigenomegateway.wustl.edu/browser/. The most up-to-date instructions and tools for preparing methylC track are available at http://epigenomegateway.wustl.edu/+/cmtk.Contacttwang@genetics.wustl.eduSupplementary informationSupplementary data are available at Bioinformatics online.

Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

We have extended our understanding of the molecular biology that underlies adult glioblastoma over many years. By contrast, high-grade gliomas in children and adolescents have remained a relatively under-investigated disease. The latest large-scale genomic and epigenomic profiling studies have yielded an unprecedented abundance of novel data and provided deeper insights into gliomagenesis across all age groups, which has highlighted key distinctions but also some commonalities. As we are on the verge of dissecting glioblastomas into meaningful biological subgroups, this Review summarizes the hallmark genetic alterations that are associated with distinct epigenetic features and patient characteristics in both paediatric and adult disease, and examines the complex interplay between the glioblastoma genome and epigenome.

Tumor recurrence is a leading cause of cancer mortality. Therapies for recurrent disease may fail, at least in part, because the genomic alterations driving the growth of recurrences are distinct from those in the initial tumor. To explore this hypothesis, we sequenced the exomes of 23 initial low-grade gliomas and recurrent tumors resected from the same patients. In 43% of cases, at least half of the mutations in the initial tumor were undetected at recurrence, including driver mutations in TP53, ATRX, SMARCA4, and BRAF; this suggests that recurrent tumors are often seeded by cells derived from the initial tumor at a very early stage of their evolution. Notably, tumors from 6 of 10 patients treated with the chemotherapeutic drug temozolomide (TMZ) followed an alternative evolutionary path to high-grade glioma. At recurrence, these tumors were hypermutated and harbored driver mutations in the RB (retinoblastoma) and Akt-mTOR (mammalian target of rapamycin) pathways that bore the signature of TMZ-induced mutagenesis.

Diffuse gliomas consist of both low- and high-grade varieties, each with distinct morphological and biological features. The often extended periods of relative indolence exhibited by low-grade gliomas (LGG; WHO grade II) differ sharply from the aggressive, rapidly fatal clinical course of primary glioblastoma (GBM; WHO grade IV). Nevertheless, until recently, the molecular foundations underlying this stark biological contrast between glioma variants remained largely unknown. The discoveries of distinctive and highly recurrent genomic and epigenomic abnormalities in LGG have both informed a more accurate classification scheme and pointed to viable avenues for therapeutic development. As such, the field of neuro-oncology now seems poised to capitalize on these gains to achieve significant benefit for LGG patients. This report will briefly recount the proceedings of a workshop held in January 2013 and hosted by Accelerate Brain Cancer Cure (ABC(2)) on the subject of LGG. While much of the meeting covered recent insights into LGG biology, its focus remained on how best to advance the clinical management, whether by improved preclinical modeling, more effective targeted therapeutics and clinical trial design, or innovative imaging technology.

Recent advancements in sequencing-based DNA methylation profiling methods provide an unprecedented opportunity to map complete DNA methylomes. These include whole-genome bisulfite sequencing (WGBS, MethylC-seq, or BS-seq), reduced-representation bisulfite sequencing (RRBS), and enrichment-based methods such as MeDIP-seq, MBD-seq, and MRE-seq. These methods yield largely comparable results but differ significantly in extent of genomic CpG coverage, resolution, quantitative accuracy, and cost, at least while using current algorithms to interrogate the data. None of these existing methods provides single-CpG resolution, comprehensive genome-wide coverage, and cost feasibility for a typical laboratory. We introduce methylCRF, a novel conditional random fields-based algorithm that integrates methylated DNA immunoprecipitation (MeDIP-seq) and methylation-sensitive restriction enzyme (MRE-seq) sequencing data to predict DNA methylation levels at single-CpG resolution. Our method is a combined computational and experimental strategy to produce DNA methylomes of all 28 million CpGs in the human genome for a fraction (

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.

Accurate detection of epimutations in tumor cells is crucial for -understanding the molecular pathogenesis of cancer. Alterations in DNA methylation in cancer are functionally important and clinically relevant, but even this well-studied area is continually re-evaluated in light of unanticipated results, such as the strong association between aberrant DNA methylation in adult tumors and polycomb group profiles in embryonic stem cells, cancer-associated genetic mutations in epigenetic regulators such as DNMT3A and TET family genes, and the discovery of altered 5-hydroxymethylcytosine, a product of TET proteins acting on 5-methylcytosine, in human tumors with TET mutations. The abundance and distribution of covalent histone modifications in primary cancer tissues relative to normal cells is an important but largely uncharted area, although there is good evidence for a mechanistic role of cancer-specific alterations in histone modifications in tumor etiology, drug response, and tumor progression. Meanwhile, the discovery of new epigenetic marks continues, and there are many useful methods for epigenome analysis applicable to primary tumor samples, in addition to cancer cell lines. For DNA methylation and hydroxymethylation, next-generation sequencing allows increasingly inexpensive and quantitative whole-genome profiling. Similarly, the refinement and maturation of chromatin immunoprecipitation with next-generation sequencing (ChIP-seq) has made possible genome-wide mapping of histone modifications, open chromatin, and transcription factor binding sites. Computational tools have been developed apace with these epigenome methods to better enable accurate interpretation of the profiling data.

Biologists possess the detailed knowledge critical for extracting biological insight from genome-wide data resources, and yet they are increasingly faced with nontrivial computational analysis challenges posed by genome-scale methodologies. To lower this computational barrier, particularly in the early data exploration phases, we have developed an interactive pattern discovery and visualization approach, Spark, designed with epigenomic data in mind. Here we demonstrate Spark's ability to reveal both known and novel epigenetic signatures, including a previously unappreciated binding association between the YY1 transcription factor and the corepressor CTBP2 in human embryonic stem cells.

The H3K9me3 histone modification is often found at promoter regions, where it functions to repress transcription. However, we have previously shown that 3' exons of zinc finger genes (ZNFs) are marked by high levels of H3K9me3. We have now further investigated this unusual location for H3K9me3 in ZNF genes. Neither bioinformatic nor experimental approaches support the hypothesis that the 3' exons of ZNFs are promoters. We further characterized the histone modifications at the 3' ZNF exons and found that these regions also contain H3K36me3, a mark of transcriptional elongation. A genome-wide analysis of ChIP-seq data revealed that ZNFs constitute the majority of genes that have high levels of both H3K9me3 and H3K36me3. These results suggested the possibility that the ZNF genes may be imprinted, with one allele transcribed and one allele repressed. To test the hypothesis that the contradictory modifications are due to imprinting, we used a SNP analysis of RNA-seq data to demonstrate that both alleles of certain ZNF genes having H3K9me3 and H3K36me3 are transcribed. We next analyzed isolated ZNF 3' exons using stably integrated episomes. We found that although the H3K36me3 mark was lost when the 3' ZNF exon was removed from its natural genomic location, the isolated ZNF 3' exons retained the H3K9me3 mark. Thus, the H3K9me3 mark at ZNF 3' exons does not impede transcription and it is regulated independently of the H3K36me3 mark. Finally, we demonstrate a strong relationship between the number of tandemly repeated domains in the 3' exons and the H3K9me3 mark. We suggest that the H3K9me3 at ZNF 3' exons may function to protect the genome from inappropriate recombination rather than to regulate transcription.

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.

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.

This article discusses research being done on the human epigenome. It references a study by R. Lister et al., published in the journal "Nature" in 2009. The study provided complete maps of deoxyribonucleic acid (DNA) methylation in two human cell types and integrated the maps with ribonucleic acid (RNA) transcript sequences and expression levels. It identified human embryonic stem (ES) cells as the most interesting cell type for comprehensive methylome mapping because of their inherent developmental potential and pluripotency.

The article offers information on the Roadmap Epigenome Browser which is based on the Human Epigenome Browser from the Washington University and integrates data from the U.S. National Institute of Health Roadmap Epigenomics Consortium and ENCODE list of functional elements in the human genome. It informs that Roadmap Epigenome Brower enables researchers to explore the tissue-specific regulatory roles of genetic variants in the context of disease. It discusses single nucleotide polymorphism.

The authors focus on a research related to the genome interactions in eukaryotes performed using the WashU Epigenome Browser located at http://epigenomegateway.wustl.edu. They state that they have tried to explore the histone modification profile of human cells and long-range interaction data of human genome. They also mention that human genomes can show similar interaction pattern despite they have different histone modifications.