Materials and Methods
Animals
Male Sprague Dawley rats (initial weight 250–300 g; Charles River, Raleigh, NC, USA) were single-housed under a 12-h reverse light/dark cycle (ZT0 = 19:00) and had access to food and water ad libitum. Animals were housed in a temperature- and humidity-controlled animal facility with AAALAC accreditation. All experiments were approved by the Institutional Animal Care and Use Committee at Texas A&M University and conducted according to specifications of the National Institutes of Health as outlined in the Guide for the Care and Use of Laboratory Animals.
Viruses
The transduction and monosynaptic spread of glycoprotein (G)-deleted rabies virus (RV) with EnvA pseudotyping (EnvA-RV-EGFP) is limited to neurons expressing TVA receptor and rabies glycoprotein (RG) in a Cre-dependent manner (Wickersham et al., 2007a, b; Watabe-Uchida et al., 2012). We targeted specific PL neurons by injecting either NAc core or RMTg with retrograde Cre-expressing virus (CAV2-Cre) and injecting PL with Cre-dependent helper AAV viruses (TVA and RG) and EnvA-RV-EGFP.
EnvA-RV-EGFP (titer of 4.3 × 10e8) was provided by Ed Callaway via the Gene Transfer, Targeting, and Therapeutics Facility at the Salk Institute (Wickersham et al., 2007a). AAV-TVA (AAV1-EF1a-FLEX-TVAmCherry, a titer of 4 × 10e12) and AAV-RG (AAV1-CA-FLEX-RG, a titer of 4 × 10e12) were obtained through the viral vector core at the University of North Carolina (developed by Watabe-Uchida et al., 2012). CAV2-Cre (CAV2-CMV-Cre recombinase) was provided by Eric J. Kremer at the Institut de Genetique Moleculaire de Montpellier (Kremer et al., 2000). We diluted CAV2-Cre 1:10 in 10% glycerol in PBS for a final titer of 10.3 × 10e11 ppml, based on preliminary experiments indicating cell toxicity at higher titers. All viral procedures were approved by the Institutional Biosafety Committee at Texas A&M University and conducted according to specifications of the National Institutes of Health.
Stereotaxic Surgery
Rats were anesthetized with vaporized isoflurane (induced at 5%, maintained at 1–2%), given the analgesic ketoprofen (2 mg/kg, s.c.), and placed in a stereotaxic frame (Kopf, Tujunga, CA, USA). A skin incision was made over the skull and holes were drilled in the skull over the target sites. Intracranial injections were made using a pulled glass micropipette and Nanoject 2010 injector (World Precision Instruments, Sarasota, FL, USA). After the injection, the pipette was slowly raised from the brain and the skin was stapled closed.
For conventional retrograde tracing, rats were given unilateral microinjections of 100 nl of 2% Fluoro-Gold (FG, Fluorochrome, Denver, CO, USA) into unilateral NAc core (AP +1.8, ML +2.6, DV –7.3 from bregma, 6°) and 100 nl of 0.2% cholera toxin subunit b (CTB, List Biological Laboratories, Campbell, CA, USA) into RMTg (AP –7.2, ML +1.9, DV –7.6 from dura, 10°). FG was used only in NAc because it tended to cause toxicity in RMTg, evident by circling of animals after unilateral injection. Rats were sacrificed at least 1 week after injections.
For monosynaptic retrograde tracing, rats were given two surgeries. During the first stereotaxic surgery, CAV2-Cre (1,000 nl) plus biotin dextran (10 K MW, 0.2% final volume, for visualization of injection site) was injected unilaterally into either NAc core (AP +1.8, ML +2.6, DV –7.3 from bregma, 6°) or RMTg (AP –7.2, ML +1.8, DV –7.6 from dura, 10°), and a mixture of AAV-RG and AAV-TVA (1,000 nl of a 1:1 mix) was injected into ipsilateral PL (AP +3.1, ML +1.2, DV –3.8 from skull, 12°). Three weeks later during a second stereotaxic surgery, EnvA-RV-EGFP (1,000 nl) was injected into the same location in ipsilateral PL. Rats were sacrificed 7 days after rabies injections.
Tissue Processing
Rats were deeply anesthetized with isoflurane and then transcardially perfused with 0.9% NaCl followed by 10% neutral-buffered formalin via a peristaltic pump (Cole Parmer, Vernon Hills, IL, USA). Brains were collected, post-fixed in 10% formalin overnight at 4°C, placed in 20% sucrose in phosphate-buffered saline (PBS) with 0.02% sodium-azide at 4°C for ≥2 days, frozen in dry ice, cut into 40-μm thick coronal sections on a cryostat, and collected into PBS-azide for storage before immunohistochemistry (IHC). For IHC, all incubations and rinses took place at room temperature on a shaker. Sections were rinsed in PBS three times between steps.
For FG and CTB tracing, free-floating sections were incubated overnight in goat anti-CTB (1:50 K, List Biological Labs Cat# 703, RRID: AB_10013220) in PBS with 0.25% Triton X-100 (PBST, Sigma–Aldrich), and then incubated for 1 h in Alexa Fluor 594-conjugated donkey anti-goat (1:500 in PBST, Jackson Immunoresearch Labs cat# 705-585-147, RRID:AB_2340433). Sections were mounted onto glass slides, coverslipped with ProLong Diamond Antifade (Thermo Fisher Scientific, Fair Lawn, NJ, USA), sealed with clear nail polish, and stored at 4°C.
For RV tracing, one well (every 12th section) was used for mCherry/EGFP fluorescence (for counting starter cells), and one well (every 12th section) was used for IHC to visualize mCherry/EGFP with DAB (for counting inputs to PL). For fluorescence, sections were mounted onto glass slides, coverslipped with ProLong Diamond Antifade, sealed with clear nail polish, and stored at 4°C. For DAB, free-floating sections were incubated for 15 min in 0.3% H2O2 in PBS, overnight in mouse anti-DsRed (1:2 K in PBST, Takara Bio Cat# 632392, RRID: AB_2801258), 1 h in biotinylated donkey anti-mouse (1:500 in PBST, Jackson Immunoresearch Labs Cat# 715–065–151, RRID: AB_2340785), and 45 min in ABC (Vector Elite Kit, 1:500 in PBST, Vector). The reaction was visualized via incubation for 10 min in 0.025% 3,3′-diaminobenzidine (DAB), 0.05% nickel ammonium sulfate, and 0.015% H2O2 (to yield black DAB color). Sections were then incubated overnight in rabbit anti-GFP (1:50 K in PBST, Abcam Cat# ab290, RRID:AB_303395), 1 h in biotinylated donkey anti-rabbit (1:500 in PBST, Jackson Immunoresearch Labs Cat# 711-065-152, RRID: AB_2340593), 45 min in ABC (Vector Elite Kit, 1:500 in PBST, Vector), and visualized via the same DAB reaction but without nickel (to yield brown DAB color). Sections were mounted onto Superfrost Plus slides, dried, counterstained with Methyl Green (0.5% in sodium acetate buffer, Sigma–Aldrich), and coverslipped with Permount (Thermo Fisher Scientific, Fair Lawn, NJ, USA).
Data Analysis
Fluorescent and brightfield images were acquired using an Olympus BX51 microscope. For CTB and FG tracing, neurons were counted in PL on two sections per rat. Four rats were excluded from CTB/FG analysis due to small/missing RMTg injections. One rat was removed from CTB/FG analysis after being identified as a significant outlier (Grubb's method) due to a higher percentage of double labeling. The remaining data were analyzed with a paired t-test. For rabies tracing, all EGFP-labeled neurons were counted on all sections stained with DAB (every 12th section) and categorized into brain regions according to boundaries identified by Paxinos and Watson (2007). RMTg boundaries were defined (Figures 1B, 2C) according to Smith et al. (2019). A neuron was only counted if a soma was clearly visible. Fluorescent sections were used to count all double-labeled neurons (EGFP + mCherry) on every 12th section. Expression of EGFP alone indicated input cells, and expression of both mCherry and EGFP indicated starter cells. Inputs were normalized within each animal and calculated as a percentage of total inputs per rat. Brain regions were only included in analysis and figures if normalized inputs were ≥0.2% in ≥2 rats. Data were analyzed using two-way ANOVAs (with repeated measures when appropriate) and Sidak post hoc analyses.
Results
Conventional Retrograde Tracing in PL Targets
Sparse double labeling was seen in PL (Figure 1A) when conventional retrograde tracers (FG and CTB) were injected into ipsilateral NAc core and RMTg (Figure 1B; n = 11). We observed NAc-projecting PL neurons in layers II/III and upper layer V, whereas RMTg-projecting neurons were located in deeper layer V (Figure 1A). Double-labeled neurons (projecting to both NAc core and RMTg) represented only 0.5% ± 0.1 of PL neurons projecting to NAc core and 1.1% ± 0.3 of PL neurons projecting to RMTg (Figure 1C; t(10) = 3.004, p = 0.013). These data indicate that PL neurons projecting to NAc core vs. RMTg are predominantly separate subpopulations.
Monosynaptic Retrograde Tracing in PL Subpopulations
To determine whether these PL subpopulations differ in terms of input specificity, we conducted monosynaptic retrograde tracing in each projection subpopulation via retrograde Cre delivery (Figure 2A). CAV2-Cre was injected into either NAc core (n = 4) or RMTg (n = 4; Figures 2B,C), and Cre-dependent helper AAVs driving expression of TVA-mCherry and RG were injected into PL. Three weeks later, EnvA-RV-EGFP was injected into PL. Rabies only infects and monosynaptically spreads from neurons expressing both TVA and RG.
We mapped and counted the number of starter cells (mCherry + EGFP) using fluorescent labeling (representative images in Figures 2D–G). We counted the number of input cells (EGFP) across different brain regions using immunohistochemistry with DAB labeling (representative images in Figure 4). Starter cells projecting to either NAc core or RMTg were primarily in PL, but also in neighboring infralimbic cortex (IL; Figure 3A). Although the average number of starter cells differed for NAc vs. RMTg (t(6) = 5.50, p = 0.0015), the number of transsynaptically labeled cells was correlated with the number of starter cells (Figure 3B; r = 0.91, p = 0.0019).
Inputs were normalized within each animal and calculated as a percentage of total inputs per rat (Figure 3C). The total ipsilateral inputs were significantly greater than total contralateral inputs (main effect for contra/ipsi: F(1,6) = 37,354, p < 0.0001), and the two PL subpopulations differed in terms of ipsilateral vs. contralateral inputs (main effect for PL subpopulation: F(1,6) = 14.24, p = 0.009; interaction: F(1,6) = 83.43, p < 0.0001). Therefore, we separately analyzed contralateral and ipsilateral inputs when comparing afferents for the two PL subpopulations. PL neurons projecting to NAc core vs. RMTg were significantly different in terms of ipsilateral input (two-way ANOVA, overall effect of PL subpopulation: F(1,6) = 53.44, p = 0.0003; overall effect of afferents: F(29,174) = 86.11, p < 0.0001; interaction: F(29,174) = 5.05; p < 0.0001). We found significant differences for ipsilateral PL and IL (p's < 0.0001), and mediodorsal thalamic nucleus (MD; p = 0.036). PL neurons projecting to NAc core vs. RMTg were also significantly different in terms of contralateral input (overall effect of PL subpopulation: F(1,6) = 119.9, p < 0.0001; overall effect of afferents: F(6,36) = 60.46, p < 0.0001; interaction: F(6,36) = 34.06; p < 0.0001), with significant differences for contralateral PL (p < 0.0001).
For both populations, the primary source of afferent information was neighboring prefrontal neurons in PL and IL (Figures 4A,B). However, RMTg-projecting PL neurons received a significantly higher percentage of inputs from local ipsilateral PL and IL neurons (Figure 3C; PL: 39.4% ± 0.2 and IL: 20.1% ± 3.8), as compared to NAc-projecting PL neurons (PL: 26.9% ± 4.0 and IL: 11.8% ± 1.3). We also observed that RMTg-projecting PL neurons received a higher percentage of input from the medial orbitofrontal cortex (OFC), while NAc-projecting PL neurons received a higher percentage of input from the dorsal tenia tecta, lateral OFC, and agranular insular cortex. The percentage of inputs was similar for NAc-projecting and RMTg-projecting PL neurons for secondary motor cortex, anterior olfactory nucleus, ventral OFC, dorsal peduncular cortex, cingulate cortex areas 1 and 2, entorhinal/ectorhinal/perirhinal cortex, and piriform cortex. Interestingly, NAc-projecting PL neurons received a significantly greater percentage of contralateral cortical input from PL (2.92% ± 0.30) as compared to RMTg-projecting PL neurons (0.43% ± 0.15). NAc-projecting PL neurons also received a greater percentage of contralateral inputs from IL, secondary motor cortex, all areas of OFC (medial, ventral, and lateral), and cingulate cortex area 1.
Both PL subpopulations received prominent afferent information from the thalamus (Figures 4C,D), including subregions within the midline thalamic nuclei (anterior/posterior parts of paraventricular, paratenial, rhomboid, reuniens, ventral reuniens), anterior thalamic nuclei (anteromedial, ventrolateral/dorsomedial parts of anteroventral, interanteromedial), reticular thalamic nucleus, MD (mediodorsal, lateral/central/medial parts of mediodorsal, intermediodorsal), lateral thalamic nuclei (dorsomedial part of lateroventral, mediorostral part of lateral posterior), intralaminar thalamic nuclei (paracentral, central medial, centrolateral, and parafascicular), and ventral thalamic nuclei (ventral anterior, ventromedial, ventral posteromedial, ventrolateral; Figure 3C). NAc-projecting PL neurons received a higher percentage of inputs from all thalamic regions, as compared to RMTg-projecting PL neurons, with a significant difference observed for MD (6.6% ± 0.7 vs. 2.1% ± 0.6).