Abstract
Brainstem median raphe (MR) neurons expressing the serotonergic regulator gene Pet1 send collateralized projections to forebrain regions to modulate affective, memory-related, and circadian behaviors. Some Pet1 neurons express a surprisingly incomplete battery of serotonin pathway genes, with somata lacking transcripts for tryptophan hydroxylase 2 (Tph2) encoding the rate-limiting enzyme for serotonin [5-hydroxytryptamine (5-HT)] synthesis, but abundant for vesicular glutamate transporter type 3 (Vglut3) encoding a synaptic vesicle-associated glutamate transporter. Genetic fate maps show these nonclassical, putatively glutamatergic Pet1 neurons in the MR arise embryonically from the same progenitor cell compartment—hindbrain rhombomere 2 (r2)—as serotonergic TPH2+ MR Pet1 neurons. Well established is the distribution of efferents en masse from r2-derived, Pet1-neurons; unknown is the relationship between these efferent targets and the specific constituent source-neuron subgroups identified as r2-Pet1Tph2-high versus r2-Pet1Vglut3-high. Using male and female mice, we found r2-Pet1 axonal boutons segregated anatomically largely by serotonin+ versus VGLUT3+ identity. The former present in the suprachiasmatic nucleus, paraventricular nucleus of the thalamus, and olfactory bulb; the latter are found in the hippocampus, cortex, and septum. Thus r2-Pet1Tph2-high and r2-Pet1Vglut3-high neurons likely regulate distinct brain regions and behaviors. Some r2-Pet1 boutons encased interneuron somata, forming specialized presynaptic “baskets” of VGLUT3+ or VGLUT3+/5-HT+ identity; this suggests that some r2-Pet1Vglut3-high neurons may regulate local networks, perhaps with differential kinetics via glutamate versus serotonin signaling. Fibers from other Pet1 neurons (non-r2-derived) were observed in many of these same baskets, suggesting multifaceted regulation. Collectively, these findings inform brain organization and new circuit nodes for therapeutic considerations.
SIGNIFICANCE STATEMENT Our findings match axonal bouton neurochemical identity with distant cell bodies in the brainstem raphe. The results are significant because they suggest that disparate neuronal subsystems derive from Pet1+ precursor cells of the embryonic progenitor compartment rhombomere 2 (r2). Of these r2-Pet1 neuronal subsystems, one appears largely serotonergic, as expected given expression of the serotonergic regulator PET1, and projects to the olfactory bulb, thalamus, and suprachiasmatic nucleus. Another expresses VGLUT3, suggesting principally glutamate transmission, and projects to the hippocampus, septum, and cortex. Some r2-Pet1 boutons—those that are VGLUT3+ or VGLUT3+/5-HT+ co-positive—comprise “baskets” encasing interneurons, suggesting that they control local networks perhaps with differential kinetics via glutamate versus serotonin signaling. Results inform brain organization and circuit nodes for therapeutic consideration.
Introduction
Neurons of the brainstem median raphe (MR) that express the serotonergic regulator gene Pet1 (alias Fev) project to numerous forebrain sites (Kosofsky and Molliver, 1987; Kiyasova and Gaspar, 2011; Bang et al., 2012; Gaspar and Lillesaar, 2012) to modulate diverse processes including emotional behavior (Teissier et al., 2015; Abela et al., 2020), sensorimotor gating (Okaty et al., 2015), memory (Ohmura et al., 2010; Baskin et al., 2020), and sleep and circadian rhythms (Meyer-Bernstein and Morin, 1996; Glass et al., 2003; Iwasaki et al., 2018). Recent work shows that some of these MR Pet1 neurons, even with their expression of Pet1, lack detectable transcript levels for various serotonin pathway genes (Pelosi et al., 2014; Okaty et al., 2015, 2019; Sos et al., 2017; Ren et al., 2019), but demonstrate high levels of Slc17a8 (alias Vglut3) mRNA, which encodes the type 3 transporter [vesicular glutamate transporter type 3 (VGLUT3)] for loading glutamate into synaptic vesicles. Genetic fate maps coupled with single-cell transcriptomic analyses show that these Vglut3+, tryptophan hydroxylase 2 (Tph2)–, putative glutamatergic MR Pet1 neurons arise embryonically from the same progenitor cell compartment—hindbrain rhombomere 2 (r2)—as do more classically serotonergic Tph2+, serotonin [5-hydroxytryptamine (5-HT)+] MR Pet1 neurons (Okaty et al., 2015). Though the brain regions innervated by the collective r2-derived Pet1 neuron population have been mapped (Bang et al., 2012), the degree to which the nonclassical r2-Pet1Vglut3-high neuron subgroup and the classical r2-Pet1Tph2-high neuron subgroup differentially contribute to these projections is unknown. Here we present progress toward resolving this circuitry.
One possibility is that the r2-Pet1Tph2-high and r2-Pet1Vglut3-high neuron subgroups exhibit converging circuitry, innervating the same target regions and potentially exerting regulatory effects of different directions and/or timescales. An alternative extreme is that they exhibit divergent circuitry, each innervating different brain regions to control different functions. Also possible is a combination of both as a function of target site. Because neurochemical identity (serotonin and TPH2 vs VGLUT3 as a proxy for glutamate) largely distinguishes the r2-Pet1Tph2-high from r2-Pet1Vglut3-high neuron somata, we reasoned it might also distinguish their boutons and thus target brain regions, which in turn could help discern from among these circuitry possibilities.
In this study, we used dual recombinase-based genetic labeling strategies (Awatramani et al., 2003; Bang et al., 2012), retrograde viral tracing, immunohistology, and single-cell in situ mRNA detection combined with confocal microscopy to describe the neurochemical identity and distribution of axonal boutons arising from the r2-Pet1 neuron population. We found anatomic segregation of r2-Pet1 axonal boutons by serotonin+ versus VGLUT3+ identity, suggesting that the r2-Pet1Tph2-high and r2-Pet1Vglut3-high neuron subgroups innervate distinct brain regions. Further, we found that some VGLUT3+ r2-Pet1 boutons comprised a specialized presynaptic structure that has been called a pericellular basket (Köhler et al., 1982; Hornung and Celio, 1992; Dinopoulos et al., 1993; Riedel et al., 2008), here found encasing the soma and dendrites of individual hippocampal, cortical, and septal GABAergic neurons of specific subclasses. This finding suggests that some r2-Pet1Vglut3-high neurons may have a privileged capacity to control local network excitability. In the septum, we observed co-positive VGLUT3+/5-HT+ r2-Pet1 boutons, and they were enriched in these pericellular baskets, suggesting that this may be a specialized site of 5-HT/glutamate cotransmission. Further, we observed that the septal cells targeted by r2-Pet1 baskets commonly also received serotonin+ fibers from other (non-r2) Pet1 neurons, suggesting even additional regulatory complexity. This work provides a foundational framework for pursuit of downstream circuits and cells relevant to specific behavioral modulation by r2-Pet1 neurons, and further extends our knowledge of the functional organization of the median raphe nucleus.
Materials and Methods
Animals
All experimental procedures were approved by the Harvard Institutional Animal Care and Use Committee. Mice were housed in a colony maintained on a 12 h light/dark cycle. All transgenic strains described below have been previously backcrossed to the C57BL/6J inbred strain (stock #000664, The Jackson Laboratory) for at least nine generations. In all experiments, adult mice were used [at least postnatal day 60 (P60)], as previous research has shown that the postnatal maturation of serotonin neuron fiber morphology is finalized before this age (Maddaloni et al., 2017).
For projection mapping studies, r2-Pet1 neuron boutons were visualized by cell-autonomous expression of a synaptophysin-GFP (synGFP) fusion protein using intersectional genetics. These studies used triple-transgenic mice harboring the r2-specific Cre driver Hoxa2-cre (alias r2Hoxa2-cre; Awatramani et al., 2003) and the Flp driver Pet1-Flpe (Jensen et al., 2008) in combination with the dual Cre/Flp-responsive allele Gt(ROSA)26Sortm10(CAG-Syp/EGFP+,-tdTomato)Dym (termed RC-FPSit; stock #030206, The Jackson Laboratory; Niederkofler et al., 2016). Triple-transgenic r2Hoxa2-cre, Pet1-Flpe, RC-FPSit+/− mice were generated by crossing double-transgenic Pet1-Flpe, RC-FPSit−/− males to r2Hoxa2-cre females. To fluorescently label the somata and fibers of r2-Pet1 neurons with EGFP (called the intersectional cell population) and all other Pet1 neurons with mCherry (called the subtractive cell population), we generated triple-transgenic mice bearing the drivers r2Hoxa2-cre and Pet1-Flpe, together with the dual Cre/Flp-responsive ROSA26 allele Gt(ROSA)26Sortm8(CAG-mCherry,-EGFP)Dym (termed RC-FrePe; stock #029486, The Jackson Laboratory; Dymecki et al., 2010; Brust et al., 2014; Fig. 1). Triple-transgenic r2Hoxa2-cre, Pet1-Flpe, RC-FrePe+/− mice were generated by crossing double-transgenic Pet1-Flpe hemizygous, RC-FrePe−/− homozygous males to hemizygous r2Hoxa2-cre females. We also used a similar intersectional line Gt(ROSA)26Sortm1.3(CAG-tdTomato,-EGFP)Pjen (termed RC-FLTG; Plummer et al., 2015) and comparable breeding strategy to visualize r2-Pet1 and other Pet1 subtractive fibers (stock #026932, The Jackson Laboratory). To fluorescently mark Piezo2-Pet1 neurons, we generated triple-transgenic mice bearing the drivers Piezo2-GFP-IRES-cre (Woo et al., 2014) and Pet1-Flpe, combined with the conditional allele Gt(ROSA)26Sortm65.1(CAG-tdTomato)Hze (termed RC-Ai65; stock #021875, The Jackson Laboratory; Madisen et al., 2012, 2015). Triple-transgenic r2Hoxa2-cre, Pet1-Flpe, RC-Ai65+/− mice were generated by crossing double-transgenic Pet1-Flpe hemizygous, RC-Ai65−/− homozygous mice with heterozygous Piezo2-GFP-IRES-cre mice.
Tissue preparation for immunohistochemistry and fluorescence in situ hybridization
All experiments used brain sections from mice perfused transcardially with 0.1 m PBS, followed by 4% paraformaldehyde (PFA) as previously published (Brust et al., 2014). Brains were fixed overnight in 4% PFA/PBS for 24 h at 4°C. For bouton neurochemical phenotyping or RNAscope experiments, tissue was then cryoprotected in 30% sucrose/PBS for 48–72 h at 4°C and then embedded in tissue freezing medium (Triangle Biomedical Services) for storage at –80°C. Tissue was sectioned into free-floating [20 µm for single-molecule fluorescence in situ hybridization (FISH) and 40 µm for immunohistochemistry (IHC)] coronal sections in series of six or eight. For retrograde tracing, sections at 50 µm thickness were cut on a vibratome (model VT1000S, Leica) and stored in 28% sucrose 30% ethylene glycol cryoprotectant at −30°C.
Immunohistochemistry
We performed IHC using the following antibodies and dilutions: anti-calbindin (mouse monoclonal; 1:2000; catalog #C9848, Sigma-Aldrich), anti-calbindin D-28k (mouse monoclonal; 1:3000; catalog #300, Swant), anti-DsRed (rabbit polyclonal; 1:1000; catalog #632496, Takara), anti-GABA (rabbit polyclonal; 1:250; catalog #ABN131, Millipore), anti-GFP (chicken polyclonal; 1:1500; catalog #GFP-1020, Aves Labs), anti-Mouse IgG biotinylated (donkey polyclonal; 1:500; catalog #A16021, Thermo Fisher Scientific), anti-neurokinin B (NKB; rabbit polyclonal; 1:500; catalog #NB300-201SS, Novus Biologicals), anti-red fluorescent protein (rat monoclonal; 1:1000; catalog #5f8, Chromotek), anti-Rabbit IgG biotinylated (donkey polyclonal; 1:200; catalog #711–065-152, Jackson ImmunoResearch), anti-Reelin (goat polyclonal; 1:500; catalog #AF3820, Novus Biologicals), anti-5-HT (goat polyclonal; 1:500; catalog #ab66047, Abcam), anti-5-HT2C (goat polyclonal; 1:100; catalog #sc-15 081, Santa Cruz Biotechnology), anti-5-HT3A (rabbit polyclonal; 1:500; catalog #ASR-031, Alomone Labs), anti-somatostatin (SST; rabbit polyclonal; 1:500; catalog #HPA019472, Sigma-Aldrich), anti-TPH2 (rabbit polyclonal; 1:1000; catalog #NB100-74 555, Novus Biologicals), anti-VGLUT3 (rabbit polyclonal; 1:500; catalog #135–203, Synaptic Systems), anti-VGLUT3 (guinea pig polyclonal; 1:500; Synaptic Systems, catalog #135–204), and anti-VIP (vasoactive intestinal peptide; rabbit polyclonal; 1:500; catalog #20077, Immunostar).
Floating sections were incubated for 1 h at room temperature in a blocking solution of Tris-buffered saline (TBS) containing 3% normal donkey serum (NDS; Jackson ImmunoResearch) and 0.3% Triton X-100 (Sigma-Aldrich). Sections were then transferred to a solution containing primary antibodies diluted in TBS with 1% NDS and 0.3% Triton X-100 and incubated for 24 h at 4°C on a nutator mixer (BD Clay Adams). After three TBS washes, sections were transferred to a solution containing secondary antibodies diluted in TBS with 1% NDS and 0.3% Triton X-100, and incubated in a light-protected chamber for 2 h at room temperature. Detection of 5-HT3A receptor followed previously established protocols (Chang et al., 2016), with ethanol dehydration before blocking. For detection of subtractive Pet1 fibers (non-r2-Pet1 fibers) in composite basket phenotyping experiments and for detection of calbindin-expressing interneurons in the basket target quantification experiments, we used a modified protocol. Briefly, amplification of signal was achieved by modifying the second day protocol using a biotinylated secondary antibody to primary antibody host IgG (see above) for 2 h at room temperature. Following manufacturer recommendations, the calbindin D-28k antibody was dissolved in TBS with 10% NDS. For DsRed, we dissolved the biotinylated antibody in TBS with 2% NDS and 0.1% Triton X-100. To reduce nonspecific signal, we also used an avidin-biotin blocking kit (catalog #SP-2001, Vector Laboratories) on IHC day 1 following manufacturer instructions. In general, for detection of primary antibodies, species-matched Alexa Fluor 488-, Alexa Fluor 568-, Cy5-, and Alexa Fluor 647-conjugated secondary antibodies (Jackson ImmunoResearch) were used. For visualizing biotinylated antibodies, we incubated sections with Alexa Fluor 555-conjugated streptavidin (1:500; catalog #S32355, Thermo Fisher Scientific) during the same step as secondary antibodies. Sections were then washed three times in fresh TBS at room temperature as before, with 1 µg/ml DAPI (Life Technologies) added to the last wash step lasting 15 min. Sections were mounted onto no. 1.5 coverslips (Electron Microscopy Sciences) and then attached to SuperFrost Plus Histologic Slides (Thermo Fisher Scientific) using Aqua-Poly/Mount medium (Polysciences) or ProLong Glass anti-fade mountant (Thermo Fisher Scientific).
Multiplexed fluorescent in situ hybridization
Single-molecule mRNA FISH (smFISH) was performed on fresh frozen or fixed tissue sections according to the manufacturer instructions (RNAscope, Advanced Cell Diagnostics; Wang et al., 2012). Sections were mounted onto SuperFrost Plus slides (Thermo Fisher Scientific) and dried on a slide warmer at 45°C for 30 min to promote section adhesion before proceeding with the RNAscope protocol recommended by suppliers (Advanced Cell Diagnostics). The RNAscope Manual Assay 2.5 was performed using the Multiplex Fluorescence Detection Kit (catalog #320851, ACD) and probes for eGFP (catalog #400281, ACD), Fev (Pet1; catalog #413241-C3, ACD), Tph2 (catalog #Mm-318691-C2, CD), and Vglut3 (catalog #Mm-431261-C1, ACD). For combined smFISH and IHC, following FISH, sections were then immunostained for EGFP as above and counterstained with DAPI. Slides were mounted using ProLong Gold Antifade Mounting media (Thermo Fisher Scientific) and no. 1.5 thickness glass coverslips.
Retrograde tracing experiments
To label r2-Pet1 cells that project to the hippocampus (HP), we injected three adult (two male, one female) r2Hoxa2-cre, Pet1-Flpe, RC-FrePe–/+ mice with AAV2-retro-CAG-tdTomato-WPRE virus (catalog #59462-AAVrg, Addgene) bilaterally into the dorsal hippocampus [coordinates: mediolateral (ML), ±1.4 mm; anteroposterior (AP), −2.06 mm; dorsoventral (DV), −1.75 mm]. The retrograde serotype of this virus allows infection of projection neurons that innervate the region targeted by viral injections and results in virus-induced marking of infected neurons by tdTomato. Three weeks postsurgery, mice were perfused and tissue processed as above. Sections were stained for EGFP and VGLUT3 or 5-HT using a modified immunohistochemistry protocol. Briefly, sections were rinsed 3× in PBS, then 1× in PBS with 0.1% Triton X-100 (PBST; Sigma-Aldrich) before blocking in PBST with 5% NDS for 2 h at room temperature. After incubating with primary antibodies for 48 h at 4°C, sections are rinsed 3× in PBST and then incubated with secondary antibodies in PBST with 2% NDS for 2 h at room temperature. After rinsing 3× in PBST (last rinse for 10 min with 1.5 μg/ml DAPI), sections were mounted on #1.5 thickness coverslips and then adhered to slides with Prolong Glass mounting media.
Confocal and wide-field microscopy.
Confocal imaging was performed using a Nikon Ti Inverted Spinning Disk Microscope equipped with a Yokogawa CSU-W1 Spinning disk scanhead with 50 µm pinhole disk, a PI Z piezo stage insert, a TOPTICA iChrome MLE laser launch, and a Plan Apo λ 60× 1.4 numerical aperture (NA) oil objective for bouton neurochemical phenotyping and RNAscope experiments. For retrograde tracing experiments, a 20× Plan Apo λ 20×/0.75 NA Air differential interference contrast (DIC) objective was used to generate tiled images. We also used this objective to quantify cell types targeted by r2-Pet1 baskets and to quantify the proportion of baskets that formed composites with other Pet1 fibers. Images were captured using an Andor Zyla 4.2 Plus sCMOS monochrome camera and Nikon Elements Acquisition Software AR 5.02. In all colocalization experiments, single-label controls were used to determine imaging settings to eliminate channel bleedthrough. Signal from the different channels was acquired sequentially using a Semrock Di01-T405/488/568/647 multi-bandpass dichroic mirror and bandpass emission filters for green (Chroma ET525/36m), red (Chroma ET 605/52m), and far red (Chroma ET 705/72m) channels. The z-stack acquisition was optimized to minimize axial shift between channels.
Several pericellular basket example images were taken using a laser-scanning confocal Zeiss Axio Observer Z1 equipped with LSM780 Scan Head, motorized stage, and with either a Plan Apo 63×/1.4 Oil DIC III or Plan Apo 40×/1.3 Oil DIC III objective. These images were acquired using 405, 488, 561, and 633 nm laser lines, and a signal from each channel was acquired sequentially using Quasar PMT and GAaSP 32 Channel Spectral detectors. The range of wavelengths collected for each channel was determined using the Zeiss Zen Black (2011) software Smart Setup settings and compared with negative control slides. Voxel size was set to 0.21 × 0.21 × 1.13 μm.
Wide-field fluorescent micrographs were collected using a Zeiss Axioplan2 Upright Microscope with manual stage and an Axiocam 506 monochrome camera. Samples were illuminated with a Lumencor SOLA LED light engine. Signal from the different channels was acquired sequentially with FITC, DsRed, and Cy5 filter sets and used either a 10×/0.45 NA Air Plan Apochromat or 20×/0.8 Air Plan Apochromat objective. For all images, when brightness/contrast settings were adjusted using ImageJ, identical settings were used for each channel applied uniformly across images within a given experiment. Images are presented as maximum intensity z-stack projections generated using Fiji (NIH; Schneider et al., 2012).
Bouton neurochemical immunophenotyping analysis and quantification.
Histologic studies for bouton neurochemical phenotyping were performed using tissue collected from P50 to P100 male (n = 5) and female (n = 6) triple-transgenic r2Hoxa2-cre, Pet1-Flpe, RC-FPSit mice from multiple litters. To determine whether r2-Pet1 neuron boutons contained VGLUT3, 5-HT, or both, tissue was stained for EGFP, VGLUT3, and 5-HT (see above for antibodies used). We collected four images of each ROI per animal when possible (mean number of images per region per animal = 3.82) using a Nikon Ti spinning disk microscope with 1.4 NA 60× objective (see above). Colocalization of synGFP boutons with VGLUT3 and/or 5-HT was performed in a semiautomated manner in Imaris [Bitplane version 9.1.2 (January 29, 2018) Build 45902]. Briefly, each channel (synGFP, 5-HT, and VGLUT3) was converted into a spots object based on user-defined thresholds. If present, visual artifacts [e.g., highly autofluorescent cells of the olfactory bulb (OB)] were manually removed. The MATLAB (version 2017b) plugin “Colocalize Spots” was used to test for colocalization based on distance between synGFP centroids and 5-HT or VGLUT3 centroids. A distance threshold of 1 µm was chosen empirically based on positive control (staining for synGFP with two secondary antibodies conjugated to different fluorophores) and negative control (randomly shuffling the synGFP channels of test images of a particular region to mismatch synGFP and VGLUT3 or synGFP and 5-HT signals) experiments. In these positive and negative control experiments, these analysis parameters resulted in a 1.03% false-negative rate and 2.16% false-positive rate (Fig. 2A–C). The same approach was used to neurochemically phenotype r2-Pet1 and other Pet1 boutons contributing to baskets in the lateral septum using images from adult r2Hoxa2-cre, Pet1-Flpe, RC-FrePe mice (n = 3).
Single-molecule mRNA FISH transcript quantification.
For semiautomated analysis of mRNA transcript number in r2-Pet1 neurons, three-color smFISH was performed to detect eGFP, Pet1, and Tph2 in triple-transgenic r2Hoxa2-cre, Pet1-Flpe, RC-FrePe mice from multiple litters. ImageJ was used to perform a 50 pixel rolling-ball background subtraction, to adjust window/level settings, and to generate threshold-transformed masks of each fluorescence channel. For each eGFP-marked r2-Pet1 cell, Tph2 and Pet1 transcripts were automatically quantified using the Analyze Particles function in Fiji (Schindelin et al., 2012) with a particle size of 1–10 pixels. Alternatively, a custom ImageJ analysis script (Okaty et al., 2020) was used to count Vglut3 and Tph2 mRNA puncta within EGFP+ cell somata using the Default or MaxEntropy thresholds to segment cell somata and the Find Maxima function with a noise tolerance of 1000 to identify mRNA puncta. Analysis was performed on maximum intensity projections of z-stacks as this method reduced computational load and time, and produced results that were highly similar to a full 3D counting approach in Imaris (data not shown).
Quantification of cell somata in retrograde tracing experiments.
For retrograde tracing experiments, images from virus-injected brains were analyzed by manual counting of EGFP+ and tdTomato+ cells using the Cell Counter plugin for Fiji (NIH). First, we selected display settings using single-label controls. We quantified the percentage of cells expressing VGLUT3 or 5-HT for tdTomato+/EGFP+ cells as well as for all tdTomato+ cells and all EGFP+ cells.
Quantification of Pet1 composite baskets in septum.
To quantify the proportion of septal cells targeted by r2-Pet1 baskets that also received boutons from non-r2 Pet1 cells, we analyzed triple-transgenic n = 4 (2 male, 2 female) r2Hoxa2-cre, Pet1-Flpe, RC-FrePe mice and 1 r2Hoxa2-cre, Pet1-Flpe, RC-FLTG (male) mouse. First, we manually identified EGFP+ r2-Pet1 baskets using only the green channel; then, each encased cell was manually reviewed for whether it was also targeted by mCherry+ non-r2 Pet1 fibers. We quantified the proportion of r2-Pet1-targeted somata receiving (1) a non-r2 Pet1 basket (thus receiving a composite basket), (2) non-r2 Pet1 boutons that did not meet basket criteria, or (3) receiving no boutons from other Pet1 cells.
Quantification of cell somata targeted by r2-Pet1 baskets.
To quantify the proportion of VGLUT3+ and calbindin+ somata encased by r2-Pet1 baskets in the cortex, hippocampus, and septum, we executed the following. For VGLUT3+ somata and for calbindin+ somata in hippocampus and septum, a user manually segmented cells in Fiji using the magic wand tool on a median-filtered image (radius: 5 for calbindin and 10 for VGLUT3). For calbindin+ cells in the cortex, cell somata were amenable to automated segmentation using a custom Fiji macro. Briefly, the macro first corrected for uneven illumination by subtracting from the calbindin channel a Gaussian blurred image (radius, 100) of that channel. This was followed by a minor Gaussian blur (radius, 2) and autothresholding using the “Triangle” method. Resulting cell objects were filtered using the Analyze Particles function to limit cells to sizes between 40 and 300 µm2 area and 0.4–1.00 circularity; clumped cells were split using the Adjustable Watershed plugin with a tolerance of 1. Finally, these images were manually reviewed, and any autofluorescent vasculature incorrectly counted by the macro as cells were removed. To examine cells in specific subregions of each image, another custom script sorted cells into user-defined regions of interest (ROIs). After segmentation, another custom script was used to quantify the area of r2-Pet1 boutons overlapping the quadrants of each cell body. Each cell was divided into quadrants with a gap at the center to prevent a bouton located at the center to be counted in all quadrants. The script then segmented r2-Pet1 fibers following the same automated method above (but without filtering for cell-sized objects) and quantified the bouton area in each quadrant. Cells were counted as receiving a basket if they had innervation in three or more quadrants of at least 1 µm2 each and an overall innervation of >12 µm2. These parameters produced data that corresponded well to manually counting basket. All code is publicly available on Github in the repository “auto-basket-detector-2D” by user rsenft1.
Statistical analyses
Statistical tests were performed in Prism 8 (GraphPad Software). To explore the relationship between Tph2 and Pet1 expression, we performed a Spearman correlation on cells with countable Tph2 levels. To compare the neurochemical phenotypes of boutons in baskets against those outside baskets, we performed a within-animal two-tailed paired t test. To compare the percentage of hippocampal-projecting r2-Pet1 cells stained for 5-HT versus the percentage stained for VGLUT3, we performed a paired t test. We also used paired t tests to compare the neurochemical phenotypes of r2-Pet1 versus other Pet1 boutons contributing to composite baskets and to compare r2-Pet1 basket targets across regions. Results are reported as the mean ± SEM, unless otherwise stated.
Results
r2-Pet1 somata stratify by Tph2 versus Vglut3 mRNA levels even with abundant Pet1 transcripts
Previous molecular analyses of the r2-Pet1 neuron population suggested two categories of constituent cells based on neurotransmitter-related transcript levels: one cell subgroup was abundant in Tph2 mRNA; and the other was abundant in Vglut3 mRNA but demonstrating low or no detectable Tph2 transcripts (Okaty et al., 2015). This inverse relationship between Tph2 and Vglut3 mRNA abundance characterizing these two Pet1-expressing cell subgroups, referred to as r2-Pet1Tph2-high and r2-Pet1Vglut3-high, contrasts with most other regions of the Pet1+ brainstem raphe. More typical among raphe Pet1+ cells are abundant levels of both Tph2 and Vglut3, or solely Tph2, with these serotonergic (TPH2+) cells intermingled among various non-Pet1 raphe neurons, many of which are themselves glutamatergic. The distinctive, largely “one-or-the-other” neurochemical stratification distinguishing r2-Pet1 cell bodies suggests that the cognate axonal boutons might also show exclusivity for neurochemical phenotype (5-HT vs VGLUT3 indicative of glutamate release; Amilhon et al., 2010) and thus provide a means to resolve the efferent circuitry for each subgroup from within the broad innervation profile mapped for the r2-Pet1 population en masse (Bang et al., 2012). Before embarking on bouton phenotyping, however, we first set out to re-examine Tph2, Vglut3, and Pet1 transcript levels in r2-Pet1 neuron somata using the sensitive and specific approach of multiplexed FISH for mRNA detection.
We first performed single-molecule transcript analyses in tissue sections from triple-transgenic (r2Hoxa2-cre, Pet1-Flpe, RC-FrePe) mice in which EGFP protein and mRNA expression selectively marked the r2-Pet1 neuron population (Fig. 1A–D), on top of which Pet1 and Tph2 mRNA could be visualized (Fig. 1D). As expected, given the use of the Pet1-Flpe driver, all eGFP transcript-marked cells showed Pet1 transcripts (fluorescent puncta), with a mean absolute count of 59.4 ± 4.6 Pet1 mRNA puncta per cell, corresponding to an absolute range of 10–145 Pet1 transcripts/cell (n = 54 cells across three mice; Fig. 1E). This finding validates that r2-Pet1 neurons, indeed constituents of both subgroups, continue to express Pet1 in adulthood. Quantitative analyses of puncta representing individual Tph2 mRNA transcripts in these cells normalized to cell area showed that ∼30% of EGFP-marked r2-Pet1 neurons (16 of 54 cells) expressed high levels of Tph2 (Tph2 transcripts too numerous to count), while the remainder of cells (38) had few or no Tph2 puncta (Fig. 1E′). For cells with countable levels of Tph2 puncta, area-normalized Tph2 transcript levels exhibited a positive correlation with Pet1 transcript counts (Spearman correlation, r = 0.3413; n = 38 cells, p = 0.036). Overall, these findings are consistent with our prior r2-Pet1 single-cell RNA sequencing (scRNAseq) results (Okaty et al., 2015).
A, The developmentally specified r2-Pet1 neuron lineage comprises two transcriptomically polarized subtypes, one appearing more serotonergic (Tph2-high, Vglut3-low) and the other more glutamatergic (Tph2-low, Vglut3-high). r2-Pet1 cells are derived from progenitor cells in the r2 domain of the embryonic hindbrain that postmitotically express Pet1 and are intermingled with other lineages in the adult brainstem median raphe (B8 and B5) and B9 nuclei. B, C, We genetically access r2-Pet1 cells using ROSA26 (R26) knock-in dual Cre- and Flp-responsive intersectional alleles partnered with the r2- and serotonin lineage-specific drivers r2Hoxa2-cre and Pet1-Flpe, respectively. D–E′, By fluorescence in situ hybridization performed on brain sections (D), all EGFP-labeled r2-Pet1 cells expressed detectable levels of Pet1 transcripts (E), and Pet1 mRNA levels in r2-Pet1 cells weakly correlated with Tph2 transcript abundance (r = 0.3413; n = 38 cells, p = 0.036, Spearman correlation; E′). F–H, r2-Pet1 cells comprise transcriptomically distinct subtypes, one appearing more serotonergic and expressing high levels of Tph2 and low levels of Vglut3 (r2-Pet1Tph2-high), and one appearing more glutamatergic with the inverse pattern of gene expression (r2-Pet1Vglut3-high) as seen in RNA FISH-IHC fluorescent photomicrographs of tissue sections in which EGFP labeling marks r2-Pet1 neurons (F) and quantification of puncta counts per r2-Pet1 cell (G, H). Coronal brain schematics for all figures adapted from Franklin and Paxinos (2008) with permission (copyright Academic Press-Elsevier).
Next, we explored the relationship between Tph2 and Vglut3 transcript abundance in the r2-Pet1 cells, using RNAscope FISH on tissue sections from a triple transgenic (r2Hoxa2-cre, Pet1-Flpe, RC-FrePe) mouse, but now visualizing the r2-Pet1 cells by immunofluorescent detection of EGFP protein (Fig. 1F,F′). We counted Tph2 and Vglut3 puncta per EGFP+ cell (n = 252 r2-Pet1 cells) and calculated the base 2 logarithm of the Vglut3/Tph2 transcript ratio. If ≥1, we classified those cells as belonging to the r2-Pet1Vglut3-high group, and if −1 or less, to the r2-Pet1Tph2-high group (Fig. 1G,H), meaning that a ratio of at least 2:1 in transcript count is necessary to be classified as either cell type. By these criteria, we found 172 (68.3%) cells as r2-Pet1Vglut3-high, and 65 (25.8%) cells as r2-Pet1Tph2-high. We also examined transcript counts normalized to cell area for cells classified into these groups, finding that r2-Pet1Tph2-high cells exhibited high Tph2 expression (mean, 0.9438 transcripts/µm2; 95% CI, 0.86–1.03) and low Vglut3 expression (mean, 0.1096 transcripts/µm2; 95% CI, 0.084–0.14), while r2-Pet1Vglut3-high cells expressed high Vglut3 (mean, 0.6308 transcripts/µm2; 95% CI, 0.58–0.69) and low Tph2 (mean, 0.063 transcripts/µm2; 95% CI, 0.051–0.075). A small subset of r2-Pet1 cells (5.9%) contained an intermediate Vglut3/Tph2 transcript ratio or vice versa (i.e., a <2:1 comparative enrichment of either Vglut3 or Tph2. For simplicity, the present study approximates the molecular heterogeneity of r2-Pet1 neurons using the neuron subgroups descriptors of r2-Pet1Vglut3-high and r2-Pet1Tph2-high cells, with the caveat that a small minority of r2-Pet1 neurons likely harbor an intermediate phenotype. Collectively, these RNAscope in situ hybridization findings are in line with the previously reported scRNAseq datasets (Okaty et al., 2015), confirming that the r2-Pet1 neuron group shows a segregation of somata into two subgroups based on expression of neurochemical pathway genes.
Segregation anatomically of r2-Pet1 boutons by serotonin+ versus VGLUT3+ neurochemical identity
To visualize and subsequently explore the anatomic distribution and neurochemical identity of r2-Pet1 neuron axonal boutons in the adult brain, we used the intersectional RC-FPSit allele (Niederkofler et al., 2016), which drives the expression of a synaptophysin–GFP (synGFP) fusion protein selectively in cells with a history of both Cre and Flp expression. Synaptophysin, a synaptic vesicle protein, enriches the fused GFP signal in both axon terminal and en passant boutons (De Paola et al., 2003; Kelsch et al., 2008; Li et al., 2010; Tripodi et al., 2011); both types of boutons have been shown to participate in vesicle-mediated neurotransmission (Beaudet and Descarries, 1981; Agnati et al., 1995; De-Miguel and Trueta, 2005; Trueta and De-Miguel, 2012). In triple-transgenic (r2Hoxa2-cre, Pet1-Flpe, RC-FPSit) mice, the synGFP robustly and selectively illuminated r2-Pet1 neuron axonal boutons compared with interbouton (i.e., intervaricose) axon segments. We observed the densest r2-Pet1 bouton distribution in the following brain areas: OB (glomerular layer), prelimbic cortex (superficial layers 1–3), suprachiasmatic nucleus of the hypothalamus (SCN), paraventricular nucleus of the thalamus (PVT), HP [present throughout many layers and highest at the stratum radiatum-stratum lacunosum moleculare (sr/slm) border], medial and lateral septal nuclei (MS and LS, respectively), and diagonal band of Broca. This pattern of regions corroborates previous reports of fiber distribution from r2-Pet1 neurons (Bang et al., 2012).
Next, we probed synGFP-marked r2-Pet1 boutons for coimmunolocalization with 5-HT and VGLUT3. Brain sections (every eighth 40 µm serial section from five male and six female P60 r2Hoxa2-cre, Pet1-Flpe, RC-FPSit mice) were immunostained and imaged on a confocal microscope, and the resulting images were analyzed using Imaris (Bitplane) to calculate the fraction of synGFP+ boutons that were also immunopositive for 5-HT, VGLUT3, both, or neither, for four images per region per mouse (see Materials and Methods). We analyzed 66,435 synGFP+ boutons across 13 forebrain subregions, for an average of 5,110 boutons per region and 465 boutons per region per mouse. Generally, we observed a segregation of bouton neurochemical identity by anatomic target site. In the glomerular layer of the olfactory bulb, the PVT, and the SCN, the r2-Pet1 boutons were largely immunopositive for 5-HT but not VGLUT3 (Fig. 2D–F). Furthermore, in the PVT, approximately half of the GFP-marked boutons lacked immunodetectable levels of even 5-HT. By contrast, dorsal cortical and hippocampal regions were largely innervated by VGLUT3-only r2-Pet1 boutons (86.3 ± 0.94%, across the seven regions; Fig. 3A–C). Few 5-HT+ r2-Pet1 boutons were observed (1.35 ± 0.51% of all boutons across these regions), and a modest proportion lacked staining for either (9.38 ± 0.68% across these regions). These regions harbored synGFP–, 5-HT+ fibers, suggesting that these areas are also innervated by other (non-r2-derived) subtypes of serotonergic neurons, which is not surprising given that r1- and r3-derived Pet1 serotonergic neuron cell bodies are also present in the MR (Jensen et al., 2008), and that dorsal raphe (DR) efferents also innervate these target regions (Azmitia and Segal, 1978). These results show that r2-Pet1 cells are in part responsible for the largely VGLUT3+ innervation of the hippocampus and cortex previously observed as arising from the median raphe (Jackson et al., 2008; Szőnyi et al., 2016). In the dorsolateral septum (DLS) and medial septum (MS) and caudal stratum radiatum of the CA3 hippocampal subfield (cCA3-SR), the r2-Pet1 boutons were largely VGLUT3 only (66 ± 8.6%), but ∼20% stained positive for both 5-HT and VGLUT3 (20.4 ± 5.2%; Fig. 4A–C). Furthermore, in the DLS, we also found 5-HT-only r2-Pet1 boutons (8.04 ± 1.6%), making the DLS unique among examined regions in harboring all four types of r2-Pet1 boutons—5-HT-only, VGLUT3-only, co-positive, and negative 5-HT and VGLUT3. The DLS, and the septum more generally, also exhibited extensive non-r2-Pet1 5-HT+ and VGLUT3+ innervation, with many of these fibers immunopositive for both 5-HT and VGLUT3, suggesting another lineage of Pet1 cells also innervates the septum and may be able to cotransmit 5-HT and glutamate.
The OB, SCN, and PVT are innervated by 5-HT+ r2-Pet1 boutons. A, B, Example photomicrograph of a colocalized bouton (A) and criteria for object-based colocalization (B). C–C′′, Positive control labeling of EGFP with multiple fluorophores (C) and a negative shuffled channel control (C′) were used to assess false-negative and false-positive rates, which are quantified in C′′. D, Fluorescent photomicrograph from the olfactory bulb glomerular layer showing r2-Pet1 boutons labeled by 5-HT (arrowheads demarcate two of many) but not VGLUT3. E, F, Coronal schematics with regions demonstrating this innervation pattern in blue (E) and quantification of colocalization (F).
The cortex, HP, and diagonal band are innervated by VGLUT3+ r2-Pet1 boutons. A, Fluorescent photomicrograph from dorsal hippocampal CA1 sr/slm showing r2-Pet1 boutons labeled by VGLUT3 but not 5-HT (arrowheads show examples). B, C, Coronal schematics with regions demonstrating this innervation pattern shown in red (B) and quantification of colocalization (C).
The medial and dorsolateral septum and caudal hippocampus are innervated by r2-Pet1 boutons that are either labeled by VGLUT3 alone or are 5-HT and VGLUT3 co-positive. A, Fluorescent photomicrograph from the medial septum showing VGLUT3+ r2-Pet1 boutons (open arrowheads as examples) and 5-HT+/VGLUT3+ r2-Pet1 boutons (filled arrowhead). B, C, Coronal schematics with regions demonstrating this innervation pattern in magenta (B) and quantification of colocalization (C).
Driver Piezo2-cre coupled with Pet1-Flpe captured MR boutons positive for serotonin but not VGLUT3
Given that the r2-Pet1 boutons found in the OB, SCN, and PVT were largely of the 5-HT-only class, we hypothesized that the source neuron subset was the r2-Pet1Tph2-high neuron subpopulation and not the r2-Pet1Vglut3-high subgroup. Alternatively, it is possible that a single r2-Pet1 neuron projects both serotonergic and glutamatergic boutons, for example, through differential axonal trafficking of mRNA and/or protein. To test our hypothesis, we devised means to intersectionally label the r2-Pet1Tph2-high neuron subpopulation. Single-cell transcriptomic analysis of Pet1 MR neurons showed Piezo2 transcripts selectively enriched in r2-Pet1Tph2-high neurons compared with other Pet1 neurons (Okaty et al., 2015). This was further supported in the study by Ren et al. (2019) with the characterization of MR cluster 4, a Piezo2+ Tph2-high MR population expressing similarly modest levels of Vglut3 as the r2-Pet1Tph2-high group. Thus, we generated triple-transgenic Piezo2-GFP-IRES-Cre, Pet1-Flpe, RC-Ai65 mice to illuminate somata and axons of r2-Pet1Tph2-high neurons selectively, which would also reveal the efferent target regions and allow neurochemical phenotyping for comparison with our generated bouton neurochemical maps for r2-Pet1 cells.
Cell bodies labeled in Piezo2-GFP-IRES-Cre, Pet1-Flpe, RC-Ai65 animals were restricted to the MR and were largely immunopositive for 5-HT (95.2 ± 4.76%), with a small fraction of cells copositive for VGLUT3 and 5-HT (4.76 ± 4.76%) and no cells immunopositive for VGLUT3 in the absence of 5-HT (0.00 ± 0%; Fig. 5A). We observed fewer Piezo2-Pet1 neurons than would be expected of the entire r2-Pet1Tph2-high population, and cell bodies were restricted to the midline, suggesting that they are a subset of the r2-Pet1Tph2-high group, or at least a subset accessible using this recombinase driver combination. We observed dense, 5-HT+ Piezo2-Pet1 fibers within the SCN (73.2 ± 3.27% 5-HT+) and the PVT (60.5 ± 3.16% 5-HT+; Fig. 5B,C). Piezo2-Pet1 fibers were present, but less dense, in other regions where we observed 5-HT+ r2-Pet1 innervation such as OB, LS, and MS (Fig. 5D–F). In these regions, the majority of Piezo2-Pet1 fibers were immunopositive for 5-HT+ (OB, 74.7 ± 6.01% 5-HT+; LS, 56.2 ± 4.97%; MS, 71.6 ± 1.90%) but negative for VGLUT3 (OB, 0.250 ± 0.126% VGLUT3+; LS, 1.04 ± 0.536%; MS, 1.342 ± 0.514%). Across all regions, a small percentage of Piezo2-Pet1 boutons were immunopositive for both VGLUT3 and 5-HT, highest in the OB (2.98 ± 0.201% VGLUT3+/5-HT+) and SCN (2.87 ± 0.437% VGLUT3+/5-HT+). Upon exploring the target areas harboring VGLUT3-only boutons from r2-Pet1 neurons, we failed to detect innervation by the Piezo2-Pet1 neuron group (dorsal hippocampus shown; Fig. 5G). Together, these findings suggest that Piezo2-Pet1 neurons represent a r2-Pet1Tph2-high population with projections largely restricted to the 5-HT-only SCN, PVT, and OB.
Piezo2-Pet1 neurons are of the r2-Pet1Tph2-high subtype and harbor projections selectively to regions innervated by 5-HT+ r2-Pet1 boutons. We previously identified Piezo2 as a marker of the r2-Pet1Tph2-high subtype (Okaty et al., 2015). A, Piezo2-Pet1 cell bodies in the MR are labeled by 5-HT but not VGLUT3. B–F, Piezo2-Pet1 neurons give rise to 5-HT+ projections (arrowheads) in the SCN (B), PVT (C), lateral septum (D), medial septum (E), and olfactory bulb glomerular layer (F). G, By contrast, Piezo2-Pet1 cells do not innervate the dorsal hippocampus, a region shown to have predominately VGLUT3+ r2-Pet1 innervation.
Retrograde viral tracing from the dorsal hippocampus highlighted r2-Pet1 somata positive for VGLUT3 but not serotonin
Given that the r2-Pet1 boutons found in the dorsal hippocampus were of the VGLUT3-only class, we hypothesized that the source neuron subset was the r2-Pet1Vglut3-high neuron subpopulation and not the r2-Pet1Tph2-high subgroup. To test this, we injected bilaterally into the dorsal hippocampus (AP, −2.06; ML, ±2.06; DV, −1.75 mm) of r2Hoxa2-cre, Pet1-Flpe, RC-FrePe mice (n = 3) a retrogradely transported, tdTomato-expressing adeno-associated virus (AAV) and analyzed the neurochemical phenotype (staining for either 5-HT or VGLUT3) of the tdTomato-labeled (and thus dorsal hippocampus-projecting) subset of r2-Pet1 (EGFP+) somata (Fig. 6). Worth noting, VGLUT3, in contrast to other VGLUTs, localizes both to axonal boutons and the somatodendritic compartment as well (Fremeau et al., 2004; Herzog et al., 2004), which we exploited here. To identify VGLUT3+ cell bodies, we used a knock-out mouse-verified antibody (Schallier et al., 2011) that was previously used to label VGLUT3+ somata in the raphe (Sos et al., 2017; Okaty et al., 2020) and found a distribution of VGLUT3-immunolabeled cell bodies consistent with prior anatomic studies describing the colocalization of Vglut3 mRNA and its protein in MR neurons (Herzog et al., 2004). We examined 1804 neurons, with 1263 EGFP+ neurons, 629 tdTomato+ neurons, and 88 neurons colabeled for tdTomato and EGFP, with the latter indicating that our retrograde tracing conditions labeled a fraction of the r2-Pet1 cell group. Of the dual-labeled (EGFP+, tdTomato+) r2-Pet1 somata, 82.3 ± 1.46% stained immunopositive for VGLUT3, but only 2.0 ± 2.0% stained positive for 5-HT. Thus, the hippocampal-projecting r2-Pet1 somata were significantly more likely to be VGLUT3+ than 5-HT+ (t = 106.6, df = 2, p < 0.0001, paired t test; Fig. 6E). This finding is in line with the r2-Pet1 bouton immunophenotype of VGLUT3-only mapped to the dorsal hippocampus. Of all retrogradely labeled tdTomato+ cells in the MR region, most were VGLUT3+ (76.8 ± 6.6%) and few were 5-HT+ (9.6 ± 3.86%), consistent with previous literature reporting a predominant VGLUT3+ innervation of the hippocampus from the MR (Jackson et al., 2008; Szőnyi et al., 2016; Fig. 6E′). In contrast, when sampling all r2-Pet1 cells (EGFP labeled), a smaller fraction of cells was immunoreactive for VGLUT3 (52.0 ± 2.19%) and 28.74 ± 0.94% were 5-HT+ (Fig. 6E′′). These results demonstrate that the r2-Pet1 cells that project to the dorsal hippocampus overwhelmingly express VGLUT3 in the soma, consistent with the neurochemical phenotype of the mapped r2-Pet1 boutons, suggesting that the source neurons are the r2-Pet1Vglut3-high.
Hippocampal retrograde tracing shows that VGLUT3+ boutons in this region preferentially arise from r2-Pet1Vglut3-high neurons. A, Schematic showing injection, in the dorsal hippocampus (bilaterally) at the CA1 sr/slm border, of a retrograde serotype AAV expressing tdTomato under a CAG promoter performed in r2Hoxa2-cre, Pet1-Flpe, RC-FrePe mice expressing EGFP in r2-Pet1 cells. B, Fluorescent photomicrograph showing the general localization of retrogradely labeled cells (tdTomato+) and r2-Pet1 cells (EGFP+) in the MR. C, D, r2-Pet1 cells that project to the hippocampus (EGFP+/tdTomato+) have somata immunopositive for VGLUT3 (C) but not 5-HT (D). Hippocampus-projecting r2-Pet1 cells were more likely to be VGLUT3+ than 5-HT+ (E, t = 106.6, df = 2, p < 0.0001, paired t-test). We also quantified VGLUT3 and 5-HT immunopositivity for all tdTomato+ cells in the MR (E′) or the entire population of r2-Pet1 cells (E′′). Asterisks indicate statistically significant differences.
r2-Pet1 VGLUT3+ boutons comprise pericellular baskets in certain brain regions
In brain regions with r2-Pet1 VGLUT3+ innervation, including the cortex (retrosplenial, cingulate, prelimbic, somatosensory, and motor cortex), hippocampus, and septal nuclei, we found r2-Pet1 boutons in dense arrays surrounding a target cell body and proximal dendrites, akin to previously described “pericellular baskets” (Fig. 7A–D; Köhler et al., 1982; Freund et al., 1990; Hornung et al., 1990; DeFelipe et al., 1991; Vu and Törk, 1992; Dinopoulos et al., 1993; Hioki et al., 2004; Riedel et al., 2008). Operationally, we defined r2-Pet1 pericellular baskets when r2-Pet1 axons wrapped around more than half of the soma perimeter (or when no cell stain was present, the nucleus perimeter) or when boutons were closely apposing the soma on multiple sides such that dividing the cell into quadrants would result in boutons apposing the cell in all quadrants. Such pericellular basket structures rarely characterized r2-Pet1 boutons found in the OB and PVT and not at all in the SCN—all target regions harboring r2-Pet1 boutons of the 5-HT-only class. Additionally, Piezo2-Pet1 fibers representing a subset of the r2-Pet1Tph2-high subtype were not observed to make pericellular baskets. Thus, this specialized presynaptic basket structure may be feature specific to r2-Pet1Vglut3-high cells, but not r2-Pet1Tph2-high cells.
r2-Pet1 neuron boutons contribute to pericellular basket arrays in the cortex, hippocampus, and septum. A, Fluorescent photomicrograph showing a basket comprising r2-Pet1 boutons surrounding the soma and proximal dendrites of a hippocampal CA3 cell retrogradely labeled by injection of AAV2-retro-CAG-tdTomato-WPRE into CA1. B, C, Fluorescent photomicrographs showing pericellular basket arrays in the motor cortex (B) and hippocampal subgranular zone (SGZ) (C) comprising VGLUT3+ r2-Pet1 boutons. D, By contrast, r2-Pet1 basket arrays in the medial septum comprise boutons that are 5-HT–/VGLUT3+ as well as copositive 5-HT+/VGLUT3+. D′, In the medial septum, co-positive 5-HT+/VGLUT3+ boutons are enriched in pericellular baskets, (t = 3.140, df = 6, p = 0.020, paired t test). Operationally, we defined r2-Pet1 pericellular baskets as occurring when r2-Pet1 axons surrounded >50% of the perimeter of the cell body or nucleus and/or r2-Pet1 boutons closely apposed the target cell body in all four surrounding quadrants. Arrowheads are placed above the basket-targeted cell body.
We show in Figure 4 that r2-Pet1 boutons in the septum are either VGLUT3+ or VGLUT3+ and 5-HT+. We examined a subset of 11 medial septum images from r2Hoxa2-cre, Pet1-Flpe, RC-FPSit mice (n = 7), quantified in Figures 2–4, and examined the neurochemical phenotypes of boutons belonging to baskets (manually defined as above) versus other boutons in the image innervating the same region but not contributing to baskets. Of the r2-Pet1 boutons comprising pericellular baskets (473 boutons examined across all mice), 31.98% were co-positive, versus only 15.9% of boutons in nonbasket configurations (1300 nonbasket boutons examined across all mice; Fig. 7D′), showing that r2-Pet1 boutons that contribute to baskets in the medial septum were significantly more likely to be co-positive for both VGLUT3 and 5-HT than nonbasket boutons in the same region (t = 3.140, df = 6, p = 0.020, paired t test). These results indicate that in the septum, boutons that contribute to baskets may be more likely to cotransmit 5-HT and glutamate.
Some pericellular baskets are comprised of boutons from multiple Pet1 neuron subtypes exhibiting distinct neurochemical profiles
In some cases, our 5-HT and synGFP staining revealed synGFP–, 5-HT+ fibers converging on the same target cell as the synGFP+ r2-Pet1 boutons comprising a basket (Fig. 8A). This suggested that the former might arise from a non-r2-derived Pet1 serotonergic neuron (i.e., a different subtype of serotonergic neuron). Such composite baskets (r2-Pet1 VGLUT3+ boutons and other 5-HT+ boutons) were found to be most prevalent in the stratum radiatum of the hippocampal CA3 subfield (Fig. 8A) and in the medial and lateral septal nuclei (Fig. 8B). In line with the observation that the non-r2-Pet1 fibers in these composite baskets were found to be 5-HT+ (or SERT+), we discovered that some of the postsynaptic target cells in the cortex and dorsal hippocampus expressed 5-HT3A (Fig. 8C), the ionotropic excitatory 5-HT receptor, and thus could be capable of receiving the basket-derived 5-HT signal. Moreover, the nature of the 5-HT3A receptor suggests that 5-HT may evoke fast excitatory synaptic firing (Lee et al., 2010) in these specific target cells. In the lateral septum, some basket-targeted cells were immunopositive for the slower signaling G-protein-coupled receptor 5-HT2C (Fig. 8D); none were immunopositive for 5-HT3A (data not shown). Together, these findings suggest that 5-HT released from pericellular basket-associated boutons may act on distinct temporal scales across different brain regions.
Cells targeted by r2-Pet1 pericellular baskets also receive 5-HT input. A, B, Fluorescent photomicrographs in hippocampus (A) and dorsolateral septum (B) showing baskets composed of both synGFP+ r2-Pet1 boutons and non-synGFP-labeled 5-HT+ boutons. C, D, Consistent with the non-r2-Pet1 serotonergic innervation, some basket target cells expressed the 5-HT receptors 5-HT3A (C) and 5-HT2C (D). Arrowheads are placed above nuclei targeted by the r2-Pet1 pericellular basket.
To assess whether multiple Pet1 neuron subtypes converge to form pericellular baskets around the same downstream cell, we labeled simultaneously r2-Pet1 fibers versus other (non-r2) Pet1 fibers with different fluorophores and examined basket composition in the septum, a region abundant in baskets. Triple-transgenic r2Hoxa2-cre, Pet1-Flpe, RC-FrePe (n = 4, 2 males, 2 females) mice and one r2Hoxa2-cre, Pet1-Flpe, RC-FLTG (n = 1, male) mouse equally enabled this analysis, as each strain selectively marked r2-Pet1 somata and fibers intersectionally with EGFP while simultaneously labeling other Pet1 cells with mCherry or tdTomato, respectively, as the subtractive population (Fig. 9A,A′). We examined 163 r2-Pet1 pericellular baskets, finding many with nearby subtractive Pet1 boutons closely apposing the same cell (64.5 ± 3.95% of r2-Pet1 baskets). Further, for 35.1 ± 3.10% of r2-Pet1 baskets, this subtractive Pet1 contribution also formed a pericellular basket (Fig. 9B). For an additional third of r2-Pet1 baskets, there was no other Pet1 innervation closely apposing the targeted cell nucleus (35.5 ± 3.95% of baskets). These results demonstrate that r2-Pet1 pericellular baskets are a heterogeneous group with a large portion converging with other Pet1 neurons to form composite pericellular baskets on the same downstream cell.
Multiple Pet1 neuron lineages contribute to composite pericellular baskets in the septum. A, A′, Fluorescent photomicrographs showing composite baskets in the lateral septum in which r2-Pet1 axons (green) and other non-r2 Pet1 axons (magenta) converge onto the same targeted cell body, indicated by arrowheads. B, Quantification of the proportion of r2-Pet1 baskets that were composite with other Pet1 axons.
We next queried whether r2-Pet1 boutons in composite pericellular baskets differed neurochemically from other Pet1 boutons contributing to the same baskets. We sampled 132 composite baskets (3,294 EGFP boutons and 2,672 mCherry boutons) of the lateral septum from adult r2Hoxa2-cre, Pet1-Flpe, RC-FrePe (n = 3, 2 male, 1 female) mice. In separate experiments, we immunodetected EGFP (r2-Pet1 boutons), mCherry (other Pet1 boutons), and either 5-HT (Fig. 10A) or VGLUT3 (Fig. 10B). Composite baskets showed no differences in amounts of r2-Pet1 EGFP boutons (25.6 ± 1.40) versus other Pet1 mCherry boutons (19.8 ± 1.52). Baskets also contained additional boutons immunopositive for 5-HT (31.4 ± 2.01) or VGLUT3 (32.2 ± 1.05) not captured using our transgenic fluorescent labeling approach. These boutons could be from Pet1 neurons missed by our drivers or immunohistochemical methods or alternatively, could reflect boutons from non-Pet1 neurons that form baskets around these same target cells (Fig. 10C,D). With respect to the Pet1-captured basket-forming fibers, we found the non-r2 Pet1 boutons were more likely to be 5-HT+ (70.6 ± 2.67%) than the r2-Pet1 boutons (45.9 ± 3.42%; t = 12.7, df = 2, p = 0.0062, paired t test; Fig. 10E). Both populations contributing to composite baskets were similarly likely to be VGLUT3+ (r2-Pet1, 65.1 ± 1.80%; other Pet1, 58.6 ± 4.61%; t = 1.02, df = 2, p = 0.415, paired t test; Fig. 10F). Altogether, these findings demonstrate that (1) r2-Pet1 and other Pet1 neurons commonly form composite pericellular baskets around the same downstream target cells; and (2) in the lateral septum composite baskets are neurochemically diverse, with r2-Pet1 boutons less likely to be 5-HT+, but are similarly likely to be VGLUT3+ as other Pet1 fibers contributing to the same baskets.
In pericellular baskets, r2-Pet1 boutons (EGFP+) and other Pet1 boutons (mCherry+) are neurochemically distinct. A, B, Composite baskets with r2-Pet1 boutons and other Pet1 boutons in the lateral septum with colabeling for 5-HT (A) and VGLUT3 (B). Arrowheads demarcate the location of the nucleus targeted by the pericellular basket. C, D, We found similar numbers of EGFP and mCherry boutons contributing to composite Pet1 lateral septum pericellular baskets. E, F, However, r2-Pet1 boutons were less likely to be 5-HT+ (t = 12.7, df = 2, p = 0.0062, paired t test; E), but were similarly likely to be VGLUT3+ (t = 1.02, df = 2, p = 0.415, paired t test; F). The orange, purple, and blue lines indicate individual animal means, and individual measurements of baskets are indicated by differently shaped points in C and D, and by lightly shaded orange, purple, and blue lines in E and F. Asterisks indicate statistically significant differences.
r2-Pet1 pericellular baskets target specific GABAergic neuron subtypes
Next, we applied various cell subclass markers to further explore the identity of cells targeted by r2-Pet1 pericellular baskets. We found that r2-Pet1 pericellular baskets targeted GABA+ neurons in both cortical regions and in the septum (Fig. 11A,B), a finding consistent with previous literature on the pericellular baskets formed by median raphe 5-HT neurons (Halasy et al., 1992; Hornung and Celio, 1992). Further, we found that target cells in the cortex and hippocampus compared with the septum could be differentiated by labeling with markers of specific interneuron classes. In the cortex and hippocampus, we observed r2-Pet1 baskets around Reelin+ and VIP+ interneurons (Fig. 11C,D) as well as VGLUT3+ interneurons that were sometimes positive for NKB (Fig. 11E). The latter group is likely to be the well described cholecystokinin (CCK)+/VGLUT3+ basket cell type, which in turn forms pericellular baskets targeting cortical pyramidal neurons (Somogyi et al., 2004; Fasano et al., 2017; Booker and Vida, 2018; Harris et al., 2018). In contrast to the hippocampus and cortex, the targets of r2-Pet1 pericellular baskets in the septum remain largely mysterious (Riedel et al., 2008), though some of these targeted cells express SST. In all regions, r2-Pet1 fibers do not form baskets around parvalbumin+ interneurons (data not shown), consistent with previous studies from the MR (Freund et al., 1990; Hornung and Celio, 1992; Papp et al., 1999) and suggesting the r2-Pet1 pericellular basket regulates specific subtypes of interneurons and possibly GABAergic projection neurons in the septum.
r2-Pet1 pericellular baskets target interneurons. A, B, Targets of pericellular baskets were overwhelmingly GABA+ in dorsal hippocampus (A) and medial septum (B). C–F, Cortical and hippocampal basket targets were often Reelin+ (C), calbindin+, VIP+ (D), or VGLUT3+/NKB+ (E), whereas subcortical septal basket targets were sometimes SST+ (F). Arrowheads demarcate the cell body targeted by the r2-Pet1 pericellular basket.
We then quantified the prevalence of r2-Pet1 baskets on two interneuron classes in cortex, hippocampus, and septum: calbindin+ and VGLUT3+ somata. We examined sections containing dorsal hippocampus and septum in adult r2Hoxa2-cre, Pet1-Flpe, RC-FPSit mice (n = 6, 3 male and 3 female). Using a combination of manual and automatic segmentation combined with automatic detection of baskets (see Materials and Methods), we quantified the proportion of calbindin+ and VGLUT3+ somata receiving a r2-Pet1 basket. We observed that VGLUT3+ somata were most prevalent in layer 2 and layer 6 in cortex, and in the sr and sr/slm border region in hippocampus, which is in line with previous studies (Fremeau et al., 2004; Fig. 12A). Calbindin+ somata were highly abundant in cortical layer 2/3 and sporadically present in deeper layers, were distributed through the sr and sr/slm of the dorsal hippocampus (Fig. 12B), and were enriched in the lateral septum (Fig. 12C). We quantified the proportion of each cell type targeted by r2-Pet1 baskets and found it varied as a function of region. We observed a higher proportion of cortical VGLUT3+ somata (31.2 ± 3.18%, 343 total cells examined) were targeted by an r2-Pet1 pericellular basket relative to hippocampal VGLUT3+ somata (16.6 ± 4.12%; 329 cells examined; t = 8.60, df = 5, p = 0.0004, paired t test; Fig. 12D). Examining subregions of cortex and hippocampus in which we observed an average of ≥10 cells per animal, we found that the greatest percentage of VGLUT3+ somata was targeted in secondary motor cortex (M2; 39.9 ± 4.29%) compared with other cortical subregions [primary motor cortex (M1), 26.6 ± 6.29; parietal association cortex (PtA), 31.2 ± 7.93%] and without appreciable difference between the two hippocampal subregions analyzed (CA1 sr/slm, 16.3 ± 4.87%; CA3 sr/slm, 19.4 ± 3.90%; Fig. 12D′). Regional targeting of calbindin+ interneurons exhibited the opposite pattern: cortical calbindin+ neurons were rarer targets (0.771 ± 0.127%, 15,331 total cells examined) compared with dorsal hippocampal calbindin+ neurons (21.2 ± 2.71%, 225 total cells examined; t = 7.65, df = 5, p = 0.0006, paired t-test; Fig. 12E). In the septum, there were no VGLUT3+ somata, and calbindin+ somata were rarely targeted by r2-Pet1 baskets (0.229 ± 0.170% of cells targeted, 2919 total cells examined). These results demonstrate target specificity of r2-Pet1 baskets wherein different neuron types were preferentially targeted in different regions.
r2-Pet1 pericellular baskets exhibit regional target cell-type specificity. A, B, Examples of r2-Pet1 pericellular baskets targeting VGLUT3+ somata in cortex (A) and calbindin+ somata in dorsal hippocampus (B). C, Rarely, calbindin+ cells are targeted by r2-Pet1 baskets in the lateral septum. Arrowheads indicate the cell body targeted by the r2-Pet1 basket. D, D′, VGLUT3+ somata were more likely to be targeted by r2-Pet1 baskets in cortex than hippocampus (dHP; t = 8.60, df = 5, p = 0.0004, paired t-test; D), with some variation across primary motor cortex (M1), secondary motor cortex (M2), and parietal association cortex (PtA), and hippocampal subfields CA1 and CA3 sr/slm (D′). E, In contrast, a higher proportion of calbindin+ neurons were targeted by r2-Pet1 baskets in the dHP than the cortex or septum (t = 7.65, df = 5, p = 0.0006, paired t test; E). Asterisks indicate statistically significant differences.
Discussion
Strategy
The r2-Pet1 neuron group is hodologically, functionally, and molecularly distinct from other Pet1+ neurons and itself is composed of two subgroups—the r2-Pet1Tph2-high group being serotonergic and the r2-Pet1Vglut3-high subgroup being largely VGLUT3+ and glutamatergic. We proposed that the cognate axonal boutons would likely show a matching neurochemical phenotype (5-HT+ or VGLUT3+) and thus might be used to resolve downstream circuitry specific to each of these subtypes of r2-Pet1 neurons.
Main findings
Significant findings include the following. (1) r2-Pet1 boutons were found distributed anatomically by neurochemical phenotype. 5-HT+/VGLUT3– boutons were found in the SCN PVT, and OB. The reciprocal bouton identity of VGLUT3+, 5-HT– was found in the hippocampus, prelimbic cortex, and septum. Rarely observed were co-positive boutons or a mix of each single-positive bouton type, the exception being in the septum. (2) Piezo2-Pet1 neurons reside within the r2-Pet1Tph2-high group and exhibited axonal projections to the SCN, PVT, and OB. (3) Retrograde labeling from the hippocampus illuminated somata of r2-Pet1Vglut3-high, but not r2-Pet1Tph2-high cells. Findings 2 and 3 substantiate the bouton-to-neuron subtype matchup of finding 1 and suggest that the r2-Pet1 neuronal group is composed of two distinct neuronal subsystems that modulate different brain regions using different neurotransmitters. They likewise suggest that the majority of r2-Pet1Vglut3-high neurons are excitatory, potentially releasing glutamate as their sole transmitter, even with the expression of Pet1 and having arisen from the same progenitor cell compartment as the classically serotonergic Pet1Tph2-high cells. (4) Some VGLUT3+ r2-Pet1 boutons composed pericellular baskets encasing the postsynaptic neuron soma. (5) Some basket-forming boutons in the septum were both 5-HT+ and VGLUT3+, suggestive of possible cotransmission of 5-HT and glutamate (Amilhon et al., 2010). (6) The r2-Pet1 baskets often encased GABAergic interneurons of the Reelin+, VIP+, or VGLUT3+ classes, in some cases as composites containing fibers from different Pet1 neuron subtypes. (7) Some basket-encased, postsynaptic targets in the septum expressed the metabotropic 5-HT2C receptor, and others in the cortex and hippocampus expressed the ionotropic 5-HT3A receptor, suggesting that the basket-released 5-HT may act on distinct temporal scales in different target brain regions.
Separable neuronal subsystems comprise the r2-Pet1 neuron population of the median raphe
A simple interpretation of our r2-Pet1 bouton neurochemical mapping is a divergent circuitry model whereby the r2-Pet1Tph2-high and r2-Pet1Vglut3-high neuron subgroups innervate different brain regions and deploy different neurotransmitters (Fig. 13). While the specific functions of each subgroup have yet to be determined, this divergent circuitry model supports certain predictions. The first is that the midline subcortical projections characterizing the r2-Pet1Tph2-high subgroup may regulate sensory modulation of circadian behavior. The SCN and PVT targets are each implicated in regulating circadian rhythm and themselves are reciprocally connected structures (Weber et al., 1998; Morin and Meyer-Bernstein, 1999). Loss of MR 5-HT input to the SCN has been shown to alter circadian rhythm and increase the sensitivity to light-induced circadian phase shifting (Bradbury et al., 1997; Morin and Meyer-Bernstein, 1999) Similarly, lesions of the PVT suggest that it too may mediate circadian rhythm shifts in response to light (Salazar-Juárez et al., 2002). Thus, r2-Pet1Tph2-high projections to these regions may suppress changes in circadian rhythm by modulating regional responsiveness to deviations from the expected pattern of illumination.
A, Summary of findings. r2-Pet1 VGLUT3+ and 5-HT+ boutons are spatially segregated across the forebrain, with cortical and hippocampal areas innervated by VGLUT3+ r2-Pet1 boutons; midline structures such as the SCN, OB, and PVT innervated by 5-HT+ r2-Pet1 boutons; and septal areas innervated by a mix of VGLUT3+ and copositive 5-HT+/VGLUT3+ boutons. B, These results suggest functionally separable subtypes of r2-Pet1 cells. C, In the septum, cortex, and hippocampus, regions targeted by the r2-Pet1Vglut3-high subtype, r2-Pet1 boutons formed pericellular baskets targeting inhibitory interneurons, suggesting additional functional specialization of this r2-Pet1 subset.
VGLUT3+ r2-Pet1 boutons, by contrast, were found in limbic regions such as the septum, hippocampus, and prelimbic cortex, suggesting that r2-Pet1Vglut3-high neurons may modulate memory and sensorimotor gating. The medial septum and hippocampus are involved in generating theta rhythm, an oscillating pattern of brain activity that contributes to spatial and episodic memory formation (Buzsáki, 2002; Hasselmo, 2005; Munn et al., 2015). Stimulation of the MR has been shown to desynchronize theta (Vertes, 1981; Vinogradova et al., 1999; Jackson et al., 2008; Wang et al., 2015). Though this effect has previously been attributed to 5-HT release (Assaf and Miller, 1978; Kinney et al., 1996), recent work has suggested there also exists a nonserotonergic mechanism (Jackson et al., 2008). If r2-Pet1Vglut3-high neurons desynchronize theta, their activity may perturb memory formation (Vertes, 2005). Indeed, in line with this hypothesis, recent experiments acutely inhibiting r2-Pet1 neurons en masse during conditioned learning led to a later failure in the normal extinction of that conditioned place preference (Baskin et al., 2020). Finally, chronic silencing of r2-Pet1 neurons has been shown to enhance prepulse inhibition (PPI), a measure of sensorimotor gating (Okaty et al., 2015). Given the extensive VGLUT3+ r2-Pet1 innervation of regions that modulate PPI, such as the prelimbic cortex and hippocampus (Bakshi and Geyer, 1998; Rohleder et al., 2016), and the sensitivity of VGLUT3+ MR neurons to sensory stimulation (Domonkos et al., 2016), we suggest that this effect may have been driven by a loss of neurotransmission from the r2-Pet1Vglut3-high neuron subtype. Thus, we predict r2-Pet1Vglut3-high neurons may regulate sensorimotor gating and modulate theta to affect memory.
Pericellular baskets characterize some r2-Pet1Vglut3-high axonal collaterals
We observed VGLUT3+ r2-Pet1 boutons in pericellular baskets in hippocampus, cortex, and septum, with baskets of the latter region also containing 5-HT+/VGLUT3+ co-positive r2-Pet1 boutons. We propose that these co-positive boutons derive from r2-Pet1Vglut3-high collaterals because our analyses of Piezo2-Pet1 fibers from the r2-Pet1Tph2-high subgroup failed to reveal baskets. This would further suggest that a fraction of r2-Pet1Vglut3-high cells express sufficient, although low, levels of TPH2 or the 5-HT reuptake transporter SLC6A4 to generate 5-HT+ boutons. Pericellular baskets formed by 5-HT+ fibers (Köhler et al., 1982; Freund et al., 1990; Hornung et al., 1990; Dinopoulos et al., 1993; Leger et al., 2001) and VGLUT3+ fibers (Hioki et al., 2004; Riedel et al., 2008; Amilhon et al., 2010) have been reported previously; we suggest that fibers from r2-Pet1Vglut3-high cells may be reflected, in part, in the latter.
Also revealed here are composite baskets, where axons from multiple Pet1 neuron subtypes converge on a single targeted cell in cortex, hippocampus, and septum. In lateral septum composite baskets, r2-Pet1 varicosities were less likely to be 5-HT+, while both r2-Pet1 and other Pet1 varicosities were similarly likely to be VGLUT3+, suggesting that these subtypes of Pet1 neurons may release different neurotransmitters onto the basket-targeted cell. The commonality of VGLUT3 in baskets suggests that it may play a role in basket function, possibly by synergizing cotransmission of 5-HT and glutamate (Amilhon et al., 2010). Given that previous retrograde tracing studies identified inputs from both the DR and MR to the LS (Köhler et al., 1982; Deng et al., 2019); we suggest that the other source of non-r2 Pet1 5-HT+ input to pericellular baskets is from En1-Pet1 neurons arising embryonically from r1 and/or r0 (Jensen et al., 2008). These results suggest that pericellular baskets are complex, potentially releasing either 5-HT, glutamate, or both depending on which set of Pet1 neurons and associated axonal terminals are excited.
In the cortex, we observed r2-Pet1 VGLUT3+ pericellular baskets around subtypes of GABAergic interneurons. A particularly intriguing target was a set of inhibitory neurons that themselves were also VGLUT3+ and which likely correspond to a subset of CCK-expressing interneurons referred to as “basket cells” (Somogyi et al., 2004; Fasano et al., 2017; Harris et al., 2018). Previous work proposed that the source of these presynaptic basket structures encasing these VGLUT3+ basket interneurons was intracortical rather than long-range from the raphe, because the baskets persisted after pharmacological ablation of serotonin neurons with 5,7-DHT (Hioki et al., 2004). It is possible that r2-Pet1Vglut3-high neurons are poorly killed by 5,7-DHT, because of their reduced expression of serotonin transporter SLC6A4, monoamine oxidase A (MAOA), and/or MAOB (Okaty et al., 2015), which may be important for the uptake and toxicity of 5,7-DHT (Finnegan et al., 1989; Paterak and Stefański, 2014). Additionally, there are likely multiple sources for these VGLUT3+ presynaptic baskets, such that the loss of those composed of fibers from the MR goes undetected. Our work then clarifies the existence of a long-range pericellular basket-forming projection from Pet1 MR cells to a specific subgroup of cortical interneurons.
An analogous projection from r2-Pet1 neurons to VGLUT3+ putative basket cell interneurons in the hippocampus was also observed. As in the cortex, many of the postsynaptic cells were immunopositive for 5-HT3A, the ionotropic 5-HT receptor. It has been suggested that an MR glutamatergic input to 5-HT3A-expressing hippocampal interneurons “prepotentiates” them to be responsive to subsequent presentation of 5-HT (Jackson et al., 2009). We suggest that the r2-Pet1Vglut3-high subtype is well positioned for this role, as it targets 5-HT3A-expressing interneurons and, by forming pericellular baskets around these cells, could deliver strong, temporally precise excitation immediately before 5-HT release. This circuit could allow r2-Pet1Vglut3-high neurons to control the activity of basket cell interneurons that then each have highly collateralized axons extending to many principal cells, giving the r2-Pet1 neuron a privileged position for tuning overall cortical and hippocampal network activity.
Conclusion
We describe neurochemically divergent groups of r2-Pet1 neurons, one largely serotonergic and one largely glutamatergic, each arising from the same embryonic progenitor zone, rhombomere 2. We propose that the anatomic segregation of their collective ascending boutons by neurochemical phenotype reflects two separate, largely non-overlapping neuronal subsystems that may each mediate distinct functions of the median raphe. The work described here uncovers a new level of structural organization within the median raphe.
Footnotes
This work was supported by National Institutes of Health Grants F31-NS-108406 (to R.A.S.), T32-HL-007901 (to N.S.), and R21-DA-036056 (to S.M.D. and M.E.F.). The experiments presented in this article were made possible by the use of microscopes in the MicRoN Microscopy Core (Harvard Medical School), use of Imaris software in the Neurobiology Imaging Facility (Harvard Medical School; funded by National Institute of Neurological Disorders and Stroke P30 Core Center Grant NS-072030) and critical conversations with Dr. Paula Montero-Llopis and Ryan Stephansky of MicRoN. We thank Drs. Plummer, Jensen, and colleagues for sharing the RC-FLTG mouse line. We also thank Drs. Benjamin Okaty, Ryan Dosumu-Johnson, and Benjamin Rood for discussions and feedback; and Jia Jia Mai for technical support.
The authors declare no competing financial interests.
- Correspondence should be addressed to Susan M. Dymecki at dymecki{at}genetics.med.harvard.edu