Data availability
All data associated with this project have been uploaded to Figshare (https://figshare.com/authors/David_Marcus/23255861).
Code availability
All code associated with this project has been uploaded to Figshare (https://figshare.com/authors/David_Marcus/23255861).
References
Devane, W. A., Dysarz, F. A., Johnson, M. R., Melvin, L. S. & Howlett, A. C. Determination and characterization of a cannabinoid receptor in rat brain. Mol. Pharmacol. 34, 605–613 (1988).
Article CAS PubMed Google Scholar
Herkenham, M. et al. Cannabinoid receptor localization in brain. Proc. Natl Acad. Sci. USA 87, 1932–1936 (1990).
Article ADS CAS PubMed Central PubMed Google Scholar
Devane, W. A. et al. Isolation and structure of a brain constituent that binds to the cannabinoid receptor. Science 258, 1946–1949 (1992).
Article ADS CAS PubMed Google Scholar
Stella, N., Schweitzer, P. & Piomelli, D. A second endogenous cannabinoid that modulates long-term potentiation. Nature 388, 773–778 (1997).
Article ADS CAS PubMed Google Scholar
Munro, S., Thomas, K. L. & Abu-Shaar, M. Molecular characterization of a peripheral receptor for cannabinoids. Nature 365, 61–65 (1993).
Article ADS CAS PubMed Google Scholar
Wilson, R. I. & Nicoll, R. A. Endogenous cannabinoids mediate retrograde signalling at hippocampal synapses. Nature 410, 588–592 (2001).
Article ADS CAS PubMed Google Scholar
Kano, M., Ohno-Shosaku, T., Hashimotodani, Y., Uchigashima, M. & Watanabe, M. Endocannabinoid-mediated control of synaptic transmission. Physiol. Rev. 89, 309–380 (2009).
Article CAS PubMed Google Scholar
Kano, M. Control of synaptic function by endocannabinoid-mediated retrograde signaling. Proc. Jpn Acad. Ser. B Phys. Biol. Sci. 90, 235–250 (2014).
Article CAS PubMed Central PubMed Google Scholar
Kreitzer, A. C. & Regehr, W. G. Retrograde inhibition of presynaptic calcium influx by endogenous cannabinoids at excitatory synapses onto Purkinje cells. Neuron 29, 717–727 (2001).
Article CAS PubMed Google Scholar
Horne, S. J., Topp, T. E. & Quigley, L. Depression and the willingness to expend cognitive and physical effort for rewards: a systematic review. Clin. Psychol. Rev. 88, 102065 (2021).
Article PubMed Google Scholar
Volkow, N. D. & Morales, M. The brain on drugs: from reward to addiction. Cell 162, 712–725 (2015).
Article CAS PubMed Google Scholar
Robinson, T. E. & Berridge, K. C. The psychology and neurobiology of addiction: an incentive-sensitization view. Addiction 95, S91–S117 (2000).
Article PubMed Google Scholar
Robinson, T. E. & Berridge, K. C. The neural basis of drug craving: an incentive-sensitization theory of addiction. Brain Res. Rev. 18, 247–291 (1993).
Article CAS PubMed Google Scholar
Purcell, J. R. et al. A review of risky decision-making in psychosis-spectrum disorders. Clin. Psychol. Rev. 91, 102112 (2022).
Article PubMed Google Scholar
Hikida, T., Morita, M. & Macpherson, T. Neural mechanisms of the nucleus accumbens circuit in reward and aversive learning. Neurosci. Res. 108, 1–5 (2016).
Article PubMed Google Scholar
Scofield, M. D. et al. The nucleus accumbens: mechanisms of addiction across drug classes reflect the importance of glutamate homeostasis. Pharmacol. Rev. 68, 816–871 (2016).
Article CAS PubMed Central PubMed Google Scholar
Castro, D. C. & Bruchas, M. R. A motivational and neuropeptidergic hub: anatomical and functional diversity within the nucleus accumbens shell. Neuron 102, 529–552 (2019).
Article CAS PubMed Central PubMed Google Scholar
Al-Hasani, R. et al. Distinct subpopulations of nucleus accumbens dynorphin neurons drive aversion and reward. Neuron 87, 1063–1077 (2015).
Article CAS PubMed Central PubMed Google Scholar
Cooper, S., Robison, A. J. & Mazei-Robison, M. S. Reward circuitry in addiction. Neurotherapeutics 14, 687–697 (2017).
Article CAS PubMed Central PubMed Google Scholar
Nieh, E. H., Kim, S.-Y., Namburi, P. & Tye, K. M. Optogenetic dissection of neural circuits underlying emotional valence and motivated behaviors. Brain Res. 1511, 73–92 (2013).
Article CAS PubMed Google Scholar
Floresco, S. B. The nucleus accumbens: an interface between cognition, emotion, and action. Annu. Rev. Psychol. 66, 25–52 (2015).
Article PubMed Google Scholar
Richard, J. M., Castro, D. C., DiFeliceantonio, A. G., Robinson, M. J. F. & Berridge, K. C. Mapping brain circuits of reward and motivation: in the footsteps of Ann Kelley. Neurosci. Biobehav. Rev. 37, 1919–1931 (2013).
Article PubMed Google Scholar
Christoffel, D. J. et al. Input-specific modulation of murine nucleus accumbens differentially regulates hedonic feeding. Nat. Commun. 12, 2135 (2021).
Article ADS CAS PubMed Central PubMed Google Scholar
Reed, S. J. et al. Coordinated reductions in excitatory input to the nucleus accumbens underlie food consumption. Neuron 99, 1260–1273.e4 (2018).
Article CAS PubMed Google Scholar
Lafferty, C. K., Yang, A. K., Mendoza, J. A. & Britt, J. P. Nucleus accumbens cell type- and input-specific suppression of unproductive reward seeking. Cell Rep. 30, 3729–3742.e3 (2020).
Article CAS PubMed Google Scholar
Deroche, M. A., Lassalle, O., Castell, L., Valjent, E. & Manzoni, O. J. Cell-type- and endocannabinoid-specific synapse connectivity in the adult nucleus accumbens core. J. Neurosci. 40, 1028–1041 (2020).
Article PubMed Google Scholar
Folkes, O. M. et al. An endocannabinoid-regulated basolateral amygdala–nucleus accumbens circuit modulates sociability. J. Clin. Invest. 130, 1728–1742 (2020).
Article CAS PubMed Central PubMed Google Scholar
Mateo, Y. et al. Endocannabinoid actions on cortical terminals orchestrate local modulation of dopamine release in the nucleus accumbens. Neuron 96, 1112–1126.e5 (2017).
Article CAS PubMed Central PubMed Google Scholar
Wenzel, J. M. et al. Phasic dopamine signals in the nucleus accumbens that cause active avoidance require endocannabinoid mobilization in the midbrain. Curr. Biol. 28, 1392–1404.e5 (2018).
Article CAS PubMed Central PubMed Google Scholar
Kondev, V. et al. Synaptic and cellular endocannabinoid signaling mechanisms regulate stress-induced plasticity of nucleus accumbens somatostatin neurons. Proc. Natl Acad. Sci. USA 120, e2300585120 (2023).
Article CAS PubMed Central PubMed Google Scholar
Penzo, M. A. & Gao, C. The paraventricular nucleus of the thalamus: an integrative node underlying homeostatic behavior. Trends Neurosci. 44, 538–549 (2021).
Article CAS PubMed Central PubMed Google Scholar
Zhou, K. & Zhu, Y. The paraventricular thalamic nucleus: a key hub of neural circuits underlying drug addiction. Pharmacol. Res. 142, 70–76 (2019).
Article PubMed Google Scholar
Dong, A. et al. A fluorescent sensor for spatiotemporally resolved imaging of endocannabinoid dynamics in vivo. Nat. Biotech. 40, 787–798 (2022).
Beas, S. et al. Dissociable encoding of motivated behavior by parallel thalamo-striatal projections. Curr. Biol. 34, 1549–1560.e3 (2024).
Article CAS PubMed Central PubMed Google Scholar
Gao, C. et al. Two genetically, anatomically and functionally distinct cell types segregate across anteroposterior axis of paraventricular thalamus. Nat. Neurosci. 23, 217–228 (2020).
Article CAS PubMed Central PubMed Google Scholar
Do-Monte, F. H., Minier-Toribio, A., Quiñones-Laracuente, K., Medina-Colón, E. M. & Quirk, G. J. Thalamic regulation of sucrose seeking during unexpected reward omission. Neuron 94, 388–400.e4 (2017).
Article CAS PubMed Central PubMed Google Scholar
Paniccia, J. E. et al. Restoration of a paraventricular thalamo-accumbal behavioral suppression circuit prevents reinstatement of heroin seeking. Neuron https://doi.org/10.1016/j.neuron.2023.11.024 (2023).
Yao, Z. et al. A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain. Nature 624, 317–332 (2023).
Article ADS CAS PubMed Central PubMed Google Scholar
Lein, E. S. et al. Genome-wide atlas of gene expression in the adult mouse brain. Nature 445, 168–176 (2007).
Article ADS CAS PubMed Google Scholar
Petreanu, L., Huber, D., Sobczyk, A. & Svoboda, K. Channelrhodopsin-2-assisted circuit mapping of long-range callosal projections. Nat. Neurosci. 10, 663–668 (2007).
Article CAS PubMed Google Scholar
Gunduz-Cinar, O. et al. A cortico-amygdala neural substrate for endocannabinoid modulation of fear extinction. Neuron 111, 3053–3067.e10 (2023).
Article CAS PubMed Central PubMed Google Scholar
Hunker, A. C. et al. Conditional single vector CRISPR/SaCas9 viruses for efficient mutagenesis in the adult mouse nervous system. Cell Rep. 30, 4303–4316.e6 (2020).
Article CAS PubMed Central PubMed Google Scholar
Schiffmann, S. N. & Vanderhaeghen, J. J. Adenosine A2 receptors regulate the gene expression of striatopallidal and striatonigral neurons. J. Neurosci. 13, 1080–1087 (1993).
Article CAS PubMed Central PubMed Google Scholar
Pereira, T. D. et al. SLEAP: a deep learning system for multi-animal pose tracking. Nat. Methods 19, 486–495 (2022).
Article CAS PubMed Central PubMed Google Scholar
Lobo, M. K. et al. Cell type-specific loss of BDNF signaling mimics optogenetic control of cocaine reward. Science 330, 385–390 (2010).
Article ADS CAS PubMed Central PubMed Google Scholar
Guillaumin, M. C. C., Viskaitis, P., Bracey, E., Burdakov, D. & Peleg-Raibstein, D. Disentangling the role of NAc D1 and D2 cells in hedonic eating. Mol. Psychiatry 28, 3531–3547 (2023).
Article PubMed Central PubMed Google Scholar
Walle, R. et al. Nucleus accumbens D1- and D2-expressing neurons control the balance between feeding and activity-mediated energy expenditure. Nat. Commun. 15, 2543 (2024).
Article ADS CAS PubMed Central PubMed Google Scholar
Domingues, A. V. et al. Dynamic representation of appetitive and aversive stimuli in nucleus accumbens shell D1- and D2-medium spiny neurons. Nat. Commun. 16, 59 (2025).
Article ADS PubMed Central PubMed Google Scholar
Pedersen, C. E. et al. Medial accumbens shell spiny projection neurons encode relative reward preference. Preprint at bioRxiv https://doi.org/10.1101/2022.09.18.508426 (2024).
Zingg, B., Dong, H.-W., Tao, H. W. & Zhang, L. I. Application of AAV1 for anterograde transsynaptic circuit mapping and input-dependent neuronal cataloging. Curr. Protoc. 2, e339 (2022).
Article PubMed Central PubMed Google Scholar
Zingg, B., Peng, B., Huang, J., Tao, H. W. & Zhang, L. I. Synaptic specificity and application of anterograde transsynaptic AAV for probing neural circuitry. J. Neurosci. 40, 3250–3267 (2020).
Article CAS PubMed Central PubMed Google Scholar
Xiao, X. et al. A genetically defined compartmentalized striatal direct pathway for negative reinforcement. Cell 183, 211–227.e20 (2020).
Article CAS PubMed Central PubMed Google Scholar
Li, H. et al. Neurotensin orchestrates valence assignment in the amygdala. Nature 608, 586–592 (2022).
Article ADS CAS PubMed Central PubMed Google Scholar
Alsammani, A., Stacey, W. C. & Gliske, S. V. Estimation of circular statistics in the presence of measurement bias. IEEE J. Biomed. Health Inform. 28, 1089–1100 (2024).
Article PubMed Central PubMed Google Scholar
Arski, O. N. et al. Epilepsy disrupts hippocampal phase precision and impairs working memory. Epilepsia 63, 2583–2596 (2022).
Article CAS PubMed Google Scholar
Vinck, M., Battaglia, F. P., Womelsdorf, T. & Pennartz, C. Improved measures of phase-coupling between spikes and the local field potential. J. Comput. Neurosci. 33, 53–75 (2012).
Article MathSciNet PubMed Google Scholar
Valentino, R. J. & Volkow, N. D. Cannabis and cannabinoid signaling: research gaps and opportunities. J. Pharmacol. Exp. Ther. 391, 154–158 (2024).
Article CAS PubMed Central PubMed Google Scholar
Lutz, B. Neurobiology of cannabinoid receptor signaling. Dialogues Clin. Neurosci. 22, 207–222 (2020).
Article PubMed Central PubMed Google Scholar
Dudok, B. et al. Retrograde endocannabinoid signaling at inhibitory synapses in vivo. Science 383, 967–970 (2024).
Article ADS CAS PubMed Central PubMed Google Scholar
Zhao, Z. et al. Cannabinoids regulate an insula circuit controlling water intake. Curr. Biol. 34, 1918–1929.e5 (2024).
Article CAS PubMed Google Scholar
Planert, H., Berger, T. K. & Silberberg, G. Membrane properties of striatal direct and indirect pathway neurons in mouse and rat slices and their modulation by dopamine. PLoS ONE 8, e57054 (2013).
Article ADS CAS PubMed Central PubMed Google Scholar
Kawaguchi, Y. Physiological, morphological, and histochemical characterization of three classes of interneurons in rat neostriatum. J. Neurosci. 13, 4908–4923 (1993).
Article CAS PubMed Central PubMed Google Scholar
Parker, K. E. et al. A paranigral VTA nociceptin circuit that constrains motivation for reward. Cell 178, 653–671.e19 (2019).
Article CAS PubMed Central PubMed Google Scholar
Machado, A. S., Darmohray, D. M., Fayad, J., Marques, H. G. & Carey, M. R. A quantitative framework for whole-body coordination reveals specific deficits in freely walking ataxic mice. eLife 4, e07892 (2015).
Article PubMed Central PubMed Google Scholar
Zhou, P. et al. Efficient and accurate extraction of in vivo calcium signals from microendoscopic video data. eLife 7, e28728 (2018).
Article PubMed Central PubMed Google Scholar
Resendez, S. L. et al. Visualization of cortical, subcortical and deep brain neural circuit dynamics during naturalistic mammalian behavior with head-mounted microscopes and chronically implanted lenses. Nat. Protoc. 11, 566–597 (2016).
Article CAS PubMed Central PubMed Google Scholar
Sheintuch, L. et al. Tracking the same neurons across multiple days in Ca2+ imaging data. Cell Rep. 21, 1102–1115 (2017).
Article CAS PubMed Central PubMed Google Scholar
Download references
Acknowledgements
We thank the Molecular Genetics Resource Core for the Center in Neurobiology of Addiction, Pain, and Emotion and its director, S. Schattauer, for generating the GRABeCB2.0 and CRISPR viruses used in this study; A. Suko for laboratory management and organization; T. Hobbs, C. Pizzano and V. Lau for colony management; and the entire Bruchas laboratory as well as other members of the NAPE Center at the University of Washington for resources and critical feedback.
Funding
This work was supported by the National Institute on Drug Abuse: F32 DA054709 (to D.J.M.), K99/R00 DA059617 (to D.J.M.), R37 DA033396 (to M.R.B.) and R21s DA056816 and DA057186 (to M.R.B. and N.S.); the National Institute of Mental Health R01 MH112355 (to M.R.B.); and the National Center for Complementary and Integrative Health RO1 AT011524 (to B.B.L.). Further support was provided by the UW Addictions, Drug, and Alcohol Institute research grant (to D.J.M.) and the Scan Design Foundation Innovative Pain Research Grant (to D.J.M.).
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks Anna Beyeler, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 2-AG is released in the NAc following rewarding stimuli.
(a) Schematic for fiber photometry recordings of GRABeCB2.0 in the NAc during Pavlovian reward conditioning. (b) Fiber placements for GRABeCB2.0 cohort displayed in Fig. 1. (c) Photometry trace of cue-aligned GRABeCB2.0 signal in SR141716A (10 mg/kg) treated animals (N = 11). (d-f) Photometry trace of lick-aligned GRABeCB2.0 signal in vehicle, SR141716A (10 mg/kg), and DO34 (20 mg/kg) treated animals (N = 11). (g) Comparison of cue-aligned GRABeCB2.0 signal (0-26 s post cue) showing attenuation of GRABeCB2.0 signal (0-20 s post lick) in SR141716A (p = 0.0131) and DO34 (p = 0.0079) treated animals compared to vehicle (N = 11). (h) Comparison of lick-aligned (0-20 s) GRABeCB2.0 signal in SR141716A (p = 0.0343) and DO34 (p = 0.0469) treated animals compared to vehicle (N = 11). (i-k) Photometry traces of cue-aligned GRABeCB2.0 signal in vehicle, JZL184 (5 mg/kg), and URB597 3 mg/kg) treated animals (N = 8). (l) Comparison of cue-aligned (0-26 s) GRABeCB2.0 signal in JZL184 (p = 0.0914) and URB597 (p = 0.6983) treated animals compared to vehicle (N = 8) (m) Schematic for Pavlovian reward conditioning with different cue modalities (n,o) Photometry traces of light and tone cue aligned GRABeCB2.0 signal (N = 7). (p) Comparison of light and tone aligned GRABeCB2.0 signal (N = 7). (q) Comparison of number of licks in light and tone Pavlovian reward conditioning sessions (N = 7). (r) Schematic for reward omission assay. (s,t) Photometry trace and quantification of cue-aligned (0-26 s) GRABeCB2.0 signal during rewarded trials in the reward omission assay (N = 11). (u,v) Photometry trace and quantification of cue-aligned (0-26 s) GRABeCB2.0 signal during non-rewarded trials in the reward omission assay (N = 11). (w) Schematic for Pavlovian fear conditioning/extinction. (x,y) Photometry trace of GRABeCB2.0 signal during Pavlovian fear conditioning and extinction with representative heatmaps (N = 11). (z,aa) Comparison of GRABeCB2.0 signal during foot shock (0-5 s, p = 0.0049) and cued recall (0-20 s, p = 0.8987) (N = 11). (bb) Summary schematic of elicited eCB release in the NAc. Error bars represent ± SEM (p,q,t,v,z,aa), or min/max values (g,h,l); box plots represent interquartile range with line at median; N represents number of mice. p values reported from one-way ANOVA (g,h,l) and two-tailed paired t test (p,q,t,v,z,aa). *p < 0.05, **p < 0.01.
Extended Data Fig. 2 Histological analysis of CB1R expressing inputs to the NAc and histological/electrophysiological characterization of excitatory aPVTNTS projections.
(a-c) (Top) TdTomato positive NAc projecting neurons in the medial prefrontal cortex (mPFC) Basolateral Amygdala (BLA) and ventral Hippocampus (vHip) and (bottom) RNAscope image and quantification of Cre and CB1R transcript levels in the mPFC BLA and vHip (N = 3,3,3). (d-f) Optically evoked input-output curve from D1R+ and D1R- neurons in vehicle (D1R + n = 19, D1- n = 20, 7.36 · 10−5), SR141617A (D1R+ n = 14 D1R- n = 12, p = 0.0050), and JZL184 (D1+ n = 17, D1R- n = 15, p = 0.9698) treated slices. (g,h) 1μM JZL184 treatment significantly reduces the maximum oEPSC amplitude (p = 0.0193) and increases the Paired Pulse Ratio (p = 0.0355) in D1(+) neurons compared to vehicle treated slices (Veh: n = 19, SR141716A: n = 14, JZL184: n = 16). (i,j) 1μM SR141716A significantly increases the maximum oEPSC amplitude (p = 0.0016) and reduces the Paired Pulse Ratio (p = 0.0092) in D1(−) neurons (Veh: n = 20, SR141716A: n = 11, JZL184: n = 15). (k-n) Electrophysiological comparison of sag current amplitude (p = 0.3871) afterhyperpolzarization voltage (p = 0.3404), holding current (p = 0.6077), and membrane resistance (p = 0.6166) between D1(+) (n = 14) and D1(−) (n = 13) neurons. (o) Schematic for electrophysiological characterization of aPVT output to NAc, BNST, CeA, and ZI. (p-s) Optically evoked input/output curves for aPVT input to the NAc (n = 21), ZI (n = 8) BNST (n = 4), and CeA (n = 4). (t) Comparison of maximal aPVT evoked oEPSC in the NAc, ZI, BNST, and CeA. (u-w) Schematics, representative images, and quantification of collateralization of aPVT projections to the NAc and the ZI, BNST, and CeA, respectively (n = 3,3,3). (x-z) Representative images and quantification of collateralization of pPVT projections to the NAc and the ZI, BNST, and CeA, respectively. Error bars represent ± SEM (d-f,p-s), or min/max values (g-n,t); box plots represent interquartile range with line at median; n represent number of neurons. p values reported from two-tailed unpaired t-test (k-n), one-way ANOVA (g-j), and two-way ANOVA with Holm-Sidak post-hoc correction (d-f). *p < 0.05, **p < 0.01, ****p < 0.0001.
Extended Data Fig. 3 Inhibition of aPVTNTS-NAc terminals is not dependent on reward receipt and is not associated with locomotor activity.
(a) Schematic for fiber photometry recordings from GCaMP6s expressing aPVTNTS-NAc terminals during ad libitum sucrose consumption and Pavlovian reward conditioning. (b) Fiber placements for aPVTNTS-NAc GCaMP cohort displayed in Fig. 2. (c) Photometry trace of GCaMP6s signal aligned to first and last licks in a licking bout during ad libitum consumption (N = 10). (d,e) Transition diagram showing probabilities of transitions between approach, engage, and null/exit behaviors on day 1 and day 5 of conditioning during the post-sipper time window, with size of circles corresponding to number of events. (f,g) Photometry traces of GCaMP6s signal aligned to first lick and last lick on day 1 of reward conditioning. (h,i) Photometry traces of GCaMP6s signal aligned to first lick and last lick on day 5 of reward conditioning. (j,k) Quantification of first lick (0-20 s, p = 0.8736, N = 10) and last lick (0-20 s, p = 0.0004, N = 10) aligned GCaMP6s signals on day 1 and day 5 of reward conditioning. (l) Schematic for fiber photometry recordings from GCaMP6s expressing aPVTNTS-NAc terminals during reward omission. (m,n) Photometry traces of GCaMP6s signal during rewarded and non-rewarded trials during the reward omission test, demonstrating aPVTNTS-NAc terminal inhibition in the absence of reward receipt. (o) Schematic for fiber photometry recordings from GCaMP6s expressing aPVTNTS-NAc terminals during reward conditioning. (p) Photometry trace of GCaMP6s signal on day 5 of reward conditioning during ignored cue/reward trials (no reward engagement), demonstrating a lack of inhibition when the cue/reward is ignored. (q) Schematic for fiber photometry recordings from GCaMP6s expressing aPVTNTS-NAc terminals during Fixed Ratio 1 (FR1) training and Progressive Ratio (PR) training. (r,s) Photometry traces of GCaMP6s signal aligned to cue on day 1 and day 3 of FR1 conditioning, demonstrating inhibition during active operant responding for reward. (t,u) Photometry traces of GCaMP6s signal aligned to cue and nose poke during PR conditioning. (v) Schematic for fiber photometry recordings of GCaMP6s from aPVTNTS-NAc terminals during the open field assay. (w,x) Photometry traces of GCaMP6s signal aligned to walk and rear events during the open field assay. All error bars represent ± SEM; n represent number of neurons. p values reported from two-tailed paired t-test (j,k). ***p < 0.001.
Extended Data Fig. 4 aPVTNTS-NAc terminal activity is negatively correlated with engagement in defensive freezing behaviors.
(a) Schematic for fiber photometry recordings from GCaMP6s expressing aPVTNTS-NAc terminals during Pavlovian fear conditioning and extinction. (b) Photometry trace of GCaMP6s signal during Pavlovian fear conditioning. (c) Quantification of GCaMP6s signal during foot shock (0-5 s, N = 11, p = 0.0002). (d) Photometry trace of GCaMP6s signal aligned to freeze end on extinction day 1. (e) Quantification of GCaMP6s signal aligned to freeze end (0-20 s, N = 11, p = 0.0024). (f) Correlation between freezing time and tone-aligned photometry Z-score (0-20 s) across 5 days of extinction (N = 11, R2 = 0.13, p = 0.0068). (g) Schematic for optogenetic manipulation of aPVT-NAc terminals using DIO-ChR2 and DIO-PPO. (h) 20hz optical stimulation ChR2 (N = 9) expressing aPVTNTS-NAc terminals during the 20 second cue presentation decreases freezing time compared to eYFP (N = 9) controls (p = 0.0053). (i) 10hz optical stimulation of PPO (N = 14) expression aPVTNTS-NAc terminals does not alter freezing behavior compared to eYFP (N = 11) controls (p = 0.1470). Error bars represent ± SEM (c,e), or min/max values (h,i); box plots represent interquartile range with line at median; N represent number of mice. p values reported from paired two-tailed t-test (c,e), unpaired two-tailed t-test (h,i) and simple linear regression (f). **p < 0.01, ***p < 0.001.
Extended Data Fig. 5 Neither activation nor inhibition of aPVTNTS-NAc terminals alters locomotion, drives real-time place preference (RTPP), or drive intracranial self-stimulation (ICSS).
(a) Schematic for optogenetic manipulation of aPVT-NAc terminals using DIO-ChR2 and DIO-PPO. (b) 20hz photo-stimulation of ChR2 expressing aPVTNTS-NAc terminals during cue presentation reduces engagement in consummatory behaviors (ChR2: N = 9 p = 0.0427, eYFP: N = 9), p = 0.8190). (c,d) Neither 20hz photo-stimulation (2 mins off, 2 mins on, 2 mins off) of ChR2 expressing aPVTNTS-NAc terminals (N = 10, p = 0.3655) nor eYFP expressing aPVTNTS-NAc terminals (N = 8, p = 0.8947) alters locomotor behavior. (e) 20hz photo-activation of aPVTNTS-NAc terminals does not support RTPP (ChR2: N = 10 p = 0.8141, eYFP: N = 9 p = 0.7191). (f) 20hz photo-activation (2 second burst per nose poke) of aPVTNTS-NAc terminals does not support ICSS (ChR2: N = 10 p = 0.0968, eYFP: N = 8, p = 0.1670). (g) 10hz photo-inhibition of aPVTNTS-NAc terminals does not support RTPP (PPO: N = 7 p = 0.7492, eYFP: N = 6 p = 0.3216). (h) 10hz photo-inhibition of aPVTNTS-NAc terminals does not support ICSS (PPO: N = 13 p = 0.2565, eYFP: N = 11 p = 0.4527). (i) Fiber placements for optogenetics cohort displayed in Fig. 2. Error bars represent ± SEM (b,e-h), or min/max values (c,d); box plots represent interquartile range with line at median; N represent number of mice. p values reported from paired two-tailed t-test (b,e-h), and one-way ANOVA (c,d). *p < 0.05.
Extended Data Fig. 6 Effect of CB1R antagonism on inhibition of aPVTNTS-NAc terminals and characterization of AAV1-FLEX-sgCNR1 deletion of CB1R from aPVT NTS neurons.
(a) Schematic for fiber photometry recordings from GCaMP6s expressing aPVTNTS-NAc terminals during Pavlovian reward conditioning. (b,c) Photometry trace with representative lick raster and heatmap of cue-aligned GCaMP6s signal in control mice or mice treated with 10 mg/kg SR141716A (N = 10). (d,e) Effect SR141716A on cue-aligned photometry Z-score (p = 0.0064) and total number of licks (p = 0.0270) throughout the session (N = 10). (f,g,h) Diagram showing % of neurons expressing NTS, CB1R, or NTS and CB1R in the aPVT from NTS-Cre mice injected with sgROSA (N = 6) or sgCNR1 (N = 6). (i) Quantification of % of total neurons with CB1R transcript (p = 0.0064). (j) % of NTS neurons CB1R transcript (p = 0.0333). (k) % of NTS neurons expressing CB1R (p = 0.0055). (l)% of total neurons expressing NTS (p = 0.6507). (m) % of total neurons that are NTS positive and CB1R negative (p = 0.0354). (n) % of total neurons that are CB1R positive and NTS negative (p = 0.1512). Error bars represent ± SEM (d,e), or min/max values (i-n); box plots represent interquartile range with line at median; N represent number of mice. p values reported from paired (d,e) and unpaired (i-n) two-tailed t-test. *p < 0.05, **p < 0.01, ***p < 0.001.
Extended Data Fig. 7 Effect of AAV1-FLEX-sgCNR1 deletion of CB1R from aPVT NTS neurons and BLA Vglut1 neurons on reward consumption and Pavlovian reward conditioning.
(a) Schematic for CB1R deletion from the aPVT and photometry recordings ofGCaMP6s aPVTNTS-NAc terminals during adlib sucrose consumption and Pavlovian reward conditioning. (b) Fiber placements for aPVT CB1R deletion/aPVTNTS-NAc cohort in Fig. 3. (c,d) Photometry traces of GCaMP6s signal aligned first lick during adlib sucrose consumption from sgROSA (N = 7) and sgCNR1 (N = 7) injected mice. (e) Comparison of first lick-aligned (0-20 s) GCaMP6s signal between sgROSA (N = 7) and sgCNR1 (N = 7) injected mice (p = 0.1741) during adlib sucrose consumption. (f-h) Comparison of total lick number (p = 0.3705), bout length (p = 0.0393, and bout number (p = 0.2828), between sgROSA (N = 7) and sgCNR1 (N = 7) injected mice during adlib sucrose consumption. (i) Comparison GCaMP6s signal (0-20 s) during rewarded trials in vehicle treated sgROSA (N = 6) and sgCNR1 (N = 7) injected mice during Pavlovian reward conditioning (p = 0.0373). (j) Comparison of # of licks showing a main effect of DO34 treatment on reducing sucrose consumption (Main effects: Treatment p = 0.0025, genotype p = 0.0146, interaction p = 0.5607). (k,l) Photometry traces of GCaMP6s signal aligned to footshock during Pavlovian fear conditioning. (m,n) Photometry traces of GCaMP6s signal aligned to tone presentation of day 1 of fear extinction. (o) Schematic for CB1R deletion from the BLA and photometry recordings of GCaMP6s BLAVGlut-NAc terminals during adlib sucrose consumption and Pavlovian reward conditioning. (p) Fiber placements for BLA CB1R deletion/BLAVGlut-NAc photometry. (q,r) Photometry traces of GCaMP6s signal during adlib sucrose consumption from sgROSA (n = 7) and sgCNR1 (n = 7) mice. (s) # of licks during adlib sucrose consumption (sgROSA N = 7, sgCNR1 N = 7). (t) Comparison of lick-aligned GCaMP6s signal during adlib sucrose consumption. (u-x) Photometry traces of GCaMP6s signal from vehicle and DO34 treated sgROSA (N = 7) and sgCNR1 (N = 7) mice. (y) Comparison of engagement time demonstrating a main effect of both treatment and genotype (Main effects: Treatment p = 0.0025, genotype p = 0.0146, interaction p = 0.5607, (sgROSA N = 7, sgCNR1 N = 7). (z) Comparison of GCaMP6s signal (0-26 s post-cue) demonstrating a main effect of both treatment and genotype (Main effects: Treatment p = 0.0395, genotype p = 0.0066, interaction p = 0.6757, (sgROSA N = 7, sgCNR1 N = 7). (aa) Comparison of total lick number demonstrating a trend toward a main effect of treatment (Main effects: Treatment p = 0.0109, genotype p = 0.4333, interaction p = 0.4333, (sgROSA N = 7, sgCNR1 N = 7). All error bars represent min/max values; box plots represent interquartile range with line at median; N represent number of mice. p values reported from unpaired two-tailed t-test (e-i, s,t) or Two-way ANOVA (j, y-aa). *p < 0.05, **p < 0.01, ***p < 0.001.
Extended Data Fig. 8 AAV1-FLEX-sgCNR1 deletion of CB1R from aPVT NTS neurons does not affect fear conditioning or extinction and foot shock induced eCB release and binding to aPVT terminals in the NAc.
(a) Schematic for CB1R deletion from the BLA and photometry recordings of GCaMP6s BLAVGlut -NAc terminals during Pavlovian fear conditioning and extinction. (b,c) Photometry traces of GCaMP6s signal aligned to foot shock from sgROSA (N = 7) and sgCNR1 (N = 7) injected mice during Pavlovian fear conditioning. (d,e) Comparison of freezing behavior (p = 0.1880) and GCaMP6s signal (0-5 s, p = 0.4953) between sgROSA (N = 7) and sgCNR1 (N = 7) injected animals during fear conditioning. (f,g) Photometry traces of GCaMP6s signal aligned to cue presentation from sgROSA (N = 7) and sgCNR1 (N = 7) injected mice during fear extinction. (h,i) Comparison of freezing behavior (p = 0.9355) and GCaMP6s signal (0-20 s, p = 0.5312) from sgROSA (N = 7) and sgCNR1 (N = 7) injected mice. (j) Schematic for photometry recordings of GRABeCB2.0 or GRABeCBMUT during fear conditioning and fear extinction. (k) Photometry traces of GRABeCB2.0 (n = 8) or GRABeCBMUT (n = 4) signal during fear conditioning. (l) Photometry trace of GRABeCB2.0 (n = 8) or GRABeCBMUT (n = 4) signal during fear extinction. (m) Fiber placements for GRABeCB mice used in j–m and Fig. 4j–n. All error bars represent min/max values; box plots represent interquartile range with line at median; N represent number of mice. p values reported from unpaired two-tailed t-test (d,e,h,i).
Extended Data Fig. 9 Characterization of transsynaptic labeling, specificity of closed-loop optogenetic stimulation, and cluster classification.
(a) Schematic for validation of transsynaptic labeling approach. Transsynaptic AAV1-DIO-FLP was injected into the aPVT, AAV5-DIO-tdTomato was injected into the dorsolateral striatum (DLS), and a cocktail of AAV5-DIO-tdTomato and AAV5-fDIO-eYFP was injected in the NAc of Penk-Cre mice. The DLS, which expresses Penk but does not receive aPVT input, showed expression of AAV5-DIO-tdTomato but not AAV5-DIO-eYFP. The NAc, which expresses Penk and does receives aPVT input, showed expression of both AAV5-DIO-tdTomato and AAV5-DIO-eYFP. (b) Fiber placements for transsynaptic optogenetics cohort in Fig. 4 and representative image. (c-f) Closed-loop optogenetic stimulation during sipper port engagement does not affect % time disengaging, walking, rearing, or grooming. (g) Schematic for closed-loop optogenetic activation of aPVT-NAcPENK neurons during walk events. (h,i) Effect of walk elicited closed-loop activation of aPVT-NAcPENK neurons on % engagement time in ChRimson (N = 7) and eYFP (N = 5) expressing animals. (j,k) Effect of walk elicited closed-loop activation of aPVT-NAcPENK neurons on % walking time in ChRimson (N = 7) and eYFP (N = 5) expressing animals. (l) Lens placements for transsynaptic 1-photon imaging cohort in Fig. 5. (m) Dendrogram plot of hierarchical clustering of 206 tracked aPVT-NAcPenk neurons based on Principal Component Analysis (PCA) of the activity of each neuron in the 60 second window following cue-onset. (n) Percent of variance explained by the top 20 PCs. (o) Support Vector Machine decoding of cluster identity, trained on PavD5 cluster activity.
Extended Data Fig. 10 Validation of generalized linear model and Hilbert and Rayleigh analysis of cluster entrainment to aPVTNTS-NAc GCaMP6s and aPVT-NAc GRABeCB2.0 signal.
(a-e) Cue-aligned traces of % time engaging with sipper port, disengaging from sipper port, walking, rearing, and grooming on PavD1 and PavD5 (N = 10) from GRIN lens implanted animals. (f-g) Heatmap of % of neurons within each cluster encoding each behavioral state on PavD5, with a β coefficient threshold of +/− 0.5. (h,i) Heatmap of % of neurons within each cluster encoding each behavioral state on PavD1, with a β coefficient threshold of +/− 0.5. (j,k) Effect of individually dropping each predictor on model accuracy for PavD5 and PavD1. (l-p) Contribution of effect sizes for each neuron in each cluster to the generalized linear model (T values). (q-u) Overlay of Pavlovian reward conditioning day 5 aPVT-NAcPenk cluster GCaMP6s traces with aPVTNTS-NAc terminal GCaMP6s traces and Rayleigh plots from cluster 1, 2, 3, 4, and 5. (v-z) Overlay of Pavlovian reward conditioning day 5 aPVT-NAcPenk cluster GCaMP6s traces with aPVT-NAc terminal GRABeCB2.0 traces and Rayleigh plots from cluster 1, 2, 3, 4, and 5.
Supplementary information
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
Reprints and permissions
About this article
Cite this article
Marcus, D.J., English, A.E., Chun, G. et al. Endocannabinoids facilitate reward engagement through retrograde gain control. Nature (2026). https://doi.org/10.1038/s41586-026-10967-w
Download citation
Received:
Accepted:
Published:
Version of record:
DOI: https://doi.org/10.1038/s41586-026-10967-w