Cell-Type-Specific Afferent Innervation of the Nucleus Accumbens Core and Shell

Introduction

The nucleus accumbens (NAc), a brain region located in the ventral aspect of the basal ganglia, is widely recognized to act as a limbic-motor interface (Floresco, 2015). The NAc has been implicated in various important functions including reward-related processes (Ambroggi et al., 2008; Carlezon and Thomas, 2009; Day et al., 2011; Russo and Nestler, 2013), spatial navigation (Floresco et al., 1997; Ito et al., 2008), feeding (Krause et al., 2010; Kenny, 2011; Urstadt et al., 2013), and sexual motivation (Everitt, 1990; Beny-Shefer et al., 2017). In addition, dysfunction of this region is thought to be involved in numerous neurological and psychiatric disorders such as depression (Russo and Nestler, 2013; Francis and Lobo, 2017), anxiety disorders (Fu et al., 2017), Alzheimer’s disease (Ferrier et al., 1983; Nie et al., 2017), Parkinson’s disease (Albin et al., 1989; DeLong, 1990), and drug addiction (Russo et al., 2010; Lüscher, 2016). Previous studies have established that the NAc integrates inputs from many brain areas, such as the ventral midbrain (Fallon and Moore, 1978; Phillipson and Griffiths, 1985), basal amygdaloid complex (Wright et al., 1996; Groenewegen et al., 1999), hippocampus (Kelley and Domesick, 1982; Lopes da Silva et al., 1984; Yang and Mogenson, 1984), thalamus (Phillipson and Griffiths, 1985; Berendse and Groenewegen, 1990; Brog et al., 1993), and prelimbic and prefrontal cortex (Beckstead, 1979; Phillipson and Griffiths, 1985; Montaron et al., 1996; Gorelova and Yang, 1997). However, how the diverse functions of the NAc are regulated by these complex anatomical connections remains poorly understood.

The NAc, however, is an inhomogeneous tissue that is often divided into two primary segments: the shell and the core (Zahm and Brog, 1992). The functions of these two subregions differ significantly in many behavioral paradigms. For example, approach toward reward-related stimuli is mediated by the NAc core but not the shell (Parkinson et al., 2000; Di Ciano et al., 2001; Salamone and Correa, 2012), whereas the opposite holds true for learning about the irrelevance of stimuli (Gal et al., 2005; Floresco et al., 2006). Additionally, the projection neurons within the NAc are GABAergic medium spiny neurons (MSNs), which are dichotomized based on their expression of either D1 or D2 dopamine receptors (Gerfen et al., 1990; Kawaguchi, 1997). These two MSN subtypes, which together make up >95% of all NAc neurons, exert balanced but antagonistic roles in reinforcement learning. In visual cue-based reward learning, NAc (core + shell) D1R-MSNs are critical for learning acquisition, whereas D2R-MSNs are needed for the flexibility to the learning switch (Yawata et al., 2012). Opposing roles for these two MSN subtypes in learning is further supported by the finding that optogenetic activation of D2R-MSNs suppresses cocaine reward, with opposite effects induced by activation of D1R-MSNs (Lobo et al., 2010).

Overall, to dissect the diverse functions of NAc, it is crucial to determine the direct synaptic inputs to different anatomical subregions in a cell-type-specific manner. Previous studies using traditional tracers have provided valuable information about the connections of the NAc (Kelley and Domesick, 1982; Berendse and Groenewegen, 1990; Brog et al., 1993; Wright et al., 1996; Gorelova and Yang, 1997). However, technical limitations of these traditional methods make it difficult to map the whole-brain inputs to specific neuron types. Furthermore, traditional tracers may be taken up by fibers that pass by the NAc, resulting in nonspecific labeling. Recently, the development of monosynaptic retrograde transsynaptic tracing based on modified rabies virus (RV) has greatly facilitated the mapping of whole-brain inputs to genetically defined cell types (Wickersham et al., 2007; Watabe-Uchida et al., 2012; Pollak Dorocic et al., 2014; Weissbourd et al., 2014; Do et al., 2016).

In the present study, we utilized such a virus-based labeling technique to identify monosynaptic inputs to D1R- and D2R-MSNs within different NAc subregions. We then compared the input distribution patterns among four groups defined by cell type and location (i.e., D1R-MSNs in the NAc core, D1R-MSNs in the NAc shell, D2R-MSNs in the NAc core and D2R-MSNs in the NAc shell). Although the input distributions of different cell types were highly similar, different NAc subregions had significantly different biased inputs from numerous brain areas. These findings provide an anatomical foundation for future studies of the neural circuits that underlie the diverse functions of the NAc.

Materials and Methods

Animals

All the experimental procedures were approved by the Hubei Provincial Animal Care and Use Committee and the experimental guidelines of the Animal Experimentation Ethics Committee of Huazhong University of Science and Technology. BAC-transgenic D1R-Cre [MMRRC Tg (Drd1acre) EY262Gsat], and D2R-Cre [MMRRC Tg (Drd2-cre) ER44Gsat] mouse lines were acquired from the Mutant Mouse Regional Resource Center (MMRRC; Davis, California). These mice were backcrossed with C57BL/6J mice. C57BL/6J mice were purchased from Beijing Vital River (Beijing). All animals were housed with their littermates in a dedicated housing room under a 12/12 h light/dark cycle, and food and water were available ad libitum. For all experiments, adult male mice (2–3 months old) were used.

Virus Information

All viral tools used in this study were packaged and provided by BrainVTA (BrainVTA Co., Ltd., Wuhan, China). The detailed production and concentration procedures for modified RV were previously described (Osakada et al., 2011). The final titer of RV-EnvA-ΔG-dsRed was 2 × 10⁸ infecting units per milliliter. For adeno-associated viruses (AAVs), AAV9-EF1a-FLEX-EGFP-2a-TVA and AAV9-EF1a-FLEX-RG were packaged into 2/9 serotypes with final titers at 1–5 × 10¹² genome copies per milliliter.

Surgery and Viral Injections

Before virus injection, experimental mice were anesthetized with 2% chloral hydrate, 10% ethylurethane, and 1.7 mg/mL xylazine in a 0.9% NaCl mixture (0.09 ml/10 g body weight, i.p.) and mounted in a stereotaxic holder (item: 68030, RWD, Shenzhen, China). After the skulls of the experimental mice were adjusted to be parallel to the reference panel, a craniotomy (∼0.5 mm diameter) was made above the targeted areas with a dental drill. For cell-type-specific tracing, 25–30 nl of helper AAV (AAV9-EF1a-FLEX-EGFP-2a-TVA and AAV9-EF1a-FLEX-RG mixed at a 1:2 ratio of viral particles) was stereotactically injected into the NAc core (coordinates in mm: AP + 1.34, ML 1.10, DV:-4.3) or the NAc shell (coordinates in mm: AP + 1.34, ML 0.55, DV:-4.7) using a glass micropipette connected to a syringe pump (Item: 53311, Quintessential stereotaxic injector, Stoelting, United States). The micropipette was held in place for 10 min after the injection before being slowly retracted from the brain. Two weeks after helper AAV injection, 150–200 nl of RV-EnvA-ΔG-dsRed was injected into the same location under biosafety level 2 conditions through the same procedure mentioned above.

Histology and Immunostaining

One week after injection of RV, mice were perfused intracardially with 0.01 M phosphate-buffered saline (PBS), followed by 4% paraformaldehyde (PFA) in 0.01 M PBS. Mouse brains were carefully removed and placed in 4% PFA solution overnight for postfixation. The procedures for agarose embedding were generally performed as reported previously (Sallee and Russell, 1993; Jiang et al., 2017). Briefly, after postfixation, each intact brain was washed with PBS for 4 h and dehydrated in 10% sucrose in PBS overnight. The brains were slightly dried and embedded in melted agarose (Sigma-Aldrich, United States) using a silicone mold in a 55°C water bath for 0.5 h. The orientations of the brains were carefully adjusted in the water bath. The mold was then placed at room temperature until the agarose solidified. Subsequently, the brains were separated from the mold and stored in PBS at 4°C until sectioning. Coronal brain slices (50 μm) were prepared using a vibratome (Leica VT1000, Leica Microsystems). Every second section was selected for whole-brain mapping and visualized with a slide scanning microscope (Nikon).

To characterize the starter cells, sections near the injection site were selected for immunofluorescent staining. These sections were blocked with 3% bovine serum albumin (BSA) in PBS-0.3% Triton X-100 for 1 h at room temperature, incubated with the primary rabbit anti-DARPP-32 monoclonal antibody (1:200, Abcam ab40801) overnight at 4°C, washed with PBS (three times for 10 min each at room temperature), incubated with Alexa-Fluor 647 goat anti-rabbit (1:500, Invitrogen A32733) for 2 h at room temperature, incubated with DAPI (1 ng/ml) for 10 min, and finally washed with PBS (five times for 7 min each at room temperature) prior to mounting the slides. We imaged the sections using an LSM 710 inverted confocal microscope (Zeiss).

Image Analysis

The locations of labeled neurons were manually registered, and the number of cells in each region from both hemispheres was quantified using Fiji (NIH) according to a standard mouse brain atlas (Paxinos and Franklin, 2001). The input from each region was normalized by dividing the number of labeled neurons found in that region by the total number of labeled neurons (excluding the NAc) in each brain to obtain the percentage of total inputs. For statistical analyses, two-tailed unpaired t-tests, the Mann–Whitney U test or one-way ANOVA followed by the Bonferroni correction were performed using SPSS (version 13.0). To quantify the similarity in input patterns, we calculated Pearson’s correlation coefficients. Error bars indicate the SEM.

Results

Strategies for Tracing Monosynaptic Inputs to D1R- and D2R-MSNs Within Different NAc Subregions

Our aim is to trace, determine and compare the direct inputs to D1R- and D2R- MSNs in the NAc core and shell (hereafter called D1R-core, D1R-shell, D2R-core, and D2R-shell neurons). We genetically targeted different projection neuron subpopulations in the NAc by using D1R- and D2R-Cre mouse lines, which have been shown to label D1R- and D2R-expressing neuron populations with high specificity (Gong et al., 2007). To demonstrate the whole-brain monosynaptic inputs for each cell type, we used a retrograde transsynaptic tracing system based on a modified RV, SAD-ΔG-dsRed (EnvA) (Wickersham et al., 2007). First, we coexpressed the avian receptor TVA and rabies glycoprotein G (RG) specifically in each cell type, which was achieved through injection of two Cre-inducible AAV vectors (AAV9-EF1a-DIO-eGFP-TVA and AAV9-EF1a-DIO-RG) into D1R- or D2R-Cre mice. Injections were targeted to either the NAc core or the NAc shell. Two weeks later, we injected SAD-ΔG-dsRed(EnvA) into the same area (Figures 1A,B). This modified RV can infect only those cells expressing TVA and requires RG to spread retrogradely to presynaptic cells. Mice were perfused 1 week later; the whole brain was sectioned at 50 μm, and every second section was selected for further analysis.

FIGURE 1

FIGURE 1. Experimental strategy for identification of monosynaptic inputs to D1R- and D2R-MSNs in different NAc subregions. (A) Recombinant AAV strains and RV. (B) Experimental design. (C) Representative coronal brain sections near the injection sites. Blue indicates cell nuclear staining with DAPI. Scale bar, 1000 μm. Inset, enlarged view of the area in white the box showing starter cells (yellow, expressing both eGFP and dsRed, indicated by white arrowheads). Scale bar, 30 μm. AcbC, nucleus accumbens, core (NAc core); aca, anterior commissure, anterior part; AcbSh, nucleus accumbens, shell (NAc shell). (D) Numbers of starter neurons in the individual animals. (E) Proportions of labeled starter neurons in injection regions. LS, lateral septal nucleus; BST, bed nucleus of the stria terminalis. (F) Numbers of transsynaptically labeled neurons, i.e., input neurons in the individual animals. (G) A linear relationship was detected between the number of starter and input neurons.

The starter cells were identified based on the coexpression of TVA-eGFP and dsRed. We observed a large number of double-labeled neurons near the injection sites (Figures 1C,D). Combined with DARPP-32 immunostaining, a reliable marker of both populations of MSNs (Ouimet et al., 1984; Matamales et al., 2009), we found that almost all eGFP- and dsRed-double-positive neurons (starter cells) were also DARPP-32-positive, confirming that they are MSNs, in both D1R-Cre and D2R-Cre mice (Supplementary Figure S1). For each group, the location of the vast majority of starter cells was restricted to the injected subregion of the NAc (Figure 1C), although we found a small number of eGFP- and dsRed-double-positive neurons in neighboring nuclei: the lateral septal nucleus (LS), bed nucleus of the stria terminalis (BST), anterior olfactory nucleus (AON), and caudate putamen (CPu). However, these neurons made up a small fraction of total starter neurons in most animals (Figure 1E). We found 7336–78192 dsRed-labeled presynaptic neurons in each brain (n = 24) (Figure 1F). Variation of the cell numbers across animals was partly due to different injection volumes. Nevertheless, the numbers of transsynaptically labeled neurons had a linear relationship with the numbers of starter neurons (Figure 1G).

A number of control experiments were conducted to define the specificity of our tracing approach. We first injected SAD-ΔG-dsRed (EnvA) into the NAc core of wild-type animals (i.e., Cre-negative) without prior AAV injection. This treatment resulted in no dsRed-labeled neurons (Figure 2A), indicating the dependence of the RV infection on AAV-induced expression of TVA. Next, we repeated the retrograde transsynaptic tracing strategy in wild-type animals. Consistent with previous studies using similar methods (Watabe-Uchida et al., 2012; Weissbourd et al., 2014; Do et al., 2016), slight local background infection by RV was observed (Figure 2B), most likely due to the leaky expression of a low level of TVA. Note that all these non-specifically labeled neurons were found near the injection sites, within the NAc. Thus, to rule out the impact of this nonspecifically labeling on the mapping of long-range inputs, the following analysis excluded data originating from the NAc. Furthermore, to verify that the transsynaptic spread is under the tight control of RG expression, we injected only AAV9-EF1a-DIO-eGFP-TVA and SAD-ΔG-dsRed (EnvA) into the NAc core of D1R-Cre mice. A large number of labeled neurons were found near the injection site, and almost all dsRed-positive neurons coexpressed eGFP (Figures 2C₁,C₂). In contrast, we did not detect any labeled neurons outside the NAc (Figure 2C₁). Together, these results demonstrated that our tracing approach was suitable for mapping the long-range monosynaptic inputs to the NAc and that labeled neurons outside the injection sites represent monosynaptic inputs to projection neurons in the NAc.

FIGURE 2

FIGURE 2. Control experiments for rabies-mediated transsynaptic tracing. (A) Injection of RV without prior AAV injection resulted in no dsRed-labeled neurons, indicating the dependence of the RV infection on AAV-induced expression of TVA. Scale bar, 500 μm. (B) Injecting AAV-DIO-EGFP-TVA, AAV-DIO-RG, and RV into the wild-type mice led to a low level of nonspecific labeled RV neurons that expressed dsRed at the injection site. Scale bar, 500 μm. (C₁) Injection of AAV-DIO-EGFP-TVA and RV without prior AAV-DIO-RG injection in the D1R-Cre mice resulted in no DsRed-labeled input neurons. Scale bar, 500 μm. (C₂) Enlarged view of the boxed region in panel (C₁); all dsRed-positive neurons coexpressed eGFP. Scale bar, 50 μm.

Long-Range Cell-Type-Specific Inputs to the NAc Core and Shell

Representative coronal images from the transsynaptic tracing of inputs to each group show the whole-brain distribution of dsRed-labeled neurons (Figure 3A). We found that NAc (including core and shell) neurons integrate inputs from diverse brain regions, from the anterior neocortex to specific midbrain areas (Figure 3A). The dsRed-labeled neurons were predominantly found ipsilateral to the injection site, although sparser labeling was also observed in the contralateral hemisphere (data not shown). To quantify the distributions of long-range monosynaptic inputs, we counted the number of dsRed-labeled neurons within individual regions of each brain. Since the total number of labeled neurons varied across brain samples, to directly compare between experiments, we normalized the data in each region by the total number of labeled neurons in each brain (Supplementary Table S1).

FIGURE 3

FIGURE 3. Overview of whole brain input to D1R-MSNs and D2R-MSNs in different NAc subregions. (A) Representative coronal sections showing labeling of monosynaptic inputs to D1R-core, D1R-shell, D2R-core, and D2R-shell neurons. Only the side ipsilateral to the injection site is shown. Scale bar, 1 mm. (B) Proportions of total inputs from 11 brain areas. Mean ± SEM (n = 6 mice each for the D1R-core, D1R-shell, D2R-core, and D2R- shell groups). ^∗∗∗p < 0.001, ^∗∗p < 0.01, and ^∗p < 0.05. Only significant differences between the same cell type in different subregions or between different cell types in the same subregions are marked; one-way ANOVA with Bonferroni correction.

When the brain was divided into 11 major brain areas (Figure 3B), the hippocampus contributed most of the long-range inputs to the NAc, followed by the cortex, thalamus, olfactory areas, amygdala, pallidum, hypothalamus, septum/diagonal band of Broca, striatum (excluding the NAc), midbrain, and hindbrain (Figure 2B). Even at this gross anatomical level, the pattern of inputs targeting the NAc core versus NAc shell shows a notable difference, particularly in the cortex and hippocampus (Figure 2B): the cortex preferentially projects to the NAc core, whereas the hippocampus preferentially projects to the NAc shell.

To investigate input tracing in more detail, we further divided the brain into 84 regions. Figure 4 displays the proportion of input neurons in specific regions. Despite the functional differences of D1R-MSNs and D2R-MSNs in the NAc, the overall input distribution patterns of D1R-MSNs (Figure 4, left) and D2R-MSNs (Figure 4, right) were similar regardless of whether they were located in the NAc core or NAc shell. All regions projecting to D1R-MSNs also provided direct inputs to D2R-MSNs and vice versa. For both D1R-MSNs and D2R-MSNs, most of the inputs to the NAc shell were concentrated in a few brain regions (Figure 4, blue), and conversely, inputs to the NAc core were broadly scattered (Figure 4, red). For example, the subiculum (S), paraventricular thalamic nucleus (PV), hippocampal CA1 region and ventral pallidum (VP) were the four major presynaptic regions that projected to the NAc shell, and inputs from these four regions accounted for ∼50% of total inputs to the NAc shell (42.3 ± 7.6% for D1R-MSNs; 57.5 ± 13.0% for D2R-MSNs). These results suggested that both types of MSNs in the NAc might be regulated by many common upstream regions, whereas the NAc core might integrate information from a broader brain area than the NAc shell.

FIGURE 4

FIGURE 4. Quantitative analysis of the proportions of whole-brain input to D1R-MSNs and D2R-MSNs in different NAc subregions. (Left) Monosynaptic inputs to D1R-MSNs in the NAc core (red) and NAc shell (blue). (Right) Monosynaptic inputs to D2R-MSNs in the NAc core (red) and NAc shell (blue). Mean ± SEM (n = 6 mice for each group).

Multiple Nuclei Provide Significant Contrasting Inputs to the NAc Core and Shell

As mentioned above, the overall input distribution patterns of D1R-MSNs and D2R-MSNs in both NAc subregions show a high degree of similarity; therefore, we will consider these two types of neurons in each subregion together in the next analysis.

When the virus injection was targeted to the NAc core, the labeled neurons were widely distributed across cortical areas. In contrast, only sparse labeling was found in cortical areas when we targeted the injection to the NAc shell (Figures 5A,B, 6A,B; also see Figure 4). This difference was especially pronounced in the anterior cortex. For example, we found dense labeling in the medial and ventral orbital cortex (MO/VO), agranular insular cortex (AI) and prelimbic cortex (PrL) when tracing inputs to the NAc core, with significantly fewer labeled neurons when tracing inputs to the NAc shell. In addition, other anterior cortical areas, including the lateral and dorsolateral orbital cortex (LO/DLO), cingulate cortex (Cg) and secondary motor cortex (M2), also preferentially project to the NAc core (Figures 5A,B, 6A,B).

FIGURE 5

FIGURE 5. Monosynaptic inputs to D1R-MSNs from the cortex and thalamus. (A) Distributions of anterior cortex inputs to D1R-MSNs in the NAc core (upper panels) and shell (lower panels). Scale bar, 200 μm. (B) Proportion of total inputs to D1R-MSNs in the NAc core (red) and shell (blue) from 11 cortical areas that contained relatively large numbers of input neurons (>1% in at least one experimental group). Mean ± SEM (n = 6 mice for each group). ^∗∗∗p < 0.001, ^∗∗p < 0.01, and ^∗p < 0.05, two-tailed unpaired t-test or Mann–Whitney U test. (C) Distributions of thalamic inputs to D1R-MSNs in the NAc core (right) and shell (left). Inset, enlarged view of the area in the white box showing spatially separated labeling. Scale bar, 200 μm.

FIGURE 6

FIGURE 6. Monosynaptic inputs to D2R-MSNs from the cortex and thalamus. (A) Distributions of anterior cortex inputs to D2R-MSNs in the NAc core (upper panels) and shell (lower panels). Scale bar, 200 μm. (B) Proportion of total inputs to D2R-MSNs in the NAc core (red) and shell (blue) from 11 cortical areas that contained relatively large numbers of input neurons (>1% in at least one experimental group). Mean ± SEM (n = 6 mice for each group). ^∗∗∗p < 0.001, ^∗∗p < 0.01, and ^∗p < 0.05, two-tailed unpaired t-test or Mann–Whitney U test. (C) Distributions of thalamic inputs to D2R-MSNs in the NAc core (right) and shell (left). Inset, enlarged view of the area in the white box showing spatially separated labeling. Scale bar, 200 μm.

Among the rostral forebrain areas, NAc core injection led to dense labeling in the AON and piriform cortex (Pir), whereas NAc shell injection caused the densest labeling in the tenia tecta (DTT/VTT) (Figure 4). Despite being near the NAc (the injection site), the dorsal striatum, which prominently includes the CPu, receives only a small fraction of the total inputs in both NAc core- and shell-targeted cases. In fact, labeled neurons in this structure may be “spill-over” from the NAc (Figures 3A, 4). That is, labeled neurons were often present at the border between the CPu and NAc, whereas the center of the CPu showed very sparse or no labeling. Conventional tracing studies have identified a strong reciprocal connection between the NAc and VP (Groenewegen et al., 1999). In accordance with these studies, we found that the VP provides many inputs to both the NAc core and NAc shell, with a slight preference for the NAc shell (Figure 4).

In the amygdala, two subdivisions account for nearly all amygdala projections to the NAc: the basolateral amygdaloid nucleus (BL, including the anterior, posterior and ventral parts of the basolateral nucleus, BLA/BLP/BLV) and basomedial amygdaloid nucleus (BM, including the anterior and posterior parts of the basomedial nucleus, BMA/BMP) (Figures 3A, 4). While the BL makes up a similar proportion of the total inputs to the NAc core and shell, the BM makes up a significantly larger proportion of inputs to the NAc core (Figure 4). We found a sparser distribution of labeled neurons in the central nucleus (Ce) and medial nucleus (Me) for both NAc core- and shell-targeted cases (Figure 4). Interestingly, the lateral amygdaloid nucleus (La) shows a notable preference to project to the NAc core (Figure 4).

Consistent with conventional retrograde tracing studies (Brog et al., 1993), we found that the LS, particularly the ventral and dorsal divisions (LSV/LSD), sends prominent monosynaptic inputs to the NAc shell but has very sparse connectivity with the NAc core (Figures 3A, 4). The connectivity between the hippocampus and the NAc has attracted great interest as a circuit involved in spatial navigation and recognition of novelty (Floresco et al., 1997; Ito et al., 2008). Here, we found many inputs from the hippocampus. Indeed, the S contained the largest numbers of total labeled neurons for both NAc core- and shell-targeted cases, although with significantly more projections to the NAc shell (Figures 3A, 4). The NAc shell also received dense inputs from the CA1 region of the hippocampus, whereas the NAc core received relatively few inputs from this region (Figures 3A, 4).

In the thalamus, the projections to the NAc core and shell were roughly spatially segregated. For example, more dorsal and medial structures, such as the PV, particularly the anterior dorsal part, project preferentially to the NAc shell. In contrast, more ventral and lateral structures, such as the paratenial thalamic nucleus (PT) and mediodorsal thalamic nucleus (MD), show a notable trend toward projecting to the NAc core (Figures 4, 5C, 6C). These results correspond well to previous experiments using traditional retrograde tracers to label inputs to different subregions of NAc (Brog et al., 1993). Other thalamic nuclei, including the intermediodorsal (IMD), central medial (CM), paracentral (PC), and centrolateral thalamic nuclei (CL) also send more inputs to the NAc core than to the NAc shell (Figure 4). Additionally, we detected a large number of labeled neurons in the parafascicular nucleus (PF) for both NAc core- and shell-targeted cases (Figures 3A, 4).

Hypothalamic nuclei play a critical role in the regulation of feeding and circadian rhythms (Saper, 2006). Many hypothalamic areas provide scattered inputs to both the NAc core and NAc shell with similar patterns (Figure 4). We found that there was a slight preference for hypothalamic nuclei, such as the lateral preoptic area (LPO), medial preoptic area (MPA), and substantia innominata (SI), to project to the NAc core. In these areas, the lateral hypothalamic area (LH) is unique in that its inputs to the NAc are strongly biased toward the shell. Interestingly, our cell-type-specific tracing revealed that the LH preferentially targets D1R-MSNs over D2R-MSNs in the NAc shell (see section “Discussion”).

There were very few inputs from the midbrain and hindbrain to either the NAc core or the NAc shell (Figure 4). The ventral tegmental area (VTA) and substantia nigra pars compacta (SNc), known to contain a high density of dopaminergic neurons, provided the major midbrain inputs to both NAc subregions. Some labeling was also sparsely distributed in the dorsal raphe nucleus (DR), median raphe nucleus (MnR), rostral linear nucleus of the raphe (RLi), periaqueductal gray (PAG), pedunculopontine tegmental nucleus (PPTg), retrorubral field (RRF), locus coeruleus (LC), and parabrachial nucleus (PB), suggesting rather weak inputs from these structures to the NAc.

In summary, these results show that although inputs to the NAc core and shell originate from largely overlapping areas, many cortical and subcortical nuclei provide significant contrasting inputs to these two subregions.

Comparison Between Inputs to D1R- and D2R-MSNs in Different NAc Subregions

In our analysis of long-range tracing, we found that although the input patterns to D1R-MSNs and D2R-MSNs are highly similar, the input patterns to different subregions of the NAc differ significantly. To validate our previous findings, we quantified the similarities and differences in the distribution patterns of inputs to the four groups (D1R-core, D1R-shell, D2R-core, and D2R-shell neurons). Figures 7A–D show the correlation coefficients for all pairs (Ogawa et al., 2014). Each open circle in the scatter plots represents one brain region (significant differences in red, p < 0.05), and the diagonal represents the same input proportion for each pair. Comparing the input patterns of D1R- and D2R-MSNs, we found that the vast majority of the brain regions, represented by the open circles centered around the diagonals, showed that they provided similar numbers of inputs to D1R-MSNs and D2R-MSNs in both NAc subregions (Figures 7A,B). In the NAc core- and shell-targeted cases, we found that only 2 and 24 of 84 brain regions, respectively, provided significantly different inputs to D1R-MSNs and D2R-MSNs. The correlation coefficients between inputs to D1R-MSNs and D2R-MSNs in both the NAc core and the NAc shell were extremely high (correlation coefficient, r = 0.987 for D1R-core versus D2R-core; r = 0.967 for D1R-shell versus D2R-shell, p < 0.001; Figures 7A,B), confirming their high degree of similarity as observed earlier. On the other hand, when we computed the correlation coefficients between inputs to the same cell type in different NAc subregions, we found relatively low correlation coefficients (r = 0.607 for D1R-core versus D1R-shell; r = 0.602 for D2R-core versus D2R-shell, p < 0.001; Figures 7C,D). The open circles representing brain regions were relatively scattered from the diagonal. The D1R-MSNs in the NAc core and shell received significantly different inputs from 57 of 84 brain regions (Figure 7C). There were 52 of 84 brain regions that provided significantly different inputs to D2-MSNs in different NAc subregions (Figure 7D). The results are presented in Table 1.

FIGURE 7

FIGURE 7. Comparisons of input distributions across the four groups. (A) Comparison between inputs to D1R-core and D2R-core neurons. (B) Comparison between inputs to D1R-shell and D2R-shell neurons. (C) Comparison between inputs to D1R-core and D1R-shell neurons. (D) Comparison between inputs to D2R-core and D2R-shell neurons. Values are the means of the percent of total inputs from each region. Red circles indicate significant differences (p < 0.05, two-tailed unpaired t-test or Mann–Whitney U test). r, Pearson’s correlation coefficient. (E) Summary of selected prominent monosynaptic inputs to the NAc core (top, red) and shell (bottom, blue), considering D1R-MSNs and D2R-MSNs together. Line thickness represents the number of inputs.

TABLE 1

TABLE 1. Brain areas providing significantly different inputs to each comparison pair.

Discussion

Using rabies-mediated transsynaptic tracing, we mapped the whole-brain monosynaptic inputs to the two main projection cell types (D1R-MSNs and D2R-MSNs) in different NAc subregions. Our experiments confirmed many previously demonstrated connections, but with cell-type specificity and quantitative analyses. Our data showed an extremely high degree of overall similarity between inputs to D1R-MSNs and D2R-MSNs in both NAc subregions. However, the input patterns to different NAc subregions differ significantly (Figure 7E). Both D1R- and D2R-MSNs in the NAc core are more likely than those in the shell to receive monosynaptic inputs from the anterior cortex, while D1R-MSNs in the NAc shell receive disproportionally higher input from the LH. Through quantitative comparisons, we identified a list of brain nuclei that provided significantly different inputs to each comparison pair (Table 1).

Accumulated evidence has indicated that D1R-MSNs and D2R-MSNs in the NAc (core + shell) play differential and often opposing roles in reinforcement learning (Lobo et al., 2010; Yawata et al., 2012). A recent study combining an AAV-based anterograde tracing approach and optogenetics demonstrated the classical separation of D1R-MSNs in the direct pathway and D2R-MSNs in the indirect pathway did not necessarily apply in NAc (Kupchik et al., 2015). However, the organizing principles of afferents to these two MSN subtypes in the NAc remain unclear. Here, we found that D1R-MSNs and D2R-MSNs in both NAc subregions received quantitatively similar patterns of inputs, suggesting that these two MSN subtypes integrate information from largely overlapping brain areas. Interestingly, our result is similar to the results of recent transsynaptic tracing studies in the DR (Weissbourd et al., 2014), the VTA (Beier et al., 2015), the basal forebrain (BF) (Do et al., 2016), the arcuate nucleus (ARC) (Wang et al., 2015), the dorsal striatum (Guo et al., 2015), and the medial part of the olfactory tubercle (OT) (Zhang et al., 2017), where different cell types within the same brain structures receive similar overall inputs from diverse brain regions. Given the extensive interconnection between local axonal collaterals of D1R- and D2R-MSNs (Kawaguchi et al., 1989; Taverna et al., 2008; Chuhma et al., 2011), the coinnervation of these two MSN subtypes may reflect complex computations by striatal microcircuits and a need to modulate MSN activity through collateral inhibition (Taverna et al., 2008). It has been reported that D1 and D2 receptors exhibit a marked difference in their dopamine binding affinity (Maeno, 1982; Richfield et al., 1989); thus, phasic and tonic firing of dopamine neurons differentially modulate D1 and D2 receptors, respectively (Grace et al., 2007; Hikosaka, 2007). These characteristic features of dopamine transmission may result in the distinct functions mediated by D1R- and D2R-MSNs despite their shared upstream dopamine pathway (Hikida et al., 2010). However, given the vast differences in the expression profiles of functional proteins and intracellular signaling molecules between D1R- and D2R-MSNs (Gerfen et al., 1990; Grace et al., 2007; Hikosaka, 2007; Heiman et al., 2008), it remains to be determined whether these two MSN subtypes in the NAc are connected by separate subpopulations within the overlapping upstream brain areas. Future research will require the development of new technologies that allow us to label the presynaptic partners of separate cell types with different colors within the same animal.

The transsynaptic tracing strategy allowed us to quantitatively compare inputs to different cell types of different NAc subregions. Although inputs to D1R-MSNs and D2R-MSNs in the NAc show a high degree of similarity, we identified multiple nuclei that provide significant biased inputs to different NAc subregions. The difference observed in the anatomical connections of the NAc core and shell likely underlies their distinct behavioral functions. For example, we found that many anterior cortical areas (including the VO/MO, LO/DLO, AI, PrL, M2, and Cg) were preferentially labeled when starter cells were located in the NAc core but not the shell. Interestingly, most of these structures have long been implicated in value-based decision making (Hauber and Sommer, 2009; Sul et al., 2010, 2011). Thus, our results may provide explanations for the prominent role of the NAc core rather than the shell in instigating approach behavior toward motivationally relevant stimuli.

The NAc shell-to-LH pathway is well documented to play an important role in feeding control (Kelley et al., 2005; Urstadt et al., 2013; O’Connor et al., 2015). For example, pharmacological inhibition of the NAc shell triggers increased feeding that can be prevented by concomitant infusion of a gamma-aminobutyric acid type A (GABAA) receptor agonist into the LH (Urstadt et al., 2013). A recent study showed that D1R-MSNs in the NAc shell provide the predominant source of accumbal input to the LH (O’Connor et al., 2015). In our tracing study, we found that the LH also projects preferentially to the NAc shell, consistent with previous findings (Brog et al., 1993). On the other hand, although the input from the LH is strongly biased toward the NAc shell as opposed to the core, we extend the previous findings by showing that the LH projects preferentially to D1R-MSNs, over D2R-MSNs in the NAc shell. Although we have not yet determined the functional significance of this reverse pathway, our findings provide an anatomical basis that may help elucidate the computational mechanisms underlying the regulation of energy homeostasis.

It is worth mentioning that we detected only a small proportion of total inputs from the VTA and SNc, which contain the major dopaminergic populations, as well as from the DR and MnR, the sites of the major serotonergic populations. These results appear to differ from some traditional tracing studies indicating that the striatum (the CPu and NAc) receives substantial dopaminergic (Fallon and Moore, 1978; Phillipson and Griffiths, 1985) and serotonergic innervation (Vertes, 1991; Waselus et al., 2006). In contrast, our data are consistent with certain previous studies of the dorsal striatum (Wall et al., 2013; Guo et al., 2015), where rabies-mediated transsynaptic tracing led to weak labeling in the VTA and SNc (Wall et al., 2013), and in the DR and MnR (Guo et al., 2015). It was speculated that this weak transsynaptic RV labeling results from the nature of peculiar synapses between dopaminergic neurons and MSNs (Wall et al., 2013). Similar to dopaminergic projections, many serotonergic projections often form varicosities that lack the junctional complexes classically considered to be the morphological substrate for chemical transmission in the central nervous system (Beaudet and Descarries, 1981; Descarries et al., 2010). Our results, together with previous findings, raise the possibility that these specialized synaptic structures may prevent the RV from spreading from the striatal neurons to the presynaptic partners.

Rabies-mediated transsynaptic tracing has been widely used for elucidating neuronal connectivity. As previously published studies have done before us, we must emphasize again the limitations of this technique (Callaway and Luo, 2015), including the fact that it can label only a fraction of inputs to starter cells and that the number of labeled input neurons does not necessarily reflect functional connectivity strength. For example, one input neuron may form synapses with several MSNs, and one starter cell may also be connected to multiple input neurons. It has been established that the spread of RV may be exclusively transsynaptic (restricted to synapses) (Ugolini, 1995; Callaway, 2008; Ugolini, 2008). However, it is not entirely clear how RV crosses synapses and whether transsynaptic spread of RV is biased toward certain connections. Additionally, although we revealed the input neural circuitries of different NAc subregions, it should be noted that the quantification of inputs depends largely on the location of the virus injection and the spread of the virus. It is difficult to cover the entire NAc subregion while leaving adjacent areas uninfected by the viral injection. To lower the possibility of nonspecific infection, we injected a small volume of virus and used a very slow injection rate to limit the spread of the virus to a small range (see Figure 1C). It is possible that other regions that were not targeted in our experiments could show different input patterns.

Another potential caveat arises when we consider that some MSNs coexpress D1 and D2 receptors (D1/D2-MSNs) (Bertran-Gonzalez et al., 2008), posing the possibility that the results of our tracing study using D1R- and D2R-Cre mouse lines contain an overlapping population of starter cells. These D1/D2-MSNs were estimated to account for only 5% of MSNs in the NAc core but up to 17% of MSNs in the NAc shell of mice (Bertran-Gonzalez et al., 2008; Perreault et al., 2011; Gangarossa et al., 2013; Gagnon et al., 2017). Recent studies have shown that this phenomenon of the coexpression of D1 and D2 receptors also exists in rats and nonhuman primates (Perreault et al., 2016; Rico et al., 2017). Furthermore, cholinergic interneurons (CINs), accounting for <5% of the total neuronal population, coexpress D2 receptors in both the striatum and NAc (Le Moine et al., 1990; Pisani et al., 2000; Gagnon et al., 2017), posing the possibility that CINs actually account for a small fraction of the labeled D2-MSNs in our tracing study. It would be intriguing to determine whether the patterns of input organization of these D1/D2-MSNs and CINs are similar to those of D1R-MSNs and D2R-MSNs in the NAc. According to these caveats and limitations, we suggest that our study of inputs to D1R-MSNs and D2R-MSNs in the NAc has revealed trends that likely underestimate the true specificity within NAc circuits. Future studies combining new genetic and viral approaches are necessary to target diverse cell subtypes in the NAc with higher specificity.

In summary, we mapped the organization of inputs to D1R-MSNs and D2R-MSNs in different NAc subregions (core and shell) using transsynaptic retrograde tracing with a modified RV. Although the input distributions of different cell types were highly similar, different NAc subregions had significantly different biased inputs from numerous brain areas. The similarities and differences observed in our study may provide new insight into the diverse functions of the NAc.

Data Availability

All data generated or analyzed in this study are included in the manuscript.

Author Contributions

ZL and TX conceived and designed the project. ZL performed the tracing experiments. ZL, GF, AL, and JY performed the whole-brain data acquisition and analyzed the data. ZL and ZC generated the figures. ZL and TX wrote the paper.

Funding

This work was supported by NSFC Projects Nos. 61721092, 81671374, and 91749209, the director fund of the WNLO, and the Independent Innovation Foundation of HUST (Grant No. 2018KFYXMPT018).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We appreciate the MOST group members of the Britton Chance Center for Biomedical Photonics for assistance with the experiments.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnana.2018.00084/full#supplementary-material

FIGURE S1 | Confocal images showing almost all eGFP- and dsRed-double-positive neurons (starter cells) were also DARPP-32-positive in D1R-Cre mice (top, arrowheads) and D2R-Cre mice (bottom, arrowheads). Scale bar, 10 μm.

TABLE S1 | Raw data used for quantification.

References

Albin, R. L., Young, A. B., and Penney, J. B. (1989). The functional anatomy of basal ganglia disorders. Trends Neurosci. 12, 366–375. doi: 10.1016/0166-2236(89)90074-X

CrossRef Full Text | Google Scholar

Ambroggi, F., Ishikawa, A., Fields, H. L., and Nicola, S. M. (2008). Basolateral amygdala neurons facilitate reward-seeking behavior by exciting nucleus accumbens neurons. Neuron 59, 648–661. doi: 10.1016/j.neuron.2008.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Beaudet, A., and Descarries, L. (1981). The fine structure of central serotonin neurons. J. Physiol. (Paris) 77, 193–203.

Google Scholar

Beckstead, R. M. (1979). An autoradiographic examination of corticocortical and subcortical projections of the mediodorsal-projection (prefrontal) cortex in the rat. J. Comp. Neurol. 184, 43–62. doi: 10.1002/cne.901840104

PubMed Abstract | CrossRef Full Text | Google Scholar

Beier, K. T., Steinberg, E. E., DeLoach, K. E., Xie, S., Miyamichi, K., Schwarz, L., et al. (2015). Circuit architecture of VTA dopamine neurons revealed by systematic input-output mapping. Cell 162, 622–634. doi: 10.1016/j.cell.2015.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Beny-Shefer, Y., Zilkha, N., Lavi-Avnon, Y., Bezalel, N., Rogachev, I., Brandis, A., et al. (2017). Nucleus accumbens dopamine signaling regulates sexual preference for females in male mice. Cell Rep. 21, 3079–3088. doi: 10.1016/j.celrep.2017.11.062

PubMed Abstract | CrossRef Full Text | Google Scholar

Berendse, H. W., and Groenewegen, H. J. (1990). Organization of the thalamostriatal projections in the rat, with special emphasis on the ventral striatum. J. Comp. Neurol. 299, 187–228. doi: 10.1002/cne.902990206

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertran-Gonzalez, J., Bosch, C., Maroteaux, M., Matamales, M., Herve, D., Valjent, E., et al. (2008). Opposing patterns of signaling activation in dopamine D1 and D2 receptor-expressing striatal neurons in response to cocaine and haloperidol. J. Neurosci. 28, 5671–5685. doi: 10.1523/jneurosci.1039-08.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Brog, J. S., Salyapongse, A., Deutch, A. Y., and Zahm, D. S. (1993). The patterns of afferent innervation of the core and shell in the “accumbens” part of the rat ventral striatum: immunohistochemical detection of retrogradely transported fluoro-gold. J. Comp. Neurol. 338, 255–278. doi: 10.1002/cne.903380209

PubMed Abstract | CrossRef Full Text | Google Scholar

Callaway, E. M. (2008). Transneuronal circuit tracing with neurotropic viruses. Curr. Opin. Neurobiol. 18, 617–623. doi: 10.1016/j.conb.2009.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Callaway, E. M., and Luo, L. (2015). monosynaptic circuit tracing with glycoprotein-deleted rabies viruses. J. Neurosci. 35, 8979–8985. doi: 10.1523/jneurosci.0409-15.2015

CrossRef Full Text | Google Scholar

Carlezon, W. A. Jr., and Thomas, M. J. (2009). Biological substrates of reward and aversion: a nucleus accumbens activity hypothesis. Neuropharmacology 56(Suppl. 1), 122–132. doi: 10.1016/j.neuropharm.2008.06.075

PubMed Abstract | CrossRef Full Text | Google Scholar

Chuhma, N., Tanaka, K. F., Hen, R., and Rayport, S. (2011). Functional connectome of the striatal medium spiny neuron. J. Neurosci. 31, 1183–1192. doi: 10.1523/jneurosci.3833-10.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Day, J. J., Jones, J. L., and Carelli, R. M. (2011). Nucleus accumbens neurons encode predicted and ongoing reward costs in rats. Eur. J. Neurosci. 33, 308–321. doi: 10.1111/j.1460-9568.2010.07531.x

PubMed Abstract | CrossRef Full Text | Google Scholar

DeLong, M. R. (1990). Primate models of movement disorders of basal ganglia origin. Trends Neurosci. 13, 281–285. doi: 10.1016/0166-2236(90)90110-V

CrossRef Full Text | Google Scholar

Descarries, L., Riad, M., and Parent, M. (2010). Ultrastructure of the serotonin innervation in the mammalian central nervous system. Handb. Behav. Neurosci. 21, 65–101. doi: 10.1016/S1569-7339(10)70072-2

CrossRef Full Text | Google Scholar

Di Ciano, P., Cardinal, R. N., Cowell, R. A., Little, S. J., and Everitt, B. J. (2001). Differential involvement of NMDA, AMPA/kainate, and dopamine receptors in the nucleus accumbens core in the acquisition and performance of pavlovian approach behavior. J. Neurosci. 21, 9471–9477. doi: 10.1523/JNEUROSCI.21-23-09471.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Do, J. P., Xu, M., and Lee, S. H. (2016). Cell type-specific long-range connections of basal forebrain circuit. eLife 5:e13214. doi: 10.7554/eLife.13214

PubMed Abstract | CrossRef Full Text | Google Scholar

Everitt, B. J. (1990). Sexual motivation: a neural and behavioural analysis of the mechanisms underlying appetitive and copulatory responses of male rats. Neurosci. Biobehav. Rev. 14, 217–232. doi: 10.1016/S0149-7634(05)80222-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Fallon, J. H., and Moore, R. Y. (1978). Catecholamine innervation of the basal forebrain IV. Topography of the dopamine projection to the basal forebrain and neostriatum. J. Compar. Neurol. 180, 545–579. doi: 10.1002/cne.901800310

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrier, I., Cross, A., Johnson, J., Roberts, G., Crow, T., Corsellis, J., et al. (1983). Neuropeptides in Alzheimer type dementia. J. Neurol. Sci. 62, 159–170. doi: 10.1016/0022-510X(83)90196-X

CrossRef Full Text | Google Scholar

Floresco, S. B. (2015). The nucleus accumbens: an interface between cognition, emotion, and action. Annu. Rev. Psychol. 66, 25–52. doi: 10.1146/annurev-psych-010213-115159

PubMed Abstract | CrossRef Full Text | Google Scholar

Floresco, S. B., Ghods-Sharifi, S., Vexelman, C., and Magyar, O. (2006). Dissociable roles for the nucleus accumbens core and shell in regulating set shifting. J. Neurosci. 26, 2449–2457. doi: 10.1523/jneurosci.4431-05.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Floresco, S. B., Seamans, J. K., and Phillips, A. G. (1997). Selective roles for hippocampal, prefrontal cortical, and ventral striatal circuits in radial-arm maze tasks with or without a delay. J. Neurosci. 17, 1880–1890. doi: 10.1523/JNEUROSCI.17-05-01880.1997

CrossRef Full Text | Google Scholar

Francis, T. C., and Lobo, M. K. (2017). Emerging role for nucleus accumbens medium spiny neuron subtypes in depression. Biol. Psychiatry 81, 645–653. doi: 10.1016/j.biopsych.2016.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, K., Miyamoto, Y., Sumi, K., Saika, E., Muramatsu, S.-I., Uno, K., et al. (2017). Overexpression of transmembrane protein 168 in the mouse nucleus accumbens induces anxiety and sensorimotor gating deficit. PLoS One 12:e0189006. doi: 10.1371/journal.pone.0189006

PubMed Abstract | CrossRef Full Text | Google Scholar

Gagnon, D., Petryszyn, S., Sanchez, M. G., Bories, C., Beaulieu, J. M., De Koninck, Y., et al. (2017). Striatal neurons expressing D1 and D2 receptors are morphologically distinct and differently affected by dopamine denervation in mice. Sci. Rep. 7:41432. doi: 10.1038/srep41432

PubMed Abstract | CrossRef Full Text | Google Scholar

Gal, G., Schiller, D., and Weiner, I. (2005). Latent inhibition is disrupted by nucleus accumbens shell lesion but is abnormally persistent following entire nucleus accumbens lesion: the neural site controlling the expression and disruption of the stimulus preexposure effect. Behav. Brain Res. 162, 246–255. doi: 10.1016/j.bbr.2005.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Gangarossa, G., Espallergues, J., de Kerchove, d’Exaerde, A., El Mestikawy, S., Gerfen, C. R., et al. (2013). Distribution and compartmental organization of GABAergic medium-sized spiny neurons in the mouse nucleus accumbens. Front. Neural Circits 7:22. doi: 10.3389/fncir.2013.00022

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerfen, C. R., Engber, T. M., Mahan, L. C., Susel, Z., Chase, T. N., Monsma, F. J., et al. (1990). D1 and D2 dopamine receptor-regulated gene expression of striatonigral and striatopallidal neurons. Science 250, 1429–1432. doi: 10.1126/science.2147780

CrossRef Full Text | Google Scholar

Gong, S., Doughty, M., Harbaugh, C. R., Cummins, A., Hatten, M. E., Heintz, N., et al. (2007). Targeting Cre recombinase to specific neuron populations with bacterial artificial chromosome constructs. J. Neurosci. 27, 9817–9823. doi: 10.1523/jneurosci.2707-07.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorelova, N., and Yang, C. R. (1997). The course of neural projection from the prefrontal cortex to the nucleus accumbens in the rat. Neuroscience 76, 689–706. doi: 10.1016/S0306-4522(96)00380-6

CrossRef Full Text | Google Scholar

Grace, A. A., Floresco, S. B., Goto, Y., and Lodge, D. J. (2007). Regulation of firing of dopaminergic neurons and control of goal-directed behaviors. Trends Neurosci. 30, 220–227. doi: 10.1016/j.tins.2007.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Groenewegen, H. J., Wright, C. I., Beijer, A. V., and Voorn, P. (1999). Convergence and segregation of ventral striatal inputs and outputs. Ann. N. Y. Acad. Sci. 877, 49–63. doi: 10.1111/j.1749-6632.1999.tb09260.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Q., Wang, D., He, X., Feng, Q., Lin, R., Xu, F., et al. (2015). Whole-brain mapping of inputs to projection neurons and cholinergic interneurons in the dorsal striatum. PLoS One 10:e0123381. doi: 10.1371/journal.pone.0123381

PubMed Abstract | CrossRef Full Text | Google Scholar

Hauber, W., and Sommer, S. (2009). Prefrontostriatal circuitry regulates effort-related decision making. Cereb. Cortex 19, 2240–2247. doi: 10.1093/cercor/bhn241

PubMed Abstract | CrossRef Full Text | Google Scholar

Heiman, M., Schaefer, A., Gong, S., Peterson, J. D., Day, M., Ramsey, K. E., et al. (2008). A translational profiling approach for the molecular characterization of CNS cell types. Cell 135, 738–748. doi: 10.1016/j.cell.2008.10.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Hikida, T., Kimura, K., Wada, N., Funabiki, K., and Nakanishi, S. (2010). Distinct roles of synaptic transmission in direct and indirect striatal pathways to reward and aversive behavior. Neuron 66, 896–907. doi: 10.1016/j.neuron.2010.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Hikosaka, O. (2007). Basal ganglia mechanisms of reward-oriented eye movement. Ann. N. Y. Acad. Sci. 1104, 229–249. doi: 10.1196/annals.1390.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, R., Robbins, T. W., Pennartz, C. M., and Everitt, B. J. (2008). Functional interaction between the hippocampus and nucleus accumbens shell is necessary for the acquisition of appetitive spatial context conditioning. J. Neurosci. 28, 6950–6959. doi: 10.1523/jneurosci.1615-08.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, T., Long, B., Gong, H., Xu, T., Li, X., Duan, Z., et al. (2017). A platform for efficient identification of molecular phenotypes of brain-wide neural circuits. Sci. Rep. 7:13891. doi: 10.1038/s41598-017-14360-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawaguchi, Y. (1997). Neostriatal cell subtypes and their functional roles. Neurosci. Res. 27, 1–8. doi: 10.1016/S0168-0102(96)01134-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawaguchi, Y., Wilson, C. J., and Emson, P. C. (1989). Intracellular recording of identified neostriatal patch and matrix spiny cells in a slice preparation preserving cortical inputs. J. Neurophysiol. 62, 1052–1068. doi: 10.1152/jn.1989.62.5.1052

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelley, A. E., Baldo, B. A., Pratt, W. E., and Will, M. J. (2005). Corticostriatal-hypothalamic circuitry and food motivation: integration of energy, action and reward. Physiol. Behav. 86, 773–795. doi: 10.1016/j.physbeh.2005.08.066

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelley, A. E., and Domesick, V. B. (1982). The distribution of the projection from the hippocampal formation to the nucleus accumbens in the rat: an anterograde- and retrograde-horseradish peroxidase study. Neuroscience 7, 2321–2335. doi: 10.1016/0306-4522(82)90198-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kenny, P. J. (2011). Reward mechanisms in obesity: new insights and future directions. Neuron 69, 664–679. doi: 10.1016/j.neuron.2011.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Krause, M., German, P. W., Taha, S. A., and Fields, H. L. (2010). A pause in nucleus accumbens neuron firing is required to initiate and maintain feeding. J. Neurosci. 30, 4746–4756. doi: 10.1523/jneurosci.0197-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kupchik, Y. M., Brown, R. M., Heinsbroek, J. A., Lobo, M. K., Schwartz, D. J., and Kalivas, P. W. (2015). Coding the direct/indirect pathways by D1 and D2 receptors is not valid for accumbens projections. Nat. Neurosci. 18, 1230–1232. doi: 10.1038/nn.4068

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Moine, C., Tison, F., and Bloch, B. (1990). D2 dopamine receptor gene expression by cholinergic neurons in the rat striatum. Neurosci. Lett. 117, 248–252. doi: 10.1016/0304-3940(90)90671-U

CrossRef Full Text | Google Scholar

Lobo, M. K., Covington, H. E. III, Chaudhury, D., Friedman, A. K., Sun, H., Damez-Werno, D., et al. (2010). Cell type-specific loss of BDNF signaling mimics optogenetic control of cocaine reward. Science 330, 385–390. doi: 10.1126/science.1188472

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopes, da Silva, F. H., Arnolds, D. E., and Neijt, H. C. (1984). A functional link between the limbic cortex and ventral striatum: physiology of the subiculum accumbens pathway. Exp. Brain Res. 55, 205–214. doi: 10.1007/BF00237271

PubMed Abstract | CrossRef Full Text | Google Scholar

Lüscher, C. (2016). The emergence of a circuit model for addiction. Annu. Rev. Neurosci. 39, 257–276. doi: 10.1146/annurev-neuro-070815-013920

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeno, H. (1982). Dopamine receptors in canine caudate nucleus. Mol. Cell. Biochem. 43, 65–80. doi: 10.1007/BF00423094

CrossRef Full Text | Google Scholar

Matamales, M., Bertran-Gonzalez, J., Salomon, L., Degos, B., Deniau, J. M., Valjent, E., et al. (2009). Striatal medium-sized spiny neurons: identification by nuclear staining and study of neuronal subpopulations in BAC transgenic mice. PLoS One 4:e4770. doi: 10.1371/journal.pone.0004770

PubMed Abstract | CrossRef Full Text | Google Scholar

Montaron, M. F., Deniau, J. M., Menetrey, A., Glowinski, J., and Thierry, A. M. (1996). Prefrontal cortex inputs of the nucleus accumbens-nigro-thalamic circuit. Neuroscience 71, 371–382. doi: 10.1016/0306-4522(95)00455-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Nie, X., Sun, Y., Wan, S., Zhao, H., Liu, R., Li, X., et al. (2017). Subregional structural alterations in hippocampus and nucleus accumbens correlate with the clinical impairment in patients with Alzheimer’s disease clinical spectrum: parallel combining volume and vertex-based approach. Front. Neurol. 8:399. doi: 10.3389/fneur.2017.00399

CrossRef Full Text | Google Scholar

O’Connor, E. C., Kremer, Y., Lefort, S., Harada, M., Pascoli, V., Rohner, C., et al. (2015). Accumbal D1R neurons projecting to lateral hypothalamus authorize feeding. Neuron 88, 553–564. doi: 10.1016/j.neuron.2015.09.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogawa, S. K., Cohen, J. Y., Hwang, D., Uchida, N., and Watabe-Uchida, M. (2014). Organization of monosynaptic inputs to the serotonin and dopamine neuromodulatory systems. Cell Rep. 8, 1105–1118. doi: 10.1016/j.celrep.2014.06.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Osakada, F., Mori, T., Cetin, A. H., Marshel, J. H., Virgen, B., and Callaway, E. M. (2011). New rabies virus variants for monitoring and manipulating activity and gene expression in defined neural circuits. Neuron 71, 617–631. doi: 10.1016/j.neuron.2011.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Ouimet, C. C., Miller, P. E., Hemmings, H. C. Jr., Walaas, S. I., and Greengard, P. (1984). DARPP-32, a dopamine- and adenosine 3’:5’-monophosphate-regulated phosphoprotein enriched in dopamine-innervated brain regions. III. Immunocytochemical localization. J. Neurosci. 4, 111–124. doi: 10.1523/JNEUROSCI.04-01-00111.1984

CrossRef Full Text | Google Scholar

Parkinson, J. A., Willoughby, P. J., Robbins, T. W., and Everitt, B. J. (2000). Disconnection of the anterior cingulate cortex and nucleus accumbens core impairs Pavlovian approach behavior: further evidence for limbic cortical-ventral striatopallidal systems. Behav. Neurosci. 114, 42–63. doi: 10.1037/0735-7044.114.1.42

PubMed Abstract | CrossRef Full Text | Google Scholar

Paxinos, G., and Franklin, K. J. (2001). The Mouse Brain in Stereotaxic Coordinates. Amsterdam: Elsevier Science.

Google Scholar

Perreault, M. L., Hasbi, A., O’Dowd, B. F., and George, S. R. (2011). The dopamine d1-d2 receptor heteromer in striatal medium spiny neurons: evidence for a third distinct neuronal pathway in Basal Ganglia. Front. Neuroanat. 5:31. doi: 10.3389/fnana.2011.00031

PubMed Abstract | CrossRef Full Text | Google Scholar

Perreault, M. L., Hasbi, A., Shen, M. Y. F., Fan, T., Navarro, G., Fletcher, P. J., et al. (2016). Disruption of a dopamine receptor complex amplifies the actions of cocaine. Eur. Neuropsychopharmacol. 26, 1366–1377. doi: 10.1016/j.euroneuro.2016.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillipson, O., and Griffiths, A. (1985). The topographic order of inputs to nucleus accumbens in the rat. Neuroscience 16, 275–296. doi: 10.1016/0306-4522(85)90002-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Pisani, A., Bonsi, P., Centonze, D., Calabresi, P., and Bernardi, G. (2000). Activation of D2-like dopamine receptors reduces synaptic inputs to striatal cholinergic interneurons. J. Neurosci. 20:Rc69. doi: 10.1523/JNEUROSCI.20-07-j0003.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Pollak Dorocic, I., Furth, D., Xuan, Y., Johansson, Y., Pozzi, L., Silberberg, G., et al. (2014). A whole-brain atlas of inputs to serotonergic neurons of the dorsal and median raphe nuclei. Neuron 83, 663–678. doi: 10.1016/j.neuron.2014.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Richfield, E. K., Penney, J. B., and Young, A. B. (1989). Anatomical and affinity state comparisons between dopamine D1 and D2 receptors in the rat central nervous system. Neuroscience 30, 767–777. doi: 10.1016/0306-4522(89)90168-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Rico, A. J., Dopeso-Reyes, I. G., Martinez-Pinilla, E., Sucunza, D., Pignataro, D., Roda, E., et al. (2017). Neurochemical evidence supporting dopamine D1-D2 receptor heteromers in the striatum of the long-tailed macaque: changes following dopaminergic manipulation. Brain Struct. Funct. 222, 1767–1784. doi: 10.1007/s00429-016-1306-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, S. J., Dietz, D. M., Dumitriu, D., Morrison, J. H., Malenka, R. C., and Nestler, E. J. (2010). The addicted synapse: mechanisms of synaptic and structural plasticity in nucleus accumbens. Trends Neurosci. 33, 267–276. doi: 10.1016/j.tins.2010.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, S. J., and Nestler, E. J. (2013). The brain reward circuitry in mood disorders. Nat. Rev. Neurosci. 14, 609–625. doi: 10.1038/nrn3381

PubMed Abstract | CrossRef Full Text | Google Scholar

Salamone, J. D., and Correa, M. (2012). The mysterious motivational functions of mesolimbic dopamine. Neuron 76, 470–485. doi: 10.1016/j.neuron.2012.10.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Sallee, C. J., and Russell, D. F. (1993). Embedding of neural tissue in agarose or glyoxyl agarose for vibratome sectioning. Biotech. Histochem. 68, 360–368. doi: 10.3109/10520299309105642

CrossRef Full Text | Google Scholar

Saper, C. B. (2006). Staying awake for dinner: hypothalamic integration of sleep, feeding, and circadian rhythms. Prog. Brain Res. 153, 243–252. doi: 10.1016/s0079-6123(06)53014-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sul, J. H., Jo, S., Lee, D., and Jung, M. W. (2011). Role of rodent secondary motor cortex in value-based action selection. Nat. Neurosci. 14, 1202–1208. doi: 10.1038/nn.2881

PubMed Abstract | CrossRef Full Text | Google Scholar

Sul, J. H., Kim, H., Huh, N., Lee, D., and Jung, M. W. (2010). Distinct roles of rodent orbitofrontal and medial prefrontal cortex in decision making. Neuron 66, 449–460. doi: 10.1016/j.neuron.2010.03.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Taverna, S., Ilijic, E., and Surmeier, D. J. (2008). Recurrent collateral connections of striatal medium spiny neurons are disrupted in models of Parkinson’s disease. J. Neurosci. 28, 5504–5512. doi: 10.1523/jneurosci.5493-07.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ugolini, G. (1995). Specificity of rabies virus as a transneuronal tracer of motor networks: transfer from hypoglossal motoneurons to connected second-order and higher order central nervous system cell groups. J. Comp. Neurol. 356, 457–480. doi: 10.1002/cne.903560312

PubMed Abstract | CrossRef Full Text | Google Scholar

Ugolini, G. (2008). Use of rabies virus as a transneuronal tracer of neuronal connections: implications for the understanding of rabies pathogenesis. Dev. Biol. (Basel) 131, 493–506.

PubMed Abstract | Google Scholar

Urstadt, K. R., Kally, P., Zaidi, S. F., and Stanley, B. G. (2013). Ipsilateral feeding-specific circuits between the nucleus accumbens shell and the lateral hypothalamus: regulation by glutamate and GABA receptor subtypes. Neuropharmacology 67, 176–182. doi: 10.1016/j.neuropharm.2012.10.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Vertes, R. P. (1991). A PHA-L analysis of ascending projections of the dorsal raphe nucleus in the rat. J. Comp. Neurol. 313, 643–668. doi: 10.1002/cne.903130409

PubMed Abstract | CrossRef Full Text | Google Scholar

Wall, N. R., De La Parra, M., Callaway, E. M., and Kreitzer, A. C. (2013). Differential innervation of direct- and indirect-pathway striatal projection neurons. Neuron 79, 347–360. doi: 10.1016/j.neuron.2013.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., He, X., Zhao, Z., Feng, Q., Lin, R., Sun, Y., et al. (2015). Whole-brain mapping of the direct inputs and axonal projections of POMC and AgRP neurons. Front. Neuroanat. 9:40. doi: 10.3389/fnana.2015.00040

PubMed Abstract | CrossRef Full Text | Google Scholar

Waselus, M., Galvez, J. P., Valentino, R. J., and Van Bockstaele, E. J. (2006). Differential projections of dorsal raphe nucleus neurons to the lateral septum and striatum. J. Chem. Neuroanat. 31, 233–242. doi: 10.1016/j.jchemneu.2006.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Watabe-Uchida, M., Zhu, L., Ogawa, S. K., Vamanrao, A., and Uchida, N. (2012). Whole-brain mapping of direct inputs to midbrain dopamine neurons. Neuron 74, 858–873. doi: 10.1016/j.neuron.2012.03.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Weissbourd, B., Ren, J., DeLoach, K. E., Guenthner, C. J., Miyamichi, K., and Luo, L. (2014). Presynaptic partners of dorsal raphe serotonergic and GABAergic neurons. Neuron 83, 645–662. doi: 10.1016/j.neuron.2014.06.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickersham, I. R., Lyon, D. C., Barnard, R. J., Mori, T., Finke, S., Conzelmann, K. K., et al. (2007). Monosynaptic restriction of transsynaptic tracing from single, genetically targeted neurons. Neuron 53, 639–647. doi: 10.1016/j.neuron.2007.01.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, C. I., Beijer, A. V., and Groenewegen, H. J. (1996). Basal amygdaloid complex afferents to the rat nucleus accumbens are compartmentally organized. J. Neurosci. 16, 1877–1893. doi: 10.1523/JNEUROSCI.16-05-01877.1996

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C. R., and Mogenson, G. J. (1984). Electrophysiological responses of neurones in the nucleus accumbens to hippocampal stimulation and the attenuation of the excitatory responses by the mesolimbic dopaminergic system. Brain Res. 324, 69–84. doi: 10.1016/0006-8993(84)90623-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Yawata, S., Yamaguchi, T., Danjo, T., Hikida, T., and Nakanishi, S. (2012). Pathway-specific control of reward learning and its flexibility via selective dopamine receptors in the nucleus accumbens. Proc. Natl. Acad. Sci. U.S.A. 109, 12764–12769. doi: 10.1073/pnas.1210797109

PubMed Abstract | CrossRef Full Text | Google Scholar

Zahm, D. S., and Brog, J. S. (1992). On the significance of subterritories in the “accumbens” part of the rat ventral striatum. Neuroscience 50, 751–767. doi: 10.1016/0306-4522(92)90202-D

CrossRef Full Text | Google Scholar

Zhang, Z., Zhang, H., Wen, P., Zhu, X., Wang, L., Liu, Q., et al. (2017). Whole-brain mapping of the inputs and outputs of the medial part of the olfactory tubercle. Front. Neural Circuits 11:52. doi: 10.3389/fncir.2017.00052

PubMed Abstract | CrossRef Full Text | Google Scholar