Alterations in DNA 5-hydroxymethylation patterns in the hippocampus of an experimental model of chronic epilepsy.

Temporal lobe epilepsy (TLE) is a type of focal epilepsy characterized by spontaneous recurrent seizures originating from the hippocampus. The epigenetic reprogramming hypothesis of epileptogenesis suggests that the development of TLE is associated with alterations in gene transcription changes resulting in a hyperexcitable network in TLE. DNA 5-methylcytosine (5-mC) is an epigenetic mechanism that has been associated with chronic epilepsy. However, the contribution of 5-hydroxymethylcytosine (5-hmC), a product of 5-mC demethylation by the Ten-Eleven Translocation (TET) family proteins in chronic TLE is poorly understood. 5-hmC is abundant in the brain and acts as a stable epigenetic mark altering gene expression through several mechanisms. Here, we found that the levels of bulk DNA 5-hmC but not 5-mC were significantly reduced in the hippocampus of human TLE patients and in the kainic acid (KA) TLE rat model. Using 5-hmC hMeDIP-sequencing, we characterized 5-hmC distribution across the genome and found bidirectional regulation of 5-hmC at intergenic regions within gene bodies. We found that hypohydroxymethylated 5-hmC intergenic regions were associated with several epilepsy-related genes, including Gal, SV2, and Kcnj11 and hyperdroxymethylation 5-hmC intergenic regions were associated with Gad65, TLR4, and Bdnf gene expression. Mechanistically, Tet1 knockdown in the hippocampus was sufficient to decrease 5-hmC levels and increase seizure susceptibility following KA administration. In contrast, Tet1 overexpression in the hippocampus resulted in increased 5-hmC levels associated with improved seizure resiliency in response to KA. These findings suggest an important role for 5-hmC as an epigenetic regulator of epilepsy that can be manipulated to influence seizure outcomes.

Distinct roles of Bdnf I and Bdnf IV transcript variant expression in hippocampal neurons.

Brain-derived neurotrophic factor (Bdnf) plays a critical role in brain development, dendritic growth, synaptic plasticity, as well as learning and memory. The rodent Bdnf gene contains nine 5' non-coding exons (I-IXa), which are spliced to a common 3' coding exon (IX). Transcription of individual Bdnf variants, which all encode the same BDNF protein, is initiated at unique promoters upstream of each non-coding exon, enabling precise spatiotemporal and activity-dependent regulation of Bdnf expression. Although prior evidence suggests that Bdnf transcripts containing exon I (Bdnf I) or exon IV (Bdnf IV) are uniquely regulated by neuronal activity, the functional significance of different Bdnf transcript variants remains unclear. To investigate functional roles of activity-dependent Bdnf I and IV transcripts, we used a CRISPR activation system in which catalytically dead Cas9 fused to a transcriptional activator (VPR) is targeted to individual Bdnf promoters with single guide RNAs, resulting in transcript-specific Bdnf upregulation. Bdnf I upregulation is associated with gene expression changes linked to dendritic growth, while Bdnf IV upregulation is associated with genes that regulate protein catabolism. Upregulation of Bdnf I, but not Bdnf IV, increased mushroom spine density, volume, length, and head diameter, and also produced more complex dendritic arbors in cultured rat hippocampal neurons. In contrast, upregulation of Bdnf IV, but not Bdnf I, in the rat hippocampus attenuated contextual fear expression. Our data suggest that while Bdnf I and IV are both activity-dependent, BDNF produced from these promoters may serve unique cellular, synaptic, and behavioral functions.

Ten simple rules for using public biological data for your research.

With an increasing amount of biological data available publicly, there is a need for a guide on how to successfully download and use this data. The 10 simple rules for using public biological data are: (1) use public data purposefully in your research; (2) evaluate data for your use case; (3) check data reuse requirements and embargoes; (4) be aware of ethics for data reuse; (5) plan for data storage and compute requirements; (6) know what you are downloading; (7) download programmatically and verify integrity; (8) properly cite data; (9) make reprocessed data and models Findable, Accessible, Interoperable, and Reusable (FAIR) and share; and (10) make pipelines and code FAIR and share. These rules are intended as a guide for researchers wanting to make use of available data and to increase data reuse and reproducibility.

Distinct subpopulations of D1 medium spiny neurons exhibit unique transcriptional responsiveness to cocaine.

Drugs of abuse increase extracellular concentrations of dopamine in the nucleus accumbens (NAc), resulting in transcriptional alterations that drive long-lasting cellular and behavioral adaptations. While decades of research have focused on the transcriptional mechanisms by which drugs of abuse influence neuronal physiology and function, few studies have comprehensively defined NAc cell type heterogeneity in transcriptional responses to drugs of abuse. Here, we used single nucleus RNA-seq (snRNA-seq) to characterize the transcriptome of over 39,000 NAc cells from male and female adult Sprague-Dawley rats following acute or repeated cocaine experience. This dataset identified 16 transcriptionally distinct cell populations, including two populations of medium spiny neurons (MSNs) that express the Drd1 dopamine receptor (D1-MSNs). Critically, while both populations expressed classic marker genes of D1-MSNs, only one population exhibited a robust transcriptional response to cocaine. Validation of population-selective transcripts using RNA in situ hybridization revealed distinct spatial compartmentalization of these D1-MSN populations within the NAc. Finally, analysis of published NAc snRNA-seq datasets from non-human primates and humans demonstrated conservation of MSN subtypes across rat and higher order mammals, and further highlighted cell type-specific transcriptional differences across the NAc and broader striatum. These results highlight the utility in using snRNA-seq to characterize both cell type heterogeneity and cell type-specific responses to cocaine and provides a useful resource for cross-species comparisons of NAc cell composition.

DOCK3 is a dosage-sensitive regulator of skeletal muscle and Duchenne muscular dystrophy-associated pathologies.

DOCK3 is a member of the DOCK family of guanine nucleotide exchange factors that regulate cell migration, fusion and viability. Previously, we identified a dysregulated miR-486/DOCK3 signaling cascade in dystrophin-deficient muscle, which resulted in the overexpression of DOCK3; however, little is known about the role of DOCK3 in muscle. Here, we characterize the functional role of DOCK3 in normal and dystrophic skeletal muscle. Utilizing Dock3 global knockout (Dock3 KO) mice, we found that the haploinsufficiency of Dock3 in Duchenne muscular dystrophy mice improved dystrophic muscle pathologies; however, complete loss of Dock3 worsened muscle function. Adult Dock3 KO mice have impaired muscle function and Dock3 KO myoblasts are defective for myogenic differentiation. Transcriptomic analyses of Dock3 KO muscles reveal a decrease in myogenic factors and pathways involved in muscle differentiation. These studies identify DOCK3 as a novel modulator of muscle health and may yield therapeutic targets for treating dystrophic muscle symptoms.

Enhancer RNAs predict enhancer-gene regulatory links and are critical for enhancer function in neuronal systems.

Genomic enhancer elements regulate gene expression programs important for neuronal fate and function and are implicated in brain disease states. Enhancers undergo bidirectional transcription to generate non-coding enhancer RNAs (eRNAs). However, eRNA function remains controversial. Here, we combined Assay for Transposase-Accessible Chromatin using Sequencing (ATAC-Seq) and RNA-Seq datasets from three distinct neuronal culture systems in two activity states, enabling genome-wide enhancer identification and prediction of putative enhancer-gene pairs based on correlation of transcriptional output. Notably, stimulus-dependent enhancer transcription preceded mRNA induction, and CRISPR-based activation of eRNA synthesis increased mRNA at paired genes, functionally validating enhancer-gene predictions. Focusing on enhancers surrounding the Fos gene, we report that targeted eRNA manipulation bidirectionally modulates Fos mRNA, and that Fos eRNAs directly interact with the histone acetyltransferase domain of the enhancer-linked transcriptional co-activator CREB-binding protein (CBP). Together, these results highlight the unique role of eRNAs in neuronal gene regulation and demonstrate that eRNAs can be used to identify putative target genes.

Human and rodent temporal lobe epilepsy is characterized by changes in O-GlcNAc homeostasis that can be reversed to dampen epileptiform activity.

Temporal Lobe Epilepsy (TLE) is frequently associated with changes in protein composition and post-translational modifications (PTM) that exacerbate the disorder. O-linked-ß-N-acetyl glucosamine (O-GlcNAc) is a PTM occurring at serine/threonine residues that is derived from and closely associated with metabolic substrates. The enzymes O-GlcNActransferase (OGT) and O-GlcNAcase (OGA) mediate the addition and removal, respectively, of the O-GlcNAc modification. The goal of this study was to characterize OGT/OGA and protein O-GlcNAcylation in the epileptic hippocampus and to determine and whether direct manipulation of these proteins and PTM's alter epileptiform activity. We observed reduced global and protein specific O-GlcNAcylation and OGT expression in the kainate rat model of TLE and in human TLE hippocampal tissue. Inhibiting OGA with Thiamet-G elevated protein O-GlcNAcylation, and decreased both seizure duration and epileptic spike events, suggesting that OGA may be a therapeutic target for seizure control. These findings suggest that loss of O-GlcNAc homeostasis in the kainate model and in human TLE can be reversed via targeting of O-GlcNAc related pathways.

Autosomal dominant retinitis pigmentosa rhodopsin mutant Q344X drives specific alterations in chromatin complex gene transcription.

Purpose: Epigenetic and transcriptional mechanisms have been shown to contribute to long-lasting functional changes in adult neurons. The purpose of this study was to identify any such modifications in diseased retinal tissues from a mouse model of rhodopsin mutation-associated autosomal dominant retinitis pigmentosa (ADRP), Q344X, relative to age-matched wild-type (WT) controls. Methods: We performed RNA sequencing (RNA-seq) at poly(A) selected RNA to profile the transcriptional patterns in 3-week-old ADRP mouse model rhodopsin Q344X compared to WT controls. Differentially expressed genes were determined by DESeq2 using the Benjamini & Hochberg p value adjustment and an absolute log2 fold change cutoff. Quantitative western blots were conducted to evaluate protein expression levels of histone H3 phosphorylated at serine 10 and histone H4. qRT-PCR was performed to validate the expression patterns of differentially expressed genes. Results: We observed significant differential expression in 2151 genes in the retina of Q344X mice compared to WT controls, including downregulation in the potassium channel gene, Kcnv2, and differential expression of histone genes, including the H1 family histone member, H1foo; the H3 histone family 3B, H3f3b; and the histone deacetylase 9, Hdac9. Quantitative western blots revealed statistically significant decreased protein expression of both histone H3 phosphorylated at serine 10 and histone H4 in 3-week-old Q344X retinas. Furthermore, qRT-PCR performed on select differentially expressed genes based on our RNA-seq results revealed matched expression patterns of up or downregulation. Conclusions: These findings provide evidence that transcriptomic alterations occur in the ADRP mouse model rhodopsin Q344X retina and that these processes may contribute to the dysfunction and neurodegeneration seen in this animal model.

Acute ocular hypertension in the living human eye: Model description and initial cellular responses to elevated intraocular pressure.

This initial methods study presents the initial immunohistochemical and transcriptomic changes in the optic nerve head and retina from three research-consented brain-dead organ donors following prolonged and transient intraocular pressure (IOP) elevation. In this initial study, research-consented brain-dead organ donors were exposed to unilateral elevation of IOP for 7.5 h (Donor 1), 30 h (Donor 2), and 1 h (Donor 3) prior to organ procurement. Optic nerve tissue and retinal tissue was obtained following organ procurement for immunohistological and transcriptomic analysis. Optic nerve sections in Donor 1 exposed to 7.5-hours of unilateral sub-ischemic IOP elevation demonstrated higher levels of protein expression of the astrocytic marker, glial fibrillary acidic protein (GFAP), within the lamina cribrosa with greatest expression inferior temporally in the treated eye compared to control. Spatial transcriptomic analysis performed on optic nerve head tissues from Donor 2 exposed to 30 h of unilateral IOP elevation demonstrated differential transcription of mRNA across laminar and scleral regions. Immunohistochemistry of retinal sections from Donor 2 exhibited higher GFAP and IBA1 expression in the treated eye compared with control, but this was not observed in Donor 3, which was exposed to only 1-hour of IOP elevation. While there were no differences in GFAP protein expression in the retina following the 1-hour IOP elevation in Donor 3, there were higher levels of transcription of GFAP in the inner nuclear layer, and CD44 in the retinal ganglion cell layer, indicative of astrocytic and Müller glial reactivity as well as an early inflammatory response, respectively. We found that transcriptomic differences can be observed across treated and control eyes following unilateral elevation of IOP in brain dead organ donors. The continued development of this model affords the unique opportunity to define the acute mechanotranscriptomic response of the optic nerve head, evaluate the injury and repair mechanisms in the retina in response to IOP elevation, and enable correlation of in vivo imaging and functional testing with ex vivo cellular responses for the first time in the living human eye.

Reelin marks cocaine-activated striatal ensembles, promotes neuronal excitability, and regulates cocaine reward.

Drugs of abuse activate defined neuronal ensembles in brain reward structures such as the nucleus accumbens (NAc), which are thought to promote the enduring synaptic, circuit, and behavioral consequences of drug exposure. While the molecular and cellular effects arising from experience with drugs like cocaine are increasingly well understood, the mechanisms that sculpt NAc ensemble participation are largely unknown. Here, we leveraged unbiased single-nucleus transcriptional profiling to identify expression of the secreted glycoprotein Reelin (encoded by the Reln gene) as a marker of cocaine-activated neuronal ensembles within the rat NAc. Multiplexed in situ detection confirmed selective expression of the immediate early gene Fos in Reln+ neurons after cocaine experience, and also revealed enrichment of Reln mRNA in Drd1 + medium spiny neurons (MSNs) in both the rat and human brain. Using a novel CRISPR interference strategy enabling selective Reln knockdown in the adult NAc, we observed altered expression of genes linked to calcium signaling, emergence of a transcriptional trajectory consistent with loss of cocaine sensitivity, and a striking decrease in MSN intrinsic excitability. At the behavioral level, loss of Reln prevented cocaine locomotor sensitization, abolished cocaine place preference memory, and decreased cocaine self-administration behavior. Together, these results identify Reelin as a critical mechanistic link between ensemble participation and cocaine-induced behavioral adaptations.

Temporally specific gene expression and chromatin remodeling programs regulate a conserved Pdyn enhancer.

Neuronal and behavioral adaptations to novel stimuli are regulated by temporally dynamic waves of transcriptional activity, which shape neuronal function and guide enduring plasticity. Neuronal activation promotes expression of an immediate early gene (IEG) program comprised primarily of activity-dependent transcription factors, which are thought to regulate a second set of late response genes (LRGs). However, while the mechanisms governing IEG activation have been well studied, the molecular interplay between IEGs and LRGs remain poorly characterized. Here, we used transcriptomic and chromatin accessibility profiling to define activity-driven responses in rat striatal neurons. As expected, neuronal depolarization generated robust changes in gene expression, with early changes (1 hr) enriched for inducible transcription factors and later changes (4 hr) enriched for neuropeptides, synaptic proteins, and ion channels. Remarkably, while depolarization did not induce chromatin remodeling after 1 hr, we found broad increases in chromatin accessibility at thousands of sites in the genome at 4 hr after neuronal stimulation. These putative regulatory elements were found almost exclusively at non-coding regions of the genome, and harbored consensus motifs for numerous activity-dependent transcription factors such as AP-1. Furthermore, blocking protein synthesis prevented activity-dependent chromatin remodeling, suggesting that IEG proteins are required for this process. Targeted analysis of LRG loci identified a putative enhancer upstream of Pdyn (prodynorphin), a gene encoding an opioid neuropeptide implicated in motivated behavior and neuropsychiatric disease states. CRISPR-based functional assays demonstrated that this enhancer is both necessary and sufficient for Pdyn transcription. This regulatory element is also conserved at the human PDYN locus, where its activation is sufficient to drive PDYN transcription in human cells. These results suggest that IEGs participate in chromatin remodeling at enhancers and identify a conserved enhancer that may act as a therapeutic target for brain disorders involving dysregulation of Pdyn.

Distinct subpopulations of D1 medium spiny neurons exhibit unique transcriptional responsiveness to cocaine.

Drugs of abuse increase extracellular concentrations of dopamine in the nucleus accumbens (NAc), resulting in transcriptional alterations that drive long-lasting cellular and behavioral adaptations. While decades of research have focused on the transcriptional mechanisms by which drugs of abuse influence neuronal physiology and function, few studies have comprehensively defined NAc cell type heterogeneity in transcriptional responses to drugs of abuse. Here, we used single nucleus RNA-seq (snRNA-seq) to characterize the transcriptome of over 39,000 NAc cells from male and female adult Sprague-Dawley rats following acute or repeated cocaine experience. This dataset identified 16 transcriptionally distinct cell populations, including two populations of medium spiny neurons (MSNs) that express the Drd1 dopamine receptor (D1-MSNs). Critically, while both populations expressed classic marker genes of D1-MSNs, only one population exhibited a robust transcriptional response to cocaine. Validation of population-selective transcripts using RNA in situ hybridization revealed distinct spatial compartmentalization of these D1-MSN populations within the NAc. Finally, analysis of published NAc snRNA-seq datasets from non-human primates and humans demonstrated conservation of MSN subtypes across rat and higher order mammals, and further highlighted cell type-specific transcriptional differences across the NAc and broader striatum. These results highlight the utility in using snRNA-seq to characterize both cell type heterogeneity and cell type-specific responses to cocaine and provides a useful resource for cross-species comparisons of NAc cell composition.

Temporally specific gene expression and chromatin remodeling programs regulate a conserved Pdyn enhancer.

Neuronal and behavioral adaptations to novel stimuli are regulated by temporally dynamic waves of transcriptional activity, which shape neuronal function and guide enduring plasticity. Neuronal activation promotes expression of an immediate early gene (IEG) program comprised primarily of activity-dependent transcription factors, which are thought to regulate a second set of late response genes (LRGs). However, while the mechanisms governing IEG activation have been well studied, the molecular interplay between IEGs and LRGs remain poorly characterized. Here, we used transcriptomic and chromatin accessibility profiling to define activity-driven responses in rat striatal neurons. As expected, neuronal depolarization generated robust changes in gene expression, with early changes (1 h) enriched for inducible transcription factors and later changes (4 h) enriched for neuropeptides, synaptic proteins, and ion channels. Remarkably, while depolarization did not induce chromatin remodeling after 1 h, we found broad increases in chromatin accessibility at thousands of sites in the genome at 4 h after neuronal stimulation. These putative regulatory elements were found almost exclusively at non-coding regions of the genome, and harbored consensus motifs for numerous activity-dependent transcription factors such as AP-1. Furthermore, blocking protein synthesis prevented activity-dependent chromatin remodeling, suggesting that IEG proteins are required for this process. Targeted analysis of LRG loci identified a putative enhancer upstream of Pdyn (prodynorphin), a gene encoding an opioid neuropeptide implicated in motivated behavior and neuropsychiatric disease states. CRISPR-based functional assays demonstrated that this enhancer is both necessary and sufficient for Pdyn transcription. This regulatory element is also conserved at the human PDYN locus, where its activation is sufficient to drive PDYN transcription in human cells. These results suggest that IEGs participate in chromatin remodeling at enhancers and identify a conserved enhancer that may act as a therapeutic target for brain disorders involving dysregulation of Pdyn.

Alterations in DNA 5-hydroxymethylation Patterns in the Hippocampus of an Experimental Model of Refractory Epilepsy.

Temporal lobe epilepsy (TLE) is a type of focal epilepsy characterized by spontaneous recurrent seizures originating from the hippocampus. The epigenetic reprogramming hypothesis of epileptogenesis suggests that the development of TLE is associated with alterations in gene transcription changes resulting in a hyperexcitable network in TLE. DNA 5-methylcytosine (5-mC) is an epigenetic mechanism that has been associated with chronic epilepsy. However, the contribution of 5-hydroxymethylcytosine (5-hmC), a product of 5-mC demethylation by the Ten-Eleven Translocation (TET) family proteins in chronic TLE is poorly understood. 5-hmC is abundant in the brain and acts as a stable epigenetic mark altering gene expression through several mechanisms. Here, we found that the levels of bulk DNA 5-hmC but not 5-mC were significantly reduced in the hippocampus of human TLE patients and in the kainic acid (KA) TLE rat model. Using 5-hmC hMeDIP-sequencing, we characterized 5-hmC distribution across the genome and found bidirectional regulation of 5-hmC at intergenic regions within gene bodies. We found that hypohydroxymethylated 5-hmC intergenic regions were associated with several epilepsy-related genes, including Gal , SV2, and Kcnj11 and hyperdroxymethylation 5-hmC intergenic regions were associated with Gad65 , TLR4 , and Bdnf gene expression. Mechanistically, Tet1 knockdown in the hippocampus was sufficient to decrease 5-hmC levels and increase seizure susceptibility following KA administration. In contrast, Tet1 overexpression in the hippocampus resulted in increased 5-hmC levels associated with improved seizure resiliency in response to KA. These findings suggest an important role for 5-hmC as an epigenetic regulator of epilepsy that can be manipulated to influence seizure outcomes.

Distinct roles of Bdnf I and Bdnf IV transcript variant expression in hippocampal neurons.

Brain-derived neurotrophic factor (Bdnf) plays a critical role in brain development, dendritic growth, synaptic plasticity, as well as learning and memory. The rodent Bdnf gene contains nine 5' non-coding exons (I-IXa), which are spliced to a common 3' coding exon (IX). Transcription of individual Bdnf variants, which all encode the same BDNF protein, is initiated at unique promoters upstream of each non-coding exon, enabling precise spatiotemporal and activity-dependent regulation of Bdnf expression. Although prior evidence suggests that Bdnf transcripts containing exon I (Bdnf I) or exon IV (Bdnf IV) are uniquely regulated by neuronal activity, the functional significance of different Bdnf transcript variants remains unclear. To investigate functional roles of activity-dependent Bdnf I and IV transcripts, we used a CRISPR activation (CRISPRa) system in which catalytically-dead Cas9 (dCas9) fused to a transcriptional activator (VPR) is targeted to individual Bdnf promoters with single guide RNAs (sgRNAs), resulting in transcript-specific Bdnf upregulation. Bdnf I upregulation is associated with gene expression changes linked to dendritic growth, while Bdnf IV upregulation is associated with genes that regulate protein catabolism. Upregulation of Bdnf I, but not Bdnf IV, increased mushroom spine density, volume, length, and head diameter, and also produced more complex dendritic arbors in cultured rat hippocampal neurons. In contrast, upregulation of Bdnf IV, but not Bdnf I, in the rat hippocampus attenuated contextual fear expression. Our data suggest that while Bdnf I and IV are both activity-dependent, BDNF produced from these promoters may serve unique cellular, synaptic, and behavioral functions.

Development of adaptive anoikis resistance promotes metastasis that can be overcome by CDK8/19 Mediator kinase inhibition.

Anoikis resistance or evasion of cell death triggered by cell detachment into suspension is a hallmark of cancer that is concurrent with cell survival and metastasis. The effects of frequent matrix detachment encounters on the development of anoikis resistance in cancer remains poorly defined. Here we show using a panel of ovarian cancer models, that repeated exposure to suspension stress in vitro followed by attached recovery growth leads to the development of anoikis resistance paralleling in vivo development of anoikis resistance in ovarian cancer ascites. This resistance is concurrent with enhanced invasion, chemoresistance and the ability of anoikis adapted cells to metastasize to distant sites. Adapted anoikis resistant cells show a heightened dependency on oxidative phosphorylation and can also evade immune surveillance. We find that such acquired anoikis resistance is not genetic, as acquired resistance persists for a finite duration in the absence of suspension stress. Transcriptional reprogramming is however essential to this process, as acquisition of adaptive anoikis resistance in vitro and in vivo is exquisitely sensitive to inhibition of CDK8/19 Mediator kinase, a pleiotropic regulator of transcriptional reprogramming. Our data demonstrate that growth after recovery from repeated exposure to suspension stress is a direct contributor to metastasis and that inhibition of CDK8/19 Mediator kinase during such adaptation provides a therapeutic opportunity to prevent both local and distant metastasis in cancer.

An atlas of transcriptionally defined cell populations in the rat ventral tegmental area.

The ventral tegmental area (VTA) is a complex brain region that is essential for reward function and frequently implicated in neuropsychiatric disease. While decades of research on VTA function have focused on dopamine neurons, recent evidence has identified critical roles for GABAergic and glutamatergic neurons in reward processes. Additionally, although subsets of VTA neurons express genes involved in the synthesis and transport of multiple neurotransmitters, characterization of these combinatorial populations has largely relied on low-throughput methods. To comprehensively define the molecular architecture of the VTA, we performed single-nucleus RNA sequencing on 21,600 cells from the rat VTA. Analysis of neuronal subclusters identifies selective markers for dopamine and combinatorial neurons, reveals expression profiles for receptors targeted by drugs of abuse, and demonstrates population-specific enrichment of gene sets linked to brain disorders. These results highlight the heterogeneity of the VTA and provide a resource for further exploration of VTA gene expression.

Histone macroH2A1 is a stronger regulator of hippocampal transcription and memory than macroH2A2 in mice.

Histone variants H2A.Z and H3.3 are epigenetic regulators of memory, but roles of other variants are not well characterized. macroH2A (mH2A) is a structurally unique histone that contains a globular macrodomain connected to the histone region by an unstructured linker. Here we assessed if mH2A regulates memory and if this role varies for the two mH2A-encoding genes, H2afy (mH2A1) and H2afy2 (mH2A2). We show that fear memory is impaired in mH2A1, but not in mH2A2-deficient mice, whereas both groups were impaired in a non-aversive spatial memory task. However, impairment was larger for mH2A1- deficient mice, indicating a preferential role for mH2A1 over mH2A2 in memory. Accordingly, mH2A1 depletion in the mouse hippocampus resulted in more extensive transcriptional de-repression compared to mH2A2 depletion. mH2A1-depleted mice failed to induce a normal transcriptional response to fear conditioning, suggesting that mH2A1 depletion impairs memory by altering transcription. Using chromatin immunoprecipitation (ChIP) sequencing, we found that both mH2A proteins are enriched on transcriptionally repressed genes, but only mH2A1 occupancy was dynamically modified during learning, displaying reduced occupancy on upregulated genes after training. These data identify mH2A as a regulator of memory and suggest that mH2A1 supports memory by repressing spurious transcription and promoting learning-induced transcriptional activation.

Regulation of dopamine-dependent transcription and cocaine action by Gadd45b.

Exposure to drugs of abuse produces robust transcriptional and epigenetic reorganization within brain reward circuits that outlives the direct effects of the drug and may contribute to addiction. DNA methylation is a covalent epigenetic modification that is altered following stimulant exposure and is critical for behavioral and physiological adaptations to drugs of abuse. Although activity-related loss of DNA methylation requires the Gadd45 (Growth arrest and DNA-damage-inducible) gene family, very little is known about how this family regulates activity within the nucleus accumbens or behavioral responses to drugs of abuse. Here, we combined genome-wide transcriptional profiling, pharmacological manipulations, electrophysiological measurements, and CRISPR tools with traditional knockout and behavioral approaches in rodent model systems to dissect the role of Gadd45b in dopamine-dependent epigenetic regulation and cocaine reward. We show that acute cocaine administration induces rapid upregulation of Gadd45b mRNA in the rat nucleus accumbens, and that knockout or site-specific CRISPR/Cas9 gene knockdown of Gadd45b blocks cocaine conditioned place preference. In vitro, dopamine treatment in primary striatal neurons increases Gadd45b mRNA expression through a dopamine receptor type 1 (DRD1)-dependent mechanism. Moreover, shRNA-induced Gadd45b knockdown decreases expression of genes involved in psychostimulant addiction, blocks induction of immediate early genes by DRD1 stimulation, and prevents DRD1-mediated changes in DNA methylation. Finally, we demonstrate that Gadd45b knockdown decreases striatal neuron action potential burst duration in vitro, without altering other electrophysiological characteristics. These results suggest that striatal Gadd45b functions as a dopamine-induced gene that is necessary for cocaine reward memory and DRD1-mediated transcriptional activity.

Profiling and quantification of pluripotency reprogramming reveal that WNT pathways and cell morphology have to be reprogramed extensively.

Pluripotent state can be established via reprogramming of somatic nuclei by factors within an oocyte or by ectopic expression of a few transgenes. Considered as being extensive and intensive, the full complement of genes to be reprogrammed, however, has never been defined, nor has the degree of reprogramming been determined quantitatively. Here, we propose a new concept of reprogramome, which is defined as the full complement of genes to be reprogrammed to the expression levels found in pluripotent stem cells (PSCs). This concept in combination with RNA-seq enables us to precisely profile reprogramome and sub-reprogramomes, and study the reprogramming process with the help of other available tools such as GO analyses. With reprogramming of human fibroblasts into PSCs as an example, we have defined the full complement of the human fibroblast-to-PSC reprogramome. Furthermore, our analyses of the reprogramome revealed that WNT pathways and genes with roles in cellular morphogenesis should be extensively and intensely reprogrammed for the establishment of pluripotency. We further developed a new mathematical model to quantitate the overall reprogramming, as well as reprogramming in a specific cellular feature such as WNT signaling pathways and genes regulating cellular morphogenesis. We anticipate that our concept and mathematical model may be applied to study and quantitate other reprogramming (pluripotency reprogramming from other somatic cells, and lineage reprogramming), as well as transcriptional and epigenetic differences between any two types of cells including cancer cells and their normal counterparts.