Introduction

Dopamine neurons of the ventral tegmental area (VTADA) and their release of dopamine in the nucleus accumbens (NAc) play key roles in the encoding of primary reward (Cohen et al., 2012; Roitman et al., 2008; Romo and Schultz, 1989), reward-related learning (Schultz, 2016; Schultz et al., 1997; Watabe-Uchida et al., 2017) and motivation (Berridge and Robinson, 1998; Bromberg-Martin et al., 2010; Wise, 2004). Brief, phasic dopamine responses are evoked by primary rewarding stimuli and develop to their predictors (Konanur et al., 2023, for example). Pauses in dopamine activity are evoked by reward omission (Sugam et al., 2012; Tobler et al., 2003). Collectively, these findings support a role for dopamine in signaling reward prediction errors (RPEs) (Schultz, 2016, 1998; Schultz et al., 1997; Steinberg et al., 2013). This value-learning account has been more recently challenged to include the learning from, and updating of, the value of reward-directed actions (Coddington et al., 2023; Jeong et al., 2022). Yet there is little debate that dopamine increases to reward and reward-related cues. A role for dopamine in aversion and aversive conditioning is more controversial (Morales and Margolis, 2017). Some recordings of VTADA activity (Bromberg-Martin et al., 2010; Matsumoto and Hikosaka, 2009) or dopamine release (Mikhailova et al., 2019; Wenzel et al., 2015) support increased signaling in aversion and aversive conditioning while others support decreased signaling (Mileykovskiy and Morales, 2011; Roitman et al., 2008; Ungless et al., 2004; Wheeler et al., 2011; Zhuo et al., 2024; see Morales and Margolis, 2017 for review).

The unsettled nature of dopamine signaling in aversion could be due to several factors. First, like reward (Berridge and Robinson, 2003), aversion is a multi-dimensional construct that ranges from involuntary responses to noxious primary stimuli to passive or active avoidance in response to predictors of aversion. Aversive stimuli have included foot/tail shocks (de Jong et al., 2019; Mileykovskiy and Morales, 2011), air puffs (Matsumoto and Hikosaka, 2009; Mirenowicz and Schultz, 1996; Zhuo et al., 2024), painful pinches (Romo and Schultz, 1989), and white noise (Goedhoop et al., 2022). These are transduced along different sensory pathways (i.e., pain, somatosensory, auditory) from reward (typically gustatory; e.g. sucrose). Even when administering rewarding and aversive stimuli that utilize the same sensory pathway, as with taste, rewarding sweet and aversive bitter solutions activate different taste receptors (Schier and Spector, 2019). Complicating the picture further, mesolimbic dopamine responses to noxious stimuli have been assayed in anesthetized animals (Brischoux et al., 2009; Budygin et al., 2012; Romo and Schultz, 1989) – making behavioral assessment of aversion impossible. In awake, behaving recordings, behavioral paradigms often permit the avoidance of aversive stimuli (Kutlu et al., 2021; Oleson et al., 2012) which can ultimately be rewarding. To facilitate the sampling of aversive stimuli, water or food deprivation protocols are often used (Glover et al., 2016; Gordon-Fennell et al., 2023; Hurley et al., 2023; López et al., 2023; Miranda et al., 2023). However, dopamine responses are modulated by physiological state (i.e., hunger, thirst, sodium appetite (Cone et al., 2015; Fortin and Roitman, 2018; Hsu et al., 2020)), presenting a further confound in interpreting mesolimbic responses.

To overcome these limitations, we measured VTADA neural activity and dopamine release in the lateral shell of the NAc using fiber photometry as we varied the valence of intraoral sucrose infusions. In Paired rats, the valence of sucrose was changed by subsequently administering lithium chloride (LiCl) – which induces visceral malaise (Bernstein et al., 1992) and reliably conditions a taste aversion (CTA; Garcia et al., 1955; Garcia and Kimeldorf, 1957; Nachman and Ashe, 1973; Nolan et al., 1997; Kim et al., 2010; Swank and Bernstein, 1994; Thiele et al., 1996). In Unpaired rats, sucrose valence was unchanged by injecting LiCl 24hr after intraoral sucrose (Smith and Roll, 1967). Changes in valence were strengthened through multiple pairings of sucrose and LiCl or weakened through extinction. Throughout, rats were fed and watered ad libitum to avoid physiological need as a confound. To assay the stereotypical appetitive and aversive behavioral responses to taste stimuli (Breslin et al., 1992; Grill and Norgren, 1978a), we used a deep learning algorithm (DeepLabCut, (Mathis et al., 2018)) and home cage sucrose preference task to assess the valence of intraoral sucrose. Across testing, VTADA activity and dopamine release in the NAc differentially responded to sucrose taste based on its affective value ascertained from behavioral reactivity and sucrose preference. Thus, we conclude that the mesolimbic dopamine system differentially encodes valence – reward versus aversion – and flexibly for the same stimulus.

Results

NAc dopamine differentially responds to primary taste stimuli and correlates with differential behavioral reactivity

We expressed the dopamine fluorescent sensor AAV1.Syn.Flex.GRAB_DA2h (GRAB_DA2h) in the lateral shell of the NAc and recorded real-time fluorescence via an indwelling fiber optic (Figure 1A-B, Supplementary Figure 1). To characterize dopamine and behavioral responses to innately appetitive and aversive taste stimuli, we made intraoral infusions of sucrose and quinine in naïve rats. On average, intraoral infusions of sucrose and quinine evoked different dopamine release responses (t(5) = 3.61, p<0.05; Fig. 1C), which is consistent with prior work from our group (Hsu et al., 2020; Roitman et al., 2008). Across individual trials, responses to sucrose infusions were higher than those to quinine (Fig. 1D). Using Receiver Operating Characteristic (ROC) analysis (Cone et al., 2015; Green and Swets, 1974), we found dopamine responses to sucrose and quinine delivery on individual trials to be highly discriminable (area under the ROC curve (AUCROC)=0.81; Fig. 1E).

In vivo fiber photometry in the mesolimbic dopamine system captures phasic dopamine responses to primary taste stimuli.

(A) Representative images of GRABDA_2h and Cre-dependent GCamp6f expression in the NAc and VTA, respectively. Top Row: Dopamine release recordings from the lateral subregion of the NAc shell are confirmed via viral expression of dopamine-sensor, GRABDA2h (green), probed against DAPI (blue) to visualize sensor location targeted to the lateral shell of the NAc (NAcLS), which borders the NAc core (NAcC). Bottom row: Dopamine cellular activity was recorded in TH Cre+ rats in the VTA (paranigral nucleus (PN), the parabrachial pigmented area; (PBP)). TH+ (red) colocalized with intracellular calcium-sensor, GCaMP6f (green) to demonstrate isolation of the VTADA population. (B) Real-time dopamine release in the NAc across processing steps from a representative rat receiving 5s 0.3M sucrose intraoral infusions (light green bars). Fluorescence excited by the 465nm (Ca2+- and GRABDA2h-dependent, green), and 405nm LED (Ca2+-independent, purple) was captured. 465 and 405nm traces were scaled and subtracted to remove motion artifacts and photobleaching (blue). Fluorescence was then normalized to the whole recording session and represented as a z-score (black). (C) Top: heat maps for NAc dopamine release for 30 trials/session in which 200µL of 0.3M sucrose (green) or 0.001M quinine (orange) was delivered over 5s. Bottom: Average dopamine release aligned to the onset of intraoral delivery. Dotted line represents onset of infusion and gray shading reflects infusion duration and time window for statistical analysis. Bottom inset: Average z-score during infusion period. (D) Relative frequency histogram of dopamine responses (mean z-score) to sucrose and quinine from every trial reported. (E) ROC of D determined a discriminable difference between dopamine responses. Dark lines and shading in C are represented as means; *p<0.05 in paired t-test.

Appetitive and aversive tastes arouse well-characterized, stereotypical responses in rats (Grill and Norgren, 1978a). Specifically, intraoral infusions of appetitive taste stimuli evoke mouth movements and tongue protrusions, but the head and body are relatively still. In contrast, aversive taste stimuli cause gapes, headshakes, forelimb flails, and chin rubs (Grill and Norgren, 1978a) – and therefore greater head and whole-body movement. We hypothesized that these responses could be well captured by measuring movement of the nose and forepaws (Fig. 2A). Here, primary taste stimuli produced distinct nose and forepaw responses to intraoral infusions relative to a 5s baseline period just before infusion onset (behavioral reactivity, Fig. 2B). Quinine produced greater average nose (t(5) = 3.69, p<0.05, Fig. 2C); and forepaw (t(5) = 3.82, p<0.05; Fig. 2D) movement relative to sucrose in the same rats. Importantly, these behavioral measures were negatively correlated to dopamine responses from the same rats where greater dopamine responses were associated with less movement evoked by the stimulus (nose movement to dopamine: r2=0.44, slope=-0.041; p<0.05); forepaw moment to dopamine: r2=0.44, slope =-0.029; p<0.05; Fig. 2E & 2F).

Aversive taste stimuli are linked to increased movement of nose and forepaws.

(A) Representative image of nose, forepaws, tail base, and chamber legs tracking of a rat in a cylindrical chamber from a below chamber perspective. Positional coordinates of selected features were obtained using model created via DeepLabCut, an open-source deep-learning pose estimation program. Custom MATLAB scripts were used to analyze movement from the positional data. (B) Representative movement of nose and forepaws tracked by a DeepLabCut model during 5s pre-infusion, sucrose or quinine infusion, and post-infusion periods. (C&D) Behavioral reactivity was measured as the average change in nose movement or forepaw movement from baseline to infusion period. Intraoral infusion of quinine produces a greater behavioral reactivity. (E&F) Relationship between average change in behavioral reactivity and mean z-score of NAc dopamine during infusion averaged by session. Data in C & D are represented as means; *p<0.05 in paired t-tests. Line in E & F denote the linear relationship between parameters with dotted lines as 95% confidence intervals. P-value of linear regressions indicate slope’s deviation from zero.

Dopamine responses differentially encode the same taste stimulus based on behavioral reactivity

CTA robustly shifts behavior to the same taste stimulus from ingestion to aversion (Grill and Norgren, 1978b). Here, we measured either dopamine release or VTADA activity (Supplementary Fig. 2) and behavioral responses from initially naïve rats (Conditioning Day, CD) to intraoral infusions of 0.3M sucrose. We then administered either Saline (Unpaired) or malaise-inducing LiCl (Paired). The subsequent day, rats did not receive intraoral infusions but did receive the counterbalanced injection in their home cage and were untreated the following day (Fig. 3A). All rats then received another session of sucrose intraoral infusions (Test Day, TD). Thus, by TD, all rats had equal exposure to all stimuli with the lone difference being that Paired rats had sucrose and LiCl administered in close temporal proximity whereas Unpaired rats had sucrose and LiCl administered at least 24hr apart. On TD, we measured each subjects’ dopamine and behavioral responses to intraoral sucrose as per CD (Fig. 3A). In Unpaired rats, dopamine release to intraoral sucrose was unchanged between days (Unpaired: t(9) = 0.83, p>0.05; Figure 3B). In sharp contrast, in Paired rats NAc dopamine release evoked by intraoral sucrose was significantly suppressed on TD relative to CD (Paired: t(10) = 3.97, p<0.01; Fig. 3B). Using an ROC to test for between subject differences, we found that dopamine release to intraoral sucrose on individual trials from Paired and Unpaired rats could not be discriminated on CD (AUCROC=0.50) but could be well discriminated on TD (AUCROC=0.80; Fig. 3C-F).

Pairing of LiCl-induced malaise to sucrose suppresses phasic dopamine responses to intraoral sucrose delivery.

(A) Schematic of Single-pairing CTA timeline. Subjects were first habituated to 30 brief intraoral infusions (200µL/trial) of water at varying intertrial intervals (35-55s) daily for two days. On Conditioning Day (CD), rats received intraoral infusions of 0.3M sucrose and were then injected i.p. with saline (Unpaired) or 0.15M LiCl (Paired). Rats received the counterbalanced injection in the home cage the next day and were untreated the following day. On Test Day (TD), rats received intraoral infusions (parameters identical to CD). (B) NAc dopamine across trials and sessions before and after intraoral sucrose delivery onset. Top: NAc dopamine release (in z-score) across 30 trials/session during both CD and TD. Bottom: Average traces of dopamine release aligned to the onset of intraoral delivery of sucrose on CD and TD. Insets: Mean z-score quantified during 5s intraoral delivery of sucrose. (C-D) Relative frequency histogram of dopamine release responses to sucrose on CD and TD for every trial reported as mean z-score for both Unpaired and Paired subjects. (E) ROC of relative frequency distributions of mean z-score acquired from each trial of on CD and TD between treatment groups. (F) Plotted area under the curve (AUCROC) values of E. (G-K) Recordings of VTADA activity from Unpaired and Paired rats reported with same conventions as B-F. Data in B & G are represented as means; *p<0.05, ** p<0.01, paired t-test.

Similar to dopamine release in the NAc, VTADA activity was unchanged between days in Unpaired rats (Unpaired: t(5) = 0.60, p>0.05; Fig 3G) but suppressed on TD relative to CD in Paired rats (Paired: t(6) = 2.627, p<0.05; Fig 3G). Likewise, individual VTADA responses between groups could not be discriminated on CD (AUCROC=0.55) but were well discriminated on TD (AUCROC=0.76; Fig. 3H-K).

We analyzed behavioral reactivity to sucrose infusions on CD and TD in Unpaired and Paired rats during which dopamine release or VTADA activity were also measured. For rats in which NAc dopamine release was measured, intraoral infusions of sucrose elicited different nose movement responses across testing days (Treatment (Unpaired vs Paired) x Day (CD vs TD): F(1,18)=6.61, p<0.05; Fig. 4A). In the Unpaired rats, nose movement reactions were comparable between days (p>0.05; Fig. 4A) but Paired subjects showed increased nose movement to intraoral infusions from TD relative to CD (p=0.0009; Fig. 4A). Analysis of forepaw movement also indicated differences between Treatment groups (main effect of Treatment: F(1,18)=8.81, p<0.01; Fig. 4B). Forepaw movement did not change from CD to TD in Unpaired rats (p>0.05; Fig. 4B) but increased in Paired rats (p<0.05; Fig. 4B). Further, Paired rats showed greater forepaw reactivity than Unpaired rats on TD (p<0.01; Fig. 4B). Examining the relationship between dopamine release and behavioral reactivity revealed negative correlations between dopamine release and nose movement (r2=0.26, slope =-0.025, p<0.001; Fig. 4C) and between dopamine release and forepaw moment (r2=0.22, slope=-0.021, p<0.005; Fig. 4D).

Suppressed dopamine responses correlate with enhanced behavioral reactivity to intraoral sucrose delivery after CTA formation.

(A&B) Behavioral reactivity was quantified as the change in movement from baseline to infusion period. In Unpaired subjects, the average behavioral reactivity of nose (left) or forepaw (right) movement did not change from CD to TD. In Paired subjects, behavioral reactivity increased from CD to TD. (C) Relationship between behavioral reactivity of nose movement and mean z-score of NAc dopamine during infusion averaged by session. (D) Relationship between average change in behavioral reactivity of forepaw movement and mean z-score of VTADA activity responses during infusion averaged by session. Data in A, B, E & G are represented as means; *p<0.05, **p<0.01; ***p<0.005, Two-Way RM ANOVA with Uncorrected Fisher’s LSD post-hoc. Line in C, D, F & H denote the linear relationship between parameters with dotted lines as 95% confidence intervals. P-value of linear regressions indicate slope’s deviation from zero.

Similar results were observed in the separate group of rats from which VTADA activity recordings were made. Intraoral infusions of sucrose elicited different nose movement responses between subjects across testing days (Treatment x Day: F(1,10)=6.60, p<0.05; Fig. 4E). Specifically, nose movement was comparable between CD and TD in Unpaired rats (p>0.05; Fig. 4E) but significantly greater on TD relative to CD in Paired rats (p<0.005; Fig. 4E). Forepaw movement responses between subjects and across testing days were similar (Treatment x Day: F(1,10)=6.34, p<0.05; Fig. 4F). Forepaw movement did not change across testing days in Unpaired rats (p>0.05; Fig. 4F) but increased in Paired rats (p<0.005; Fig. 4F). Further, dopamine activity was also negatively correlated to nose movement (r2=0.48, slope =-0.057, p<0.0005; Fig. 4G) and forepaw moment (r2=0.51, slope=-0.044, p<0.0001; Fig. 4H).

Repeated ‘safe’ re-exposure to sucrose reduces the suppression of sucrose-evoked phasic VTADA activity in parallel with behavioral extinction

A CTA can be extinguished when the conditioned stimulus (CS (i.e., intraoral sucrose)) is repeatedly presented in the absence of the unconditioned stimulus (US (i.e., malaise induced by LiCl injection) (Hadamitzky et al., 2015)); however, it is unclear if CTA-suppressed dopamine responses to intraoral sucrose recover with extinction. Here, VTADA activity was recorded (Supplementary Fig. 3) in a separate cohort of subjects that underwent Single-pairing CTA followed by five consecutive Extinction sessions, where intraoral sucrose was presented without subsequent injection (Fig. 5A). To examine voluntary sucrose preference or avoidance, rats were given access to sucrose and water for two hours following each intraoral sucrose session (Fig. 5A). In a within-subjects comparison, intraoral sucrose infusions produced comparable VTADA activity responses in Unpaired rats across CD and all Extinction days (E1-E5, Unpaired: main effect of session: F(2.60, 15.59)=2.06, p>0.05; Fig. 5B & 5C). In contrast, VTADA responses to intraoral sucrose were influenced by CTA formation and extinction (Paired: main effect of Day: F(2.18, 15.27)=10.32, p>0.005; Fig. 5B & 5C). CTA formation suppressed sucrose-driven VTADA responses (Paired: C1 vs. E1, p<0.005; Fig. 5B & 5C) but VTADA responses were recovered to CD levels during with re-exposure (Paired: C1 vs. E2, p>0.05; Fig. 5B & 5C). To assess between subject differences across Days, we used an ROC to compare the distribution of VTADA responses on individual trials between Treatments (Unpaired vs Paired). VTADA responses to intraoral sucrose were discriminable on E1 (AUCROC = 0.76) and E2 (AUCROC = 0.77; Fig. 5E-F) but not during E3-E5 (AUCROC ≤ 0.63).

Sucrose exposure under extinction conditions ameliorates CTA’s suppression of phasic dopamine responses to intraoral sucrose.

(A) Schematic of the CTA paradigm consisting of a single-pairing of LiCl or saline pairing to sucrose followed by five consecutive sessions intraoral sucrose not subject to additional US (malaise) exposure. CTA training as in Figure 3 was conducted followed by five Extinction sessions (E1-E5), in which rats received intraoral infusions (parameters identical to TD) without any additional injections. (B) VTADA across trials and sessions before and after intraoral sucrose delivery onset. Top: VTADA activity (in z-score) across 30 trials of the CD and E1-E5 sessions. Bottom: Average traces of VTADA activity aligned to the onset of intraoral delivery of sucrose. (C) Mean z-score during 5s intraoral delivery of sucrose. (D) Relative frequency histogram of VTADA responses to sucrose on all test days for every trial reported as mean z-score for both Unpaired and Paired subjects. (E) ROC of relative frequency distributions of mean z-score acquired from each trial of on all test sessions between treatment groups. (F) Plotted area under the curve (AUCROC) values of E. Data in C are represented as means; ** p<0.01, One-Way RM ANOVA with Dunnett’s multiple comparisons test post hoc.

Preference scores also reflected changing responses to sucrose. Paired rats showed an avoidance of sucrose relative to Unpaired rats (main effect of group F(1,13)=35.24, p<0.0001; Fig. 6A) and sucrose preference scores of Paired rats changed across Extinction days (Treatment X Day: F(4,52)=33.14, p<0.0001; Fig. 6A). Relative to Unpaired rats, Paired subjects had a significantly lower preference for sucrose on E1 (p<0.0001; Fig. 6A) and E2 (p<0.0001; Fig. 6A), which recovered to Unpaired levels by E3 (p>0.05; Fig. 6A). Importantly, sucrose preference was predicted by the earlier recording of VTADA responses to intraoral sucrose (r2=0.38, slope=0.33, p<0.0001; Fig. 6B).

Suppressed dopamine responses to sucrose predict conditioned taste avoidance.

(A) Each Extinction session was followed by a two-hour Two-Bottle Preference test with access to both sucrose and water. Sucrose preference scores were calculated as the percent of sucrose solution consumed from the sum of sucrose and water consumed. Paired (red, closed circles) rats showed decreased sucrose preference on E1 and E2 relative to Unpaired rats (blue, open circles). (B) Average VTADA responses during intraoral sucrose delivery were positively correlated to sucrose preference scores calculated from the Two-Bottle Preference test. Data in A are represented as means; ****p<0.001, Two-Way RM ANOVA with Šidák multiple comparisons test post hoc. Line in B denotes the linear relationship between parameters with dotted lines as 95% confidence intervals. P-value of linear regression indicates slope’s deviation from zero.

To confirm the recovery of VTADA responses after CTA extinction was due to presentations of intraoral sucrose without subsequent US exposure (malaise) and not simply the passage of time, we recorded VTADA activity from a separate group of rats using the Single-pairing CTA paradigm but with a delayed TD (Supplementary Fig. 4). Here, rats were re-exposed to intraoral sucrose seven days after CD (Fig. 7A); thus, the timing of TD was equated with E5 from the Single-pairing paradigm with Extinction (Fig. 5A). Intraoral sucrose evoked comparable VTADA responses in Unpaired rats between CD and TD (Unpaired: t(4) = 0.97, p>0.05; Fig. 7B). In contrast, dopamine responses were significantly suppressed in Paired rats on the delayed TD relative to CD (Paired: t(4) = 10.24, p=0.0005; Fig. 7B). Further, while individual VTADA responses to intraoral sucrose in Unpaired and Paired rats could not be discriminated on CD (AUCROC = 0.59), they were well discriminated on TD (AUCROC = 0.84, Fig. 7C-F).

Delayed testing after sucrose pairing to LiCl-induced malaise suppresses dopamine response to sucrose.

(A) Schematic of the Delayed-test CTA paradigm. CTA training as in Figure 3 was conducted. To match the timeline of Single-pairing CTA with Extinction, rats had five additional untreated days followed by Test Day (TD), which corresponded to the delay from CD to E5. On TD, rats received intraoral infusions (parameters identical to E1). (B) VTADA across trials and sessions before and after intraoral sucrose delivery onset. Top: VTADA activity (in z-score) across 30 trials during CD and TD sessions. Bottom: Average traces of VTADA activity aligned to the onset of intraoral delivery of sucrose. Inset: Average z-score quantified during 5s intraoral delivery of sucrose. (C-D) Relative frequency histogram of dopamine release responses to sucrose on CD and TD for every trial reported as mean z-score for both Unpaired and Paired subjects. (E) ROC of relative frequency distributions of mean z-score acquired from each trial of on CD and TD between treatment groups. (F) Plotted area under the curve (AUCROC) values of E. Data in B are represented as means; ***p<0.005, paired t-test.

Suppression of phasic dopamine responses in CTA scales with conditioning and extinction

To determine if phasic dopamine responses scale with the strength of a taste aversion, VTADA activity was monitored in a separate cohort of rats undergoing a Repeated-pairing CTA with Extinction paradigm (Supplementary Fig. 5). Here, three CD cycles with accompanying noncontingent injection and ‘off’’ days (C1-C3) were administered to Unpaired and Paired rats followed immediately by Extinction days E1-E8. In addition, sucrose preference tests were conducted after each Extinction day (Fig. 7A). In Unpaired rats, average VTADA activity was comparable across all Conditioning (C1-C3) and Extinction days (E1-E8, Unpaired: main effect of day: F(3.88, 23.30)=1.99, p>0.05; Fig. 8B-D). In contrast, Paired rats showed significant changes in VTADA responses to sucrose across the paradigm (Paired: main effect of day: F(2.86, 19.98)=16.07, p<0.0001; Fig. 8B-D). In Paired rats, VTADA responses were suppressed on C3 through E4 (Paired: Relative to C1: C3-E4 p<0.05; Fig. 8B-D) but not subsequent Extinction days (Paired: Relative to C1: E5-E8 p>0.05; Fig. 8B-D). Between subject differences were also evident; responses to intraoral sucrose were discriminable between Unpaired and Paired groups on C3 (AUCROC = 0.88) and E1-E6 (AUCROC >0.7; see Fig. 8E-G for individual AUCROCS).

Phasic dopamine responses to intraoral sucrose scale to the strength of the CTA.

(A) Schematic of the Repeated-pairing CTA paradigm. The three-day CTA conditioning process (as in Figure 3) was repeated three times (C1-C3). Following conditioning, subjects received eight daily Extinction sessions (E1-E8) with intraoral infusions of sucrose followed by no US (malaise) exposure. After the behavioral sessions, rats had access to sucrose and water consumption for two-hours (Two-Bottle Preference test). (B) VTADA across trials and sessions before and after intraoral sucrose delivery onset. Top: Heat maps for VTADA activity (in z-score) for 30 trials in each session during the three conditioning sessions (C1-C3). Bottom: Associated average traces of VTADA activity aligned to the onset of intraoral delivery of sucrose. (C) Conventions as in B but for all eight non-reinforced sucrose sessions (E1-E8). (D) Average z-score quantified during 5s intraoral delivery of sucrose for all conditioning and subsequent test sessions. (E) Relative frequency histogram of VTADA responses to sucrose on all test days for every trial reported as mean z-score for both Unpaired and Paired subjects. (F) ROC of relative frequency distributions of mean z-score acquired from each trial of on all conditioning (left) and extinction (right) days between treatment groups. (G) Plotted area under the curve (AUCROC) values of F. Data in D are represented as means; * p<0.05, One-Way RM ANOVA with Dunnett’s multiple comparisons test post hoc.

Behavioral reactivity was simultaneously recorded with VTADA activity. Intraoral sucrose delivery produced comparable behavioral responses in Unpaired rats across all days (Unpaired: main effect of day: F(3.19, 19.12)=0.74, p>0.05; Fig. 9A). In contrast, Paired subjects showed shifts in behavioral responses across conditioning and extinction days (Paired: main effect of day: F(4.02, 28.15)=7.70, p=0.0003; Fig. 9A). Relative to the first conditioning day, head movement was elevated from C3 through E3 and decreased with subsequent extinction sessions (Paired: Relative to C1: C3, E2, E3 p<0.05, E1 p<0.005; Fig. 9A). VTADA activity and behavioral reactivity were negatively correlated across all rats and days (r2=0.43, slope=-9.52, p<0.0001; Fig. 9B).

Intraoral sucrose driven dopamine responses are negatively correlated to behavioral reactivity and predict sucrose preference.

(A) Head movement to intraoral sucrose averaged across trials. Top: Unpaired rats showed comparable average change in head movement behavioral reactivity across testing sessions. Bottom: Paired subjects showed an increase in average head movement responses to intraoral sucrose from C1 to C2-E3. (B) Relationship between average head movement behavioral reactivity and VTADA responses (in mean z-score) during infusion averaged by session. (C) After each E1-E8 session, rats were administered the Two-Bottle Preference test. Paired (red, closed circles) rats showed decreased sucrose preference on E1-E5 relative to Unpaired rats (blue, open circles). (D) Average VTADA responses during intraoral sucrose delivery were positively correlated to sucrose preference scores calculated from the Two-Bottle Preference test. Data in A are represented as means; *p<0.05, **p<0.01, One-Way RM ANOVA with Dunnett’s multiple comparisons test post hoc. Data in C are represented as means; *p<0.05, **p<0.01, ***p<0.005, ****p<0.001, Two-Way RM ANOVA with Šidák multiple comparisons test post hoc. Line in B & D denotes the linear relationship between parameters with dotted lines as 95% confidence intervals. P-value of linear regression indicates slope’s deviation from zero.

After each extinction session, we conducted a sucrose preference test. Paired rats avoided sucrose compared to Unpaired rats (main effect of group: F(1,13)=59.78, p<0.0001; Fig. 9C). Specifically, Paired rats had a lower sucrose preference on E1 through E5 (E1&E2 p<0.0001, E3 p<0.005, E4 p=0.0001, E5 p<0.05; Fig. 9C). Similar to behavioral reactivity but with more predictive power, VTADA responses were tightly correlated with post-session sucrose preference (r2=0.59, slope=0.48, p<0.0001; Fig. 9D).

Discussion

Unexpected rewards and reward-predictive cues evoke phasic dopamine activity and NAc dopamine release (Schultz, 1998). Responses to aversion remain more controversial (McCutcheon et al., 2012; Morales and Margolis, 2017) and may differ based on sensory transduction pathways (McCutcheon et al., 2012) and whether the aversive outcome can be avoided (Goedhoop et al., 2022; Oleson et al., 2012) – which contribute to competing hypotheses that phasic dopamine responses signal either valence (reward versus aversion) or salience (reward and aversion (Kutlu et al., 2021)). Here, we circumvent common obstacles to determining responses to aversive stimuli by capturing real-time dopamine responses to the same taste stimulus – sucrose – and modulating its valence through conditioning. Intraoral delivery of sucrose eliminated avoidance as a potential confound. Measurement of behavior ensured that rats were responding to the intraoral infusion, eliminating differences in stimulus salience as a contributing factor. A key aspect of our paradigms is that all rats had identical exposure to sucrose and illness with the temporal proximity of sucrose and illness being the lone difference. Ultimately, we show that dopamine responses to sucrose scale with valence across acquisition and extinction of a CTA.

Phasic dopamine activity differentially encodes the valence of primary appetitive and aversive taste stimuli

Dopamine responses to primary aversive stimuli are mixed. Noxious stimuli have been reported to increase the firing rate of dopamine cell bodies in the VTA (Anstrom et al., 2009; Anstrom and Woodward, 2005) and dopamine release in the NAc (Budygin et al., 2012; Kiyatkin, 1995; Mikhailova et al., 2019; Young, 2004; Young et al., 1993). Conversely, others have reported that noxious stimuli inhibit VTADA firing (Brischoux et al., 2009; Ungless et al., 2004) and NAc dopamine release (Mantz et al., 1989; Wenzel et al., 2015) or that responses vary depending on the VTADA subpopulation or dopamine release site (Brischoux et al., 2009; de Jong et al., 2019; Zweifel et al., 2011). Aversive stimuli can drive excitatory inputs to the VTA (Amo et al., 2024; Faget et al., 2024). With respect to taste stimuli specifically, some reports indicate decreased (Roitman et al. 2008; Twining et al., 2015), while others report increased (Bassareo et al., 2002; Kutlu et al., 2021) dopamine responses to bitter quinine. Here, we used intraoral delivery of appetitive and aversive taste stimuli to circumvent confounds including active or passive avoidance, and food or water deprivation. Our findings are similar to prior work from our group (Hsu et al., 2020; McCutcheon et al., 2012; Roitman et al., 2008): intraoral sucrose increased NAc dopamine release while the response in the same rats to quinine was significantly lower. Neuroanatomical data support projections to the VTA from relatively early nodes of the gustatory pathway (e.g. parabrachial nucleus; PBN) and specifically sucrose and quinine responsive neurons in the PBN (Boughter et al., 2019). It remains to be determined if sucrose- and quinine- responsive PBN neurons synapse onto different VTA cell types (e.g. dopamine neurons for the former and GABA neurons for the latter).

Differences in stimulus salience could contribute to the dopamine responses observed here. Indeed, recent work has supported dopamine alignment with the encoding of stimulus salience (Chen and Bruchas, 2021; Kutlu et al., 2021), to contribute to prioritizing attentional processes towards stimuli of high intensity. Further, human VTA activation is greater for novel, more highly salient stimuli than to familiar stimuli (Bunzeck and Düzel, 2006). Here, behavioral responses to quinine were greater than those for sucrose supporting quinine as a highly salient stimulus. Nonetheless, the strong inverse relationship between behavioral reactivity and dopamine release observed here, and particularly the lack of a dopamine response to quinine refutes the idea that phasic dopamine encodes the salience of the stimulus.

Phasic dopamine activity differentially encodes appetitive versus aversive sucrose

The differential dopamine response to appetitive and aversive tastes supports affective encoding by the mesolimbic system. However, sucrose and quinine act on different taste receptors (Schier and Spector, 2019) and may access the VTA via parallel pathways (Boughter et al., 2019). We therefore employed a CTA paradigm to measure dopamine responses to the same stimulus when the stimulus was appetitive versus aversive. While dopamine responses were identical across groups on CD (when rats were initially naïve to intraoral sucrose), pairing intraoral sucrose infusions with LiCl injection caused changed, suppressed dopamine responses to intraoral sucrose when it was next administered (TD). Importantly, the temporal relationship between illness and initial sucrose exposure, as in seminal studies (Smith and Roll, 1967), determined both the formation of a CTA (as evidence in behavioral reactivity measures) and suppressed dopamine activity and release. Indeed, LiCl injection administered 24hr after intraoral sucrose (Unpaired) failed to modulate dopamine signaling. Recent work has suggested that dopamine release from terminals within the NAc could result from intra-NAc processes that are independent of cell body activity (Cachope and Cheer, 2014; Mohebi et al., 2023; Threlfell et al., 2012). However, we found similar modulation of dopamine responses in Paired rats when recordings captured NAc dopamine release or VTADA cell body activity. Our data are more consistent with the good agreement between dopamine cell body activity, release, and post-synaptic signaling reported by others (Lee et al., 2020). Local GABA neurons innervate and suppress VTADA neurons through neurotransmission at GABAA receptors, which can pause VTADA activity and downstream dopamine release in the NAc (Mathon et al., 2005; Tan et al., 2012; Van Zessen et al., 2012). Understanding local GABA dynamics will be an intriguing direction to identify mechanisms for suppression of dopamine signaling to innately aversive tastes as well as those that have acquired aversion through conditioning.

Throughout our studies, we found an inverse relationship between dopamine release/activity and behavioral reactivity to intraoral infusions. Decreases in dopamine activity can lead to an uncoupling of dopamine and high affinity dopamine D2-like receptors expressed on medium spiny neurons (Dreyer et al., 2010). In turn, this could contribute to the different patterns of NAc activity reported for rewarding versus aversive primary (Roitman et al., 2005) and conditioned tastes (Roitman et al., 2010; Wheeler et al., 2008). As manipulations of NAc activity can influence taste reactivity (Morales and Berridge, 2020), differential dopamine responses to sucrose in rats with versus without CTA likely applies a critical filter for NAc processing in the service of oromotor output.

In our purely classical conditioning CTA paradigm, taste (i.e., intraoral sucrose) served as a conditioned stimulus for visceral malaise. Other aversive conditioning paradigms have shown development of decreased dopamine responses to cues predictive of primary aversive stimuli (Goedhoop et al., 2022; Zhuo et al., 2024). There are unique aspects of CTA worth emphasizing. As we’ve shown here, CTA can form with just one taste-illness pairing. In addition, CTA tolerates relatively long delays between the conditioned (taste) and unconditioned stimulus (LiCl). While LiCl was injected immediately after the conditioning session, this interval exceeded 30 minutes from the initial taste exposure. Moreover, indices of the malaise induced by LiCl do not emerge for many minutes, e.g. ∼5min for lying on belly (Aguilar-Rivera et al., 2020), ∼20min for pica (Aguilar-Rivera et al., 2020), and ∼20min for decrease in core body temperature (Bernstein et al., 1992; Guimaraes et al., 2015). Thus, taste and illness memory traces must be integrated over many minutes. The PBN is essential for CTA formation (Spector et al., 1992) and contains neurons that process gustatory and illness information (Carter et al., 2015). The amygdala is also critical for CTA acquisition (Gao et al., 2023; Inui et al., 2019; Morin et al., 2021; Morris et al., 1999; St. Andre and Reilly, 2007). The same neurons within the basolateral nucleus of the amygdala are responsive to both novel taste and LiCl (Barot et al., 2008; Zimmerman et al., 2024). The rostral tegmental nucleus (RMTg) is a robust source of GABAergic input onto VTADA neurons and suppresses their firing (Balcita-Pedicino et al., 2011; Bourdy et al., 2014). Primary aversive stimuli, including LiCl, increase RMTg neuronal firing (Li et al., 2019). Further, re-exposure to LiCl-paired saccharin increases cFos expression in the RMTg (Glover et al., 2016). Likewise, LiCl injections suppress electrically-evoked NAc dopamine release (Fortin et al. 2016). How taste and illness traces and their integration ultimately modulate the mesolimbic system is an important ongoing direction.

Suppression of dopamine signaling scales with CTA expression and extinction

Multiple conditioning trials caused greater suppression in dopamine activity evoked by the sucrose CS. This has been previously reported for audiovisual cues that predict discrete, noxious stimuli (e.g. air puff (Zhuo et al., 2024) and foot shock (Wilkinson et al., 1998)). Here, greater suppression occurred with each taste-illness pairing over just three conditioning trials spaced days apart. Taste, relative to other sensory stimuli, serves as a uniquely strong predictor for illness (Garcia and Koelling, 1966). After establishing CTA, repeated re-exposure to the conditioned stimulus without subsequent malaise drives CTA extinction (Cantora et al., 2006; Mickley et al., 2004; Nolan et al., 1997), allowing for new, contextual knowledge to conflict with the original association between the taste and negative post-ingestive outcome (Bouton and Bolles, 1979). We found that the suppressed dopamine response to sucrose gradually returned to conditioning day levels after repeated extinction sessions. The number of extinction sessions needed for the restoration of dopamine signaling scaled with the number of conditioning sessions (i.e., two extinction sessions for one pairing and five extinction sessions for three pairings). The return of dopamine responses to conditioning day levels was not simply due to the passage of time but rather required CS re-exposure. As with all experiments herein, dopamine responses were negatively correlated with behavioral reactivity across acquisition and extinction. Importantly, dopamine responses to intraoral sucrose were also predictive of home cage sucrose preference across extinction. The infralimbic (IL) cortex plays a role in extinction learning (Quirk et al., 2006). Stimulation of the IL suppresses aversive taste reactivity to sucrose after it had been paired with LiCl (Hurley and Carelli, 2020). Thus, while initial taste-LiCl associations may be due to subcortical (e.g. PBN, RMTg, amygdala) influences on sucrose-evoked VTADA activity and NAc dopamine release, extinction of the suppressed dopamine response to sucrose may well engage top-down cortical processes.

Our studies support a role for dopamine in differentially encoding taste valence. It is important to recognize that there is considerable heterogeneity among dopamine neurons (Phillips et al., 2022). Indeed, VTA neurons that release both glutamate and dopamine are important for aversion, and it is the dopamine signal from these neurons that is critical for this process (Warlow et al., 2024). The technical approaches used here to measure VTADA activity and dopamine release in the NAc average responses across an unknown number of neurons. Thus, while our results suggest that most dopamine neurons decrease their activity and release of dopamine in the NAc to primary and conditioned aversive taste, subpopulations likely play different and critical roles in appetitive and aversive responses. Combining molecular markers and defining input/output relationships for VTA dopamine neurons will be critical for further advancing an understanding of dopamine and aversion. A range of clinically used drugs – from anti-obesity treatments to chemotherapeutics – induce malaise and can lead to the development of CTAs and a reduction in motivation. Thus, further understanding of how these associative processes influence dopamine signaling is important to yield therapeutic approaches with greater efficacy and improved adherence.

Materials & methods

Subjects

Male and female (intact, naturally cycling) Long Evans rats (>250g) were bred from females heterozygous for expressing Cre recombinase under the control of the tyrosine hydroxylase promoter (TH:Cre) (Rat Research Resource Center, RRRC No. 659; Witten et al., 2011) and male wildtype rats (Charles River Laboratories). Offspring were genotyped for Cre (Transnetyx, Inc.) via tissue from ear punches. TH:Cre+ (n=31 males, 22 females) rats were used for experiments where VTADA activity was measured. TH:Cre-(n=13 males, 14 females) rats were used for experiments where NAc dopamine release was measured. Rats were individually-housed in a temperature- and humidity-controlled room and on a 12:12h light:dark schedule (lights on 0700h). All experiments were conducted during the light cycle. Food and water were provided ad libitum. All studies were conducted in accordance with the National Institutes for Health Guide for the Care and Use of Laboratory Animals and approved by the Animal Care Committee at the University of Illinois Chicago.

Viruses

Adeno-associated viruses (AAVs) packaged with fluorescent sensors were used for in vivo fiber photometry. AAV1.hSyn.Flex.GCaMP6f.WPRE.SV40 (5 x 1012GC/mL, Addgene) was used to detect calcium (Ca2+, Chen et al., 2013). AAV9.hSyn.GRAB_DA2h (1 x 1013GC/mL, Addgene) was used to detect dopamine (Sun et al., 2020).

Surgeries

Rats were anesthetized with isoflurane, injected with an analgesic (1mg/kg meloxicam, subcutaneous), and placed in a stereotaxic instrument. To record VTADA activity, 1µL of AAV1.hSyn.Flex.GCaMP6f.WPRE.SV40 was injected into the VTA (Anterior/Posterior (AP): - 5.40mm, Medial/Lateral (ML): +0.70mm, Dorsal/Ventral (DV): -8.15mm, relative to Bregma) of TH:Cre+ rats at a rate of 0.1μL/min. A 5-min post-injection period to allow for diffusion followed. A fiber optic (flat 400-μm core, 0.48 numerical aperture, Doric Lenses Inc.) was then implanted just dorsal (DV: -8.00mm) to the virus injection. To measure dopamine release, 1µL of AAV9-hsyn-GRAB_DA2h was injected into the NAc lateral shell (AP: +1.5mm, ML: +2.5mm, DV: -8.0mm relative to Bregma) of TH:Cre-rats. A fiber optic was implanted just dorsal to the virus injection (DV: -7.9mm). For all surgeries, an intraoral catheter was also implanted. Catheters, made from ∼6cm length of PE6 tubing (Scientific Commodities, Inc.) threaded through a Teflon washer, were inserted via hypodermic needle just lateral to the first maxillary molar. The needle was guided subcutaneously and exteriorized out of the incision at the top of the head. Implants were cemented to skull screws using Metabond (Parkell, Inc.) and dental acrylic. To allow for recovery and construct expression, experiments began at minimum 21-days after surgery (Fig. 1A).

Fiber photometry recording

Light-emitting diodes (LEDs; Doric Lenses) emitted 465nm (Ca2+ or dopamine-dependent) and 405nm (Ca2+-independent) wavelengths. At 30µW emission power, frequencies of 465nm and 405nm light were sinusoidally modulated at 210Hz and 330Hz, respectively. Light was coupled to a filter cube (FMC4, Doric Lenses) and converged into an optical fiber patch cord mated to the fiber optic implant of the rat. Fluorescence was collected by the same fiber/patch cord and focused onto a photoreceiver (Visible Femtowatt Photoreceiver Model 2151, Newport). A lock-in amplifier and data acquisition system (RZ10X; Tucker Davis Technologies), was used to demodulate the fluorescence due to 465nm and 405nm excitation. Transistor-transistor logic (TTL) signals from the system operating behavioral events (Med-Associates, Inc.) were sent to the data acquisition system, timestamped, and recorded along with fluorescence using software (Synapse Suite, Tucker Davis Technologies).

After acquisition, data were Fourier-transformed and the 405nM excitation signal was subtracted from the 465nM excitation signal to account for movement artifacts and photobleaching (ΔF/F). The subtracted signal was smoothed using a custom fifth order bandpass Butterworth filter (cutoff frequencies: 0.05Hz, 2.25Hz) and returned to the time domain using custom scripts written in MATLAB (version R2022b, Mathworks; code available on Github). To compare changes in fluorescence across recording sessions and across rats in different experimental groups, the processed signal for each recording session was normalized to the session’s average fluorescence and converted to z-scores. The normalized signal was then aligned to events of interest (Fig. 1B).

Intraoral Infusion Protocols

Before each session, intraoral catheters were flushed with distilled water to ensure patency. All training and experimental sessions took place in standard chambers (ENV-009A-CT, Med Associates Inc.) or custom-made cylindrical chambers for behavioral reactivity recordings. At the start of a session, a fluid line from a syringe containing a test solution was gravity-fed into a 2-way solenoid valve (The Lee Company) and the fluid line exiting the solenoid was connected to the intraoral catheter. A TTL from Med Associates software and hardware triggered the opening and closing of the solenoid valve. Rats were habituated to receive 30 brief (5s, 200µL) intraoral infusions of water repeated at varying intertrial intervals (35-55s) for two days. Intraoral delivery ensured that rats tasted the stimulus but were not obligated to consume it. This session structure was used throughout all experimental conditions detailed below.

Primary taste

Rats received intraoral infusions of either 0.3M sucrose or 0.001M quinine and received the other solution the following day with solution order counterbalanced across rats (Fig. 1E).

Single-pairing CTA

After habituation, all rats underwent a three-day training cycle. The first day was Conditioning Day (CD), where rats received intraoral infusions of 0.3M sucrose and were randomly assigned to either the Paired or Unpaired group. Immediately after the CD session, rats were injected intraperitoneally (i.p.; 20mL/kg) with either 0.15M LiCl for Paired rats or 0.15M NaCl for Unpaired rats. On the second day, rats did not receive intraoral infusions. Instead, they were injected with the counterbalanced treatment and returned to the home cage. On the third day, rats were left undisturbed in their home cage. This three-day CTA training cycle was followed by the Test Day (TD), in which rats received a session of intraoral infusions of sucrose with no subsequent injections.

Single-pairing CTA with Extinction

Treatment of Paired and Unpaired animals was the same as Single-pairing CTA above except that TD served as the first (E1) of five extinction sessions (E1-E5). For each extinction day, a session of intraoral infusions was administered as on TD. Extinction sessions were administered on consecutive days.

Delayed-test CTA

The same procedure as the Single-pairing CTA paradigm was used. However, TD was not performed until seven days after CD, to match the interval from CD to E5 from the Single-pairing CTA with Extinction condition.

Repeated-pairing CTA with Extinction

Paired and Unpaired animals underwent three conditioning cycles (C1-C3) of the procedure described in Single-pairing CTA. Following the conclusion of the third cycle, extinctions sessions were administered for 8 consecutive days (E1-E8).

Two-Bottle Preference Test

After each Extinction session, rats were returned to the home cage and given 2-hour access to two sipper bottles, one containing 0.3M sucrose and the other water. Bottle locations were switched after 1-hour to avoid a side bias. Change in bottle weight from before to after 2-hour access was used to calculate sucrose preference as sucrose consumed (g) / total fluid consumed (g).

Behavioral Recordings

Behavior in all sessions was video captured at 10 fps for later analysis. In the Primary taste and Single-pairing CTA conditions, rats were recorded in a cylindrical chamber with a clear floor using an adjacent video camera (Teledyne FLIR, Model: BFS-U3-13Y3M-C) that captured behavioral reactivity via the reflection from an angled mirror positioned under the chamber (below view, Grill and Norgren, 1978a). In the studies of Single- and Repeated-pairing CTA with Extinction and Delayed-test CTA, rats were recorded in standard Med Associates chambers with a camera (Wo-We Webcam 720P) positioned above the chamber (above view). Video was monitored using recording software (below view: Teledyne FLIR, Spinnaker SDK; above view: Synapse, Tucker Davis Technologies) and stored for off-line analysis. A custom Med Associates program, linked to the camera’s operation, sent TTLs to time stamp the onset and offset of the video recordings on the fiber photometry software (Synapse, Tucker Davis Technologies) to later sync video to photometry recordings.

Behavioral reactivity was calculated as the change in movement from the 5s pre-infusion period to the 5s intraoral infusion. Positional coordinates for body parts of the rat were obtained using the open-source deep-neural network toolbox, DeepLabCut (Mathis et al., 2018). A custom DeepLabCut model was developed and implemented to track the nose, forepaws, and chamber legs (for scale) in the below view recordings. For above view recordings, a different DeepLabCut model was developed to track the ears and chamber corners. Custom MATLAB codes were employed to calculate the midpoint between the ears (reported as head position) and track distance moved during specified time points, using the coordinates from stationary landmarks (i.e., chamber legs and corners) to normalize distances from pixels to metric units. Any coordinates with model likelihood estimates below 95% were excluded (≤9% of timepoints) and replaced with the linear interpolation of coordinates at neighboring high-confidence timepoints.

Immunohistochemistry

Following completion of experiments, rats were deeply anesthetized with isoflurane and transcardially perfused with 0.9% NaCl followed by 10% buffered formalin solution (HT501320, Sigma Aldrich). Brains were removed and stored in formalin and switched to 20% sucrose the following morning. All brains were sectioned at 40μm on a freezing stage microtome (SM2010R, Leica Biosystems). Sections were collected and processed to fluorescently tag GFP (as an indicator of GCaMP6f or GRAB_DA2h expression) and/or TH via immunohistochemistry. Primary antibodies were incubated at 4°C (washes and other steps at room temperature). Tissues were washed with 1X Potassium Phosphate Buffered Saline (KPBS) six times for ten-minutes, permeabilized in 0.3% Triton-X 100 for 30-minutes and blocked in 2% normal donkey serum for 30-minutes. Sections were incubated in rabbit anti-TH (AB152, Sigma Aldrich) and/or chicken anti-GFP (AB13907, Abcam) antibodies overnight (∼18-hours). After six ten-minute 1xKPBS washes, secondary antibody (Cy3 conjugated donkey anti-rabbit and AF488 conjugated donkey anti-chicken; Jackson Immunoresearch) was applied and sections were incubated for 2-hours at room-temp followed by a single wash. Sections were then mounted onto glass slides, air dried, and coverslipped with Fluoroshield with DAPI (F6057, Sigma Aldrich). Data from rats with GFP expression and fiber placements within the borders of the VTA or NAc (Paxinos and Charles Watson, 2007) were included in analyses.

Data Analyses

To quantify results from in vivo fiber photometry experiments, the mean z-score during the 5s infusion period was measured on each trial and averaged across trials for each session. One- or two-way repeated measures analysis of variance (ANOVA) or paired t-tests were used for statistical comparisons. When group main effects were found with two or more treatments, Tukey’s (for one-way ANOVAs) and Sidak’s (for two-way ANOVAs) post hoc tests were employed. Linear regression was used to calculate p values, r2 goodness-of-fit, 95% confidence bands of the best-fit line, and linear equations for average dopamine response vs. average behavioral reactivity or sucrose preference score. A p value of < 0.05 was used to determine statistical significance. These statistical analyses were performed using GraphPad Prism 10.0 Software (GraphPad Software Inc.).

Finally, we performed ROC analysis (Green and Swets, 1966; Cone et al., 2015) to quantify how discriminable neural responses were to different test conditions via R Studio using the pROC package (package version 1.18.5, RStudio Team (2020), code available on GitHub (Robin et al., 2011)). For each trial, we measured the mean z-score during the 5s intraoral infusion period. The analysis shows the ability of an ideal observer to classify a neural response as driven by sucrose/Unpaired infusion (compared with quinine/Paired). We plotted the rate of hits as a function of false alarm rate across a range of thresholds determined as the mean of the difference between consecutive values of observed z-scores. We then computed the area under the ROC curve (AUCROC) and used a threshold of AUCROC > 0.7 to indicate that an ideal observer could reliably discriminate between conditions.

Acknowledgements

We gratefully acknowledge the help of Drs. Ted Hsu and Vaibhav Konanur for technical support in fiber photometry training and analysis.

Additional information

Author Contributions

Designed research: MKL, SJH, MFR; Performed research: SJH, MKL, PB, & RD; Data Analysis and Visualization: MKL, RD, ATK, JDR, MFR; Wrote the paper: MKL & MFR.

Ethics Declarations

The authors declare no competing interests.

Funding Sources

NIH T32DK128782 (MKL) and R01DA025634 (MFR) were awarded and used to fund experimental procedures. Funding sources were not involved in the study design, data collection or interpretation, or influence journal choice.

NAc lateral shell recording placements from Primary Taste experiment.

The fiber optic implant locations in the NAc lateral shell were verified via post-experimental histological reconstruction. Measurement values denote anterior-posterior coronal plane in relation to Bregma. Reconstructed images do not represent the specific hemisphere where fiber optics were implanted.

NAc lateral shell and VTA recording placements from Single Paring CTA experiment.

The fiber optic implant locations in the NAc lateral shell and VTA were verified via post-experimental histological reconstruction. Measurement values denote anterior-posterior coronal plane in relation to Bregma. Treatment group is indicated by outline color: Unpaired (blue outline) and Paired (red outline). Reconstructed images do not represent the specific hemisphere where fiber optics were implanted.

VTA recording placements from Delayed Test CTA experiment.

The fiber optic implant locations in the VTA were verified via post-experimental histological reconstruction. Measurement values denote anterior-posterior coronal plane in relation to Bregma. Treatment group is indicated by outline color: Unpaired (blue outline) and Paired (red outline). Reconstructed images do not represent the specific hemisphere where fiber optics were implanted.

VTA recording placements from Single Paring CTA with Extinction experiment.

The fiber optic implant locations in the VTA were verified via post-experimental histological reconstruction. Measurement values denote anterior-posterior coronal plane in relation to Bregma. Treatment group is indicated by outline color: Unpaired (blue outline) and Paired (red outline). Reconstructed images do not represent the specific hemisphere where fiber optics were implanted.

VTA recording placements from Repeated Paring CTA with Extinction experiment.

The fiber optic implant locations in the VTA were verified via post-experimental histological reconstruction. Measurement values denote anterior-posterior coronal plane in relation to Bregma. Treatment group is indicated by outline color: Unpaired (blue outline) and Paired (red outline). Reconstructed images do not represent the specific hemisphere where fiber optics were implanted.