The paraventricular thalamus provides a polysynaptic brake on limbic CRF neurons to sex-dependently blunt binge alcohol drinking and avoidance behavior in mice

Bed nucleus of the stria terminalis (BNST) neurons that synthesize corticotropin-releasing factor (CRF) drive binge alcohol drinking and anxiety. Here, we found that female C57BL/6J mice binge drink more than males and have greater basal BNSTCRF neuron excitability and synaptic excitation. We identified a dense VGLUT2 + synaptic input from the paraventricular thalamus (PVT) that releases glutamate directly onto BNSTCRF neurons but also engages a large BNST interneuron population to ultimately inhibit BNSTCRF neurons, and this polysynaptic PVTVGLUT2-BNSTCRF circuit is more robust in females than males. Chemogenetic inhibition of the PVTBNST projection promoted binge alcohol drinking only in female mice, while activation reduced avoidance behavior in both sexes. Lastly, repeated binge drinking produced a female-like phenotype in the male PVT-BNSTCRF excitatory synapse without altering the function of PVTBNST neurons per se. Our data describe a complex, feedforward inhibitory PVTVGLUT2-BNSTCRF circuit that is sex-dependent in its function, behavioral roles, and alcohol-induced plasticity.

A lcohol use disorder (AUD) is highly co-expressed with other neuropsychiatric diseases including anxiety disorders, with women having increased susceptibility to this comorbidity compared to men (61% vs. 35%, respectively) 1 . Binge alcohol drinking is a primary risk factor for the development of these conditions, and females across mammalian species display greater binge drinking and transition from first alcohol use to disease states more quickly than males [2][3][4][5][6] . These initial differences in alcohol sensitivity suggest that the mechanisms underlying early drinking contribute to disease susceptibility and are important targets for intervention 2 . However, there is a relative dearth of studies investigating the mechanisms controlling binge alcohol drinking in females and the sex differences in the expression of this behavior and consequent vulnerability to disease. The bed nucleus of the stria terminalis (BNST) is a hub in the brain circuits underlying anxiety and alcohol/substance use disorders in humans and is highly sexually dimorphic in mammals 7,8 . The BNST is enriched with neurons that synthesize and release corticotropin-releasing factor (CRF), a stress neuropeptide involved in the development and maintenance of anxiety and addictive disorders, and activation of BNST CRF neurons drives binge drinking behavior and produces anxiety 9,10 ; however, the identity and organization of upstream excitatory circuits controlling BNST CRF neuron function and its role in alcohol drinking and anxiety are poorly understood. The BNST is anatomically connected to the paraventricular nucleus of the thalamus (PVT), a region likewise implicated in the etiology of alcohol and substance use disorders and associated behaviors including anxiety [11][12][13][14] . Human neuroimaging studies have demonstrated a functional connection between the thalamus and BNST that is denser in females than males 7,8 and decreased thalamic projection strength in young individuals with alcohol abuse 15 . These converging lines of evidence suggest that the PVT-BNST projection might contribute to sex-dependent behavior via modulation of BNST CRF neurons. Here, we systematically examined in both sexes the anatomical and functional architecture of the PVT-BNST CRF circuit and its role in binge alcohol drinking and anxiety behaviors, as well as alcohol-induced plasticity.

Results
Greater BNST CRF neuron excitation and alcohol drinking in females. First, we examined the relationship between sexdependent binge alcohol consumption and BNST CRF neuron excitation. We showed that female C57BL/6J mice consistently consume more alcohol than males using the Drinking in the Dark (DID) 16,17 model of binge drinking (Fig. 1a). Using whole-cell patch-clamp slice electrophysiological recordings in BNST CRF neurons (Fig. 1b), we found that the proportion of BNST CRF neurons active in their basal state was more than twice as high in females compared to males (60% in females vs. 24% in males; Fig. 1c). Further evaluation of BNST CRF neurons in slice showed that the frequency of spontaneous excitatory postsynaptic currents (sEPSCs), including activity-independent "miniature" EPSCs (mEPSCs), was higher in females than males, but the frequency of inhibitory postsynaptic currents (sIPSCs and mIPSCs) was not (Fig. 1d, e; Supplementary Fig. 1). These results suggest that spontaneous glutamate, but not GABA, release onto BNST CRF neurons is greater in females than males, contributing to an overall increased synaptic drive onto BNST CRF neurons in females (biased toward excitation) compared to males (Fig. 1f).
To determine whether increased excitability and excitation of BNST CRF neurons may be related to potentiated alcohol intake in females, we assessed the ability of chemogenetic inhibition of BNST CRF neurons to suppress binge alcohol consumption in males and females. By comprehensively analyzing the effects of sex and agonist dose in a broader dataset from a mixed-sex cohort of mice we previously published 18 , we found that females with the hM4D Gi-coupled designer receptor exclusively activated by designer drug (Gi-DREADD) in BNST CRF neurons required a higher dose of CNO (10 mg/kg) to attenuate binge drinking than their male counterparts (3 mg/kg; Supplementary Fig. 2). Together, these data suggest that the level of activity in the BNST CRF neuron population was sufficiently higher in females to necessitate more robust functional inhibition to affect binge drinking behavior.
Robust projection from the paraventricular thalamus modulating BNST CRF neurons. The sex differences in the underlying physiology and excitatory synaptic input to BNST CRF neurons led us to examine whether there are sex differences in the anatomical and/or functional density of specific glutamatergic inputs to the BNST that may modulate BNST CRF neurons. Using a viral retrograde tracing approach in VGLUT2-ires-Cre (VGLUT2-Cre) mice to label VGLUT2-positive (VGLUT2+) and negative (VGLUT2-) BNST-projecting neurons (Fig. 2a), we identified several known sources of excitatory input to the BNST, including a well-characterized input from the basolateral amygdala (BLA) shown to reduce avoidance behavior in males 19,20 . We found that the brain region with the densest projection to the BNST was the PVT, which was similarly robust in both sexes across the anterior-posterior extent of the PVT and enriched in the anteriormid PVT (Fig. 2b, [21][22][23][24] . To determine whether PVT VGLUT2 neurons play a role in binge drinking behavior, we injected a Cre-dependent virus expressing the kappa opioid receptor (KOR)-based Gi-DREADD (Gi-KORD) 25 or control (CON) virus into the PVT of VGLUT2-Cre mice ( Supplementary  Fig. 4a). Following administration of the Gi-KORD ligand Salvinorin B (SalB; 17 mg/kg, s.c.), male and female Gi-KORD mice but not CONs displayed blunted binge drinking behavior (Supplementary Fig. 4b). In contrast, SalB activation of the Gi-KORD did not alter sucrose drinking or avoidance behavior in the open field test (OF; Supplementary Fig. 4c-e), suggesting that PVT VGLUT2 neurons drive alcohol binge drinking behavior specifically without altering reward-seeking behavior or anxiety more generally.
We next evaluated sex differences in the function of BNSTprojecting PVT (PVT BNST ) neurons. Whole-cell patch-clamp recordings of PVT BNST neurons indicated no sex differences in the intrinsic or synaptic excitability of this neuron population apart from decreased sIPSC amplitude in females ( Supplementary  Fig. 5). To establish whether there are sex differences in the function of the PVT BNST circuit and postsynaptic responses from BNST CRF neurons, we injected a CaMKIIα-driven channelrhodopsin (ChR2) virus into the PVT of CRF-reporter mice and recorded from BNST CRF neurons (Fig. 3a, b). We found that ChR2 activation of the PVT synaptic input (2 ms pulse of 490 nm LED) elicited an optically-evoked monosynaptic EPSC and polysynaptic IPSC (oEPSC and oIPSC, respectively) in BNST CRF neurons, as both currents could be abolished by bath application of the voltagegated sodium channel blocker tetrodotoxin (TTX, 500 nM) but only the oEPSC was restored by addition of the voltage-gated potassium channel blocker 4-aminopyridine (4-AP, 100 μM; Fig. 3c,d). Investigation of potential sex differences in PVT-BNST CRF excitatory synapses revealed that the paired-pulse ratios of pharmacologically-isolated PVT-evoked oEPSCs in BNST CRF neurons separated by 50 ms were higher in males than females, suggesting a greater probability of evoked presynaptic glutamate release from PVT terminals in the BNST in females (Fig. 3e). In contrast, AMPA/NMDA ratios of oEPSCs were similar between males and females, indicating no postsynaptic sex difference (Fig. 3f). Intriguingly, PPR was below 1.0 in females (Fig. 3e), suggesting that these PVT terminals require relatively little excitation to elicit glutamate release, and thus PVT-BNST CRF synapses may serve as low-pass synaptic filters in females. Further, quantification of oEPSC and oIPSC amplitudes within individual neurons across increasing LED power showed a more robust response in BNST CRF neurons of females compared to males (Fig. 4a, b, Supplementary Fig. 6a, b), providing converging evidence for increased evoked glutamate release from PVT terminals in the BNST. In conjunction with increased spontaneous glutamate release onto BNST CRF neurons in females (Fig. 1, Supplementary Fig. 1), these results suggest higher PVT glutamate tone in the BNST of females than males.
However, oIPSCs were larger than oEPSCs in most individual BNST CRF neurons in both sexes, resulting in a net effect of synaptic inhibition by PVT afferent activation in both males and females (Fig. 4c), similar to what has been described for other limbic projections of the PVT [22][23][24] . In addition, the latency between optical excitation of PVT terminals and initiation of oPSCs was significantly longer for oIPSCs than oEPSCs within individual BNST CRF neurons in both sexes (Fig. 4d), providing converging evidence that the PVT-evoked inhibition onto BNST CRF neurons was polysynaptic in nature. These results suggest that increased glutamate release from PVT inputs in females also results in increased polysynaptic inhibition of BNST CRF neurons. We confirmed that activating PVT BNST glutamatergic inputs results in a net hyperpolarization of the membrane potential of BNST CRF neurons using whole-cell current-clamp recordings. LED excitation of ChR2 from PVT terminals produced an optically-evoked postsynaptic potential (oPSP) that was negative (hyperpolarizing) in most BNST CRF neurons in both males and females, regardless of their resting membrane potential (Fig. 4e,f). Evaluation of oPSPs in response to PVT ChR2 excitation across a range of frequencies ( Supplementary Fig. 6c,d) showed that 20 Hz stimulation (the same frequency used to evaluate PPR) elicited a greater magnitude second oPSP than the first oPSP in BNST CRF neurons of males but no difference at any frequency in females, consistent with higher PPR in males than females. Thus, while female PVT inputs readily release glutamate, those in males have a lower probability of release but may be more prone to synaptic facilitation upon repeated stimulation, indicating an increased capacity for plasticity.  PVT BNST circuit regulation of alcohol drinking and anxiety behaviors. Given the inhibitory effect of ex vivo PVT BNST afferent activation on BNST CRF neuron excitability, we next evaluated the role of the PVT BNST circuit in binge drinking and anxiety-like behavior. We used a multiplexed chemogenetic strategy 25 to bidirectionally and independently manipulate the PVT BNST circuit during behavior by injecting either a cocktail of a Cre-dependent excitatory (hM3D) Gq-DREADD + Credependent inhibitory Gi-KORD or a Cre-dependent CON vector in retrogradely Cre-labeled PVT BNST neurons (Fig. 5a, b; Supplementary Fig. 7). We found that CNO administration (5 mg/kg, i.p.) to activate the Gq-DREADD prior to alcohol access during Day 4 DID did not reliably alter binge consumption of 20% alcohol compared to vehicle baseline (Fig. 5c); however, there was a trend for a decrease in consumption in DREADD males (see Supplementary Fig. 8a), suggesting that activation of the PVT BNST pathway may be sufficient to blunt binge alcohol drinking in a subset of males. In contrast, Salvinorin B (SalB) administration (17 mg/kg, s.c.) to activate the Gi-KORD robustly increased binge alcohol drinking in DREADD females but had no effect in DREADD males or either CON group (Fig. 5d, Supplementary Fig. 8b), suggesting that tonic activity of the PVT BNST circuit ( Supplementary Fig. 5c) engaging interneurons is necessary for active suppression of alcohol drinking behavior in females. As BNST CRF neurons show larger postsynaptic responses to the PVT BNST input in females, reducing the activity of the PVT BNST input may result in more robust disinhibition of BNST CRF neurons to further enhance binge drinking. Intriguingly, neither chemogenetic manipulation affected 1% sucrose consumption in a similar DID paradigm ( Supplementary Fig. 8c, d). These results suggest that the PVT BNST projection plays a unique inhibitory role in binge drinking behavior that is not dependent on modulation of the general rewarding or aversive aspects of drug intake and consummatory behavior previously shown to be mediated by other major limbic outputs of the PVT such as the nucleus accumbens 23,26,27 . We also probed whether the PVT BNST circuit plays a role in anxiety-like behavior, given our previous work showing that activity of the BNST CRF neuron population is anxiogenic 18,28 and that the PVT regulates emotional behaviors 14,22 . CNO administration to activate the Gq-DREADD decreased avoidance of the open arms on the elevated plus maze (EPM) in DREADD mice Fig. 2 The PVT provides a dense glutamatergic projection to the BNST. a Viral unilateral retrograde tracing strategy in the BNST of VGLUT2-Cre mice (N's = 4 M, 4 F). b % cells in the basolateral amygdala (BLA) and paraventricular nucleus of the thalamus (PVT) sections that project to the BNST (the densest sources of glutamatergic input identified), calculated as (total # of retrogradely-labeled cells/total DAPI+) × 100 and presented as mean values ± SEM. Mixed-effects model: main effect of subregion (F 4,23 = 32.01, P < 0.0001, not indicated) but no effect of sex (P = 0.967) or interaction (P = 0.524); post-hoc two-tailed t-tests with Holm-Sidak corrections between subregions: BLA has fewer BNST-projecting neurons than the PVT across all A/P coordinates (all $$$$ Ps < 0 . 0001 ) , and that the amPVT has more than the mPVT (t 23 = 2.94, *P = 0.037) and pPVT (t 23 = 3.15, *P = 0.027), with no other differences: aPVT vs. amPVT (P = 0.085); aPVT vs. mPVT (P = 0.873); aPVT vs. pPVT (P = 0.873); mPVT vs. pPVT (P = 0.873). aPVT, anterior PVT; amPVT anterior-mid PVT; mPVT, mid PVT: pPVT, posterior PVT. c Proportion of VGLUT2-Cre+ (EGFP-labeled) and VGLUT2-Cre− (tdTomatolabeled) BNST-projecting PVT (PVT BNST ) neurons in males and females, showing that almost all PVT BNST neurons are VGLUT2+ in both males and females with no difference between sexes (two-tailed unpaired t-test: t 6 = 1.24, P = 0.261). d Representative images of coronal brain slices from a virusinjected mouse illustrating the expression of DAPI (blue) and all BNST-projecting cells (both VGLUT2-Cre + and -in green) in the BLA and across the A/P extent of the PVT. e Density heat maps illustrating the average number of BNST-projecting cell bodies within a 50 µM radius of each identified projection neuron for samples from each sex, matched and scaled similarly to representative images in d (top: males, bottom: females). compared to CON mice ( Fig. 5e-g), an anxiolytic effect of Gq-DREADD activation that we replicated in the open field test (OF; Supplementary Fig. 8e,f). In contrast, SalB administration to activate the Gi-KORD did not affect avoidance of the center of the open field ( Fig. 5h-j). Together, these results suggest that in an anxiogenic context when BNST CRF neuron activity is high 29,30 , activation of the PVT BNST synaptic brake is sufficient to decrease anxiety-like behavior.
Repeated alcohol use produces a female-like phenotype in BNST CRF synapses in males. Finally, we examined the sexdependent plasticity in this circuit following repeated binge alcohol drinking to understand whether there are changes in the PVT VGLUT2 -BNST CRF circuit during withdrawal that could precipitate or contribute to increased disease vulnerability. We found that one day after three cycles of EtOH DID (Fig. 6a), the proportion of BNST CRF neurons in an active state was increased following alcohol exposure in males (Fig. 6b); intriguingly, this proportion was unchanged by a history of alcohol exposure in females, suggesting that basal BNST CRF neuron population-level activity may be near a maximum in naïve females ( Fig. 1b;   Fig. 6b) in the absence of a discrete salient stimulus (for example, during an anxiety assay such as EPM here). Alternatively, the population level excitability of BNST CRF neurons in females may be relatively impervious to the effects of extended alcohol exposure. Notably, a range of neural activity patterns was observed in all groups, illustrating the heterogeneity of the BNST CRF population in both sexes as previously described in males 31,32 . Most other measures of intrinsic excitability and current-injected firing were neither different between males and females nor altered following EtOH DID exposure when measured either at their resting membrane potential (RMP; Supplementary Fig. 9) or at a common hyperpolarized potential of −70 mV ( Supplementary  Fig. 10). This suggests that the primary effect of repeated alcohol exposure on BNST CRF neuron excitability is on the populationlevel activity, with males developing a female-like phenotype. Examination of synaptic transmission showed that sEPSC frequency was increased in EtOH males but unchanged in females (Fig. 6c), while sEPSC amplitude was unchanged in both sexes (Fig. 6d). This led to an increase in excitatory synaptic drive onto BNST CRF neurons in males but not females (Fig. 6e). In contrast, there were no changes in the frequency, amplitude, or synaptic drive of sIPSCs in either sex (Fig. 6f-h). Assessment of the kinetics of sPSCs revealed increased half-width and weighted tau of sEPSCs in EtOH males, consistent with an increase in glutamatergic transmission ( Supplementary Fig. 11).
We further examined plasticity in PVT-evoked postsynaptic responses from BNST CRF neurons and found no significant effect of EtOH on oEPSC raw amplitude across increasing LED power (Fig. 7a). However, when oEPSC amplitude was normalized to maximum within each cell, there was a leftward shift in the power-response curve in males but not females ( Supplementary  Fig. 12a, b), demonstrating that lower optical power was able to elicit larger oEPSCs in EtOH males compared to controls. Further, analysis of the slope of the initial increase in raw oEPSC amplitude across 5-25% LED power confirmed this alcoholinduced increase in the excitatory postsynaptic responses of BNST CRF neurons in males ( Supplementary Fig. 12c, d), mirroring the alcohol exposure effect on sEPSCs in males (Fig. 6c) and suggesting that PVT-BNST CRF synapses comprise a large contingent of the glutamate synapses vulnerable to alcoholinduced plasticity. oIPSCs were also somewhat enhanced by alcohol exposure, particularly in males, as analysis of raw oIPSCs across LED power showed a sex × EtOH interaction (Fig. 7b), however, normalized oIPSC amplitude was not significantly altered in BNST CRF neurons of EtOH mice in either sex ( Supplementary Fig. 12e, f); the initial slope was increased in EtOH mice, but this was not significantly different specifically within either sex ( Supplementary Fig. 12g, h). Following voluntary alcohol exposure, the E/I ratio of oPSCs remained significantly below 1.0 in males but not females (Fig. 7c), suggesting that while the overall E/I balance toward inhibition was maintained in males (that is, increased PVT-mediated excitation was met with increased polysynaptic inhibition), the net inhibitory effect of PVT afferent activation at BNST CRF neurons was attenuated in females. However, oEPSC/oIPSC ratios were not different between CON and EtOH females, suggesting that this statistical lack of E/I ratio below one may be due to increased variability in the responses of individual neurons/ synapses in females rather than a robust loss of polysynaptic inhibition, Notably, we observed very few sex differences or alcohol-induced changes in the electrophysiological properties, excitability, and synaptic transmission of PVT BNST neurons themselves ( Supplementary Figs. 13 and 14), suggesting that sexdependent and alcohol-dependent differences in the PVT-BNST CRF synapse strength are not due to differences in PVT BNST neuronal excitability per se. Together, our data suggest that the effects of sex and repeated binge alcohol drinking are on the tonic activity of the BNST CRF neuron population, due to basal differences in and differential plasticity in PVT VGLUT2 -BNST CRF synapses and BNST CRF neuron responses to the largely tonicallyactive PVT BNST afferents.

Discussion
Altogether, we found that female mice display more robust binge drinking behavior than males and that BNST CRF neurons are both more tonically active and receive increased excitatory synaptic input in females relative to males. We further demonstrated that the PVT provides a dense excitatory input to the BNST in both sexes 33 that is functionally more robust in females than males, consistent with human neuroimaging studies on thalamic outputs 7,8 . We report that while PVT glutamate neurons directly synapse onto BNST CRF neurons as recently shown anatomically in males 33 , they also provide robust feedforward inhibition of this population by recruiting interneurons, providing a synaptic "brake" on BNST CRF neuron excitability; such feedforward inhibition has also been reported for other PVT projections, such as those to the nucleus accumbens (NAc) and central  amygdala (CeA) 23,24 , to modulate appetitive, drug, and fearrelated behaviors in males [22][23][24]27,34 . For example, the PVT regulates appetitive learning and behavioral responses for the availability of sucrose, food, and water, particularly through its projection to the NAc [35][36][37] . Intriguingly, we found that inhibition of the entire PVT glutamate neuron population suppressed voluntary, predictable binge drinking behavior (without an effect on palatable binge sucrose consumption using a parallel paradigm; Supplementary Fig. 4), consistent with many of these previously reported complex roles of the PVT in motivated behaviors. In contrast, we found that the PVT BNST circuit plays a sex-dependent but opposite role in binge alcohol drinking, such that removal of the PVT BNST brake disinhibits alcohol consumption in females without affecting sucrose intake. In contrast, PVT BNST activation may inhibit alcohol consumption in a subset of males, pointing to a difference in the set point for engagement of the PVT-BNST CRF circuit (and thus modulation) between sexes. Our findings provide insight into the mechanisms underlying binge alcohol drinking in females. They also highlight underlying circuit features that are different between males and females, conferring complex sex-dependent control and circuit modulation of binge alcohol consumption. They further point to a unique role for the PVT BNST circuit in alcohol drinking that is independent of the role of the PVT in general appetitive reward.
In addition, we found that activation of the PVT BNST pathway in aversive contexts, when BNST CRF neuron activity is high 29,30 , is sufficient to reduce avoidance behavior in both sexes, similar to the role of the BLA input in males 19 . These results suggest a critical sex-dependent role of this thalamo-limbic circuit in the control of and relationship between alcohol drinking and anxiety states. Further, repeated binge drinking produced a female-like phenotype in PVT-BNST CRF synapses and BNST CRF neurons in males, characterized by increased excitation observed in naïve females and associated with greater binge drinking behavior. Notably, while the overall excitability phenotype becomes more similar to a female state, the mechanism(s) driving the alcoholinduced increase in excitability in males may be distinct from the underlying mechanisms of sex differences in excitability in the naïve population. However, our results provide a solid basis for future examination of the molecular similarities and differences between BNST CRF neurons in males and females, as well as the effects of chronic alcohol. In addition, other glutamatergic inputs  > 0.05). c oEPSC/ oIPSC ratio in BNST CRF neurons of EtOH mice. Data are presented as mean values ± SEM. EtOH M have E/I ratios below 1 (one-sample t-test: t 13 = 3.18, ## P = 0.007), which does not differ from that of CON M (two-tailed unpaired t-test: t 23 = 1.01, P = 0.321). EtOH F does not have an E/I ratio below 1.0 (one-sample t-test: t 6 = 1.71, P = 0.137), however, this ratio does not significantly differ from CON F (two-tailed unpaired t-test: t 20 = 1.84, P = 0.081). and the BNST CRF neuron population itself interact highly with other BNST neuron subpopulations, including PKCδ 33 , neuropeptide Y 18 , and dynorphin 20 , among others. However, the literature is not yet comprehensive and provides some directly conflicting evidence regarding the independence and interactions between some BNST subpopulations. For example, one recent report suggests that the CRF and PKCδ populations in the male BNST receive direct synaptic input from the PVT and other regions as independent subpopulations that are inhibitory upon one another and modulate anxiety-like behavior in opposing manners (with CRF neurons being anxiogenic) 33 . In contrast, another recent study shows that these are highly overlapping populations and that acute stress recruits additional PKCδ expression in CRF neurons in females but not males 38 , suggesting an important and perhaps sex-dependent function of the overlapping population of neurons. While both studies find that stress activates both types of neurons, they come to different conclusions about the organization and role(s) of these neurons in behavior. Further, BNST CRF neurons are themselves complex, serving as both projection and interneurons 28 ; as such, some may laterally inhibit others to participate in the polysynaptic inhibition we observed here. Altogether, our results and others point to the need for future studies to disentangle the independence, overlap, and interactions between BNST subpopulations, with special attention paid to the effects of sex and various physiological stressors (including alcohol and other drugs of abuse). Understanding the microcircuit organization and molecular makeup of these neurons is critical to defining the broader circuit architecture containing the PVT VGLUT2 -BNST CRF circuit described here.

Methods
Subjects. All experimental mice were male and female adult mice on a C57BL/6J background strain. Wild-type C57BL/6J mice were purchased as adults from Jackson Laboratory, and all transgenic lines were bred in our animal facility. CRFires-Cre (CRF-Cre) 18,39 and VGLUT2-ires-Cre (VGLUT2-Cre) 40 mice were bred with WT C57BL/6J mice, and hemizygous CRF-Cre mice were bred with homozygous floxed Ai9-tdTomato or floxed EGFP-L10a mice purchased from Jackson Laboratory (stocks 007909 and 024750) to produce CRF-Cre-reporter mice. Mice were group-housed with ad libitum access to food and water in colony room on a 12:12 h reverse light cycle, with lights off at 7:30 a.m. Mice were singly housed for one week prior to the onset of behavioral experiments and remained singly housed thereafter. Experiments began approximately 3 h into the dark phase of the light cycle. All experimental procedures were approved by the Institutional Animal Care and Use Committees at Weill Cornell Medicine and the University of North Carolina-Chapel Hill.
Behavior assays. The standard Drinking in the Dark (DID) binge alcohol drinking paradigm was used in mice to model human binge consumption behavior 17 . For each cycle of EtOH DID, three hours into the dark cycle, the home cage water bottle was replaced with a bottle containing 20% (v/v) alcohol (EtOH) for two hours on Days 1-3 and four hours on Day 4, followed by three days of forced abstinence between consecutive cycles. A similar access schedule was used to evaluate binge sucrose consumption, except that home cage water bottles were replaced with 1% (w/v) sucrose. For all drinking experiments, empty "dummy" cages on the same rack as behavior mice received the same EtOH or sucrose bottle replacement, and consumption was adjusted for a leak from dummy bottles and then normalized to bodyweight. The open-field test (OF) was used to evaluate avoidance and locomotor behavior as previously performed 18 . Mice were placed in the 50 × 50 cm arena for 60 min, and Ethovision video tracking (Noldus, Wageningen, Netherlands) was used to quantify raw locomotor and location data used to calculate measures including distance traveled and time spent in each compartment (center vs. periphery, total). The elevated plus maze (EPM) was also used to assess anxiety-like behaviors and was conducted in a plexiglass maze with two open and two closed arms (35 cm length × 5.5 cm width, with walls 15 cm tall over the closed arms). Mice were placed in the center of the EPM for five-minute trials and movement and time spent in each compartment were tracked using SMART 2.0 and Ethovision 10 and 11. Total time and percent time spent in each arm were quantified.
Stereotaxic surgeries. For experiments requiring site-directed administration of viral vectors or retrobeads, mice were anesthetized with 2% isoflurane (VetEquip, Livermore, CA) in 0.8% oxygen in an induction chamber (VetEquip, Livermore, CA) then placed in an Angle Two mouse stereotaxic frame (Leica Biosystems, Wetzlar, Germany) and secured with ear bars into a nose cone delivering isoflurane to maintain anesthesia. Mice were given a subcutaneous injection of meloxicam (2 mg/kg) for preemptive analgesia and 0.1 mL of 0.25% Marcaine around the incision site. A Neuros 7000 series 1 µL Hamilton syringe with a 33-gauge needle (Reno, NV) connected to a remote automated microinfusion pump (KD Scientific, Holliston, MA) was used for construct delivery at a rate of 50-100 nL/min to the PVT  Fig. 5a, b). Approximately one week following surgery, mice started the DID procedure with a baseline cycle followed by a cycle in which they received 0.9% sterile saline vehicle (10 ml/kg, i.p.) injections on Days 2 and 4 forty min prior to alcohol access to habituate to the injection procedure. Chemogenetic manipulations began the next cycle with saline on Day 2 and clozapine-n-oxide (CNO, 5 mg/kg in 0.9% saline) on Day 4. The Gi-KORD was subsequently evaluated similarly with a DMSO vehicle (1 ml/kg, s.c.) injection cycle and manipulation cycle with Salvinorin B (SalB, 17 mg/kg in DMSO). Mice displaying injection stress defined by >1 g/kg reduction in and <1 g/kg vehicle baseline drinking were excluded from the analysis. For the OF and EPM, half of the animals of each sex received the vehicle and half received the DREADD activator (CNO or SalB) 40 min prior to the assay. For sucrose DID, the same SalB or CNO drug administration procedure was used as that for EtOH DID.
Brain extraction and fluorescence immunohistochemistry. Following behavior procedures, mice were deeply anesthetized with pentobarbital (100 mg/kg, i.p.) and transcardially perfused with sterile phosphate-buffered saline (PBS) followed by 4% paraformaldehyde (PFA). Brains were extracted, post-fixed overnight in 4% PFA, and then placed in PBS until they were sliced on the coronal plane in 45 μm sections on a VT1000S vibratome (Leica Biosystems) to check injection placements and viral expression (hit maps of these expression data are presented in Supplementary Fig. 7a, b). To amplify the expression of fluorophore tags, coronal slices containing the PVT and BNST from DREADD and CON mouse brains underwent immunofluorescence staining. Slices were washed twice in PBS followed by 0.2% triton (Fisher Bioreagents, Hampton, NH) for 10 min each and then blocked in 5% normal donkey serum (NDS) for 30 min. Tissue was incubated in primary antibody (DsRed rabbit polyclonal 1:400, Takara, Kusatsu, Japan) overnight at room temperature. The next day, slices were washed 3 times in 0.2% triton for 10 min each and then blocked in 5% NDS for 30 min. Tissue was then incubated in secondary antibody (Alexafluor-568 donkey anti-rabbit, 1:250, Invitrogen, Carlsbad, CA) for 2 h and subsequently washed twice in 0.2% Triton and then in PBS for 10 min each. Slices were counterstained with DAPI and mounted on slices, and coverslipped with Vectashield hard mount antifade mounting medium (Vector Labs, Burlingame, CA), and stored in the dark at 4°C until imaged (as described below) to verify surgical placements and viral expression.
Retrograde neuronal tracing. For anatomical circuit tracing of excitatory inputs to the BNST, VGLUT2-Cre mice received stereotaxic injections of a retrograde virus (AAVrg-Ef1a-DO-DIO-TdTomato_EGFP-WPRE-pA) unilaterally into the BNST, which was is retrogradely trafficked to ultimately express GFP in Cre+ cells and tdTomato in Cre− cells that project to the BNST. Following three weeks to allow for optimal viral expression, mice were sacrificed and their brains harvested for quantification.
Image acquisition and analysis. Coronal brain slices were collected and imaged for all experiments to confirm and quantify viral expression and immunolabeling. Images were acquired on a Zeiss LSM 880 Laser Scanning Confocal microscope (Carl Zeiss, Oberkochen, Germany). For the retrograde tracing experiment, images of the PVT and BLA were quantified using ImageJ 1.50 (US National Institute of Health) to count the total number of DAPI-stained nuclei, GFP+, and tdTomato+ cells. The coordinates of each cell were analyzed using a custom MATLAB (MathWorks, Natick, MA) program and normalized to the most dorsomedial point of the PVT. Heatmaps of the density of BNST-projecting neurons were generated as described elsewhere 41 , overlaying data from four mice per sex. The proportion of PVT projectors that are VGLUT2+ was calculated as GFP+/(GFP + plus tdTo-mato+). Coordinates analyzed (mm from Bregma): BLA (−1.55), anterior (aPVT, −0.46 mm), anterior-mid (am-PVT, −0.70), mid (mPVT, −1.22), and posterior (pPVT, −1.91). ex vivo slice electrophysiology and calcium imaging. Slice electrophysiology experiments were performed as previously conducted in our laboratory 9,18 . Mice were decapitated under isoflurane anesthesia and their brains were rapidly extracted. Coronal BNST slices (300 µm) were prepared on a VT1200 vibratome (Leica Biosystems) in ice-cold, oxygenated (95% O 2 /5% CO 2 ) sucrose artificial cerebrospinal fluid (aCSF) containing (in mM): 194 sucrose, 20 NaCl, 4.4 KCl, 2 CaCl2, 1 Mg Cl2, 1.2 NaH2PO4, 10 glucose, and 26 NaHCO3 (pH 7.3 and 30 mOsm). Slices were transferred to a holding chamber with 30°C oxygenated normal aCSF (in mM):124 NaCl, 4.4 KCl, 2 CaCl2, 1.2 MgSO4, 1 NaH2PO4, 10 glucose, and 26 NaHCO3 (pH 7.3 and 30 mOsm) and allowed to equilibrate for at least one hour. For electrophysiological recordings, slices were transferred to a submerged recording chamber (Warner Instruments, Hamden, CT) and perfused at a rate of 2 mL/min with 30°C oxygenated normal aCSF.
BNST CRF neurons were identified for recording with their tdTomato tag using a 580 nm LED under 40x objective (Olympus, Tokyo, Japan). Signals were acquired using a Multiclamp 700B amplifier (Molecular Devices), digitized, and analyzed via pClamp 10 software (Molecular Devices). Input resistance and access resistance were continuously monitored throughout experiments, and cells in which access resistance changed by more than 20% were not included in data analysis. Excitability experiments were performed in current-clamp configuration using a potassium gluconate-based intracellular recording solution containing (in mM): 135 KGluc, 5 NaCl, 2 MgCl 2 −6H2O, 10 HEPES, 0.6 EGTA, 4 Na-ATP and 0. Prior to slice electrophysiology experiments, mice underwent three cycles of EtOH DID (EtOH group) or a water control DID procedure in which the replacement bottle contained water instead of EtOH (CON group) to allow for investigation of basic sex differences in function, as well as for direct comparison to EtOH mice. Twenty-four hours following the onset of the last EtOH or water bottle access, mice were sacrificed for slice electrophysiology experiments as described above. For experiments investigating PVT BNST neurons, mice received bilateral intra-BNST injections of green or red retrobeads (250 nL, Lumafluor) prior to DID to label this population for identification during recordings. For experiments characterizing postsynaptic responses in BNST CRF neurons to PVT glutamate inputs, mice received an intra-PVT stereotaxic injection of AAV5-CamKIIα-hChR2(H134R)-eYFP.WPRE.hGH (200 nL, Penn Vector Core and Addgene) three weeks prior to DID procedures. During recordings, 1 ms 490 nm LED stimulation was used to optically stimulate ChR2 + PVT cell bodies to elicit action potentials to confirm sufficient PVT expression and ChR2 fidelity at 1, 2, 5, 10, 20, and 50 Hz. One 2 ms stimulation every 10 s was used in the BNST to optically-evoke glutamate release from PVT terminals while recording postsynaptic responses (oPSCs) from BNST CRF neurons in voltage-clamp, and 2 ms stimulation at 1, 2, 5, 10, and 20 Hz was used while measuring postsynaptic potentials (oPSPs) in BNST CRF neurons in current-clamp. Polysynaptic and monosynaptic oPSCs were assessed with the bath application of TTX (500 nM) to block all neurotransmission followed by the addition of 4-aminopyridine (100 µM) to reinstate monosynaptic transmission.
Following behavior, a subset of PVT BNST DREADD mice were sacrificed for slice electrophysiology confirmation of the approach using bath application of CNO (10 µM) during current-clamp recordings of DREADD + PVT cell bodies in the presence of TTX (Supplementary Fig. 7c). To provide confirmation of the approach at the BNST CRF neuron population level, CRF-Cre mice received bilateral intra-BNST injections of the calcium sensor GCaMP6s (AAV4-Syn-FLEX-GCaMP6s; UPenn Vector Core; 500 nL) and an intra-PVT injection of the excitatory hM3D DREADD virus (AAV2-CaMKIIα-hM3D(Gq)-mCherry, 200 nL) or control virus (AAV8-CaMKIIα-mCherry; Supplementary Fig. 7d). Four weeks later, fresh brain slices were acutely prepared as described for slice electrophysiology recordings above. GCaMP6s was excited with 35% 470 nm LED (CoolLED) at a frequency of 1 Hz for 10 s every min to minimize photobleaching, and videos were acquired at a frame rate of 10 Hz with an optiMOS monochrome camera (QImaging, Surrey, British Columbia, Canada) and MicroManager 1.4 software, across the entire experiment including a five min baseline, 10 min bath application of CNO (10 μM), and 10 min washout period. The experimenter maintained an objective to focus on the BNST z-plane of interest containing several CRF GCaMP6s neurons throughout the experiment. Custom code (MATLAB 9.4) was used to analyze changes in fluorescence intensity in individual CRF GCaMP6 neurons (ΔF) compared to background fluorescence within the frame (intensity of the entire field of view, F) throughout the video. CRF GCaMP6s cells from CON virus mice were used to quantify the inherent linear decay in fluorescence of GCaMP6s, and fluorescence in hM3D CRF GCaMP6s cell fluorescence was normalized to this decay for statistical analysis.
Statistical analyses. Statistical analyses were performed in GraphPad Prism 8. Data for all dependent measures were examined for their distributions in normal and log space, outliers, and equality of variance across groups using Q-Q plots. Electrophysiological properties (except I h and membrane potential) and synaptic transmission data were lognormally distributed, analyzed using log-transformed values, and presented as raw values on a log2 scale in figures; all other data were normally distributed, analyzed in raw space, and presented on a linear scale in figures. Outliers according to Q-Q plots were excluded (however this was rare and reported). Data are presented as mean ± SEM, and raw data points are included in all figures except those with more than two repeated measures when there were too many raw data points to be clearly represented.
Two-way analysis of variance (ANOVA) and unpaired t-tests were used to evaluate the effects of sex and alcohol on synaptic transmission and excitability data. Repeated measures ANOVAs (RM-ANOVAs) were used to examine the effects of treatment, cycles of binge drinking, different anatomical subregions, etc., within individual animals across experimental groups; mixed-effects models were used when one or more matched data points was unavailable for an individual animal. To prevent false-positive results and overinterpretation of RM-ANOVAs on behavioral data, sphericity was not assumed and the Geisser and Greenhouse correction of degrees of freedom was employed for the repeated measure. For all ANOVAs, significant effects were further probed with appropriate post hoc paired or unpaired t-tests with Holm-Sidak correction for multiple comparisons, and multiplicity-adjusted P values are reported. Differences in proportions between groups were assessed with Fisher's exact tests, and differences within the group from a null hypothesis value were evaluated using one-sample t-tests. Statistical comparisons were always performed with an alpha level of 0.05 and using twotailed analyses.