Single-nucleus RNA sequencing of pre-malignant liver reveals disease-associated hepatocyte state with HCC prognostic potential.

Current approaches to staging chronic liver diseases have limited utility for predicting liver cancer risk. Here, we employed single-nucleus RNA sequencing (snRNA-seq) to characterize the cellular microenvironment of healthy and pre-malignant livers using two distinct mouse models. Downstream analyses unraveled a previously uncharacterized disease-associated hepatocyte (daHep) transcriptional state. These cells were absent in healthy livers but increasingly prevalent as chronic liver disease progressed. Copy number variation (CNV) analysis of microdissected tissue demonstrated that daHep-enriched regions are riddled with structural variants, suggesting these cells represent a pre-malignant intermediary. Integrated analysis of three recent human snRNA-seq datasets confirmed the presence of a similar phenotype in human chronic liver disease and further supported its enhanced mutational burden. Importantly, we show that high daHep levels precede carcinogenesis and predict a higher risk of hepatocellular carcinoma development. These findings may change the way chronic liver disease patients are staged, surveilled, and risk stratified.

Synthetic Epigenetic Reprogramming of Mesenchymal to Epithelial States Using the CRISPR/dCas9 Platform in Triple Negative Breast Cancer.

Epithelial-mesenchymal transition (EMT) is a reversible transcriptional program invoked by cancer cells to drive cancer progression. Transcription factor ZEB1 is a master regulator of EMT, driving disease recurrence in poor-outcome triple negative breast cancers (TNBCs). Here, this work silences ZEB1 in TNBC models by CRISPR/dCas9-mediated epigenetic editing, resulting in highly-specific and nearly complete suppression of ZEB1 in vivo, accompanied by long-lasting tumor inhibition. Integrated "omic" changes promoted by dCas9 linked to the KRAB domain (dCas9-KRAB) enabled the discovery of a ZEB1-dependent-signature of 26 genes differentially-expressed and -methylated, including the reactivation and enhanced chromatin accessibility in cell adhesion loci, outlining epigenetic reprogramming toward a more epithelial state. In the ZEB1 locus transcriptional silencing is associated with induction of locally-spread heterochromatin, significant changes in DNA methylation at specific CpGs, gain of H3K9me3, and a near complete erasure of H3K4me3 in the ZEB1 promoter. Epigenetic shifts induced by ZEB1-silencing are enriched in a subset of human breast tumors, illuminating a clinically-relevant hybrid-like state. Thus, the synthetic epi-silencing of ZEB1 induces stable "lock-in" epigenetic reprogramming of mesenchymal tumors associated with a distinct and stable epigenetic landscape. This work outlines epigenome-engineering approaches for reversing EMT and customizable precision molecular oncology approaches for targeting poor outcome breast cancers.

Retrotransposon instability dominates the acquired mutation landscape of mouse induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) can in principle differentiate into any cell of the body, and have revolutionized biomedical research and regenerative medicine. Unlike their human counterparts, mouse iPSCs (miPSCs) are reported to silence transposable elements and prevent transposable element-mediated mutagenesis. Here we apply short-read or Oxford Nanopore Technologies long-read genome sequencing to 38 bulk miPSC lines reprogrammed from 10 parental cell types, and 18 single-cell miPSC clones. While single nucleotide variants and structural variants restricted to miPSCs are rare, we find 83 de novo transposable element insertions, including examples intronic to Brca1 and Dmd. LINE-1 retrotransposons are profoundly hypomethylated in miPSCs, beyond other transposable elements and the genome overall, and harbor alternative protein-coding gene promoters. We show that treatment with the LINE-1 inhibitor lamivudine does not hinder reprogramming and efficiently blocks endogenous retrotransposition, as detected by long-read genome sequencing. These experiments reveal the complete spectrum and potential significance of mutations acquired by miPSCs.

Genomic Targeting of TET Activity for Targeted Demethylation Using CRISPR/Cas9.

Methylation of DNA at cytosine bases is an important DNA modification underlying normal development and disease states. Despite decades of research into the biological function of DNA methylation, most of the observations so far have relied primarily on associative data between observed changes in DNA methylation states and local changes in transcriptional activity or chromatin state processes. This is primarily due to the lack of molecular tools to precisely modify DNA methylation in the genome. Recent advances in genome editing technologies have allowed repurposing the CRISPR-Cas9 system for epigenome editing by fusing the catalytically dead Cas9 (dCas9) to epigenome modifying enzymes. Moreover, methods of recruiting multiple protein domains, including the SunTag system, have increased the efficacy of epigenome editing at target sites. Here, we describe an end-to-end protocol for efficient targeted removal of DNA methylation by recruiting multiple catalytic domain of TET1 enzymes to the target sites with the dCas9-SunTag system, including sgRNA design, molecular cloning, delivery of plasmid into mammalian cells, and targeted DNA methylation analysis.

Cerebrospinal fluid liquid biopsy for detecting somatic mosaicism in brain.

Brain somatic mutations are an increasingly recognized cause of epilepsy, brain malformations and autism spectrum disorders and may be a hidden cause of other neurodevelopmental and neurodegenerative disorders. At present, brain mosaicism can be detected only in the rare situations of autopsy or brain biopsy. Liquid biopsy using cell-free DNA derived from cerebrospinal fluid has detected somatic mutations in malignant brain tumours. Here, we asked if cerebrospinal fluid liquid biopsy can be used to detect somatic mosaicism in non-malignant brain diseases. First, we reliably quantified cerebrospinal fluid cell-free DNA in 28 patients with focal epilepsy and 28 controls using droplet digital PCR. Then, in three patients we identified somatic mutations in cerebrospinal fluid: in one patient with subcortical band heterotopia the LIS1 p. Lys64* variant at 9.4% frequency; in a second patient with focal cortical dysplasia the TSC1 p. Phe581His*6 variant at 7.8% frequency; and in a third patient with ganglioglioma the BRAF p. Val600Glu variant at 3.2% frequency. To determine if cerebrospinal fluid cell-free DNA was brain-derived, whole-genome bisulphite sequencing was performed and brain-specific DNA methylation patterns were found to be significantly enriched (P = 0.03). Our proof of principle study shows that cerebrospinal fluid liquid biopsy is valuable in investigating mosaic neurological disorders where brain tissue is unavailable.

Retention of paternal DNA methylome in the developing zebrafish germline.

Two waves of DNA methylation reprogramming occur during mammalian embryogenesis; during preimplantation development and during primordial germ cell (PGC) formation. However, it is currently unclear how evolutionarily conserved these processes are. Here we characterise the DNA methylomes of zebrafish PGCs at four developmental stages and identify retention of paternal epigenetic memory, in stark contrast to the findings in mammals. Gene expression profiling of zebrafish PGCs at the same developmental stages revealed that the embryonic germline is defined by a small number of markers that display strong developmental stage-specificity and that are independent of DNA methylation-mediated regulation. We identified promoters that are specifically targeted by DNA methylation in somatic and germline tissues during vertebrate embryogenesis and that are frequently misregulated in human cancers. Together, these detailed methylome and transcriptome maps of the zebrafish germline provide insight into vertebrate DNA methylation reprogramming and enhance our understanding of the relationships between germline fate acquisition and oncogenesis.

Epigenomic Landscapes of hESC-Derived Neural Rosettes: Modeling Neural Tube Formation and Diseases.

We currently lack a comprehensive understanding of the mechanisms underlying neural tube formation and their contributions to neural tube defects (NTDs). Developing a model to study such a complex morphogenetic process, especially one that models human-specific aspects, is critical. Three-dimensional, human embryonic stem cell (hESC)-derived neural rosettes (NRs) provide a powerful resource for in vitro modeling of human neural tube formation. Epigenomic maps reveal enhancer elements unique to NRs relative to 2D systems. A master regulatory network illustrates that key NR properties are related to their epigenomic landscapes. We found that folate-associated DNA methylation changes were enriched within NR regulatory elements near genes involved in neural tube formation and metabolism. Our comprehensive regulatory maps offer insights into the mechanisms by which folate may prevent NTDs. Lastly, our distal regulatory maps provide a better understanding of the potential role of neurological-disorder-associated SNPs.

Synthetically controlling dendrimer flexibility improves delivery of large plasmid DNA.

Tools for editing the genome and epigenome have revolutionised the field of molecular biology and represent a new frontier in targeted therapeutic intervention. Although efficiencies and specificities of genome editing technologies have improved with the development of TALEs and CRISPR platforms, intracellular delivery of these larger constructs still remains a challenge using existing delivery agents. Viral vectors, including lentiviruses and adeno-associated viruses, as well as some non-viral strategies, such as cationic polymers and liposomes, are limited by packaging capacity, poor delivery, toxicity, and immunogenicity. We report a highly controlled synthetic strategy to engineer a flexible dendritic polymer using click chemistry to overcome the aforementioned delivery challenges associated with genome engineering technologies. Using a systematic approach, we demonstrate that high transfection efficiencies and packaging capacity can be achieved using this non-viral delivery methodology to deliver zinc fingers, TALEs and CRISPR/dCas9 platforms.

Epigenomic Signatures of Neuronal Diversity in the Mammalian Brain.

Neuronal diversity is essential for mammalian brain function but poses a challenge to molecular profiling. To address the need for tools that facilitate cell-type-specific epigenomic studies, we developed the first affinity purification approach to isolate nuclei from genetically defined cell types in a mammal. We combine this technique with next-generation sequencing to show that three subtypes of neocortical neurons have highly distinctive epigenomic landscapes. Over 200,000 regions differ in chromatin accessibility and DNA methylation signatures characteristic of gene regulatory regions. By footprinting and motif analyses, these regions are predicted to bind distinct cohorts of neuron subtype-specific transcription factors. Neuronal epigenomes reflect both past and present gene expression, with DNA hyper-methylation at developmentally critical genes appearing as a novel epigenomic signature in mature neurons. Taken together, our findings link the functional and transcriptional complexity of neurons to their underlying epigenomic diversity.

Analysis of an artificial zinc finger epigenetic modulator: widespread binding but limited regulation.

Artificial transcription factors (ATFs) and genomic nucleases based on a DNA binding platform consisting of multiple zinc finger domains are currently being developed for clinical applications. However, no genome-wide investigations into their binding specificity have been performed. We have created six-finger ATFs to target two different 18 nt regions of the human SOX2 promoter; each ATF is constructed such that it contains or lacks a super KRAB domain (SKD) that interacts with a complex containing repressive histone methyltransferases. ChIP-seq analysis of the effector-free ATFs in MCF7 breast cancer cells identified thousands of binding sites, mostly in promoter regions; the addition of an SKD domain increased the number of binding sites ∼ 5-fold, with a majority of the new sites located outside of promoters. De novo motif analyses suggest that the lack of binding specificity is due to subsets of the finger domains being used for genomic interactions. Although the ATFs display widespread binding, few genes showed expression differences; genes repressed by the ATF-SKD have stronger binding sites and are more enriched for a 12 nt motif. Interestingly, epigenetic analyses indicate that the transcriptional repression caused by the ATF-SKD is not due to changes in active histone modifications.

Epigenomic analysis of multilineage differentiation of human embryonic stem cells.

Epigenetic mechanisms have been proposed to play crucial roles in mammalian development, but their precise functions are only partially understood. To investigate epigenetic regulation of embryonic development, we differentiated human embryonic stem cells into mesendoderm, neural progenitor cells, trophoblast-like cells, and mesenchymal stem cells and systematically characterized DNA methylation, chromatin modifications, and the transcriptome in each lineage. We found that promoters that are active in early developmental stages tend to be CG rich and mainly engage H3K27me3 upon silencing in nonexpressing lineages. By contrast, promoters for genes expressed preferentially at later stages are often CG poor and primarily employ DNA methylation upon repression. Interestingly, the early developmental regulatory genes are often located in large genomic domains that are generally devoid of DNA methylation in most lineages, which we termed DNA methylation valleys (DMVs). Our results suggest that distinct epigenetic mechanisms regulate early and late stages of ES cell differentiation.

Global DNA hypomethylation coupled to repressive chromatin domain formation and gene silencing in breast cancer.

While genetic mutation is a hallmark of cancer, many cancers also acquire epigenetic alterations during tumorigenesis including aberrant DNA hypermethylation of tumor suppressors, as well as changes in chromatin modifications as caused by genetic mutations of the chromatin-modifying machinery. However, the extent of epigenetic alterations in cancer cells has not been fully characterized. Here, we describe complete methylome maps at single nucleotide resolution of a low-passage breast cancer cell line and primary human mammary epithelial cells. We find widespread DNA hypomethylation in the cancer cell, primarily at partially methylated domains (PMDs) in normal breast cells. Unexpectedly, genes within these regions are largely silenced in cancer cells. The loss of DNA methylation in these regions is accompanied by formation of repressive chromatin, with a significant fraction displaying allelic DNA methylation where one allele is DNA methylated while the other allele is occupied by histone modifications H3K9me3 or H3K27me3. Our results show a mutually exclusive relationship between DNA methylation and H3K9me3 or H3K27me3. These results suggest that global DNA hypomethylation in breast cancer is tightly linked to the formation of repressive chromatin domains and gene silencing, thus identifying a potential epigenetic pathway for gene regulation in cancer cells.

Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells.

5-hydroxymethylcytosine (5hmC) is a modified base present at low levels in diverse cell types in mammals. 5hmC is generated by the TET family of Fe(II) and 2-oxoglutarate-dependent enzymes through oxidation of 5-methylcytosine (5mC). 5hmC and TET proteins have been implicated in stem cell biology and cancer, but information on the genome-wide distribution of 5hmC is limited. Here we describe two novel and specific approaches to profile the genomic localization of 5hmC. The first approach, termed GLIB (glucosylation, periodate oxidation, biotinylation) uses a combination of enzymatic and chemical steps to isolate DNA fragments containing as few as a single 5hmC. The second approach involves conversion of 5hmC to cytosine 5-methylenesulphonate (CMS) by treatment of genomic DNA with sodium bisulphite, followed by immunoprecipitation of CMS-containing DNA with a specific antiserum to CMS. High-throughput sequencing of 5hmC-containing DNA from mouse embryonic stem (ES) cells showed strong enrichment within exons and near transcriptional start sites. 5hmC was especially enriched at the start sites of genes whose promoters bear dual histone 3 lysine 27 trimethylation (H3K27me3) and histone 3 lysine 4 trimethylation (H3K4me3) marks. Our results indicate that 5hmC has a probable role in transcriptional regulation, and suggest a model in which 5hmC contributes to the 'poised' chromatin signature found at developmentally-regulated genes in ES cells.

The Arabidopsis glutathione transferase gene family displays complex stress regulation and co-silencing multiple genes results in altered metabolic sensitivity to oxidative stress.

Plant glutathione transferases (GSTs) are induced by diverse biotic and abiotic stimuli, and are important for protecting plants against oxidative damage. We have studied the primary transcriptional stress response of the entire Arabidopsis GST family to seven stresses, including both biotic and abiotic stimuli, with a focus on early changes in gene expression. Our results indicate that individual GST genes are highly specific in their induction patterns. Furthermore, we have been able to link individual GSTs to particular stress stimuli. Using RNAi, we successfully co-silenced a group of four phi GSTs that represent some of the most highly expressed GST genes. Despite a marked reduction in total phi GST protein levels, the transgenic plants showed no reduction in GST activity as measured using the model substrate 1-chloro-2,4-dinitrobenzene (CDNB), and appeared to have surprisingly robust physical phenotypes during stress. However, analysis of metabolite pools showed oxidation of the glutathione pool in the RNAi lines, and we observed alterations in carbon and nitrogen compounds following salicylic acid and hydrogen peroxide stress treatments, indicative of oxidative modification of primary metabolism. Thus, there appears to be a high degree of functional redundancy within the Arabidopsis GST family, with extensive disruption being required to reveal the roles of phi GSTs in protection against oxidative stress.