Viral Tracing Confirms Paranigral Ventral Tegmental Area Dopaminergic Inputs to the Interpeduncular Nucleus Where Dopamine Release Encodes Motivated Exploration

Abstract Midbrain dopaminergic (DAergic) neurons of the ventral tegmental area (VTA) are engaged by rewarding stimuli and encode reward prediction error to update goal-directed learning. However, recent data indicate that VTA DAergic neurons are functionally heterogeneous with emerging roles in aversive signaling, salience, and novelty, based in part on anatomic location and projection, highlighting a need to functionally characterize the repertoire of VTA DAergic efferents in motivated behavior. Previous work identifying a mesointerpeduncular circuit consisting of VTA DAergic neurons projecting to the interpeduncular nucleus (IPN), a midbrain area implicated in aversion, anxiety-like behavior, and familiarity, has recently come into question. To verify the existence of this circuit, we combined presynaptic targeted and retrograde viral tracing in the dopamine transporter-Cre mouse line. Consistent with previous reports, synaptic tracing revealed that axon terminals from the VTA innervate the caudal IPN; whereas, retrograde tracing revealed DAergic VTA neurons, predominantly in the paranigral region, project to the nucleus accumbens shell, as well as the IPN. To test whether functional DAergic neurotransmission exists in the IPN, we expressed the genetically encoded DA sensor, dLight 1.2, in the IPN of C57BL/6J mice and measured IPN DA signals in vivo during social and anxiety-like behavior using fiber photometry. We observed an increase in IPN DA signal during social investigation of a novel but not familiar conspecific and during exploration of the anxiogenic open arms of the elevated plus maze. Together, these data confirm VTA DAergic neuron projections to the IPN and implicate this circuit in encoding motivated exploration.


Introduction
The modulatory neurotransmitter dopamine (DA) plays critical roles in reward, learning, motivation, and action selection (Schultz et al., 1997;Berridge and Robinson, 1998;Floresco, 2015;Arber and Costa, 2022). Despite decades of intense research, the precise regulation and the circuitry architecture of DAergic neurotransmission still remain unclear (Berke, 2018). Growing evidence demonstrate that midbrain DAergic systems are integrated by a spectrum of molecularly, anatomically, and functionally distinct neuron subtypes. In addition, single-cell gene expression profiling (Tiklová et al., 2019;Poulin et al., 2020;Phillips et al., 2022), together with projection-specific functional mapping, support the hypothesis that heterogeneous DA neuronal clusters can influence individual behavioral readouts Morales and Margolis, 2017;Poulin et al., 2018).
Midbrain DA neurons in the ventral tegmental area (VTA) respond to reward (Mirenowicz and Schultz, 1996), reward-predictive cues (Flagel et al., 2011), associative learning (Saunders et al., 2018), as well as salient stimuli, such as novel social investigations (Gunaydin et al., 2014;Solié et al., 2022). In addition, some VTA DA neurons are engaged by aversive stimuli (Matsumoto and Hikosaka, 2009) or during anxiety-related and fear-related behaviors (Zweifel et al., 2011). Most VTA DA neurons send abundant projection-specific outputs to the ventral striatum nucleus accumbens (NAc) region, where they regulate reward-related and aversive processing (Lammel et al., 2012;de Jong et al., 2019), encode saliency (Kutlu et al., 2021), or promote social behaviors (Gunaydin et al., 2014), but whether the same neurons send functional projections to additional areas and how they control emotional and motivational behaviors are not fully understood.
Recent data described a mesointerpeduncular pathway consisting of VTA DAergic neurons that innervate the IPN (Zhao-Shea et al., 2015), a circuit that mediates anxietylike behavior through unique IPN microcircuitry (DeGroot et al., 2020) and that controls the motivational component of familiar social investigations (Molas et al., 2017b). Such cross talk between two adjacent midbrain structures with apparent opposing roles in regulating behavior (Wills et al., 2022) could have important implications for balancing motivational and affective behaviors. However, a recent study excluded the existence of an anatomic connection from the VTA to the IPN (Nasirova et al., 2021). Thus, a comprehensive analysis clarifying VTA DAergic neuron connections to the IPN and elucidating internal signals that trigger DA release in this brain area, would provide valuable insight into VTA DA neuron architecture, as well as intrinsic midbrain DA circuitry function.

Animals
All animal experiments were conducted in accordance with the Guide for the Care and Use of Laboratory Animals provided by the National Research Council, and with approved animal protocols from the Institutional Animal Care and Use Committee of the Institution. C57BL/6J [stock #000664, The Jackson Laboratory (https://www.jax.org/ strain/000664)] and DAT-Cre [stock #006660, The Jackson Laboratory (https://www.jax.org/strain/006660)] mice were bred in the institution animal facility. Cre lines were crossed with C57BL/6J mice, and only heterozygous animals were used for the experiments. Mice of both sexes were used in all experiments. For social experiments, juvenile stimuli always consisted of C57BL/6J mice (4-7 weeks old). Mice were group housed with a maximum of five per cage and were kept on a standard 12 h light/dark cycle (lights on at 7:00 A.M.) with ad libitum access to food and water. Following viral brain injections and recovery overnight, mice for behavior experiments were transferred to the reverse light/dark cycle room (lights on at 7:00 P.M.) for 3 weeks before fiber brain implantations or additional experiments. Mice were single housed for at least 1 week before behavior testing, which was conducted during the dark cycle (8:00 A.M. to 5:00 P.M.).

Stereotaxic surgeries
Briefly, mice (6 weeks old) were deeply anesthetized with a mixture of 100 mg/kg ketamine and 10 mg/kg xylazine (VEDCO) by intraperitoneal injection. Ophthalmic ointment was applied to maintain eye lubrication. The skin of the skull was shaved and disinfected with iodine. Mice were then placed on a heating pad and in a stereotaxic frame (Stoelting), and the skull was exposed by making a small incision with a scalpel blade. Using bregma and l as landmarks, the skull was leveled along the coronal and sagittal planes. A 0.4 mm drill was used for craniotomies at the target bregma coordinates. Microinjections were made by using a gas-tight 33 ga Hamilton 10 ml neurosyringe (catalog #1701RN, Hamilton) and a microsyringe pump (Stoelting). The following coordinates [in mm, from bregma) were used: for NAc: anteroposterior (AP), 1.0; mediolateral (ML), 60.5; dorsoventral (DV), À4.0; for VTA: AP, À3.51; ML, 60.2; DV, À4.2; and for IPN: AP, À3.51; ML, À1; DV, À4.81; and 12°angle. Viral volumes for injections were 300 nl, delivered at a constant flow rate of 30 nl/min. After injection, the needle was left unmoved for 10 min before being slowly retracted. The incision was then closed and held together with Vetbond.
After 3 weeks of recovery from virus injection, mice underwent surgery as described above for the implantation of optic fibers. Optic fiber [core diameter, 200 mm; numerical aperture (N.A.), 0.48; Doric Lenses] was placed targeting the IPN (AP, À3.8 mm; ML, À1 mm; DV, À4.61 mm; at 12.5°) and was held in place with adhesive luting cement (C&B Metabond, Parkell) followed by dental cement (Cerebond, PlasticsOne). Mice were allowed to recover for 5-7 d in the reverse light/dark cycle room before behavior tests. Injection sites and viral expression were confirmed for all animals by experimenters blinded to behavioral outcome, as previously described (Molas et al., 2017b). Animals showing no viral or off-target site viral expression or incorrect optic fiber placement (,10%) were excluded from analysis.

Fiber photometry and data analysis
Fluorescent signals from biosensors were recorded with a Fiber Photometry System (Doric Instruments). An LED driver was used to deliver excitation light from LEDs at 465 nm (output, ;8.5 mW) and at 405 nm (output, ;5 mW), which was used as an isosbestic wavelength for the indicator (Doric Instruments). The light was reflected into a 200 mm, 0.48 N.A. optic fiber patch cord via the Dual Fluorescence Minicube (Doric Instruments). Emissions were detected with a femtowatt photoreceiver (model 2151, Newport) and were amplified by transimpedance amplification to give an output voltage readout. Sampling (12 kHz) and lock-in demodulation of the fluorescence signals were controlled by Doric Neuroscience Studio software with a decimation factor of 50. A Doric Instruments behavior camera was connected to the Doric Neuroscience Studio software using a USB 3.0 Vision interface to synchronize the photometry recordings with time-locked behavioral tracking systems. All mice were habituated to the patch cord plugged to the optic fiber implant for 10 min in their home cages before the start of the experiment. For social novelty tests, recordings began with the animal in the home cage for 1 min and then placed by the experimenter to the center of the behavioral apparatus. Behavioral events were tallied from the videos in a blinded fashion, and analysis was done using the time-locked photometry recording.
Fiber photometry data analysis was performed using custom-written MATLAB and Python scripts. A low-pass filter (3 Hz) was applied to the demodulated fluorescence signals before the 405 nm channel was scaled to the 465 nm by applying a least mean squares linear regression. Scaled signals were used to calculate the DF/F 0 , where DF/F 0 = (465 nm signalfitted 405 nm signal)/fitted 405 nm signal. The z-scores were calculated using as a baseline the average DF/F 0 values from the À1.0 s before the onset of each behavioral event (considered as time 0, t = 0). For random sampling, two sets of 10 start time stamps were randomly generated, one set within the first 5-149.99 s and the other within 150-294.99 s of the 5 min recording trace. For the 20 random time stamps, DF/F 0 were extracted from À1 to 3 s, and the z scores of each event were estimated using as baseline the À1.0 s before the time stamp.

Behavioral assays
Animals were acclimated to the testing room for 30 min before any experimental assay, and all testing was performed under dim red-light conditions.

Social behavior
Social behavior experiments were performed in wildtype C57BL/6J mice expressing the dLight1.2 biosensor in the IPN. Both male and female mice were used, which interacted with a same-sex C57BL/6J juvenile conspecific. Animals were tested in a Plexiglas apparatus (42 Â 64 Â 30 cm) containing two plastic grid cylinders (diameter, 25 Â 10 cm) located at opposite corners of a rectangular maze. Subject mice were first habituated to the apparatus and the empty cylinders for a 5 min period. Following habituation, a juvenile unfamiliar C57BL/6J conspecific (4-7 weeks of age) was placed inside one of the two cylinders (counterbalanced), reducing social investigations led by the subject animals. The subject mouse was then positioned in the central zone and allowed to freely explore the social and nonsocial cylinders for 5 min. This testing phase was repeated for 24 h, on day 2, using the same juvenile conspecific located in the same compartment, which became familiar. The apparatus and cylinders were cleaned with Micro-90 Solution (International Products Corporation) to eliminate olfactory traces after each session. All sessions were video recorded and synchronized to activity dynamics. Exploration of the social and nonsocial cylinders in videos of the trials were labeled frame by frame by experimenters blind to group conditions. The onset of each behavioral exploratory event (considered as t = 0) was defined whenever the subject mouse directed its nose toward the cylinders at a distance of ,2 cm and initiated a sniffing investigation. Sitting or resting next to the cylinder or objects was not considered exploration.

Elevated plus maze
The elevated plus maze (EPM) apparatus consisted of a central junction (5 Â 5 cm) and had four arms elevated 45 cm above the floor with each arm positioned at 90°relative to the adjacent arms. Two closed arms were enclosed by high walls (30 Â 5 Â 15 cm) and the open arms were exposed (30 Â 5 Â 0.25 cm). A 60 W red fluorescent light was positioned 100 cm above the maze and was used as the illumination source. Both male and female C57BL/6J mice expressing the DA biosensor dLight1.2 in the IPN were used. The optic fiber implant was connected to the recording patch cord, and then mice were placed on the junction part of the maze facing one of the open arms. All mice were given 5 min of free exploration while their behavior was video recorded and synchronized to the dLight1.2 signals via the Doric Instruments fiber photometry system, as described above.

Immunostaining and microscopy
Mice were euthanized by injection of sodium pentobarbital (200 mg/kg, i.p.) and transcardially perfused with ice-cold 0.1 M PBS, pH 7.4, followed by 10 ml of cold 4% (w/v) paraformaldehyde (PFA) in 0.1 M PBS. Brains were postfixed in 4% PFA for 2 h and then submerged in 30% sucrose. Brains were sliced to coronal sections (25 mm) by using a freezing microtome (model HM430, Thermo Fisher Scientific). For virus expression and fiber implant verification, after washes in 0.1 M PBS, sections were mounted, air dried, and coverslipped with Vectashield Mounting Medium (Vector Laboratories). Slices were imaged using a fluorescence microscope (MicroImmagine, Carl Zeiss) connected to computer-associated image analyzer software (release 4.6.1, AxioVision). For immunohistochemical staining, brain sections were permeabilized with 0.2% Triton X-100 in 0.1 M PBS for 5 min, blocked with 2% BSA in 0.1 M PBS for 30 min, and then incubated overnight with the corresponding primary antibodies in 2% BSA at 4°C. The following primary antibodies were used: mouse anti-tyrosine hydroxylase [TH; 1:500; catalog #MAB318, Millipore (https://www.emdmillipore.com/ US/en/product/Anti-Tyrosine-Hydroxylase-Antibody-clone-LNC1,MM_NF-MAB318)]; and guinea pig anti-synaptophysin [1:300; catalog #AGP-144, Alomone Labs (https:// www.alomone.com/p/guinea-pig-anti-synaptophysinantibody/ANR-013-GP]. Slices were subsequently washed in 0.1 M PBS, blocked with 2% donkey (or goat) serum (Sigma-Aldrich) for 30 min and then incubated in secondary antibodies for 1 h [1:800; donkey anti-mouse 647 (catalog #A31571), goat anti-guinea pig 594 (catalog #A11076), Thermo Fisher Scientific]. After washes in 0.1 M PBS, sections were mounted, air dried, and coverslipped with Vectashield medium with DAPI (Vector Laboratories). Images were obtained using a confocal microscope (LSM 700, Zeiss) at 10Â or at 10Â with a 1.5 zoom. Images were analyzed using ImageJ Fiji to create a zoomed-in inset (the red line square on the images, with a 1.5 or 2 zoom factor). The ImageJ JAcoP method was used for colocalization analysis between VTA DA ! IPN fibers expressing AxonGCaMP and synaptophysin staining. Briefly, each image threshold was set automatically for analysis before Mander's coefficient was applied to obtain the fraction of synaptophysin (red) overlapping with VTA DA ! IPN terminals (green) and vice versa. For quantification of fluorescently labeled axons from VTA DA neurons innervating the IPN, the Digital Enhancement of Fibers with Noise Elimination (DEFiNE) method was applied (Powell et al., 2019; available for download at: https://figshare.com/s/1be5a1e77c4d4431769a). Axons were quantified in confocal images that were not processed through the clean images function, but each input image was a single-channel maximum intensity projection. Quantification was performed in ROIs (0.3 Â 0.4 mm) randomly allocated within the anterior (bregma, À3.4 mm) and posterior (bregma, À3.8 mm) IPN.

Statistical analysis
Statistical analyses for fiber photometry were performed using parametric tests on z-scored data after testing for normality. One-way or two-way repeatedmeasures (RM) ANOVA with Dunnett's multiple-comparisons tests or Bonferroni's post hoc tests was conducted for the analyses involving the comparison of group means, as indicated. The z-scores are presented as the mean 6 SEM of all events for transitions between open arms (included junction) and closed arms, and for social approach behaviors. Comparisons of z scores were made using the calculated average for each animal. All analyses were performed using Prism 9 (GraphPad). Statistical significance was accepted at p , 0.05 (see Table 1 for statistics summary).

Data availability
The code used for fiber photometry data analysis is freely available on GitHub (https://github.com/TapperLab/ TapperLab) and also as Extended Data 1.

Results
We used a genetic strategy to target putative DA neuron subtypes and rigorously investigate DAergic projections from the VTA to the neighboring IPN. To this aim, we specifically selected a knock-in genetic mouse line that expresses Cre recombinase under the transcriptional control of the endogenous DA transporter (DAT) promoter. In this mouse line, Cre recombinase expression is driven from the 39 untranslated region of the endogenous DAT gene by means of an internal ribosome entry sequence (IRES) to reduce interference with DAT function (Bäckman et al., 2006). Some neurons in the IPN express Th mRNA, which can lead to recombination in Th-IRES-Cre mice, although these neurons have low/ undetectable TH protein in the adult brain (Poulin et al., 2018). These Th 1 IPN neurons are not related to midbrain DA neurons, as they are not derived from the midbrain floor plane, and they lack the expression of typical DAergic neuronal markers such as DAT, NURR1, FOXA2, or PITX3 (Poulin et al., 2018); therefore, using DAT-Cre mice restricts and minimizes expression to midbrain DA neurons. Previous work expressed a Cre-dependent virus in the VTA of DAT-Cre animals and detected neuronal projections innervating mainly the caudal part of the IPN (cIPN; Molas et al., 2017b;DeGroot et al., 2020). To verify that VTA DA ! IPN projections are indeed axonal terminals and not simply DA dendritic elements extending into the IPN, we injected adeno-associated viruses (AAVs) containing the hSyn.Flex.mGFP.2A.synaptophysin.mRuby construct into the VTA of DAT-Cre mice (Fig. 1A). Following Cre recombination, synaptophysin fused to the mRuby red fluorophore is selectively transported into the axonal compartments of the transfected neurons ( Fig. 1A; Zhang et al., 2016b). TH immunostaining demonstrated efficient recombination restricted to DA neurons in the midbrain (Fig. 1B). Furthermore, via circuit mapping, abundant axon terminals were detected in the NAc region, the principal output target of VTA DA neurons (Fig. 1C). These VTA DA ! NAc axon terminals intensely expressed synaptophysin-mRuby fused protein (Fig. 1C), altogether validating the viral-mediated genetic strategy. To delineate the VTA DA ! IPN circuit, we used a group of six mice, with comparable results. All injected animals reliably exhibited VTA DA synaptophysin-mRuby axon terminals innervating the IPR region of the cIPN (Fig. 1D). Additional VTA DA axonal varicosities were also detected targeting the cIPN IPDM/IPDL subregion (Fig. 1D), consistent with previous data (Molas et al., 2017b;DeGroot et al., 2020).
To reassure that VTA DA neurons send neuronal projections innervating the neighboring IPN and that these are active presynaptic axons, DAT-Cre mice received an injection of Cre-dependent AxonGCaMP in the VTA expressed via AAV5-mediated gene delivery ( Fig. 2A). This genetically encoded calcium indicator is uniformly enriched in axons, allowing for structure-specific labeling of presynaptic terminals . Similarly, as described above, presynaptic terminals from VTA DA neuronal inputs were observed in the IPR and IPDM regions of the cIPN (Fig. 2B-D). Moreover, immunostaining against synaptophysin protein revealed robust colocalization between the GFP 1 (AxonGCaMP) and synaptophysin (Fig. 2B-D), confirming active presynaptic structures.
Distinct VTA DA projection populations regulate reward associations and motivation via specific NAc inputs (Heymann et al., 2020). To elucidate the projection specificity of VTA DA that innervates the cIPN, DAT-Cre mice received a coinjection of AAV2-hsyn-mCherry (localization marker) together with AAVrg-hsyn-DIO-eGFP in the NAc region (Fig. 3A). Imaging of the target injection site confirmed viral-mediated gene delivery restricted mainly to the shell area of the NAc (Fig. 3B). In addition, to verify the retrolabeled neurons detected in the VTA were positive for DAergic markers, brain slices of the injected animals were immunostained against TH protein. All the experimental animals (n = 6 mice) exhibited abundant terminal projections from retrolabeled VTA ! NAc projecting neurons that innervated the IPR region of the cIPN (Fig. 3C). The cell bodies from VTA ! IPN projecting neurons mostly localized in the paranigral (PN) area of the VTA (Fig. 3C) and were indeed DAergic, as shown by colocalization with TH staining (Fig. 3C). For visualization enhancement and quantification of the fluorescently labeled axons, we used the DEFiNE method. Axonal fibers innervating the IPN from retrolabeled VTA DA ! NAc projecting neurons were highly enriched in posterior regions of the IPN (i.e., cIPN) compared with anterior IPN bregma (Fig.  4A,B). Similarly, quantification of axonal fibers originating from direct infusion of the synaptophysin-mRuby construct in VTA DA neurons revealed increased axon terminals at more posterior IPN bregma compared with anterior (Fig. 4C,D). Noticeably, in anterior IPN bregma, the number of axonal fibers was higher when VTA DA neurons were directly transfected with the synaptophysin-mRuby construct as opposed to retrolabeled VTA DA ! NAc projecting neurons (Fig. 4E). In contrast, these two viralmediated VTA DA neuron-labeling strategies resulted in a similar number of axonal fibers at posterior IPN bregma (Fig. 4F).
Previous work suggested that the VTA DA ! IPN circuit is engaged during anxiety-like behaviors (DeGroot et al., 2020) and when mice encounter unfamiliar conspecifics (Molas et al., 2017b). Although DA signals have been detected in acute mouse IPN slices (DeGroot et al., 2020), the real-time dynamics of in vivo IPN DAergic neurotransmission have never been reported. To this aim, here we recorded IPN DA dynamics in freely behaving mice using the genetically encoded DA sensor dLight1.2 (Patriarchi et al., 2018). Fluctuations in IPN DA signals were recorded during the three-chamber sociability task, when mice encountered a new juvenile conspecific ( Fig. 5A; Materials and Methods). On the following day, subject mice were presented to the same juvenile conspecific in the same location, which became familiar ( Fig. 5A; Materials and Methods). To this aim, we virally expressed dLight1.2 in the IPN of C57BL/6J mice, enabling ultrafast optical DA recordings, and at 3 weeks post-viral transduction we implanted an optic fiber targeting the injection site (Fig.  5B). IPN DA dynamics were time locked to when animals approached and initiated a sniffing investigation of conspecific stimuli (Fig. 5C). Demodulated fluorescence signals were obtained from the 465 and 405 nm channels in a 5 min trial (Fig. 5D). The 405 nm channel was scaled to the 465 nm by applying a least mean squares linear regression (Fig. 5E). Scaled signals were used to calculate the DF/F 0 where DF/F 0 = (465 nm signalfitted 405 nm signal)/fitted 405 nm signal (Fig. 5F). On day 1 of the sociability test, sniffing investigation of a novel conspecific significantly increased the release of DA in the IPN (Fig. 5G-I). However, IPN DA signals rapidly habituated on the next session, as the conspecific became familiar (Fig. 5J-L). Random sampling of IPN DA signals without being time locked to social sniffing investigations did not result in apparent changes in activity either when mice interacted with a novel conspecific (Fig. 5M-O) or when this became familiar (Fig. 5P-R).
To further investigate IPN DA signals triggered by additional behaviors, we recorded IPN DA dynamics in mice tested in the EPM (Fig. 6A), a well established  White arrows indicate presynaptic puncta colocalization. Inset, Magnified view of colocalization between eGFP and the synaptophysin marker. Scale bar, 100 mm. C, Quantification of the colocalization coefficient between eGFP and synaptophysin staining from single-plane confocal images containing the cIPN (n = 6 mice; 4 males, 2 females). D, Top, AxonGCaMP expression in the VTA of DAT-Cre mice (eGFP, green), synaptophysin immunostaining (red) and colocalization of the two channels (merge) in brain slices paradigm to measure anxiety-like behaviors in rodents (Walf and Frye, 2007). As mice investigated the open arms of the EPM, the release of DA in the IPN significantly increased (Fig. 6B-E). Conversely, the transition from the open to the closed EPM compartments led to reductions in IPN DA signals (Fig. 6B,F-H). Timelocked IPN DA signals when mice entered the open arms were higher compared with when entering the closed arms of the EPM or with non-time-locked random sampling signals (Fig. 6I-K). All the recorded animals were verified for correct viral expression and fiber placement within the cIPN (Fig. 7).

Discussion
DA dysfunction has been implicated in numerous brain diseases, including addiction, depression, schizophrenia, Parkinson's disease, and anxiety disorders (Horga and Abi-Dargham, 2019;Nestler and Lüscher, 2019;Ruitenberg et al., 2021;Taylor et al., 2022;Zalachoras et al., 2022). A Figure 3. VTA DA neurons from the PN send projections to the IPN. A, Schematic of viral strategy used. DAT-Cre mice were injected with a viral mixture of AAV-hSyn-DIO-eGFP (retrograde) and AAV2-hSyn-mCherry (location marker; 1:1) into the NAc. B, Representative image showing the virus injection site targeting the NAc shell area (AcbSh). Scale bar, 100 mm. C, Example of injected animals, cases 1-6 (4 males, 2 females), all showing retrolabeled eGFP 1 neurons in the VTA colabeled with TH staining. For each case: top, TH immunostaining (magenta), retrolabeled eGFP 1 neurons from the NAc (green), and overlay of the two channels (merge) in brain slices containing the cIPN. Scale bar, 100 mm. Insets: right, a magnified view enclosing the PN and IPR in the merged channel (red square, 2Â zoom in); bottom, enlarged view of the PN and IPR region from the top images with a right inset image of the merge channel demonstrating AcbSh-projecting neurons in the PN region are DAergic (TH 1 ) and also send efferents to the IPR in the cIPN (red square, 2Â zoom in). Scale bar, 100 mm. continued containing the cIPN. Scale bar, 100 mm. Bottom, Enlarged view of the IPR region from the top images (gray square). White arrows denote VTA DA eGFP 1 presynaptic projections in the IPR colocalized with synaptophysin puncta. Scale bar, 100 mm.
comprehensive understanding of the circuit architecture and the functional mapping of DA neurons is imperative to gain insights into inherent regulation of DA neurotransmission in health and disease. The present study confirms the existence of a mesointerpeduncular pathway that connects the VTA with the IPN, thereby modulating behavioral states with implications in overall midbrain DA circuitry function.
the DAT-Cre knock-in mouse line, in which Cre mimics the expression pattern of the plasma membrane dopamine transporter (Bäckman et al., 2006;Lammel et al., 2015) and, therefore, demonstrates higher specificity targeting putative midbrain DA neurons (Poulin et al., 2018). VTA DA axons preferentially innervated the IPR region of the cIPN, as previously reported. These axon terminals were detected in most injected animals across multiple experimental cohorts and appeared to be more obvious in those mice where viral expression extended to the PN region of the VTA. In addition, synaptically targeted markers localized in terminal projections from the VTA DA ! IPN circuit, whereby protein immunostaining revealed active presynaptic terminals rather than passing fibers. Interestingly, the cIPN is highly enriched in neurons expressing the D 1 receptor (Molas et al., 2017b), but also in serotonergic cell bodies (Groenewegen et al., 1986). Serotonergic IPN neurons innervate the ventral hippocampus (vHipp) to mediate active stress coping and natural reward (Sherafat et al., 2020). Considering that cIPN neurons can amplify VTA DA signals through a microcircuit that spans to additional IPN subregions (DeGroot et al., 2020), if Figure 5. Novel social encounters trigger IPN DA signals. A, Schematic of the experimental approach used to measure IPN DA activity during interactions with novel and familiar social stimuli. Subject mice were exposed to the same juvenile C57BL/6 conspecific on days 1 (novel) and 2 (familiar), while IPN DA signals were recorded using the dLight1.2 biosensor. B, Schematic of AAV-dLight viral injection strategy in the IPN of C57BL/6 mice (left) and representative pictograph of DA sensor dLight1.2 (green) expression with optic probe location targeting the cIPN (right). Scale bar, 100 mm. C, Illustration of a social sniffing investigation. D, Example of raw signals (in volts) corresponding to the 465 and 405 nm channels recording during a 5 min interaction with a new social stimulus. E, The 405 nm channel is scaled to the 465 nm by applying a least mean squares linear regression. F, Scaled signals are used to calculate the DF/F 0 value, where DF/F 0 = (465 nm signalfitted 405 nm signal)/fitted 405 nm signal. G, DF/F 0 values time locked to IPN DA signals relative to the initiation of a social sniffing investigation (red line) on day 1, when mice interact with a novel conspecific. H, Heatmap representations (top) and z-score values (bottom) of the time-locked IPN DA signals relative to social novelty explorations. I, Average z score per second compared with the baseline signal from 1 s before the onset of each social sniffing event (preonset, gray). Statistical comparisons were made using an average z score per animal (n = 10 mice; 6 males, 4 females). Significant increases in IPN DA activity were observed 2;3 s after onset of novel social sniffing investigations. One-way RM ANOVA (F (3,39) = 21.80, p , 0.0001). Dunnett's multiple-comparisons test: **p , 0.01, ***p , 0.001. J, DF/F 0 values time locked to IPN DA signals relative to the time initiating a social sniffing investigation (red line) on day 2, when mice interact with a familiar conspecific.  Anatomical and functional connectivity of midbrain DA neurons has been broadly investigated across animal species (Swanson, 1982;Morales and Margolis, 2017). Numerous studies identified the source of synaptic input to DA neurons (Lammel et al., 2012;Watabe-Uchida et al., 2012;Beier et al., 2015), as well as output targets (Lammel et al., 2011(Lammel et al., , 2008Poulin et al., 2018;Heymann et al., 2020). While consistent data indicate that the NAc is the major target of VTA DA neurons, additional structures such as the amygdala, cortex, hippocampus, ventral pallidum, septum, periaqueductal gray, bed nucleus of stria terminalis, olfactory tubercle, and locus coeruleus, among others, also receive DAergic inputs from the VTA. Noticeably, most of the circuit-tracing studies traditionally focus on those regions with the highest abundance of DA terminal projections, neglecting target-specific sites that receive sparse DAergic inputs. For instance, VTA neurons send local, topographically organized axonal connections that innervate the VTA itself (Adell and Artigas, 2004;Ferreira et al., 2008;Aransay et al., 2015), which overall have received less attention. Of note, Aransay et al. (2015) also reported VTA innervation to the IPN, which, although less frequent, nevertheless supports a direct anatomic link between the VTA and IPN.
The anatomic location of DA neuron synaptic output can be a critical factor determining its intrinsic properties and behavioral outcomes (Lammel et al., 2011;de Jong et al., 2019). Our data show that a subpopulation of NAc shell-projecting VTA DA neurons in the PN region may preferentially project into the IPN to innervate cIPN, as reported previously (DeGroot et al., 2020). Emerging evidence suggests that subpopulations of VTA DAergic neurons can innervate more than one brain structure (Aransay et al., 2015). Specifically, medial shell NAc-projecting DA neurons send significant collaterals outside the striatum, including the septum and ventral pallidum, indicating that this DA subpopulation is capable of simultaneously influencing neural activity in multiple brain regions . Since the same DA neurons presumably innervate the IPN, the data together position the IPN as an integral member within specific VTA DAergic subcircuitries.
Our photometry results demonstrate that innate DA signals in the IPN are triggered with motivated exploration, when mice investigate novel conspecific individuals and when they explore the anxiogenic arms of the EPM. These results affirm that social interactions bear rewarding aspects and recruit neural circuits of motivation (Chevallier et al., 2012), including DAergic systems (Gunaydin et al., 2014;Hung et al., 2017;Bariselli et al., 2018;Solié et al., 2022). Given that the NAc shell represents a storage site for social memories (Okuyama et al., 2016), one possibility could be that innate IPN signals contribute to social novelty and familiarity responses, supporting previous findings (Molas et al., 2017b). On the other hand, NAc shell-projecting VTA DA neurons are recruited by aversive stimuli and cues that predict them (de Jong et al., 2019). Increased IPN DA signals with the exploration of anxiogenic environments would result from activation of a neural network that strengthens responses to aversive stimuli to modulate anxiety-like behavior.
A recent study excluded the existence of an anatomic connection from the VTA to the IPN (Nasirova et al., 2021). One possible explanation for the discrepancy in the results may be that most of the viral-mediated circuit tracing in the study by Nasirova et al. (2021) was done in a Cre mouse line that only targets IPN neurons expressing the a5 nAChR subunit. Although neurons in continued multiple-comparisons test: *p , 0.05, ***p , 0.001. I, Example of IPN DA DF/F 0 values time locked to the transition to the open or closed arms of the EPM compared with DF/F 0 values obtained with random sampling across the 5 min recording session. J, The zscore values of I. K, Mean z-score values of the baseline and the 3 s open and closed EPM arm exploratory event for the true signal compared with random sampling signal. Two-way RM ANOVA, significant time Â z score interaction: F (2,36) = 4.14, p = 0.024; p = 0.0001, Bonferroni's post hoc test: **p , 0.001, ***p , 0.001. All data represent the mean 6 SEM. the IPN are highly enriched in a5*-nAChRs (Ables et al., 2017), some subpopulations do not express the a5-encoding gene. Thus, limiting IPN circuit tracing to an a5-expressing neuronal subtype does not accurately reflect total IPN connectivity. In addition, for the viral-mediated retrograde tracing analysis, the authors selected IPN brain slices with a maximum IPN caudal bregma coordinate of À3.6 mm according to the atlas of Paxinos and Franklin (2001;Nasirova et al., 2021). As mentioned above, VTA DA neurons that project to the IPN localize more caudal, at coordinates À3.63 to À4.03 mm from bregma, which were likely missed in the analysis. Noticeably, previous work using rabies tracing from overall IPN neurons did detect sparse cell bodies localized in caudal VTA (Lima et al., 2017). Nasirova et al. (2021) used the Allen Connectivity Atlas to reinforce their negative data. However, the few Allen examples performed in the Slc6a3-Cre (DAT-Cre) line lack viral expression transfecting caudal VTA PN neurons, thereby precluding the detection of any putative VTA DA innervation to the IPN. Additionally, Nasirova et al. (2021) included examples of VTA Cre-mediated anterograde tracing in DAT-Cre mice, but, for this experiment, the authors used a nonvalidated Cre-dependent synaptically targeted GFP marker, which presented strong labeling of cell bodies in the medial mamillary nucleus and also the IPN itself (Nasirova et al., 2021), two brain regions lacking DA neurons, thus raising questions regarding the specificity of the virus and therefore the validity of the results. Surprisingly, the article by Nasirova et al. (2021) failed to cite, consider, or discuss the study by DeGroot et al. (2020), which used a multidisciplinary approach and specifically demonstrated the following: (1) DA detection in IPN slices using a genetically encoded DA sensor; (2) optogenetic activation of VTA DA IPN inputs elicits a postsynaptic response that is blocked by a D 1 receptor antagonist; (3) retrograde Credependent AAV-eGFP injection into the medial nucleus accumbens shell labels VTA neurons that clearly project into the IPN of DAT-Cre mice (a result that was repeated here with the addition of TH staining to label DAergic neurons); and (4) optogenetic activation or silencing the DAergic IPN input decreases and increases anxiety-like behavior, respectively.
In summary, the present study was able to confirm the existence of a mesointerpeduncular pathway that connects the VTA with the IPN, replicating previous findings (Aransay et al., 2015;Zhao-Shea et al., 2015;Molas et al., 2017b;DeGroot et al., 2020). These results may significantly influence the prevailing models of intrinsic midbrain DA circuitry as well as of IPN function. Considering that VTA DAergic neurons also send projections to the mHb  (Phillipson and Pycock, 1982;Beier et al., 2015), the data together suggest a complex direct dopaminergic modulation of the habenulointerpeduncular tract that may have strong impact on reward-related, aversive/affective motivated behaviors. Finally, beyond the VTA-IPN axis, and bearing in mind that the activation of small subsets of neuronal ensembles can lead to selective widespread activation of neural networks with concomitant behavioral outcome (Marshel et al., 2019;Dalgleish et al., 2020), the present work emphasizes the need of investigating sparse, functionally relevant neglected circuits that may serve as signal amplification to computationally process motivational information.