Target-specific M1 inputs to infragranular S1 pyramidal neurons

Transcription

1 J Neurophysiol 116: , 216. First published June 22, 216; doi:1.1152/jn Target-specific M1 inputs to infragranular S1 pyramidal neurons Amanda K. Kinnischtzke, Erika E. Fanselow, and Daniel J. Simons Center for Neuroscience, University of Pittsburgh, Pittsburgh, Pennsylvania; and Department of Neurobiology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania Submitted 16 November 215; accepted in final form 16 June 216 Kinnischtzke AK, Fanselow EE, Simons DJ. Target-specific M1 inputs to infragranular S1 pyramidal neurons. J Neurophysiol 116: , 216. First published June 22, 216; doi:1.1152/jn The functional role of input from the primary motor cortex (M1) to primary somatosensory cortex (S1) is unclear; one key to understanding this pathway may lie in elucidating the cell-type specific microcircuits that connect S1 and M1. Recently, we discovered that a subset of pyramidal neurons in the infragranular layers of S1 receive especially strong input from M1 (Kinnischtzke AK, Simons DJ, Fanselow EE. Cereb Cortex 24: , 214), suggesting that M1 may affect specific classes of pyramidal neurons differently. Here, using combined optogenetic and retrograde labeling approaches in the mouse, we examined the strengths of M1 inputs to five classes of infragranular S1 neurons categorized by their projections to particular cortical and subcortical targets. We found that the magnitude of M1 synaptic input to S1 pyramidal neurons varies greatly depending on the projection target of the postsynaptic neuron. Of the populations examined, M1-projecting corticocortical neurons in L6 received the strongest M1 inputs, whereas ventral posterior medial nucleus-projecting corticothalamic neurons, also located in L6, received the weakest. Each population also possessed distinct intrinsic properties. The results suggest that M1 differentially engages specific classes of S1 projection neurons, thereby regulating the motor-related influence S1 exerts over subcortical structures. barrel cortex; sensorimotor integration; cortical circuitry NEW & NOTEWORTHY We examined the specificity of primary motor cortex (M1) inputs to specific classes of primary somatosensorty cortex (S1) pyramidal neurons and found that M1 axons synapse on infragranular S1 neurons with greater selectivity than previously appreciated. These results suggest M1 can differentially regulate particular circuits in S1 under conditions of high M1 activity (e.g., whisking). Additionally, each type of pyramidal neuron has a unique set of intrinsic properties that allow them function differently within the circuit. DURING ACTIVE EXPLORATION, sensory and motor systems act in concert in order for sensory feedback to inform ongoing and future movements. Although different regions within the central nervous system process sensory and motor information, the structures are highly interconnected, forming large-scale sensorimotor integration systems (Kleinfeld et al. 26; Petersen 27). At the level of the telencephalon, primary somatosensory (S1) and primary motor (M1) cortices are reciprocally connected (Cauller et al. 1998; Miyashita et al. 1994; Porter and White 1983; White and DeAmicis 1977). Both anatomic Address for reprint requests and other correspondence: A. Kinnischtzke, Dept. of Neuroscience, Columbia Univ., 53 Hammer Health Science Bldg., 71 W. 168th St., New York, NY 132 ( ak3573@cumc.columbia.edu). and physiological studies have demonstrated that S1 sends a robust projection to superficial layers in M1, in particular contacting neurons that project back to S1 (Aronoff et al. 21; Mao et al. 211). Recent studies have determined that M1 densely innervates and provides the strongest input to the infragranular layers in S1 (Kinnischtzke et al. 214; Veinante and Deschenes 23; Zagha et al. 213; Zhang and Deschenes 1998), yet the specificity of its connections with excitatory cells in the deep layers remains poorly understood. The somatic and dendritic morphology of pyramidal neurons is relatively uniform, at least compared with inhibitory interneurons, and consequently they are often considered as a single population. However, accumulating evidence demonstrates that pyramidal neurons differ functionally in terms of their intrinsic properties and their local and extrinsic projections. Recent studies have found that excitatory connections among local S1 pyramidal neurons are not random but instead form sub-circuits based on similarity of physiological properties (Otsuka and Kawaguchi, 28) receptive field structure (Kampa et al. 211; Cossell et al. 215), and/or the projection target of the neurons (Brown and Hestrin, 29; Otsuka and Kawaguchi, 211). Our previous study (Kinnischtzke et al. 214) demonstrated that M1 engages pyramidal cells distributed throughout the cortical layers, consistent with the view that M1 sends a general modulatory signal to S1. On the other hand, we also observed significant variation in the strengths of M1 input, with a subset of neurons in the superficial aspect of L6 receiving especially strong inputs. If M1 synapses specifically form stronger connections with pyramidal neurons projecting to a common target, it suggests M1 may regulate S1 communication with long-range cortical and/or subcortical structures in a more precise manner than previously suggested. Layer (L)5 and L6 of S1 contain a diversity of pyramidal neuron populations that project to a number of cortical and subcortical sites. Here, we investigated synaptic input from M1 onto several types of infragranular S1 pyramidal neurons: pyramidal cells projecting back to M1 (M1 projecting), corticothalamic neurons projecting to either the posteromedial thalamus (PoM projecting) or ventral posterior medial nucleus (VPM projecting), and corticofugal neurons projecting to the spinal trigeminal nucleus (Sp5 projecting). Our results demonstrate that both the intrinsic electrophysiological properties as well as strength and effectiveness of the M1 inputs the neuron receives are strongly related to its projection target. METHODS All experiments were carried out in compliance with University of Pittsburgh School of Medicine animal use policies and were approved by the University of Pittsburgh Institutional Animal Care and Use Committee /16 Copyright 216 the American Physiological Society 1261

2 1262 Channelrhodopsin and retrograde labeling procedures. Viral constructs carrying channelrhodopsin (ChR)2 and mcherry genes were injected into the mouse primary motor cortex as previously described (Kinnischtzke et al. 214). Briefly, mice (1 16 postnatal days of age) were anesthetized with isoflurane (1 2% in oxygen), and a craniotomy was made over the vibrissal area within M1 (1. mm anterior to and.8 mm lateral from the bregma). The virus was pressure injected at depths of.9 and.4 mm (corresponding to the deep and superficial layers in M1) within a single penetration. A volume of.1.2 l was injected at each depth. We used the adeno-associated virus (AAV) construct AAV2/5.CamKII.hChR2(H134R)-mCherry.WPRE.SV4 (University of Pennsylvania Vector Core, permission from Dr. Karl Deisseroth) for all experiments, as this construct results in strong anterograde expression but virtually no retrograde transport (Kinnischtzke et al. 214). During the same surgery, mice were also injected with red and green fluorescent retrograde tracers (RetroBeads, Lumafluor). All mice had green retrobeads injected in M1; the undiluted beads were typically coinjected simultaneously with the ChR2 virus (at 1:4 with the virus to achieve 25 5 nl of retrobead injection along with 1 2 nl of virus injection). Mice were then injected with red fluorescent retrobeads in one of the following locations: VPM, PoM, or Sp5. For injections in the thalamus, a small craniotomy was performed dorsal to VPM (1.4 mm posterior and 1.7 mm lateral to the bregma) or PoM (1.6 mm posterior and 1.4 mm lateral to the bregma). Beads were pressure injected using a picospritzer (5 1 nl; depth of 3.2 mm for VPM and 3.1 mm for PoM). The pipette was withdrawn slowly 5 1 min after injection to limit inadvertent leakage in the cortex. For Sp5 injections, the brain stem was exposed directly posterior to the cerebellum, and the pipette was advanced at a 4 angle. Retrobeads were pressure injected using a picospritzer at depths of.8 and.6 mm (.1 l at each depth). Experiments were conducted using both male and female mice from the GIN transgenic mouse line (Oliva et al. 2) to be consistent with our previous study. Slice preparation. Electrophysiological experiments were performed 4 8 wk after the injection surgery to allow for transport and full expression of the virus. Mice were between postnatal days 44 and 74 at the time of experiments. This age is not fully adult yet past the stage of major cortical reorganization; therefore, we expect these data to be similar to an adult state. The animal was anesthetized with isoflurane, and the brain was removed and placed in ice-cold artificial cerebrospinal fluid (ACSF) containing (in mm) 126 NaCl, 3 KCl, 1.25 NaH 2 PO 4, 2 MgSO 4, 26 NaHCO 3, 1 dextrose, and 2 CaCl 2 saturated with 95% O 2-5% CO 2. The tissue was cut into 4- m slices in the coronal plane using a vibratome. Slices were incubated at 32 C for 3 45 min and then maintained at room temperature until used for recording. Slices containing S1 barrel cortex were identified by the presence of layer IV barrels and a barrel-related pattern of ChR2- mcherry fluorescence, as labeling was dense in deep layers, L5 and L6, and then ascended between the barrels in and L4 (see Fig. 4A). Confirmation of the injection targets was confirmed by examination of fluorescence at the injection site as well as by presence of the expected pattern of laminar cortical labeling. If animals injected with retrobeads in VPM exhibited significant L5B cortical labeling, the injection was considered PoM contamination and the animal was not included in the study. We recorded primarily from the larger, more medially situated barrels (rows D and E), as this is where the ChR2-mCherry labeling was typically strongest. We recorded from pyramidal neurons in 1 3 adjacent barrel-related columns/slice. Recording procedures. All recorded cells were positive for either red or green retrobeads, which allowed for unambiguous classification of the projection target. Whole cell recordings were performed using glass micropipettes (4 1 M ) filled with internal solution containing (in mm) 135 K-gluconate, 4 KCl, 2 NaCl, 1 HEPES,.2 EGTA, 4 ATP-Mg,.3 GTP-Tris, and 14 phosphocreatine-tris (ph 7.25, mosm). Biocytin (.5%) was added to the internal solution in most experiments. Membrane potentials reported here were not corrected for liquid junction potential. Recordings were conducted at 32 C. When patching, cell-attached seal resistances were 1G, and series resistance after achieving whole cell configuration was 5 2 M. After establishing whole cell configuration, a series of current steps was presented in current clamp for use in analysis of the cell s electrophysiological characteristics. Current steps were 6 ms in duration and presented in 2-pA steps ranging from 1 to 3 pa. Steps were presented 5 s apart. In voltage-clamp experiments, series resistances ranged from 1 to 4 M and were compensated for up to 7%. For excitatory postsynaptic current (EPSC) measurements, all cells were held at a potential of 8 mv. Data were collected using a Multiclamp 7B amplifier and pclamp1 software (Molecular Devices) at a sampling rate of 2 khz. To reduce between-experiment variability related to differences in ChR expression levels, neurons were recorded in a paired manner. For each slice, at least one M1-projecting neuron was recorded in addition to cells having identified projection targets (e.g., Sp5, PoM, or VPM). The amplitude of the M1 input to cells recorded across different animals could then be normalized by the M1 input to a standard population, M1-projecting cells. Laminar definitions. Measures of laminar position are reported as normalized distance from the top of (value of.) and the bottom of L6 (value of 1.). The L4/L5 boundary is distinct, and the normalization procedure accounted for slight variability between experiments in cortical depth due, for example, to section plane. Laminar boundaries were based on Nissl-stained coronal sections in S1. The actual, non-normalized, laminar boundary measurements were similar to previously reported values (Hooks et al. 211). Borders were identified on the basis of changes in cell densities indicated by Nissl stain. The L4/ boundary was distinguishable by an abrupt loss of barrel structure and a decrease in the number and density of cells. The bottom of was typically 1 m below the L4/ border. The boundary between L5B and L6, visible as a notable increase in cell density near the top of L6, was typically located approximately halfway between the top of and the subcortical white matter. Optical stimulation procedures. All experimental procedures were carried out as previously described (Kinnischtzke et al. 214). Fullfield blue light (47 nm, OptoLED, Cairn Research) was delivered through a 4 objective; light intensity at the surface of the slice was 2 mw/mm 2. Each cell recorded was centered in the field of view to maximally activate the synapses onto that particular neuron. The spatial diameter of light through the 4 objective was 25 m, thereby covering the width of approximately one barrel-related column. A single pulse (duration: 1. ms) was delivered per trial, and individual trials were separated by 8 s. Pulse duration (1 ms) was long enough to reliably elicit synaptic responses from postsynaptic neurons yet short enough to limit recruitment of polysynaptic activity. Data analysis. Data analyses were performed using custom-written MATLAB code (The MathWorks, Natick, MA; A. K. Kinnischtzke). Intrinsic properties were derived from a series of direct current steps that were presented to the cell (see above). Input resistance was calculated as the slope of the voltage-current relationship for subthreshold voltage deflections that were within 3 mv of the resting membrane potential. Using hyperpolarizing current steps, sag current magnitude (in mv) was measured as the difference in the minimum voltage during the initial 1 ms and the mean voltage during the last 5 ms in the step. Because sag currents exhibit a voltage dependence, only current steps where the minimum voltage value was between 9 and 75 mv were used. Action potential (AP) features were calculated from the first AP in the first current step that elicited at least two APs at 1 Hz or greater. The AP threshold was measured as the maximum of the second derivative in the voltage 1 ms preceding the AP peak. AP half-width was calculated as the width halfway between the AP threshold and the AP peak. To measure afterhyperpolarization (AHP) depth, the mini-

3 1263 mum voltage (V min ) within 5 ms after the AP threshold was subtracted from the AP threshold (see Fig. 3D). To determine the presence or absence of a depolarizing afterpotential (DAP), the local maximum (V max ) was found between V min and the voltage during the following 1 ms. The voltage 17 ms after V min was also calculated (V end ). If V max V end, the DAP magnitude was calculated as the difference between V max and V min (see Fig. 3D). If V max V end, the DAP magnitude was. Rheobase current was the first current step amplitude that elicited at least one AP. All spike train properties were calculated from the first current step that resulted in an average of 1 Hz firing or greater. This was typically 3 pa (range: 8 pa) above the rheobase current. A burst was defined as two or more APs with an interspike interval (ISI) of 12 ms or less. An initial adaptation ratio was defined as the ISI between the first and second APs divided by the ISI between the second and third APs. A steady-state adaptation ratio was calculated as the ISI between the third and fourth spikes divided by the ISI of the last two spikes in the train. For the analysis of M1 input, individual trials were averaged together because we observed little trial-to-trial variability in the responses to optical stimulation. We did not use methods for limiting polysynaptic activity (Petreanu et al. 29). Although these conditions are ideal for unambiguous detection of the presence of monosynaptic connections, the addition of TTX and 4-aminopyridine to block Na and K channels, respectively, alters the normal synaptic strength and properties of synaptic connections (Cruikshank et al. 21) making comparison of synaptic strengths difficult. In addition, blockade of APs with the addition of TTX would not allow for an examination of the intrinsic properties of the cells. Instead, we limited our analysis of postsynaptic responses to an initial time window after the light pulse (1 ms) to minimize the contribution of polysynaptic activity. Activity in this early time window (first 1 ms) also showed very little trial-to-trial variability, consistent with it being monosynaptic in nature; however, it is possible that some of this activity could include polysynaptic events that were reliably evoked by the optical stimulation. For the analysis of EPSCs, a cell was considered to have received input if the peak of the EPSC exceeded 15 times the SD of the baseline preceding the response. This detected EPSC responses of 5 pa and larger. To capture the start of the EPSC, we measured the first of 2 consecutive data points that exceeded 1 SD below the baseline current (i.e., exceeded the noise level). This value was taken as the response onset. EPSC amplitude was calculated as the EPSC peak minus the value at the EPSC onset. In current clamp, spike probability was measured as the number of trials that elicited a spike divided by the total number of trials. Statistical comparisons of multiple groups were performed using ANOVA followed by multiple pairwise comparison tests that were corrected using the Tukey-Kramer method. Statistical comparison of proportions or probabilities between groups was tested using a Fisher s exact test. Comparison of EPSC amplitudes for within-slice recordings was conducted using a paired t-test. RESULTS Anatomic distribution of L5 and L6 pyramidal neurons. To test for M1 input to pyramidal neurons projecting to a specific cortical or subcortical target, we injected fluorescent retrograde tracers in combination with an AAV containing the gene for ChR2 (see METHODS). In all animals, a green fluorescent retrograde tracer was coinjected with the ChR2 virus into M1 and a red fluorescent retrograde tracer was injected into either the VPM, PoM, or Sp5 in the brain stem (Fig. 1). As previously described, M1 axons are found throughout L1 L6 in S1 but are especially dense in L1, L5, and L6 (Fig. 1A) (Kinnischtzke et al. 214). Although M1 axons diffusely and uniformly terminated in the deep layers (Fig. 1A), pyramidal neurons projecting to specific downstream targets instead tended to cluster in specific laminar locations within L5 or L6 (Fig. 1B). For each retrograde injection target (i.e., M1, PoM, Sp5, or VPM), the position of retrogradely labeled neurons was consistent with known laminar patterns in rodent S1 (Alloway et al. 24; Chmielowska et al. 1989; Deschenes et al. 1994; Hattox and Nelson 27; Killackey et al. 1989; Mao et al. 211; Wise and Jones 1977). M1-projecting neurons were found throughout most layers in S1 but tended to be concentrated in bands within L2/L3,, superficial L6, and deep L6. Neurons projecting to Sp5 were located exclusively in L5B. Corticothalamic neurons projecting to PoM were found in L5B, in a similar position as Sp5-projecting neurons, as well as in deep L6. Finally, corticothalamic neurons projecting to VPM were only found in L6 and tended to be concentrated in the upper half of L6, as previously described (Chmielowska et al. 1989; Zhang and Deschenes 1997). The retrogradely labeled neurons recorded in this study were grouped into the following five populations for further study: M1-projecting cells in ( M1-projecting), M1-projecting neurons at the very top of L6 (L6 M1-projecting), Sp5- projecting neurons, L5B PoM-projecting neurons (PoM-projecting), and VPM-projecting neurons. Because only two L6 PoM-projecting neurons were recorded, they were not included as an additional population (see DISCUSSION). Depth measurements, normalized to the L4/ border, revealed that each of the studied populations of pyramidal neurons were largely segregated within a layer, e.g., a sublayer, with only a small degree of overlap among populations residing within that layer (Fig. 1C). The exception was that PoM-projecting neurons in L5B and Sp5-projecting neurons were completely overlapping in their depth location (Fig. 1C). Coexpression of red and green fluorescent retrobeads was observed in only two neurons, both from animals injected in M1 and Sp5 (n 2; data not shown). These cells were located in the deep part of L5B and were included in both Sp5- and L6 M1-projecting data sets. No other cases of colabeling of green (M1-projecting) and red (PoM-, Sp5-, or VPM-projecting) fluorescent retrobeads were observed in this study. Examination of the labeled pyramidal neurons revealed morphological differences among them consistent with those of similarly identified populations of S1-projecting neurons (Chagnac-Amitai et al. 199; Deschenes et al. 1994; Hattox and Nelson 27; Le Be et al. 27; Zhang and Deschenes 1997); these anatomic features were therefore were not investigated further. Intrinsic properties. We compared the electrophysiological characteristics of M1-projecting neurons (n 23), L6 M1-projecting neurons (n 19), Sp5-projecing neurons (n 24), L5B PoM-projecting neurons (n 15), and VPM-projecting neurons (n 19). As seen in juvenile mice (Hattox and Nelson 27), many subthreshold and AP characteristics were strongly dependent on the projection target of the pyramidal neuron (Table 1 and Figs. 2 and 3). All statistical relationships are shown in Table 1. Comparisons of passive properties demonstrated that VPMprojecting neurons had the highest input resistance ( M ) and PoM- and Sp5-projecting neurons had the lowest input resistance ( and M, respectively;

4 1264 A B1 2 Sp5-proj. cells M1-proj. cells L2/3 L2/3 L2/3 L4 L4 L4 L5B L5B L5B L6 L6 L6 C 3 Norm. depth from top of 4 VPM-proj. cells PoM-proj. cells.2 L5B.4 L2/3.8 L5B L5B Number of cells M1-proj. Sp5-proj. PoM-proj. L6 M1-proj. L4 L4.6 L6 L2/3 L6 L6 VPM-proj. Fig. 1. Laminar positions of primary somatosensory cortex (S1) pyramidal neurons in layer (L)5 and L6 are organized by projection target. A: confocal image of channelrhodopsin (ChR)2-primary motor cortex (M1) axons within S1 labeled with mcherry showing the typical pattern of innervation within S1. B: retrobead injections into M1 (1) resulted in retrogradely labeled pyramidal neurons throughout L2 L6. Labeled cells were especially numerous in, the top of L6, and the bottom of L6. Red fluorescent retrobeads injected into the spinal trigeminal nucleus (Sp5; 2) exclusively labeled S1 neurons in L5B. Injection of retrobeads into the posteromedial thalamus (PoM; 3) labeled neurons in both L5B and L6, whereas retrobeads in the ventral posterior medial nucleus (VPM; 4) only labeled neurons in L6 of S1. In 1 4, the insets illustrate the injection target (top; ) and an example injection for each target location (bottom). C: laminar depth measurements for each recorded cell were normalized to the top of. Histograms illustrate the distribution of each population within L5 and L6. Open circles and error bars (right) indicate the means and SDs, respectively, for each population. Fig. 3A). These differences are inversely related to the size of the neurons, as Sp5- and PoM-projecting neurons are among the largest neurons in S1 and VPM-projecting neurons are the smallest. M1- and L6 M1-projecting neurons had significantly more hyperpolarized resting membrane potentials ( and mv, respectively) than the other three pyramidal neuron types. We also examined the presence of sag potentials, which in pyramidal neurons is mediated by hyperpolarization-activated current (Ih) (Sheets et al. 211). PoM- and Sp5-projecting neurons exhibited the most sag potential ( and mv, respectively), and sag potentials in PoM-projecting cells were significantly larger than those of M1-, VPM-, and L6 M1-projecting neurons (Fig. 3, B and C). L6 M1projecting neurons, in particular, had notably small sag potentials ( mv). Each population of pyramidal neurons exhibited distinct AP features as well. Sp5- and PoM-projecting neurons had the narrowest APs (AP half-width of and mv, respectively) as well as the most hyperpolarized threshold for AP initiation ( and mv, respectively). By comparison, APs of VPM-projecting neurons were the broadest (AP half-width of ms) and had the most depolarized AP threshold ( mv). and L6 M1-projecting neurons had broader APs than Sp5- or PoM-projecting neurons but narrower than VPM-projecting neurons ( M1-projecting neurons: ms and L6 M1-projecting neurons: ms). AP thresholds of and L6 neurons were more hyperpolarized than those of VPM-projecting neurons but were more depolarized than those of Sp5- or PoM-projecting neurons ( M1-projecting neurons: mv and L6 M1-projecting neurons: mv). As described below, AP thresholds were important determinants of whether the different classes of neurons tended to fire spikes in response to M1 inputs. Striking differences between pyramidal neuron subtypes were observed in the voltage trajectory after the AP. VPMprojecting neurons were distinct in exhibiting a deep, fast AHP that was not seen in the other cell types (Figs. 2B and 3, D and E). Sp5- and PoM-projecting neurons, on the other hand,

5 1265 Table 1. Table of intrinsic properties Property Projection Target Significance Level M1 (1) Sp5 (2) PoM (3) L6M1(4) VPM (5) P.5 P.1 P.5 Number Passive properties Input resistance, M :4 1:2 5:1,2,3,4 Resting membrane potential, mv :5 4:2,3,5 Sag current, mv :5 4:1,2,3 3:1,4 4:5 AP features AP half-width, ms :2,3 4:2,3 5:2,3 AP threshold, mv :2,3,5 5:2,3,4 AP peak, mv NS NS NS Afterhyperpolarization depth, mv :3 5:1,2,4 Percentage with depolarizing afterpotential 59 (13/22) 92 (22/24) 87 (13/15) 94 (15/16) 5 (1/19) 1:2 1:4 5:1,2,3,4 Spiking properties Rheobase, pa NS NS NS Initial adaptation ratio :4 5:1,2,4 Steady-state adaptation ratio :1,4,5 Percentage with initial burst 4 (1/23) 33 (8/24) 27 (4/15) 58 (11/19) (/19) 1:2 1:4 4:2,3 5:2,4 3:5 Percentage with steady-state bursts (/23) 17 (4/24) 13 (2/15) 5 (1/19) (/19) NS NS NS Values are means SE or percentages with numbers in parentheses. Primary somatosensory cortex (S1) pyramidal neurons exhibit diverse intrinsic properties based on the projection target. Statistical relationships are shown on the right. See METHODS for a description of each property. The right side indicates statistical relationships between each pyramidal neuron type. Group numbers (1 5) are indicated at the top and as follows: layer (L5)A primary motor cortex (M1)-projecting neurons (1), spinal trigeminal nucleus (Sp5)-projecting neurons (2), posteromedial thalamus (PoM)-projecting neurons (3), L6 M1-projecting neurons (4), and ventral posterior medial nucleus (VPM)-projecting neurons (5). AP, action potential; NS, not significant. displayed a small DAP after a spike (Figs. 2, C and D, and 3F). DAPs were observed in many L6 M1-projecting and some M1-projecting neurons as well (Fig. 3F). Interestingly, in Sp5-, PoM-, and some L6 M1-projecting neurons, a DAP was present after every spike (Fig. 2, C, D, and E, right), whereas for M1-projecting and the remaining L6 M1-projecting neurons, a DAP, if present at all, was only found after the first AP in a train (Figs. 2, A and E, left, and 3F). A B C M1-proj. VPM-proj. Sp5-proj mv 2 ms D PoM-proj. E L6 M1-proj Fig. 2. Examples of firing properties for S1 pyramidal neuron populations. A: M1-projecting neurons exhibited mostly regular spiking behavior, and the first interspike interval (ISI) tended to be shorter than subsequent ISIs. B: VPM-projecting neurons were also regular spiking with large afterhyperpolarizations (AHPs) after each action potential (AP). Sp5-projecting neurons (C) and PoM-projecting neurons (D) fired either single spikes that were followed by a depolarizing afterpotential (DAP; C and D, left) or a mix of bursting and regular spiking (C, right) or showed intrinsic bursting dynamics (D, right). E: L6 M1-projecting neurons often fired an initial spike doublet, and in some cells, subsequent spikes were followed by a DAP (left), whereas others were not (right). In A E, the first spike in the train was expanded and shown to the right. Scale bars for expanded spike in A E are 1 ms and 1 mv.

6 1266 A B C 7 Input res. (M ) ms 2mV Sag magnitude (mv) M1 Sp5 PoM L6-M1 VPM -M1 Sp5 PoM L6-M1 VPM D E 18 F (1) VPM-proj. PoM-proj. 1 mv 5ms AHP DAP 1 mv (2) 5ms 5mV 5ms AHP depth (mv) % cells with DAP st AP 2nd AP -M1 Sp5 PoM L6-M1 VPM -M1 Sp5 PoM L6-M1 VPM % of cells RS cells B/RS cells IB cells G H I -M1 Sp5 PoM L6-M1 VPM AHP ratio M1 Sp5 PoM L6-M1 VPM Steady-state adapt. ratio Initial adapt. ratio Fig. 3. Comparison of intrinsic properties between S1 pyramidal neurons populations. A: population average of input resistance for each type of pyramidal neuron. B: illustration of sag currents. Traces are averaged voltage responses for all neurons in each population. The dashed box in the inset illustrates the portion of the current step shown in the averaged traces. C: average sag currents were largest in PoM-projecting neurons and smallest in L6 M1-projecting neurons. D: first spike waveforms (left), averaged across neurons, highlight the differences between populations. Inset 1 shows measurements for AHP depth (left) and DAP magnitude (right). Inset 2 shows the expansion of voltage trajectory in the dashed box and highlights differences in the voltage trajectories for each cell type after an AP. The color scheme is as in B. E: population average of AHP depth. F: percentage of neurons in each population demonstrating DAPs after the first spike in the current step (solid bars) and second spike in the current step (open bars). G: number of cells in each population that showed regular spiking (RS; open bars), initial bursts followed by regular spiking (B/RS; gray bars), or intrinsic bursting (IB; solid bars). H: ratio of the first to second AHP depth in a spike train. I: VPM-projecting neurons (blue dots) displaying distinct adaptation dynamics, as the initial adaptation ratio was higher than the steady-state adaptation ratio. The color scheme for all images with color is as follows: M1-projecting (gray), Sp5-projecting (magenta), PoM-projecting (green), L6 M1-projecting (red), and VPM-projecting (blue). In A I, statistical relationships are not shown and instead are shown in Table 1. Single dots indicate individual data points, and all bars and error bars indicate means and SEs. DAPs may underlie the generation of bursting behaviors observed in some pyramidal neurons (Llano and Sherman 29). We therefore compared the proportions of pyramidal neurons that were either completely regular spiking (RS), displayed an initial burst followed by single spikes (B/RS), or exhibited bursting throughout the current step [intrinsic bursting (IB); Fig. 3G]. Consistent with DAPs relating to bursting dynamics, pyramidal neuron types that were most likely to exhibit DAPs, particularly after the second and later APs (e.g., Sp5-, PoM-, and L6M1-projecting neurons; Fig. 3F, solid bars), were also more likely to exhibit initial and steady-state bursting behavior (Fig. 3I). VPM-projecting cells were always RS (19/19 cells; Fig. 2B). -M1 projecting cells were almost always RS cells (19/23 cells; Fig. 2A), with a few cells exhibiting an initial burst followed by regular spiking behavior (B/RS; Fig. 3I). Neurons projecting to Sp5 could be RS, B/RS,

7 or IB cells (RS: 12/24 cells, B/RS: 7/24 cells, and IB: 5/24 cells; Fig. 2C); the same was true for PoM-projecting neurons (RS: 8/15 cells; B/RS: 5/15 cells; and IB: 2/15 cells; Fig. 2D). Interestingly, most L6 M1-projecting neurons were B/RS (14/19 cells) as they often exhibited an initial spike doublet (Fig. 2E). Most of the remaining L6 M1-projecting neurons were RS (4/19 cells), although one was IB (1/19 cells). Finally, we examined changes in spike properties across the spike train. Again, VPM-projecting neurons were distinctive as AHP depth decreased from the first spike to subsequent spikes, resulting in an AHP ratio that was larger than 1 (Fig. 3H). In the other cell types, the first AHP was shallower than subsequent AHPs and AHP ratios were 1 (Fig. 3H). VPM-projecting neurons also displayed different adaptation dynamics, with the ISI gradually increasing across the spike train (Fig. 2B). This resulted in initial adaptation ratios larger than the steadystate adaptation ratios (.9.5 vs..79.3; Fig. 3I). For each of the other cell types, the initial adaptation ratio was typically smaller than the steady-state adaptation ratio. This A DIC was particularly evident for the neurons that exhibited an initial burst (Fig. 3I, open circles). IB cells were not included in analysis of adaptation dynamics as the irregular spike patterns within and between bursts yielded highly variable adaptation values. Motor cortex inputs. To account for possible differences in ChR2 expression between slices, all animals were injected with ChR2 and green retrobeads in M1 as well as red retrobeads in a second target (either VPM, PoM, or Sp5). We were therefore able to use a paired design in which M1-projecting neurons were recorded in the same slice as cells from one of the other pyramidal neuron populations: L6 M1-, Sp5-, PoM-, or VPM-projecting neurons (Figs. 4 and 5). Paired comparisons were used to determine the strength of M1 inputs relative to a standard population ( M1-projecting neurons). Comparison between M1-projecting neurons and unlabeled neurons confirmed that this population was representative of the amount of M1 input received by all neurons (data not shown). B fluorescence 1267 M1-proj. Sp5-proj. 2 mv 4x 1 ms x -67 mv -64 mv 1 pa 4 mv 1 ms C M1-proj. D PoM-proj. E n.s. n.s EPSC amp. (pa) 1 ms EPSC amp. (pa) 2 mv mv 5 pa 1 ms M1 Sp5 M1 PoM Fig. 4. M1 input to M1-, Sp5-, and PoM-projecting neurons is equivalent. A: example of a S1 slice at 4 magnification under DIC illumination (top left) and epifluorescence (top right) showing bright, retrogradely labeled neurons in L5B that project to Sp5. Diffuse fluorescent labeling was also observed, indicating the presence of ChR2-expressing M1 axon terminals. An example of a recorded pyramidal neuron under DIC optics (bottom left) that expressed fluorescent retrobeads (bottom right) is shown. B: example M1-projecting (left) and Sp5-projecting (right) neurons recorded in the same slice that showed characteristic current step responses (top) and received similar excitatory input from M1 when recorded in current clamp (middle) or voltage clamp (bottom). For current-clamp recordings, multiple traces were overlaid, and the average excitatory postsynaptic potential is shown in black. Excitatory postsynaptic current (EPSC) traces show only the trial-averaged response. C: example M1-projecting and L5B PoM-projecting neurons recorded in the same slice. Conventions are as in B. Cells were held at 8 mv for voltage-clamp recordings in B and C. Population-averaged EPSC responses to M1 stimulation showed that Sp5-projecting neurons (D) and PoM-projecting neurons (E) receive a statistically equivalent magnitude of M1 input as M1-projecting neurons ( M1-projecting vs. Sp5-projecting neurons: P.31; M1-projecting vs. PoM-projecting neurons: P.27). In D and E, the closed circles, open squares, and denote RS cells, B/RS cells, and IB cells, respectively. The open circles and error bars to the left and right of the paired data are the population means and SEs.

8 1268 The recording from a M1-projecting neuron and a Sp5-projecting neuron within the same slice (Fig. 4A) demonstrated that optical stimulation of ChR2-M1 axons evoked a short-latency EPSC that was similar in amplitude for both cells (Fig. 4B). Comparisons between the populations of recorded pairs (n 12 pairs) indicated no significant difference in M1 input between these populations ( M1-projecting neurons: A M1-proj. L6 M1-proj. 2 mv 2 ms mv -62 mv 2 mv 5 pa 1 ms B M1-proj. VPM-proj. 2 mv 2 ms mv 5 pa 1 ms EPSC amp. (pa) D 15 EPSC amp. (pa) C -M1 L6-M pa and Sp5-projecting neurons: pa, P.31; Fig. 4D). Similar results were found for M1- and PoM-projecting neurons (n 9 pairs), indicating that M1-projecting and PoM-projecting neurons also receive inputs of equivalent strength from M1 ( M1projecting neurons: pa and PoM-projecting neurons: pa, P.27; Fig. 4, C and E). Sp5and PoM-projecting cells that exhibited either initial bursts (Fig. 4, D and E, open circles) or were intrinsically bursting (Fig. 4, D and E, ) received a comparable magnitude of M1 input as the cells that were regular spiking (Fig. 4, D and E, closed circles). Thus, Sp5- and PoM-projecting neurons, which, as noted above, have similar sublaminar locations and intrinsic properties, were similar with respect to the strengths of the M1 inputs they receive. Unlike M1-, Sp5-, and PoM-projecting neurons, many L6 M1-projecting neurons had large-amplitude M1-evoked responses (Fig. 5A). Indeed, M1 input was often strong enough to evoke APs in L6 M1-projecting neurons; in contrast, APs were never observed in M1-projecting neurons (Fig. 5A; see also below). Across the population (n 9 pairs), L6 M1-projecting neurons on average received significantly stronger input from M1 than M1-projecting neurons ( vs pa, P.1; Fig. 5C). Many of the L6 M1-projecting neurons that received strong M1 input exhibited an initial spike doublet (Fig. 5C, open squares); large-amplitude M1 inputs were observed in regular spiking cells as well (Fig. 5C, closed circles). M1 input to VPM-projecting neurons was notably small, and these cells received on average the smallest M1 inputs (n 16 pairs; Fig. 5, B and D). In the sample of paired cells, the amplitude of M1 input to VPM-projecting neurons was significantly smaller than of M1-projecting cells recorded in the same slice ( vs pa, P.8). The two L6 PoM-projecting cells we recorded had similar intrinsic properties and received comparable amounts of M1 input as VPM-projecting neurons (data not shown). Comparisons normalized to within-slice M1 inputs confirmed that L6 M1-projecting cells received significantly stronger and VPM-projecting cells received significantly weaker synaptic input from M1 than M1-projecting neurons (P.5; Fig. 6A). Because each sample of paired M1projecting neurons was similar (P.43 by ANOVA), we pooled the data to compare all five groups directly (Fig. 6B). The pooled data showed similar results as the data from individual pairs, confirming that synaptic inputs from M1 differentially targeted populations of pyramidal neurons based on their axonal projections (P.1 by ANOVA; Fig. 6B), with L6 M1-projecting neurons receiving significantly -M1 VPM Fig. 5. M1 provides strong input to L6 M1-projecting neurons and weak input to VPM-projecting neurons. A: example M1-projecting (left) and L6 M1-projecting (right) neuron responses to current steps (top). In this pair, the L6 M1-projecting neuron received larger M1 input, as shown in current clamp (middle) and voltage clamp (bottom). Note that the L6 M1-projecting neuron fired spikes on several trials (one trial is shown in the dashed box inset). Cells were recorded in the same slice. B: same as in A but for a M1-projecting neuron (left) and a VPM-projecting neuron (right) recorded in the same slice. C: on average, L6 M1-projecting neurons received significantly larger input than paired M1-projecting cells (P.5). D: VPM-projecting neurons received significantly less M1 input than paired M1-projecting neurons (P.5). Conventions for C and D are the same as in Fig. 4, D and E.

9 1269 A Normalized M1 input M1 Sp5 PoM M1 L6 VPM B EPSC amp (pa) M1 Sp5 PoM L6 M1 VPM C Spike probability Spike probability M1-proj PoM-proj EPSC amp (pa) Sp5-proj VPM-proj EPSC amp (pa) D 1 E 4 F G Spike probability M1 Sp5 PoM L6 M1 VPM AP thresh - V rest 2 Sp5 PoM L6 VPM M1 M1-59 Sp5-proj. 1 mv 1 ms 5 pa -72 L6 M1-proj. Spike probability R 2 =.52 p = AP thresh - V rest Fig. 6. Optical stimulation of ChR2-M1 terminals evoked suprathreshold responses in multiple types of S1 pyramidal neurons. A: paired recordings were normalized to the magnitude of M1 input to the M1-projecting cell. M1 input was significantly greater to L6 M1-projecting cells (P.5) and less to VPM-projecting cells (P.5) compared with the M1-projecting population. B: population average of EPSC amplitude evoked by M1 stimulation showed that L6 M1-projecting neurons received significantly larger M1 input than the other pyramidal neuron populations (P.5). C: pyramidal neurons in each population with the largest M1-evoked EPSC often showed nonzero spike probabilities when recorded in current clamp. Note that the top left image shows both M1-projecting (gray circles) and L6 M1-projecting (solid circles) neurons. D: average spike probability in L6 M1-projecting neurons was significantly greater than in M1- and VPM-projecting neurons (P.5) but not Sp5- and PoM-projecting neurons (P.5). E: the average voltage difference from the resting membrane potential to AP threshold was significantly smaller in Sp5- and PoM-projecting neurons than other pyramidal neuron classes (P.5). F: example Sp5-projecting neuron (left) that spiked on each trial (black) and a L6 M1-projecting neuron (right) that fired on some trials (black) but not others (gray). G: scatterplot of the distance between rest and AP threshold and the M1-evoked spike probability for all neurons that spiked. Symbol conventions are as follows: Sp5-projecting neurons (solid circles), PoM-projecting neurons (open triangles), L6 M1-projecting neurons (open circles), and VPM-projecting neurons (open square). The dashed line shows the linear least-squares fit. larger inputs than any other population and VPM-projecting neurons receiving the smallest (Fig. 6B). M1-evoked suprathreshold responses. Because these pyramidal neurons display a range of intrinsic properties, we also recorded in current-clamp to provide a more physiological measure of the cell s response to activation of M1 given the intrinsic properties of the neuron and amount of excitatory and inhibitory synaptic input produced by M1 activation. Analysis of the M1-evoked EPSP amplitude, rather than EPSC, illustrated a similar effect as described above, in that L6 M1- projecting neurons received significantly more M1 input than the other populations we recorded (data not shown). Many pyramidal neurons responded by spiking, usually those in populations that received the largest amplitude of excitatory input from M1 (Fig. 6C). For example, the L6 M1-projecting population had the highest proportion of neurons exhibiting suprathreshold responses (n 9/19). Interestingly, although M1-, Sp5-, and PoM-projecting neurons received similar mean amplitudes of input from M1 (Fig. 6, A and B), M1-projecting neurons never spiked (n /23), whereas several Sp5- and PoM-projecting neurons did (n 6/24 and n 3/15, respectively). VPM-projecting neurons received notably small excitatory inputs from M1 (Fig. 6, A and B), and, despite having the largest input resistance (Fig. 3A), they rarely responded by firing APs (n 1/2). Average spiking probability was highest in L6 M1-projecting neurons (.3.9; Fig. 6D); however, spiking probabilities in Sp5- (.23.8) and PoM-projecting (.19.1) neurons were statistically equivalent to L6 M1-projecting cells even though they received significantly smaller excitatory currents from M1. For example, Sp5- and PoM-projecting neurons reliably fired APs for EPSC amplitudes around 5 pa, as did a few L6 M1-projecting neurons. However, other L6 M1-projecting neurons spiked only occasionally despite receiving excitatory inputs of 1 pa or greater (Fig. 6C). Of the populations where two or more cells had a nonzero spike probability, the

10 127 amplitude of M1 input alone was not a reliable predictor of whether or not the cell would spike in response to photostimulation of M1 axons (L6 M1: R 2.29, Sp5: R 2.59, and PoM: R 2.74, P.5 for all three). Therefore, we further examined spiking probability as it related to both the magnitude of M1 input currents and the voltage difference between resting potential and AP threshold. The voltage difference from rest to AP threshold was smallest for Sp5- and PoM-projecting pyramidal neurons (P.5; Fig. 6E), reflecting the fact that Sp5- and PoM-projecting neurons have the lowest voltage thresholds for AP generation (Table 1). M1- and L6 M1-projecting pyramidal neurons exhibited both more hyperpolarized resting membrane potentials and higher thresholds for spike initiation (Fig. 6E and Table 1). We hypothesized that the smaller voltage differential in Sp5- and PoM-projecting neurons compared with and L6 M1-projecting neurons accounts for differences in spiking probabilities to M1 input (Fig. 6F). We confirmed this by comparing the voltage difference between resting and AP threshold with spike probability output for all neurons having suprathreshold responses. Neurons with the smallest voltage differences, which tended to be Sp5- and PoM-projecting neurons, had the highest spike probability, whereas those with larger voltage differences had a lower probability of spiking (Fig. 6G), despite many of them receiving larger-amplitude EPSCs. Relationship between the project target, M1 input, and laminar depth. Our findings demonstrate both laminar- and projection target-specific selectivity of M1 projections onto infragranular pyramidal neurons. ANOVA with three factors demonstrated a significant effect of projection target (P.3), a weak effect at trend level of layer (L5, L6: P.8), and no effect of normalized distance from the L4/L5 border (P.99). To disambiguate further the laminar and projection target effects, we compared the present data involving targetidentified neurons with data from our previous study (Kinnischtzke et al. 214) in which unlabeled infragranular pyramidal cells were recorded in unbiased fashion, without projection target indicators. The sample of cells projecting to unknown A Normalized depth from EPSC amp (pa) B EPS Camp (pa) targets overall had a similar distribution of M1 input amplitude by depth as cells with a known target (Fig. 7A). We compared amplitudes of M1 inputs to projection-identified neurons with unidentified neurons located within the same depth range (Fig. 7B). In some cases, it is possible the unlabeled cells are of the same projection type as those they are being compared with; the unlabeled population as a group, however, serves to provide an indication of the amount of M1 input that depth typically receives from M1. For example, M1-projecting neurons, located in at a normalized depth of.17 (5 95% percentile), received the same low-amplitude inputs as cells of an unknown projection target located within the same range of depths ( vs pa, P.56; Fig. 7D). Similarly, M1 inputs to L6 M1- projecting cells, which received on average the strongest M1 input, were equivalent to those of unidentified cells at the same depths ( vs pa, P.86). On the other hand, PoM-, Sp5-, and VPM-projecting populations each received markedly less M1 input than nearby cells having unidentified targets (P.4, P.8, and P.3, respectively). These findings indicate that the strength of M1 input to reciprocally connected M1-projecting infragranular neurons is governed by their laminar (or sublaminar) location. In contrast, the amount of M1 input to at least three classes of subcortically projecting cells (Sp5, PoM, and VPM) is determined based on the projection target. DISCUSSION Here, we combined optogenetic techniques with retrograde tracers to investigate M1 inputs to selected classes of S1 infragranular pyramidal neurons. Retrograde tracers allowed identification of the projection target for each recorded neuron while injection of an AAV containing the gene for ChR2 permitted selective activation of synaptic inputs originating in M1. Determining synaptic connectivity using optogenetic activation of ChR2-expressing terminals is ideal for targeting specifically the synapses of interest, in this case, M1 axons within S1. This approach does result in near-synchronous n.s. n.s. M1 Sp5 PoM L6 M1 VPM Proj. target [-.17] [ ] [ ] [ ] [.5-.75] Norm. depth Fig. 7. M1 inputs contact S1 neurons based on the laminar position or projection target. A: population graph of the amplitude of the M1-evoked EPSC as a function of the depth for each neuron indicating that M1 inputs are particularly strong around the border between L5 and L6 for neurons with known projection targets ( M1-projecting neurons: gray open circles, Sp5-projecting neurons: open diamonds, PoM-projecting neurons: open triangles, L6 M1-projecting neurons: black open circles, and VPM-projecting neurons: open squares) and unknown projection targets (black filled circles). B: comparison of neurons with known axonal projection (gray bars) with cells at the same laminar position but unknown projection (solid bars) demonstrated that Sp5-, PoM-, and VPM-projecting neurons all received significantly less input than neighboring cells at the same depths (P.5, P.1, and P.5 by a two-sample t-test).