Abstract
The motor cortex comprises the primary descending circuits for flexible control of voluntary movements and is critically involved in motor skill learning. Motor skill learning is impaired in patients with Parkinson's disease, but the precise mechanisms of motor control and skill learning are still not well understood. Here we have used transgenic mice, electrophysiology, in situ hybridization, and neural tract-tracing methods to target genetically defined cell types expressing D1 and D2 dopamine receptors in the motor cortex. We observed that putative D1 and D2 dopamine receptor-expressing neurons (D1+ and D2+, respectively) are organized in highly segregated, nonoverlapping populations. Moreover, based on ex vivo patch-clamp recordings, we showed that D1+ and D2+ cells have distinct morphological and electrophysiological properties. Finally, we observed that chemogenetic inhibition of D2+, but not D1+, neurons disrupts skilled forelimb reaching in adult mice. Overall, these results demonstrate that dopamine receptor-expressing cells in the motor cortex are highly segregated and play a specialized role in manual dexterity.
Significance Statement
The primary motor cortex (M1) is a command center governing voluntary motor control and performance of fine motor skills. Dopamine (DA) signaling in the M1 is required for skill learning and the underlying synaptic plasticity, although little is known about DA-recipient neurons in the M1. Here, we show that neurons in mouse M1 that express D1 and D2 receptors are arranged into distinct, nonoverlapping populations. Furthermore, we show that D2-expressing neurons in M1 control skilled forelimb reaching in adult mice. Overall, our findings reveal distinct circuits for DA receptor-expressing cells in the M1 as well as a unique function for D2-expressing neurons in manual dexterity.
Introduction
The brain's dopamine (DA) system is critical for movement control and motor learning, as evidenced by its role in Parkinson's disease (PD), a degenerative, neurological disorder characterized by motor deficits caused by progressive DA depletion (Ehringer and Hornykiewicz, 1960; Kish et al., 1988). Primary motor symptoms of PD are tremor, rigidity, and bradykinesia, but they are often accompanied by upper limb motor impairments, particularly in fine motor skills, including reaching and grasping (Jankovic, 2008; Quinn et al., 2013; Parma et al., 2014; Fasano et al., 2022). Movement abnormalities in PD arise from disordered neural activity in the corticobasal ganglia circuits resulting from the loss of DA signaling (Albin et al., 1989; DeLong, 1990; McGregor and Nelson, 2019). DA action is mediated by Gs-coupled D1 and Gi-coupled D2 receptors, and neurons expressing these receptors are highly segregated in the basal ganglia circuits (Gerfen et al., 1990; Bertran-Gonzalez, 2010; Wei et al., 2018), which suggests that they have specialized roles in the regulation of motor function. In line with this hypothesis, animal studies utilizing optogenetics and chemogenetics show that cell type-specific modifications of activity of D1 and D2 receptor-expressing neurons in the basal ganglia modulate specific parkinsonian motor behaviors (Kravitz et al., 2010; Alcacer et al., 2017).
Despite extensive research on the function of D1 and D2 receptor-expressing projection neurons in the basal ganglia, still little is known about the role of DA receptor-expressing neurons in cortical circuits. In this study, we focused on the primary motor cortex (M1), a command center governing voluntary motor control and performance of fine motor skills (Guo et al., 2015a). In rodents, this cortical region is strongly innervated by DA fibers and is enriched in both D1 and D2 receptors (Hosp et al., 2011; Kunori et al., 2014; Vitrac et al., 2014). Importantly, DA signaling in the M1 is necessary for skill learning and synaptic plasticity, as these processes are disrupted in animal models of PD (Molina-Luna et al., 2009; Guo et al., 2015b; Li et al., 2017; Aeed et al., 2021). Nevertheless, due to the complexity of M1 circuits, precise mechanisms of fine motor control are still not well understood. Therefore, using transgenic mice, electrophysiology, in situ hybridization, neural tract-tracing, and chemogenetics, we aimed to determine laminar organization of D1 and D2 receptor-expressing neurons in the M1, their properties, and their role in skilled forelimb reaching in adult mice.
Materials and Methods
Study approval
Procedures were approved by the 2nd Local Institutional Animal Care and Use Committee in Krakow (approval number 65/2019) and conducted in accordance with the directive 2010/63/EU of the European Parliament and of the Council of Sept. 22, 2010, on the protection of animals used for scientific purposes and with the Polish Act on the Protection of Animals Used for Scientific or Educational Purposes of Jan. 15, 2015.
Animals
Animals were housed 2–5 per cage in an animal facility room with a controlled temperature (22 ± 2°C) and humidity (40–60% RH), under a 12 h light/dark cycle. Unless otherwise specified, mice had ad libitum access to water and laboratory chow (RM1A, Special Diet Services). Drd1aCre mice (Lemberger et al., 2007) were obtained from the German Cancer Research Center, Heidelberg and Drd2Cre (Gong et al., 2007) from the University of California, Davis (MMRRC_032108-UCD). Drd1a-tdTomato line 6 (Ade et al., 2011) and Ai14 (tdTomato) Cre reporter line (Madisen et al., 2010) were purchased from The Jackson Laboratory (IMSR_JAX:016204 and IMSR_JAX:007914). All mice were congenic with the C57BL/6N background (>8 generations of backcrosses prior to initiation of the study). For the purpose of the project, the Drd2Cre strain was crossed with Drd1a-tdTomato and Ai14 (tdTomato) strains to obtain double transgenic animals. Genotyping was performed using a standard PCR assay according to previously described protocols and genotyping protocols available in the JAX database. The age of animals at the onset of the experiments was 8–12 weeks, with the exception of the patch-clamp experiments where 6–8-week-old animals were used. Each experimental group consisted of mice from at least two litters. Both sexes were used.
Stereotaxic injections
Animals were anesthetized with a mixture of ketamine (100 mg/ml) and xylazine (20 mg/ml) and positioned in a stereotaxic frame (ASI Instruments). Body temperature was maintained at 37°C throughout the procedure. A small craniotomy was made above the caudal forelimb area of the motor cortex using coordinates relative to Bregma: 0.3 ± 0.2 mm anterior, 1.5 ± 0.2 mm lateral, 1.5 ± 0.2 mm ventral. Borosilicate glass pipettes with 30–40 µm tip diameters were backfilled with oil and Cre-dependent viral vectors and connected to a microliter glass Hamilton syringe via Tygon tubing. A 200 nl of solution was pressure injected in batches of 20 nl, with 60 s intervals between injections. The pipette was left in the tissue for 5 min to ensure effective diffusion and was then slowly withdrawn. After the surgery, subcutaneous injections of anti-inflammatory drug Tolfedine (4%, Vetoquinol) and 5% glucose solution were given to alleviate pain and prevent dehydration. Similar procedure was followed in case of striatal injections, where Cre-dependent hM4Di viral vector was injected into the dorsal striatum using coordinates relative to Bregma: 0.1 ± 0.1 mm anterior, 2.0 ± 0.1 mm lateral, 3.0 ± 0.2 mm ventral (site 1); 0.4 ± 0.1 mm anterior, 2.0 ± 0.1 mm lateral, 3.0 ± 0.2 mm ventral (site 2); 0.7 ± 0.1 mm anterior, 2.0 ± 0.1 mm lateral, 3.0 ± 0.2 mm ventral (site 3). For the purpose of tracing experiments, Drd1aCre and Drd2Cre mice were injected with rAAV2-Ef1a-DIO-mCherry, and Drd2Cre::Drd1a-tdTomato mice were injected with rAAV2-Ef1a-DIO-EYFP. For behavioral experiments and multielectrode array (MEA) recordings, Drd1aCre mice, Drd2Cre mice, and their wild-type littermates received injections of pAAV-hSyn-DIO-hM4D(Gi)-mCherry. Viral vectors were obtained from Addgene and UNC Vector Core. Animals were returned to their home cages for 2 weeks before being used in behavioral and MEA experiments or 4 weeks in case of anatomy and tracing studies.
Histology and fluorescence microscopy
Mice were killed via intraperitoneal injection of Morbital (Biowet; sodium pentobarbital 133.3 mg/ml + pentobarbital 26.7 mg/ml) and perfused with ice-cold phosphate-buffered saline (PBS), pH 7.4, followed by 4% formaldehyde, pH 7.4. Dissected brains were fixed in 4% formaldehyde for 12 h at 4°C. Coronal slices (50 µm) containing the motor cortex and striatum, thalamus, midbrain, pons, and medulla were cut on a VT1000S vibrating blade microtome (Leica Biosystems). In the case of Drd1a-tdTomato and Drd2Cre::Drd1a-tdTomato mice, the tdTomato signal was enhanced with antibody against red fluorescent protein (RFP, Rockland). For tdTomato labeling, slices were first blocked in PBS-T containing 10% normal donkey serum (NDS, Jackson ImmunoResearch) and 0.6% Triton X-100 (Sigma-Aldrich) for 1 h at room temperature; followed by incubation with primary antibody anti-RFP 1:1,000 (Rockland) in PBS-T (2% NDS, 0.3% Triton X-100) overnight at 4°C; and secondary antibody Cy3 1:400 (Jackson ImmunoResearch) overnight at 4°C. Finally, rinsed slices were mounted on glass slides and coverslipped with Fluoroshield containing DAPI (Sigma-Aldrich). Slices obtained from Drd1aCre, Drd2Cre, or Drd2Cre::Ai14 mice were mounted on glass slides and coverslipped. Images were acquired using Axio Imager.M2 fluorescence microscope (Zeiss) equipped with Axiocam 503 mono camera. Whole brain images were acquired with 5× objective, and z-stack images of regions of interest were acquired with 20× objective. ZEN software (Zeiss) was used for initial image acquisition and preprocessing. ImageJ software (Schneider et al., 2012) was used for subsequent image processing that involved cropping and adjustment of brightness and contrast. Cell counting was performed in a ∼700 µm wide region of interest across the total depth of M1, using the ImageJ Cell Counter plugin. Landmarks indicating the borders between layers were identified based on the distance from the pia to the white matter, as described previously (Shepherd, 2009; Hooks et al., 2011).
Multiplex fluorescent in situ hybridization (RNAscope)
In situ hybridization using the RNAscope HiPlex12 v2 Assay for AF488, Atto550, and Atto647 detection (Advanced Cell Diagnostics, ACD) was performed on fresh frozen 16 µm brain sections obtained from Drd1a-tdTomato and Drd2Cre::Ai14 mice as described previously (Szlaga et al., 2022). During the primary hybridization step, the following probes (ACD) were used for round 1 (R1); Drd1R (Mm-Drd1-T1, catalog #461901-T1), tdTomato (tdTomato-T2, catalog #317041-T2), Drd2 (Mm-Drd2-T3, catalog #406501-T3), and round 2 (R2); vGlut1 (Mm-Slc17a7-T4, catalog #416631-T4), vGAT1 (Mm-Slc32a1-T6, catalog #319191-T6). Prior to hybridization with Fluoro T1-T3 and T4-T6, the FEPE Reagent was applied to minimize tissue autofluorescence. Images of the same M1 area for R1 and R2 were acquired with an Axio Imager M2 fluorescent microscope (Zeiss) with an automatic stage and Axiocam 503 mono camera (Zeiss) and processed using Zen (Zeiss), CorelDraw 2020 (Corel Corporation), ImageJ (Schneider et al., 2012), and HiPlex Image Registration Software v1.0 (ACD). Cells expressing the mRNAs tested were counted, using the ImageJ Cell Counter plugin, in three regions of interest: one in layer II/III and two in layers V–VI. The presence of tdTomato mRNA was used as a marker for neurons expressing either Drd1a or Drd2 receptors in Drd1a-tdTomato and Drd2Cre::Ai14 mice, respectively. Cells were identified by the presence of a cell-like distribution of fluorescent mRNA dots and/or a DAPI-stained nucleus. mRNA was assessed as present if at least two unambiguous dots were observed. The area-proportional Euler diagrams representing different cell groups and their relationships in specific layers of the cortex were generated using the open-source software Edeap (Wybrow et al., 2021).
Whole-cell patch-clamp recordings
Animals were deeply anesthetized with isoflurane (Baxter) and decapitated. Brains were collected in the ice-cold artificial cerebrospinal fluid (ASCF) containing the following (in mM): 92 NaCl, 2.5 KCl, 1.25 NaH2PO4, 25 NaHCO3, 20 HEPES, 5 Na+ ascorbate, 3 Na+ pyruvate, 2 thiourea, 10 glucose, 10 MgSO4, 0.5 CaCl2, pH 7.4 (osmolarity 290–300 mOsm/kg), saturated with carbogen (95% O2 and 5% CO2). Brains were cut into 200-µm-thick coronal slices containing M1 (AP: 0.00–0.60 mm from Bregma) on VT1000S vibrating blade microtome (Leica). Slices were immediately transferred to the incubation chamber containing carbogenated, warm (32°C) ACSF containing the following (in mM): 118 NaCl, 25 NaHCO3, 3 KCl, 1.2 NaH2PO4, 2 CaCl2, 1.3 MgSO4, and 10 glucose. After a recovery period (90–120 min, room temperature), slices were placed in a recording chamber, where the tissue was perfused (2 ml/min) with carbogenated, warm (32°C) ACSF of the same composition. Borosilicate glass pipettes (Sutter Instrument), 6–8 Ω tip resistance, were filled with a solution containing the following (in mM): 145 potassium gluconate, 2 MgCl2, 4 Na2ATP, 0.4 Na3GTP, 1 EGTA, 10 HEPES, pH 7.3 (osmolarity: 290–300 mOsm/kg), and 0.05% biocytin for subsequent immunofluorescent identification. All reagents were purchased from Sigma-Aldrich, except biocytin (Tocris Bioscience). The calculated liquid junction potential was +15 mV and analyzed data were corrected for this value.
Whole-cell current- and voltage-clamp recordings were obtained from pyramidal cells located in the motor cortex deep output layer V. Cells were identified with the Axio Examiner.A1 microscope (Zeiss) using video-enhanced infrared-differential interference contrast. The neuronal identity of dopaminoceptive neurons was confirmed by the presence of tdTomato excited at 530 nm with a Colibri LED light fluorescent system (Zeiss). Signal recordings were made with a SEC-05X amplifier (NPI Electronic) and Micro3 1401 converter (Cambridge Electronic Design, CED). The recorded signal was low-pass filtered at 3 kHz and digitized at 20 kHz. Electrophysiological properties of examined neurons were assessed based on their responses to a series of stimulation protocols in voltage- and current-clamp mode. To measure the I–V relationships of the steady-state current, voltage steps stimulations ranging from −120 to −50 mV (10 mV change, pulse duration 500 ms) were delivered from a holding potential of −75 mV. Neuronal excitability (input–output relationship) was measured as the number of spikes elicited by applied current pulses from −250 to +500 pA (50 pA increment, pulse duration 500 ms). Linear regression for excitability and steady-state current I–V relationship was calculated in GraphPad Prism software. The voltage sag was measured from the voltage response to a −250 pA hyperpolarizing current step (500 ms). The rheobase (the current value at the moment of reaching the AP threshold) was determined based on the voltage response to current ramp stimulation (0–1 nA, 1 s). Passive membrane properties (resistance, capacitance, time constant) were calculated from the voltage responses to a −50 pA current pulse (500 ms). Action potential (AP) properties were calculated from single AP evoked from a membrane potential of −75 mV, by a 0.5 ms depolarizing current pulse. Custom-written scripts in Signal and Spike2 software (CED) were used for data analysis.
Morphological reconstruction
Slices from patch-clamp experiments were fixed overnight at 4°C in 4% formaldehyde in PBS; blocked for 24 h at 4°C in PBS containing 10% NDS (Jackson ImmunoResearch) and 0.6% Triton X-100 (Sigma-Aldrich); followed by incubation with ExtrAvidin-Cy3 1:200 (Sigma-Aldrich) in PBS (2% NDS, 0.3% Triton X-100) for another 48 h at 4°C. Rinsed slices were mounted on glass slides and coverslipped with Fluoroshield containing DAPI (Sigma-Aldrich). To compare the dendritic morphology of recorded neurons, cells well-filled with biocytin that had clearly visible dendritic trees and no major truncations were further imaged using a LSM 780 META laser scanning confocal microscope on AxioObserver Z1 (Zeiss) under 20× magnification, and additional z-stack images were taken, using LSM 710 META on AxioObserver Z1 (Zeiss), under 40× magnification. Dendrite tracing and 3D reconstruction (10 µm steps) were performed in ImageJ using the Simple Neurite Tracer plugin (Longair et al., 2011). The width of the apical dendritic shaft was measured 5 µm above the soma, but due to the truncation of apical shafts in the majority of imaged cells, only the basal part of the dendritic tree was traced. L-Measure software (Scorcioni et al., 2008) was further used, to acquire morphological parameters: total dendritic length, maximal branch order, number of primary dendrites, branches, bifurcations, and dendritic tips.
MEA recordings
The multielectrode extracellular ex vivo experiments were performed as previously described (Chrobok et al., 2021). Two weeks after the viral vector transfection (described above) mice were decapitated, and the slices were prepared as for the patch-clamp recordings. The brains were collected in carbogenated, ice-cold ACSF, comprising the following (in mM): 25 NaHCO3, 3 KCl, 1.4 Na2HPO4, 2 CaCl2, 10 MgCl2·6H2O, 10 glucose, 185 sucrose; and the 250-µm-thick coronal slices were cut using a vibrating blade microtome (Leica). Sections containing the striatum were transferred to an incubation chamber filled with carbogenated, preheated (32°C) ACSF, with the same composition as for the patch-clamp recordings, with the addition of 0.01 g L−1 phenol red. The slices were incubated for 90 min and transferred to the recording wells of the MEA2100-Systems (Multi Channel Systems). The striatum was positioned upon the recording electrodes of the 8 × 8 perforated MEA (200 µm spacing, 60PedotpMEA200-30iR-Au, Multi Channel Systems). Carbogenated, warm (32°C) ACSF was used for continuous perfusion of the slices (2 ml/min). Before the recording, slices were left to settle for an hour. The clozapine N-oxide (CNO; Sigma-Aldrich) was freshly diluted in the recording ACSF (10 µM; 10 ml) and delivered by bath perfusion. The raw signal was sampled at 20 kHz and recorded using Multi Channel Experimenter (Multi Channel Systems). The data files were converted to HDF5 and CED-64 formats using Multi Channel Data Manager (Multi Channel Systems). Using custom-written Matlab scripts and the KiloSort algorithm (Pachitariu et al., 2016), the initial automatic spike sorting was performed. Then, the signal was transferred into the bandpass filtered (Butterworth filter; fourth order; cut off: 0.3–7.5 kHz) CED-64 file, and single units were manually refined using principal components analysis and autocorrelation. Afterward, the data were binned (30 s) and analyzed using NeuroExplorer 5 and the temporal heatmaps encoding the single-unit activity were prepared using a custom-written Matlab script to illustrate the drug-induced changes. A response to CNO application was considered significant if it differed by one standard deviation from the baseline mean of the recorded signal.
Motor skill training
To evaluate the function of cortical DA receptor-expressing neurons, a single pellet-reaching task was performed. The design of the Plexiglas training chamber and the experimental protocol was based on the previously described protocol (Farr and Whishaw, 2002), with modifications. To ensure motivation during training, animals were food deprived to ∼85% of their initial body weight. To habituate mice to food reward, food pellets (20 mg, dustless precision pellets, Bioserv) were given in the home cage 3 d before training. To habituate mice to the testing environment, a day before training each mouse was individually placed for 20 min in the Plexiglas testing chamber (20 cm long, 8 cm wide, 20 cm high, with a 0.5 cm wide vertical slit on the front narrow wall) and was allowed to explore the apparatus, consume pellets laying on the floor, or reach for multiple pellets located outside the box. During the subsequent training (days 1–7), mice were trained to extend their right forelimb (contralateral to hM4Di injected hemisphere) through a narrow slit to grasp and retrieve a single food pellet located on an elevated platform (1.5 cm high, 1 cm away from the opening, and centered 0.5 cm to the left). Following this initial training, mice performance was assessed for additional three sessions (days 8–10) during which baseline success rate was measured (pretest). On the subsequent test days (days 11–13), animals received intraperitoneal injections of CNO (2 mg/kg) 30 min before being placed in the testing apparatus. Reaching events during pretest and CNO treatment sessions were recorded at 100 frames-per-second and 640 × 480 pixels by a color video camera (Basler acA1440-220uc). Each session consisted of 40 trials that started when the animal approached the pellet and the reach attempt was recorded for 5 s. The success rate was calculated as a percentage of trials in which the pellet was successfully retrieved (regardless of the number of reaching attempts). Animals that failed to maintain at least 30% of success rate during the pretest were excluded from the analysis. Individual reach outcomes (from all reaching attempts) were analyzed and reach failures were classified as “no-grab” (pellet was touched but not correctly grasped), “miss” (pellet was not touched), or “drop” (pellet was retrieved and dropped before being placed in the mouth). Attempts to reach with the tongue or left forelimb (ipsilateral to hM4Di injected hemisphere) were omitted from the analysis. To account for day-to-day variability in behavior, data from the respective pretest and CNO treatment test days were averaged.
Kinematic analysis
To quantify reaching kinematics, the first reach attempt (irrespective of the outcome) extracted from trials recorded in the pretest and test sessions was analyzed. The position of the paw and the pellet were tracked with DeepLabCut (Mathis et al., 2018) from the moment the paw was lifted from the ground until it touched the pellet. Trials in which the starting position of the paw was not on the ground or in which the entire reaching motion was not clearly visible were excluded from analysis. In a small fraction of frames, the deep neural network failed to correctly recognize the position of the paw. In such cases, paw coordinates were linearly extrapolated. Paw distance to pellet was calculated as sqrt(x2 + y2) normalized to the number of pixels by 1 mm, and velocity was calculated as a derivative of this distance.
Open field
To assess the effects of CNO on basic locomotor activity, animals were tested in the open field (35.5 cm long × 25.5 cm wide × 22 cm high). Animals received an intraperitoneal injection of CNO (2 mg/kg) 30 min prior to the placement in the arena filled with the home cage bedding. A single 10 min session was recorded at 60 frames-per-second and 640 × 480 pixels by a monochrome video camera (Creative HD 720p). The position of the animal was tracked with DeepLabCut. The total distance traveled and average movement speed were measured.
Experimental design and statistical analyses
Results are presented as mean ± SEM. Statistical analysis was based on the assumption that the samples follow a Gaussian distribution. Student's t test was applied for statistical comparisons between two groups and one-way or two-way ANOVA followed by post hoc analysis (Bonferroni's multiple-comparisons test or Fisher's least significant difference test) was used for analysis with multiple groups, and repeated measures were incorporated when appropriate. p < 0.05 was considered statistically significant. All statistical analyses were conducted using GraphPad Prism software.
Data availability
Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.
Results
Primary motor cortex D1+ and D2+ neurons are organized in differently distributed, highly separated populations
To determine the distribution of putative D1 and D2 receptor-expressing neurons (D1+ and D2+, respectively) throughout the layers of the motor cortex, we used DA receptor-specific reporter lines: Drd1aCre mice injected with DIO-mCherry virus and double transgenic Drd2Cre::Ai14 mice (Fig. 1A). We counted fluorescently labeled cells throughout the superficial (I, II/III) and deep output (V–VI) layers. We found that D1+ cells were mainly distributed throughout the deep output layers V–VI, while D2+ positive neurons were primarily distributed in superficial layer II/III (Fig. 1B,C; Table 1). To further verify if M1 neurons coexpress D1 and D2 receptors, we used double transgenic Drd1a-tdTomato::Drd2Cre mice injected with DIO-EYFP virus (Fig. 1D). We counted fluorescently labeled D1+ (tdTomato-positive) and D2+ (EYFP-positive) cells and found that the majority of studied neurons were positive for only one fluorescent reporter, with only 4.15% positive for both tdTomato and EYFP (Fig. 1E). Overall, these results indicate that D1+ and D2+ neurons in M1 are organized in differently distributed, highly separated populations.
Distinct neurochemical profiles define the layer-specific distribution of D1+ and D2+ M1 neurons
Cortical neurons can be broadly classified into glutamatergic pyramidal cells or GABAergic interneurons (BRAIN Initiative Cell Census Network (BICCN), 2021). To determine whether D1+ and D2+ neurons in M1 belong to one or the other group, RNAscope multiplex in situ hybridization was performed, using brain sections from Drd1a-tdTomato and Drd2Cre::Ai14 mice (Fig. 2A) and specific probes for tdTomato, vesicular glutamate (vGlut1) and GABA (vGAT1) transporters mRNA, as well as D1 and D2 receptors mRNA (Fig. 2B). In both mouse lines, the observed tdTomato mRNA expression pattern corresponded to the distribution of tdTomato-immunolabeled D1+ and D2+ neurons described in the previous section (Fig. 1A). Furthermore, the majority of tdTomato mRNA-expressing (tdTomato+) cells in slices from Drd1a-tdTomato mice were located within deep cortical layers V–VI, while tdTomato+ neurons in slices from Drd2Cre::Ai14 mice were primarily present in layer II/III (Fig. 2A, Table 2). Most of tdTomato+ cells in both Drd1a-tdTomato and Drd2Cre::Ai14 mice coexpressed vGlut1 mRNA (vGlut1+), indicating that they are primarily glutamatergic pyramidal neurons, with the remaining smaller proportion being vGAT1 mRNA-expressing (vGAT1+) GABAergic neurons, and a marginal fraction of cells coexpressing mRNA for both transporters (vGlut1+/vGAT1+; Fig. 2C, Table 1). Interestingly, the observed neurochemical phenotypes of tdTomato+ neurons in Drd1a-tdTomato and Drd2Cre::Ai14 animals depended on their laminar distribution (Fig. 2D, Table 1). In deep cortical layers V–VI (Fig. 2D, right; Table 1), most tdTomato+ neurons in Drd1a-tdTomato mice were vGlut1+ and constituted a vast majority of all vGlut1+ neurons identified in this area, while the remaining tdTomato+ cells belonged to a smaller vGAT1+ population. In contrast, the fraction of tdTomato+ neurons coexpressing vGlut1 and vGAT1 mRNA in Drd2Cre::Ai14 mice in the deep layers was almost equal. This pattern was reversed in the superficial layer II/III (Fig. 2D, left; Table 1), where populations of tdTomato+ cells coexpressing vGlut1 and vGAT1 mRNA in Drd1a-tdTomato mice were more evenly distributed, while the majority of tdTomato+ cells in Drd2Cre::Ai14 mice coexpressed vGlut1 mRNA. Consistent with the results shown in Figure 1D and E, we observed minimal Drd2 and Drd1a mRNA colocalization in tdTomato+ cells in Drd1a-tdTomato and Drd2Cre::Ai14 mice, respectively (Table 2). These results provide further evidence supporting the segregation of D1+ and D2+ neurons in M1 into distinct subpopulations, with D1+ cells located mainly in the deep layers and displaying predominantly glutamatergic phenotype and D2+ neurons being GABAergic, while the opposite is true in the superficial layers.
D1+ and D2+ pyramidal neurons have distinct electrophysiological and morphological properties
Previous studies have shown that D1+ and D2+ layer V pyramidal cell subclasses in the prefrontal cortex (PFC) have unique electrophysiological and morphological properties (Gee et al., 2012; Seong and Carter, 2012; Clarkson et al., 2017). To test if similar diversity applies to the motor cortex, whole-cell patch-clamp recordings and subsequent morphological analysis were performed from fluorescently labeled D1+ and D2+ layer V pyramidal neurons in brain slices obtained from Drd1a-tdTomato and Drd2Cre::Ai14 mice, respectively (Fig. 3A). To test for possible differences in the voltage dependence of the current flowing through the membrane of D1+ and D2+ neurons, their steady-state current–voltage (I–V) relationship was characterized and compared using the linear regression curve slope (Fig. 3C, Table 1). The performed analysis showed that examined slopes differed significantly and subsequent post hoc tests revealed differences between tested groups at low voltage step values, indicating differences in the types or properties of ion channels expressed by D1+ and D2+ neurons. We did not observe differences between D1+ and D2+ neurons in their resting membrane potential or passive membrane properties (resistance, capacitance, time constant; Fig. 3D, Table 1). To characterize and compare the excitability of D1+ and D2+ neurons, the relationship between the injected current and recorded firing rate was evaluated using the linear regression curve slope (Fig. 3E, Table 1). The performed analysis revealed that D2+ neuron excitability curve slope was significantly steeper than D1+ neurons, indicating higher excitability (stronger neuronal gain) of D2+ cells.
Previous studies showed that D1+ and D2+ PFC pyramidal neurons can be identified based on the characteristics of the voltage sag induced by hyperpolarization-activated cation currents (Ih currents). Specifically, D1+ PFC neurons were characterized by the absence of voltage sag in response to hyperpolarizing current (Seong and Carter, 2012), while D2+ neurons were characterized by a large sag amplitude (Gee et al., 2012; Clarkson et al., 2017). In line with these findings, we observed significantly larger voltage sag amplitude in D2+ than that in in D1+ M1 pyramidal neurons (Fig. 3D, Table 1). Moreover, analysis of the shape of evoked APs revealed that D2+ M1 pyramidal neurons have smaller AP half width and lower value of afterhyperpolarization (AHP) trough when compared with D1+ cells (Fig. 3F,G; Table 1). Importantly, the observed higher sag amplitude and shorter AP duration in D2+ neurons may contribute to greater excitability of D2+ relative to D1+ neurons (Hogan and Poroli, 2008). At the same time, ramp current stimulations revealed that rheobase (the minimum current necessary to elicit an AP) is higher in D2+ compared with those in D1+ neurons (Fig. 3D, Table 1). Given the lack of differences in the AP threshold (Fig. 3G, Table 1), a higher rheobase value in D2+ neurons may arise from the summation of small (below significance threshold) differences in membrane parameters influencing the rheobase, such as membrane potential, time constant, resistance, and/or capacitance.
During patch-clamp recordings, cells were filled with biocytin, allowing post hoc morphological reconstruction of their dendritic trees (Fig. 4A,B). Previous data suggested possible differences in the width of the apical dendrite shafts (measured 5 µm above the soma) between D1+ and D2+ in the PFC (Gee et al., 2012); therefore, we performed such measurements and analysis but have not found similar differences between D1+ and D2+ M1 pyramidal neurons (Fig. 4C,D; Table 1). Due to the relatively small slice thickness (∼200 µm, which was required to detect tdTomato fluorescence with our recording setup) and resulting partial truncation of the apical dendrite tufts, only the basal dendritic trees were compared. The Sholl intersection profile revealed that D2+ neurons have significantly more complex basal dendritic trees in the distal parts of the basal dendrites (Fig. 4E, Table 1). We also found that D2+ neurons, as compared with D1+ cells, have greater total dendritic length, a larger number of primary dendrites, and a larger number of dendritic tips (Fig. 4F, Table 1). Taken together, it was revealed that layer V D2+ pyramidal neurons are characterized by increased intrinsic excitability and higher dendritic tree complexity.
D1+ and D2+ motor cortex neurons have distinct projection targets across the brain
Motor cortex pyramidal neurons comprise three broad classes of output projection neurons: (1) intratelencephalic (IT), targeting cortical and striatal regions with somas distributed across layers II–VI; (2) pyramidal tract (PT), innervating the brainstem and spinal cord with somas located in layer V; and (3) corticothalamic (CT), projecting to the thalamus with somas in layer VI (Muñoz-Castañeda et al., 2021). We used anterograde tracing to trace axonal projections of D1+ and D2+ neurons in Drd1aCre and Drd2Cre mice and observed that both D1+ and D2+ neurons originating in M1 have extensive projections to the cortical and subcortical regions of the brain. Their fibers are present in contralateral M1 (Fig. 5A), ipsilateral somatosensory cortex (Fig. 5B), cerebral peduncle (Fig. 5C), pontine nuclei (Fig. 5D), medullary pyramids (Fig. 5E), and their decussation (Fig. 5F). In addition, we found that D1+ neurons innervate the dorsolateral striatum (Fig. 5A) and thalamic nuclei (Fig. 5B). These results show that D1+ cells constitute the majority of output projection neurons from the motor cortex and that this is a heterogeneous group that overlaps with the three main classes of output neurons. In turn, D2+ neurons are a relatively smaller fraction of output projection neurons that primarily constitute fibers of corticocortical and to a lesser extent corticopontine and corticospinal tracts. Overall, the output patterns of D1+ and D2+ neurons indicate that both classes of projection neurons are well positioned to mediate the control of skilled motor behavior.
Inhibition of D2+ neurons disrupts skilled forelimb grasping
The motor cortex is involved in the control of skilled forelimb movements, particularly reaching and grasping (Guo et al., 2015a). To test the functional role of newly identified classes of D1+ and D2+ neurons in controlling the execution of skilled forelimb movements, we used a chemogenetic approach (Roth, 2016) and the pellet-reaching task (Farr and Whishaw, 2002; Fig. 6A,B). Drd1aCre mice, Drd2Cre mice, and their wild-type (WT) littermates were injected with the DIO-hM4Di virus, inducing Cre-dependent expression of Gi-coupled designer receptors exclusively activated by designer drugs (DREADDs), in the left M1, contralateral to the reaching forelimb (Fig. 6C). This allowed for cell type-specific inhibition of neural activity by injections of DREADDs ligand, CNO (Extended Data Fig. 6-1).Mice were trained to extend their right forelimb through a narrow slit to grasp and retrieve a single food pellet located on an elevated platform. Once their performance was stable across 3 consecutive days of the pretest, mice were intraperitoneally injected with CNO (2 mg/kg). While the performance of Drd1aCre mice remained intact, we observed a reduction of successful reaches after the activation of inhibitory DREADDs in Drd2Cre mice (Fig. 6D, Table 1). A detailed analysis of individual reach outcomes revealed that the CNO-induced reduction of success rate in Drd2Cre mice was associated with degraded performance during the grasp phase, as evidenced by increased “no-grab” incidence (Fig. 6E, Table 1). We further tested whether chemogenetic inhibition of D1+ or D2+ M1 neurons influences reach kinematics, but no significant effect on movement distance or speed was observed in any of the groups tested (Fig. 6F, Table 1). In addition, no effects of CNO treatment on general locomotor activity were observed, when hM4Di injected mice were tested in the open field (Fig. 6G, Table 1). Overall, these results show that D2+ neurons in M1 play a specialized role in manual dexterity, rather than gross forelimb kinematics or locomotor behavior.
Figure 6-1
Validation of chemogenetic inhibition of neuronal activity in Drd1aCre and Drd2Cre mice using striatal neurons as an example. (A) Coronal section obtained from a wild-type (left) and Drd1aCre mouse (right) showing slice placement upon the multi-electrode array (black dots indicate single recording electrodes) and hM4Di-mCherry expression in the striatum after recording. Scale bars, 200 μm. (B) Table summarizing recorded neurons and their responses to CNO application. (C) Temporal heatmaps encoding single-unit activity (SUA) of all recorded neurons and their response to CNO application (10 μM, 10 ml; indicated by a black line). Each row indicates a single unit. White dotted horizontal lines classify units to either Drd1aCre mice (n = 46 units), Drd2Cre mice (n = 34), or wild-type mice (n = 35). The cell activity was normalized from 0 to 1 and was sorted by the defined window corresponding to the CNO response. Bin, 30 s. (D) Examples of neurons recorded during CNO administration in Drd1aCre, Drd2Cre and wild-type mice. Top panels show separated spikes of a single neuron, bottom panels show a corresponding frequency histogram. The orange-shaded rectangle indicates the duration of CNO action. Bin, 60 s. (B-D) Slices were obtained from n = 1 hM4Di-mCherry injected animal for each group, n = 4 slices per animal. Download Figure 6-1, TIF file.
Discussion
Here, we have identified two differentially distributed, largely nonoverlapping populations of putative DA receptor-expressing neurons in the primary motor cortex. Specifically, we observed that neurons located in the superficial layer II/III predominantly expressed D2 receptors, while the majority of the cells located in the deep output layers V and VI expressed D1 receptors. Our observation corresponds with previous studies demonstrating the segregated expression of D1 and D2 receptors in striatal and prefrontal neurons (Gerfen et al., 1990; Bertran-Gonzalez, 2010; Wei et al., 2018), suggesting that this segregation is well conserved throughout the central nervous system and can reflect the existence of distinct populations of neurons with different functional roles. The majority of DA receptor-expressing cells in M1 were found to be glutamatergic pyramidal neurons with extensive corticocortical and subcortical projections. Interestingly, within the same output layer V, we observed significant differences in the morphological and electrophysiological properties between D2+ and D1+ pyramidal cells, with the former exhibiting increased intrinsic excitability and greater complexity of the dendritic tree. These findings are in line with previous reports showing differences in the morphology and excitability of DA receptor-expressing pyramidal cells located in the prefrontal and motor cortices (Gee et al., 2012; Seong and Carter, 2012; Clarkson et al., 2017; Swanson et al., 2021). Collectively, these reports and our present results strongly suggest that the length and complexity of the dendritic tree are unique features that distinguish between layer V D1+ and D2+ pyramidal neurons, across prefrontal and motor areas of the cerebral cortex.
Previous studies have shown that motor skill learning is accompanied by functional and structural reorganization in the upper layers of the M1 forelimb region, contralateral to the trained forelimb (Rioult-Pedotti et al., 1998, 2000; Harms et al., 2008). This reorganization process has been shown to be dependent on DA signaling (Molina-Luna et al., 2009; Rioult-Pedotti et al., 2015) and to coincide with the emergence of population neuronal activity during learning (Peters et al., 2014). Our study revealed that inhibition of D2+ neurons in the motor cortex disrupts skilled forelimb reaching in adult mice that have learned the task. This observation aligns with previous reports demonstrating that blocking D2 receptors impairs skill learning and reduces plasticity within M1 (Rioult-Pedotti et al., 2015). As previously discussed, D2+ neurons exhibit increased intrinsic excitability, a hallmark of the formation of neuronal ensembles (Alejandre-García et al., 2022). Therefore, D2+ neurons may be predisposed for recruitment during skill learning and subsequent reactivation during learned movement execution. Nevertheless, the selectivity of this effect was surprising, as the D1+ neurons constitute a strong representation of layer V PT neurons that also undergo rapid plastic changes during motor skill learning (Xu et al., 2009) and can directly mediate spinal circuits for skilled reaching and grasping (Alstermark and Isa, 2012). However, a recent study showed that inhibition of PT neurons does not affect reaching and grasping in mice, while inhibition of IT neurons severely disrupts performance of the skilled reach-to-grasp task (Park et al., 2022). In line with this, D2+ neurons largely comprise the IT-type subclass of neurons located in superficial layer II/III. Considering the top-down organization of M1 (Weiler et al., 2008; Anderson et al., 2010), where layer II/III neurons provide output to neurons located in lower layers, the role of D2+ neurons in motor performance may be of greater significance. Moreover, recent studies show that layer II/III neurons exhibit more prevalent direction-selective activity than layer V neurons (Galiñanes et al., 2018; Currie et al., 2022) and encode previous reach outcomes independent of kinematics and reward (Levy et al., 2020), which suggests that they play a superior role in motor skill learning and subsequent movement execution.
Taken together, our findings have important implications for understanding the pathophysiology of movement disorders, in particular PD. The M1 is the most important structure for volitional motor control, and DA depletion influences the firing rate and synchronization of motor cortex neurons (Underwood and Parr-Brownlie, 2021). Although we did not directly investigate the mechanisms underlying the effect of experimental PD or the intricate effects of DA receptor activation on cortical neural activity, we did observe that inhibiting D2 receptor-expressing neurons in M1 results in impaired motor performance in a forelimb reaching task. This observation allows us to speculate that this neuronal population can contribute to upper limb motor impairments observed in PD patients. Experimental data showing that D2 receptor blockade-induced bradykinesia is associated with reduced motor cortex activity (Parr-Brownlie, 2005), as well as reduced neuroplasticity in the human M1 (Monte-Silva et al., 2011), support this hypothesis. Moreover, as a growing body of evidence suggests, D2 receptor activation increases the excitability of pyramidal neurons in PFC and M1 (Wang and Goldman-Rakic, 2004; Gee et al., 2012; Vitrac et al., 2014). Therefore, we propose that the therapeutic profile of neurostimulation methods aimed at modulating motor cortex activity, which have already demonstrated efficacy in improving motor symptoms of PD (Khedr et al., 2003; Fregni et al., 2006; Yang et al., 2018), could potentially be significantly augmented by cell type-specific modulation of cortical dopaminoceptive circuits, through pharmacological activation of DA receptors.
Footnotes
The authors declare no competing financial interests.
We thank Monika Baginska for excellent technical assistance in genotyping and maintaining the colony of transgenic animals.
This study was supported by The National Science Centre (Poland) grant SONATINA 2018/28/C/NZ4/00102 (to P.E.C.).
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.