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.

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.

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.