Epigenetic modulation through BET bromodomain inhibitors as a novel therapeutic strategy for progranulin-deficient frontotemporal dementia

Frontotemporal dementia (FTD) is a debilitating neurodegenerative disorder with currently no disease-modifying treatment options available. Mutations in GRN are one of the most common genetic causes of FTD, near ubiquitously resulting in progranulin (PGRN) haploinsufficiency. Small molecules that can restore PGRN protein to healthy levels in individuals bearing a heterozygous GRN mutation may thus have therapeutic value. Here, we show that epigenetic modulation through bromodomain and extra-terminal domain (BET) inhibitors (BETi) potently enhance PGRN protein levels, both intracellularly and secreted forms, in human central nervous system (CNS)-relevant cell types, including in microglia-like cells. In terms of potential for disease modification, we show BETi treatment effectively restores PGRN levels in neural cells with a GRN mutation known to cause PGRN haploinsufficiency and FTD. We demonstrate that BETi can rapidly and durably enhance PGRN in neural progenitor cells (NPCs) in a manner dependent upon BET protein expression, suggesting a gain-of-function mechanism. We further describe a CNS-optimized BETi chemotype that potently engages endogenous BRD4 and enhances PGRN expression in neuronal cells. Our results reveal a new epigenetic target for treating PGRN-deficient forms of FTD and provide mechanistic insight to aid in translating this discovery into therapeutics.

Maternal dietary fat during lactation shapes single nucleus transcriptomic profile of postnatal offspring hypothalamus in a sexually dimorphic manner in mice

Maternal overnutrition during lactation predisposes offspring to develop metabolic diseases and exacerbates the relevant syndromes in males more than females in later life. The hypothalamus is a heterogenous brain region that regulates energy balance. Here we combined metabolic trait quantification of mother and offspring mice under low and high fat diet (HFD) feeding during lactation, with single nucleus transcriptomic profiling of their offspring hypothalamus at peak lacation to understand the cellular and molecular alterations in response to maternal dietary pertubation. We found significant expansion in neuronal subpopulations including histaminergic (Hdc), arginine vasopressin/retinoic acid receptor-related orphan receptor β (Avp/Rorb) and agouti-related peptide/neuropeptide Y (AgRP/Npy) in male offspring when their mothers were fed HFD, and increased Npy-astrocyte interactions in offspring responding to maternal overnutrition. Our study provides a comprehensive offspring hypothalamus map at the peak lactation and reveals how the cellular subpopulations respond to maternal dietary fat in a sex-specific manner during development.

Kdm1a safeguards the topological boundaries of PRC2-repressed genes and prevents aging-related euchromatinization in neurons

Kdm1a is a histone demethylase linked to intellectual disability with essential roles during gastrulation and the terminal differentiation of specialized cell types, including neurons, that remains highly expressed in the adult brain. To explore Kdm1a's function in adult neurons, we develop inducible and forebrain-restricted Kdm1a knockouts. By applying multi-omic transcriptome, epigenome and chromatin conformation data, combined with super-resolution microscopy, we find that Kdm1a elimination causes the neuronal activation of nonneuronal genes that are silenced by the polycomb repressor complex and interspersed with active genes. Functional assays demonstrate that the N-terminus of Kdm1a contains an intrinsically disordered region that is essential to segregate Kdm1a-repressed genes from the neighboring active chromatin environment. Finally, we show that the segregation of Kdm1a-target genes is weakened in neurons during natural aging, underscoring the role of Kdm1a safeguarding neuronal genome organization and gene silencing throughout life.

Integrative omics analysis reveals epigenomic and transcriptomic signatures underlying brain structural deficits in major depressive disorder

Several lines of evidence support the involvement of transcriptomic and epigenetic mechanisms in the brain structural deficits of major depressive disorder (MDD) separately. However, research in these two areas has remained isolated. In this study, we proposed an integrative strategy that combined neuroimaging, brain-wide gene expression, and peripheral DNA methylation data to investigate the genetic basis of gray matter abnormalities in MDD. The MRI T1-weighted images and Illumina 850 K DNA methylation microarrays were obtained from 269 patients and 416 healthy controls, and brain-wide transcriptomic data were collected from Allen Human Brain Atlas. The between-group differences in gray matter volume (GMV) and differentially methylated CpG positions (DMPs) were examined. The genes with their expression patterns spatially related to GMV changes and genes with DMPs were overlapped and selected. Using principal component regression, the associations between DMPs in overlapped genes and GMV across individual patients were investigated, and the region-specific correlations between methylation status and gene expression were examined. We found significant associations between the decreased GMV and DMPs methylation status in the anterior cingulate cortex, inferior frontal cortex, and fusiform face cortex regions. These DMPs genes were primarily enriched in the neurodevelopmental and synaptic transmission process. There was a significant negative correlation between DNA methylation and gene expression in genes associated with GMV changes of the frontal cortex in MDD. Our findings suggest that GMV abnormalities in MDD may have a transcriptomic and epigenetic basis. This imaging-transcriptomic-epigenetic integrative analysis provides spatial and biological links between cortical morphological deficits and peripheral epigenetic signatures in MDD.

Chronic oxytocin administration stimulates the oxytocinergic system in children with autism

Clinical efficacy of intranasal administration of oxytocin is increasingly explored in autism spectrum disorder, but to date, the biological effects of chronic administration regimes on endogenous oxytocinergic function are largely unknown. Here exploratory biological assessments from a completed randomized, placebo-controlled trial showed that children with autism (n = 79, 16 females) receiving intranasal oxytocin for four weeks (12 IU, twice daily) displayed significantly higher salivary oxytocin levels 24 hours after the last oxytocin nasal spray administration, but no longer at a four-week follow up session. Regarding salivary oxytocin receptor gene (OXTR) epigenetics (DNA-methylation), oxytocin-induced reductions in OXTR DNA-methylation were observed, suggesting a facilitation of oxytocin receptor expression in the oxytocin compared to the placebo group. Notably, heightened oxytocin levels post-treatment were significantly associated with reduced OXTR DNA-methylation and improved feelings of secure attachment. These findings indicate that four weeks of chronic oxytocin administration stimulated the endogenous oxytocinergic system in children with autism.

USP7/Maged1-mediated H2A monoubiquitination in the paraventricular thalamus: an epigenetic mechanism involved in cocaine use disorder

The risk of developing drug addiction is strongly influenced by the epigenetic landscape and chromatin remodeling. While histone modifications such as methylation and acetylation have been studied in the ventral tegmental area and nucleus accumbens (NAc), the role of H2A monoubiquitination remains unknown. Our investigations, initially focused on the scaffold protein melanoma-associated antigen D1 (Maged1), reveal that H2A monoubiquitination in the paraventricular thalamus (PVT) significantly contributes to cocaine-adaptive behaviors and transcriptional repression induced by cocaine. Chronic cocaine use increases H2A monoubiquitination, regulated by Maged1 and its partner USP7. Accordingly, Maged1 specific inactivation in thalamic Vglut2 neurons, or USP7 inhibition, blocks cocaine-evoked H2A monoubiquitination and cocaine locomotor sensitization. Additionally, genetic variations in MAGED1 and USP7 are linked to altered susceptibility to cocaine addiction and cocaine-associated symptoms in humans. These findings unveil an epigenetic modification in a non-canonical reward pathway of the brain and a potent marker of epigenetic risk factors for drug addiction in humans.

Identification of candidate DNA methylation biomarkers related to Alzheimer's disease risk by integrating genome and blood methylome data

Alzheimer disease (AD) is a common neurodegenerative disease with a late onset. It is critical to identify novel blood-based DNA methylation biomarkers to better understand the extent of the molecular pathways affected in AD. Two sets of blood DNA methylation genetic prediction models developed using different reference panels and modelling strategies were leveraged to evaluate associations of genetically predicted DNA methylation levels with AD risk in 111,326 (46,828 proxy) cases and 677,663 controls. A total of 1,168 cytosine-phosphate-guanine (CpG) sites showed a significant association with AD risk at a false discovery rate (FDR) < 0.05. Methylation levels of 196 CpG sites were correlated with expression levels of 130 adjacent genes in blood. Overall, 52 CpG sites of 32 genes showed consistent association directions for the methylation-gene expression-AD risk, including nine genes (CNIH4, THUMPD3, SERPINB9, MTUS1, CISD1, FRAT2, CCDC88B, FES, and SSH2) firstly reported as AD risk genes. Nine of 32 genes were enriched in dementia and AD disease categories (P values ranged from 1.85 × 10-4 to 7.46 × 10-6), and 19 genes in a neurological disease network (score = 54) were also observed. Our findings improve the understanding of genetics and etiology for AD.

Sperm chromatin accessibility's involvement in the intergenerational effects of stress hormone receptor activation

Dexamethasone is a stress hormone receptor agonist used widely in clinics. We and others previously showed that paternal administration of dexamethasone in mice affects the phenotype of their offspring. The substrate of intergenerational transmission of environmentally induced effects often involves changes in sperm RNA, yet other epigenetic modifications in the germline can be affected and are also plausible candidates. First, we tested the involvement of altered sperm RNAs in the transmission of dexamethasone induced phenotypes across generations. We did this by injecting sperm RNA into naïve fertilized oocytes, before performing metabolic and behavioral phenotyping of the offspring. We observed phenotypic changes in discordance with those found in offspring generated by in vitro fertilization using sperm from dexamethasone exposed males. Second, we investigated the effect of dexamethasone on chromatin accessibility using ATAC sequencing and found significant changes at specific genomic features and gene regulatory loci. Employing q-RT-PCR, we show altered expression of a gene in the tissue of offspring affected by accessibility changes in sperm. Third, we establish a correlation between specific DNA modifications and stress hormone receptor activity as a likely contributing factor influencing sperm accessibility. Finally, we independently investigated this dependency by genetically reducing thymine-DNA glycosylase levels and observing concomitant changes at the level of chromatin accessibility and stress hormone receptor activity.

CBLL1 is hypomethylated and correlates with cortical thickness in transgender men before gender affirming hormone treatment

Gender identity refers to the consciousness of being a man, a woman or other condition. Although it is generally congruent with the sex assigned at birth, for some people it is not. If the incongruity is distressing, it is defined as gender dysphoria (GD). Here, we measured whole-genome DNA methylation by the Illumina © Infinium Human Methylation 850k array and reported its correlation with cortical thickness (CTh) in 22 transgender men (TM) experiencing GD versus 25 cisgender men (CM) and 28 cisgender women (CW). With respect to the methylation analysis, TM vs. CW showed significant differences in 35 CpGs, while 2155 CpGs were found when TM vs. CM were compared. With respect to correlation analysis, TM showed differences in methylation of CBLL1 and DLG1 genes that correlated with global and left hemisphere CTh. Both genes were hypomethylated in TM compared to the cisgender groups. Early onset TM showed a positive correlation between CBLL1 and several cortical regions in the frontal (left caudal middle frontal), temporal (right inferior temporal, left fusiform) and parietal cortices (left supramarginal and right paracentral). This is the first study relating CBLL1 methylation with CTh in transgender persons and supports a neurodevelopmental hypothesis of gender identity.

Brain-wide correspondence of neuronal epigenomics and distant projections

Single-cell analyses parse the brain's billions of neurons into thousands of 'cell-type' clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Single-cell DNA methylome and 3D multi-omic atlas of the adult mouse brain

Cytosine DNA methylation is essential in brain development and is implicated in various neurological disorders. Understanding DNA methylation diversity across the entire brain in a spatial context is fundamental for a complete molecular atlas of brain cell types and their gene regulatory landscapes. Here we used single-nucleus methylome sequencing (snmC-seq3) and multi-omic sequencing (snm3C-seq)1 technologies to generate 301,626 methylomes and 176,003 chromatin conformation-methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell taxonomy with 4,673 cell groups and 274 cross-modality-annotated subclasses. We identified 2.6 million differentially methylated regions across the genome that represent potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide spatial transcriptomics data validated the association of spatial epigenetic diversity with transcription and improved the anatomical mapping of our epigenetic datasets. Furthermore, chromatin conformation diversities occurred in important neuronal genes and were highly associated with DNA methylation and transcription changes. Brain-wide cell-type comparisons enabled the construction of regulatory networks that incorporate transcription factors, regulatory elements and their potential downstream gene targets. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a whole-brain SMART-seq2 dataset. Our study establishes a brain-wide, single-cell DNA methylome and 3D multi-omic atlas and provides a valuable resource for comprehending the cellular-spatial and regulatory genome diversity of the mouse brain.

Divergent single cell transcriptome and epigenome alterations in ALS and FTD patients with C9orf72 mutation

A repeat expansion in the C9orf72 (C9) gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here we investigate single nucleus transcriptomics (snRNA-seq) and epigenomics (snATAC-seq) in postmortem motor and frontal cortices from C9-ALS, C9-FTD, and control donors. C9-ALS donors present pervasive alterations of gene expression with concordant changes in chromatin accessibility and histone modifications. The greatest alterations occur in upper and deep layer excitatory neurons, as well as in astrocytes. In neurons, the changes imply an increase in proteostasis, metabolism, and protein expression pathways, alongside a decrease in neuronal function. In astrocytes, the alterations suggest activation and structural remodeling. Conversely, C9-FTD donors have fewer high-quality neuronal nuclei in the frontal cortex and numerous gene expression changes in glial cells. These findings highlight a context-dependent molecular disruption in C9-ALS and C9-FTD, indicating unique effects across cell types, brain regions, and diseases.

Epigenetic regulation in major depression and other stress-related disorders: molecular mechanisms, clinical relevance and therapeutic potential

Major depressive disorder (MDD) is a chronic, generally episodic and debilitating disease that affects an estimated 300 million people worldwide, but its pathogenesis is poorly understood. The heritability estimate of MDD is 30-40%, suggesting that genetics alone do not account for most of the risk of major depression. Another factor known to associate with MDD involves environmental stressors such as childhood adversity and recent life stress. Recent studies have emerged to show that the biological impact of environmental factors in MDD and other stress-related disorders is mediated by a variety of epigenetic modifications. These epigenetic modification alterations contribute to abnormal neuroendocrine responses, neuroplasticity impairment, neurotransmission and neuroglia dysfunction, which are involved in the pathophysiology of MDD. Furthermore, epigenetic marks have been associated with the diagnosis and treatment of MDD. The evaluation of epigenetic modifications holds promise for further understanding of the heterogeneous etiology and complex phenotypes of MDD, and may identify new therapeutic targets. Here, we review preclinical and clinical epigenetic findings, including DNA methylation, histone modification, noncoding RNA, RNA modification, and chromatin remodeling factor in MDD. In addition, we elaborate on the contribution of these epigenetic mechanisms to the pathological trait variability in depression and discuss how such mechanisms can be exploited for therapeutic purposes.

Inheritance of associative memories and acquired cellular changes in C. elegans

Experiences have been shown to modulate behavior and physiology of future generations in some contexts, but there is limited evidence for inheritance of associative memory in different species. Here, we trained C. elegans nematodes to associate an attractive odorant with stressful starvation conditions and revealed that this associative memory was transmitted to the F1 progeny who showed odor-evoked avoidance behavior. Moreover, the F1 and the F2 descendants of trained animals exhibited odor-evoked cellular stress responses, manifested by the translocation of DAF-16/FOXO to cells' nuclei. Sperm, but not oocytes, transmitted these odor-evoked cellular stress responses which involved H3K9 and H3K36 methylations, the small RNA pathway machinery, and intact neuropeptide secretion. Activation of a single chemosensory neuron sufficed to induce a serotonin-mediated systemic stress response in both the parental trained generation and in its progeny. Moreover, inheritance of the cellular stress responses increased survival chances of the progeny as exposure to the training odorant allowed the animals to prepare in advance for an impending adversity. These findings suggest that in C. elegans associative memories and cellular changes may be transferred across generations.

Mir125b-2 imprinted in human but not mouse brain regulates hippocampal function and circuit in mice

Genomic imprinting predominantly occurs in the placenta and brain. Few imprinted microRNAs have been identified in the brain, and their functional roles in the brain are not clear. Here we show paternal, but not maternal, expression of MIR125B2 in human but not mouse brain. Moreover, Mir125b-2^m-/p- mice showed impaired learning and memory, and anxiety, whose functions were hippocampus-dependent. Hippocampal granule cells from Mir125b-2^m-/p- mice displayed increased neuronal excitability, increased excitatory synaptic transmission, and decreased inhibitory synaptic transmission. Glutamate ionotropic receptor NMDA type subunit 2A (Grin2a), a key regulator of synaptic plasticity, was physically bound by miR-125b-2 and upregulated in the hippocampus of Mir125b-2^m-/p- mice. Taken together, our findings demonstrate MIR125B2 imprinted in human but not mouse brain, mediated learning, memory, and anxiety, regulated excitability and synaptic transmission in hippocampal granule cells, and affected hippocampal expression of Grin2a. Our work provides functional mechanisms of a species-specific imprinted microRNA in the brain.

Noncanonical regulation of imprinted gene Igf2 by amyloid-beta 1-42 in Alzheimer's disease

Reduced insulin-like growth factor 2 (IGF2) levels in Alzheimer's disease (AD) may be the mechanism relating age-related metabolic disorders to dementia. Since Igf2 is an imprinted gene, we examined age and sex differences in the relationship between amyloid-beta 1-42 (Aβ₄₂) accumulation and epigenetic regulation of the Igf2/H19 gene cluster in cerebrum, liver, and plasma of young and old male and female 5xFAD mice, in frontal cortex of male and female AD and non-AD patients, and in HEK293 cell cultures. We show IGF2 levels, Igf2 expression, histone acetylation, and H19 ICR methylation are lower in females than males. However, elevated Aβ₄₂ levels are associated with Aβ₄₂ binding to Igf2 DMR2, increased DNA and histone methylation, and a reduction in Igf2 expression and IGF2 levels in 5xFAD mice and AD patients, independent of H19 ICR methylation. Cell culture results confirmed the binding of Aβ₄₂ to Igf2 DMR2 increased DNA and histone methylation, and reduced Igf2 expression. These results indicate an age- and sex-related causal relationship among Aβ₄₂ levels, epigenomic state, and Igf2 expression in AD and provide a potential mechanism for Igf2 regulation in normal and pathological conditions, suggesting IGF2 levels may be a useful diagnostic biomarker for Aβ₄₂ targeted AD therapies.

The autism risk factor CHD8 is a chromatin activator in human neurons and functionally dependent on the ERK-MAPK pathway effector ELK1

The chromodomain helicase DNA-binding protein CHD8 is the most frequently mutated gene in autism spectrum disorder. Despite its prominent disease involvement, little is known about its molecular function in the human brain. CHD8 is a chromatin regulator which binds to the promoters of actively transcribed genes through genomic targeting mechanisms which have yet to be fully defined. By generating a conditional loss-of-function and an endogenously tagged allele in human pluripotent stem cells, we investigated the molecular function and the interaction of CHD8 with chromatin in human neurons. Chromatin accessibility analysis and transcriptional profiling revealed that CHD8 functions as a transcriptional activator at its target genes in human neurons. Furthermore, we found that CHD8 chromatin targeting is cell context-dependent. In human neurons, CHD8 preferentially binds at ETS motif-enriched promoters. This enrichment is particularly prominent on the promoters of genes whose expression significantly changes upon the loss of CHD8. Indeed, among the ETS transcription factors, we identified ELK1 as being most highly correlated with CHD8 expression in primary human fetal and adult cortical neurons and most highly expressed in our stem cell-derived neurons. Remarkably, ELK1 was necessary to recruit CHD8 specifically to ETS motif-containing sites. These findings imply that ELK1 and CHD8 functionally cooperate to regulate gene expression and chromatin states at MAPK/ERK target genes in human neurons. Our results suggest that the MAPK/ERK/ELK1 axis potentially contributes to the pathogenesis caused by CHD8 mutations in human neurodevelopmental disorders.

Cell-type specific profiling of histone post-translational modifications in the adult mouse striatum

Epigenetic gene regulation in the heterogeneous brain remains challenging to decipher with current strategies. Bulk tissue analysis from pooled subjects reflects the average of cell-type specific changes across cell-types and individuals, which obscures causal relationships between epigenetic modifications, regulation of gene expression, and complex pathology. To address these limitations, we optimized a hybrid protocol, ICuRuS, for the isolation of nuclei tagged in specific cell-types and histone post translational modification profiling from the striatum of a single mouse. We combined affinity-based isolation of the medium spiny neuron subtypes, Adenosine 2a Receptor or Dopamine Receptor D1, with cleavage of histone-DNA complexes using an antibody-targeted micrococcal nuclease to release DNA complexes for paired end sequencing. Unlike fluorescence activated cell sorting paired with chromatin immunoprecipitation, ICuRuS allowed for robust epigenetic profiling at cell-type specific resolution. Our analysis provides a framework to understand combinatorial relationships between neuronal-subtype-specific epigenetic modifications and gene expression.

High-throughput robust single-cell DNA methylation profiling with sciMETv2

DNA methylation is a key epigenetic property that drives gene regulatory programs in development and disease. Current single-cell methods that produce high quality methylomes are expensive and low throughput without the aid of extensive automation. We previously described a proof-of-principle technique that enabled high cell throughput; however, it produced only low-coverage profiles and was a difficult protocol that required custom sequencing primers and recipes and frequently produced libraries with excessive adapter contamination. Here, we describe a greatly improved version that generates high-coverage profiles (~15-fold increase) using a robust protocol that does not require custom sequencing capabilities, includes multiple stopping points, and exhibits minimal adapter contamination. We demonstrate two versions of sciMETv2 on primary human cortex, a high coverage and rapid version, identifying distinct cell types using CH methylation patterns. These datasets are able to be directly integrated with one another as well as with existing snmC-seq2 datasets with little discernible bias. Finally, we demonstrate the ability to determine cell types using CG methylation alone, which is the dominant context for DNA methylation in most cell types other than neurons and the most applicable analysis outside of brain tissue.

Methylation pattern and mRNA expression of synapse-relevant genes in the MAM model of schizophrenia in the time-course of adolescence

Schizophrenia is highly heritable and aggregating in families, but genetics alone does not exclusively explain the pathogenesis. Many risk factors, including childhood trauma, viral infections, migration, and the use of cannabis, are associated with schizophrenia. Adolescence seems to be the critical period where symptoms of the disease manifest. This work focuses on studying an epigenetic regulatory mechanism (the role of DNA methylation) and its interaction with mRNA expression during development, with a particular emphasis on adolescence. The presumptions regarding the role of aberrant neurodevelopment in schizophrenia were tested in the Methyl-Azoxy-Methanol (MAM) animal model. MAM treatment induces neurodevelopmental disruptions and behavioral deficits in off-springs of the treated animals reminiscent of those observed in schizophrenia and is thus considered a promising model for studying this pathology. On a gestational day-17, adult pregnant rats were treated with the antimitotic agent MAM. Experimental animals were divided into groups and subgroups according to substance treatment (MAM and vehicle agent [Sham]) and age of analysis (pre-adolescent and post-adolescent). Methylation and mRNA expression analysis of four candidate genes, which are often implicated in schizophrenia, with special emphasis on the Dopamine hypothesis i.e., Dopamine receptor D₂ (Drd2), and the "co-factors" Disrupted in schizophrenia 1 (DISC1), Synaptophysin (Syp), and Dystrobrevin-binding protein 1 (Dtnbp1), was performed in the Gyrus cingulum (CING) and prefrontal cortex (PFC). Data were analyzed to observe the effect of substance treatment between groups and the impact of adolescence within-group. We found reduced pre-adolescent expression levels of Drd2 in both brain areas under the application of MAM. The "co-factor genes" did not show high deviations in mRNA expression levels but high alterations of methylation rates under the application of MAM (up to ~20%), which diminished in the further time course, reaching a comparable level like in Sham control animals after adolescence. The pre-adolescent reduction in DRD2 expression might be interpreted as downregulation of the receptor due to hyperdopaminergic signaling from the ventral tegmental area (VTA), eventually even to both investigated brain regions. The notable alterations of methylation rates in the three analyzed co-factor genes might be interpreted as attempt to compensate for the altered dopaminergic neurotransmission.

DNA methylation and hydroxymethylation characterize the identity of D1 and D2 striatal projection neurons

Neuronal DNA modifications differ from those in other cells, including methylation outside CpG context and abundant 5-hydroxymethylation whose relevance for neuronal identities are unclear. Striatal projection neurons expressing D1 or D2 dopamine receptors allow addressing this question, as they share many characteristics but differ in their gene expression profiles, connections, and functional roles. We compare translating mRNAs and DNA modifications in these two populations. DNA methylation differences occur predominantly in large genomic clusters including differentially expressed genes, potentially important for D1 and D2 neurons. Decreased gene body methylation is associated with higher gene expression. Hydroxymethylation differences are more scattered and affect transcription factor binding sites, which can influence gene expression. We also find a strong genome-wide hydroxymethylation asymmetry between the two DNA strands, particularly pronounced at expressed genes and retrotransposons. These results identify novel properties of neuronal DNA modifications and unveil epigenetic characteristics of striatal projection neurons heterogeneity.

Dementia with Lewy bodies post-mortem brains reveal differentially methylated CpG sites with biomarker potential

Dementia with Lewy bodies (DLB) is a common form of dementia with known genetic and environmental interactions. However, the underlying epigenetic mechanisms which reflect these gene-environment interactions are poorly studied. Herein, we measure genome-wide DNA methylation profiles of post-mortem brain tissue (Broadmann area 7) from 15 pathologically confirmed DLB brains and compare them with 16 cognitively normal controls using Illumina MethylationEPIC arrays. We identify 17 significantly differentially methylated CpGs (DMCs) and 17 differentially methylated regions (DMRs) between the groups. The DMCs are mainly located at the CpG islands, promoter and first exon regions. Genes associated with the DMCs are linked to "Parkinson's disease" and "metabolic pathway", as well as the diseases of "severe intellectual disability" and "mood disorders". Overall, our study highlights previously unreported DMCs offering insights into DLB pathogenesis with the possibility that some of these could be used as biomarkers of DLB in the future.

Epigenome-wide DNA methylation analysis of whole blood cells derived from patients with GAD and OCD in the Chinese Han population

Generalized anxiety disorder (GAD) and obsessive-compulsive disorder (OCD) had high comorbidity and affected more than 44 million people around the world leading to a huge burden on health and economy. Here, we conducted an epigenome-wide DNA methylation study employing 93 patients with GAD, 65 patients with OCD, and 302 health controls, to explore epigenetic alterations associated with the onset and differences of GAD and OCD. We identified multiple differentially methylated positions (DMPs) and regions (DMRs): three DMP genes included RIOK3 (cg21515243, p = 8.00 × 10^-10), DNASE2 (cg09379601, p = 1.10 × 10^-9), and PSMB4 (cg01334186, p = 3.70 × 10^-7) and two DMR genes USP6NL (p = 4.50 × 10^-4) and CPLX1 (p = 6.95 × 10^-4) were associated with the onset of GAD and OCD; three DMPs genes included LDLRAP1 (cg21400344, p = 4.40 × 10^-12), ACIN1 (cg23712970, p = 2.98×10^-11), and SCRT1 (cg25472897, p = 5.60 × 10^-11) and three DMR genes WDR19 (p = 3.39 × 10^-3), SYCP1 (p = 6.41 × 10^-3), and FAM172A (p = 5.74 × 10^-3) were associated with the differences between GAD and OCD. Investigation of epigenetic age and chronological age revealed a different epigenetic development trajectory of GAD and OCD. Conclusively, our findings which yielded robust models may aid in distinguishing patients from healthy controls (AUC = 0.90-0.99) or classifying patients with GAD and OCD (AUC = 0.89-0.99), and may power the precision medicine for them.

Meta-analysis of epigenome-wide association studies of major depressive disorder

Epigenetic mechanisms have been hypothesized to play a role in the etiology of major depressive disorder (MDD). In this study, we performed a meta-analysis between two case-control MDD cohorts to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) in MDD. Using samples from two Cohorts (a total of 298 MDD cases and 63 controls with repeated samples, on average ~ 1.8 samples/subject), we performed an EWAS meta-analysis. Multiple cytosine-phosphate-guanine sites annotated to TNNT3 were associated with MDD reaching study-wide significance, including cg08337959 (p = 2.3 × 10^-11). Among DMPs with association p values less than 0.0001, pathways from REACTOME such as Ras activation upon Ca²⁺ influx through the NMDA receptor (p = 0.0001, p-adjusted = 0.05) and long-term potentiation (p = 0.0002, p-adjusted = 0.05) were enriched in this study. A total of 127 DMRs with Sidak-corrected p value < 0.05 were identified from the meta-analysis, including DMRs annotated to TNNT3 (chr11: 1948933 to 1949130 [6 probes], Sidak corrected P value = 4.32 × 10^-41), S100A13 (chr1: 153599479 to 153600972 [22 probes], Sidak corrected P value = 5.32 × 10^-18), NRXN1 (chr2: 50201413 to 50201505 [4 probes], Sidak corrected P value = 1.19 × 10^-11), IL17RA (chr22: 17564750 to 17565149, Sidak corrected P value = 9.31 × 10^-8), and NPFFR2 (chr4: 72897565 to 72898212, Sidak corrected P value = 8.19 × 10^-7). Using 2 Cohorts of depression case-control samples, we identified DMPs and DMRs associated with MDD. The molecular pathways implicated by these data include mechanisms involved in neuronal synaptic plasticity, calcium signaling, and inflammation, consistent with reports from previous genetic and protein biomarker studies indicating that these mechanisms are involved in the neurobiology of depression.

Two kinds of transcription factors mediate chronic morphine-induced decrease in miR-105 in medial prefrontal cortex of rats

Chronic morphine administration alters gene expression in different brain regions, an effect which may contribute to plastic changes associated with addictive behavior. This change in gene expression is most possibly mediated by addictive drug-induced epigenetic remodeling of gene expression programs. Our previous studies showed that chronic morphine-induced decrease of miR-105 in the medial prefrontal cortex (mPFC) contributed to context-induced retrieval of morphine withdrawal memory. However, how chronic morphine treatment decreases miR-105 in the mPFC still remains unknown. The present study shows that chronic morphine induces addiction-related change in miR-105 in the mPFC via two kinds of transcription factors: the first transcription factor is CREB activated by mu receptors-ERK-p90RSK signaling pathway and the second transcription factor is glucocorticoid receptor (GR), which as a negative transcription factor, mediates chronic morphine-induced decrease in miR-105 in the mPFC of rats.

NRSF/REST lies at the intersection between epigenetic regulation, miRNA-mediated gene control and neurodevelopmental pathways associated with Intellectual disability (ID) and Schizophrenia

Genetic evidence indicates disrupted epigenetic regulation as a major risk factor for psychiatric disorders, but the molecular mechanisms that drive this association remain to be determined. EHMT1 is an epigenetic repressor that is causal for Kleefstra Syndrome (KS), a genetic disorder linked with neurodevelopmental disorders and associated with schizophrenia. Here, we show that reduced EHMT1 activity decreases NRSF/REST protein leading to abnormal neuronal gene expression and progression of neurodevelopment in human iPSC. We further show that EHMT1 regulates NRSF/REST indirectly via repression of miRNA and leads to aberrant neuronal gene regulation and neurodevelopment timing. Expression of a NRSF/REST mRNA that lacks the miRNA-binding sites restores neuronal gene regulation to EHMT1 deficient cells. Significantly, the EHMT1-regulated miRNA gene set not only controls NRSF/REST but is enriched for association for Intellectual Disability (ID) and schizophrenia. This reveals a broad molecular interaction between H3K9 demethylation, NSRF/REST regulation and risk for ID and Schizophrenia.

Dimensions of childhood adversity differentially affect biological aging in major depression

Adverse childhood experiences have been consistently linked with physical and mental health disorders in adulthood that may be mediated, in part, via the effects of such exposures on biological aging. Using recently developed "epigenetic clocks", which provide an estimate of biological age, several studies have demonstrated a link between the cumulative exposure to childhood adversities and accelerated epigenetic aging. However, not all childhood adversities are equivalent and less is known about how distinct dimensions of childhood adversity relate to epigenetic aging metrics. Using two measures of childhood adversity exposure, we assess how the dimensions of Maltreatment and Household Dysfunction relate to epigenetic aging using two "second-generation" clocks, GrimAge and PhenoAge, in a cohort of unmedicated somatically healthy adults with moderate to severe major depression (n = 82). Our results demonstrate that the dimension of Maltreatment is associated with epigenetic age acceleration (EAA) using the PhenoAge but not the GrimAge clock. This association was observed using both the Childhood Trauma questionnaire (CTQ; β = 0.272, p = 0.013) and the Adverse Childhood Experiences (ACEs) questionnaire (β = 0.307, p = 0.005) and remained significant when adjusting for exposure to the dimension of Household Dysfunction (β = 0.322, p = 0.009). In contrast, the dimension of Household Dysfunction is associated with epigenetic age deceleration (β = -0.194, p = 0.083) which achieved significance after adjusting for exposure to the dimension of Maltreatment (β = -0.304, p = 0.022). This study is the first to investigate these effects among individuals with Major Depressive Disorder and suggests that these dimensions of adversity may be associated with disease via distinct biological mechanisms.

Cocaine regulation of Nr4a1 chromatin bivalency and mRNA in male and female mice

Cocaine epigenetically regulates gene expression via changes in histone post-translational modifications (HPTMs). We previously found that the immediate early gene Nr4a1 is epigenetically activated by cocaine in mouse brain reward regions. However, few studies have examined multiple HPTMs at a single gene. Bivalent gene promoters are simultaneously enriched in both activating (H3K4me3 (K4)) and repressive (H3K27me3 (K27)) HPTMs. As such, bivalent genes are lowly expressed but poised for activity-dependent gene regulation. In this study, we identified K4&K27 bivalency at Nr4a1 following investigator-administered cocaine in male and female mice. We applied sequential chromatin immunoprecipitation and qPCR to define Nr4a1 bivalency and expression in striatum (STR), prefrontal cortex (PFC), and hippocampus (HPC). We used Pearson's correlation to quantify relationships within each brain region across treatment conditions for each sex. In female STR, cocaine increased Nr4a1 mRNA while maintaining Nr4a1 K4&K27 bivalency. In male STR, cocaine enriched repressive H3K27me3 and K4&K27 bivalency at Nr4a1 and maintained Nr4a1 mRNA. Furthermore, cocaine epigenetically regulated a putative NR4A1 target, Cartpt, in male PFC. This study defined the epigenetic regulation of Nr4a1 in reward brain regions in male and female mice following cocaine, and, thus, shed light on the biological relevance of sex to cocaine use disorder.

Regulation of human cortical interneuron development by the chromatin remodeling protein CHD2

Mutations in the chromodomain helicase DNA binding protein 2 (CHD2) gene are associated with neurodevelopmental disorders. However, mechanisms by which CHD2 regulates human brain development remain largely uncharacterized. Here, we used a human embryonic stem cell model of cortical interneuron (hcIN) development to elucidate its roles in this process. We identified genome-wide CHD2 binding profiles during hcIN differentiation, defining direct CHD2 targets related to neurogenesis in hcIN progenitors and to neuronal function in hcINs. CHD2 bound sites were frequently coenriched with histone H3 lysine 27 acetylation (H3K27ac) and associated with high gene expression, indicating roles for CHD2 in promoting gene expression during hcIN development. Binding sites for different classes of transcription factors were enriched at CHD2 bound regions during differentiation, suggesting transcription factors that may cooperatively regulate stage-specific gene expression with CHD2. We also demonstrated that CHD2 haploinsufficiency altered CHD2 and H3K27ac coenrichment on chromatin and expression of associated genes, decreasing acetylation and expression of cell cycle genes while increasing acetylation and expression of neuronal genes, to cause precocious differentiation. Together, these data describe CHD2 direct targets and mechanisms by which CHD2 prevents precocious hcIN differentiation, which are likely to be disrupted by pathogenic CHD2 mutation to cause neurodevelopmental disorders.

Spatial profiling of chromatin accessibility in mouse and human tissues

Cellular function in tissue is dependent on the local environment, requiring new methods for spatial mapping of biomolecules and cells in the tissue context¹. The emergence of spatial transcriptomics has enabled genome-scale gene expression mapping^2-5, but the ability to capture spatial epigenetic information of tissue at the cellular level and genome scale is lacking. Here we describe a method for spatially resolved chromatin accessibility profiling of tissue sections using next-generation sequencing (spatial-ATAC-seq) by combining in situ Tn5 transposition chemistry⁶ and microfluidic deterministic barcoding⁵. Profiling mouse embryos using spatial-ATAC-seq delineated tissue-region-specific epigenetic landscapes and identified gene regulators involved in the development of the central nervous system. Mapping the accessible genome in the mouse and human brain revealed the intricate arealization of brain regions. Applying spatial-ATAC-seq to tonsil tissue resolved the spatially distinct organization of immune cell types and states in lymphoid follicles and extrafollicular zones. This technology progresses spatial biology by enabling spatially resolved chromatin accessibility profiling to improve our understanding of cell identity, cell state and cell fate decision in relation to epigenetic underpinnings in development and disease.

Differential microRNA expression analyses across two brain regions in Alzheimer's disease

Dysregulation of microRNAs (miRNAs) is involved in the pathogenesis of neurodegenerative diseases, including Alzheimer's disease (AD). Hitherto, sample sizes from differential miRNA expression studies in AD are exceedingly small aggravating any biological inference. To overcome this limitation, we investigated six candidate miRNAs in a large collection of brain samples. Brain tissue was derived from superior temporal gyrus (STG) and entorhinal cortex (EC) from 99 AD patients and 91 controls. MiRNA expression was examined by qPCR (STG) or small RNA sequencing (EC). Brain region-dependent differential miRNA expression was investigated in a transgenic AD mouse model using qPCR and FISH. Total RNA sequencing was used to assess differential expression of miRNA target genes. MiR-129-5p, miR-132-5p, and miR-138-5p were significantly downregulated in AD vs. controls both in STG and EC, while miR-125b-5p and miR-501-3p showed no evidence for differential expression in this dataset. In addition, miR-195-5p was significantly upregulated in EC but not STG in AD patients. The brain region-specific pattern of miR-195-5p expression was corroborated in vivo in transgenic AD mice. Total RNA sequencing identified several novel and functionally interesting target genes of these miRNAs involved in synaptic transmission (GABRB1), the immune-system response (HCFC2) or AD-associated differential methylation (SLC16A3). Using two different methods (qPCR and small RNA-seq) in two separate brain regions in 190 individuals we more than doubled the available sample size for most miRNAs tested. Differential gene expression analyses confirm the likely involvement of miR-129-5p, miR-132-5p, miR-138-5p, and miR-195-5p in AD pathogenesis and highlight several novel potentially relevant target mRNAs.

Cell-type-specific epigenetic effects of early life stress on the brain

Early life stress (ELS) induces long-term phenotypic adaptations that contribute to increased vulnerability to a host of neuropsychiatric disorders. Epigenetic mechanisms, including DNA methylation, histone modifications and non-coding RNA, are a proposed link between environmental stressors, alterations in gene expression, and phenotypes. Epigenetic modifications play a primary role in shaping functional differences between cell types and can be modified by environmental perturbations, especially in early development. Together with contributions from genetic variation, epigenetic mechanisms orchestrate patterns of gene expression within specific cell types that contribute to phenotypic variation between individuals. To date, many studies have provided insights into epigenetic changes resulting from ELS. However, most of these studies have examined heterogenous brain tissue, despite evidence of cell-type-specific epigenetic modifications in phenotypes associated with ELS. In this review, we focus on rodent and human studies that have examined epigenetic modifications induced by ELS in select cell types isolated from the brain or associated with genes that have cell-type-restricted expression in neurons, microglia, astrocytes, and oligodendrocytes. Although significant challenges remain, future studies using these approaches can enable important mechanistic insight into the role of epigenetic variation in the effects of ELS on brain function.

Perinatal exposure to fluoxetine and maternal adversity affect myelin-related gene expression and epigenetic regulation in the corticolimbic circuit of juvenile rats

Many pregnant women experience symptoms of depression, and are often treated with selective serotonin reuptake inhibitor (SSRI) antidepressants, such as fluoxetine. In utero exposure to SSRIs and maternal depressive symptoms is associated with sex-specific effects on the brain and behavior. However, knowledge about the neurobiological mechanisms underlying these sex differences is limited. In addition, most animal research into developmental SSRI exposure neglects the influence of maternal adversity. Therefore, we used a rat model relevant to depression to investigate the molecular effects of perinatal fluoxetine exposure in male and female juvenile offspring. We performed RNA sequencing and targeted DNA methylation analyses on the prefrontal cortex and basolateral amygdala; key regions of the corticolimbic circuit. Perinatal fluoxetine enhanced myelin-related gene expression in the prefrontal cortex, while inhibiting it in the basolateral amygdala. SSRI exposure and maternal adversity interacted to affect expression of genes such as myelin-associated glycoprotein (Mag) and myelin basic protein (Mbp). We speculate that altered myelination reflects altered brain maturation. In addition, these effects are stronger in males than in females, resembling known behavioral outcomes. Finally, Mag and Mbp expression correlated with DNA methylation, highlighting epigenetic regulation as a potential mechanism for developmental fluoxetine-induced changes in myelination.

Early life adversity shapes neural circuit function during sensitive postnatal developmental periods

Early life adversity (ELA) is a major risk factor for mental illness, but the neurobiological mechanisms by which ELA increases the risk for future psychopathology are still poorly understood. Brain development is particularly malleable during prenatal and early postnatal life, when complex neural circuits are being formed and refined through an interplay of excitatory and inhibitory neural input, synaptogenesis, synaptic pruning, myelination, and neurogenesis. Adversity that influences these processes during sensitive periods of development can thus have long-lasting and pervasive effects on neural circuit maturation. In this review, we will discuss clinical and preclinical evidence for the impact of ELA on neural circuit formation with a focus on the early postnatal period, and how long-lasting impairments in these circuits can affect future behavior. We provide converging evidence from human and animal studies on how ELA alters the functional development of brain regions, neural circuits, and neurotransmitter systems that are crucial for cognition and affective behavior, including the hippocampus, the hypothalamus-pituitary-adrenal (HPA) axis, neural networks of fear responses and cognition, and the serotonin (5-HT) system. We also discuss how gene-by-environment (GxE) interactions can determine individual differences in susceptibility and resilience to ELA, as well as molecular pathways by which ELA regulates neural circuit development, for which we emphasize epigenetic mechanisms. Understanding the molecular and neurobiological mechanisms underlying ELA effects on brain function and psychopathology during early postnatal sensitive periods may have great potential to advance strategies to better treat or prevent psychiatric disorders that have their origin early in life.

Accelerated epigenetic aging in suicide attempters uninfluenced by high intent-to-die and choice of lethal methods

Suicide attempts (SA) are associated with excess non-suicidal mortality, putatively mediated in part by premature cellular senescence. Epigenetic age (EA) estimators of biological age have been previously demonstrated to strongly predict physiological dysregulation and mortality risk. Herein, we investigate if violent SA with high intent-to-die is predictive of epigenetics-derived estimates of biological aging. The genome-wide methylation pattern was measured using the Illumina Infinium Methylation EPIC BeadChip in whole blood of 88 suicide attempters. Subjects were stratified into two groups based on the putative risk of later committed suicide (low- [n = 58] and high-risk [n = 30]) in dependency of SA method (violent or non-violent) and/or intent-to-die (high/low). Estimators of intrinsic and extrinsic EA acceleration, one marker optimized to predict physiological dysregulation (DNAmPhenoAge/AgeAccelPheno) and one optimized to predict lifespan (DNAmGrimAge/AgeAccelGrim) were investigated for associations to severity of SA, by univariate and multivariate analyses. The study was adequately powered to detect differences of 2.2 years in AgeAccelGrim in relation to SA severity. Baseline DNAmGrimAge exceeded chronological age by 7.3 years on average across all samples, conferring a mean 24.6% increase in relation to actual age. No individual EA acceleration marker was differentiated by suicidal risk group (p > 0.1). Thus, SA per se but not severity of SA is related to EA, implicating that excess non-suicidal mortality in SA is unrelated to risk of committed suicide. Preventative healthcare efforts aimed at curtailing excess mortality after SA may benefit from acting equally powerful to recognize somatic comorbidities irrespective of the severity inherent in the act itself.

Gene regulation by gonadal hormone receptors underlies brain sex differences

Oestradiol establishes neural sex differences in many vertebrates1-3 and modulates mood, behaviour and energy balance in adulthood4-8. In the canonical pathway, oestradiol exerts its effects through the transcription factor oestrogen receptor-α (ERα)9. Although ERα has been extensively characterized in breast cancer, the neuronal targets of ERα, and their involvement in brain sex differences, remain largely unknown. Here we generate a comprehensive map of genomic ERα-binding sites in a sexually dimorphic neural circuit that mediates social behaviours. We conclude that ERα orchestrates sexual differentiation of the mouse brain through two mechanisms: establishing two male-biased neuron types and activating a sustained male-biased gene expression program. Collectively, our findings reveal that sex differences in gene expression are defined by hormonal activation of neuronal steroid receptors. The molecular targets we identify may underlie the effects of oestradiol on brain development, behaviour and disease.

Epigenome-wide DNA methylation in obsessive-compulsive disorder

In adult patients with obsessive-compulsive disorder (OCD), altered DNA methylation has been discerned in several candidate genes, while DNA methylation on an epigenome-wide level has been investigated in only one Chinese study so far. Here, an epigenome-wide association study (EWAS) was performed in a sample of 76 OCD patients of European ancestry (37 women, age ± SD: 33.51 ± 10.92 years) and 76 sex- and age-matched healthy controls for the first time using the Illumina MethylationEPIC BeadChip. After quality control, nine epigenome-wide significant quantitative trait methylation sites (QTMs) and 21 suggestive hits were discerned in the final sample of 68 patients and 68 controls. The top hit (cg24159721) and four other significant QTMs (cg11894324, cg01070250, cg11330075, cg15174812) map to the region of the microRNA 12136 gene (MIR12136). Two additional significant CpG sites (cg05740793, cg20450977) are located in the flanking region of the MT-RNR2 (humanin) like 8 gene (MT-RNRL8), while two further QTMs (cg16267121, cg15890734) map to the regions of the MT-RNR2 (humanin) like 3 (MT-RNRL3) and MT-RNR2 (humanin) like 2 (MT-RNRL2) genes. Provided replication of the present findings in larger samples, the identified QTMs might provide more biological insight into the pathogenesis of OCD and thereby could in the future serve as peripheral epigenetic markers of OCD risk with the potential to inform targeted preventive and therapeutic efforts.

An epigenetic association analysis of childhood trauma in psychosis reveals possible overlap with methylation changes associated with PTSD

Patients with a severe mental disorder report significantly higher levels of childhood trauma (CT) than healthy individuals. Studies have suggested that CT may affect brain plasticity through epigenetic mechanisms and contribute to developing various psychiatric disorders. We performed a blood-based epigenome-wide association study using the Childhood Trauma Questionnaire-short form in 602 patients with a current severe mental illness, investigating DNA methylation association separately for five trauma subtypes and the total trauma score. The median trauma score was set as the predefined cutoff for determining whether the trauma was present or not. Additionally, we compared our genome-wide results with methylation probes annotated to candidate genes previously associated with CT. Of the patients, 83.2% reported CT above the cutoff in one or more trauma subtypes, and emotional neglect was the trauma subtype most frequently reported. We identified one significant differently methylated position associated with the gene TANGO6 for physical neglect. Seventeen differentially methylated regions (DMRs) were associated with different trauma categories. Several of these DMRs were annotated to genes previously associated with neuropsychiatric disorders such as post-traumatic stress disorder and cognitive impairments. Our results support a biomolecular association between CT and severe mental disorders. Genes that were previously identified as differentially methylated in CT-exposed subjects with and without psychosis did not show methylation differences in our analysis. We discuss this inconsistency, the relevance of our findings, and the limitations of our study.

Brain microRNAs are associated with variation in cognitive trajectory in advanced age

In advancing age, some individuals maintain a stable cognitive performance over time, while others experience a rapid decline. Such variation in cognitive trajectory is only partially explained by common neurodegenerative pathologies. Hence, we aimed to identify new molecular processes underlying variation in cognitive trajectory using brain microRNA profile followed by an integrative analysis with brain transcriptome and proteome. Individual cognitive trajectories were derived from longitudinally assessed cognitive-test scores of older-adult brain donors from four longitudinal cohorts. Postmortem brain microRNA profiles, transcriptomes, and proteomes were derived from the dorsolateral prefrontal cortex. The global microRNA association study of cognitive trajectory was performed in a discovery (n = 454) and replication cohort (n = 134), followed by a meta-analysis that identified 6 microRNAs. Among these, miR-132-3p and miR-29a-3p were most significantly associated with cognitive trajectory. They explain 18.2% and 2.0% of the variance of cognitive trajectory, respectively, and act independently of the eight measured neurodegenerative pathologies. Furthermore, integrative transcriptomic and proteomic analyses revealed that miR-132-3p was significantly associated with 24 of the 47 modules of co-expressed genes of the transcriptome, miR-29a-3p with 3 modules, and identified 84 and 214 downstream targets of miR-132-3p and miR-29a-3p, respectively, in cognitive trajectory. This is the first global microRNA study of cognitive trajectory to our knowledge. We identified miR-29a-3p and miR-132-3p as novel and robust contributors to cognitive trajectory independently of the eight known cerebral pathologies. Our findings lay a foundation for future studies investigating mechanisms and developing interventions to enhance cognitive stability in advanced age.

Identification of a modular super-enhancer in murine retinal development

Super-enhancers are expansive regions of genomic DNA comprised of multiple putative enhancers that contribute to the dynamic gene expression patterns during development. This is particularly important in neurogenesis because many essential transcription factors have complex developmental stage- and cell-type specific expression patterns across the central nervous system. In the developing retina, Vsx2 is expressed in retinal progenitor cells and is maintained in differentiated bipolar neurons and Müller glia. A single super-enhancer controls this complex and dynamic pattern of expression. Here we show that deletion of one region disrupts retinal progenitor cell proliferation but does not affect cell fate specification. The deletion of another region has no effect on retinal progenitor cell proliferation but instead leads to a complete loss of bipolar neurons. This prototypical super-enhancer may serve as a model for dissecting the complex gene expression patterns for neurogenic transcription factors during development. Moreover, it provides a unique opportunity to alter expression of individual transcription factors in particular cell types at specific stages of development. This provides a deeper understanding of function that cannot be achieved with traditional knockout mouse approaches.

Corticosterone-mediated regulation and functions of miR-218-5p in rat brain

Chronic stress is one of the key precipitating factors in major depressive disorder (MDD). Stress associated studies have underscored the mechanistic role of epigenetic master players like microRNAs (miRNAs) in depression pathophysiology at both preclinical and clinical levels. Previously, we had reported changes in miR-218-5p expression in response to corticosterone (CORT) induced chronic stress. MiR-218-5p was one of the most significantly induced miRNAs in the prefrontal cortex (PFC) of rats under chronic stress. In the present report, we have investigated how chronic CORT exposure mechanistically affected miR-218-5p expression in the rat brain and how miR-218 could trigger molecular changes on its downstream regulatory pathways. Elevated expression of miR-218-5p was found in the PFC of CORT-treated rats. A glucocorticoid receptor (GR) targeted Chromatin-Immunoprecipitation (ChIP) assay revealed high GR occupancy on the promoter region of Slit3 gene hosting miR-218-2 in its 3rd intron. RNA-sequencing data based on RNA Induced silencing Complex Immunoprecipitation (RISC-IP) with AGO2 in SH-SY5Y cells detected six consistent target genes of miR-218-5p (APOL4, DTWD1, BNIP1, METTL22, SNAPC1, and HDAC6). The expression of all five genes, except APOL4, was successfully validated with qPCR in CORT-treated rat PFC. Further, Hdac6-based ChIP-seq experiment helped in mapping major genomic loci enriched for intergenic regions in the PFC of CORT-treated rat. A proximity-based gene ontology (GO) analysis revealed a majority of the intergenic sites to be part of key genes implicated in central nervous system functions, notably synapse organization, neuron projection morphogenesis, and axonogenesis. Our results suggest that the upregulation of miR-218-5p in PFC of CORT-treated rats possibly resulted from GR biding in the promoter region of Slit3 gene. Interestingly, Hdac6 was one of the consistent target genes potentially found to regulate CNS related genes by chromatin modification. Collectively, these findings establish the role of miR-218-5p in chronic stress and the epigenetic function it plays to induce chromatin-based transcriptional changes of several CNS genes in triggering stress-induced disorders, including depression. This also opens up the scope to understand the role of miR-218-5p as a potential target for noncoding RNA therapeutics in clinical depression.

Differential brain ADRA2A and ADRA2C gene expression and epigenetic regulation in schizophrenia. Effect of antipsychotic drug treatment

Postsynaptic α_2A-adrenoceptor density is enhanced in the dorsolateral prefrontal cortex (DLPFC) of antipsychotic-treated schizophrenia subjects. This alteration might be due to transcriptional activation, and could be regulated by epigenetic mechanisms such as histone posttranslational modifications (PTMs). The aim of this study was to evaluate ADRA2A and ADRA2C gene expression (codifying for α₂-adrenoceptor subtypes), and permissive and repressive histone PTMs at gene promoter regions in the DLPFC of subjects with schizophrenia and matched controls (n = 24 pairs). We studied the effect of antipsychotic (AP) treatment in AP-free (n = 12) and AP-treated (n = 12) subgroups of schizophrenia subjects and in rats acutely and chronically treated with typical and atypical antipsychotics. ADRA2A mRNA expression was selectively upregulated in AP-treated schizophrenia subjects (+93%) whereas ADRA2C mRNA expression was upregulated in all schizophrenia subjects (+53%) regardless of antipsychotic treatment. Acute and chronic clozapine treatment in rats did not alter brain cortex Adra2a mRNA expression but increased Adra2c mRNA expression. Both ADRA2A and ADRA2C promoter regions showed epigenetic modification by histone methylation and acetylation in human DLPFC. The upregulation of ADRA2A expression in AP-treated schizophrenia subjects might be related to observed bivalent chromatin at ADRA2A promoter region in schizophrenia (depicted by increased permissive H3K4me3 and repressive H3K27me3) and could be triggered by the enhanced H4K16ac at ADRA2A promoter. In conclusion, epigenetic predisposition differentially modulated ADRA2A and ADRA2C mRNA expression in DLPFC of schizophrenia subjects.

Early-life nicotine or cotinine exposure produces long-lasting sleep alterations and downregulation of hippocampal corticosteroid receptors in adult mice

Early-life exposure to environmental toxins like tobacco can permanently re-program body structure and function. Here, we investigated the long-term effects on mouse adult sleep phenotype exerted by early-life exposure to nicotine or to its principal metabolite, cotinine. Moreover, we investigated whether these effects occurred together with a reprogramming of the activity of the hippocampus, a key structure to coordinate the hormonal stress response. Adult male mice born from dams subjected to nicotine (NIC), cotinine (COT) or vehicle (CTRL) treatment in drinking water were implanted with electrodes for sleep recordings. NIC and COT mice spent significantly more time awake than CTRL mice at the transition between the rest (light) and the activity (dark) period. NIC and COT mice showed hippocampal glucocorticoid receptor (GR) downregulation compared to CTRL mice, and NIC mice also showed hippocampal mineralocorticoid receptor downregulation. Hippocampal GR expression significantly and inversely correlated with the amount of wakefulness at the light-to-dark transition, while no changes in DNA methylation were found. We demonstrated that early-life exposure to nicotine (and cotinine) concomitantly entails long-lasting reprogramming of hippocampal activity and sleep phenotype suggesting that the adult sleep phenotype may be modulated by events that occurred during that critical period of life.

An atlas of gene regulatory elements in adult mouse cerebrum

The mammalian cerebrum performs high-level sensory perception, motor control and cognitive functions through highly specialized cortical and subcortical structures1. Recent surveys of mouse and human brains with single-cell transcriptomics2-6 and high-throughput imaging technologies7,8 have uncovered hundreds of neural cell types distributed in different brain regions, but the transcriptional regulatory programs that are responsible for the unique identity and function of each cell type remain unknown. Here we probe the accessible chromatin in more than 800,000 individual nuclei from 45 regions that span the adult mouse isocortex, olfactory bulb, hippocampus and cerebral nuclei, and use the resulting data to map the state of 491,818 candidate cis-regulatory DNA elements in 160 distinct cell types. We find high specificity of spatial distribution for not only excitatory neurons, but also most classes of inhibitory neurons and a subset of glial cell types. We characterize the gene regulatory sequences associated with the regional specificity within these cell types. We further link a considerable fraction of the cis-regulatory elements to putative target genes expressed in diverse cerebral cell types and predict transcriptional regulators that are involved in a broad spectrum of molecular and cellular pathways in different neuronal and glial cell populations. Our results provide a foundation for comprehensive analysis of gene regulatory programs of the mammalian brain and assist in the interpretation of noncoding risk variants associated with various neurological diseases and traits in humans.

DNA methylation atlas of the mouse brain at single-cell resolution

Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing1,2 to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data3 enabled prediction of high-confidence enhancer-gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments4. By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum.

Epigenomic diversity of cortical projection neurons in the mouse brain

Neuronal cell types are classically defined by their molecular properties, anatomy and functions. Although recent advances in single-cell genomics have led to high-resolution molecular characterization of cell type diversity in the brain1, neuronal cell types are often studied out of the context of their anatomical properties. To improve our understanding of the relationship between molecular and anatomical features that define cortical neurons, here we combined retrograde labelling with single-nucleus DNA methylation sequencing to link neural epigenomic properties to projections. We examined 11,827 single neocortical neurons from 63 cortico-cortical and cortico-subcortical long-distance projections. Our results showed unique epigenetic signatures of projection neurons that correspond to their laminar and regional location and projection patterns. On the basis of their epigenomes, intra-telencephalic cells that project to different cortical targets could be further distinguished, and some layer 5 neurons that project to extra-telencephalic targets (L5 ET) formed separate clusters that aligned with their axonal projections. Such separation varied between cortical areas, which suggests that there are area-specific differences in L5 ET subtypes, which were further validated by anatomical studies. Notably, a population of cortico-cortical projection neurons clustered with L5 ET rather than intra-telencephalic neurons, which suggests that a population of L5 ET cortical neurons projects to both targets. We verified the existence of these neurons by dual retrograde labelling and anterograde tracing of cortico-cortical projection neurons, which revealed axon terminals in extra-telencephalic targets including the thalamus, superior colliculus and pons. These findings highlight the power of single-cell epigenomic approaches to connect the molecular properties of neurons with their anatomical and projection properties.