The neuronal stimulation-transcription coupling map

Excitation-transcription coupling, neuronal gene expression and synaptic plasticity

Excitation-transcription coupling (E-TC) links synaptic and cellular activity to nuclear gene transcription. It is generally accepted that E-TC makes a crucial contribution to learning and memory through its role in underpinning long-lasting synaptic enhancement in late-phase long-term potentiation and has more recently been linked to late-phase long-term depression: both processes require de novo gene transcription, mRNA translation and protein synthesis. E-TC begins with the activation of glutamate-gated N-methyl-D-aspartate-type receptors and voltage-gated L-type Ca2+ channels at the membrane and culminates in the activation of transcription factors in the nucleus. These receptors and ion channels mediate E-TC through mechanisms that include long-range signalling from the synapse to the nucleus and local interactions within dendritic spines, among other possibilities. Growing experimental evidence links these E-TC mechanisms to late-phase long-term potentiation and learning and memory. These advances in our understanding of the molecular mechanisms of E-TC mean that future efforts can focus on understanding its mesoscale functions and how it regulates neuronal network activity and behaviour in physiological and pathological conditions.

Seizure enhances SUMOylation and zinc-finger transcriptional repression in neuronal nuclei

A single episode of pilocarpine-induced status epilepticus can trigger the development of spontaneous recurrent seizures in a rodent model for epilepsy. The initial seizure-induced events in neuronal nuclei that lead to long-term changes in gene expression and cellular responses likely contribute toward epileptogenesis. Using a transgenic mouse model to specifically isolate excitatory neuronal nuclei, we profiled the seizure-induced nuclear proteome via tandem mass tag mass spectrometry and observed robust enrichment of nuclear proteins associated with the SUMOylation pathway. In parallel with nuclear proteome, we characterized nuclear gene expression by RNA sequencing which provided insights into seizure-driven transcriptional regulation and dynamics. Strikingly, we saw widespread downregulation of zinc-finger transcription factors, specifically proteins that harbor Krüppel-associated box (KRAB) domains. Our results provide a detailed snapshot of nuclear events induced by seizure activity and demonstrate a robust method for cell-type-specific nuclear profiling that can be applied to other cell types and models.

Activity-regulated gene expression across cell types of the mouse hippocampus

Activity-regulated gene (ARG) expression patterns in the hippocampus (HPC) regulate synaptic plasticity, learning, and memory, and are linked to both risk and treatment responses for many neuropsychiatric disorders. The HPC contains discrete classes of neurons with specialized functions, but cell type-specific activity-regulated transcriptional programs are not well characterized. Here, we used single-nucleus RNA-sequencing (snRNA-seq) in a mouse model of acute electroconvulsive seizures (ECS) to identify cell type-specific molecular signatures associated with induced activity in HPC neurons. We used unsupervised clustering and a priori marker genes to computationally annotate 15,990 high-quality HPC neuronal nuclei from N = 4 mice across all major HPC subregions and neuron types. Activity-induced transcriptomic responses were divergent across neuron populations, with dentate granule cells being particularly responsive to activity. Differential expression analysis identified both upregulated and downregulated cell type-specific gene sets in neurons following ECS. Within these gene sets, we identified enrichment of pathways associated with varying biological processes such as synapse organization, cellular signaling, and transcriptional regulation. Finally, we used matrix factorization to reveal continuous gene expression patterns differentially associated with cell type, ECS, and biological processes. This work provides a rich resource for interrogating activity-regulated transcriptional responses in HPC neurons at single-nuclei resolution in the context of ECS, which can provide biological insight into the roles of defined neuronal subtypes in HPC function.

Transcriptional changes in specific subsets of Drosophila neurons following inhibition of the serotonin transporter

The transcriptional effects of SSRIs and other serotonergic drugs remain unclear, in part due to the heterogeneity of postsynaptic cells, which may respond differently to changes in serotonergic signaling. Relatively simple model systems such as Drosophila afford more tractable microcircuits in which to investigate these changes in specific cell types. Here, we focus on the mushroom body, an insect brain structure heavily innervated by serotonin and comprised of multiple different but related subtypes of Kenyon cells. We use fluorescence-activated cell sorting of Kenyon cells, followed by either bulk or single-cell RNA sequencing to explore the transcriptomic response of these cells to SERT inhibition. We compared the effects of two different Drosophila Serotonin Transporter (dSERT) mutant alleles as well as feeding the SSRI citalopram to adult flies. We find that the genetic architecture associated with one of the mutants contributed to significant artefactual changes in expression. Comparison of differential expression caused by loss of SERT during development versus aged, adult flies, suggests that changes in serotonergic signaling may have relatively stronger effects during development, consistent with behavioral studies in mice. Overall, our experiments revealed limited transcriptomic changes in Kenyon cells, but suggest that different subtypes may respond differently to SERT loss-of-function. Further work exploring the effects of SERT loss-of-function in other circuits may be used help to elucidate how SSRIs differentially affect a variety of different neuronal subtypes both during development and in adults.

Modulation of chromatin architecture influences the neuronal nucleus through activity-regulated gene expression

The disruption of chromatin-regulating genes is associated with many neurocognitive syndromes. While most of these genes are ubiquitously expressed across various cell-types, many chromatin regulators act upon activity regulated genes (ARGs) that play central roles in synaptic development and plasticity. Recent literature suggests a link between ARG expression disruption in neurons with the human phenotypes observed in various neurocognitive syndromes. Advances in chromatin biology have demonstrated how chromatin structure, from nucleosome occupancy to higher-order structures such as topologically associated domains, impacts the kinetics of transcription. This review discusses the dynamics of these various levels of chromatin structure and their influence on the expression of ARGs.

Molecular encoding of stimulus features in a single sensory neuron type enables neuronal and behavioral plasticity

Neurons modify their transcriptomes in response to an animal's experience. How specific experiences are transduced to modulate gene expression and precisely tune neuronal functions are not fully defined. Here, we describe the molecular profile of a thermosensory neuron pair in C. elegans experiencing different temperature stimuli. We find that distinct salient features of the temperature stimulus, including its duration, magnitude of change, and absolute value, are encoded in the gene expression program in this single neuron type, and we identify a novel transmembrane protein and a transcription factor whose specific transcriptional dynamics are essential to drive neuronal, behavioral, and developmental plasticity. Expression changes are driven by broadly expressed activity-dependent transcription factors and corresponding cis-regulatory elements that nevertheless direct neuron- and stimulus-specific gene expression programs. Our results indicate that coupling of defined stimulus characteristics to the gene regulatory logic in individual specialized neuron types can customize neuronal properties to drive precise behavioral adaptation.

Transcriptional changes in specific subsets of Drosophila neurons following inhibition of the serotonin transporter

The transcriptional effects of SSRIs and other serotonergic drugs remain unclear, in part due to the heterogeneity of postsynaptic cells, which may respond differently to changes in serotonergic signaling. Relatively simple model systems such as Drosophila afford more tractable microcircuits in which to investigate these changes in specific cell types. Here, we focus on the mushroom body, an insect brain structure heavily innervated by serotonin and comprised of multiple different but related subtypes of Kenyon cells. We use fluorescence activated cell sorting of Kenyon cells, followed by either or bulk or single cell RNA sequencing to explore the transcriptomic response of these cells to SERT inhibition. We compared the effects of two different Drosophila Serotonin Transporter (dSERT) mutant alleles as well as feeding the SSRI citalapram to adult flies. We find that the genetic architecture associated with one of the mutants contributed to significant artefactual changes in expression. Comparison of differential expression caused by loss of SERT during development versus aged, adult flies, suggests that changes in serotonergic signaling may have relatively stronger effects during development, consistent with behavioral studies in mice. Overall, our experiments revealed limited transcriptomic changes in Kenyon cells, but suggest that different subtypes may respond differently to SERT loss-of-function. Further work exploring the effects of SERT loss-of-function in other Drosophila circuits may be used help to elucidate how SSRIs differentially affect a variety of different neuronal subtypes both during development and in adults.

Molecular encoding of stimulus features in a single sensory neuron type enables neuronal and behavioral plasticity

Neurons modify their transcriptomes in response to an animal’s experience. How specific experiences are transduced to modulate gene expression and precisely tune neuronal functions are not fully defined. Here, we describe the molecular profile of a thermosensory neuron pair in C. elegans experiencing different temperature stimuli. We find that distinct salient features of the temperature stimulus including its duration, magnitude of change, and absolute value are encoded in the gene expression program in this single neuron, and identify a novel transmembrane protein and a transcription factor whose specific transcriptional dynamics are essential to drive neuronal, behavioral, and developmental plasticity. Expression changes are driven by broadly expressed activity-dependent transcription factors and corresponding cis -regulatory elements that nevertheless direct neuron- and stimulus-specific gene expression programs. Our results indicate that coupling of defined stimulus characteristics to the gene regulatory logic in individual specialized neuron types can customize neuronal properties to drive precise behavioral adaptation.

Acute Fasting Modulates Food-Seeking Behavior and Neural Signaling in the Piriform Cortex

It is well known that the state of hunger can modulate hormones and hypothalamic neural circuits to drive food-seeking behavior and consumption. However, the role the sensory cortex plays in regulating foraging is much less explored. Here, we investigated whether acute fasting in mice can alter an odor-guided foraging behavior and how it can alter neurons and synapses in the (olfactory) piriform cortex (PC). Acute hunger enhances the motivation of a mouse to search for food pellets and increases food intake. The foraging behavior strongly activates the PC, as revealed by c-Fos immunostaining. The activation of PC is accompanied by an increase in excitation-inhibition ratio of synaptic density. Fasting also enhances the phosphorylation of AMP kinase, a biochemical energy regulator. Taken together, our results uncover a new regulatory brain region and implicate the PC in controlling foraging behavior.

Stimulus-specific remodeling of the neuronal transcriptome through nuclear intron-retaining transcripts

The nuclear envelope has long been considered primarily a physical barrier separating nuclear and cytosolic contents. More recently, nuclear compartmentalization has been shown to have additional regulatory functions in controlling gene expression. A sizeable proportion of protein-coding mRNAs is more prevalent in the nucleus than in the cytosol, suggesting regulated mRNA trafficking to the cytosol, but the mechanisms underlying controlled nuclear mRNA retention remain unclear. Here, we provide a comprehensive map of the subcellular localization of mRNAs in mature mouse cortical neurons, and reveal that transcripts retained in the nucleus comprise the majority of stable intron-retaining mRNAs. Systematically probing the fate of nuclear transcripts upon neuronal stimulation, we found opposite effects on sub-populations of transcripts: while some are targeted for degradation, others complete splicing to generate fully mature mRNAs that are exported to the cytosol and mediate rapid increases in protein levels. Finally, different forms of stimulation mobilize distinct groups of intron-retaining transcripts, with this selectivity arising from the activation of specific signaling pathways. Overall, our findings uncover a cue-specific control of intron retention as a major regulator of acute remodeling of the neuronal transcriptome.

Nuclear RNA binding regulates TDP-43 nuclear localization and passive nuclear export

Nuclear clearance of the RNA-binding protein TDP-43 is a hallmark of neurodegeneration and an important therapeutic target. Our current understanding of TDP-43 nucleocytoplasmic transport does not fully explain its predominantly nuclear localization or mislocalization in disease. Here, we show that TDP-43 exits nuclei by passive diffusion, independent of facilitated mRNA export. RNA polymerase II blockade and RNase treatment induce TDP-43 nuclear efflux, suggesting that nuclear RNAs sequester TDP-43 in nuclei and limit its availability for passive export. Induction of TDP-43 nuclear efflux by short, GU-rich oligomers (presumably by outcompeting TDP-43 binding to endogenous nuclear RNAs), and nuclear retention conferred by splicing inhibition, demonstrate that nuclear TDP-43 localization depends on binding to GU-rich nuclear RNAs. Indeed, RNA-binding domain mutations markedly reduce TDP-43 nuclear localization and abolish transcription blockade-induced nuclear efflux. Thus, the nuclear abundance of GU-RNAs, dictated by the balance of transcription, pre-mRNA processing, and RNA export, regulates TDP-43 nuclear localization.

Epigenetic and Transcriptional Regulation of Spontaneous and Sensory Activity Dependent Programs During Neuronal Circuit Development

Spontaneous activity generated before the onset of sensory transduction has a key role in wiring developing sensory circuits. From axonal targeting, to synapse formation and elimination, to the balanced integration of neurons into developing circuits, this type of activity is implicated in a variety of cellular processes. However, little is known about its molecular mechanisms of action, especially at the level of genome regulation. Conversely, sensory experience-dependent activity implements well-characterized transcriptional and epigenetic chromatin programs that underlie heterogeneous but specific genomic responses that shape both postnatal circuit development and neuroplasticity in the adult. In this review, we focus on our knowledge of the developmental processes regulated by spontaneous activity and the underlying transcriptional mechanisms. We also review novel findings on how chromatin regulates the specificity and developmental induction of the experience-dependent program, and speculate their relevance for our understanding of how spontaneous activity may act at the genomic level to instruct circuit assembly and prepare developing neurons for sensory-dependent connectivity refinement and processing.

A transcriptional rheostat couples past activity to future sensory responses

Animals traversing different environments encounter both stable background stimuli and novel cues, which are thought to be detected by primary sensory neurons and then distinguished by downstream brain circuits. Here, we show that each of the ∼1,000 olfactory sensory neuron (OSN) subtypes in the mouse harbors a distinct transcriptome whose content is precisely determined by interactions between its odorant receptor and the environment. This transcriptional variation is systematically organized to support sensory adaptation: expression levels of more than 70 genes relevant to transforming odors into spikes continuously vary across OSN subtypes, dynamically adjust to new environments over hours, and accurately predict acute OSN-specific odor responses. The sensory periphery therefore separates salient signals from predictable background via a transcriptional rheostat whose moment-to-moment state reflects the past and constrains the future; these findings suggest a general model in which structured transcriptional variation within a cell type reflects individual experience.

RyR-mediated Ca2+ release elicited by neuronal activity induces nuclear Ca2+ signals, CREB phosphorylation, and Npas4/RyR2 expression

The expression of several hippocampal genes implicated in learning and memory processes requires that Ca²⁺ signals generated in dendritic spines, dendrites, or the soma in response to neuronal stimulation reach the nucleus. The diffusion of Ca²⁺ in the cytoplasm is highly restricted, so neurons must use other mechanisms to propagate Ca²⁺ signals to the nucleus. Here, we present evidence showing that Ca²⁺ release mediated by the ryanodine receptor (RyR) channel type-2 isoform (RyR2) contributes to the generation of nuclear Ca²⁺ signals induced by gabazine (GBZ) addition, glutamate uncaging in the dendrites, or high-frequency field stimulation of primary hippocampal neurons. Additionally, GBZ treatment significantly increased cyclic adenosine monophosphate response element binding protein (CREB) phosphorylation-a key event in synaptic plasticity and hippocampal memory-and enhanced the expression of Neuronal Per Arnt Sim domain protein 4 (Npas4) and RyR2, two central regulators of these processes. Suppression of RyR-mediated Ca²⁺ release with ryanodine significantly reduced the increase in CREB phosphorylation and the enhanced Npas4 and RyR2 expression induced by GBZ. We propose that RyR-mediated Ca²⁺ release induced by neuronal activity, through its contribution to the sequential generation of nuclear Ca²⁺ signals, CREB phosphorylation, Npas4, and RyR2 up-regulation, plays a central role in hippocampal synaptic plasticity and memory processes.

Sleep Loss Drives Brain Region-Specific and Cell Type-Specific Alterations in Ribosome-Associated Transcripts Involved in Synaptic Plasticity and Cellular Timekeeping

Sleep and sleep loss are thought to impact synaptic plasticity, and recent studies have shown that sleep and sleep deprivation (SD) differentially affect gene transcription and protein translation in the mammalian forebrain. However, much less is known regarding how sleep and SD affect these processes in different microcircuit elements within the hippocampus and neocortex, for example, in inhibitory versus excitatory neurons. Here, we use translating ribosome affinity purification (TRAP) and in situ hybridization to characterize the effects of sleep versus SD on abundance of ribosome-associated transcripts in Camk2a-expressing (Camk2a+) pyramidal neurons and parvalbumin-expressing (PV+) interneurons in the hippocampus and neocortex of male mice. We find that while both Camk2a+ neurons and PV+ interneurons in neocortex show concurrent SD-driven increases in ribosome-associated transcripts for activity-regulated effectors of plasticity and transcriptional regulation, these transcripts are minimally affected by SD in hippocampus. Similarly, we find that while SD alters several ribosome-associated transcripts involved in cellular timekeeping in neocortical Camk2a+ and PV+ neurons, effects on circadian clock transcripts in hippocampus are minimal, and restricted to Camk2a+ neurons. Taken together, our results indicate that SD effects on transcripts associated with translating ribosomes are both cell type-specific and brain region-specific, and that these effects are substantially more pronounced in the neocortex than the hippocampus. We conclude that SD-driven alterations in the strength of synapses, excitatory-inhibitory (E-I) balance, and cellular timekeeping are likely more heterogeneous than previously appreciated.SIGNIFICANCE STATEMENT Sleep loss-driven changes in transcript and protein abundance have been used as a means to better understand the function of sleep for the brain. Here, we use translating ribosome affinity purification (TRAP) to characterize changes in abundance of ribosome-associated transcripts in excitatory and inhibitory neurons in mouse hippocampus and neocortex after a brief period of sleep or sleep loss. We show that these changes are not uniform, but are generally more pronounced in excitatory neurons than inhibitory neurons, and more pronounced in neocortex than in hippocampus.

The m6A-epitranscriptome in brain plasticity, learning and memory

Activity-dependent gene expression and protein translation underlie the ability of neurons to dynamically adjust their synaptic strength in response to sensory experience and during learning. The emerging field of epitranscriptomics (RNA modifications) has rapidly shifted our views on the mechanisms that regulate gene expression. Among hundreds of biochemical modifications on RNA, N⁶-methyladenosine (m⁶A) is the most abundant reversible mRNA modification in the brain. Its dynamic nature and ability to regulate all aspects of mRNA processing have positioned m⁶A as an important and versatile regulator of nervous system functions, including neuronal plasticity, learning and memory. In this review, we summarise recent experimental evidence that supports the role of m⁶A signalling in learning and memory, as well as providing an overview of the underlying molecular mechanisms in neurons. We also discuss the consequences of perturbed m⁶A signalling and/or its regulatory networks which are increasingly being linked to various cognitive disorders in humans.

Egr2 induction in spiny projection neurons of the ventrolateral striatum contributes to cocaine place preference in mice

Drug addiction develops due to brain-wide plasticity within neuronal ensembles, mediated by dynamic gene expression. Though the most common approach to identify such ensembles relies on immediate early gene expression, little is known of how the activity of these genes is linked to modified behavior observed following repeated drug exposure. To address this gap, we present a broad-to-specific approach, beginning with a comprehensive investigation of brain-wide cocaine-driven gene expression, through the description of dynamic spatial patterns of gene induction in subregions of the striatum, and finally address functionality of region-specific gene induction in the development of cocaine preference. Our findings reveal differential cell-type specific dynamic transcriptional recruitment patterns within two subdomains of the dorsal striatum following repeated cocaine exposure. Furthermore, we demonstrate that induction of the IEG Egr2 in the ventrolateral striatum, as well as the cells within which it is expressed, are required for the development of cocaine seeking.

The human brain is ever changing, constantly rewiring itself in response to new experiences, knowledge or information from the environment. Addictive drugs such as cocaine can hijack the genetic mechanisms responsible for this plasticity, creating dangerous, obsessive drug-seeking and consuming behaviors.

Cocaine-induced plasticity is difficult to apprehend, however, as brain regions or even cell populations can react differently to the compound. For instance, sub-regions in the striatum – the brain area that responds to rewards and helps to plan movement – show distinct responses during progressive exposure to cocaine. And while researchers know that the drug immediately changes how neurons switch certain genes on and off, it is still unclear how these genetic modifications later affect behavior.

Mukherjee, Gonzales et al. explored these questions at different scales, first focusing on how progressive cocaine exposure changed the way various gene programs were activated across the entire brain. This revealed that programs in the striatum were the most affected by the drug.

Examining this region more closely showed that cocaine switches on genes in specific ‘spiny projection’ neuron populations, depending on where these cells are located and the drug history of the mouse. Finally, Mukherjee, Gonzales et al. used genetically modified mice to piece together cocaine exposure, genetic changes and modifications in behavior. These experiments revealed that the drive to seek cocaine depended on activation of the Egr2 gene in populations of spiny projection neurons in a specific sub-region of the striatum. The gene, which codes for a protein that regulates how genes are switched on and off, was itself strongly activated by cocaine intake.

Cocaine addiction can have devastating consequences for individuals. Grasping how this drug alters the brain could pave the way for new treatments, while also providing information on the basic mechanisms underlying brain plasticity.

Activity labeling in vivo using CaMPARI2 reveals intrinsic and synaptic differences between neurons with high and low firing rate set points

Neocortical pyramidal neurons regulate firing around a stable mean firing rate (FR) that can differ by orders of magnitude between neurons, but the factors that determine where individual neurons sit within this broad FR distribution are not understood. To access low- and high-FR neurons for ex vivo analysis, we used Ca2+- and UV-dependent photoconversion of CaMPARI2 in vivo to permanently label neurons according to mean FR. CaMPARI2 photoconversion was correlated with immediate early gene expression and higher FRs ex vivo and tracked the drop and rebound in ensemble mean FR induced by prolonged monocular deprivation. High-activity L4 pyramidal neurons had greater intrinsic excitability and recurrent excitatory synaptic strength, while E/I ratio, local output strength, and local connection probability were not different. Thus, in L4 pyramidal neurons (considered a single transcriptional cell type), a broad mean FR distribution is achieved through graded differences in both intrinsic and synaptic properties.

Transcriptome-Wide Analysis of Interplay between mRNA Stability, Translation and Small RNAs in Response to Neuronal Membrane Depolarization

Experience-dependent changes to neural circuitry are shaped by spatially-restricted activity-dependent mRNA translation. Although the complexity of mRNA translation in neuronal cells is widely appreciated, translational profiles associated with neuronal excitation remain largely uncharacterized, and the associated regulatory mechanisms are poorly understood. Here, we employed ribosome profiling, mRNA sequencing and small RNA sequencing to profile transcriptome-wide changes in mRNA translation after whole cell depolarization of differentiated neuroblast cultures, and investigate the contribution of sequence-specific regulatory mechanisms. Immediately after depolarization, a functional partition between transcriptional and translational responses was uncovered, in which many mRNAs were subjected to significant changes in abundance or ribosomal occupancy, but not both. After an extended (2 h) post-stimulus rest phase, however, these changes became synchronized, suggesting that there are different layers of post-transcriptional regulation which are temporally separated but become coordinated over time. Globally, changes in mRNA abundance and translation were found to be associated with a number of intrinsic mRNA features, including mRNA length, GC% and secondary structures; however, the effect of these factors differed between both post-depolarization time-points. Furthermore, small RNA sequencing revealed that miRNAs and tRNA-derived small RNA fragments were subjected to peak changes in expression immediately after stimulation, during which these molecules were predominantly associated with fluctuations in mRNA abundance, consistent with known regulatory mechanisms. These data suggest that excitation-associated neuronal translation is subjected to extensive temporal coordination, with substantial contributions from a number of sequence-dependent regulatory mechanisms.

MEF2C and HDAC5 regulate Egr1 and Arc genes to increase dendritic spine density and complexity in early enriched environment

We investigated the effects of environmental enrichment during critical period of early postnatal life and how it interplays with the epigenome to affect experience-dependent visual cortical plasticity. Mice raised in an EE from birth to during CP have increased spine density and dendritic complexity in the visual cortex. EE upregulates synaptic plasticity genes, Arc and Egr1, and a transcription factor MEF2C. We also observed an increase in MEF2C binding to the promoters of Arc and Egr1. In addition, pups raised in EE show a reduction in HDAC5 and its binding to promoters of Mef2c, Arc and Egr1 genes. With an overexpression of Mef2c, neurite outgrowth increased in complexity. Our results suggest a possible underlying molecular mechanism of EE, acting through MEF2C and HDAC5, which drive Arc and Egr1. This could lead to the observed increased dendritic spine density and complexity induced by early EE.

Subregion-specific rules govern the distribution of neuronal immediate-early gene induction

The induction of immediate-early gene (IEG) expression in brain nuclei in response to an experience is necessary for the formation of long-term memories. Additionally, the rapid dynamics of IEG induction and decay motivates the common use of IEG expression as markers for identification of neuronal assemblies ("ensembles") encoding recent experience. However, major gaps remain in understanding the rules governing the distribution of IEGs within neuronal assemblies. Thus, the extent of correlation between coexpressed IEGs, the cell specificity of IEG expression, and the spatial distribution of IEG expression have not been comprehensively studied. To address these gaps, we utilized quantitative multiplexed single-molecule fluorescence in situ hybridization (smFISH) and measured the expression of IEGs (Arc, Egr2, and Nr4a1) within spiny projection neurons (SPNs) in the dorsal striatum of mice following acute exposure to cocaine. Exploring the relevance of our observations to other brain structures and stimuli, we also analyzed data from a study of single-cell RNA sequencing of mouse cortical neurons. We found that while IEG expression is graded, the expression of multiple IEGs is tightly correlated at the level of individual neurons. Interestingly, we observed that region-specific rules govern the induction of IEGs in SPN subtypes within striatal subdomains. We further observed that IEG-expressing assemblies form spatially defined clusters within which the extent of IEG expression correlates with cluster size. Together, our results suggest the existence of IEG-expressing neuronal "superensembles," which are associated in spatial clusters and characterized by coherent and robust expression of multiple IEGs.