EphB1 controls long-range cortical axon guidance through a cell non-autonomous role in GABAergic cells.

EphB1 is required for proper guidance of cortical axon projections during brain development, but how EphB1 regulates this process remains unclear. We show here that EphB1 conditional knockout (cKO) in GABAergic cells (Vgat-Cre), but not in cortical excitatory neurons (Emx1-Cre), reproduced the cortical axon guidance defects observed in global EphB1 KO mice. Interestingly, in EphB1 cKOVgat mice, the misguided axon bundles contained co-mingled striatal GABAergic and somatosensory cortical glutamatergic axons. In wild-type mice, somatosensory axons also co-fasciculated with striatal axons, notably in the globus pallidus, suggesting that a subset of glutamatergic cortical axons normally follows long-range GABAergic axons to reach their targets. Surprisingly, the ectopic axons in EphB1 KO mice were juxtaposed to major blood vessels. However, conditional loss of EphB1 in endothelial cells (Tie2-Cre) did not produce the axon guidance defects, suggesting that EphB1 in GABAergic neurons normally promotes avoidance of these ectopic axons from the developing brain vasculature. Together, our data reveal a new role for EphB1 in GABAergic neurons to influence proper cortical glutamatergic axon guidance during brain development.

PRMT5-mediated arginine methylation of FXR1 is essential for RNA binding in cancer cells.

Emerging evidence indicates that arginine methylation promotes the stability of arginine-glycine-rich (RGG) motif-containing RNA-binding proteins (RBPs) and regulates gene expression. Here, we report that post-translational modification of FXR1 enhances the binding with mRNAs and is involved in cancer cell growth and proliferation. Independent point mutations in arginine residues of FXR1's nuclear export signal (R386 and R388) and RGG (R453, R455 and R459) domains prevent it from binding to RNAs that form G-quadruplex (G4) RNA structures. Disruption of G4-RNA structures by lithium chloride failed to bind with FXR1, indicating its preference for G4-RNA structure containing mRNAs. Furthermore, loss-of-function of PRMT5 inhibited FXR1 methylation both in vivo and in vitro, affecting FXR1 protein stability, inhibiting RNA-binding activity and cancer cell growth and proliferation. Finally, the enhanced crosslinking and immunoprecipitation (eCLIP) analyses reveal that FXR1 binds with the G4-enriched mRNA targets such as AHNAK, MAP1B, AHNAK2, HUWE1, DYNC1H1 and UBR4 and controls its mRNA expression in cancer cells. Our findings suggest that PRMT5-mediated FXR1 methylation is required for RNA/G4-RNA binding, which promotes gene expression in cancer cells. Thus, FXR1's structural characteristics and affinity for RNAs preferentially G4 regions provide new insights into the molecular mechanism of FXR1 in oral cancer cells.

Epigenetic function during heroin self-administration controls future relapse-associated behavior in a cell type-specific manner.

Opioid use produces enduring associations between drug reinforcement/euphoria and discreet or diffuse cues in the drug-taking environment. These powerful associations can trigger relapse in individuals recovering from opioid use disorder (OUD). Here, we sought to determine whether the epigenetic enzyme, histone deacetylase 5 (HDAC5), regulates relapse-associated behavior in an animal model of OUD. We examined the effects of nucleus accumbens (NAc) HDAC5 on both heroin- and sucrose-seeking behaviors using operant self-administration paradigms. We utilized cre-dependent viral-mediated approaches to investigate the cell-type-specific effects of HDAC5 on heroin-seeking behavior, gene expression, and medium spiny neuron (MSN) cell and synaptic physiology. We found that NAc HDAC5 functions during the acquisition phase of heroin self-administration to limit future relapse-associated behavior. Moreover, overexpressing HDAC5 in the NAc suppressed context-associated and reinstated heroin-seeking behaviors, but it did not alter sucrose seeking. We also found that HDAC5 functions within dopamine D1 receptor-expressing MSNs to suppress cue-induced heroin seeking, and within dopamine D2 receptor-expressing MSNs to suppress drug-primed heroin seeking. Assessing cell-type-specific transcriptomics, we found that HDAC5 reduced expression of multiple ion transport genes in both D1- and D2-MSNs. Consistent with this observation, HDAC5 also produced firing rate depression in both MSN classes. These findings revealed roles for HDAC5 during active heroin use in both D1- and D2-MSNs to limit distinct triggers of drug-seeking behavior. Together, our results suggest that HDAC5 might limit relapse vulnerability through regulation of ion channel gene expression and suppression of MSN firing rates during active heroin use.

scToppR: a coding-friendly R interface to ToppGene.

MOTIVATION: The scToppR package provides a ToppGene interface from R programs/scripts to fully access/control the database for functional enrichment without the need for active interaction on its Web site (https://toppgene.cchmc.org/). RESULTS: The library facilitates the functional enrichment analysis and visualization by interacting with ToppGene, downloading the functional enrichment dataframes, and using R environment to visualize the final results. AVAILABILITY: Code and documentation are currently available at https://github.com/BioinformaticsMUSC/scToppR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Human-unique brain cell clusters are associated with learning disorders and human episodic memory activity.

The advanced evolution of the human cerebral cortex forms the basis for our high-level cognitive functions. Through a comparative analysis of single-nucleus transcriptome data from the human neocortex and that of chimpanzees, macaques, and marmosets, we discovered 20 subgroups of cell types unique to the human brain, which include 11 types of excitatory neurons. Many of these human-unique cell clusters exhibit significant overexpression of genes regulated by human-specific enhancers. Notably, these specific cell clusters also express genes associated with disease risk, particularly those related to brain dysfunctions like learning disorders. Furthermore, genes linked to cortical thickness and human episodic memory encoding activities show heightened expression within these cell subgroups. These findings underscore the critical role of human brain-unique cell clusters in the evolution of human brain functions.

Novel transcriptional networks regulated by CLOCK in human neurons.

The molecular mechanisms underlying human brain evolution are not fully understood; however, previous work suggested that expression of the transcription factor CLOCK in the human cortex might be relevant to human cognition and disease. In this study, we investigated this novel transcriptional role for CLOCK in human neurons by performing chromatin immunoprecipitation sequencing for endogenous CLOCK in adult neocortices and RNA sequencing following CLOCK knockdown in differentiated human neurons in vitro. These data suggested that CLOCK regulates the expression of genes involved in neuronal migration, and a functional assay showed that CLOCK knockdown increased neuronal migratory distance. Furthermore, dysregulation of CLOCK disrupts coexpressed networks of genes implicated in neuropsychiatric disorders, and the expression of these networks is driven by hub genes with human-specific patterns of expression. These data support a role for CLOCK-regulated transcriptional cascades involved in human brain evolution and function.

Role of macrophages in regression of myocardial fibrosis following alleviation of left ventricular pressure overload.

Sustained hemodynamic pressure overload (PO) produced by murine transverse aortic constriction (TAC) causes myocardial fibrosis; removal of TAC (unTAC) returns left ventricle (LV) hemodynamic load to normal and results in significant, but incomplete regression of myocardial fibrosis. However, the cellular mechanisms that result in these outcomes have not been defined. The objective was to determine temporal changes in myocardial macrophage phenotype in TAC and unTAC and determine whether macrophage depletion alters collagen degradation after unTAC. Myocardial macrophage abundance and phenotype were assessed by immunohistochemistry, flow cytometry, and gene expression by RT-PCR in control (non-TAC), 2 wk, 4 wk TAC, and 2 wk, 4 wk, and 6 wk unTAC. Myocardial cytokine profiles and collagen-degrading enzymes were determined by immunoassay and immunoblots. Initial collagen degradation was detected with collagen-hybridizing peptide (CHP). At unTAC, macrophages were depleted with clodronate liposomes, and endpoints were measured at 2 wk unTAC. Macrophage number had a defined temporal pattern: increased in 2 wk and 4 wk TAC, followed by increases at 2 wk unTAC (over 4 wk TAC) that then decreased at 4 wk and 6 wk unTAC. At 2 wk unTAC, macrophage area was significantly increased and was regionally associated with CHP reactivity. Cytokine profiles in unTAC reflected a proinflammatory milieu versus the TAC-induced profibrotic milieu. Single-cell sequencing analysis of 2 wk TAC versus 2 and 6 wk unTAC revealed distinct macrophage gene expression profiles at each time point demonstrating unique macrophage populations in unTAC versus TAC myocardium. Clodronate liposome depletion at unTAC reduced CHP reactivity and decreased cathepsin K and proMMP2. We conclude that temporal changes in number and phenotype of macrophages play a critical role in both TAC-induced development and unTAC-mediated partial, but incomplete, regression of myocardial fibrosis.NEW & NOTEWORTHY Our novel findings highlight the dynamic changes in myocardial macrophage populations that occur in response to PO and after alleviation of PO. Our data demonstrated, for the first time, a potential benefit of macrophages in contributing to collagen degradation and the partial regression of interstitial fibrosis following normalization of hemodynamic load.

Optimization of a mouse model of pancreatic cancer to simulate the human phenotypes of metastasis and cachexia.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) presents with a high mortality rate. Two important features of PDAC contribute to this poor outcome. The first is metastasis which occurs in ~ 80% of PDAC patients. The second is cachexia, which compromises treatment tolerance for patients and reduces their quality of life. Although various mouse models of PDAC exist, recapitulating both metastatic and cachectic features have been challenging. METHODS: Here, we optimize an orthotopic mouse model of PDAC by altering several conditions, including the subcloning of parental murine PDAC cells, implantation site, number of transplanted cells, and age of recipient mice. We perform spatial profiling to compare primary and metastatic immune microenvironments and RNA sequencing to gain insight into the mechanisms of muscle wasting in PDAC-induced cachexia, comparing non-metastatic to metastatic conditions. RESULTS: These modifications extend the time course of the disease and concurrently increase the rate of metastasis to approximately 70%. Furthermore, reliable cachexia endpoints are achieved in both PDAC mice with and without metastases, which is reminiscent of patients. We also find that cachectic muscles from PDAC mice with metastasis exhibit a similar transcriptional profile to muscles derived from mice and patients without metastasis. CONCLUSION: Together, this model is likely to be advantageous in both advancing our understanding of the mechanism of PDAC cachexia, as well as in the evaluation of novel therapeutics.

Adaptor protein complex 2 in the orbitofrontal cortex predicts alcohol use disorder.

Alcohol use disorder (AUD) is a life-threatening disease characterized by compulsive drinking, cognitive deficits, and social impairment that continue despite negative consequences. The inability of individuals with AUD to regulate drinking may involve functional deficits in cortical areas that normally balance actions that have aspects of both reward and risk. Among these, the orbitofrontal cortex (OFC) is critically involved in goal-directed behavior and is thought to maintain a representation of reward value that guides decision making. In the present study, we analyzed post-mortem OFC brain samples collected from age- and sex-matched control subjects and those with AUD using proteomics, bioinformatics, machine learning, and reverse genetics approaches. Of the 4,500+ total unique proteins identified in the proteomics screen, there were 47 proteins that differed significantly by sex that were enriched in processes regulating extracellular matrix and axonal structure. Gene ontology enrichment analysis revealed that proteins differentially expressed in AUD cases were involved in synaptic and mitochondrial function, as well as transmembrane transporter activity. Alcohol-sensitive OFC proteins also mapped to abnormal social behaviors and social interactions. Machine learning analysis of the post-mortem OFC proteome revealed dysregulation of presynaptic (e.g., AP2A1) and mitochondrial proteins that predicted the occurrence and severity of AUD. Using a reverse genetics approach to validate a target protein, we found that prefrontal Ap2a1 expression significantly correlated with voluntary alcohol drinking in male and female genetically diverse mouse strains. Moreover, recombinant inbred strains that inherited the C57BL/6J allele at the Ap2a1 interval consumed higher amounts of alcohol than those that inherited the DBA/2J allele. Together, these findings highlight the impact of excessive alcohol consumption on the human OFC proteome and identify important cross-species cortical mechanisms and proteins that control drinking in individuals with AUD.

Big tau aggregation disrupts microtubule tyrosination and causes myocardial diastolic dysfunction: from discovery to therapy.

BACKGROUND: Amyloid plaques and neurofibrillary tangles, the molecular lesions that characterize Alzheimer's disease (AD) and other forms of dementia, are emerging as determinants of proteinopathies 'beyond the brain'. This study aims to establish tau's putative pathophysiological mechanistic roles and potential future therapeutic targeting of tau in heart failure (HF). METHODS AND RESULTS: A mouse model of tauopathy and human myocardial and brain tissue from patients with HF, AD, and controls was employed in this study. Tau protein expression was examined together with its distribution, and in vitro tau-related pathophysiological mechanisms were identified using a variety of biochemical, imaging, and functional approaches. A novel tau-targeting immunotherapy was tested to explore tau-targeted therapeutic potential in HF. Tau is expressed in normal and diseased human hearts, in contradistinction to the current oft-cited observation that tau is expressed specifically in the brain. Notably, the main cardiac isoform is high-molecular-weight (HMW) tau (also known as big tau), and hyperphosphorylated tau segregates in aggregates in HF and AD hearts. As previously described for amyloid-beta, the tauopathy phenotype in human myocardium is of diastolic dysfunction. Perturbation in the tubulin code, specifically a loss of tyrosinated microtubules, emerged as a potential mechanism of myocardial tauopathy. Monoclonal anti-tau antibody therapy improved myocardial function and clearance of toxic aggregates in mice, supporting tau as a potential target for novel HF immunotherapy. CONCLUSION: The study presents new mechanistic evidence and potential treatment for the brain-heart tauopathy axis in myocardial and brain degenerative diseases and ageing.

FoxP1 orchestration of ASD-relevant signaling pathways in the striatum.

Mutations in the transcription factor Forkhead box p1 (FOXP1) are causative for neurodevelopmental disorders such as autism. However, the function of FOXP1 within the brain remains largely uncharacterized. Here, we identify the gene expression program regulated by FoxP1 in both human neural cells and patient-relevant heterozygous Foxp1 mouse brains. We demonstrate a role for FoxP1 in the transcriptional regulation of autism-related pathways as well as genes involved in neuronal activity. We show that Foxp1 regulates the excitability of striatal medium spiny neurons and that reduction of Foxp1 correlates with defects in ultrasonic vocalizations. Finally, we demonstrate that FoxP1 has an evolutionarily conserved role in regulating pathways involved in striatal neuron identity through gene expression studies in human neural progenitors with altered FOXP1 levels. These data support an integral role for FoxP1 in regulating signaling pathways vulnerable in autism and the specific regulation of striatal pathways important for vocal communication.

Genomics at cellular resolution: insights into cognitive disorders and their evolution.

High-throughput genomic sequencing approaches have held the promise of understanding and ultimately leading to treatments for cognitive disorders such as autism spectrum disorders, schizophrenia and Alzheimer's disease. Although significant progress has been made into identifying genetic variants associated with these diseases, these studies have also uncovered that these disorders are mostly genetically complex and thus challenging to model in non-human systems. Improvements in such models might benefit from understanding the evolution of the human genome and how such modifications have affected brain development and function. The intersection of genome-wide variant information with cell-type-specific expression and epigenetic information will further assist in resolving the contribution of particular cell types in evolution or disease. For example, the role of non-neuronal cells in brain evolution and cognitive disorders has gone mostly underappreciated until the recent availability of single-cell transcriptomic approaches. In this review, we discuss recent studies that carry out cell-type-specific assessments of gene expression in brain tissue across primates and between healthy and disease populations. The emerging results from these studies are beginning to elucidate how specific cell types in the evolved human brain are contributing to cognitive disorders.

Hippocampal subfield transcriptome analysis in schizophrenia psychosis.

We have previously demonstrated functional and molecular changes in hippocampal subfields in individuals with schizophrenia (SZ) psychosis associated with hippocampal excitability. In this study, we use RNA-seq and assess global transcriptome changes in the hippocampal subfields, DG, CA3, and CA1 from individuals with SZ psychosis and controls to elucidate subfield-relevant molecular changes. We also examine changes in gene expression due to antipsychotic medication in the hippocampal subfields from our SZ ON- and OFF-antipsychotic medication cohort. We identify unique subfield-specific molecular profiles in schizophrenia postmortem samples compared with controls, implicating astrocytes in DG, immune mechanisms in CA3, and synaptic scaling in CA1. We show a unique pattern of subfield-specific effects by antipsychotic medication on gene expression levels with scant overlap of genes differentially expressed by SZ disease effect versus medication effect. These hippocampal subfield changes serve to confirm and extend our previous model of SZ and can explain the lack of full efficacy of conventional antipsychotic medication on SZ symptomatology. With future characterization using single-cell studies, the identified distinct molecular profiles of the DG, CA3, and CA1 in SZ psychosis may serve to identify further potential hippocampal-based therapeutic targets.

Accelerated evolution of oligodendrocytes in the human brain.

Recent discussions of human brain evolution have largely focused on increased neuron numbers and changes in their connectivity and expression. However, it is increasingly appreciated that oligodendrocytes play important roles in cognitive function and disease. Whether both cell types follow similar or distinctive evolutionary trajectories is not known. We examined the transcriptomes of neurons and oligodendrocytes in the frontal cortex of humans, chimpanzees, and rhesus macaques. We identified human-specific trajectories of gene expression in neurons and oligodendrocytes and show that both cell types exhibit human-specific up-regulation. Moreover, oligodendrocytes have undergone more pronounced accelerated gene expression evolution in the human lineage compared to neurons. We highlighted human-specific coexpression networks with specific functions. Our data suggest that oligodendrocyte human-specific networks are enriched for alternative splicing and transcriptional regulation. Oligodendrocyte networks are also enriched for variants associated with schizophrenia and other neuropsychiatric disorders. Such enrichments were not found in neuronal networks. These results offer a glimpse into the molecular mechanisms of oligodendrocytes during evolution and how such mechanisms are associated with neuropsychiatric disorders.

Human Genomic Signatures of Brain Oscillations During Memory Encoding.

Memory encoding is an essential step for all learning. However, the genetic and molecular mechanisms underlying human memory encoding remain poorly understood, and how this molecular framework permits the emergence of specific patterns of brain oscillations observed during mnemonic processing is unknown. Here, we directly compare intracranial electroencephalography recordings from the neocortex in individuals performing an episodic memory task with human gene expression from the same areas. We identify genes correlated with oscillatory memory effects across 6 frequency bands. These genes are enriched for autism-related genes and have preferential expression in neurons, in particular genes encoding synaptic proteins and ion channels, supporting the idea that the genes regulating voltage gradients are involved in the modulation of oscillatory patterns during successful memory encoding across brain areas. Memory-related genes are distinct from those correlated with other forms of cognitive processing and resting state fMRI. These data are the first to identify correlations between gene expression and active human brain states as well as provide a molecular window into memory encoding oscillations in the human brain.

ELAVL2-regulated transcriptional and splicing networks in human neurons link neurodevelopment and autism.

The role of post-transcriptional gene regulation in human brain development and neurodevelopmental disorders remains mostly uncharacterized. ELAV-like RNA-binding proteins (RNAbps) are a family of proteins that regulate several aspects of neuronal function including neuronal excitability and synaptic transmission, both critical to the normal function of the brain in cognition and behavior. Here, we identify the downstream neuronal transcriptional and splicing networks of ELAVL2, an RNAbp with previously unknown function in the brain. Expression of ELAVL2 was reduced in human neurons and RNA-sequencing was utilized to identify networks of differentially expressed and alternatively spliced genes resulting from haploinsufficient levels of ELAVL2. These networks contain a number of autism-relevant genes as well as previously identified targets of other important RNAbps implicated in autism spectrum disorder (ASD) including RBFOX1 and FMRP. ELAVL2-regulated co-expression networks are also enriched for neurodevelopmental and synaptic genes, and include genes with human-specific patterns of expression in the frontal pole. Together, these data suggest that ELAVL2 regulation of transcript expression is critical for neuronal function and clinically relevant to ASD.

MEF2C Hypofunction in GABAergic Cells Alters Sociability and Prefrontal Cortex Inhibitory Synaptic Transmission in a Sex-Dependent Manner.

Background: Heterozygous mutations or deletions of MEF2C cause a neurodevelopmental disorder termed MEF2C haploinsufficiency syndrome (MCHS), characterized by autism spectrum disorder and neurological symptoms. In mice, global Mef2c heterozygosity has produced multiple MCHS-like phenotypes. MEF2C is highly expressed in multiple cell types of the developing brain, including GABAergic (gamma-aminobutyric acidergic) inhibitory neurons, but the influence of MEF2C hypofunction in GABAergic neurons on MCHS-like phenotypes remains unclear. Methods: We employed GABAergic cell type-specific manipulations to study mouse Mef2c heterozygosity in a battery of MCHS-like behaviors. We also performed electroencephalography, single-cell transcriptomics, and patch-clamp electrophysiology and optogenetics to assess the impact of Mef2c haploinsufficiency on gene expression and prefrontal cortex microcircuits. Results: Mef2c heterozygosity in developing GABAergic cells produced female-specific deficits in social preference and altered approach-avoidance behavior. In female, but not male, mice, we observed that Mef2c heterozygosity in developing GABAergic cells produced 1) differentially expressed genes in multiple cell types, including parvalbumin-expressing GABAergic neurons, 2) baseline and social-related frontocortical network activity alterations, and 3) reductions in parvalbumin cell intrinsic excitability and inhibitory synaptic transmission onto deep-layer pyramidal neurons. Conclusions: MEF2C hypofunction in female, but not male, developing GABAergic cells is important for typical sociability and approach-avoidance behaviors and normal parvalbumin inhibitory neuron function in the prefrontal cortex of mice. While there is no apparent sex bias in autism spectrum disorder symptoms of MCHS, our findings suggest that GABAergic cell-specific dysfunction in females with MCHS may contribute disproportionately to sociability symptoms.

High affinity chimeric antigen receptor signaling induces an inflammatory program in human regulatory T cells.

Regulatory T cells (Tregs) are promising cellular therapies to induce immune tolerance in organ transplantation and autoimmune disease. The success of chimeric antigen receptor (CAR) T-cell therapy for cancer has sparked interest in using CARs to generate antigen-specific Tregs. Here, we compared CAR with endogenous T cell receptor (TCR)/CD28 activation in human Tregs. Strikingly, CAR Tregs displayed increased cytotoxicity and diminished suppression of antigen-presenting cells and effector T (Teff) cells compared with TCR/CD28 activated Tregs. RNA sequencing revealed that CAR Tregs activate Teff cell gene programs. Indeed, CAR Tregs secreted high levels of inflammatory cytokines, with a subset of FOXP3+ CAR Tregs uniquely acquiring CD40L surface expression and producing IFNγ. Interestingly, decreasing CAR antigen affinity reduced Teff cell gene expression and inflammatory cytokine production by CAR Tregs. Our findings showcase the impact of engineered receptor activation on Treg biology and support tailoring CAR constructs to Tregs for maximal therapeutic efficacy.

Histone deacetylase 5 in prelimbic prefrontal cortex limits context-associated cocaine seeking.

Background: Repeated cocaine use produces neuroadaptations that support drug craving and relapse in substance use disorders (SUDs). Powerful associations formed with drug-use environments can promote a return to active drug use in SUD patients, but the molecular mechanisms that control the formation of these prepotent drug-context associations remain unclear. Methods: In the rat intravenous cocaine self-administration (SA) model, we examined the role and regulation of histone deacetylase 5 (HDAC5) in the prelimbic (PrL) and infralimbic (IL) cortices in context-associated drug seeking. To this end, we employed viral molecular tools, chemogenetics, RNA-sequencing, electrophysiology, and immunohistochemistry. Results: In the PrL, reduction of endogenous HDAC5 augmented context-associated, but not cue-or drug prime-reinstated cocaine seeking, whereas overexpression of HDAC5 in PrL, but not IL, reduced context-associated cocaine seeking, but it had no effects on sucrose seeking. In contrast, PrL HDAC5 overexpression following acquisition of cocaine SA had no effects on future cocaine seeking. We found that HDAC5 and cocaine SA altered the expression of numerous PrL genes, including many synapse-associated genes. HDAC5 significantly increased inhibitory synaptic transmission onto PrL deep-layer pyramidal neurons, and it reduced the induction of FOS-positive neurons in the cocaine SA environment. Conclusions: Our findings reveal an essential and selective role for PrL HDAC5 to limit associations formed in cocaine, but not sucrose, SA environments, and that it alters the PrL excitatory/inhibitory balance, possibly through epigenetic regulation of synaptic genes. These results further position HDAC5 as a key factor regulating reward-circuit neuroadaptations that underlie common relapse triggers in SUD.

NPAS4 supports cocaine-conditioned cues in rodents by controlling the cell type-specific activation balance in the nucleus accumbens.

Powerful associations that link drugs of abuse with cues in the drug-paired environment often serve as prepotent relapse triggers. Drug-associated contexts and cues activate ensembles of nucleus accumbens (NAc) neurons, including D1-class medium spiny neurons (MSNs) that typically promote, and D2-class MSNs that typically oppose, drug seeking. We found that in mice, cocaine conditioning upregulated transiently the activity-regulated transcription factor, Neuronal PAS Domain Protein 4 (NPAS4), in a small subset of NAc neurons. The NPAS4+ NAc ensemble was required for cocaine conditioned place preference. We also observed that NPAS4 functions within NAc D2-, but not D1-, MSNs to support cocaine-context associations and cue-induced cocaine, but not sucrose, seeking. Together, our data show that the NPAS4+ ensemble of NAc neurons is essential for cocaine-context associations in mice, and that NPAS4 itself functions in NAc D2-MSNs to support cocaine-context associations by suppressing drug-induced counteradaptations that oppose relapse-related behaviour.