Genome-wide profiling of 5-formylcytosine reveals its roles in epigenetic priming.

TET proteins oxidize 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). 5fC and 5caC are excised by mammalian DNA glycosylase TDG, implicating 5mC oxidation in DNA demethylation. Here, we show that the genomic locations of 5fC can be determined by coupling chemical reduction with biotin tagging. Genome-wide mapping of 5fC in mouse embryonic stem cells (mESCs) reveals that 5fC preferentially occurs at poised enhancers among other gene regulatory elements. Application to Tdg null mESCs further suggests that 5fC production coordinates with p300 in remodeling epigenetic states of enhancers. This process, which is not influenced by 5hmC, appears to be associated with further oxidation of 5hmC and commitment to demethylation through 5fC. Finally, we resolved 5fC at base resolution by hydroxylamine-based protection from bisulfite-mediated deamination, thereby confirming sites of 5fC accumulation. Our results reveal roles of active 5mC/5hmC oxidation and TDG-mediated demethylation in epigenetic tuning at regulatory elements.

Base-resolution analysis of 5-hydroxymethylcytosine in the mammalian genome.

The study of 5-hydroxylmethylcytosines (5hmC) has been hampered by the lack of a method to map it at single-base resolution on a genome-wide scale. Affinity purification-based methods cannot precisely locate 5hmC nor accurately determine its relative abundance at each modified site. We here present a genome-wide approach, Tet-assisted bisulfite sequencing (TAB-Seq), that when combined with traditional bisulfite sequencing can be used for mapping 5hmC at base resolution and quantifying the relative abundance of 5hmC as well as 5mC. Application of this method to embryonic stem cells not only confirms widespread distribution of 5hmC in the mammalian genome but also reveals sequence bias and strand asymmetry at 5hmC sites. We observe high levels of 5hmC and reciprocally low levels of 5mC near but not on transcription factor-binding sites. Additionally, the relative abundance of 5hmC varies significantly among distinct functional sequence elements, suggesting different mechanisms for 5hmC deposition and maintenance.

Structural dynamics control the MicroRNA maturation pathway.

MicroRNAs (miRNAs) are crucial gene expression regulators and first-order suspects in the development and progression of many diseases. Comparative analysis of cancer cell expression data highlights many deregulated miRNAs. Low expression of miR-125a was related to poor breast cancer prognosis. Interestingly, a single nucleotide polymorphism (SNP) in miR-125a was located within a minor allele expressed by breast cancer patients. The SNP is not predicted to affect the ground state structure of the primary transcript or precursor, but neither the precursor nor mature product is detected by RT-qPCR. How this SNP modulates the maturation of miR-125a is poorly understood. Here, building upon a model of RNA dynamics derived from nuclear magnetic resonance studies, we developed a quantitative model enabling the visualization and comparison of networks of transient structures. We observed a high correlation between the distances between networks of variants with that of their respective wild types and their relative degrees of maturation to the latter, suggesting an important role of transient structures in miRNA homeostasis. We classified the human miRNAs according to pairwise distances between their networks of transient structures.

Base-resolution methylation patterns accurately predict transcription factor bindings in vivo.

Detecting in vivo transcription factor (TF) binding is important for understanding gene regulatory circuitries. ChIP-seq is a powerful technique to empirically define TF binding in vivo. However, the multitude of distinct TFs makes genome-wide profiling for them all labor-intensive and costly. Algorithms for in silico prediction of TF binding have been developed, based mostly on histone modification or DNase I hypersensitivity data in conjunction with DNA motif and other genomic features. However, technical limitations of these methods prevent them from being applied broadly, especially in clinical settings. We conducted a comprehensive survey involving multiple cell lines, TFs, and methylation types and found that there are intimate relationships between TF binding and methylation level changes around the binding sites. Exploiting the connection between DNA methylation and TF binding, we proposed a novel supervised learning approach to predict TF-DNA interaction using data from base-resolution whole-genome methylation sequencing experiments. We devised beta-binomial models to characterize methylation data around TF binding sites and the background. Along with other static genomic features, we adopted a random forest framework to predict TF-DNA interaction. After conducting comprehensive tests, we saw that the proposed method accurately predicts TF binding and performs favorably versus competing methods.

Regulation of DNA methylation turnover at LTR retrotransposons and imprinted loci by the histone methyltransferase Setdb1.

During mammalian development, DNA methylation patterns need to be reset in primordial germ cells (PGCs) and preimplantation embryos. However, many LTR retrotransposons and imprinted genes are impervious to such global epigenetic reprogramming via hitherto undefined mechanisms. Here, we report that a subset of such genomic regions are resistant to widespread erasure of DNA methylation in mouse embryonic stem cells (mESCs) lacking the de novo DNA methyltransferases (Dnmts) Dnmt3a and Dnmt3b. Intriguingly, these loci are enriched for H3K9me3 in mESCs, implicating this mark in DNA methylation homeostasis. Indeed, deletion of the H3K9 methyltransferase SET domain bifurcated 1 (Setdb1) results in reduced H3K9me3 and DNA methylation levels at specific loci, concomitant with increased 5-hydroxymethylation (5hmC) and ten-eleven translocation 1 binding. Taken together, these data reveal that Setdb1 promotes the persistence of DNA methylation in mESCs, likely reflecting one mechanism by which DNA methylation is maintained at LTR retrotransposons and imprinted genes during developmental stages when DNA methylation is reprogrammed.

Integrating DNA methylation dynamics into a framework for understanding epigenetic codes.

Genomic function is dictated by a combination of DNA sequence and the molecular mechanisms controlling access to genetic information. Access to DNA can be determined by the interpretation of covalent modifications that influence the packaging of DNA into chromatin, including DNA methylation and histone modifications. These modifications are believed to be forms of "epigenetic codes" that exist in discernable combinations that reflect cellular phenotype. Although DNA methylation is known to play important roles in gene regulation and genomic function, its contribution to the encoding of epigenetic information is just beginning to emerge. Here we discuss paradigms associated with the various components of DNA methylation/demethylation and recent advances in the understanding of its dynamic regulation in the genome, integrating these mechanisms into a framework to explain how DNA methylation could contribute to epigenetic codes.

TET1 plays an essential oncogenic role in MLL-rearranged leukemia.

The ten-eleven translocation 1 (TET1) gene is the founding member of the TET family of enzymes (TET1/2/3) that convert 5-methylcytosine to 5-hydroxymethylcytosine. Although TET1 was first identified as a fusion partner of the mixed lineage leukemia (MLL) gene in acute myeloid leukemia carrying t(10,11), its definitive role in leukemia is unclear. In contrast to the frequent down-regulation (or loss-of-function mutations) and critical tumor-suppressor roles of the three TET genes observed in various types of cancers, here we show that TET1 is a direct target of MLL-fusion proteins and is significantly up-regulated in MLL-rearranged leukemia, leading to a global increase of 5-hydroxymethylcytosine level. Furthermore, our both in vitro and in vivo functional studies demonstrate that Tet1 plays an indispensable oncogenic role in the development of MLL-rearranged leukemia, through coordination with MLL-fusion proteins in regulating their critical cotargets, including homeobox A9 (Hoxa9)/myeloid ecotropic viral integration 1 (Meis1)/pre-B-cell leukemia homeobox 3 (Pbx3) genes. Collectively, our data delineate an MLL-fusion/Tet1/Hoxa9/Meis1/Pbx3 signaling axis in MLL-rearranged leukemia and highlight TET1 as a potential therapeutic target in treating this presently therapy-resistant disease.

U1 small nuclear ribonucleoprotein complex and RNA splicing alterations in Alzheimer's disease.

Deposition of insoluble protein aggregates is a hallmark of neurodegenerative diseases. The universal presence of ß-amyloid and tau in Alzheimer's disease (AD) has facilitated advancement of the amyloid cascade and tau hypotheses that have dominated AD pathogenesis research and therapeutic development. However, the underlying etiology of the disease remains to be fully elucidated. Here we report a comprehensive study of the human brain-insoluble proteome in AD by mass spectrometry. We identify 4,216 proteins, among which 36 proteins accumulate in the disease, including U1-70K and other U1 small nuclear ribonucleoprotein (U1 snRNP) spliceosome components. Similar accumulations in mild cognitive impairment cases indicate that spliceosome changes occur in early stages of AD. Multiple U1 snRNP subunits form cytoplasmic tangle-like structures in AD but not in other examined neurodegenerative disorders, including Parkinson disease and frontotemporal lobar degeneration. Comparison of RNA from AD and control brains reveals dysregulated RNA processing with accumulation of unspliced RNA species in AD, including myc box-dependent-interacting protein 1, clusterin, and presenilin-1. U1-70K knockdown or antisense oligonucleotide inhibition of U1 snRNP increases the protein level of amyloid precursor protein. Thus, our results demonstrate unique U1 snRNP pathology and implicate abnormal RNA splicing in AD pathogenesis.

Sjogren syndrome antigen B (SSB)/La promotes global microRNA expression by binding microRNA precursors through stem-loop recognition.

MicroRNAs (miRNA) control numerous physiological and pathological processes. Typically, the primary miRNA (pri-miRNA) transcripts are processed by nuclear Drosha complex into ~70-nucleotide stem-loop precursor miRNAs (pre-miRNA), which are further cleaved by cytoplasmic Dicer complex into ~21-nucleotide mature miRNAs. However, it is unclear how nascent pre-miRNAs are protected from ribonucleases, such as MCPIP1, that degrade pre-miRNAs to abort miRNA production. Here, we identify Sjögren syndrome antigen B (SSB)/La as a pre-miRNA-binding protein that regulates miRNA processing in vitro. All three RNA-binding motifs (LAM, RRM1, and RRM2) of La/SSB are required for efficient pre-miRNA binding. Intriguingly, La/SSB recognizes the characteristic stem-loop structure of pre-miRNAs, of which the majority lack a 3' UUU terminus. Moreover, La/SSB associates with endogenous pri-/pre-miRNAs and promotes miRNA biogenesis by stabilizing pre-miRNAs from nuclease (e.g. MCPIP1)-mediated decay in mammalian cells. Accordingly, we observed positive correlations between the expression status of La/SSB and Dicer in human cancer transcriptome and prognosis. These studies identify an important function of La/SSB as a global regulator of miRNA expression, and implicate stem-loop recognition as a major mechanism that mediates association between La/SSB and diverse RNA molecules.

Genome-wide DNA hydroxymethylation changes are associated with neurodevelopmental genes in the developing human cerebellum.

5-Hydroxymethylcytosine (5-hmC) is a newly discovered modified form of cytosine that has been suspected to be an important epigenetic modification in neurodevelopment. While DNA methylation dynamics have already been implicated during neurodevelopment, little is known about hydroxymethylation in this process. Here, we report DNA hydroxymethylation dynamics during cerebellum development in the human brain. Overall, we find a positive correlation between 5-hmC levels and cerebellum development. Genome-wide profiling reveals that 5-hmC is highly enriched on specific gene regions including exons and especially the untranslated regions (UTRs), but it is depleted on introns and intergenic regions. Furthermore, we have identified fetus-specific and adult-specific differentially hydroxymethylated regions (DhMRs), most of which overlap with genes and CpG island shores. Surprisingly, during development, DhMRs are highly enriched in genes encoding mRNAs that can be regulated by fragile X mental retardation protein (FMRP), some of which are disrupted in autism, as well as in many known autism genes. Our results suggest that 5-hmC-mediated epigenetic regulation may broadly impact the development of the human brain, and its dysregulation could contribute to the molecular pathogenesis of neurodevelopmental disorders. Accession number: Sequencing data have been deposited to GEO with accession number GSE40539.

Integrating 5-hydroxymethylcytosine into the epigenomic landscape of human embryonic stem cells.

Covalent modification of DNA distinguishes cellular identities and is crucial for regulating the pluripotency and differentiation of embryonic stem (ES) cells. The recent demonstration that 5-methylcytosine (5-mC) may be further modified to 5-hydroxymethylcytosine (5-hmC) in ES cells has revealed a novel regulatory paradigm to modulate the epigenetic landscape of pluripotency. To understand the role of 5-hmC in the epigenomic landscape of pluripotent cells, here we profile the genome-wide 5-hmC distribution and correlate it with the genomic profiles of 11 diverse histone modifications and six transcription factors in human ES cells. By integrating genomic 5-hmC signals with maps of histone enrichment, we link particular pluripotency-associated chromatin contexts with 5-hmC. Intriguingly, through additional correlations with defined chromatin signatures at promoter and enhancer subtypes, we show distinct enrichment of 5-hmC at enhancers marked with H3K4me1 and H3K27ac. These results suggest potential role(s) for 5-hmC in the regulation of specific promoters and enhancers. In addition, our results provide a detailed epigenomic map of 5-hmC from which to pursue future functional studies on the diverse regulatory roles associated with 5-hmC.

Chemical modification-assisted bisulfite sequencing (CAB-Seq) for 5-carboxylcytosine detection in DNA.

5-Methylcytosine (5mC) in DNA can be oxidized stepwise to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) by the TET family proteins. Thymine DNA glycosylase can further remove 5fC and 5caC, connecting 5mC oxidation with active DNA demethylation. Here, we present a chemical modification-assisted bisulfite sequencing (CAB-Seq) that can detect 5caC with single-base resolution in DNA. We optimized 1-ethyl-3-[3-dimethylaminopropyl]carbodiimide hydrochloride (EDC)-catalyzed amide bond formation between the carboxyl group of 5caC and a primary amine group. We found that the modified 5caC can survive the bisulfite treatment without deamination. Therefore, this chemical labeling coupled with bisulfite treatment provides a base-resolution detection and sequencing method for 5caC.

Fragile x mental retardation protein regulates proliferation and differentiation of adult neural stem/progenitor cells.

Fragile X syndrome (FXS), the most common form of inherited mental retardation, is caused by the loss of functional fragile X mental retardation protein (FMRP). FMRP is an RNA-binding protein that can regulate the translation of specific mRNAs. Adult neurogenesis, a process considered important for neuroplasticity and memory, is regulated at multiple molecular levels. In this study, we investigated whether Fmrp deficiency affects adult neurogenesis. We show that in a mouse model of fragile X syndrome, adult neurogenesis is indeed altered. The loss of Fmrp increases the proliferation and alters the fate specification of adult neural progenitor/stem cells (aNPCs). We demonstrate that Fmrp regulates the protein expression of several components critical for aNPC function, including CDK4 and GSK3beta. Dysregulation of GSK3beta led to reduced Wnt signaling pathway activity, which altered the expression of neurogenin1 and the fate specification of aNPCs. These data unveil a novel regulatory role for Fmrp and translational regulation in adult neurogenesis.

MicroRNA miR-137 regulates neuronal maturation by targeting ubiquitin ligase mind bomb-1.

The maturation of young neurons is regulated by complex mechanisms and dysregulation of this process is frequently found in neurodevepmental disorders. MicroRNAs have been implicated in several steps of neuronal maturation including dendritic and axonal growth, spine development, and synaptogenesis. We demonstrate that one brain-enriched microRNA, miR-137, has a significant role in regulating neuronal maturation. Overexpression of miR-137 inhibits dendritic morphogenesis, phenotypic maturation, and spine development both in brain and cultured primary neurons. On the other hand, a reduction in miR-137 had opposite effects. We further show that miR-137 targets the Mind bomb one (Mib1) protein through the conserved target site located in the 3' untranslated region of Mib1 messenger RNA. Mib1 is an ubiquitin ligase known to be important for neurodevelopment. We show that exogenously expressed Mib1 could partially rescue the phenotypes associated with miR-137 overexpression. These results demonstrate a novel miRNA-mediated mechanism involving miR-137 and Mib1 that function to regulate neuronal maturation and dendritic morphogenesis during development.

The loss of methyl-CpG binding protein 1 leads to autism-like behavioral deficits.

Methyl-CpG binding proteins (MBDs) are central components of DNA methylation-mediated epigenetic gene regulation. Alterations of epigenetic pathways are known to be associated with several neurodevelopmental disorders, particularly autism. Our previous studies showed that the loss of Mbd1 led to reduced hippocampal neurogenesis and impaired learning in mice. However, whether MBD1 regulates the autism-related cognitive functions remains unknown. Here we show that Mbd1 mutant (Mbd1(-/-)) mice exhibit several core deficits frequently associated with autism, including reduced social interaction, learning deficits, anxiety, defective sensory motor gating, depression and abnormal brain serotonin activity. Furthermore, we find that Mbd1 can directly regulate the expression of Htr2c, one of the serotonin receptors, by binding to its promoter, and the loss of Mbd1 led to elevated expression of Htr2c. Our results, therefore, demonstrate the importance of epigenetic regulation in mammalian brain development and cognitive functions. Understanding how the loss of Mbd1 could lead to autism-like behavioral phenotypes would reveal much-needed information about the molecular pathogenesis of autism.

Single-Cell Genetic Analysis Using Automated Microfluidics to Resolve Somatic Mosaicism.

Somatic mosaicism occurs throughout normal development and contributes to numerous disease etiologies, including tumorigenesis and neurological disorders. Intratumor genetic heterogeneity is inherent to many cancers, creating challenges for effective treatments. Unfortunately, analysis of bulk DNA masks subclonal phylogenetic architectures created by the acquisition and distribution of somatic mutations amongst cells. As a result, single-cell genetic analysis is becoming recognized as vital for accurately characterizing cancers. Despite this, methods for single-cell genetics are lacking. Here we present an automated microfluidic workflow enabling efficient cell capture, lysis, and whole genome amplification (WGA). We find that ~90% of the genome is accessible in single cells with improved uniformity relative to current single-cell WGA methods. Allelic dropout (ADO) rates were limited to 13.75% and variant false discovery rates (SNV FDR) were 4.11x10(-6), on average. Application to ER-/PR-/HER2+ breast cancer cells and matched normal controls identified novel mutations that arose in a subpopulation of cells and effectively resolved the segregation of known cancer-related mutations with single-cell resolution. Finally, we demonstrate effective cell classification using mutation profiles with 10X average exome coverage depth per cell. Our data demonstrate an efficient automated microfluidic platform for single-cell WGA that enables the resolution of somatic mutation patterns in single cells.

Role of Tet1 and 5-hydroxymethylcytosine in cocaine action.

Ten-eleven translocation (TET) enzymes mediate the conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), which is enriched in brain, and its ultimate DNA demethylation. However, the influence of TET and 5hmC on gene transcription in brain remains elusive. We found that ten-eleven translocation protein 1 (TET1) was downregulated in mouse nucleus accumbens (NAc), a key brain reward structure, by repeated cocaine administration, which enhanced behavioral responses to cocaine. We then identified 5hmC induction in putative enhancers and coding regions of genes that have pivotal roles in drug addiction. Such induction of 5hmC, which occurred similarly following TET1 knockdown alone, correlated with increased expression of these genes as well as with their alternative splicing in response to cocaine administration. In addition, 5hmC alterations at certain loci persisted for at least 1 month after cocaine exposure. Together, these reveal a previously unknown epigenetic mechanism of cocaine action and provide new insight into how 5hmC regulates transcription in brain in vivo.

AGO3 Slicer activity regulates mitochondria-nuage localization of Armitage and piRNA amplification.

In Drosophila melanogaster the reciprocal "Ping-Pong" cycle of PIWI-interacting RNA (piRNA)-directed RNA cleavage catalyzed by the endonuclease (or "Slicer") activities of the PIWI proteins Aubergine (Aub) and Argonaute3 (AGO3) has been proposed to expand the secondary piRNA population. However, the role of AGO3/Aub Slicer activity in piRNA amplification remains to be explored. We show that AGO3 Slicer activity is essential for piRNA amplification and that AGO3 inhibits the homotypic Aub:Aub Ping-Pong process in a Slicer-independent manner. We also find that expression of an AGO3 Slicer mutant causes ectopic accumulation of Armitage, a key component in the primary piRNA pathway, in the Drosophila melanogaster germline granules known as nuage. AGO3 also coexists and interacts with Armitage in the mitochondrial fraction. Furthermore, AGO3 acts in conjunction with the mitochondria-associated protein Zucchini to control the dynamic subcellular localization of Armitage between mitochondria and nuage in a Slicer-dependent fashion. Collectively, our findings uncover a new mechanism that couples mitochondria with nuage to regulate secondary piRNA amplification.

Genome-wide antagonism between 5-hydroxymethylcytosine and DNA methylation in the adult mouse brain.

Mounting evidence points to critical roles for DNA modifications, including 5-methylcytosine (5mC) and its oxidized forms, in the development, plasticity and disorders of the mammalian nervous system. The novel DNA base 5-hydroxymethylcytosine (5hmC) is known to be capable of initiating passive or active DNA demethylation, but whether and how extensively 5hmC functions in shaping the post-mitotic neuronal DNA methylome is unclear. Here we report the genome-wide distribution of 5hmC in dentate granule neurons from adult mouse hippocampus in vivo. 5hmC in the neuronal genome is highly enriched in gene bodies, especially in exons, and correlates with gene expression. Direct genome-wide comparison of 5hmC distribution between embryonic stem cells and neurons reveals extensive differences, reflecting the functional disparity between these two cell types. Importantly, integrative analysis of 5hmC, overall DNA methylation and gene expression profiles of dentate granule neurons in vivo reveals the genome-wide antagonism between these two states of cytosine modifications, supporting a role for 5hmC in shaping the neuronal DNA methylome by promoting active DNA demethylation.

Tet-mediated covalent labelling of 5-methylcytosine for its genome-wide detection and sequencing.

5-methylcytosine is an epigenetic mark that affects a broad range of biological functions in mammals. The chemically inert methyl group prevents direct labelling for subsequent affinity purification and detection. Therefore, most current approaches for the analysis of 5-methylcytosine still have limitations of being either density-biased, lacking in robustness and consistency, or incapable of analysing 5-methylcytosine specifically. Here we present an approach, TAmC-Seq, which selectively tags 5-methylcytosine with an azide functionality that can be further labelled with a biotin for affinity purification, detection and genome-wide mapping. Using this covalent labelling approach, we demonstrate high sensitivity and specificity for known methylated loci, as well as increased CpG dinucleotide coverage at lower sequencing depth as compared with antibody-based enrichment, providing an improved efficiency in the 5-methylcytosine enrichment and genome-wide profiling.