Abstract
The debilitating psychomotor symptoms of Huntington's disease (HD) are linked partly to degeneration of the basal ganglia indirect pathway. At early symptomatic stages, before major cell loss, indirect pathway neurons exhibit numerous cellular and synaptic changes in HD and its models. However, the impact of these alterations on circuit activity remains poorly understood. To address this gap, optogenetic- and reporter-guided electrophysiological interrogation was used in early symptomatic male and female Q175 HD mice. D2 dopamine receptor-expressing striatal projection neurons (D2-SPNs) were hypoactive during synchronous cortical slow-wave activity, consistent with known reductions in dendritic excitability and cortical input strength. Downstream prototypic parvalbumin-expressing external globus pallidus (PV+ GPe) neurons discharged at 2-3 times their normal rate, even during periods of D2-SPN inactivity, arguing that defective striatopallidal inhibition was not the only cause of their hyperactivity. Indeed, PV+ GPe neurons also exhibited abnormally elevated autonomous firing ex vivo. Optogenetic inhibition of PV+ GPe neurons in vivo partially and fully ameliorated the abnormal hypoactivity of postsynaptic subthalamic nucleus (STN) and putative PV– GPe neurons, respectively. In contrast to STN neurons whose autonomous firing is impaired in HD mice, putative PV– GPe neuron activity was unaffected ex vivo, implying that excessive inhibition was responsible for their hypoactivity in vivo. Together with previous studies, these data demonstrate that (1) indirect pathway nuclei are dysregulated in Q175 mice through changes in presynaptic activity and/or intrinsic cellular and synaptic properties; and (2) prototypic PV+ GPe neuron hyperactivity and excessive target inhibition are prominent features of early HD pathophysiology.
SIGNIFICANCE STATEMENT The early symptoms of Huntington's disease (HD) are linked to degenerative changes in the action-suppressing indirect pathway of the basal ganglia. Consistent with this linkage, the intrinsic properties of cells in this pathway exhibit complex alterations in HD and its models. However, the impact of these changes on activity is poorly understood. Using electrophysiological and optogenetic approaches, we demonstrate that the indirect pathway is highly dysregulated in early symptomatic HD mice through changes in upstream activity and/or intrinsic properties. Furthermore, we reveal that hyperactivity of external globus pallidus neurons and excessive inhibition of their targets are key features of early HD pathophysiology. Together, these findings could help to inform the development and targeting of viral-based, gene therapeutic approaches for HD.
Introduction
Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder caused by expansion of trinucleotide CAG repeats (>35) in exon 1 of the huntingtin gene (HTT) (Bates et al., 2015). The cardinal features of HD are progressive dysregulation and degeneration of the basal ganglia and cortex, and the emergence of associated motor, cognitive, and psychiatric symptoms (Bates et al., 2015). Although mutant huntingtin (mHTT) species impair multiple processes that are critical for neuronal function, the mechanisms that underlie the relative vulnerability of cortico-basal ganglia circuit function remain poorly understood (Reddy and Shirendeb, 2012; Seredenina and Luthi-Carter, 2012; Johri et al., 2013; Tong et al., 2014; Martin et al., 2015; Rosas et al., 2018). Under normal conditions, the basal ganglia promote contextually appropriate behavior, in part through the processing of functionally diverse cortical inputs (Mink and Thach, 1993; Nelson and Kreitzer, 2014; Klaus et al., 2019) by complex microcircuits, including the so-called direct and indirect pathways. The results of basal ganglia computation are then broadcast to the thalamus, midbrain, and brainstem, affecting behavior (Albin et al., 1989; Gerfen et al., 1990). D2 dopamine receptor expressing striatal projection neurons (D2-SPNs) that are the source of the so-called action-suppressing indirect pathway are susceptible early in HD (Albin et al., 1992; Richfield et al., 1995; Reiner and Deng, 2018). As a result, loss of indirect pathway function is posited to underlie chorea and impaired behavioral control, in the initial clinical phase of adult-onset HD (Bates et al., 2015; Reiner and Deng, 2018).
Examination of HD and its experimental models at early symptomatic ages has revealed that neurons in the basal ganglia and cortex exhibit complex changes in their intrinsic and synaptic properties (Cepeda et al., 2007; Raymond et al., 2011; Plotkin and Surmeier, 2015). For example, D2-SPNs exhibit loss of cortico-striatal long-term potentiation (LTP), and reductions in axospinous synapse density, miniature excitatory postsynaptic current (mEPSC) amplitude, and dendritic excitability (Plotkin and Surmeier, 2014; Plotkin et al., 2014; Sebastianutto et al., 2017; Carrillo-Reid et al., 2019). These observations suggest that D2-SPNs will be less effectively engaged by cortical excitation in HD mice. Consistent with this view, the activities of cortical neurons and SPNs in vivo are less well correlated in HD mice (Miller et al., 2008; Walker et al., 2008; Estrada-Sanchez et al., 2015). However, downstream external globus pallidus (GPe) neurons exhibit no sign of disinhibition in vivo, that is, they are not hyperactive (Beaumont et al., 2016), possibly because striatopallidal transmission has been upregulated (Perez-Rosello et al., 2019) and autonomous subthalamic nucleus (STN) activity has been downregulated (Atherton et al., 2016) in compensation. Thus, the effects of cortical drive on the indirect pathway cannot be easily inferred from ex vivo observations alone. Furthermore, most of the in vivo electrophysiological studies on HD mice are difficult to interpret because the recorded nuclei are comprised of multiple cell types that were not discriminated. To address these gaps, we used optogenetic- and reporter-guided electrophysiological interrogation of indirect pathway neurons in a well-characterized mouse model of HD at an early symptomatic stage. Specifically, we used Q175 model mice, in which exon 1 of human HTT (containing ∼190 trinucleotide repeats) is knocked into the endogenous Htt gene. This model exhibits many of the progressive molecular, neuropathological, and behavioral abnormalities seen in HD patients (Heikkinen et al., 2012; Menalled et al., 2012). In pilot experiments, we found that stereotyped patterns of cortical activity present under urethane anesthesia were similar in WT and Q175 mice, allowing us to compare firing in the indirect pathway during analogous cortical states. We also used patch-clamp recording in ex vivo brain slices to determine whether the changes in GPe neuron activity that we observed in Q175 mice in vivo were due in part to alterations in autonomous firing. Together with other studies, our data argue that mHTT profoundly dysregulates indirect pathway activity through complex changes in presynaptic activity and/or cellular and synaptic properties.
Materials and Methods
Animals
Procedures were performed in compliance with the policies of the National Institutes of Health and approved by the Institutional Animal Care and Use Committee of Northwestern University. Mice were maintained on a 14 h light/10 h dark cycle with food and water ad libitum, and monitored regularly by animal care technicians, veterinarians, and research staff.
Heterozygous Q175 mice (B6J.129S1-Htttm1.1Mfc/190ChdiJ; RRID:IMSR_JAX:029928; The Jackson Laboratory) were bred with homozygous PT-kj18-cre mice (Tg(Sim1-cre)KJ18Gsat; RRID:MMRRC_036958-UCD; MMRRC Repository), or A2A-cre mice (Tg(Adora2a-cre)KG139Gsat; RRID:MMRRC_036158-UCD; MMRRC Repository), or PV-cre mice (B6.Cg-Pvalbtm1.1(cre)Aibs/J; RRID:IMSR_JAX:012358; The Jackson Laboratory) to generate offspring that were either homozygous for Htt (WT), or heterozygous for Htt and mHtt (Q175), and heterozygous for Cre-recombinase. The following experimental mice were used (median and age range are reported): WT/PT-kj18-cre: age = 233, 225-239 d old, n = 5; Q175/PT-kj18-cre: age = 238, 228-244 d old, n = 4; WT/A2A-cre: age = 210, 210-233 d old, n = 3; Q175/A2A-cre: age = 211, 211-235 d old, n = 3; WT/PV-cre: age = 219.5, 204-232 d old, n = 4; Q175/PV-cre: age = 220, 205-234 d old, n = 5. Experimental mice were male, except for PV-cre mice, where both male (WT/PV-cre: n = 3; Q175/PV-cre: n = 3) and female (WT/PV-cre: n = 1; Q175/PV-cre: n = 2) mice were used. Data from male and female PV-cre mice were overlapping and therefore pooled.
Stereotaxic injection of viral vectors
Anesthesia was induced with vaporized 3%-4% isoflurane (Smiths Medical ASD) followed by an intraperitoneal injection of ketamine (100 mg/kg). After securing the mouse in a stereotaxic instrument (Neurostar), anesthesia was maintained with 1%-2% isoflurane. Adeno-associated viruses (AAVs) diluted in HEPES-buffered synthetic interstitial fluid (HBS SIF: 140 mm NaCl, 23 mm glucose, 15 mm HEPES, 3 mm KCl, 1.5 mm MgCl2, 1.6 mm CaCl2; pH 7.2 with NaOH; 300-310 mOsm/L) were then injected under stereotaxic guidance. AAV injection was conducted over 5-10 min at each site. An additional 5-10 min was then allowed for the injectate to diffuse before syringe retraction. ChR2(H134R)-eYFP was virally expressed in: (1) layer V pyramidal tract (PT) neurons in PT-kj18-cre mice through unilateral injection of AAV9.EF1a.DIO.hChR2(H134R)-eYFP.WPRE.hGH (RRID:Addgene_20298; 1 × 1013 genome copies/ml; AP: 0.6 mm, 1.2 mm, 1.8 mm; ML: 1.5 mm; DV: 1.0 mm; 0.5 μl per injection); and (2) D2-SPNs in A2A-cre mice through unilateral injections of AAV9.EF1a.DIO.hChR2(H134R)-eYFP.WPRE.hGH (RRID:Addgene_20298; 3 × 1011 genome copies/ml; AP: 0.4 mm, 0.9 mm; ML: 2.2 mm; DV: 3.7 mm, 2.7 mm; 0.3 μl per injection). Arch-GFP was virally expressed in parvalbumin-expressing (PV+) GPe neurons in PV-cre mice through unilateral injection of AAV9.CBA.Flex.Arch-GFP.WPRE.SV40 (RRID:Addgene_22222; 5 × 1011 genome copies/ml; AP: −0.27 mm; ML: 1.90 mm; DV: 3.95 mm, 3.45 mm; 0.25 μl per ventral injection and 0.20 μl per dorsal injection). eGFP was virally expressed in PV+ GPe neurons in PV-cre mice through unilateral injection of AAV9.Syn.DIO.eGFP.WPRE.hGH (RRID:Addgene_100043; 3 × 1011 genome copies/ml; AP: −0.27 mm; ML: 1.90 mm; DV: 3.95 mm, 3.45 mm; 0.25 μl per ventral injection and 0.20 μl per dorsal injection).
In vivo electrophysiological recording
Two to 4 weeks following surgery, anesthesia was induced with vaporized 3%-4% isoflurane followed by intraperitoneal injection of urethane (1.5 g/kg; Sigma-Aldrich). Mice were placed back into their home cage for ∼60 min until the toe-pinch withdrawal reflex was abolished. If necessary, additional urethane supplements (0.5 g/kg, i.p.) were administered every 30 min until the withdrawal reflex was eliminated. Mice were then placed into a stereotaxic instrument (David Kopf Instruments) for the duration of the recording session, with urethane supplements administered as required to maintain anesthesia. Craniotomies were drilled over the primary motor cortex (AP: 1.4 mm; ML: 1.5 mm), striatum (AP: 0.65 mm; ML: 1.95 mm), GPe (AP: −0.3 mm; ML: 2.0 mm), and/or STN (AP: −1.9 mm; ML: 1.4 mm) and irrigated with HBS. The intracranial electroencephalogram (EEG) was recorded from a peridural screw electrode (MS-51 960-1; McMaster-Carr) affixed over the ipsilateral primary motor cortex of A2A-cre and PV-cre mice. In PT-kj18-cre mice in which extracellular recordings of motor cortical neurons were made, the EEG screw was implanted over the contralateral primary motor cortex. Extracellular single-unit recordings were acquired using silicon tetrodes/optrodes (A1x4-tet-10 mm-100-121-A16 and A1x4-tet-10 mm-100-121-OA16, respectively; NeuroNexus Technologies) connected to a 64-channel Digital Lynx (Neuralynx) data acquisition system via a unity gain headstage, with a reference wire implanted adjacent to the ipsilateral temporal musculature. Signals were sampled at 40 kHz, with a gain of 14×. Online digital finite impulse response filters were applied. Single-unit activity was bandpass-filtered between 200 and 9000 Hz, and EEG and local field potential signals were bandpass-filtered between 0.1 and 400 Hz. Optogenetic stimulation was delivered using either a 473 nm diode laser (LuxX+ 473-100; Omicron-Laserage Laserprodukte) or a custom 577 nm laser system (Genesis MX STM 577-500 OPSL CW; Coherent). In order to histologically verify recording sites, silicon tetrodes/optrodes were dipped in a lipophilic florescent dye (DiI; 20 mg/ml in 50% acetone/methanol; D282; Thermo Fisher Scientific) before implantation. Sensory-evoked cortical activation (ACT) was generated by pinching the hindpaw for 5 s using a pair of fine forceps (Fine Science Tools). ChR2(H134R)-eYFP-expressing neurons were optogenetically stimulated with 473 nm light (<6 mW) for a duration of 5 ms. Stimulation was repeated at least 5 times with each trial of stimulation separated by a minimum of 2 min. Arch-GFP-expressing neurons were optogenetically inhibited through delivery of 577 nm light (<6 mW) for a duration of 5 s. Optogenetic inhibition was repeated at least 2 times with each trial being separated by a minimum of 2 min. Laser intensity was calibrated as power at the tip of the optrode before probe implantation and verified at the conclusion of each experiment.
Immunohistochemistry
After electrophysiological recording, mice were given a lethal dose of anesthetic and then perfused transcardially with ∼5-10 ml of 0.01 m phosphate buffered saline (PBS) (pH 7.4; P3813, Millipore Sigma) followed by 15-30 ml of 4% paraformaldehyde (PFA) in 0.1 m phosphate buffer (PB), pH 7.4. Each brain was then removed and postfixed overnight in 4% PFA (in 0.1 m PB, pH 7.4) before being washed in PBS, blocked, and sectioned in the coronal plane at 70 µm using a vibratome (VT1000S; Leica Microsystems). Sections were then processed for the immunohistochemical detection of NeuN, an antigen expressed by neurons that is commonly used to delineate brain structures. First, sections were washed in PBS and incubated for 48-72 h at 4°C in anti-NeuN (1:200; clone A60; Millipore Sigma; RRID:AB_2532109) in PBS with 0.5% Triton X-100 (Millipore Sigma) and 2% normal donkey serum (Jackson ImmunoResearch Laboratories). Then, sections were washed in PBS before being incubated for 90 min at room temperature in AlexaFluor-488- or 594-conjugated donkey anti-mouse IgG (1:250; Jackson ImmunoResearch Laboratories; RRID:AB_2313584; RRID:AB_2340621) in PBS with 0.5% Triton X-100 and 2% normal donkey serum. In a subset of PV-cre mice, in which PV+ GPe neurons expressed Arch-GFP or eGFP, adjacent sections of the GPe were processed for the immunohistochemical detection of PV (primary antibody: 1:1000 guinea pig anti-PV; Synaptic Systems; RRID:AB_2156476; secondary antibody: 1:250 AlexaFluor-594 donkey anti-guinea pig IgG; Jackson ImmunoResearch Laboratories; RRID:AB_2340474) or FoxP2 (primary antibody: 1:1000 rabbit anti-FoxP2; Millipore Sigma; RRID:AB_1078909; secondary antibody: 1:250 AlexaFluor-594 donkey anti-rabbit IgG; Jackson ImmunoResearch Laboratories; RRID:AB_2340621). Finally, sections were washed in PBS and mounted on glass slides with ProLong Diamond Antifade Reagent (P36965; Thermo Fisher Scientific). Mountant was allowed to cure for at least 24 h before storage at 4°C or imaging. DiI and immunofluorescent labeling were imaged using a Zeiss Axioskop 2 microscope (Carl Zeiss), equipped with an Axiocam CCD camera (426508-9901-000, Carl Zeiss), and Neurolucida software (MFB Bioscience). DiI and immunohistochemical labeling were used to map sites of recording (Fig. 1). Representative images were also acquired using confocal laser scanning microscopy (A1R; Nikon Instruments).
In vivo recording sites in the primary motor cortex, striatum, GPe, and STN. A–D, The locations of optogenetically responsive (O) and nonresponsive (X) neurons in WT (black) and Q175 (green) mice were mapped onto a standard atlas (Paxinos and Franklin, 2001) using the location of DiI-labeled recording tracks relative to cytoarchitectonic boundaries delineated by immunohistochemical labeling of NeuN. A, Recording sites in the primary motor cortex of PT-kj18-cre mice (M1, primary motor cortex; M2, secondary motor cortex). B, Recording sites in the striatum of A2A-cre mice (lv, lateral ventricle). C, Recording sites in the GPe of PV-cre mice (ic, internal capsule). D, Recording sites in the STN of PV-cre mice (zi, zona incerta).
In vivo electrophysiological analysis
Estimates of spectral power density were extracted using the Chronux data analysis toolbox (Bokil et al., 2010) for MATLAB (http://chronux.org/; The MathWorks). The EEG signal was downsampled to 1000 Hz, and spectral power was assessed at a resolution frequency of 0.061 Hz. Putative single-unit activity was discriminated with Plexon Offline Sorter software (version 3, Plexon; RRID:SCR_000012) using a combination of template matching, principal component analysis, and manual clustering. In addition, a threshold of <0.5% of interspike intervals under 2 ms was required for classification as a putative single unit (% interspike interval within 2 ms; cortical neurons: WT: 0.0, 0.0-0.0, n = 26; Q175: 0.0, 0.0-0.0, n = 43; striatal neurons: WT: 0.0, 0.0-0.0, n = 102; Q175: 0.0, 0.0-0.0, n = 65; GPe neurons: WT: 0.08, 0.02-0.29, n = 50; Q175: 0.05, 0.01-0.24, n = 55; STN neurons: WT: 0.23, 0.11-0.38, n = 30; Q175: 0.11, 0.00-0.19, n = 61; values represent median and interquartile range). Electrophysiological data were visually inspected in NeuroExplorer (Nex Technologies; RRID:SCR_001818) and then exported to MATLAB (The MathWorks; RRID:SCR_001622). Epochs with stable and robust slow-wave activity (SWA) or ACT were selected for analysis.
To determine whether neurons were responsive to optogenetic manipulation, peristimulus time histograms (PSTHs) were constructed from either 5 trials of hChR2(H134R)-eYFP stimulation or 2 or 3 trials of Arch-GFP stimulation. Cortical or striatal neurons were considered directly responsive if their activity exceeded the prestimulus mean by 2 SDs within 10 ms (bin size 2.5 ms) of hChR2(H134R)-eYFP stimulation. PV+ GPe neurons were considered directly responsive if their activity fell below 2 SDs of the prestimulus mean within 100 ms of Arch-GFP stimulation (bin size 100 ms). Putative PV– GPe neurons and STN neurons were considered to be disinhibited if their activity exceeded the prestimulus mean by 2 SDs within 100 ms of the optogenetic inhibition of PV+ GPe neurons (bin size 100 ms). To ensure that measurements were made from areas with opsin expression, analysis was restricted to recordings in which there was at least one responsive neuron on a given tetrode array. Direct or indirect responses to optogenetic stimulation were not observed when optrodes were activated in regions or mice without opsin expression (data not shown). Evoked firing was also easily discriminated from photo-electric artifacts, which were either minimal or absent.
Mean firing rates were calculated from the number of spikes divided by epoch length. The coefficient of variation (CV) of the interspike interval was used as a metric of regularity. To examine the relationship between cortical SWA and neuronal firing, phase histograms were generated in MATLAB. The EEG signal was first downsampled to 1000 Hz, and SWA was then extracted by applying a bandpass 0.1-1.5 Hz second-order Butterworth filter in the forward and reverse directions (to avoid phase shifts). The instantaneous phase of the EEG was calculated from the Hilbert transform (Le Van Quyen et al., 2001). In order to correct for the nonsinusoidal nature of slow cortical oscillations, the empirical cumulative distribution function (MATLAB) was applied (Siapas et al., 2005; Mallet et al., 2008b; Abdi et al., 2015; Kovaleski et al., 2020). Thus, each spike was assigned to a phase of the EEG from 0° to 360° (with 0°/360° and 180° corresponding to the peak active and inactive components of the EEG, respectively). Phase histograms were then constructed using 15° bins and plotted as median and interquartile range. In-phase activity was defined as activity occurring within a 180° window centered on 0°/360°. Anti-phase activity was defined as activity occurring within a 180° window centered on 180°.
Ex vivo electrophysiological recording
Mice were first lightly anesthetized with isoflurane, and then deeply anesthetized with ketamine/xylazine (87/13 mg/kg, i.p.), before being perfused transcardially with ∼10 ml of ice-cold sucrose-based ACSF (230 mm sucrose, 2.5 mm KCl, 1.25 mm NaH2PO4, 0.5 mm CaCl2, 10 mm MgSO4, 10 mm glucose, 26 mm NaHCO3, 1 mm sodium pyruvate, 5 μm L-glutathione; equilibrated with 95% O2 and 5% CO2). The brain was then removed, immersed in ice-cold 95% O2/5% CO2-equilibrated sucrose-based ACSF, and sectioned at 250 µm in the sagittal plane with a vibratome (VT1200S; Leica Microsystems). Slices were transferred to a holding chamber, immersed in ACSF (126 mm NaCl, 2.5 mm KCl, 1.25 mm NaH2PO4, 2 mm CaCl2, 2 mm MgSO4, 10 mm glucose, 26 mm NaHCO3, 1 mm sodium pyruvate, 5 μm L-glutathione; equilibrated with 95% O2 and 5% CO2), and maintained at 35°C for 30 min, and then room temperature. Next, individual slices were transferred to a recording chamber where they were perfused at 4-5 ml/min with SIF (126 mm NaCl, 3 mm KCl, 1.25 mm NaH2PO4, 1.6 mm CaCl2, 1.5 mm MgSO4, 10 mm glucose, 26 mm NaHCO3; equilibrated with 95% O2 and 5% CO2) at 35°C.
Patch-clamp recordings were made using 3-6 mΩ impedance, borosilicate glass electrodes filled with HBS SIF. Electrodes were positioned under visual guidance (Axioskop FS2, Carl Zeiss) using computer-controlled micromanipulators (Luigs and Neumann). Somatic recordings were made in the loose-seal, cell-attached configuration using an amplifier (MultiClamp 700B; Molecular Devices), and an associated digitizer (Digidata 1440A; Molecular Devices) controlled by pCLAMP 10.3 (Molecular Devices). Electrode capacitance was compensated, and signals were low-pass filtered online at 10 kHz and sampled at 25 kHz. Recordings of autonomous action potential generation were made in the presence of 20 μm DNQX, 50 μm D-AP5, 10 μm Gabazine (SR-95531), and 2 μm CGP 55845 to block synaptic transmission at AMPA, NMDA, GABAA, and GABAB receptors, respectively. All drugs used were purchased from Hello Bio and bath-applied.
PV+ GPe neurons were identified through visualization of eGFP under 460 nm light-emitting diode epifluorescent illumination (OptoLED; Cairn Research). To increase the probability that GPe neurons that did not express eGFP were PV– GPe neurons, non-eGFP-expressing GPe neurons were only recorded if eGFP-expressing PV+ GPe neurons were present in the same field of view. The frequency and regularity of eGFP-expressing PV+ GPe neuron activity and non–eGFP-expressing GPe neuron activity were calculated from 30 s recording epochs.
Experimental design and statistical analyses
Data are reported as median and interquartile range. Data are represented graphically as violin (kernel density) plots and overlaid box plots, with the median (central line), interquartile range (box), and 10%-90% range (whiskers) denoted. The number and nature of observations for each parameter are specified throughout. To ensure that the proposed research was adequately powered, sample sizes were estimated using the formulae described by Noether (1987) assuming 80% power (i.e., a 20% probability of a Type 2 error) and a two tailed α level of 0.05. For unpaired data (groups X and Y), and probabilities of X > Y (or X < Y) being 0.7, 0.8, and 0.9, the estimated sample sizes for each group are 33, 15, and 9, respectively. For paired data (where Xi and Xj are independent samples from X, reflecting effect size and sign) and the probabilities of Xi + Xj > 0 being 0.7, 0.8, and 0.9, the estimated sample sizes are 66, 30, and 17, respectively. Probabilities between 0.7 and 0.9 are representative of our historical and pilot data. To minimize assumptions concerning the distribution of data, nonparametric, two-tailed statistical comparisons were made using the Mann–Whitney U (MWU) and Wilcoxon signed-rank (WSR) tests for unpaired and paired comparisons, respectively. In addition, Fisher's exact test was used for contingency analyses. p < 0.05 was considered significant. Where appropriate, p values were adjusted for multiple comparisons using the Holm–Bonferroni method. Plots and statistical comparisons were generated in Prism (GraphPad Software; RRID:SCR_002798) and R (https://www.r-project.org/; RRID:SCR_001905).
Results
Cortical activity is similar in Q175 and WT mice during both SWA and ACT
Under anesthesia, the cortex can exhibit the following: (1) robust ∼1 Hz cortical SWA during which cortical projection neurons, including those projecting to the basal ganglia, exhibit synchronous transitions between hyperpolarized quiescent and depolarized active states, analogous to those during deep sleep (Stern et al., 1997; Steriade, 2000); and (2) spontaneous or somatosensory stimulation-triggered cortical activation ACT, during which cortical projection neurons exhibit desynchronized activity, analogous to that during arousal (Stern et al., 1997; Steriade, 2000). As a result, anesthesia has been routinely used to generate stereotyped cortical activity and probe its impact on the downstream basal ganglia in normal rodents and psychomotor disease models (Magill et al., 2001; Walters et al., 2007; Mallet et al., 2008a,b; Zold et al., 2012; Sharott et al., 2017; Aristieta et al., 2020; Kovaleski et al., 2020; Ketzef and Silberberg, 2021). Therefore, to determine whether indirect pathway activity is dysregulated in 6-month-old Q175 mice relative to WT age-matched controls, we compared neuronal activity during both cortical SWA and sensory-evoked cortical ACT under urethane anesthesia.
Cortical network activity was assessed from the intracranial EEG, which was obtained from a peridural screw “electrode” affixed over primary motor cortex. Cortical SWA was manifest in the EEG as a high-amplitude, low-frequency (∼1 Hz) oscillation on which phase-locked, low-amplitude, high-frequency oscillations were superimposed (Fig. 2A). Cortical ACT occurred spontaneously or could be triggered by somatosensory stimulation, which was signified in the EEG by diminution of the 1 Hz oscillation and persistence of higher-frequency oscillations (Fig. 2A). During cortical SWA, the powers of motor cortical oscillations in frequency bands ranging from 0 to 100 Hz were similar in Q175 and WT mice (Fig. 2B–E; Table 1). In addition, low-frequency and high-frequency cortical oscillations were attenuated and elevated, respectively, to a similar degree in Q175 and WT mice during hind paw pinch-evoked cortical ACT (Fig. 2B–E; Table 1). Together, these data suggest that motor cortical network activity is similar in 6-month-old Q175 and WT mice.
Motor cortical EEG measurements during cortical SWA and ACT in Q175 and WT micea
Motor cortical network activity is similar in Q175 and WT mice during both SWA and ACT. A–E, The motor cortical EEGs of urethane-anesthetized Q175 and WT mice were similar during SWA and ACT (A, examples; B–E, population data). Spectral power was similar for both low- (B, C) and high-frequency (D, E) EEG components. *p < 0.05.
Layer V cortical pyramidal neurons, some of which innervate the basal ganglia, are comprised of two major cell classes, PT and intratelencephalic type (IT) neurons (Harris and Shepherd, 2015). To compare the activity of PT neurons in Q175 and WT mice, an “opto-tagging” approach was used to identify the firing of PT neurons in vivo. First, ChR2(H134R)-eYFP was virally expressed in PT cortical neurons through injection of an AAV carrying a cre-dependent expression construct into the primary motor cortex of Q175 and WT mice that had been crossed with a PT neuron selective cre-driver line (PT-kj18-cre) (Fig. 3A–D). Two to 4 weeks later, the activity of cortical neurons was compared using an array of tetrodes fiber-coupled to a laser. Construction of peri-optogenetic stimulus time histograms revealed that the latency of firing evoked by optogenetic stimulation was unimodal in distribution (Fig. 3E,F). The majority of cells fired within 10 ms of the start of the 5 ms optogenetic stimulation pulse (latency; WT: 2.5, 2.5-5.0 ms, n = 15; Q175: 2.5, 2.5-3.7 ms, n = 21; values represent median and interquartile range) (Fig. 3E,F). These response latencies are consistent with the direct optogenetic stimulation of ChR2-expressing neurons (Hooks, 2018; Li et al., 2018). Therefore, neurons that fired within 10 ms of optogenetic stimulation were positively identified as PT neurons (Fig. 3D–F). Consistent with this classification, all identified PT neurons exhibited action potential properties that are typical of cortical pyramidal neurons (Mitchell et al., 2007; Kaufman et al., 2010; Takahashi et al., 2015; Lohani et al., 2019) (Fig. 3G,H). The proportion of PT neurons that exhibited short-latency responses to optogenetic activation and the latency of those evoked responses (Table 2) were similar in Q175 and WT mice. Neurons that were recorded on the same tetrode as identified PT neurons but were unresponsive to optogenetic stimulation or exhibited longer-latency presumably indirect, monosynaptic responses (WT: n = 11; Q175: n = 22) were also recorded. The majority of neurons (WT: n = 7; Q175: n = 15) that were unresponsive to optogenetic stimulation or exhibited latencies of discharge > 10 ms also exhibited action potential properties that are typical of cortical pyramidal neurons, consistent with the relative abundance of pyramidal versus interneurons in this cortical layer (Fig. 3G,H). The neurons with pyramidal-like action potential properties were termed unidentified layer V pyramidal neurons and likely comprise IT neurons and PT neurons that failed to express ChR2(H134R)-eYFP. The unidentified layer V pyramidal neurons that exhibited long-latency discharge in response to optogenetic stimulation are likely to be PT neurons because primary motor cortical PT neurons preferentially innervate each other rather than IT neurons (Kiritani et al., 2012). A minority of unidentified layer V neurons that exhibited interneuron action potential properties (Fig. 3G,H) were excluded from further description or analysis here because of their rarity.
Layer V motor cortical projection neuron activity during cortical SWA and ACT in Q175 and WT micea
Optogenetic identification of layer V motor cortical projection neurons. A, B, Cre-dependent viral expression of ChR2(H134R)-eYFP (green) in layer V PT neurons (M1, primary motor cortex; DV, dorsoventral axis; ML, mediolateral axis). B, ChR2(H134R)-eYFP expression was present in a subset of layer V, presumably PT neurons (white arrows), and absent in other layer V neurons (magenta arrows), which presumably comprise a mixture of nonexpressing PT, IT, and interneurons. Cortical neurons were also immunohistochemically labeled for the pan neuronal marker NeuN (red). C, Schematic of experimental setup illustrating placement of the EEG screw electrode (gray) and optrode (blue). D, An optogenetically stimulated layer V PT neuron and an adjacent unresponsive, unidentified (UI) neuron (blue, optogenetic stimulation; inset, spontaneous [black] and optogenetically evoked [blue] action potentials). E, F, Cortical neurons were considered to be directly responsive if their activity exceeded the prestimulus mean (solid red line) by 2 SDs (dotted red line) within 10 ms (bin size, 2.5 ms) of the start of the optogenetic stimulation pulse (E, PSTHs from examples in D; F, population latency to first spike histograms after the onset of stimulation). G, H, Action potential properties in responsive, opto-tagged PT (O) and non-responsive, UI (X) layer V cortical neurons in WT (black) and Q175 (green) mice. Spike waveforms were measured from the average extracellular waveform of each unit and classified by their peak-to-valley width [1] and valley width at its half-maximum [2]. Broad and narrow waveforms were used to discriminate pyramidal neurons from interneurons, respectively.
Consistent with previous studies, PT and unidentified layer V pyramidal neurons fired preferentially during the active component of cortical SWA (Steriade et al., 1993; Amzica and Steriade, 1998; Beltramo et al., 2013) (Fig. 4A). During cortical SWA, the frequency and regularity of PT and unidentified layer V pyramidal neuron activity were similar in Q175 and WT mice (Fig. 4A,C–E; Table 2). In addition, pinch-evoked cortical desynchronization led to a similar reduction in PT and unidentified layer V pyramidal neuron activities in both Q175 and WT mice (Fig. 4B,D,E; Table 2). Together, the EEG and single-unit data suggest that 6-month-old Q175 and WT mice exhibit similar patterns and levels of motor cortical activity under urethane anesthesia.
Layer V motor cortical projection neuron activity in Q175 and WT mice. A–E, The activities of layer V PT and unidentified pyramidal neurons in Q175 and age-matched WT mice were similar during cortical SWA and ACT. A, B, Representative examples of concurrent EEG (bandpass-filtered at 0.1-1.5 Hz) and cortical neuron activity during SWA (A) and ACT (B). The frequencies of firing of PT and unidentified pyramidal neurons were similar in both genotypes during cortical SWA (A, examples; C–E, population data). B, D, E, Hind paw pinch-evoked cortical ACT (dotted line) similarly affected PT and unidentified layer V pyramidal neuron activity in Q175 and WT mice (B, examples; D, E, population data). *p < 0.05.
D2-SPNs are hypoactive in Q175 mice during cortical SWA
The striatum is largely composed of similar numbers of direct and indirect pathway SPNs that express D1 or D2 dopamine receptors, respectively (Albin et al., 1989; Gerfen et al., 1990). As their names suggest, direct pathway D1-SPNs directly innervate basal ganglia output neurons, whereas indirect pathway D2-SPNs regulate basal ganglia output indirectly via the GPe and STN (Mink and Thach, 1993; Maurice et al., 1999; Tachibana et al., 2008). Appropriate, cortical patterning of D1- and D2-SPN activity is critical for the regulation of psychomotor function by the basal ganglia (Cui et al., 2013; Keeler et al., 2014; Nelson and Kreitzer, 2014; Tecuapetla et al., 2014; Sippy et al., 2015; Barbera et al., 2016; Lambot et al., 2016; Lemos et al., 2016; Tecuapetla et al., 2016; Klaus et al., 2019; LeBlanc et al., 2020). Although D2-SPNs are a key locus of dysregulation and degeneration in HD and its models, precisely how their activity is perturbed in vivo is unknown. To address this question, ChR2(H134R)-eYFP was virally expressed in D2-SPNs through injection of an AAV carrying a cre-dependent expression construct into the striatum of Q175 and WT mice that had been crossed with a D2-SPN selective cre-driver line (A2A-cre). Two to 3 weeks later, the activities of striatal neurons in Q175 and WT mice were compared using silicon optrodes. Consistent with the successful targeting of D2-SPNs in both Q175 and WT mice, ChR2(H134R)-eYFP was expressed in the subset of striatal neurons projecting to the GPe but not to the substantia nigra pars reticulata (Fig. 5A–D). Therefore, optogenetic activation was used to positively identify D2-SPNs (Fig. 5E–G). Construction of peri-optogenetic stimulus time histograms revealed that the latency of firing evoked by optogenetic stimulation was unimodal in distribution (Fig. 5F,G). Consistent with the direct optogenetic stimulation of ChR2 in D2-SPNs, the majority of neurons discharged within 10 ms of the start of the 5 ms optogenetic stimulation pulse (latency; WT: 5.0, 2.5-5.0 ms, n = 54; Q175: 5.0, 2.5-5.0 ms, n = 42) (Fig. 5F,G) and exhibited SPN-like rather than striatal interneuron-like action potential properties (Berke et al., 2004; Gage et al., 2010; Cayzac et al., 2011; Kim et al., 2014; Shin et al., 2018) (Fig. 5H).
Optogenetic identification of D2-SPNs. A–C, Viral-mediated, cre-dependent expression of hChR2(H134R)-eYFP (green) in D2-receptor expressing striatopallidal neurons in WT/Q175 X A2A-cre mice in the striatum (A, B, lv, lateral ventricle; DV, dorsoventral axis; ML, mediolateral axis) and GPe (C, ic, internal capsule; same orientation as in A). B, hChR2(H134R)-eYFP expression was present (white arrows) or absent (magenta arrows) in striatal neurons that were coimmunoreactive for the pan neuronal marker NeuN (red). C, Axon-terminal expression of hChR2(H134R)-eYFP in the GPe but not the substantia nigra pars reticulata (SNr), consistent with selective expression in D2- but not D1-SPNs (ic, internal capsule). D, Schematic of the experimental setup, illustrating striatal placement of the optrode. E, Example of optogenetic activation of a D2-SPN but not a neighboring unidentified (UI) striatal neuron (blue, optogenetic stimulation; inset, spontaneous [black] and optogenetically evoked [blue] action potentials). F, G, Striatal neurons were considered to be directly responsive if their activity exceeded the prestimulus mean (solid red line) by 2 SDs (dotted red line) within 10 ms (bin size, 2.5 ms) of the start of the optogenetic stimulation pulse (F, PSTHs from examples in E; G, population latency to first spike histograms after the onset of stimulation). H, Action potential properties of responsive, opto-tagged D2-SPNs and non-responsive, UI striatal neurons in WT (black) and Q175 (green) mice. Broad and narrow waveforms were used to discriminate SPNs from interneurons, respectively.
The proportion of striatal neurons that exhibited short-latency excitatory responses to optogenetic stimulation was similar in Q175 and WT mice (Table 3). Neurons that were recorded on the same tetrode as identified D2-SPNs but were unresponsive to optogenetic stimulation (WT: n = 48; Q175: n = 23) were also recorded. The majority of neurons that were unresponsive to optogenetic stimulation exhibited action potential properties that are typical of SPNs (WT: n = 45; Q175: n = 19) and likely comprise D1-SPNs and a small fraction of D2-SPNs that failed to express ChR2(H134R)-eYFP (Fig. 5H). Thus, these neurons were termed unidentified, putative D1-SPNs. A minority of unresponsive neurons exhibited interneuron-like action potential properties (Fig. 5H) and were excluded from further description or analysis here because of their rarity. Optogenetic responses were not observed when optrodes were activated in regions or mice where ChR2(H134R)-eYFP expression was absent (data not shown).
D2-SPN and unidentified, putative D1-SPN activity during cortical SWA and ACT in Q175 and WT micea
Consistent with direct cortical driving, D2-SPN and unidentified putative D1-SPN activity were phase-locked to the active component of cortical SWA in both Q175 and WT mice (Fig. 6A). The frequency of D2-SPN firing in Q175 mice was approximately one-half to two-thirds lower than that in WT control mice during cortical SWA, consistent with reductions in axospinous synapse density, mEPSC amplitude, and dendritic excitability in D2-SPNs (Fig. 6B,C; Table 3). In contrast, the frequency of unidentified, putative D1-SPN activity was similar in Q175 and WT mice during cortical SWA (Fig. 6B,C; Table 3). The phase-locking of D2-SPN and unidentified striatal neuron activity to cortical SWA was similar in WT and Q175 mice (Fig. 6D–F; Table 3).
D2-SPNs are hypoactive in Q175 mice during cortical SWA. A–C, The frequency of D2-SPN firing was lower in Q175 mice relative to age-matched WT mice during cortical SWA (A, examples; B, C, population data). The frequency of unidentified (UI), putative D1-SPN activity during SWA was similar in Q175 and age-matched WT mice (A, examples; B, C, population data). D–F, D2-SPN activity and UI, putative D1-SPN activity were similarly phase-locked to the active component of cortical SWA in Q175 and WT mice (D, population data; E, F, population spike phase histograms). G–I, Effect of pinch-evoked cortical ACT (dotted line) on both D2-SPNs and UI, putative D1-SPNs in Q175 and WT mice (G, examples; H, I, population data; an outlier data point from the WT D2-SPN group was not plotted for legibility). During cortical ACT, the frequencies of D2-SPN and UI, putative D1-SPN activity were higher in Q175 mice (H, I). *p < 0.05.
Pinch-evoked cortical ACT decreased D2-SPN firing in both genotypes (Fig. 6G,H; Table 3). However, the frequency of D2-SPN discharge was modestly but significantly elevated in Q175 relative to WT mice during cortical ACT (Fig. 6G,H; Table 3). During pinch-evoked cortical ACT, the frequency of unidentified, putative D1-SPN activity declined in WT but not in Q175 mice. As a result, the activity of unidentified, putative D1-SPN neurons was elevated in Q175 mice (Fig. 6G,I; Table 3). Together, these data demonstrate that (1) during cortical SWA, D2-SPNs are hypoactive in Q175 mice; and (2) during cortical ACT, both D2-SPNs and unidentified putative D1-SPNs are more active in Q175 mice.
Prototypic PV+ GPe neurons are hyperactive in Q175 mice
To determine whether there are changes in the activity of GPe neurons downstream of D2-SPNs, we compared the in vivo firing of GPe neurons in Q175 and age-matched WT control mice. The GABAergic projection neurons of the GPe have been divided into two major types, so-called prototypic and arkypallidal neurons (Mallet et al., 2012; Mastro et al., 2014; Abdi et al., 2015; Dodson et al., 2015; Hernandez et al., 2015). Prototypic neurons comprise approximately three-fourths of all GPe neurons, are preferentially innervated by D2-SPNs, and innervate arkypallidal and downstream basal ganglia neurons (Mastro et al., 2014; Abdi et al., 2015; Dodson et al., 2015; Hernandez et al., 2015; Aristieta et al., 2020; Ketzef and Silberberg, 2021). A subset of prototypic neurons also innervates the striatum and/or cortex (Bevan et al., 1998; Mastro et al., 2014; Saunders et al., 2016; Abecassis et al., 2020). The majority of prototypic GPe neurons express the calcium-binding protein PV (Mastro et al., 2014; Abdi et al., 2015; Dodson et al., 2015; Hernandez et al., 2015). In contrast, arkypallidal neurons comprise approximately one-fourth of GPe neurons, are preferentially innervated by D1-SPNs rather than D2-SPNs, innervate the striatum rather than prototypic GPe neurons or the downstream basal ganglia, and express the transcription factor forkhead box protein 2 (FoxP2) but not PV (Mallet et al., 2012; Mastro et al., 2014; Abdi et al., 2015; Dodson et al., 2015; Hernandez et al., 2015; Aristieta et al., 2020; Ketzef and Silberberg, 2021). Although arkypallidal neurons are more weakly innervated by D2-SPNs compared to prototypic neurons, D2-SPNs can powerfully regulate arkypallidal activity indirectly via their effects on prototypic GPe neurons (Aristieta et al., 2020; Ketzef and Silberberg, 2021).
Because the majority of GPe neurons fire tonically, we used optogenetic silencing to identify and manipulate GPe neurons in vivo. The inhibitory opsin Arch-GFP was virally expressed in prototypic PV+ GPe neurons through injection of an AAV vector carrying a cre-dependent expression construct into the GPe of Q175 and WT mice that had been crossed with a PV-cre driver line (Fig. 7A–D). Consistent with the selective expression of Arch-GFP in prototypic GPe neurons, the majority of Arch-GFP-expressing GPe neurons were immunoreactive for PV, and all were immunonegative for FoxP2 (Fig. 7B). In addition, Arch-GFP was expressed in the axon terminal fields of PV+ GPe neurons in the downstream basal ganglia, including the STN, consistent with their prototypic identity (Fig. 7C). Optogenetic inhibition was therefore used to positively identify PV+ GPe neurons. Optogenetic activation of Arch-GFP rapidly and persistently inhibited a similar proportion of neurons in the GPe of Q175 and WT mice (Fig. 7E–H; Table 4). Optogenetic responses were not observed when optrodes were activated in regions or mice where Arch-GFP expression was absent (data not shown).
PV+ GPe neuron activity during cortical SWA and ACT in Q175 and WT micea
Optogenetic identification of prototypic PV+ GPe neurons. A–C, Viral cre-dependent expression of Arch-GFP (green) in prototypic PV+ GPe neurons in Q175/WT X PV-cre mice. A, Representative confocal micrograph of the GPe (ic, internal capsule; DV, dorsoventral axis; ML, mediolateral axis). B, Arch-GFP expression was present in PV-immunoreactive GPe neurons (white arrows; top) and FoxP2-immunonegative neurons (white arrows; bottom) and absent in FoxP2-immunoreactive arkypallidal GPe neurons (magenta arrows; bottom). C, Arch-GFP-labeled axon terminal fields in the STN (top; zi, zona incerta; same orientation as in A). Arch-GFP-labeled axon terminal fields in the vicinity of NeuN-immunoreactive STN neurons (red; bottom). D, Schematic of the experimental setup illustrating optrode placement in the GPe. E–H, Optogenetic activation of Arch-GFP rapidly and persistently inhibited PV+ GPe neuron activity. E, Representative example of concurrent EEG (bandpass-filtered at 0.1-1.5 Hz) and PV+ GPe neuron activity before, during, and after activation of Arch-GFP (yellow bar). F, G, PSTHs of PV+ GPe neuron activity in the absence and presence (yellow bar) of optogenetic inhibition (bin size, 100 ms; prestimulus mean, solid red line; ±2 SDs of prestimulus mean, dotted red line; F, PSTH from neuron in E; G, population PSTH). H, Direct responses of all PV+ GPe neurons to activation of Arch-GFP.
As previously reported (Abdi et al., 2015; Kovaleski et al., 2020), PV+ GPe neurons discharged tonically in WT mice during cortical SWA (Fig. 8A). Although PV+ GPe neurons also discharged tonically in Q175 mice during cortical SWA, their median activity was more than double that in WT mice (Fig. 8B; Table 4). The regularity of PV+ GPe neuron firing, as assessed from the CV of the interspike interval, was similar in Q175 and WT mice (Fig. 8C; Table 4). PV+ GPe neuron firing was elevated during the active component of cortical SWA in Q175 mice (Fig. 8D; Table 4), consistent with the hypoactivity of upstream D2-SPNs, as described above. However, PV+ GPe neurons were also hyperactive in Q175 mice during the inactive component of cortical SWA (Fig. 8D; Table 4). Indeed, the ratio of PV+ GPe neuron activity during the inactive versus active component of cortical SWA was significantly higher in Q175 mice (Fig. 8D,E; Table 4). Given that D2-SPNs discharge only during the active phase of cortical SWA in WT and Q175 mice (Fig. 6D), the relative hyperactivity of PV+ GPe neurons in Q175 mice during the inactive phase of cortical SWA cannot be due to D2-SPN hypoactivity. During pinch-evoked cortical ACT, the firing rate of PV+ GPe neurons increased in both genotypes (Fig. 8F,G; Table 4). However, PV+ GPe neurons remained strongly hyperactive in Q175 mice relative to WT controls (Fig. 8F,G; Table 4), despite the fact that D2-SPNs are hyperactive in Q175 mice during cortical ACT (Fig. 6G,H). Together, these data demonstrate that in Q175 mice prototypic PV+ GPe neuron activity is highly elevated during both cortical SWA and ACT.
In Q175 mice, PV+ GPe neurons are hyperactive during cortical SWA and ACT. A, Representative examples of concurrent EEG (bandpass-filtered at 0.1-1.5 Hz) and PV+ GPe neuron activity in the absence and presence of optogenetic inhibition (yellow bar). A–C, During cortical SWA, the median frequency of PV+ GPe neuron firing in Q175 mice was more than twice that in WT mice, but the variability of firing was unaltered (A, representative examples; B,C, population data). D, E, In Q175 mice, the firing of PV+ GPe neurons was relatively antiphasic to cortical SWA (D, population spike phase histograms; E, population box and violin plots). F, G, Pinch-evoked cortical ACT (hind paw pinch, dotted line) increased PV+ GPe neuron firing in both genotypes. However, the frequency of PV+ GPe neuron discharge remained elevated in Q175 mice relative to WT mice (F, representative examples; G, population data). *p < 0.05.
The abnormal hypoactivity of PV– GPe neurons in Q175 mice is reversed by optogenetic inhibition of PV+ GPe neurons
We next analyzed the activity of GPe neurons that were not directly inhibited during optogenetic activation of Arch-GFP in PV-cre mice (Fig. 9; Table 5). These most likely comprise PV– prototypic GPe and arkypallidal neurons but could also include a small fraction of PV+ GPe neurons that failed to express Arch-GFP. During cortical SWA, unidentified, putative PV– GPe neurons were less active than PV+ GPe neurons in both WT and Q175 mice, consistent with previous reports in WT rodents (Fig. 9A,D; Table 5) (Abdi et al., 2015; Dodson et al., 2015; Mallet et al., 2016). Putative PV– GPe neurons also discharged more slowly in Q175 than WT mice (Fig. 9A,D,F; Table 5). In contrast, the precision of putative PV– GPe neuron activity was similar in WT and Q175 mice (Fig. 9E; Table 5). Recent studies suggest that prototypic PV+ and PV– GPe neurons inhibit each other and arkypallidal neurons through powerful local connections (Aristieta et al., 2020; Ketzef and Silberberg, 2021). In contrast, arkypallidal neuron innervation of prototypic GPe neurons is minimal (Aristieta et al., 2020; Ketzef and Silberberg, 2021). If elevated inhibition emanating from hyperactive prototypic PV+ GPe neurons is responsible for the hypoactivity of putative PV– GPe neurons in Q175 mice, then their activity should be rescued by optogenetic inhibition of PV+ GPe neurons. Consistent with this hypothesis, optogenetic silencing of PV+ GPe neurons during cortical SWA rapidly and persistently disinhibited putative PV– GPe neuron activity in both WT and Q175 mice, confirming that PV+ GPe neurons actively inhibit PV– GPe neuron activity in vivo (Fig. 9A–C,F; Table 5). Furthermore, during optogenetic inhibition of PV+ GPe neurons, the activities of putative PV– GPe neurons in WT and Q175 mice were no longer significantly different (Fig. 9A–C,F; Table 5), arguing that increased inhibition from hyperactive PV+ GPe neurons was indeed responsible for their relative hypoactivity in Q175 mice. There were no consistent differences in the phase of putative PV– GPe neuron activity (relative to cortical SWA) with or without optogenetic inhibition of PV+ GPe neurons in WT and Q175 mice (Fig. 9G,H; Table 5).
Impact of optogenetic inhibition of PV+ GPe neurons on putative PV– GPe neuron activity in Q175 and WT micea
The hypoactivity of putative PV– GPe neurons in Q175 mice is alleviated by optogenetic inhibition of hyperactive PV+ GPe neurons. A, Representative examples of concurrent EEG (bandpass-filtered at 0.1-1.5 Hz) and GPe neuronal activity in Q175/WT X PV-cre mice during optogenetic inhibition (yellow bar) of PV+ GPe neurons. A–C, Optogenetic inhibition of PV+ GPe neuron activity rapidly and persistently disinhibited putative PV– GPe neurons in Q175 and WT mice. B, C, PSTHs of putative PV– GPe neuron activity in the absence and presence of PV+ GPe neuron optogenetic inhibition (yellow bar) (bin size, 100 ms; prestimulus mean, solid red line; ±2 SDs of prestimulus mean, dotted red line; B, PSTHs from representative neurons in A; C, population PSTHs). A–F, During cortical SWA, putative PV– GPe neurons were less active than PV+ GPe neurons in both WT and Q175 mice. Furthermore, the frequency of putative PV– GPe neuron activity in Q175 mice was lower than in WT mice (A, representative examples; D, population data). E, The precision of PV– GPe neuron activity was similar in WT and Q175 mice. Optogenetic inhibition of PV+ GPe neurons disinhibited putative PV– GPe neurons in WT and Q175 mice and eliminated the difference in firing frequencies between the two genotypes, arguing that GABAergic inhibition emanating from abnormally hyperactive PV+ GPe neurons is responsible for the relative hypoactivity of PV– GPe neurons in Q175 mice (A, examples; F, population data). G, H, Relative to cortical SWA, there were no differences in the phase of putative PV– GPe neuron activity in WT and Q175 mice, with or without optogenetic inhibition of PV+ GPe neurons. *p < 0.05.
The autonomous activity of PV+ GPe neurons but not putative PV– GPe neurons is increased in brain slices from Q175 mice
The high rates of discharge of extrastriatal basal ganglia neurons in vivo are generated in part by their intrinsic autonomous activity (Wilson, 2013). Alterations in the autonomous firing of PV+ and PV– GPe neurons could therefore contribute to their abnormal in vivo activity in Q175 mice. To determine whether the autonomous firing of GPe neurons is altered in Q175 mice, an AAV vector carrying a cre-dependent eGFP-expression construct was injected into the GPe of WT or Q175 PV-cre mice. Two to 3 weeks later, brain slices were prepared and visually guided somatic patch-clamp recordings of GPe neurons were conducted. Expression of eGFP was used to identify prototypic PV+ GPe neurons. PV– GPe neurons in the vicinity of eGFP-expressing PV+ GPe neurons were identified putatively by their absence of eGFP expression. GPe neurons were recorded in the loose-seal, cell-attached, current-clamp configuration and recordings were made in the presence of AMPA, NMDA, GABAA, and GABAB receptor antagonists to minimize the impact of synaptic inputs on autonomous firing. As for previous studies (Abdi et al., 2015; Hernandez et al., 2015), we found that the frequency and precision of autonomous firing were greater and autonomously generated action potentials were briefer in prototypic PV+ GPe neurons than putative PV– GPe neurons in WT mice (Fig. 10A–F; Table 6). Thus, these data are consistent with the appropriate identification of PV+ and PV– GPe neurons rather than the failure of reporter expression in PV+ GPe neurons. In Q175 mice, the frequency and precision of autonomous PV+ GPe neuron activity were greater than in WT mice (Fig. 10A,E,F; Table 6). In contrast, there were no differences in the frequency or precision of autonomous putative PV– GPe neuron firing in Q175 and WT mice (Fig. 10B,E,F; Table 6). Together, these data argue that the autonomous firing of prototypic PV+ GPe neurons is upregulated in Q175 mice and that this alteration may contribute to the hyperactivity of these cells in vivo. In contrast, the similarity of autonomous putative PV– GPe neuron firing in Q175 and WT mice further supports the conclusion that excessive inhibition emanating from hyperactive PV+ GPe neurons is responsible for their hypoactivity in Q175 mice in vivo.
The autonomous activity of GPe neurons in ex vivo brain slices from Q175 and WT micea
The autonomous activity of PV+ GPe neurons is elevated in Q175 mice. A, B, Representative examples of autonomous GPe neuron activity recorded in the loose-seal, cell-attached, current-clamp configuration in brain slices from Q175/WT X PV-cre mice. C, D, Average spike waveforms were significantly shorter in duration in PV+ GPe neurons than putative PV– GPe neurons (C, representative examples; D, population data). A, E, F, The frequency and precision of autonomous PV+ GPe neuron activity were greater in slices from Q175 mice (A, representative examples; E, F, population data). B, E, F, In contrast, the frequency and precision of autonomous PV– GPe neuron activity were similar in Q175 and WT mice (B, representative examples; E, F, population data). *p < 0.05.
The hypoactivity of the STN in Q175 mice is partially alleviated by optogenetic inhibition of prototypic PV+ GPe neurons
The glutamatergic STN is a key component of the indirect pathway and forms a reciprocally connected network with the GPe (Mink and Thach, 1993; Maurice et al., 1999; Plenz and Kital, 1999; Nambu et al., 2002; Tachibana et al., 2008). Previous research has demonstrated that autonomous STN activity is downregulated in HD mice (Atherton et al., 2016). Furthermore, we demonstrate here that GABAergic prototypic PV+ GPe neurons, which potently inhibit STN activity through their activation of postsynaptic GABAA and GABAB receptors (Bevan et al., 2002; Hallworth and Bevan, 2005; Baufreton et al., 2009; Atherton et al., 2013; Kovaleski et al., 2020), are hyperactive in Q175 mice in vivo. Given the loss of intrinsic STN activity and increased frequency of GABAergic GPe-STN transmission, we predicted that STN neurons will be hypoactive in Q175 mice in vivo compared with WT mice. To test this prediction, we recorded the activity of STN neurons with silicon tetrodes during cortical SWA and ACT in Q175 and WT mice crossed with PV-cre mice. To isolate the effect of PV+ GPe neuron inhibition on the STN, we virally expressed Arch-GFP in PV+ GPe neurons and inhibited their activity optogenetically, as described above (Fig. 11A). As for previous studies in WT rodents, STN neurons fired in phase with cortical SWA in WT and Q175 mice (Magill et al., 2000, 2001; Mallet et al., 2008a; Callahan and Abercrombie, 2015a,b; Kovaleski et al., 2020) (Fig. 11B; Table 7). As predicted, the frequency of STN activity was significantly lower in Q175 mice (Fig. 11B,C; Table 7). The regularity of STN activity and the proportion of spikes generated during the inactive component of cortical SWA relative to those during the active component were similar in Q175 and WT mice (Fig. 11B–F; Table 7). Pinch-evoked cortical ACT elevated the frequency of STN activity in both WT and Q175 mice (Fig. 11G,H; Table 7). However, the frequency of STN activity remained lower in Q175 mice relative to that in WT mice during cortical ACT (Fig. 11G,H; Table 7).
STN neuron activity during cortical SWA and ACT in Q175 and WT micea
STN neurons are hypoactive in Q175 mice. A–D, During cortical SWA, the frequency (but not the regularity) of STN activity was lower in Q175 mice (A, schematic of experimental setup, illustrating placement of an optrode in the GPe and a tetrode array in the STN; B, representative examples of concurrent 0.1-1.5 Hz bandpass-filtered EEG and STN unit activity; C, D, population data). E, F, STN activity was similarly phase-locked to the active component of cortical SWA in Q175 and WT mice (E, population spike phase histograms; F, population data). G, H, Pinch-evoked cortical ACT (hind paw pinch, dotted line) increased STN neuron activity in both genotypes. However, during cortical ACT, STN neurons remained hypoactive in Q175 mice relative to WT (G, representative examples; H, population data). *p < 0.05.
To determine how PV+ GPe neurons regulate STN activity, the response of STN neurons to optogenetic inhibition of PV+ GPe neurons was recorded. The proportion of STN neurons that responded to inhibition of PV+ GPe neurons was similar in Q175 and WT mice (Table 8). Optogenetic inhibition of PV+ GPe neurons rapidly and persistently increased the frequency and regularity of STN activity in both WT and Q175 mice (Fig. 12A–E; Table 8). However, during optogenetic inhibition, STN neurons remained less active in Q175 mice during cortical SWA (Fig. 12A–D; Table 8) and ACT (Table 8). The ratio of anti- to in-phase STN activity was also significantly elevated by inhibition of PV+ GPe neurons in Q175 and WT mice (Fig. 12F,G; Table 8). Together, with data from this and previous studies, these findings argue that excessive inhibition arising from hyperactive PV+ GPe neurons and autonomous firing deficits contribute to the relative hypoactivity of STN neurons in Q175 mice.
Impact of optogenetic inhibition of PV+ GPe neurons on STN activity in Q175 and WT micea
STN hypoactivity is partially rescued by optogenetic inhibition of PV+ GPe neurons in Q175 mice. A, Representative examples of the EEG, bandpass-filtered at 0.1-1.5 Hz, and concurrent PV+ GPe and responsive STN neuron activity in WT and Q175 mice before, during, and after optogenetic inhibition (yellow bar) of PV+ GPe neurons. A–C, Optogenetic inhibition of PV+ GPe neuron activity rapidly and persistently disinhibited STN neuron activity. B, C, PSTHs of STN neuron activity in the absence and presence of PV+ GPe neuron optogenetic inhibition (yellow bar) (bin size, 100 ms; prestimulus mean, solid red line; ±2 SDs of prestimulus mean, dotted red line; B, representative examples from A; C, population data). A–E, Optogenetic inhibition of PV+ GPe neurons disinhibited and regularized STN activity in both genotypes. However, the relative hypoactivity of STN neurons in Q175 mice was not fully reversed by optogenetic inhibition of PV+ GPe neurons (D, E, population data). F, G, The ratio of anti- to in-phase STN activity was increased to a similar extent by inhibition of PV+ GPe neurons in WT and Q175 mice (F, population data; G, population spike phase histograms). *p < 0.05.
Discussion
Cortical network and layer V projection neuron activity
Although the cortex of symptomatic Q175 mice exhibits a variety of pathologic changes, including nuclear aggregates of mHTT (Smith et al., 2014), reduced expression of TrkB (Smith et al., 2014) and mGlu receptors (Bertoglio et al., 2018), a reduction in excitation:inhibition balance (Indersmitten et al., 2015), and cortical volume loss (Heikkinen et al., 2012), the EEG, and the firing rates and patterns of layer V PT and unidentified cortical pyramidal neurons were similar in WT and Q175 mice. The firing rate of unidentified cortical neurons is also unaltered in other HD models (Rebec, 2018).
Striatal projection neuron activity
Consistent with synaptic driving by cortical inputs, D2-SPNs and unidentified, putative D1-SPNs primarily discharged during the active component of cortical SWA. No difference in the firing of identified D2-SPNs and putative D1-SPNs was observed in WT mice during cortical SWA, consistent with some (Sharott et al., 2017) but not other (Mallet et al., 2005) studies that used a similar anesthetic regimen. Given that stimulation was not used to elicit firing in putative D1-SPNs, only neurons with spontaneous activity were detected, which may underlie the higher firing rate of D1-SPNs reported here versus the Mallet et al. (2005) study.
During cortical SWA, the frequency of D2-SPN activity in Q175 mice was half that in WT mice. In contrast, putative D1-SPN activity was similar in Q175 and WT mice. Given similar cortical activity in WT and Q175 mice, the hypoactivity of D2-SPNs in Q175 mice during SWA is consistent with reductions in cortico-striatal inputs, LTP, and strength, and dendritic hypoexcitability in D2-SPNs (Plotkin and Surmeier, 2015; Sebastianutto et al., 2017; Carrillo-Reid et al., 2019; but see Goodliffe et al., 2018). In contrast, the retention of these properties in D1-SPNs (Plotkin and Surmeier, 2015; Carrillo-Reid et al., 2019) may have contributed to their normal activity in Q175 mice. Consistent with impaired cortical driving of D2-SPNs in vivo, firing in response to cortical stimulation, burst firing, and correlated cortico-striatal activity are diminished in unidentified striatal neurons in HD mice (Beaumont et al., 2016; Rebec, 2018).
During cortical desynchronization, D2-SPN activity decreased to a lesser degree in Q175 than WT mice and unidentified/putative D1-SPN activity decreased in WT but not Q175 mice. As a result, SPN activity was elevated in Q175 mice relative to WT during ACT. While depolarization of the resting membrane potential and increased axosomatic excitability of SPNs in Q175 mice (Heikkinen et al., 2012; Indersmitten et al., 2015; Beaumont et al., 2016) do not compensate for reduced cortical driving of D2-SPNs during cortical SWA, they may contribute to elevated SPN activity during cortical desynchronization. Another possibility is that somatosensory stimulation used to trigger cortical desynchronization engages more powerful thalamo-striatal excitation in Q175 mice. Given that SPN activity is sculpted by striatal interneurons (Tepper et al., 2018), alterations in interneuron excitability and transmission (Holley et al., 2019a,b) may also contribute to aberrant SPN activity in HD mice. The spontaneous activity of unidentified striatal neurons has been reported to be elevated, reduced, or unchanged in anesthetized knock-in or awake full-length HD mice (Miller et al., 2008; Estrada-Sanchez et al., 2015; Beaumont et al., 2016). However, an absence of cell-identification approaches and/or different recording conditions in those studies make comparisons with this study problematic (Heikkinen et al., 2012; Indersmitten et al., 2015; Beaumont et al., 2016).
Prototypic PV+ GPe and putative PV– GPe neuron activity
Prototypic PV+ GPe neurons were hyperactive in Q175 mice during cortical SWA, consistent with upstream D2-SPN hypoactivity. However, additional factors must contribute because PV+ GPe neurons were also hyperactive in Q175 mice when (1) D2-SPNs were silent (in both Q175 and WT mice) during the inactive phase of cortical SWA; and (2) D2-SPNs were hyperactive in Q175 mice during cortical ACT. Furthermore, striatopallidal transmission strength increases in HD and its models through upregulation of postsynaptic GABAA receptor expression (Glass et al., 2000; Waldvogel et al., 2015; Perez-Rosello et al., 2019). Indeed, during cortical SWA, PV+ GPe neurons in Q175 mice more commonly exhibited pauses/reductions in firing that were concomitant with bouts of D2-SPN activity. Thus, PV+ GPe neuron hyperactivity may in part stem from elevated autonomous firing, as demonstrated here in Q175 mouse brain slices. Given the widespread nature and potency of prototypic GPe neuron projections in the basal ganglia (Bevan et al., 1998; Mallet et al., 2012; Mastro et al., 2014; Abdi et al., 2015), it will be critical to determine which alterations underlie their upregulated autonomous firing, and whether hyperactivity results from the cell-autonomous effects of mHTT and/or homeostatic compensation.
The hypoactivity of unidentified, putative PV– GPe neurons in Q175 mice in vivo appears to be caused by excessive inhibition because it was alleviated by optogenetic inhibition of hyperactive prototypic PV+ GPe neurons. The powerful regulation of putative PV– prototypic and arkypallidal neuron activity by prototypic PV+ GPe neurons is consistent with recent studies (Aristieta et al., 2020; Cui et al., 2021; Ketzef and Silberberg, 2021). However, optogenetic silencing of PV+ GPe neurons also disinhibited the STN, which feeds back onto GPe neurons (Aristieta et al., 2020; Pamukcu et al., 2020). Thus, STN disinhibition may have contributed to the rescue of putative PV– GPe neuron activity in Q175 mice. The finding that the autonomous firing of putative PV– GPe neurons was unaffected in Q175 mice further suggests that synaptic mechanisms underlie their hypoactivity in vivo. The low rate of putative PV– GPe neuron autonomous firing and longer duration of action potentials compared with PV+ GPe neurons in WT (and Q175) mice are consistent with previous observations in WT rodents (Abdi et al., 2015; Hernandez et al., 2015; Cui et al., 2021). Finally, our data contrast with another study that reported no change in GPe neuron activity in Q175 mice (Beaumont et al., 2016) and WT GPe activity that was substantially lower than for prototypic GPe neurons (Aristieta et al., 2020; Kovaleski et al., 2020). The identities of GPe neurons were not determined in the Beaumont (2016) study; and although a similar anesthetic was used, the depth of anesthesia was not tracked. Thus, their measurements of GPe activity were likely derived from a variety of GPe cell classes during SWA and ACT.
STN activity
In Q175 mice, STN activity was half that in WT mice during both cortical SWA and ACT. Thus, the hypoactivity of STN neurons observed here is consistent with that reported in Q175 mice in which the EEG was not measured and cortical activity state was unknown (Beaumont et al., 2016), and more robust than in YAC128 (Callahan and Abercrombie, 2015a) and R6/2 mice (Callahan and Abercrombie, 2015b), where only a subset of STN neurons exhibited diminished activity. The hypoactivity of STN neurons in Q175 mice could arise from hyperactivity of upstream prototypic PV+ GPe neurons, loss of autonomous STN activity (Atherton et al., 2016), and reduced cortical excitation (Beaumont et al., 2016). Indeed, STN hypoactivity was greatly (but not completely) alleviated in Q175 mice by optogenetic inhibition of prototypic PV+ GPe neurons, consistent with excessive inhibition from PV+ GPe neurons, but also other factors, including loss of intrinsic firing, (Atherton et al., 2016) and disinhibition of prototypic PV– GPe neurons that project to the STN (Mallet et al., 2012; Mastro et al., 2014; Hernandez et al., 2015; Aristieta et al., 2020; Cui et al., 2021). Phase-locking of STN activity to cortical SWA was unaffected in Q175 mice both before and during optogenetic inhibition of prototypic GPe neurons, arguing that excessive inhibition from prototypic GPe neurons and impaired intrinsic excitability, rather than reduced cortical drive, underlie STN hypoactivity. Because optogenetic inhibition of PV+ GPe neurons and disinhibition of the STN rescued putative PV– GPe neuron activity, STN hypoactivity may contribute to PV– GPe neuron hypoactivity in Q175 mice. In contrast, STN hypoactivity did not prevent hyperactivity in PV+ GPe neurons, presumably because their hyperactivity resulted from alterations in autonomous firing and striatopallidal transmission.
Functional implications
Together, our and other studies argue that the basal ganglia indirect pathway is profoundly dysregulated in HD mice through a complex, cell class-specific combination of state-dependent changes in presynaptic activity and aberrant intrinsic and synaptic properties. According to the classic rate model, PV+ GPe neuron hyperactivity and STN (and possibly arkypallidal) neuron hypoactivity should promote hyperkinesia (Albin et al., 1989). However, HD mice typically exhibit bradykinesia/akinesia and subtle kinematic abnormalities (Heikkinen et al., 2020). Therefore, consistent with their large number of trinucleotide repeats, HD mice may better represent juvenile-onset HD in which the initial hyperkinetic phase seen in adult-onset HD is absent (Fusilli et al., 2018; Tereshchenko et al., 2019). Because the indirect pathway is critical for normal volitional movement (Dodson et al., 2015; Tecuapetla et al., 2016; Pasquereau and Turner, 2017; Markowitz et al., 2018; Aristieta et al., 2020), dysregulation of this pathway may actively contribute to motor dysfunction in HD mice.
Footnotes
This work was supported by CHDI Foundation; National Institutes of Health, National Institute of Neurological Disorders and Stroke Grant R37 NS041280; and the Feinberg School of Medicine, Northwestern University. Confocal imaging was performed at the Northwestern University Center for Advanced Microscopy, which is supported by National Institutes of Health NCI Grant CCSG P30 CA060553 (awarded to the Robert H. Lurie Comprehensive Cancer Center). We thank Sasha Ulrich, Danielle Schowalter, and Bonnie Erjavec for the maintenance and supply of transgenic mice for this study. We also thank Drs. Vahri Beaumont, Ignacio Muñoz-Sanjuán, and Roger Cachope for their scientific input throughout this study.
The authors declare no competing financial interests.
- Correspondence should be addressed to Mark D. Bevan at m-bevan{at}northwestern.edu