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.

Developmental history shapes the epigenome and biological function of differentiated cells. Epigenomic patterns have been broadly attributed to the three embryonic germ layers. Here we investigate how developmental origin influences epigenomes. We compare key epigenomes of cell types derived from surface ectoderm (SE), including keratinocytes and breast luminal and myoepithelial cells, against neural crest-derived melanocytes and mesoderm-derived dermal fibroblasts, to identify SE differentially methylated regions (SE-DMRs). DNA methylomes of neonatal keratinocytes share many more DMRs with adult breast luminal and myoepithelial cells than with melanocytes and fibroblasts from the same neonatal skin. This suggests that SE origin contributes to DNA methylation patterning, while shared skin tissue environment has limited effect on epidermal keratinocytes. Hypomethylated SE-DMRs are in proximity to genes with SE relevant functions. They are also enriched for enhancer- and promoter-associated histone modifications in SE-derived cells, and for binding motifs of transcription factors important in keratinocyte and mammary gland biology. Thus, epigenomic analysis of cell types with common developmental origin reveals an epigenetic signature that underlies a shared gene regulatory network.

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.

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.

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.

PI5P4Ks are a class of phosphoinositide kinases that phosphorylate PI-5-P to PI-4,5-P2. Distinct localization of phosphoinositides is fundamental for a multitude of cellular functions. Here, we identify a role for peroxisomal PI-4,5-P2 generated by the PI5P4Ks in maintaining energy balance. We demonstrate that PI-4,5-P2 regulates peroxisomal fatty acid oxidation by mediating trafficking of lipid droplets to peroxisomes, which is essential for sustaining mitochondrial metabolism. Using fluorescent-tagged lipids and metabolite tracing, we show that loss of the PI5P4Ks significantly impairs lipid uptake and β-oxidation in the mitochondria. Further, loss of PI5P4Ks results in dramatic alterations in mitochondrial structural and functional integrity, which under nutrient deprivation is further exacerbated, causing cell death. Notably, inhibition of the PI5P4Ks in cancer cells and mouse tumor models leads to decreased cell viability and tumor growth, respectively. Together, these studies reveal an unexplored role for PI5P4Ks in preserving metabolic homeostasis, which is necessary for tumorigenesis.

Telomerase activation counteracts senescence and telomere erosion caused by uncontrolled proliferation. Epidermal growth factor receptor (EGFR) amplification drives proliferation while telomerase reverse transcriptase promoter (TERTp) mutations underlie telomerase reactivation through recruitment of GA-binding protein (GABP). EGFR amplification and TERTp mutations typically co-occur in glioblastoma, the most common and aggressive primary brain tumor. To determine if these two frequent alterations driving proliferation and immortality are functionally connected, we combine analyses of copy number, mRNA, and protein data from tumor tissue with pharmacologic and genetic perturbations. We demonstrate that proliferation arrest decreases TERT expression in a GABP-dependent manner and elucidate a critical proliferation-to-immortality pathway from EGFR to TERT expression selectively from the mutant TERTp through activation of AMP-mediated kinase (AMPK) and GABP upregulation. EGFR-AMPK signaling promotes telomerase activity and maintains telomere length. These results define how the tumor cell immortality mechanism keeps pace with persistent oncogene signaling and cell cycling.

TERT promoter mutations are a hallmark of glioblastoma (GBM). Accordingly, TERT and GABPB1, a subunit of the upstream mutant TERT promoter transcription factor GABP, are being considered as promising therapeutic targets in GBM. We recently reported that the expression of TERT or GABP1 modulates flux via the pentose phosphate pathway (PPP). Here, we investigated whether 13C magnetic resonance spectroscopy (MRS) of hyperpolarized (HP) δ- [1-13C]gluconolactone can serve to image the reduction in PPP flux following TERT or GABPB1 silencing. We investigated two different human GBM cell lines stably expressing shRNAs targeting TERT or GABPB1, as well as doxycycline-inducible shTERT or shGABPB1cells. MRS studies were performed on live cells and in vivo tumors, and dynamic sets of 13C MR spectra were acquired following injection of HP δ-[1-13C]gluconolactone. HP 6-phosphogluconolactone (6PG), the product of δ-[1-13C]gluconolactone via the PPP, was significantly reduced in TERT or GABPB1-silenced cells or tumors compared to controls in all our models. Furthermore, a positive correlation between TERT expression and 6PG levels was observed. Our data indicate that HP δ-[1-13C]gluconolactone, an imaging tool with translational potential, could serve to monitor TERT expression and its silencing with therapies that target either TERT or GABPB1 in mutant TERT promoter GBM patients.

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.

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

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.

Glioblastoma multiforme (GBM) is an invariably fatal central nervous system tumor despite treatment with surgery, radiation, and chemotherapy. Further insights into the molecular and cellular mechanisms that drive GBM formation are required to improve patient outcome. MicroRNAs are emerging as important regulators of cellular differentiation and proliferation, and have been implicated in the etiology of a variety of cancers, yet the role of microRNAs in GBM remains poorly understood. In this study, we investigated the role of microRNAs in regulating the differentiation and proliferation of neural stem cells and glioblastoma-multiforme tumor cells.

Epigenetic silencing involving the aberrant methylation of promoter region CpG islands is widely recognized as a tumor suppressor silencing mechanism in cancer. However, the molecular pathways underlying aberrant DNA methylation remain elusive. Recently we showed that, on a genome-wide level, CpG island loci differ in their intrinsic susceptibility to aberrant methylation and that this susceptibility can be predicted based on underlying sequence context. These data suggest that there are sequence/structural features that contribute to the protection from or susceptibility to aberrant methylation. Here we use motif elicitation coupled with classification techniques to identify DNA sequence motifs that selectively define methylation-prone or methylation-resistant CpG islands. Motifs common to 28 methylation-prone or 47 methylation-resistant CpG island-containing genomic fragments were determined using the MEME and MAST algorithms (). The five most discriminatory motifs derived from methylation-prone sequences were found to be associated with CpG islands in general and were nonrandomly distributed throughout the genome. In contrast, the eight most discriminatory motifs derived from the methylation-resistant CpG islands were randomly distributed throughout the genome. Interestingly, this latter group tended to associate with Alu and other repetitive sequences. Used together, the frequency of occurrence of these motifs successfully discriminated methylation-prone and methylation-resistant CpG island groups with an accuracy of 87% after 10-fold cross-validation. The motifs identified here are candidate methylation-targeting or methylation-protection DNA sequences.

Despite gradual legislative efforts to phase out flame retardants (FRs) from the marketplace, polybrominated diphenyl ethers (PBDEs) are still widely detected in human maternal and fetal tissues, eg, placenta, due to their continued global application in consumer goods and inherent biological persistence. Recent studies in rodents and human placental cell lines suggest that PBDEs directly cause placental toxicity. During pregnancy, trophoblasts play key roles in uterine invasion, vascular remodeling, and anchoring of the placenta-fetal unit to the mother. Thus, to study the potential consequences of PBDE exposures on human placental development, we used an in vitro model: primary villous cytotrophoblasts (CTBs). Following exposures, the endpoints that were evaluated included cytotoxicity, function (migration, invasion), the transcriptome, and the methylome. In a concentration-dependent manner, common PBDE congeners, BDE-47 and -99, significantly reduced cell viability and increased death. Upon exposures to sub-cytotoxic concentrations (≤ 5 µM), we observed BDE-47 accumulation in CTBs with limited evidence of metabolism. At a functional level, BDE-47 hindered the ability of CTBs to migrate and invade. Transcriptomic analyses of BDE-47 effects suggested concentration-dependent changes in gene expression, involving stress pathways, eg, inflammation and lipid/cholesterol metabolism as well as processes underlying trophoblast fate, eg, differentiation, migration, and vascular morphogenesis. In parallel assessments, BDE-47 induced low-level global increases in methylation of CpG islands, including a subset that were proximal to genes with roles in cell adhesion/migration. Thus, using a primary human CTB model, we showed that PBDEs induced alterations at cellular and molecular levels, which could adversely impact placental development.

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.

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.

Environmental cues provoke rapid transitions in gene expression to support growth and cellular plasticity through incompletely understood mechanisms. Lin28 RNA-binding proteins have evolutionarily conserved roles in post-transcriptional coordination of pro-growth gene expression, but signaling pathways allowing trophic stimuli to induce Lin28 have remained uncharacterized. We find that Lin28a protein exhibits rapid basal turnover in neurons and that mitogen-activated protein kinase (MAPK)-dependent phosphorylation of the RNA-silencing factor HIV TAR-RNA-binding protein (TRBP) promotes binding and stabilization of Lin28a, but not Lin28b, with an accompanying reduction in Lin28-regulated miRNAs, downstream of brain-derived neurotrophic factor (BDNF). Binding of Lin28a to TRBP in vitro is also enhanced by phospho-mimic TRBP. Further, phospho-TRBP recapitulates BDNF-induced neuronal dendritic spine growth in a Lin28a-dependent manner. Finally, we demonstrate MAPK-dependent TRBP and Lin28a induction, with physiological function in growth and survival, downstream of diverse growth factors in multiple primary cell types, supporting a broad role for this pathway in trophic responses.

Phosphatidylinositol (PI)regulating enzymes are frequently altered in cancer and have become a focus for drug development. Here, we explore the phosphatidylinositol-5-phosphate 4-kinases (PI5P4K), a family of lipid kinases that regulate pools of intracellular PI, and demonstrate that the PI5P4Kα isoform influences androgen receptor (AR) signaling, which supports prostate cancer (PCa) cell survival. The regulation of PI becomes increasingly important in the setting of metabolic stress adaptation of PCa during androgen deprivation (AD), as we show that AD influences PI abundance and enhances intracellular pools of PI-4,5-P2. We suggest that this PI5P4Kα-AR relationship is mitigated through mTORC1 dysregulation and show that PI5P4Kα colocalizes to the lysosome, the intracellular site of mTORC1 complex activation. Notably, this relationship becomes prominent in mouse prostate tissue following surgical castration. Finally, multiple PCa cell models demonstrate marked survival vulnerability following stable PI5P4Kα inhibition. These results nominate PI5P4Kα as a target to disrupt PCa metabolic adaptation to castrate resistance.

Mutations in the TERT promoter represent the genetic underpinnings of tumor cell immortality. Beyond the two most common point mutations, which selectively recruit the ETS factor GABP to activate TERT, the significance of other variants is unknown. In seven cancer types, we identify duplications of wildtype sequence within the core promoter region of TERT that have strikingly similar features including an ETS motif, the duplication length and insertion site. The duplications recruit a GABP tetramer by virtue of the native ETS motif and its precisely spaced duplicated counterpart, activate the promoter and are clonal in a TERT expressing multifocal glioblastoma. We conclude that recurrent TERT promoter duplications 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.

We present Multi-miR, a microRNA-embedded shRNA system modeled after endogenous microRNA clusters that enables simultaneous expression of up to three or four short hairpin RNAs (shRNAs) from a single promoter without loss of activity, enabling robust combinatorial RNA interference (RNAi). We further developed complementary all-in-one vectors that are over one log-scale more sensitive to doxycycline-mediated activation in vitro than previous methods and resistant to shRNA inactivation in vivo. We demonstrate the utility of this system for intracranial expression of shRNAs in a glioblastoma model. Additionally, we leverage this platform to target the redundant RAF signaling node in a mouse model of KRAS-mutant cancer and show that robust combinatorial synthetic lethality efficiently abolishes tumor growth.