Non-CG DNA methylation and MeCP2 stabilize repeated tuning of long genes that distinguish closely related neuron types.

The extraordinary diversity of neuron types in the mammalian brain is delineated at the highest resolution by subtle gene expression differences that may require specialized molecular mechanisms to be maintained. Neurons uniquely express the longest genes in the genome and utilize neuron-enriched non-CG DNA methylation (mCA) together with the Rett syndrome protein, MeCP2, to control gene expression, but the function of these unique gene structures and machinery in regulating finely resolved neuron type-specific gene programs has not been explored. Here, we employ epigenomic and spatial transcriptomic analyses to discover a major role for mCA and MeCP2 in maintaining neuron type-specific gene programs at the finest scale of cellular resolution. We uncover differential susceptibility to MeCP2 loss in neuronal populations depending on global mCA levels and dissect methylation patterns and intragenic enhancer repression that drive overlapping and distinct gene regulation between neuron types. Strikingly, we show that mCA and MeCP2 regulate genes that are repeatedly tuned to differentiate neuron types at the highest cellular resolution, including spatially resolved, vision-dependent gene programs in the visual cortex. These repeatedly tuned genes display genomic characteristics, including long length, numerous intragenic enhancers, and enrichment for mCA, that predispose them to regulation by MeCP2. Thus, long gene regulation by the MeCP2 pathway maintains differential gene expression between closely-related neurons to facilitate the exceptional cellular diversity in the complex mammalian brain.

MeCP2 represses the activity of topoisomerase IIß in long neuronal genes.

A unique signature of neurons is the high expression of the longest genes in the genome. These genes have essential neuronal functions, and disruption of their expression has been implicated in neurological disorders. DNA topoisomerases resolve DNA topological constraints and facilitate neuronal long gene expression. Conversely, the Rett syndrome protein, methyl-CpG-binding protein 2 (MeCP2), can transcriptionally repress long genes. How these factors regulate long genes is not well understood, and whether they interact is not known. Here, we identify and map a functional interaction between MeCP2 and topoisomerase IIß (TOP2ß) in mouse neurons. We profile neuronal TOP2ß activity genome wide, detecting enrichment at regulatory regions and gene bodies of long genes, including MeCP2-regulated genes. We show that loss and overexpression of MeCP2 alter TOP2ß activity at MeCP2-regulated genes. These findings uncover a mechanism of TOP2ß inhibition by MeCP2 in neurons and implicate TOP2ß dysregulation in disorders caused by MeCP2 disruption.

Distinct disease mutations in DNMT3A result in a spectrum of behavioral, epigenetic, and transcriptional deficits.

Phenotypic heterogeneity in monogenic neurodevelopmental disorders can arise from differential severity of variants underlying disease, but how distinct alleles drive variable disease presentation is not well understood. Here, we investigate missense mutations in DNA methyltransferase 3A (DNMT3A), a DNA methyltransferase associated with overgrowth, intellectual disability, and autism, to uncover molecular correlates of phenotypic heterogeneity. We generate a Dnmt3aP900L/+ mouse mimicking a mutation with mild to moderate severity and compare phenotypic and epigenomic effects with a severe R878H mutation. P900L mutants exhibit core growth and behavioral phenotypes shared across models but show subtle epigenomic changes, while R878H mutants display extensive disruptions. We identify mutation-specific dysregulated genes that may contribute to variable disease severity. Shared transcriptomic disruption identified across mutations overlaps dysregulation observed in other developmental disorder models and likely drives common phenotypes. Together, our findings define central drivers of DNMT3A disorders and illustrate how variable epigenomic disruption contributes to phenotypic heterogeneity in neurodevelopmental disease.

Towards a comprehensive regulatory map of Mammalian Genomes.

Genome mapping studies have generated a nearly complete collection of genes for the human genome, but we still lack an equivalently vetted inventory of human regulatory sequences. Cis-regulatory modules (CRMs) play important roles in controlling when, where, and how much a gene is expressed. We developed a training data-free CRM-prediction algorithm, the Mammalian Regulatory MOdule Detector (MrMOD) for accurate CRM prediction in mammalian genomes. MrMOD provides genome position-fixed CRM models similar to the fixed gene models for the mouse and human genomes using only genomic sequences as the inputs with one adjustable parameter - the significance p-value. Importantly, MrMOD predicts a comprehensive set of high-resolution CRMs in the mouse and human genomes including all types of regulatory modules not limited to any tissue, cell type, developmental stage, or condition. We computationally validated MrMOD predictions used a compendium of 21 orthogonal experimental data sets including thousands of experimentally defined CRMs and millions of putative regulatory elements derived from hundreds of different tissues, cell types, and stimulus conditions obtained from multiple databases. In ovo transgenic reporter assay demonstrates the power of our prediction in guiding experimental design. We analyzed CRMs located in the chromosome 17 using unsupervised machine learning and identified groups of CRMs with multiple lines of evidence supporting their functionality, linking CRMs with upstream binding transcription factors and downstream target genes. Our work provides a comprehensive base pair resolution annotation of the functional regulatory elements and non-functional regions in the mammalian genomes.

NSD1 deposits histone H3 lysine 36 dimethylation to pattern non-CG DNA methylation in neurons.

During postnatal development, the DNA methyltransferase DNMT3A deposits high levels of non-CG cytosine methylation in neurons. This methylation is critical for transcriptional regulation, and loss of this mark is implicated in DNMT3A-associated neurodevelopmental disorders (NDDs). Here, we show in mice that genome topology and gene expression converge to shape histone H3 lysine 36 dimethylation (H3K36me2) profiles, which in turn recruit DNMT3A and pattern neuronal non-CG methylation. We show that NSD1, an H3K36 methyltransferase mutated in NDD, is required for the patterning of megabase-scale H3K36me2 and non-CG methylation in neurons. We find that brain-specific deletion of NSD1 causes altered DNA methylation that overlaps with DNMT3A disorder models to drive convergent dysregulation of key neuronal genes that may underlie shared phenotypes in NSD1- and DNMT3A-associated NDDs. Our findings indicate that H3K36me2 deposited by NSD1 is important for neuronal non-CG DNA methylation and suggest that the H3K36me2-DNMT3A-non-CG-methylation pathway is likely disrupted in NSD1-associated NDDs.

Distinct disease mutations in DNMT3A result in a spectrum of behavioral, epigenetic, and transcriptional deficits.

Phenotypic heterogeneity is a common feature of monogenic neurodevelopmental disorders that can arise from differential severity of missense variants underlying disease, but how distinct alleles impact molecular mechanisms to drive variable disease presentation is not well understood. Here, we investigate missense mutations in the DNA methyltransferase DNMT3A associated with variable overgrowth, intellectual disability, and autism, to uncover molecular correlates of phenotypic heterogeneity in neurodevelopmental disease. We generate a DNMT3A P900L/+ mouse model mimicking a disease mutation with mild-to-moderate severity and compare phenotypic and epigenomic effects with a severe R878H mutation. We show that the P900L mutation leads to disease-relevant overgrowth, obesity, and social deficits shared across DNMT3A disorder models, while the R878H mutation causes more extensive epigenomic disruption leading to differential dysregulation of enhancers elements. We identify distinct gene sets disrupted in each mutant which may contribute to mild or severe disease, and detect shared transcriptomic disruption that likely drives common phenotypes across affected individuals. Finally, we demonstrate that core gene dysregulation detected in DNMT3A mutant mice overlaps effects in other developmental disorder models, highlighting the importance of DNMT3A-deposited methylation in neurodevelopment. Together, these findings define central drivers of DNMT3A disorders and illustrate how variable disruption of transcriptional mechanisms can drive the spectrum of phenotypes in neurodevelopmental disease.

NSD1 deposits histone H3 lysine 36 dimethylation to pattern non-CG DNA methylation in neurons.

During postnatal development the DNA methyltransferase DNMT3A deposits high levels of non-CG cytosine methylation in neurons. This unique methylation is critical for transcriptional regulation in the mature mammalian brain, and loss of this mark is implicated in DNMT3A-associated neurodevelopmental disorders (NDDs). The mechanisms determining genomic non-CG methylation profiles are not well defined however, and it is unknown if this pathway is disrupted in additional NDDs. Here we show that genome topology and gene expression converge to shape histone H3 lysine 36 dimethylation (H3K36me2) profiles, which in turn recruit DNMT3A and pattern neuronal non-CG methylation. We show that NSD1, the H3K36 methyltransferase mutated in the NDD, Sotos syndrome, is required for megabase-scale patterning of H3K36me2 and non-CG methylation in neurons. We find that brain-specific deletion of NSD1 causes alterations in DNA methylation that overlap with models of DNMT3A disorders and define convergent disruption in the expression of key neuronal genes in these models that may contribute to shared phenotypes in NSD1- and DNMT3A-associated NDD. Our findings indicate that H3K36me2 deposited by NSD1 is an important determinant of neuronal non-CG DNA methylation and implicates disruption of this methylation in Sotos syndrome. Highlights: Topology-associated DNA methylation and gene expression independently contribute to neuronal gene body and enhancer non-CG DNA methylation patterns.Topology-associated H3K36me2 patterns and local enrichment of H3K4 methylation impact deposition of non-CG methylation by DNMT3A. Disruption of NSD1 in vivo leads to alterations in H3K36me2, DNA methylation, and gene expression that overlap with models of DNMT3A disorders.

APC7 mediates ubiquitin signaling in constitutive heterochromatin in the developing mammalian brain.

Neurodevelopmental cognitive disorders provide insights into mechanisms of human brain development. Here, we report an intellectual disability syndrome caused by the loss of APC7, a core component of the E3 ubiquitin ligase anaphase promoting complex (APC). In mechanistic studies, we uncover a critical role for APC7 during the recruitment and ubiquitination of APC substrates. In proteomics analyses of the brain from mice harboring the patient-specific APC7 mutation, we identify the chromatin-associated protein Ki-67 as an APC7-dependent substrate of the APC in neurons. Conditional knockout of the APC coactivator protein Cdh1, but not Cdc20, leads to the accumulation of Ki-67 protein in neurons in vivo, suggesting that APC7 is required for the function of Cdh1-APC in the brain. Deregulated neuronal Ki-67 upon APC7 loss localizes predominantly to constitutive heterochromatin. Our findings define an essential function for APC7 and Cdh1-APC in neuronal heterochromatin regulation, with implications for understanding human brain development and disease.

A MYT1L syndrome mouse model recapitulates patient phenotypes and reveals altered brain development due to disrupted neuronal maturation.

Human genetics have defined a new neurodevelopmental syndrome caused by loss-of-function mutations in MYT1L, a transcription factor known for enabling fibroblast-to-neuron conversions. However, how MYT1L mutation causes intellectual disability, autism, ADHD, obesity, and brain anomalies is unknown. Here, we developed a Myt1l haploinsufficient mouse model that develops obesity, white-matter thinning, and microcephaly, mimicking common clinical phenotypes. During brain development we discovered disrupted gene expression, mediated in part by loss of Myt1l gene-target activation, and identified precocious neuronal differentiation as the mechanism for microcephaly. In contrast, in adults we discovered that mutation results in failure of transcriptional and chromatin maturation, echoed in disruptions in baseline physiological properties of neurons. Myt1l haploinsufficiency also results in behavioral anomalies, including hyperactivity, muscle weakness, and social alterations, with more severe phenotypes in males. Overall, our findings provide insight into the mechanistic underpinnings of this disorder and enable future preclinical studies.

CHARGE syndrome protein CHD7 regulates epigenomic activation of enhancers in granule cell precursors and gyrification of the cerebellum.

Regulation of chromatin plays fundamental roles in the development of the brain. Haploinsufficiency of the chromatin remodeling enzyme CHD7 causes CHARGE syndrome, a genetic disorder that affects the development of the cerebellum. However, how CHD7 controls chromatin states in the cerebellum remains incompletely understood. Using conditional knockout of CHD7 in granule cell precursors in the mouse cerebellum, we find that CHD7 robustly promotes chromatin accessibility, active histone modifications, and RNA polymerase recruitment at enhancers. In vivo profiling of genome architecture reveals that CHD7 concordantly regulates epigenomic modifications associated with enhancer activation and gene expression of topologically-interacting genes. Genome and gene ontology studies show that CHD7-regulated enhancers are associated with genes that control brain tissue morphogenesis. Accordingly, conditional knockout of CHD7 triggers a striking phenotype of cerebellar polymicrogyria, which we have also found in a case of CHARGE syndrome. Finally, we uncover a CHD7-dependent switch in the preferred orientation of granule cell precursor division in the developing cerebellum, providing a potential cellular basis for the cerebellar polymicrogyria phenotype upon loss of CHD7. Collectively, our findings define epigenomic regulation by CHD7 in granule cell precursors and identify abnormal cerebellar patterning upon CHD7 depletion, with potential implications for our understanding of CHARGE syndrome.

Functional and epigenetic phenotypes of humans and mice with DNMT3A Overgrowth Syndrome.

Germline pathogenic variants in DNMT3A were recently described in patients with overgrowth, obesity, behavioral, and learning difficulties (DNMT3A Overgrowth Syndrome/DOS). Somatic mutations in the DNMT3A gene are also the most common cause of clonal hematopoiesis, and can initiate acute myeloid leukemia (AML). Using whole genome bisulfite sequencing, we studied DNA methylation in peripheral blood cells of 11 DOS patients and found a focal, canonical hypomethylation phenotype, which is most severe with the dominant negative DNMT3AR882H mutation. A germline mouse model expressing the homologous Dnmt3aR878H mutation phenocopies most aspects of the human DOS syndrome, including the methylation phenotype and an increased incidence of spontaneous hematopoietic malignancies, suggesting that all aspects of this syndrome are caused by this mutation.

Deconstructing Stepwise Fate Conversion of Human Fibroblasts to Neurons by MicroRNAs.

Cell-fate conversion generally requires reprogramming effectors to both introduce fate programs of the target cell type and erase the identity of starting cell population. Here, we reveal insights into the activity of microRNAs miR-9/9∗ and miR-124 (miR-9/9∗-124) as reprogramming agents that orchestrate direct conversion of human fibroblasts into motor neurons by first eradicating fibroblast identity and promoting uniform transition to a neuronal state in sequence. We identify KLF-family transcription factors as direct target genes for miR-9/9∗-124 and show their repression is critical for erasing fibroblast fate. Subsequent gain of neuronal identity requires upregulation of a small nuclear RNA, RN7SK, which induces accessibilities of chromatin regions and neuronal gene activation to push cells to a neuronal state. Our study defines deterministic components in the microRNA-mediated reprogramming cascade.

DNMT3A Haploinsufficiency Results in Behavioral Deficits and Global Epigenomic Dysregulation Shared across Neurodevelopmental Disorders.

Mutations in DNA methyltransferase 3A (DNMT3A) have been detected in autism and related disorders, but how these mutations disrupt nervous system function is unknown. Here, we define the effects of DNMT3A mutations associated with neurodevelopmental disease. We show that diverse mutations affect different aspects of protein activity but lead to shared deficiencies in neuronal DNA methylation. Heterozygous DNMT3A knockout mice mimicking DNMT3A disruption in disease display growth and behavioral alterations consistent with human phenotypes. Strikingly, in these mice, we detect global disruption of neuron-enriched non-CG DNA methylation, a binding site for the Rett syndrome protein MeCP2. Loss of this methylation leads to enhancer and gene dysregulation that overlaps with models of Rett syndrome and autism. These findings define the effects of DNMT3A haploinsufficiency in the brain and uncover disruption of the non-CG methylation pathway as a convergence point across neurodevelopmental disorders.

Emerging Insights into the Distinctive Neuronal Methylome.

The genomes of mammalian neurons are enriched for unique forms of DNA methylation, including exceptionally high levels of non-CG methylation. Here, we review recent studies defining how non-CG methylation accumulates in neurons and is read out by the critical regulator of neuronal transcription, MeCP2. We discuss the role of gene expression and genome architecture in establishing non-CG methylation and highlight emerging mechanistic insights into how non-CG methylation and MeCP2 control transcription. Further, we describe the cell type-specific functions of this methylation and explore growing evidence that disruption of this regulatory pathway contributes to neurodevelopmental disorders. These findings uncover how the distinctive epigenome in neurons facilitates the development and function of the complex mammalian brain.

The chromatin remodeling enzyme Chd4 regulates genome architecture in the mouse brain.

The development and function of the brain require tight control of gene expression. Genome architecture is thought to play a critical regulatory role in gene expression, but the mechanisms governing genome architecture in the brain in vivo remain poorly understood. Here, we report that conditional knockout of the chromatin remodeling enzyme Chd4 in granule neurons of the mouse cerebellum increases accessibility of gene regulatory sites genome-wide in vivo. Conditional knockout of Chd4 promotes recruitment of the architectural protein complex cohesin preferentially to gene enhancers in granule neurons in vivo. Importantly, in vivo profiling of genome architecture reveals that conditional knockout of Chd4 strengthens interactions among developmentally repressed contact domains as well as genomic loops in a manner that tightly correlates with increased accessibility, enhancer activity, and cohesin occupancy at these sites. Collectively, our findings define a role for chromatin remodeling in the control of genome architecture organization in the mammalian brain.

Functions of Gtf2i and Gtf2ird1 in the developing brain: transcription, DNA binding and long-term behavioral consequences.

Gtf2ird1 and Gtf2i are two transcription factors (TFs) among the 28 genes deleted in Williams syndrome, and prior mouse models of each TF show behavioral phenotypes. Here we identify their genomic binding sites in the developing brain and test for additive effects of their mutation on transcription and behavior. GTF2IRD1 binding targets were enriched for transcriptional and chromatin regulators and mediators of ubiquitination. GTF2I targets were enriched for signal transduction proteins, including regulators of phosphorylation and WNT. Both TFs are highly enriched at promoters, strongly overlap CTCF binding and topological associating domain boundaries and moderately overlap each other, suggesting epistatic effects. Shared TF targets are enriched for reactive oxygen species-responsive genes, synaptic proteins and transcription regulators such as chromatin modifiers, including a significant number of highly constrained genes and known ASD genes. We next used single and double mutants to test whether mutating both TFs will modify transcriptional and behavioral phenotypes of single Gtf2ird1 mutants, though with the caveat that our Gtf2ird1 mutants, like others previously reported, do produce low levels of a truncated protein product. Despite little difference in DNA binding and transcriptome-wide expression, homozygous Gtf2ird1 mutation caused balance, marble burying and conditioned fear phenotypes. However, mutating Gtf2i in addition to Gtf2ird1 did not further modify transcriptomic or most behavioral phenotypes, suggesting Gtf2ird1 mutation alone was sufficient for the observed phenotypes.

MeCP2 Represses Enhancers through Chromosome Topology-Associated DNA Methylation.

The genomes of mammalian neurons contain uniquely high levels of non-CG DNA methylation that can be bound by the Rett syndrome protein, MeCP2, to regulate gene expression. How patterns of non-CG methylation are established in neurons and the mechanism by which this methylation works with MeCP2 to control gene expression is unclear. Here, we find that genes repressed by MeCP2 are often located within megabase-scale regions of high non-CG methylation that correspond with topologically associating domains of chromatin folding. MeCP2 represses enhancers found in these domains that are enriched for non-CG and CG methylation, with the strongest repression occurring for enhancers located within MeCP2-repressed genes. These alterations in enhancer activity provide a mechanism for how MeCP2 disruption in disease can lead to widespread changes in gene expression. Hence, we find that DNA topology can shape non-CG DNA methylation across the genome to dictate MeCP2-mediated enhancer regulation in the brain.

ME-Class2 reveals context dependent regulatory roles for 5-hydroxymethylcytosine.

Since the discovery of 5-hydroxymethylcytosine (5hmC) as a prominent DNA modification found in mammalian genomes, an emergent question has been what role this mark plays in gene regulation. 5hmC is hypothesized to function as an intermediate in the demethylation of 5-methylcytosine (5mC) and in the reactivation of silenced promoters and enhancers. Further, weak positive correlations are observed between gene body 5hmC and gene expression. We previously demonstrated that ME-Class is an effective tool to understand relationships between whole-genome bisulfite sequencing data and expression. In this work, we present ME-Class2, a machine-learning based tool to perform integrative 5mCG, 5hmCG and expression analysis. Using ME-Class2 we analyze whole-genome single-base resolution 5mCG and 5hmCG datasets from 20 primary tissue and cell samples to reveal relationships between 5hmCG and expression. Our analysis indicates that conversion of 5mCG to 5hmCG within 2 kb of the transcription start site associates with distinct functions depending on the summed level of 5mCG + 5hmCG. Unchanged levels of 5mCG + 5hmCG (conversion from 5mCG to stable 5hmCG) associate with repression. Meanwhile, decreases in 5mCG + 5hmCG (5hmCG-mediated demethylation) associate with gene activation. Our results demonstrate that ME-Class2 will prove invaluable to interpret genome-wide 5mC and 5hmC datasets and guide mechanistic studies into the function of 5hmCG.

The Transcriptional Regulator SnoN Promotes the Proliferation of Cerebellar Granule Neuron Precursors in the Postnatal Mouse Brain.

Control of neuronal precursor cell proliferation is essential for normal brain development, and deregulation of this fundamental developmental event contributes to brain diseases. Typically, neuronal precursor cell proliferation extends over long periods of time during brain development. However, how neuronal precursor proliferation is regulated in a temporally specific manner remains to be elucidated. Here, we report that conditional KO of the transcriptional regulator SnoN in cerebellar granule neuron precursors robustly inhibits the proliferation of these cells and promotes their cell cycle exit at later stages of cerebellar development in the postnatal male and female mouse brain. In laser capture microdissection followed by RNA-Seq, designed to profile gene expression specifically in the external granule layer of the cerebellum, we find that SnoN promotes the expression of cell proliferation genes and concomitantly represses differentiation genes in granule neuron precursors in vivo Remarkably, bioinformatics analyses reveal that SnoN-regulated genes contain binding sites for the transcription factors N-myc and Pax6, which promote the proliferation and differentiation of granule neuron precursors, respectively. Accordingly, we uncover novel physical interactions of SnoN with N-myc and Pax6 in cells. In behavior analyses, conditional KO of SnoN impairs cerebellar-dependent learning in a delayed eye-blink conditioning paradigm, suggesting that SnoN-regulation of granule neuron precursor proliferation bears functional consequences at the organismal level. Our findings define a novel function and mechanism for the major transcriptional regulator SnoN in the control of granule neuron precursor proliferation in the mammalian brain.SIGNIFICANCE STATEMENT This study reports the discovery that the transcriptional regulator SnoN plays a crucial role in the proliferation of cerebellar granule neuron precursors in the postnatal mouse brain. Conditional KO of SnoN in granule neuron precursors robustly inhibits the proliferation of these cells and promotes their cycle exit specifically at later stages of cerebellar development, with biological consequences of impaired cerebellar-dependent learning. Genomics and bioinformatics analyses reveal that SnoN promotes the expression of cell proliferation genes and concomitantly represses cell differentiation genes in vivo Although SnoN has been implicated in distinct aspects of the development of postmitotic neurons, this study identifies a novel function for SnoN in neuronal precursors in the mammalian brain.

LONGO: an R package for interactive gene length dependent analysis for neuronal identity.

Motivation: Reprogramming somatic cells into neurons holds great promise to model neuronal development and disease. The efficiency and success rate of neuronal reprogramming, however, may vary between different conversion platforms and cell types, thereby necessitating an unbiased, systematic approach to estimate neuronal identity of converted cells. Recent studies have demonstrated that long genes (>100 kb from transcription start to end) are highly enriched in neurons, which provides an opportunity to identify neurons based on the expression of these long genes. Results: We have developed a versatile R package, LONGO, to analyze gene expression based on gene length. We propose a systematic analysis of long gene expression (LGE) with a metric termed the long gene quotient (LQ) that quantifies LGE in RNA-seq or microarray data to validate neuronal identity at the single-cell and population levels. This unique feature of neurons provides an opportunity to utilize measurements of LGE in transcriptome data to quickly and easily distinguish neurons from non-neuronal cells. By combining this conceptual advancement and statistical tool in a user-friendly and interactive software package, we intend to encourage and simplify further investigation into LGE, particularly as it applies to validating and improving neuronal differentiation and reprogramming methodologies. Availability and implementation: LONGO is freely available for download at https://github.com/biohpc/longo. Supplementary information: Supplementary data are available at Bioinformatics online.