Abstract
Cell-type diversity in the brain enables the assembly of complex neural circuits, whose organization and patterns of activity give rise to brain function. However, the identification of distinct neuronal populations within a given brain region is often complicated by a lack of objective criteria to distinguish one neuronal population from another. In the external segment of the globus pallidus (GPe), neuronal populations have been defined using molecular, anatomical, and electrophysiological criteria, but these classification schemes are often not generalizable across preparations and lack consistency even within the same preparation. Here, we present a novel use of existing transgenic mouse lines, Lim homeobox 6 (Lhx6)–Cre and parvalbumin (PV)–Cre, to define genetically distinct cell populations in the GPe that differ molecularly, anatomically, and electrophysiologically. Lhx6–GPe neurons, which do not express PV, are concentrated in the medial portion of the GPe. They have lower spontaneous firing rates, narrower dynamic ranges, and make stronger projections to the striatum and substantia nigra pars compacta compared with PV–GPe neurons. In contrast, PV–GPe neurons are more concentrated in the lateral portions of the GPe. They have narrower action potentials, deeper afterhyperpolarizations, and make stronger projections to the subthalamic nucleus and parafascicular nucleus of the thalamus. These electrophysiological and anatomical differences suggest that Lhx6–GPe and PV–GPe neurons participate in different circuits with the potential to contribute to different aspects of motor function and dysfunction in disease.
Introduction
The basal ganglia are a set of subcortical nuclei that play an important role in facilitating voluntary movement in health and disease (Albin et al., 1989; DeLong, 1990). The basal ganglia are thought to control movement through the coordinated activity of the direct and indirect pathways that have opposite effects on movement (Alexander and Crutcher, 1990; Smith et al., 1998b; Kravitz et al., 2010). The external segment of the globus pallidus (GPe) is a central nucleus in the motor-suppressing indirect pathway, which receives inputs from the striatum, subthalamic nucleus (STN), and parafascicular nucleus of the thalamus (pf; Robledo and Féger, 1990; Kincaid et al., 1991; Kita and Kitai, 1991; Parent and Hazrati, 1995). In vivo recordings have shown that GPe neurons fire independently (Nini et al., 1995), with firing rates that are characterized by complex temporal patterns and low correlation during movement (DeLong et al., 1985; Wichmann et al., 1994; Jaeger et al., 1995; Bergman et al., 1998).
Overall, the specific contributions of the GPe to basal ganglia function are not well defined, but the GPe has been strongly implicated in the onset and maintenance of motor dysfunction in movement disorders, such as dystonia and Parkinson's disease (Lozano et al., 2000; Di Cristo et al., 2004; Pradeepkumar et al., 2004; Kita, 2007; Obeso et al., 2008). Most notably, the activity of GPe neurons in patients suffering from Parkinson's disease is reduced, and the synchronization of GPe neurons is thought to contribute to pathological oscillations (Nini et al., 1995; Plenz and Kital, 1999; Raz et al., 2000; Bevan et al., 2002; Levy et al., 2002; Terman et al., 2002; Isaksson et al., 2004). Consistent with the hypothesis that oscillations contribute to motor dysfunction, deep brain stimulation of the GPe has been shown to relieve parkinsonian motor symptoms in humans (Vitek et al., 1998, 2004, 2012; Yelnik et al., 2000).
Progress in linking activity of GPe neurons to behavior has been complicated by the fact that the GPe contains a heterogeneous population of neurons that likely contribute differentially to motor function in healthy and diseased individuals. Anatomical studies have shown that GPe neurons project to a number of brain areas, but most of these projections are given little consideration in models of basal ganglia function (Parent et al., 1983; Smith and Bolam, 1989; Kita and Kitai, 1994; Smith et al., 1998a; Sato et al., 2000). For example, the identification of pallidostriatal neurons using tracing techniques in rodents and primates suggests a complexity of neural circuit organization in the GPe that has long been underappreciated (Kita and Kitai, 1994; Bevan et al., 1998; Kita et al., 1999; Sato et al., 2000; Mallet et al., 2012). Attempts have been made to classify cell types in the GPe, but these criteria are often difficult to generalize across preparations (Nambu and Llinaś, 1994; Cooper and Stanford, 2000; Kita, 2007; Deister et al., 2013), rely on postmortem cell-type identification (Hoover and Marshall, 1999; Sato et al., 2000; Hoover and Marshall, 2002), or rely on activity in disease states (Chattopadhyaya et al., 2004; Opalinska et al., 2004; Mallet et al., 2012).
Here, we use transgenic mouse lines to define two genetically distinct, non-overlapping populations of GPe neurons in healthy rodents. Viral expression of enhanced yellow fluorescent protein (EYFP) revealed that Lim homeobox 6 (Lhx6)–GPe and parvalbumin (PV)–GPe neurons differed in their axonal projection patterns as well as electrophysiological properties measured in slice recordings. These results provide a new set of tools to target molecularly, anatomically, and electrophysiologically distinct cell types in the GPe that will enable studies of the organization and function of GPe circuits in health and disease.
Materials and Methods
Mouse lines and tissue processing.
All procedures were performed in accordance with the guidelines from the National Institutes of Health and with approval of the Carnegie Mellon University Institutional Animal Care and Use Committee. Male and female heterozygous 4- to 8-week-old mice on a C57BL/6J background were used for experiments. Lhx6–GPe neurons were targeted using Lhx6–GFP (GENSAT, www.gensat.org) or Lhx6–Cre (Fogarty et al., 2007) mouse lines. PV–GPe neurons were targeted using a Pvalb–2A–Cre mouse line (Madisen et al., 2010).
For immunohistochemistry, mice were perfused transcardially with PBS, followed by 4% paraformaldehyde in PBS. Brains were removed and drop-fixed with the same fixative at 4°C for 24 h. After rising with PBS, brains were transferred to 30% sucrose in PBS and stored at 4°C for at least 24 h before sectioning. Tissue was sectioned in 30 μm sections, blocked with 10% normal donkey serum, and permeabilized with 0.5% Triton X-100 for 1 h. Primary antibody incubations were performed at 4°C for 24 h, using rabbit anti-GFP (1:500; Millipore), chicken anti-GFP (1:1000; Aves), mouse anti-NeuN (1:1000; Millipore), rabbit anti-TH (1:1000), or at 4°C for 48 h when using rabbit anti-PV (1:1000; Swant). Primary antibodies were detected with Alexa Fluor 488-conjugated goat anti-chicken (1:500; Vector Laboratories), Alexa Fluor 488-conjugated goat anti-rabbit (1:500; Vector Laboratories), Alexa Fluor 568-conjugated donkey anti-rabbit (1:500; Vector Laboratories), Alexa Fluor 568-conjugated donkey anti-mouse (1:500; Vector Laboratories), Alexa Fluor 647-conjugated donkey anti-rabbit (1:500; Vector Laboratories), or Alexa Fluor 647-conjugated donkey anti-mouse (1:500; Vector Laboratories), incubated for 2 h at room temperature.
For processing of Fluoro-Gold (FG) labeling, primary antibody incubations were performed at 4°C for 48 h using rabbit anti-FG (1:3000; Fluorochrome). Primary antibody was detected using Alexa Fluor 647-conjugated donkey anti-rabbit (1:500; Vector Laboratories) for 3 h at room temperature.
Animal surgery and viral injections.
Injections of adenoassociated virus 2 (AAV2)–double-floxed inverted open reading frame (DIO)–EYFP or AAV5–DIO–channelrhodopsin 2 (ChR2)–EYFP (used for electrophysiology), produced at the University of North Carolina (Vector Core Facility), were made in 4- to 5-week-old PV–Cre or Lhx6–Cre transgenic mice. Anesthesia was induced using 50 mg/ml ketamine and 15 mg/ml xylazine and maintained throughout surgery using 2% isoflurane. Mice were placed in a stereotaxic frame (David Kopf Instruments), the scalp was opened, and bilateral holes were drilled in the skull (0.27–0.30 mm anterior, 2.1–2.2 mm lateral from bregma). One hundred fifty to 200 nl of virus were injected with a Nanoject (Drummond Scientific) through a pulled glass pipette (∼30 μm tip diameter) whose tip was positioned 3.70 mm below the top of the skull. To prevent backflow of virus, the pipette was left in the brain for 5 min after completion of the injection. Animals were housed for at least 2 weeks after injection to allow time for viral expression.
Injections of FG were done using the same surgical procedures, but bilateral holes were drilled over the striatum (1.0 mm anterior, 2.3 mm lateral from bregma), and 1 μl of FG (0.5% in saline; Fluorochrome) was injected 3.7 mm below the surface of the skull using a syringe pump (0.2 μl/min; Genietouch). Mice were perfused 2 weeks after injections, and tissue was prepared as described above.
Fluorescence quantification.
Quantification of axonal projection patterns of Lhx6–GPe and PV–GPe neurons was performed in sagittal sections from hemispheres in which GPe injections met the following criteria: viral expression within the GPe was strongest within the central plane of analysis (1.56 mm lateral; Paxinos, 2004, his Fig. 114), and labeled cell bodies were observed throughout the full rostrocaudal, dorsoventral, and mediolateral extent of the structure. Analyses were performed in seven sagittal planes (lateral from bregma; in mm): 3.12, 2.76, 2.04, 1.56, 1.20, 0.84, and 0.48 (Paxinos, 2004). Hemispheres in which there was significant spillover of virus and EYFP expression in neighboring structures, such as the reticular thalamic nucleus (RT), ventral palladium, and striatum, were excluded from analysis. To aid in visualization of axonal projections from the GPe, EYFP signal was enhanced with an immunostain against GFP. Epifluorescent images were taken at 25× magnification in seven sagittal slices of each brain, and the same sectional planes across animals were used for analyses. To control for differences in viral expression or antibody retrieval, data from each animal were normalized to the fluorescence intensity in the GPe during analysis. Images were analyzed offline using the pixel intensity-measuring tool in NIH ImageJ. Data were quantified by averaging the pixel intensity within and across each brain structure. On average, three 100 × 100 pixels (75 × 75 μm; area is 5625 μm2) squares were quantified in each of the images. The mean of all pixel intensities designated as within a single structure was found and normalized to GPe expression of each animal. For example, dorsolateral striatum was identified in three planes of section (lateral from bregma; in mm): 3.12, 2.76, and 2.04. Three measures of pixel intensity were taken by randomly placing the 75 × 75 μm square throughout the image. The mean of each image was taken and then averaged across the three planes of section to identify the average pixel intensity of dorsolateral striatum. Finally, the average was normalized to the brightest GPe expression found within the animal, which was located in the central plane of analysis (1.56 mm lateral; Paxinos, 2004, his Fig. 114). Striatum, RT, internal segment of the globus pallidus (GPi), basolateral amygdala (BLA), parasubthalamic nucleus (PSTh), and pf were analyzed using the method described above.
Because of the irregular shape of the STN and the ease at which the STN borders could be identified at 25× magnification, the NIH ImageJ hand draw tool was used to measure average pixel intensity of the entire structure. Similar to other structures, the STN was identified in more than one plane of section, and the pixel intensity from each plane was averaged and then normalized to GPe expression.
In the case of the substantia nigra (SN), a tyrosine hydroxylase (TH) stain was used to define the borders between dopamine neurons of the pars compacta (SNc) and the nondopaminergic (GABAergic) neurons of the pars reticulata (SNr). GFP and TH expression was sequentially imaged in the GPe and SN using the epifluorescent microscope at 10× magnification. For NIH ImageJ analysis, zones were clearly defined using TH expression as a marker for the SNc and using the hand draw tool described above. These zones were then overlaid on the images showing GFP expression, and the average pixel intensity was found. The average across sections was taken and normalized to GPe expression. Differences in normalized fluorescence intensities were compared using a t test.
Results of FG retrograde labeling were analyzed in two to four sections each from three PV mice (PV–Cre × Ai3 or PV–Cre × Ai9) and three Lhx6 mice (Lhx6–GFP). GPe sections were magnified to 63× on an epifluorescent microscope. The number of FG-positive (FG+) neurons in an 225 × 175 μm square of tissue were counted and scored as Lhx6+, PV+, or neither. Populations of FG+ neurons that were either Lhx6 or PV were compared using a Fisher's exact test.
Three-dimensional Sholl analysis.
To quantify the selective innervation of striatal fast-spiking striatal interneurons (FSIs) by Lhx6–GPe and PV–GPe neurons, a three-dimensional Sholl analysis was performed as described previously (Chattopadhyaya et al., 2004). Briefly, tissue from mice expressing EYFP in Lhx6–GPe or PV–GPe neurons was immunostained for PV and NeuN to identify FSIs and medium spiny neurons (MSNs) in the striatum, respectively. Although NeuN does not selectively stain MSNs, MSNs make up 95% of neurons in the striatum, so for the purposes of this analysis, we assumed that NeuN+/PV-negative (PV−) neurons were most likely MSNs. Confocal z-stack images of presumed striatal FSIs and MSNs were taken at 63× magnification, and channels were sequentially imaged to avoid bleedthrough. Stacked images were merged using NIH ImageJ, and concentric circles of increasing diameter (2 μm) were drawn from the center of each identified neuron. Innervation was quantified as the number of axonal varicosities counted in each circle.
Electrophysiological recordings.
Parasagittal sections (300 μm thickness) containing the GPe or pf were prepared from brains of 5- to 8-week-old mice. Lhx6–GPe neurons were targeted in slices from Lhx6–GFP mice, Lhx6–Cre mice injected in the GPe with AAV2–DIO–EYFP, or Lhx6–Cre mice crossed to the reporter strains Ai9 (red) or Ai3 (green) (Jackson Laboratories). PV–GPe neurons were targeted in slices from PV–Cre mice injected in the GPe with AAV2–DIO–EYFP or PV–Cre crossed to Ai9 or Ai3 reporter strains.
Slices were prepared on a Leica VT1200 vibratome in an ice-cold HEPES cutting solution (in mm): 92 NaCl, 2.5 KCl, 1.2 NaH2PO4, 30 NaHCO3, 20 HEPES, 25 glucose, 5 sodium ascorbate, 2 thiourea, 3 sodium pyruvate, 10 MgCl2, and 0.5 CaCl2. Slices were allowed to recover for 15 min at 33°C in a chamber filled with N-methyl-d-glucamine (NMDG)–HEPES recovery solution (in mm): 93 NMDG, 2.5 KCl, 1.2 NaH2PO4, 30 NaHCO3, 20 HEPES, 25 glucose, 10 MgSO4, 0.5 CaCl2, 5 sodium ascorbate, 2 thiourea, and 3 sodium pyruvate. After 15 min, slices were held at room temperature for at least 1 h before recording in a holding solution that was similar to the HEPES cutting solution but with 1 mm MgCl2 and 2 mm CaCl2. Recordings were made at 33°C in carbogenated ACSF (in mm): 125 NaCl, 26 NaHCO3, 1.25 NaH2PO4, 2.5 KCl, 12.5 glucose, 1 MgCl2, and 2 CaCl2. For recordings of intrinsic excitability, 50 μm picrotoxin and 5 μm DNQX were included in the bath to block fast inhibitory and excitatory transmission, respectively.
Data were collected with a MultiClamp 700B amplifier (Molecular Devices) and ITC-18 analog-to-digital board (HEKA) using Igor Pro software (Wavemetrics) and custom acquisition routines (Recording Artist; Richard C. Gerkin, Phoenix, AZ). Current-clamp recordings were filtered at 10 kHz and digitized at 40 kHz; voltage-clamp recordings were filtered at 2 kHz and digitized at 10 kHz. Electrodes were made from borosilicate glass (pipette resistance, 2–4 m). The internal solution for current-clamp recordings consisted of the following (in mm): 130 KeMSO3, 10 NaCl, 2 MgCl2, 0.16 CaCl2, 0.5 EGTA, 10 HEPES, 2 Mg-ATP, and 0.3 Na-GTP. The internal solution for voltage-clamp recordings consisted of the following (in mm): 120 CsMeSO3, 15 CsCl, 8 NaCl, 0.5 EGTA, 10 HEPES, 2 Mg-ATP, 0.3 Na-GTP, and 5 QX-314.
Differences in electrophysiological parameters were compared using a Wilcoxon's nonparametric test. Action potential (AP) half-width was calculated as the width of the AP, measured at a point halfway between threshold and AP peak. AP threshold was defined as the voltage at which the acceleration in voltage exceeded 3 × 105 V/s. Maximum firing rate was measured as the average maximum firing rate over a 1 s step that could be sustained without entering depolarization block. Input resistance was calculated in voltage clamp as the deviation from the holding current caused by a 5 mV hyperpolarizing step. The change in current was calculated as the difference between baseline holding current and the average current during the last 10 ms of the 40 ms hyperpolarizing voltage step. Whole-cell capacitance was calculated by integrating the area under the transient after a 5 mV hyperpolarizing voltage step from the holding potential, −80 mV. The holding current was subtracted before integrating.
Results
Distinct neuronal subpopulations identified in Lhx6 and PV transgenic mouse lines
Transgenic mouse lines were used to identify two distinct subpopulations of neurons in the GPe. In the first transgenic line, Pvalb–2a–Cre (PV–Cre; Madisen et al., 2010), gene expression was driven under the promoter for PV, a calcium binding protein found in GABAergic neurons in the CNS. In the second transgenic line, Lhx6–Cre (Fogarty et al., 2007) or Lhx6–GFP (GENSAT), gene expression was driven under the promotor for Lhx6, a homeobox protein found in a subset of GABAergic neurons in the brain (Marin et al., 2000; Cobos et al., 2006). An immunostain against PV in tissue from an Lhx6–GFP mouse revealed non-overlapping populations of neurons in the GPe (Fig. 1A). Lhx6–GPe and PV–GPe neurons each accounted for approximately one-third of NeuN+ neurons in the GPe. Lhx6–GPe neurons made up 34% (1030 of 3517) of NeuN+ GPe neurons, whereas PV–GPe neurons made up 29% (1149 of 3517), with only 2% (54 of 3517) of NeuN+ GPe neurons expressing both (Fig. 1B). Neurons were scored as positive for Lhx6 or PV if pixel intensity was >10× background.
Because the fraction of PV–GPe neurons observed with immunostaining was lower than what has been reported in other species, we used a second genetic approach to quantify the fraction of PV–GPe neurons. PV–Cre mice were crossed to fluorescent reporter mice (Ai3), and fixed sections from these animals were immunostained with NeuN and PV. Consistent with our immunohistochemistry data, genetically identified PV–GPe neurons were less than half of all GPe neurons (41%, 235 of 571), and 78% (1226 of 1544) of genetically defined PV–GPe neurons were immunopositive for PV.
To quantify the spatial distribution of Lhx6 and PV neurons within the GPe, cells were counted along mediolateral, dorsoventral, and rostrocaudal axes. No gradients were observed along the rostrocaudal and dorsoventral axes, but opposing gradients were observed along the mediolateral axis (Fig. 1C–E). Lhx6–GPe neurons were most concentrated in the medial portions of the GPe and decreased in concentration from medial to lateral, whereas PV–GPe neurons were most concentrated in the lateral and less concentrated in the more medial portions of the GPe.
Lhx6–GPe and PV–GPe neurons express different baseline firing rates
To determine whether Lhx6–GPe and PV–GPe neurons correspond to electrophysiologically distinct subpopulations, neurons of each cell type were targeted for recording in acute slices from 5- to 8-week-old Lhx6 or PV transgenic mice. To determine the rates of intrinsic pacemaking in each cell type, 50 μm picrotoxin and 5 μm DNQX were bath applied to block fast inhibitory and excitatory transmission, respectively. In the cell-attached configuration, both cell types fired tonically with low interspike variability [coefficient of variation (CV)Lhx6 = 0.31 ± 0.39, n = 25 vs CVPV = 0. 22 ± 0.25, n = 20; p = 0.1; Fig. 2A]. On average, PV–GPe neurons exhibited faster extracellular firing rates than Lhx6–GPe neurons (Lhx6, 34.2 ± 33.2 Hz, range of 1.3–142 Hz vs PV, 64.4 ± 34.3 Hz, range of 9.2–131 Hz; p = 0.001; Fig. 2B). To ensure that extracellular firing rates were not influenced by the voltage-clamp recording configuration or the gigaohm seal, we recorded extracellular firing rates in a second population of neurons using loose patch in cell-attached recording mode. Once again, PV–GPe neurons were found to have significantly higher extracellular firing rates than Lhx6–GPe neurons (Lhx6, 28.9 ± 17.7 Hz, range of 3.9–55.2 Hz, n = 20 vs PV, 66.1 ± 39.0 Hz, range of 14.8–140 Hz, n = 17; p < 0.0001). In the whole cell-recording configuration, spontaneous firing persisted in 20 of 27 PV–GPe neurons but in only 18 of 38 Lhx6–GPe neurons (Fig. 2C). Spontaneous firing rates were measured within 5 min of break-in because spontaneous firing rates tended to rundown over the duration of a recording. Spontaneous firing rates recorded in whole-cell configuration were also significantly faster in PV–GPe neurons compared with Lhx6–GPe neurons (Lhx6, 7.8 ± 14.2 Hz, n = 38 vs PV, 28.4 ± 24.9 Hz, n = 27; p = 0.0003; Fig. 2D).
Because autonomous pacemaking in GPe neurons relies in part on hyperpolarization-activated, cyclic nucleotide-gated cation (HCN) channels (Chan et al., 2004), we hypothesized that differences in intrinsic firing rates between Lhx6–GPe and PV–GPe neurons might arise from differences in hyperpolarization-activated currents (Ih). To measure Ih, Lhx6–GPe and PV–GPe neurons were given 3 s hyperpolarizing current injections of increasing magnitude (Fig. 2E). In the current-clamp recording configuration, hyperpolarization resulted in an Ih-mediated rectification of the membrane potential, often referred to as a membrane potential “sag,” which increased monotonically as a function of membrane potential (Fig. 2F). Surprisingly, we observed no difference in the maximum sag recorded in Lhx6–GPe and PV–GPe neurons (Lhx6, 16.1 ± 7.0 mV, n = 35 vs PV, 13.8 ± 5.5 mV, n = 22; p = 0.22; Fig. 2G) and no difference in the slope of the monotonic relationship between sag and membrane potential (Lhx6, −0.264 ± 0.14 vs PV, −0.259 ± 0.13; p = 0.77; Fig. 2H). These results suggest that differences in Ih expression are unlikely to account for differences in intrinsic firing rates of Lhx6–GPe and PV–GPe neurons but do not rule out this possibility entirely.
Lhx6–GPe and PV–GPe neurons have different passive and active membrane properties
Variability in a number of passive and active membrane properties has been reported across GPe neurons, including resting membrane potential, firing rate adaptation, and AP shape (Nambu and Llinaś, 1994; Cooper and Stanford, 2000). To determine whether passive membrane properties vary systematically between Lhx6–GPe and PV–GPe neurons, neurons were held in voltage clamp (−80 mV) and given brief hyperpolarizing steps (−5 mV, 100 ms). Figure 3A shows the response of a representative Lhx6–GPe neuron and a representative PV–GPe neuron. Input resistance, calculated from the current flowing during the hyperpolarizing step, was significantly greater in Lhx6–GPe relative to PV–GPe neurons (3000.6 ± 168 MΩ, n = 38 vs 170 ± 88 MΩ, n = 28; p = 0.0007; Fig. 3B). Whole-cell capacitance, calculated as the integral under the transient (Fig. 3A, inset), was significantly greater in PV–GPe relative to Lhx6–GPe neurons (241 ± 113 pF, n = 28 vs 164 ± 65 pF, n = 38; p = 0.002; Fig. 3C). Because capacitance is proportional to the surface area of a neuron, these results suggest that PV–GPe neurons are larger or have more extensive dendritic arborizations than Lhx6–GPe neurons.
To compare active membrane properties of Lhx6–GPe and PV–GPe neurons, the recording configuration was switched to current clamp, and neurons were driven to fire with 1 s depolarizing current steps. The amplitude of injected current was increased incrementally until a neuron reached its maximum firing rate, defined as the average firing rate sustained over the 1 s step without entering depolarization block (Fig. 3D). The maximum firing rate of PV–GPe neurons was significantly greater than that of Lhx6–GPe neurons (347 ± 80 Hz, n = 38 vs 248 ± 80 Hz, n = 27; p < 0.0001). Maximum firing rate was well correlated with AP width (Fig. 3F), and PV–GPe neurons had narrower APs than Lhx6–GPe neurons (0.28 ± 0.09 ms, n = 27 vs 0.37 ± 0.11 ms, n = 38; p = 0.003; Fig. 3E). Although AP widths and maximum firing rates varied continuously across the population of GPe neurons, Lhx6–GPe neurons clustered more at one end of the distribution, whereas PV–GPe neurons clustered more at the other end of the distribution (Fig. 3F). Using maximum firing rate as the sole criteria enabled Lhx6–GPe and PV–GPe neurons to be distinguished with ∼75% accuracy (78% of PV–GPe neurons had maximum firing rates >300 Hz, and 72% of Lhx6–GPe neurons had maximum firing rates <300 Hz). This degree of overlap emphasizes the importance of using genetic markers to segregate cell types in the GPe. A list of additional parameters quantifying the firing properties of Lhx6–GPe and PV–GPe neurons is provided in Table 1.
Projections of Lhx6–GPe and PV–GPe neurons to the STN
To study the axonal projection patterns of Lhx6–GPe and PV–GPe neurons, we used the Cre-dependent AAV2 to express EYFP selectively in one neuronal population or the other. Virus was stereotaxically injected into the GPe of PV–Cre and Lhx6–Cre transgenic mice, and tissue was processed and analyzed 2 weeks later. To control for variations in viral expression or immunolabeling across mice, fluorescence intensities in target brain regions were normalized to fluorescence intensity in the GPe of the same animal (see Materials and Methods). Fluorescence intensities measured in the GPe ranged from 146 to 249 a.u. but were not significantly different between Lhx6 and PV mice (Lhx6, 216 ± 27 a.u., n = 4 vs PV, 202 ± 51 a.u., n = 4; p = 0.6; Fig. 4C). To confirm that fluorescence intensity varied similarly as a function of axonal density in both mouse lines, we measured the fluorescence intensity of a small region of tissue (143 μm2) and manually counted the number of axonal varicosities present. In both mouse lines, a linear relationship was observed between fluorescence intensity and number of boutons, quantified in brain areas receiving low (striatum), medium (SNr), and high (STN) amounts of axonal projections (Fig. 4D). Therefore, we were able to use the normalized fluorescence intensity to compare the density and pattern of axonal projections from Lhx6–GPe and PV–GPe neurons to many different brain regions.
The first brain region investigated was the STN, the canonical target of the GPe (Fig. 4A). Both Lhx6–GPe and PV–GPe neurons densely innervated the STN, but innervation from PV–GPe neurons was significantly greater than innervation from Lhx6–GPe neurons (Lhx6, 59.9 ± 13.6 a.u., n = 4 vs PV, 96.2 ± 11.3 a.u., n = 4; p = 0.006; Fig. 4B). This difference in innervation density was primarily attributable to the unique patterning of Lhx6 projections to the STN. Unlike PV–GPe projections, which were evenly distributed throughout the extent of the nucleus, Lhx6–GPe projections were clustered around the periphery of the STN and did not target the central portion of the nucleus (Fig. 4E).
To examine the topographical organization of projections from the GPe to the STN, viral injections were targeted to rostral and caudal portions of the GPe (Fig. 4F). No differences in axonal projection patterns to the STN were observed in the Lhx6–Cre line from different regions of the GPe. In contrast, axonal projection patterns to the STN were topographically organized in the PV–Cre line. The rostroventral region projected to the central portion of the STN and rarely showed expression in the periphery (Fig. 4G, left). In contrast, the caudal region of the GPe projected to the periphery of the STN, similar to what was observed in the Lhx6 population (Fig. 4G, right). These results suggest a topographic organization of PV projections to distinct zones within the STN.
Lhx6–GPe and PV–GPe neurons differentially innervate the SNc but similarly innervate GABAergic output nuclei of the basal ganglia
The SN is subdivided into two regions that subserve distinct purposes in basal ganglia function (Fig. 5A). The SNr is made up of GABAergic neurons that act as a major output of the basal ganglia (Samantaray et al., 2004). The SNc is a region that lies dorsal to the SNr and contains a dense collection of dopaminergic neurons that innervate the striatum and many other structures (Björklund and Dunnett, 2007). GPe projections have been shown to innervate both the SNr and SNc in rats, but these connections are still poorly understood (Paladini et al., 1999). In the SNr, there was no significant difference in density of axons arising from Lhx6–GPe versus PV–GPe neurons (Lhx6, 15.3 ± 3.85 a.u., n = 4 vs PV, 22.8 ± 14.7 a.u., n = 3; p = 0.3; Fig. 5B). Both populations of neurons projected to the SNr diffusely with no clear topographical organization (Fig. 5C, top and middle).
A markedly different scenario was observed in the SNc. As seen in Figure 5C, the distinct clustering of axons was observed in the dorsal region of the SN. To resolve the borders of the SNc and SNr, tissue was stained with antibodies against TH to label dopamine neurons. Axonal projections to the SNc were more concentrated than those to the SNr, and the majority of SNc projections arose from Lhx6–GPe neurons (Lhx6, 41.0 ± 9.40 a.u., n = 4 vs PV, 21.4 ± 6.05 a.u., n = 3; p = 0.03; Fig. 5C).
The GPi, one of the major output structures of the basal ganglia, is nested in the internal capsule fibers between the GPe and STN (Fig. 5D). Both Lhx6–GPe and PV–GPe neurons projected to the GPi in equal amounts (Lhx6, 37.5 ± 25.7 a.u., n = 4 vs PV, 45.8 ± 11.1 a.u., n = 4; p = 0.57; Fig. 5E). As expected, axonal projections from the both populations of GPe neurons formed basket synapses onto the neurons of the GPi (Fig. 5F). Viral expression was not observed in GPi neurons.
Lhx6–GPe neurons project more densely to the dorsolateral striatum than PV–GPe neurons
The densities of axons from Lhx6–GPe or PV–GPe neurons innervating the striatum were quantified within three functionally distinct striatal zones: dorsolateral, dorsomedial, and ventral. The dorsolateral region of the striatum is important for sensorimotor processing, whereas ventral striatum deals with appetitive and rewarding aspects of behavior (Packard and Knowlton, 2002; Kelley, 2004; O'Doherty et al., 2004). The dorsomedial region is important for instrumental conditioning and behavioral flexibility and is considered the associative territory of striatum (Ragozzino et al., 2002; Yin et al., 2005). GPe neurons sent sparse projections to striatum and formed a gradient of high to low density along the dorsoventral axis and the lateromedial axis (Fig. 6B,C).
In the dorsolateral region of the striatum (DL Str), in which GPe projections were most dense, axonal projections from Lhx6–GPe neurons were significantly more dense than axonal projections from PV–GPe neurons (DL StrLhx6, 3.80 ± 1.38 a.u., n = 4 vs DL StrPV, 1.14 ± 0.27 a.u., n = 4; p = 0.009; Fig. 6B,C). In contrast, axonal projections to the medial striatum (DM Str), an associative territory of the striatum, were less dense for both cell types, and the density did not differ significantly between Lhx6–GPe and PV–GPe neurons (DM StrLhx6, 0.79 ± 1.32 a.u., n = 4 vs DM StrPV, 0.32 ± 0.29 a.u., n = 4; p = 0.5; Fig. 6B). Similarly, the ventral striatum (V Str) showed little to no axonal projections from both populations of neurons (V StrLhx6, 0.34 ± 0.23 a.u., n = 4 vs V StrPV, 0.19 ± 0.31 a.u., n = 4; p = 0.48).
Consistent with previous work, axonal projections from both Lhx6–GPe and PV–GPe neurons appeared to target GABAergic interneurons in the striatum (Bevan et al., 1998), as indicated by axonal varicosities that were clustered around PV+ interneurons, corresponding to FSIs (Fig. 6D). To quantify the selectivity of GPe innervation for striatal FSIs compared with MSNs, we counted the number of boutons surrounding each cell type using a Sholl analysis (Chattopadhyaya et al., 2004). Striatal FSIs received significantly more innervation from the GPe than MSNs (Fig. 6E,F).
Our quantification of axon densities suggests that Lhx6–GPe neurons give rise to a greater proportion of striatal projections than do PV–GPe neurons. To confirm this result, we injected the retrograde tracer FG into the dorsolateral striatum of Lhx6–GFP or PV–Cre × Ai9 mice (Fig. 6G,H). After 2 weeks, tissue was processed and analyzed to determine what percentage of retrogradely labeled neurons in the GPe were either Lhx6 or PV. As shown in Figure 6G, FG labeled significantly more Lhx6–GPe neurons than PV–GPe neurons (Lhx6, 47.5 ± 18.4%, n = 84 of 177 vs PV, 30.5 ± 11.0%, n = 59/182; p = 0.003, Fisher's exact test). Together, these data suggest that a greater proportion of GPe projections to the striatum arise from Lhx6–GPe neurons, and the focus of these projections is interneurons in the sensorimotor input region of the striatum.
Distinct projections of PV–GPe neurons to the pf
GPe projections to the thalamus have been described previously (Hazrati and Parent, 1991; Chattopadhyaya and Pal, 2004) but are typically restricted to the RT, which is located rostral to the GPe and extends farther medially in the brain (Fig. 7A). The RT is a specialized thalamic nucleus containing mostly GABAergic neurons that inhibit numerous thalamic nuclei and have no direct cortical projections (Houser et al., 1980; Ohara, 1988; Pinault and Deschênes, 1998). Consistent with previous studies, axonal projections from both populations of neurons were found in the RT with no significant difference in density (Lhx6, 28.0 ± 6.71 a.u., n = 4 vs PV, 24.4 ± 14.0 a.u., n = 4; p = 0.6; Fig. 7B,C).
In contrast to projections to the RT, Lhx6–GPe and PV–GPe neurons differed substantially in their projections to the pf. The Lhx6–GPe population exhibited virtually no axonal projections to the pf, whereas PV–GPe neurons projected densely, covering the full extent of the nucleus (Lhx6, 0.83 ± 1.3 a.u., n = 4 vs PV, 69.0 ± 19.9 a.u., n = 4; p = 0.0005; Fig. 7D,E). As seen in Figure 7E, the nuclei surrounding the pf showed no axonal projections, which would be expected if axons arose from spillover labeling of neurons in the RT. Furthermore, there was a topographic organization of projections from the GPe to the pf. On targeted injections of rostral and caudal portions of the GPe, we found that the pf was only targeted by caudal PV–GPe neurons but not by rostral PV–GPe neurons (data not shown). Injections in the rostral portion of the GPe show little to no expression within the pf. To ensure that PV–GPe axonal projections to the pf were functionally active and not just fibers of passage, a virus containing ChR2–EYFP was injected into the GPe of PV–Cre mice. Optically evoked IPSCs from a holding potential of −70 mV were observed in five of five pf neurons (565.04 ± 312.3 pA, n = 5; Fig. 7F), confirming that PV–GPe neurons form functional inhibitory synapses onto pf neurons. These results suggest the presence of a previously undescribed output projection from the GPe to the pf.
Discussion
Here, we establish the use of Lhx6 and PV transgenic mouse lines to identify genetically distinct classes of GPe neurons that differ molecularly, anatomically, and electrophysiologically. Viral-mediated expression of EYFP revealed that both Lhx6–GPe and PV–GPe neurons differed in their axonal innervation of a number of brain regions, including the striatum, SNc, STN, and the pf. These results reveal cell-type-specific segregation of GPe output and establish genetic tools for elucidating the role of distinct GPe output pathways in behavior and disease.
Lhx6 and PV transgenic mouse lines subdivide GPe neurons into two groups
Approximately two-thirds of GPe neurons were labeled in Lhx6 (one-third) or PV (one-third) transgenic mouse lines. The remaining one-third of GPe neurons could represent a third, as yet unidentified, neuronal population in the GPe or could represent incomplete labeling of Lhx6 and PV populations. Although the fraction of PV–GPe neurons in mice (one-third) was lower than that reported in rats (two-thirds; Kita, 2007), we observed a similarly low number of PV neurons using both immunohistochemical and genetic techniques. These results suggest that the number of PV–GPe neurons is lower in mice than in other species, such as rats and primates.
Lhx6 is a homeobox protein expressed in the medial ganglionic eminence (MGE) during development. The MGE is the origin of most forebrain GABAergic interneurons (Marin et al., 2000; Cobos et al., 2006), as well as GPe neurons (Flandin et al., 2010; Nóbrega-Pereira et al., 2010). It remains unclear why only a subset of GPe neurons are labeled in Lhx6 transgenic mice, but one possibility is that Lhx6–GPe neurons correspond to a fraction of MGE neurons in which Nkx2.1 expression is downregulated during embryonic stages (Nóbrega-Pereira et al., 2010). Whereas PV has been used previously to identify a subset of GPe neurons (Kita, 1994, 2007; Mallet et al., 2012), our results show for the first time that Lhx6 can be used to identify a second, non-overlapping subset of GPe neurons.
Lhx6–GPe and PV–GPe neurons were differentially distributed along the lateromedial axis of the GPe. PV–GPe neurons were more densely distributed in the lateral portion of the GPe, the sensorimotor territory of the nucleus (Alexander and Crutcher, 1990; Matsumura et al., 1995; Haber et al., 2000; François et al., 2004), in which pharmacological activation produces dyskinetic movements in primates (Matsumura et al., 1995; Grabli et al., 2004). In contrast, Lhx6–GPe neurons were more densely distributed in the medial portion of the GPe, the limbic and associative territories of the nucleus (Alexander and Crutcher, 1990; Haber et al., 2000; François et al., 2004), in which pharmacological activation produces stereotypies and hyperactivity (Grabli et al., 2004).
Lhx6 and PV neurons have different intrinsic electrophysiological properties
Differences in the electrophysiological properties of GPe neurons have been described both in vivo (Kelland et al., 1995; Mallet et al., 2008, 2012; Benhamou et al., 2012) and in slice preparations (Nambu and Llinaś, 1994; Cooper and Stanford, 2000; Kita, 2007; Deister et al., 2013), but a rigorous classification scheme has not been established. Furthermore, there is some controversy about whether or not intrinsic properties can even be used to define groups of GPe neurons because these properties vary continuously across the population and can change during the course of an experiment (Günay et al., 2008; Deister et al., 2013).
Because Lhx6–GPe and PV–GPe neurons are genetically defined, this establishes a rigorous classification scheme with which to determine how intrinsic firing properties vary across cell types in the GPe. Intracellular recordings from Lhx6–GPe and PV–GPe neurons revealed that intrinsic parameters varied continuously across the population, but a number of parameters differed significantly between these two cell types. PV–GPe neurons had narrower APs and could sustain higher maximum firing rates than Lhx6–GPe neurons, suggesting that PV–GPe neurons are better tuned to reliably respond to rapid, high-frequency inputs than Lhx6–GPe neurons. In addition, PV–GPe neurons had lower input resistances and higher whole-cell capacitances than Lhx6–GPe neurons, suggesting that PV–GPe neurons are larger.
Lhx6–GPe and PV–GPe neurons also differed with respect to spontaneous firing rates, which were higher in PV–GPe neurons. Spontaneous firing in GPe neurons is mediated by a combination of Ih and persistent Na+ currents (Chan et al., 2004, 2011; Mercer et al., 2007). Our results that sag currents were similar between Lhx6–GPe and PV–GPe neurons suggests that additional currents besides Ih may contribute to the different firing rates between these two cell populations, including persistent Na+ currents or Ca2+-sensitive K+ currents, which have also been shown to regulate firing rates in GPe neurons (Deister et al., 2009; Schultheiss et al., 2010).
Implications of genetically defined output pathways for basal ganglia function
Historically, one of the most robust classification schemes for GPe neurons has been based on their axonal projection patterns (Bevan et al., 1998; Hoover and Marshall, 1999; Kita et al., 1999; Sato et al., 2000; Hoover and Marshall, 2002; Kita, 2007; Mallet et al., 2012). Although both Lhx6–GPe and PV–GPe neurons projected to many of the same brain regions, the density of these projections differed significantly in a number of areas, including the striatum and SNc, which were more densely innervated by Lhx6–GPe neurons, and the STN and pf, which were more densely innervated by PV–GPe neurons (Fig. 8). Both cell types made similarly dense projections to a number of additional brain regions, both within and outside the basal ganglia, including the GPi, SNr, RT, BLA, and PSTh (Fig. 8). We also observed that both cell types made projections to the cortex, but the diffuse nature of these projections made quantification difficult.
Based on axonal density and retrograde labeling, we found that Lhx6–GPe neurons made up a greater portion of GPe projections to the striatum than PV–GPe neurons. This is consistent with previous studies in which most striatal-projecting neurons lack PV expression (Hoover and Marshall, 1999, 2002; Kita et al., 1999; Mallet et al., 2012). Our results suggests that Lhx6 labels a greater proportion of “pallidostriatal neurons” described previously in monkeys and rats (Bevan et al., 1998; Sato et al., 2000). Axons of pallidostriatal neurons target GABAergic interneurons in the striatum (Bevan et al., 1998), and this projection may play an important and underappreciated role in regulating striatal activity (Wilson, 2009; Gage et al., 2010).
Lhx6–GPe neurons also make stronger projections to the SNc than PV–GPe neurons. The SNc contains the majority of dopamine neurons that innervate the dorsal striatum as well as the GPe (Parent and Smith, 1987; Lavoie et al., 1989; Rommelfanger and Wichmann, 2010). Dopamine release from the SNc is critically involved in motor aspects of basal ganglia function (Gerfen et al., 1990; Björklund and Dunnett, 2007; Kreitzer and Malenka, 2008; Gerfen and Surmeier, 2011; Surmeier et al., 2011); therefore, projections from Lhx6–GPe neurons may play an underappreciated yet influential role in movement (Häusser and Yung, 1994; Paladini et al., 1999).
The axonal projections of PV–GPe neurons were significantly stronger to the STN than those of Lhx6–GPe neurons. The difference in STN innervation by PV–GPe and Lhx6–GPe neurons was most apparent in an inner region of the nucleus, essentially devoid of innervation from Lhx6–GPe neurons but strongly innervated by PV–GPe neurons located in the rostral portion of the nucleus. This result confirms the presence of a topographical projection from the GPe to the STN that has been observed in primates (François et al., 2004), but the functional role of the medial portions of the STN in rodents is not known. The GPe and STN form a reciprocally connected network that has received much attention for its role in pathological oscillations in Parkinson's disease (Bergman et al., 1994; Nini et al., 1995; Plenz and Kital, 1999; Magill et al., 2001; Bevan et al., 2002; Terman et al., 2002; Mallet et al., 2008; Holgado et al., 2010; Fan et al., 2012). Although both Lhx6–GPe and PV–GPe neurons project strongly to the STN, changes in the denser PV–GPe projections may play a more prominent role in the development or maintenance of pathological oscillations in Parkinson's disease.
The most striking difference in the projection patterns of Lhx6–GPe and PV–GPe neurons observed was in the pf, which received dense innervation from PV–GPe neurons but not Lhx6–GPe neurons. The pf plays a role in attention and behavioral flexibility, and its projections to the striatum are particularly important for the control of cholinergic interneurons (Matsumoto et al., 2001; Minamimoto and Kimura, 2002; Smith et al., 2004, 2011). The dense projections of PV–GPe neurons to the pf suggest that they are an important component of non-motor output pathways of the GPe. The pf may play a role in motor control as well, based on results that deep brain stimulation in the pf is particularly effective at alleviating tremors in patients with movement disorders (Krauss et al., 2002; Goff et al., 2009; Stefani et al., 2009).
In summary, we have shown that Lhx6 and PV transgenic mice can be used as tools to identify genetically distinct circuits within the GPe with the potential to subserve different roles in basal ganglia function. By virtue of their projections to striatal interneurons and SNc dopamine neurons, Lhx6–GPe neurons are well positioned to regulate striatal output and influence the balance of direct and indirect pathway activity to determine overall levels of motor output. In contrast, PV–GPe neurons make stronger projections to the STN and pf, brain areas with both motor and non-motor functions and areas that have been implicated specifically in the generation of pathological oscillations and tremor in Parkinson's disease, suggesting a role of PV–GPe neurons in these aspects of disease pathophysiology.
Footnotes
This work was supported by National Institutes of Health Grant R00 NS076524. We thank Richard C. Gerkin for specialized IGOR analysis and acquisition routines, Nicoletta Kessaris for providing Lhx6–Cre mice, Hongkui Zeng for providing PV–Cre mice, and Haibeng Teng of the Molecular Biosensor and Imaging Center for assistance with the confocal microscopy.
- Correspondence should be addressed to Aryn H. Gittis, Mellon Institute, Carnegie Mellon University, 4400 Fifth Avenue, Pittsburgh, PA 15213. agittis{at}cmu.edu