Convergence of case-specific epigenetic alterations identify a confluence of genetic vulnerabilities tied to opioid overdose.

Opioid use disorder is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. Opioid overdose, the most severe outcome of opioid use disorder, remains the leading cause of accidental death in the United States. We interrogated the effects of opioid overdose on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicted case-control status with high accuracy. We focused on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or "plexus": ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid use disorder and captures genetic and environmental factors perpetuating the opioid epidemic.

Blood miR-144-3p: a novel diagnostic and therapeutic tool for depression.

Major depressive disorder (MDD) is the leading cause of disability worldwide. There is an urgent need for objective biomarkers to diagnose this highly heterogeneous syndrome, assign treatment, and evaluate treatment response and prognosis. MicroRNAs (miRNAs) are short non-coding RNAs, which are detected in body fluids that have emerged as potential biomarkers of many disease conditions. The present study explored the potential use of miRNAs as biomarkers for MDD and its treatment. We profiled the expression levels of circulating blood miRNAs from mice that were collected before and after exposure to chronic social defeat stress (CSDS), an extensively validated mouse model used to study depression, as well as after either repeated imipramine or single-dose ketamine treatment. We observed robust differences in blood miRNA signatures between stress-resilient and stress-susceptible mice after an incubation period, but not immediately after exposure to the stress. Furthermore, ketamine treatment was more effective than imipramine at re-establishing baseline miRNA expression levels, but only in mice that responded behaviorally to the drug. We identified the red blood cell-specific miR-144-3p as a candidate biomarker to aid depression diagnosis and predict ketamine treatment response in stress-susceptible mice and MDD patients. Lastly, we demonstrate that systemic knockdown of miR-144-3p, via subcutaneous administration of a specific antagomir, is sufficient to reduce the depression-related phenotype in stress-susceptible mice. RNA-sequencing analysis of blood after such miR-144-3p knockdown revealed a blunted transcriptional stress signature as well. These findings identify miR-144-3p as a novel target for diagnosis of MDD as well as for antidepressant treatment, and enhance our understanding of epigenetic processes associated with depression.

Sperm transcriptional state associated with paternal transmission of stress phenotypes.

Paternal stress can induce long-lasting changes in germ cells potentially propagating heritable changes across generations. To date, no studies have investigated differences in transmission patterns between stress-resilient and -susceptible mice. We tested the hypothesis that transcriptional alterations in sperm during chronic social defeat stress (CSDS) transmit increased susceptibility to stress phenotypes to the next generation. We demonstrate differences in offspring from stressed fathers that depend upon paternal category (resilient vs susceptible) and offspring sex. Importantly, artificial insemination reveals that sperm mediates some of the behavioral phenotypes seen in offspring. Using RNA-sequencing we report substantial and distinct changes in the transcriptomic profiles of sperm following CSDS in susceptible vs resilient fathers, with alterations in long noncoding RNAs (lncRNAs) predominating especially in susceptibility. Correlation analysis revealed that these alterations were accompanied by a loss of regulation of protein-coding genes by lncRNAs in sperm of susceptible males. We also identify several co-expression gene modules that are enriched in differentially expressed genes in sperm from either resilient or susceptible fathers. Taken together, these studies advance our understanding of intergenerational epigenetic transmission of behavioral experience.SIGNIFICANCE STATEMENTThis manuscript contributes to the complex factors that influence the paternal transmission of stress phenotypes. By leveraging the segregation of males exposed to chronic social defeat stress into either resilient or susceptible categories we were able to identify the phenotypic differences in the paternal transmission of stress phenotypes across generations between the two lineages. Importantly, this work also alludes to the significance of both long noncoding RNAs and protein coding genes mediating the paternal transmission of stress. The knowledge gained from these data is of particular interest in understanding the risk for the development of psychiatric disorders such as anxiety and depression.

Complement pathway changes at age 12 are associated with psychotic experiences at age 18 in a longitudinal population-based study: evidence for a role of stress.

The complement cascade is a major component of the immune defence against infection, and there is increasing evidence for a role of dysregulated complement in major psychiatric disorders. We undertook a directed proteomic analysis of the complement signalling pathway (n = 29 proteins) using data-independent acquisition. Participants were recruited from the UK avon longitudinal study of parents and children (ALSPAC) cohort who participated in psychiatric assessment interviews at ages 12 and 18. Protein expression levels at age 12 among individuals who reported psychotic experiences (PEs) at age 18 (n = 64) were compared with age-matched controls (n = 67). Six out of the 29 targeted complement proteins or protein subcomponents were significantly upregulated following correction for multiple comparisons (VTN↑, C1RL↑, C8B↑, C8A↑, CFH↑, and C5↑). We then undertook an unbiased plasma proteomic analysis of mice exposed to chronic social stress and observed dysregulation of 11 complement proteins, including three that were altered in the same direction in individuals with PE (C1R↑, CFH↑, and C5↑). Our findings indicate that dysregulation of the complement protein pathway in blood is associated with incidence of psychotic experiences and that these changes may reflect exposure to stress.

Aberrant H3.3 dynamics in NAc promote vulnerability to depressive-like behavior.

Human major depressive disorder (MDD), along with related mood disorders, is among the world's greatest public health concerns; however, its pathophysiology remains poorly understood. Persistent changes in gene expression are known to promote physiological aberrations implicated in MDD. More recently, histone mechanisms affecting cell type- and regional-specific chromatin structures have also been shown to contribute to transcriptional programs related to depressive behaviors, as well as responses to antidepressants. Although much emphasis has been placed in recent years on roles for histone posttranslational modifications and chromatin-remodeling events in the etiology of MDD, it has become increasingly clear that replication-independent histone variants (e.g., H3.3), which differ in primary amino acid sequence from their canonical counterparts, similarly play critical roles in the regulation of activity-dependent neuronal transcription, synaptic connectivity, and behavioral plasticity. Here, we demonstrate a role for increased H3.3 dynamics in the nucleus accumbens (NAc)-a key limbic brain reward region-in the regulation of aberrant social stress-mediated gene expression and the precipitation of depressive-like behaviors in mice. We find that molecular blockade of these dynamics promotes resilience to chronic social stress and results in a partial renormalization of stress-associated transcriptional patterns in the NAc. In sum, our findings establish H3.3 dynamics as a critical, and previously undocumented, regulator of mood and suggest that future therapies aimed at modulating striatal histone dynamics may potentiate beneficial behavioral adaptations to negative emotional stimuli.

Histone arginine methylation in cocaine action in the nucleus accumbens.

Repeated cocaine exposure regulates transcriptional regulation within the nucleus accumbens (NAc), and epigenetic mechanisms-such as histone acetylation and methylation on Lys residues-have been linked to these lasting actions of cocaine. In contrast to Lys methylation, the role of histone Arg (R) methylation remains underexplored in addiction models. Here we show that protein-R-methyltransferase-6 (PRMT6) and its associated histone mark, asymmetric dimethylation of R2 on histone H3 (H3R2me2a), are decreased in the NAc of mice and rats after repeated cocaine exposure, including self-administration, and in the NAc of cocaine-addicted humans. Such PRMT6 down-regulation occurs selectively in NAc medium spiny neurons (MSNs) expressing dopamine D2 receptors (D2-MSNs), with opposite regulation occurring in D1-MSNs, and serves to protect against cocaine-induced addictive-like behavioral abnormalities. Using ChIP-seq, we identified Src kinase signaling inhibitor 1 (Srcin1; also referred to as p140Cap) as a key gene target for reduced H3R2me2a binding, and found that consequent Srcin1 induction in the NAc decreases Src signaling, cocaine reward, and the motivation to self-administer cocaine. Taken together, these findings suggest that suppression of Src signaling in NAc D2-MSNs, via PRMT6 and H3R2me2a down-regulation, functions as a homeostatic brake to restrain cocaine action, and provide novel candidates for the development of treatments for cocaine addiction.

Incubation of methamphetamine craving is associated with selective increases in expression of Bdnf and trkb, glutamate receptors, and epigenetic enzymes in cue-activated fos-expressing dorsal striatal neurons.

Cue-induced methamphetamine seeking progressively increases after withdrawal (incubation of methamphetamine craving), but the underlying mechanisms are largely unknown. We determined whether this incubation is associated with alterations in candidate genes in dorsal striatum (DS), a brain area implicated in cue- and context-induced drug relapse. We first measured mRNA expression of 24 candidate genes in whole DS extracts after short (2 d) or prolonged (1 month) withdrawal in rats following extended-access methamphetamine or saline (control condition) self-administration (9 h/d, 10 d). We found minimal changes. Next, using fluorescence-activated cell sorting, we compared gene expression in Fos-positive dorsal striatal neurons, which were activated during "incubated" cue-induced drug-seeking tests after prolonged withdrawal, with nonactivated Fos-negative neurons. We found significant increases in mRNA expression of immediate early genes (Arc, Egr1), Bdnf and its receptor (Trkb), glutamate receptor subunits (Gria1, Gria3, Grm1), and epigenetic enzymes (Hdac3, Hdac4, Hdac5, GLP, Dnmt3a, Kdm1a) in the Fos-positive neurons only. Using RNAscope to determine striatal subregion and cell-type specificity of the activated neurons, we measured colabeling of Fos with Drd1 and Drd2 in three DS subregions. Fos expression was neither subregion nor cell-type specific (52.5 and 39.2% of Fos expression colabeled with Drd1 and Drd2, respectively). Finally, we found that DS injections of SCH23390 (C17H18ClNO), a D1-family receptor antagonist known to block cue-induced Fos induction, decreased incubated cue-induced methamphetamine seeking after prolonged withdrawal. Results demonstrate a critical role of DS in incubation of methamphetamine craving and that this incubation is associated with selective gene-expression alterations in cue-activated D1- and D2-expressing DS neurons.

Threonine 149 phosphorylation enhances ΔFosB transcriptional activity to control psychomotor responses to cocaine.

Stable changes in neuronal gene expression have been studied as mediators of addicted states. Of particular interest is the transcription factor ΔFosB, a truncated and stable FosB gene product whose expression in nucleus accumbens (NAc), a key reward region, is induced by chronic exposure to virtually all drugs of abuse and regulates their psychomotor and rewarding effects. Phosphorylation at Ser(27) contributes to ΔFosB's stability and accumulation following repeated exposure to drugs, and our recent work demonstrates that the protein kinase CaMKIIα phosphorylates ΔFosB at Ser(27) and regulates its stability in vivo. Here, we identify two additional sites on ΔFosB that are phosphorylated in vitro by CaMKIIα, Thr(149) and Thr(180), and demonstrate their regulation in vivo by chronic cocaine. We show that phosphomimetic mutation of Thr(149) (T149D) dramatically increases AP-1 transcriptional activity while alanine mutation does not affect transcriptional activity when compared with wild-type (WT) ΔFosB. Using in vivo viral-mediated gene transfer of ΔFosB-T149D or ΔFosB-T149A in mouse NAc, we determined that overexpression of ΔFosB-T149D in NAc leads to greater locomotor activity in response to an initial low dose of cocaine than does WT ΔFosB, while overexpression of ΔFosB-T149A does not produce the psychomotor sensitization to chronic low-dose cocaine seen after overexpression of WT ΔFosB and abrogates the sensitization seen in control animals at higher cocaine doses. We further demonstrate that mutation of Thr(149) does not affect the stability of ΔFosB overexpressed in mouse NAc, suggesting that the behavioral effects of these mutations are driven by their altered transcriptional properties.

The Epigenetics of Anxiety Pathophysiology: A DNA Methylation and Histone Modification Focused Review.

Anxiety is one of the most common psychiatric disorders diagnosed in the United States today. Like all mental illnesses, anxiety pathology includes genetic, molecular, somatic, and behavioral characteristics. Specific brain regions implicated in anxiety include the prefrontal cortex, amygdala, hippocampus, and hypothalamus. Together, these regions regulate fear-related learning and memory processes, and are innervated by neuronal projections that use glutamate and GABA as neurotransmitters. Neurotrophic factors such as brain-derived neurotrophic factor (BDNF) are also implicated in anxiety. This review discusses the neuroepigenetics of the anxiety phenotype. While studying such changes is limited to postmortem brain studies or peripheral tissue acquisition in humans, the use of animals to model anxiety phenotypes has made epigenetic research possible. In this review, we summarize and discuss a plethora of DNA methylation, histone modification, and associated gene expression differences underscoring the anxiety phenotype. The findings we outline include expression changes of various DNA methyltransferases and changes in histone modifications that affect the hypothalamic pituitary adrenal axis and stress response as well as GABA, glutamate, and BDNF signaling in the PFC, amygdala, hypothalamus, and hippocampus. Furthermore, there have been studies showing that anxiety behaviors and biological scars from stress can be reversed using histone deacetylase inhibitors, and we discuss ideas for the future of treatment. In this review, we hope that by compiling much of the data pertaining to DNA methylation and histone modifications in vivo animal studies we are able to highlight potential avenues for future research despite existing limitations.

Sex-Specific Transcriptional Changes in Response to Adolescent Social Stress in the Brain's Reward Circuitry.

BACKGROUND: Sex differences in addiction have been described in humans and animal models. A key factor that influences addiction in both males and females is adolescent experience. Adolescence is associated with higher vulnerability to substance use disorders, and male rodents subjected to adolescent social isolation (SI) stress form stronger preferences for drugs of abuse in adulthood. However, little is known about how females respond to SI, and few studies have investigated the transcriptional changes induced by SI in the brain's reward circuitry. METHODS: We tested the hypothesis that SI alters the transcriptome in a persistent and sex-specific manner in prefrontal cortex, nucleus accumbens, and ventral tegmental area. Mice were isolated or group housed from postnatal day P22 to P42, then group housed until ∼P90. Transcriptome-wide changes were investigated by RNA sequencing after acute or chronic cocaine or saline administration. RESULTS: We found that SI disrupts sex-specific transcriptional responses to cocaine and reduces sex differences in gene expression across all three brain regions. Furthermore, SI induces gene expression profiles in males that more closely resemble group-housed females, suggesting that SI "feminizes" the male transcriptome. Coexpression analysis reveals that such disruption of sex differences in gene expression alters sex-specific gene networks and identifies potential sex-specific key drivers of these transcriptional changes. CONCLUSIONS: Together, these data show that SI has region-specific effects on sex-specific transcriptional responses to cocaine and provide a better understanding of reward-associated transcription that differs in males and females.

Crystallin Mu in Medial Amygdala Mediates the Effect of Social Experience on Cocaine Seeking in Males but Not in Females.

BACKGROUND: Social experiences influence susceptibility to substance use disorder. The adolescent period is associated with the development of social reward and is exceptionally sensitive to disruptions to reward-associated behaviors by social experiences. Social isolation (SI) during adolescence alters anxiety- and reward-related behaviors in adult males, but little is known about females. The medial amygdala (meA) is a likely candidate for the modulation of social influence on drug reward because it regulates social reward, develops during adolescence, and is sensitive to social stress. However, little is known regarding how the meA responds to drugs of abuse. METHODS: We used adolescent SI coupled with RNA sequencing to better understand the molecular mechanisms underlying meA regulation of social influence on reward. RESULTS: We show that SI in adolescence, a well-established preclinical model for addiction susceptibility, enhances preference for cocaine in male but not in female mice and alters cocaine-induced protein and transcriptional profiles within the adult meA particularly in males. To determine whether transcriptional mechanisms within the meA are important for these behavioral effects, we manipulated Crym expression, a sex-specific key driver gene identified through differential gene expression and coexpression network analyses, specifically in meA neurons. Overexpression of Crym, but not another key driver that did not meet our sex-specific criteria, recapitulated the behavioral and transcriptional effects of adolescent SI. CONCLUSIONS: These results show that the meA is essential for modulating the sex-specific effects of social experience on drug reward and establish Crym as a critical mediator of sex-specific behavioral and transcriptional plasticity.

Silent synapses dictate cocaine memory destabilization and reconsolidation.

Cocaine-associated memories are persistent, but, on retrieval, become temporarily destabilized and vulnerable to disruptions, followed by reconsolidation. To explore the synaptic underpinnings for these memory dynamics, we studied AMPA receptor (AMPAR)-silent excitatory synapses, which are generated in the nucleus accumbens by cocaine self-administration, and subsequently mature after prolonged withdrawal by recruiting AMPARs, echoing acquisition and consolidation of cocaine memories. We show that, on memory retrieval after prolonged withdrawal, the matured silent synapses become AMPAR-silent again, followed by re-maturation ~6 h later, defining the onset and termination of a destabilization window of cocaine memories. These synaptic dynamics are timed by Rac1, with decreased and increased Rac1 activities opening and closing, respectively, the silent synapse-mediated destabilization window. Preventing silent synapse re-maturation within the destabilization window decreases cue-induced cocaine seeking. Thus, cocaine-generated silent synapses constitute a discrete synaptic ensemble dictating the dynamics of cocaine-associated memories and can be targeted for memory disruption.

Fosb Induction in Nucleus Accumbens by Cocaine Is Regulated by E2F3a.

The transcription factor ΔFosB has been proposed as a molecular switch for the transition from casual, volitional drug use into a chronically addicted state, but the upstream regulatory mechanisms governing ΔFosB expression are incompletely understood. In this study, we find a novel regulatory role for the transcription factor E2F3, recently implicated in transcriptional regulation by cocaine, in controlling ΔFosB induction in the mouse nucleus accumbens (NAc) following cocaine administration. We find that an E2F consensus sequence 500 bp upstream of the Fosb transcription start site is enriched for E2F3 specifically over other E2F isoforms. We further conclude that ΔFosB expression is regulated specifically by E2F3a, not E2F3b, that E2f3a expression is specific to D1 receptor-expressing medium spiny neurons, and that E2F3a overexpression in NAc recapitulates the induction of Fosb and ΔFosb mRNA expression observed after chronic cocaine exposure. E2F3a knockdown in NAc does not abolish ΔFosb induction by cocaine, a result consistent with previously published data showing that singular knockdown of upstream regulators of ΔFosB is insufficient to block cocaine-induced expression. Finally, to elucidate potential combinatorial epigenetic mechanisms involved in E2F3a's regulation of Fosb, we explore H3K4me3 enrichment at the Fosb promoter and find that it is not enhanced by E2F3a overexpression, suggesting that it may instead be a pre-existing permissive mark allowing for E2F3a to interact with Fosb. Together, these findings support a role for E2F3a as a novel, upstream regulator of the addiction-mediating transcription factor ΔFosB in NAc.

National Institute on Drug Abuse genomics consortium white paper: Coordinating efforts between human and animal addiction studies.

The National Institute on Drug Abuse Genetics and Epigenetics Cross-Cutting Research Team convened a diverse group of researchers, clinicians, and healthcare providers on the campus of the University of California, San Diego, in June 2018. The goal was to develop strategies to integrate genetics and phenotypes across species to achieve a better understanding of substance use disorders through associations between genotypes and addictive behaviors. This conference (a) discussed progress in harmonizing large opioid genetics cohorts, (b) discussed phenotypes that are used for genetics studies in humans, (c) examined phenotypes that are used for genetics studies in animal models, (d) identified synergies and gaps in phenotypic analyses of human and animal models and (e) identified strategies to integrate genetics and genomics data with phenotypes across species. The meeting consisted of panels that focused on phenotype harmonization (Dr. Laura Bierut, Dr. Olivier George, Dr. Dan Larach and Dr. Sesh Mudumbai), translating genetic findings between species (Dr. Elissa Chesler, Dr. Gary Peltz and Dr. Abraham Palmer), interpreting and understanding allelic variations (Dr. Vanessa Troiani and Dr. Tamara Richards) and pathway conservation in animal models and human studies (Dr. Robert Hitzemann, Dr. Huda Akil and Dr. Laura Saba). There were also updates that were provided by large consortia (Dr. Susan Tapert, Dr. Danielle Dick, Dr. Howard Edenberg and Dr. Eric Johnson). Collectively, the conference was convened to discuss progress and changes in genome-wide association studies.

A novel role for E2F3b in regulating cocaine action in the prefrontal cortex.

Drug abuse is a multifaceted disorder that involves maladaptive decision making. Long-lasting changes in the addicted brain are mediated by a complex circuit of brain reward regions. The prefrontal cortex (PFC) is one region in which chronic drug exposure changes expression and function of upstream transcriptional regulators to alter drug responses and aspects of the addicted phenotype. We reported recently that the transcription factor E2F3a is a critical mediator of cocaine responses in the nucleus accumbens. E2F3a is one of two splice variants of the E2f3 gene; the other is E2F3b. Another recent study predicted E2F3 as an upstream regulator of the transcriptional response to cocaine self-administration (SA) in PFC. Based on previous findings that E2F3a and E2F3b have divergent regulatory roles, we set out to study the putative transcriptional role of these transcripts in PFC in the context of repeated I.P. cocaine exposure. We implemented viral-mediated isoform-specific gene manipulation, RNA-sequencing, advanced bioinformatics analyses, and animal behavior to determine how E2F3a and E2F3b contribute to persistent cocaine-induced transcriptional changes in PFC. We show that E2F3b, but not E2F3a, in PFC is critical for cocaine locomotor and place preference behaviors. Interestingly, RNA-seq of PFC following E2f3b overexpression or I.P. cocaine exposure showed very different effects on expression levels of differentially expressed genes. However, we found that E2F3b drives a similar transcriptomic pattern to that of cocaine SA with overlapping upstream regulators and downstream pathways predicted. These findings reveal a novel transcriptional mechanism in PFC that controls behavioral and molecular responses to cocaine.

Early life stress alters transcriptomic patterning across reward circuitry in male and female mice.

Abuse, neglect, and other forms of early life stress (ELS) significantly increase risk for psychiatric disorders including depression. In this study, we show that ELS in a postnatal sensitive period increases sensitivity to adult stress in female mice, consistent with our earlier findings in male mice. We used RNA-sequencing in the ventral tegmental area, nucleus accumbens, and prefrontal cortex of male and female mice to show that adult stress is distinctly represented in the brain's transcriptome depending on ELS history. We identify: 1) biological pathways disrupted after ELS and associated with increased behavioral stress sensitivity, 2) putative transcriptional regulators of the effect of ELS on adult stress response, and 3) subsets of primed genes specifically associated with latent behavioral changes. We also provide transcriptomic evidence that ELS increases sensitivity to future stress through enhancement of known programs of cortical plasticity.

Genome-wide transcriptional profiling of central amygdala and orbitofrontal cortex during incubation of methamphetamine craving.

Methamphetamine (Meth) seeking progressively increases after withdrawal (incubation of Meth craving), but the transcriptional mechanisms that contribute to this incubation are unknown. Here we used RNA-sequencing to analyze transcriptional profiles associated with incubation of Meth craving in central amygdala (CeA) and orbitofrontal cortex (OFC), two brain areas previously implicated in relapse to drug seeking. We trained rats to self-administer either saline (control condition) or Meth (10 days; 9 h/day, 0.1 mg/kg/infusion). Next, we collected brain tissue from CeA and OFC on withdrawal day 2 (when Meth seeking is low and non-incubated) and on day 35 (when Meth seeking is high and incubated), for subsequent RNA-sequencing. In CeA, we identified 10-fold more differentially expressed genes (DEGs) on withdrawal day 35 than day 2. These genes were enriched for several biological processes, including protein ubiquitination and histone methylation. In OFC, we identified much fewer expression changes than in CeA, with more DEGs on withdrawal day 2 than on day 35. There was a significant overlap between upregulated genes on withdrawal day 2 and downregulated genes on withdrawal day 35 in OFC. Our analyses highlight the CeA as a key region of transcriptional regulation associated with incubation of Meth seeking. In contrast, transcriptional regulation in OFC may contribute to Meth seeking during early withdrawal. Overall, these findings provide a unique resource of gene expression data for future studies examining transcriptional mechanisms in CeA that mediate Meth seeking after prolonged withdrawal.

Transcription Factor E2F3a in Nucleus Accumbens Affects Cocaine Action via Transcription and Alternative Splicing.

BACKGROUND: Lasting changes in gene expression in brain reward regions, including nucleus accumbens (NAc), contribute to persistent functional changes in the addicted brain. We and others have demonstrated that altered expression of several candidate transcription factors in NAc regulates drug responses. A recent large-scale genome-wide study from our group predicted transcription factor E2F3 (E2F3) as a prominent upstream regulator of cocaine-induced changes in gene expression and alternative splicing. METHODS: We studied expression of two E2F3 isoforms-E2F3a and E2F3b-in mouse NAc after repeated cocaine administration and assayed the effects of overexpression or depletion of E2f3 isoforms in NAc on cocaine behavioral responses. We then performed RNA sequencing to investigate the effect of E2f3a overexpression in this region on gene expression and alternative splicing and performed quantitative chromatin immunoprecipitation at downstream targets in NAc following E2f3a overexpression or repeated cocaine exposure. Sample sizes varied between experiments and are noted in the text. RESULTS: We showed that E2f3a, but not E2f3b, overexpression or knockdown in mouse NAc regulates cocaine-induced locomotor and place conditioning behavior. Furthermore, we demonstrated that E2f3a overexpression substantially recapitulates genome-wide transcriptional profiles and alternative splicing induced by cocaine. We further validated direct binding of E2F3a at key target genes following cocaine exposure. CONCLUSIONS: This study establishes E2F3a as a novel transcriptional regulator of cocaine action in NAc. The findings reveal a crucial role for E2F3a in the regulation of cocaine-elicited behavioral states. Moreover, the importance of this role is bolstered by the extensive recapitulation of cocaine's transcriptional effects in NAc by overexpression of E2f3a.

Cocaine Self-administration Alters Transcriptome-wide Responses in the Brain's Reward Circuitry.

BACKGROUND: Global changes in gene expression underlying circuit and behavioral dysregulation associated with cocaine addiction remain incompletely understood. Here, we show how a history of cocaine self-administration (SA) reprograms transcriptome-wide responses throughout the brain's reward circuitry at baseline and in response to context and/or cocaine re-exposure after prolonged withdrawal (WD). METHODS: We assigned male mice to one of six groups: saline/cocaine SA + 24-hour WD or saline/cocaine SA + 30-day WD + an acute saline/cocaine challenge within the previous drug-paired context. RNA sequencing was conducted on six interconnected brain reward regions. Using pattern analysis of gene expression and factor analysis of behavior, we identified genes that are strongly associated with addiction-related behaviors and uniquely altered by a history of cocaine SA. We then identified potential upstream regulators of these genes. RESULTS: We focused on three patterns of gene expression that reflect responses to 1) acute cocaine, 2) context re-exposure, and 3) drug + context re-exposure. These patterns revealed region-specific regulation of gene expression. Further analysis revealed that each of these gene expression patterns correlated with an addiction index-a composite score of several addiction-like behaviors during cocaine SA-in a region-specific manner. Cyclic adenosine monophosphate response element binding protein and nuclear receptor families were identified as key upstream regulators of genes associated with such behaviors. CONCLUSIONS: This comprehensive picture of transcriptome-wide regulation in the brain's reward circuitry by cocaine SA and prolonged WD provides new insight into the molecular basis of cocaine addiction, which will guide future studies of the key molecular pathways involved.

Ketamine and Imipramine Reverse Transcriptional Signatures of Susceptibility and Induce Resilience-Specific Gene Expression Profiles.

BACKGROUND: Examining transcriptional regulation by antidepressants in key neural circuits implicated in depression and understanding the relation to transcriptional mechanisms of susceptibility and natural resilience may help in the search for new therapeutic agents. Given the heterogeneity of treatment response in human populations, examining both treatment response and nonresponse is critical. METHODS: We compared the effects of a conventional monoamine-based tricyclic antidepressant, imipramine, and a rapidly acting, non-monoamine-based antidepressant, ketamine, in mice subjected to chronic social defeat stress, a validated depression model, and used RNA sequencing to analyze transcriptional profiles associated with susceptibility, resilience, and antidepressant response and nonresponse in the prefrontal cortex (PFC), nucleus accumbens, hippocampus, and amygdala. RESULTS: We identified similar numbers of responders and nonresponders after ketamine or imipramine treatment. Ketamine induced more expression changes in the hippocampus; imipramine induced more expression changes in the nucleus accumbens and amygdala. Transcriptional profiles in treatment responders were most similar in the PFC. Nonresponse reflected both the lack of response-associated gene expression changes and unique gene regulation. In responders, both drugs reversed susceptibility-associated transcriptional changes and induced resilience-associated transcription in the PFC. CONCLUSIONS: We generated a uniquely large resource of gene expression data in four interconnected limbic brain regions implicated in depression and its treatment with imipramine or ketamine. Our analyses highlight the PFC as a key site of common transcriptional regulation by antidepressant drugs and in both reversing susceptibility- and inducing resilience-associated molecular adaptations. In addition, we found region-specific effects of each drug, suggesting both common and unique effects of imipramine versus ketamine.