Cerebral Cortex Advance Access originally published online on August 24, 2005
Cerebral Cortex 2006 16(6):791-810; doi:10.1093/cercor/bhj023
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Structure of the Excitatory Receptive Fields of Infragranular Forelimb Neurons in the Rat Primary Somatosensory Cortex Responding To Touch
1 School of Biomedical Engineering, Science and Health Systems, Drexel University, Philadelphia, PA, USA, 2 Department of Neurobiology, Drexel University College of Medicine, Philadelphia, PA, USA and 3 Fundación del Hospital Nacional de Parapléjicos para la Investigación y la Integración, SESCAM, Toledo, Spain
Address correspondence to Karen A. Moxon, PhD, Drexel University, School of Biomedical Engineering, Science and Health Systems, Philadelphia, PA 19104, USA. Email: Karen.Moxon{at}drexel.edu.
| Abstract |
|---|
|
|
|---|
We quantitatively studied the excitatory receptive fields of 297 neurons recorded from the forelimb infragranular somatosensory cortex of the rat while touch stimuli were applied to discrete locations on the forelimbs. Receptive fields were highly heterogeneous, but they were regulated, on average, by an underlying spatio-temporal structure. We found the following. (i) Neurons responded with decreasing magnitude and increasing latency when the stimulus was moved from the primary location to secondary locations and to far ispilateral locations of their excitatory receptive fields, displaying smooth transitions from the primary location to secondary locations. (ii) Receptive field patterns revealed functional connectivity between the digits and ventral palm, which did not depend on whether the digits were stimulated dorsally or ventrally. (iii) The structure of the receptive fields (i.e. the neural responses to stimulation of secondary locations compared to the neural responses to stimulation of the primary location), reflected cortical (rather than body) distances. (iv) There was a functional separation between the forepaw and the rest of the forelimb. Namely: if the primary location was in the digits or palm, secondary locations were biased toward the digits and palm; if the primary location was in rest of the forelimb, secondary locations appeared equally distributed over forelimb, digits and palm. (v) More than 40% of neurons extended their receptive field to the ipsilateral forelimb, without any evident spatial organization. Overall, the stimuli evoked
3 times more spikes from secondary responses than from primary responses. These results suggest that a rich repertoire of spatio-temporal responses is available for encoding tactile information. This highly distributed receptive field structure provides the electrophysiological architecture for studying organization and plasticity of cortical somatosensory processing.
Key Words: tactile plasticity forelimb multi-electrode cortical
| Introduction |
|---|
|
|
|---|
The primary sensory cortex is organized in columns, which are subserved by a layered structure within each column (Staiger et al., 2000
Several studies in the auditory (Eggermont et al., 1981
; Suga et al., 1983
), visual (Movshon et al., 1978
; Reid and Juraska, 1991
; McLean and Waterhouse, 1994
; Ringach et al., 1997
) and trigeminal (whisker) somatosensory systems (Simons, 1985
; Simons and Carvell, 1989
; Armstrong-James et al., 1992
; Moore and Nelson, 1998
; Ghazanfar and Nicolelis, 1999
; Zhu and Connors, 1999
; Foffani and Moxon, 2004
) have shown that cells in the infragranular cortex exhibit complex spatiotemporal firing patterns and large excitatory receptive fields. These large excitatory receptive fields may be relevant for sensory exploration of natural environments (for a review, see Ghazanfar and Nicolelis, 2001
). However, little quantitative information is available about the structure of the excitatory receptive field of neurons in the infragranular somatosensory cortex of the rat, other than in the whisker area. In particular, the rat forelimb is of significant interest because rats are capable of performing complex discrimination tasks using their paws and digits, and these behaviors are more analogous to primate hands than are whisking behaviors (Whishaw et al., 1986
; Whishaw and Pellis, 1990
; Whishaw and Gorny, 1994
). In addition, the forelimb region of the somatosensory cortex, rather than the whisker region of the somatosensory cortex, is necessarily chosen in rat models studying the effects of peripheral limb damage or spinal lesion on cortical receptive field reorganization (Kaas et al., 1983
; Kaas and Collins, 2003
; Moxon and Foffani, 2003
). It is therefore important to examine quantitatively the spatial and temporal structure of the somatosensory receptive field of neurons in the forelimb somatosensory region of the cortex responding to touch.
In the present study, microelectrode arrays were implanted bilaterally into the forepaw region of the infragranular primary somatosensory cortex of the rat to record the response of ensembles of single neurons during touch stimulation of the forelimbs under light anesthesia. To describe the structure of the receptive fields consistently with studies in the whisker cortex (Ebner and Armstrong-James, 1990
; Ghazanfar and Nicolelis, 1999
; Panzeri et al., 2001
; Brecht and Sakmann, 2002
; Fox et al., 2003
), we defined the primary location as the location where a neuron exhibited the maximal response, secondary locations as all other locations on the forelimb, contralateral to the neuron, that generated an excitatory response, and far ipsilateral locations as all the locations on the ipsilateral forelimb that generated an excitatory response. Response magnitudes and latencies were used as the main quantitative parameters to compare the response of cells to stimulation of their primary location, secondary locations and far ipsilateral locations. By simultaneously recording ensembles of neurons from 32 separate electrodes equally distributed in the left and right hemispheres, and by presenting the same standardized stimuli at 20 different locations on the forelimbs we were able to objectively assess the relative strength of responses in different cortical columns/segregates, to quantify the activation of neurons from stimuli delivered to multiple forelimb locations, and to quantitatively evaluate the spatial and temporal structure of the excitatory contralateral and ipsilateral receptive fields of forelimb infragranular neurons in the rat primary somatosensory cortex responding to touch stimuli.
| Materials and Methods |
|---|
|
|
|---|
Recordings were made from seven LongEvans rats (240300 g) and all procedures were approved by the Institutional Animal Care and Use Committee at Drexel University and followed NIH Guidelines. To record from large numbers of neurons simultaneously and to minimize the effect of anesthesia on the somatosensory responses of the cells, two arrays of 16 microwire electrodes each, one for each hemisphere of the brain, were chronically implanted under Nembutal anesthesia (50 mg/kg). Animals were allowed to recover (approximately one week) and then the responses of neurons to touch stimulation of the forepaw and forelimb were quantified under light Nembutal anesthesia (induction dose: 35 mg/kg, stage III-2; Friedberg et al., 1999
Surgical Procedures to Implant Microelectrodes
Animals were anesthetized and placed in a stereotaxic apparatus for surgery (Cartesian Research, Sandy, OR). Rectangular shaped craniotomies with coordinates 0.5 mm anterior to bregma, 3.5 and 4.5 mm lateral to 1.5 posterior to bregma, and 3.0 and 4.0 mm lateral (atlas of Paxinos and Watson) (Chapin and Lin, 1984
) were performed bilaterally over the somatosensory forepaw (palm and digits) area to accommodate two electrode arrays (one per side), each consisting of two rows of eight 50 µm Teflon-coated stainless steel microwires (NB labs, Dennison, TX). The spacing between microwires before implant was 200 µm. The arrays were oriented so that rows ran from rostral to caudal. As each electrode array was implanted, neural activity was continuously monitored (for details, see Single Neuron Discrimination) and amplified through auditory speakers. The forepaw was gently tapped to elicit somatosensory responses and to ensure that electrodes were properly placed in the forepaw region. When the characteristic large amplitude of layer V neurons was recorded on the majority of electrodes, the electrode was cemented in place. The connectors were surrounded with dental cement to create an electrode cap that formed a base on which to attach a recording headstage during subsequent recording sessions.
Single Neuron Discrimination
Single neuron discrimination was done 710 days after the implantation surgery using the same methods used to describe layer V neurons in the barrel field cortex to allow direct comparison (Ghazanfar and Nicolelis, 1999
; Foffani and Moxon, 2004
; Foffani et al., 2004
). In brief, rats were anesthetized with low doses of Nembutal to minimize interference of the anesthesia on the neural recordings (Friedberg et al., 1999
) but sufficient to immobilize the rat. Stable levels of light anesthesia were maintained at different times within the same session by giving small supplements when the rat consistently responded to tail-pinch. Signals were amplified and filtered using a multi-neuron acquisition system (Plexon Inc., Dallas, TX) and the resulting signals were displayed on an oscilloscope and amplified through loudspeakers to aid in online neuronal spike sorting from all 32 channels (Wheeler, 1999
).
Receptive Field Maps
One receptive field map was performed on each animal in the following way. Since our goal was to investigate the spatial and temporal structure of the receptive field in response to touch, 10 discrete locations were chosen for stimulation on each forelimb (Fig. 1A). These locations included one spot on each of the five digits, labeled (1) D1, (2) D2, (3) D3, (4) D4 and (5) D5. D3 and D4 were stimulated on the dorsal surface while digits D1, D2 and D5 were stimulated on the ventral surface. In addition, five other arbitrary but consistent locations across all animals were stimulated. These locations included a spot on (6) the dorsal paw (PAW), (7) the ventral palm (PLM), (8) the wrist (WR), (9) distal forelimb (DFL) and (10) proximal forelimb (PFL). The wrist, distal forelimb and proximal forelimb were stimulated dorsally. Assumptions and implications of using this sparse stimulation protocol will be fully explained in the discussion. During stimulation, the paw and limb were placed on their side, digit 1 facing up, so that all locations were easily accessible by the stimulator. Each of the above locations was consecutively stimulated 100 times at 0.5 Hz with a fine tipped metal probe 1 mm in diameter. To be consistent with previous studies in the whisker system and ensure that only tactile receptors at the sight of contact were activated, the metal probe was controlled through a piezoelectric element actuated by a Grass stimulator (Model S48), which delivered squared-pulse stimuli (duration: 100 ms, frequency: 0.5 Hz), similar to previous studies (Chapin, 1986
; Foffani et al., 2004
). The tip of the metal probe moved 0.5 mm in response to the square-pulse stimuli. To control the magnitude of the stimuli at each location, the metal probe was first positioned on the skin, ensuring contact but no visual indentation under x10 magnification. The metal probe was then moved 0.5 mm away from the skin and the stimulation was started. The effect of the stimulus was viewed under x10 magnification to ensure no movement of the digits or limb. These stimulus properties and the relatively large distance between the locations stimulated make the possibility of stimulus spread across locations extremely unlikely. All locations were stimulated within the same recording session to ensure that the same neurons were recorded in response to stimulation of all locations. All 100 stimuli were given to a location and then the stimulator was moved to the next location. There was no randomization of stimuli. The frequency of stimulation of 0.5 Hz corresponds to twice the interstimulus interval previously shown not to influence subsequent responses (Chapin, 1986
). The Grass stimulator simultaneously sent pulses to the data acquisition system for precise timing of the stimulus onsets. The waveforms and action potential times of all discriminated neurons were recorded during the receptive field map and the data were stored in NeuroExplorer (Nex Technologies, Littleton, MA, version 2.66). For every location, peri-stimulus time histograms (PSTHs, 1ms binsize) of all the neurons were calculated using Nex functions and exported to Matlab (version 6.5, The Mathworks) for further analysis.
|
Quantitative Measures of Receptive Fields
In order to identify significant responses in the PSTHs, three tests were performed for every neuron and for every location (Fig. 1B): (i) a threshold was set as the average background activity of the neuron (evaluated from 100 to 5 ms before the stimulus) plus 3 SD, and the first and the last significant bin (1ms binsize) that exceeded the threshold in a window between 5 and 90 ms after the stimulus were identified; (ii) at least three bins had to be over the threshold; and (iii) the response between the first and the last significant bin had to be significantly greater than the background activity (non-paired t-test, P < 0.001). For every significant response, four parameters were extracted from the PSTH: (i) the response magnitude, defined as the integral of the PSTH between the first and the last significant bin (i.e. probability of spike per stimulus); (ii) the peak response, defined as the maximum probability of spike per bin; (iii) the first bin latency; and (iv) the peak latency, defined as the time intervals between the stimulus onset and the first significant bin or the peak, respectively. The rationale for setting the threshold as 3 SD above background was to minimize the false identification of significant responses. A posteriori, the peak responses were 18.68 SD above background on average and >4 SD above background in 97.18% of 1794 significant responses. These values suggest that the threshold choice (on average 0.017 spikes/bin) guaranteed the identification of virtually all physiological responses. For each neuron, the primary location of its excitatory receptive field was defined as the location that generated the greatest response magnitude. All the other locations where the neuron showed a significant excitatory response on the side of the body contralateral to the electrode were defined as secondary locations. The locations on the ipsilateral side of the body where the neuron showed a significant excitatory response were defined as far ipsilateral locations. An additional parameter, the normalized response magnitude, was calculated for every neuron and for every location as the ratio between the response magnitude and the response magnitude of the primary location. Therefore, the normalized response magnitude of the primary location is equal to 1.0 by definition. The discrete receptive field size was calculated as the total number of locations where a neuron exhibited a significant response. Note that this discrete definition of receptive field size (number of locations) is dependent on the stimulation protocol and is comparable with the discrete definition typically used in whisker studies (number of whiskers) but not with the continuous definition conventionally employed in non-whisker studies (mm2 of skin). To reduce errors due to the finite number of locations stimulated, neurons were included in the analyses only if the response magnitudes to stimulation of their primary location were strong enough to allow a reliable estimation of the parameters for the primary and secondary locations. Based on preliminary analyses, a threshold for the primary response magnitude was set at 0.2 spikes/stimulus. Once this minimal response magnitude for stimulation of the primary location was respected, responses to stimulation of secondary locations were considered significant according to the three criteria above, without additional thresholds (i.e. secondary response magnitudes were allowed to be <0.2 spikes/stimulus).
Comparison of Primary Location, Secondary Locations and Far Ipsilateral Locations
In order to compare the responses to stimulation of the primary location, secondary locations and far ispilateral locations of the neurons' receptive fields, and to test whether there were differences between cells whose primary location was in the glabrous surface and cells whose primary location was in the hairy surface of the forelimb, we used a two-way independent-measures analysis of variance (ANOVA) that was separately performed for each of the following variables: response magnitude, peak response, first bin latency and peak latency. Each significant response was conservatively considered as an independent sample. The first factor of the ANOVA was the position of the stimulus within the receptive field, with three levels: primary location, secondary locations or far ipsilateral locations. The second factor was the identity of the cell's primary location, with two levels: glabrous skin (i.e. grouping digits 1, 2, 5 and ventral palm together) or hairy skin (i.e. grouping digits 3, 4, dorsal paw, wrist, distal forelimb and proximal forelimb together). Sheffe's test was employed for post-hoc comparisons.
Transitions from the Primary Location to the Secondary Locations
In order to test whether receptive fields were characterized by an abrupt or smooth transition of response magnitudes and latencies when stimuli were moved from the primary location to the secondary locations, two analyses were performed. The transition of response magnitudes from the primary location to the secondary locations was examined by quantifying the discrete receptive field size as a cumulative function of the normalized response magnitude, i.e. counting the total number of locations where a neuron exhibited a significant response with normalized response magnitude greater than a given value (Nicolelis and Chapin, 1994
). The transition of response latencies from the primary location to the secondary locations was examined by quantifying the first bin latency of the responses as a function of the normalized response magnitude. The dependence between the discrete receptive field size of neurons and the location of their primary location was tested through a one-way ANOVA. The main factor was the identity of the primary location, the dependent variable was the discrete receptive field size and each neuron was considered as an independent sample. A non-paired t-test was also used when the identity of the primary location was categorized into glabrous versus hairy skin.
Receptive Field Patterns
Receptive field patterns were characterized by assigning, for each neuron, a 1 to each location eliciting a significant response when stimulated, and a 0 to the remaining locations. The receptive field pattern of a neuron was then defined by the combination of locations that elicited significant responses when stimulated. Neurons were therefore considered to possess the same receptive field pattern if they significantly responded to stimulation of the same locations.
Receptive Field Structure
To evaluate the structure of the receptive fields, neurons were grouped based on their primary location. Every secondary location for each group of neurons with the same primary location was characterized by two parameters: (i) the response likelihood, defined by counting the number of neurons responding to that secondary location divided by the number of neurons with the same primary location; and (ii) the average normalized response magnitude for that secondary location. Note that both the response likelihood and the average normalized response magnitude are numbers between 0 and 1, and that the primary location by definition corresponds to the point [1,1]. A quantitative evaluation of these spatial properties was performed by clustering the responses, as defined by response likelihood and average normalized response magnitude, based on the average Euclidean distance between pairs of objects in different clusters weighted by the number of objects in the cluster (Sneath and Sokal, 1973
; weighted pair-group average in Statistica). This is identical to the centroid method since the space defined by the system is two-dimensional with limits of 0 and 1 in both dimensions. An ANOVA was performed on each of the clusters to evaluate the most significant clustering for each primary location (i.e. the grouping of responses by stimulus location that generated the greatest F-value) using Statistica (Statsoft, Tulsa OK).
Throughout the text, results are means ± SD.
Histology
To examine the position of the electrode, post-mortem brain slices were Nissl stained to mark cell bodies. After the final recording session, rats were anesthetized and a stimulating current was placed down the electrode to mark the electrode tip (60 µA for 45 s). Each rat was transcardially perfused with 0.1 M phosphate buffer then 4% paraformaldehyde in 0.1 M phosphate buffer. The brain was removed from the skull and electrodes were extracted from the brain. The brain was stored in 4% paraformaldehyde in 0.1 M phosphate buffer overnight then transferred to a 30% sucrose solution. After 5 days, the brain was removed from solution, embedded in Tek solution and frozen. The brains were cut into 40-µm-thick coronal sections using a cryostat and collected on slides. The positions of the tips of the electrodes were easily visualized due to the hole created by the microstimulation (Fig. 1C).
| Results |
|---|
|
|
|---|
We analyzed 297 neurons that responded to touch stimulation of the forelimb, recorded from seven rats bilaterally implanted in the forepaw somatosensory cortex. Figure 2 shows typical neurophysiological responses of 6 representative neurons simultaneously recorded from one rat, with raster plots and peri-stimulus time histograms (PSTHs) corresponding to all 10 locations stimulated on the contralateral forelimb and three representative locations stimulated on the ipsilateral forelimb. The results described in the next sections are derived from the quantitative analysis of 5940 PSTHs (297 neurons x 20 locations), of which 1794 presented significant responses.
|
Distribution of Receptive Field Primary Locations
The forepaw region was targeted and, accordingly, the primary location for most neurons was in the ventral palm or the digits, with fewer neurons having primary locations in dorsal paw, wrist, distal forelimb or proximal forelimb (Table 1, Fig. 3). More precisely, 77% of the cells had their receptive field primary location on the glabrous surface of the forepaw (digits 1, 2, 5 or ventral palm), whereas 23% of the cells had their primary location on the hairy surface of the forepaw (digit 3, 4, dorsal paw, wrist, proximal forelimb and distal forelimb). The number of cells sharing the same primary location was well correlated between the two sides of the brain (Spearman: R = 0.81, P = 0.0044, n = 10; see Table 1), confirming the experimental symmetry of our implants and the functional symmetry of the distribution of receptive field primary locations. There was no positive correlation (Spearman, significance level: P < 0.01; see Table 1) between the number of cells sharing the same primary location and the number of cells sharing any other primary location, which was not surprisingly given the unpredictability of the final electrode position using microwires. The number of cells sharing the same primary location was instead well correlated with the total number of significant responses (primary or secondary) elicited at that location (Spearman: R = 0.95, P = 0.00003, n = 10; data not shown).
|
|
|
Comparison of Primary Location, Secondary Locations and Far Ipsilateral Locations
Response Magnitudes
We quantified for each neuron the response magnitude, (i.e. probability of spike per stimulus) and the peak response, (i.e. maximum probability of spike per bin at 1ms binsize) of every significant response (Table 2). As expected, there was a significant difference between responses to stimulation of the primary location, secondary locations or far ipsilateral locations, using either response magnitude or peak response (ANOVA: P < 0.000001). Neurons responded to stimulation of their primary location (n = 297 responses), secondary locations (n = 1295 responses) and far ispilateral locations (n = 202 responses) with a response magnitude of 0.74 ± 0.69, 0.47± 0.38 and 0.29 ± 0.29 spikes respectively, and with a peak response of 0.10 ± 0.07, 0.07 ± 0.04 and 0.04 ± 0.04 spikes/(bin x stimulus). No significant differences were found between the responses of neurons whose primary location was on the glabrous surface of the forelimb and the responses of neurons whose primary location was on the hairy surface of the forelimb, using either response magnitudes (P > 0.49) or peak responses (P > 0.27). There was an overall significant correlation between response magnitude and peak response to stimulation of either the primary location (Pearson: r = 0.78, P < 0.0001, n = 297), the secondary locations (r = 0.87, P < 0.0001, n = 1295), or the far ispilateral locations (r = 0.88, P < 0.001, n = 202), which suggested that these two parameters were equivalent measures of the cell's response. Because the response magnitude represented the total response probability of the cell to the stimulus, we used this measure, rather than peak response, for subsequent analyses. Overall, in the 297 neurons recorded, the stimuli evoked
3 times more spikes from secondary responses [(1295 responses x 0.47 spikes/stimulus + 202 responses x 0.29 spikes/stimulus) x 100 stimuli/response = 66723 spikes] than from primary responses (297 responses x 0.74 spikes/stimulus x 100 stimuli/response = 21978 spikes).
|
Response Latencies
To compare the temporal properties of the neural responses, for every neuron we quantified the first bin latency (i.e. latency of the first significant bin of the response) and the peak latency (i.e. latency of the response peak) of every significant response (Table 2). There was a significant difference between responses to stimulation of the primary location, secondary locations or far ipsilateral locations, using either first bin latency or peak latency (ANOVA: P < 0.00001). Neurons responded to stimulation of their primary location, secondary locations and far ipsilateral locations with a first bin latency of 14 ± 5, 17 ± 7 and 23 ± 12 ms, respectively, and with a peak latency of 22 ± 6, 23 ± 8 and 28 ± 12 ms. No significant differences in response latencies were found between cells whose primary location was on the hairy skin and cells whose primary location was on the glabrous skin when the primary location or the secondary locations were stimulated (Sheffe: P > 0.7). These results suggest that the latency of a contralateral response depended strongly on the position of the stimulus within the neuron's excitatory receptive field.
Transitions from the Primary Location to the Secondary Locations
Given our sparse stimulation protocol, one might expect that the response magnitude of a neuron would abruptly decrease and the response latency would abruptly increase when a stimulus is moved from the primary location to the secondary locations. However, this was not the case. The transition of response magnitude from the primary location to the secondary locations was quantified by examining the discrete receptive field size (i.e. the number of locations where a neuron exhibited a significant response) as a cumulative function of the normalized response magnitude. As the normalized response magnitude was allowed to decrease from 1 (i.e. the primary location) to 0 (i.e. all locations), the discrete receptive field size smoothly, not abruptly, increased (linear fit: r = 0.997, P < 0.0001, n = 11; see Fig. 4A). This occurred despite the sparseness of the locations stimulated in this study. In other words, there was a linear relationship between normalized response magnitude and discrete receptive field size. This suggests that the response magnitude of cells to stimulation of any location in their excitatory receptive field is uniformly distributed between 0 and the maximal response (i.e. the response to stimulation of the primary location). The average discrete receptive field size considering all significant responses was 6.0 ± 3.4 locations, with some cells displaying small receptive fields of one or two locations and the majority of them having large receptive fields that included several secondary locations on the contralateral forelimb and often extended to the far ipsilateral locations. The discrete receptive field size did not depend on the identity of the primary location (ANOVA, P = 0.22). In particular, the discrete receptive field size of cells whose primary location was on the glabrous skin (digits 1, 2, 5, ventral palm) was not significantly different (unpaired t-test, P = 0.94) from the discrete receptive field size of cells whose primary location was on the hairy skin (digits 3, 4, dorsal paw, wrist, distal forelimb, proximal forelimb) (Fig. 4B). On average, 40.7% of the neurons had significant responses to stimulation of at least one far ipsilateral location.
|
The transition of response latency from the primary location to the secondary locations was quantified by examining the first bin latency as a function of the normalized response magnitude (Fig. 4C). In agreement with previous studies, on average the primary location had the shortest response latency. In addition, there was a significant negative correlation between normalized response magnitude and first bin latency (Pearson: r = 0.32, P < 0.001, n = 1794; see Fig. 4C). In other words, as the normalized response magnitude was allowed to decrease from 1 to 0, the first bin latency smoothly, not abruptly, increased.
Receptive Field Patterns
Figure 5 shows the receptive fields of all neurons recorded from a representative animal (Rat 1). This animal had cells whose primary locations were predominantly on digit 5 (12 cells) and on the ventral palm (13 cells). All of the neurons whose primary location was on digit 5 also responded to stimulation of the ventral palm, and eight of the neurons whose primary location was on the palm also responded to stimulation of digit 5, suggesting a possible functional connectivity between digit 5 and the ventral palm, which is consistent with the adjacency between the two locations on the cortical map (Fig. 3). It is important to note that virtually none of the receptive fields were duplicated: of the 13 cells whose primary location was on the ventral palm and of the 12 cells whose primary location was on digit 5, every cell's receptive field had a unique pattern of response. This suggests that receptive fields were highly heterogeneous across cells, even from the same rat and even for cells with the same primary location.
|
In order to quantify this heterogeneity in our entire dataset, we defined the receptive field pattern for each neuron simply as the combination of locations that elicited significant responses when stimulated. Neurons were therefore considered to possess the same receptive field pattern if they significantly responded to stimulation of the same locations, thus without considering the variability in response magnitude. Still, the 297 neurons analyzed exhibited 121 different receptive field patterns (Fig. 6). Considering only contralateral responses, the most frequent pattern was found in 28 cells (9.4% of the sample), which responded to all contralateral locations except proximal forelimb (i.e. [19]). The subsequent most frequent patterns were small variations of the above pattern. Namely: 22 cells responded to [17], 22 other cells responded to [15, 7] and 17 cells responded to [18]. Overall, 115 neurons (39% of the sample) exhibited significant responses at least to all five digits and the ventral palm (i.e. [15, 7]). Importantly, only nine neurons responded to digits 1, 2 and 5 (stimulated ventrally) without responding to digits 3 and 4 (stimulated dorsally), and only five neurons responded to digits 3 and 4 without responding to digits 1, 2 and 5. These results (i) display the high heterogeneity of receptive field patterns across cells and (ii) suggest the presence of a general non-specific connectivity between the digits and ventral palm that does not depend on whether the digits are stimulated dorsally or ventrally.
|
Receptive Field Structure
To quantitatively examine whether the above patterns depended on the overall receptive field structure, neurons were grouped across rats based on their primary location. For each population of neurons with same primary location, two measures were employed to define the overall response at each secondary location: (i) the average normalized response magnitude, i.e. the normalized response magnitude averaged over the neurons that exhibited a significant response to stimulation of that secondary location; and (ii) the response likelihood, i.e. the fraction of neurons that exhibited a significant response to stimulation of that secondary location. Average normalized response magnitude and response likelihood were significantly correlated (Pearson: r = 0.63, P < 0.001), and were studied together with a cluster analysis. The detailed receptive fields for all 297 individual cells are shown in Figure 7, which provided the data for the quantitative cluster analysis shown in Figure 8.
|
For populations whose primary location was in the ventral palm or in the digits, the responses to stimulation of secondary locations could be separated into two significant clusters: responses to stimulation of the digits or palm clustered together and were significantly greater than the responses to the other forelimb locations stimulated (ANOVA, P < 0.05). More explicitly, if the primary location was in the ventral palm or in the digits, then both the average normalized response magnitudes and response likelihoods were quite large (typically >0.5, Fig. 8A) when the stimuli were delivered to secondary locations in the digits or ventral palm, and much smaller (typically <0.5, Fig. 8A) when the stimuli were delivered to secondary locations in the dorsal paw, wrist, distal forelimb and proximal forelimb. The receptive field structure did not appear to consistently discriminate between the median nerve innervated glabrous skin (digits 1 and 2, Figs 7A and 8A), the ulnar territory (digit 5, Figs 7B and 8A), the glabrous palm (ventral palm, Figs 7C and 8A) and the hairy skin (digits 3 and 4, Figs 7D and 8A). Of particular note, for cells whose primary location was in digits 1, 2, 5 or in the ventral palm, digits 3 and 4 which were contacted on the dorsal hairy surface clustered with the other digits and the palm which were contacted on the ventral glabrous surface and not with the other dorsal locations (Fig. 8A). Reciprocally, neurons whose primary location was in digits 3 or 4 exhibited greater responses to stimulation of the ventral palm (average normalized response magnitude = 0.58 ± 0.29 and 0.77 ± 0.13 respectively; response likelihood = 0.68 and 0.73) than to the dorsal paw (average normalized response magnitude = 0.34 ± 0.18 and 0.45 ± 0.19; response likelihood = 0.42 and 0.27). These results show that the receptive field structure does not have a direct relation with the peripheral distribution of somatosensory nerves, but does exhibit a more consistent relation with the cortical distribution of receptive field primary locations. In fact, according to the cortical map in Figure 3, digits 3 and 4 are closer to the other digits and to the ventral palm compared to the dorsal paw (in the map there is no separation between the dorsal and ventral representations of the digits). More generally, for neurons whose primary location was in the digits or in the ventral palm, the two closest secondary locations in the (average normalized response magnituderesponse likelihood) plane were, on average, the closest adjacent locations in the cortical plane (with one exception; see below) (Fig. 8A). Namely, digit 2 and the ventral palm were the two closest secondary locations for cells whose primary location was in digit 1; digit 1 and the ventral palm for cells whose primary location was in digit 2; digits 2 and 4 for cells whose primary location was in digit 3; digit 3 and the ventral palm for cells whose primary location was in digit 4; ventral palm and digit 1 (exception) for cells whose primary location was in digit 5; and digits 1 and 5 for cells whose primary location was in the ventral palm. According to the same principle, both the average normalized response magnitude and the response likelihood to stimulation of the dorsal paw, wrist, distal forelimb and proximal forelimb tended to decrease as the cortical distance from the cell's primary location increased. In other words, locations that are represented in progressively more distant regions in the brain are progressively less likely to be represented in the receptive field structure. Therefore, despite the heterogeneity of receptive field patterns, the average receptive field structure of forelimb infragranular neurons is consistent with the somatotopic organization of the primary somatosensory cortex.
For populations of cells whose primary location was in the dorsal paw, wrist or distal forelimb, the two clusters above could not be found and these neurons were equally likely to respond to all secondary locations (Fig. 8B). Only for cells whose primary location was in the dorsal paw, the two closest secondary locations were the closest adjacent locations in the cortical plane, namely digit 2 and ventral palm (Fig. 8B). It should be noted that few cells had their primary location in the dorsal paw (n = 5), wrist (n = 13) or distal forelimb (n = 7), but the two clusters were not found even when these cells were combined together (n = 25, data not shown). These results suggest two important differences when comparing secondary responses to stimuli delivered to the forepaw (digits and palm), with secondary responses to stimuli delivered to the rest of the forelimb (from the dorsal paw to the shoulder). First, more cells in the somatosensory cortex exhibit secondary responses, and second, each individual cell has a greater probability of responding when a stimulus is delivered to the forepaw compared to the rest of the forelimb.
For reasons of symmetry, one might expect that a neuron would have the highest probability of exhibiting a significant response to the far ipsilateral locations when the homologous location of the contralateral primary location is stimulated on the ispilateral side. On the contrary, the identity of the cell's primary location did not influence the organization of its response to stimulation of far ipsilateral locations. In fact, these responses were distributed across all locations stimulated (Fig. 9A) and in only 12% of the cases which is approximately what one would expect by chance the location of an ipsilateral response coincided with the location of the contralateral primary location (Fig. 9B).
|
| Discussion |
|---|
|
|
|---|
We quantitatively studied the excitatory receptive field structure of forelimb neurons in the infragranular somatosensory cortex of the rat in response to touch. Receptive fields were large and their patterns were highly heterogeneous across neurons, even for cells with the same primary location. However, the quantitative evaluation of magnitudes and latencies of the responses of neurons to stimulation of their primary location, secondary locations and far ipsilateral locations revealed an overall receptive field organization. Neurons responded with decreasing magnitude and increasing latency when the stimulus was moved away from the primary location to the secondary locations and to the far ispilateral locations. Despite the discreteness of the stimulation protocol, these transitions of response magnitudes and latencies from the primary location to the secondary locations were smooth. In addition, the identity of the primary location did not influence the discrete receptive field size of neurons, but did affect the structure of the receptive field, suggesting a functional separation between the forepaw and the rest of the forelimb. Namely: if the primary location was in the digits or in the ventral palm, the secondary locations were biased toward the digits and ventral palm, which clustered together; if the primary location was in the rest of the forelimb, the secondary locations appeared equally distributed over the entire forelimb, including digits and palm. On average, the structure of the receptive field within individual neurons was largely consistent with the overall somatotopic organization of the somatosensory cortex (i.e. the neural responses to stimulation of a secondary location in the receptive field decreased when the cortical distance between that secondary location and the primary location increased), suggesting that the receptive fields themselves are somatotopically organized. Finally, >40% of the recorded cells extended their receptive field to the ipsilateral forelimb, without any evident spatial organization. Taken together, these results suggest (i) that even a light punctuate stimulus can robustly activate large regions of the forelimb infragranular somatosensory cortex; (ii) that the majority of spikes sent by the infragranular cortex to intracortical and subcortical structures originates by secondary rather than primary responses; and (iii) that there is a rich repertoire of spatio-temporal responses distributed across the somatosensory cortex which is available for encoding information for tactile discrimination.
Studying Excitatory Receptive Fields
This work is focused on the excitatory responses to touch. Therefore, we did not explicitly consider inhibition, which is difficult to study using extracellular recordings when the background firing rates of neurons are very low. Our results are therefore exclusively related to the excitatory receptive fields of forelimb neurons in the rat infragranular somatosensory cortex and their full functional interpretation should assume the presence of spatial and temporal inhibition (Laskin and Spencer, 1979
; Simons, 1985
; Brumberg et al., 1996
; Xing and Gerstein, 1996
; DiCarlo and Johnson, 2002
; Derdikman et al., 2003
; Wirth and Luscher, 2004
).
Using a Sparse Stimulation Protocol
In this study we adopted a sparse stimulation protocol, delivering punctuate stimuli to 20 locations bilaterally distributed on the forelimbs. The rationale for this choice has three main bases: (i) to parallel as much as possible the discrete stimulation protocol necessarily employed when studying the whisker cortex (Ghazanfar and Nicolelis, 1999
; Panzeri et al., 2001
; Foffani and Moxon, 2004
); (ii) to sample the entire likely extent of the excitatory receptive field of infragranular forelimb neurons (Chapin, 1986
; Foffani et al., 2004
); and (iii) to guarantee that cortical distances (rather than body distances) between adjacent locations were qualitatively similar and comparable to cortical distances between adjacent whiskers. Importantly, these inter-location cortical distances are consistent with the distances observed between columns/segregates in the forelimb somatosensory cortex of cats and monkeys (Favorov et al., 1987
; Favorov and Whitsel, 1988a
,b
; Favorov and Diamond, 1990
), and between forepaw barrels in the somatosensory cortex of the rat (Waters et al., 1995
). With these premises, our estimates of magnitudes and latencies of responses to stimulation of the receptive field secondary locations relative to responses to stimulation of the receptive field primary location are dependent on the assumption that the identity of the primary location does not significantly change when stimuli are delivered within a skin area corresponding to the same cortical segregate, which was experimentally verified in the cat (Favorov and Diamond, 1990
). Still, in some cells we may have missed the true primary location due to the limited number of locations stimulated.
Nonetheless, our measure of discrete receptive field size should be dimensionally interpreted in the discrete cortical space (i.e. number of segregates) rather than in the continuous body space (i.e. mm2 of skin). For example, a discrete receptive field size of six locations means that a neuron is excited by afferent input that primarily project to (at least) six cortical segregates. Our result that the discrete receptive field size does not depend on the identity of the primary location, therefore, is not in contradiction with the classical concept of magnification (Mountcastle, 1957
; Wang et al., 1997
; Brown et al., 2004
) and suggests that at the infragranular level each cortical segregate receives an approximately equal amount of convergent input from the periphery, which is consistent with informational and energetic optimization principles regulating the development of cortical somatotopic maps (Durbin and Mitchison, 1990
; Montague, 1996
; Van Essen, 1997
; Plumbley, 1999
). Our result that the receptive field structure does depend on the identity of the primary location, suggests that the equal amount of convergent input can be differently distributed in different segregates, according to functional needs. The smooth transitions in response magnitudes and latencies from the primary location to the secondary locations revealed by our discrete stimulation protocol suggest that the intra-segregate heterogeneity of receptive field patterns (Favorov et al., 1987
; Favorov and Whitsel, 1988a
,b
; Favorov and Diamond, 1990
) possesses an underlying spatio-temporal organization that could be critical to guarantee spike-timing precision for distributed strategies of cortical somatosensory processing (Foffani et al., 2004
). Finally, the sparse stimulation protocol is conservative for the finding that most responses in the infragranular somatosensory cortex arise from stimulation of secondary locations, in the sense that the large inter-location distances would at most bias the data toward an underestimation of the overall contribution of secondary responses. All the above considerations support the conclusion that the entirety of the receptive field not only the primary location plays an essential role for the organization of cortical somatosensory processing.
Activating Touch Mechanoreceptors
Stimuli were optimized for activating touch mechanoreceptors. Accordingly, the majority of cells had their primary location in the glabrous skin of the forepaw and fewer neurons had their primary location in the hairy skin. This bias could be due to two reasons, which are not mutually exclusive: (i) the specific location targeted during surgery (i.e. the forepaw) and (ii) the differences in tactile mechanoreceptors density across the forelimb. Interestingly, the latter possibility is consistent with the distal-to-proximal gradient of peripheral innervation density described in cats in classical studies (Mountcastle, 1957
) and more recently by Brown and co-workers (Wang et al., 1997
; Brown et al., 2004
), and with the differences in the distribution of cutaneous mechanoreceptors observed in humans by Edin and co-workers (Edin, 1990
, 1992
; Edin and Abbs, 1991
; Essick and Edin, 1995
). The receptive field structure is therefore likely to reflect, at least in part, these differences in the distribution of cutaneous mechanoreceptors. Also the different physical characteristics of the skin area tested could have influenced our results. However, the overall excitability and discrete receptive field size of neurons did not depend on whether their primary location was on the glabrous or hairy surface of the forelimb, and in both the receptive field patterns and the receptive field structure digits 3 and 4 clustered with the glabrous forepaw locations (digits 1, 2, 5 and ventral palm) even though they were contacted on the hairy surface of the skin. From a functional point of view, this clustering is not surprising because in rodents the dorsal regions of the digits are likely to be stimulated together with the glabrous digit tips when the animal manipulates objects such as food or explores its environment with its forepaws. From the point of view of cortical anatomy, the above clustering is not surprising either. In fact, the cortical representation of the digits does not have a clear dorsalventral separation and, accordingly, in the cortical map digits 3 and 4 are closer to the other digits and to the ventral palm compared to the dorsal paw (see Fig. 3). More generally, our findings suggest that the receptive field structure of infragranular neurons in the primary somatosensory cortex reflects cortical rather than body distances, which indirectly relates our results to the distribution of cutaneous mechanoreceptors through the well-known proportionality between central magnification and peripheral innervation density (Mountcastle, 1957
; Wang et al., 1997
; Brown et al., 2004
). Because cortical distances were based on the map of the input layer (layer IV) (Chapin and Lin, 1984
), our results are also consistent with the columnar organization of the neocortex (Mountcastle, 1997
). However, it should be remarked that other kinds of stimuli and different attentional states (e.g. awake versus anesthetized) will likely produce different response patterns and that our results are specific to the infragranular layers of the rat primary somatosensory cortex. In this regard, a recent study in the whisker cortex confirmed the layer specificity in the awake state and showed that the anesthetized state is a good model for passive unexpected stimuli (Krupa et al., 2004
).
Comparison with the Trigeminal System
Despite the evident differences in the organization of peripheral receptors between the forelimb and the whisker pad, a number of similarities can be found when comparing our data with the results obtained in the infragranular whisker cortex. The distribution of surround responses, latencies and even discrete receptive field sizes are remarkably similar (Simons, 1985
; Simons and Carvell, 1989
; Armstrong-James et al., 1992
; Moore and Nelson, 1998
; Ghazanfar and Nicolelis, 1999
; Zhu and Connors, 1999
). Moreover, the ipsilateral responses we observed in the forelimb, consistent with previous studies (Zarei and Stephenson, 1996
), are also found in the whisker cortex (Shuler et al., 2002
; Manns et al., 2004
). It is not clear whether the ipsilateral responses observed in this paper originate below the thalamus (Erzurumlu and Killackey, 1980
), from thalamocortical inputs or from callosal projections, either from secondary somatosensory cortex (Wise and Jones, 1976
; Killackey and Belford, 1979
) or directly from the opposite primary somatosensory cortex through the septa (Olavarria et al., 1984
; Shuler et al., 2001
; Fox et al., 2003
). Whatever pathway is involved, our results show that the transmission of somatosensory information to the forelimb ipsilateral cortex is slower than the transmission to the contralateral cortex.
Ispilateral responses also provided a remarkable difference between the forelimb system and the whisker system. Namely, contralateral and ipsilateral receptive fields are spatially homologous in the rat whisker system (Pidoux and Verley, 1979
), whereas we did not find any spatial correspondence between the contralateral primary location of the receptive field and the locations that elicited forepaw, or even forelimb, ipsilateral responses. The specificity of different cortical layers likely plays a critical role in the organization of ipsilateral responses. Nevertheless, the fact that ipsilateral receptive fields are present in both the forelimb and whisker cortices suggests that they might play important roles in bilateral integration of somatosensory information (Shuler et al., 2001
, 2002
).
The matrix organization of the whisker pad is present at the cortical level so that external distances are geometrically transformed, although distorted, in the cortical map. The forelimb cortex has the geometrically more complex task of transforming a non-Euclidean body surface into an almost-Euclidean cortical map. As previously discussed, this transformation has been described in terms of segregates (Favorov et al., 1987
; Favorov and Whitsel, 1988a
,b
; Favorov and Diamond, 1990
) and forepaw barrels (in layer IV; Waters et al., 1995
), which define discrete transitions between receptive fields when moving across the somatosensory cortex. The whisker system and the forelimb system, therefore, share a discrete architecture of somatosensory processing at the cortical level, which represents the somatotopic reference for defining the spatio-temporal structure of the receptive fields. Neurons responding to stimulation of their primary location on the body also respond to stimuli from secondary locations, much like neurons responding to stimulation of their principal whiskers respond to stimuli from surrounding whiskers (Ghazanfar and Nicolelis, 1999
, 2001
; Petersen and Diamond, 2000
). Even more specifically, cells whose primary location is on dorsal regions of the forelimb can be recruited to respond to stimulation of the glabrous skin, suggesting that these neurons are available to assist with discrimination of objects manipulated by the paw, similarly to the recruitment of cells whose principal whisker is a caudal whiskers to assist with discrimination when the rostral whiskers are stimulated (Ghazanfar and Nicolelis, 1999
, 2001
). Overall, despite the differences in the structure of ipsilateral responses, which likely have a behavioral origin (whisker movements are mostly symmetrical), our results are consistent with the idea that the same basic principles of somatosensory processing are shared by the forelimb system and the whisker system in the infragranular primary somatosensory cortex of the rat (Foffani et al., 2004
).
Implications for Somatosensory Processing
Even though our results were obtained by stimulating one location at a time, they suggest that different neurons are tuned to different spatio-temporal stimulus patterns, which could be employed by the infragranular cortex to process more complex stimuli. In fact, neurons will be likely to respond more robustly to stimuli characterized by spatio-temporal patterns (e.g. a stimulus moving across the digits) that match the spatio-temporal structure of their excitatory receptive field. Conversely, the neural responses will be less excited or even suppressed by stimuli that do not match this structure. Even though the above hypothesis requires experimental verification, supporting evidence can be found in studies investigating the responses of somatosensory neurons to stimulation of multiple whiskers (Simons, 1985
; Ghazanfar and Nicolelis, 1997
; Shimegi et al., 1999
, 2000
; Mirabella et al., 2001
; Higley and Contreras, 2003
; Ego-Stengel et al., 2005
). Within this interpretation, the great spatio-temporal heterogeneity of receptive fields across neurons becomes functionally relevant and critical to fully uncover the information sent by the cortex to other intra-cortical structures for higher-level integration (Nicolelis et al., 1998
) and to subcortical structures for somatosensory feedback (Canedo and Aguilar, 2000
; Aguilar et al., 2003
; Alitto and Usrey, 2003
; Castro-Alamancos, 2004
; Temereanca and Simons, 2004
).
Implications for Somatosensory Plasticity
The structure of the excitatory forelimb receptive fields is also important for studying cortical reorganization and plasticity after peripheral or spinal injury. Our results suggest that the cortical responses observed after digit removal (Kelahan et al., 1981
; Merzenich et al., 1983a
,b
; Byrne and Calford, 1991
) are likely due to the previously existing responses from secondary locations of the excitatory receptive field that are providing the primary neural activity recorded after the lesion. Similarly, our results explain why the amputation of an entire limb renders silent the cortical region corresponding to the distal paw (Lane et al., 1997
; Cusick, 1996
; Stojic et al., 1998
), as it is normally unlikely to respond to the remaining stump or shoulder. The size and structure of the bilateral receptive fields will likely be important for studying long-term changes after injury and be crucial for understanding the bilateral effects of lesions and the mechanisms for plastic compensatory reorganization (Kaas et al., 1983
; McKinley et al., 1987
; Florence et al., 1998
; Jones and Pons, 1998
; Dupont et al., 2001
), especially after spinal cord injury (Jain et al., 2003
; Kaas and Collins, 2003
). In conclusion, our study suggests that even small and spatially confined stimuli in the forepaw activate large areas of the infragranular somatosensory cortex with a spatio-temporally organized receptive field structure, which represents the electrophysiological architecture for studying cortical somatosensory function and dysfunction.
| Acknowledgments |
|---|
B.T. and G.F. contributed equally to this work.This work was supported by a grant from the National Institute of Health, NIH 2P50NS24707.
| References |
|---|
|
|
|---|
Aguilar J, Rivadulla C, Soto C, Canedo C (2003) New corticocuneate cellular mechanisms underlying the modulation of cutaneous ascending transmission in anesthetized cats. J Neurophysiol 89:33283339.
Alitto HJ, Usrey WM. Corticothalamic feedback and sensory processing (2003) Curr Opin Neurobiol 13:440445.[CrossRef][Web of Science][Medline]
Armstrong-James M, Fox K, Das-Gupta A (1992) Flow of excitation within rat barrel cortex on striking a single vibrissa. J Neurophysiol 68:13451358.
Brecht M, Sakmann B (2002) Dynamic representation of whisker deflection by synaptic potentials in spiny stellate and pyramidal cells in the barrels and septa of layer 4 rat somatosensory cortex. J Physiol 15:4970.
Brown PB, Koerber HR, Millecchia R (2004) From innervation density to tactile acuity 1. Spatial representation. Brain Res 1011:1432.[CrossRef][Web of Science][Medline]
Brumberg JC, Pinto DJ, Simons DJ (1996) Spatial gradients and inhibitory summation in the rat whisker barrel system. J Neurophysiol 76:13040.
Byrne JA, Calford MB (1991) Short-term expansion of receptive fields in rat primary somatosensory cortex after hindpaw digit denervation. Brain Res 565:218224.[CrossRef][Web of Science][Medline]
Canedo A, Aguilar J (2000) Spatial and cortical influences exerted on cuneothalamic and thalamocortical neurons of the cat. Eur J Neurosci 12:25152533.[CrossRef][Web of Science][Medline]
Castro-Alamancos MA (2004) Dynamics of sensory thalamocortical synaptic networks during information processing states. Prog Neurobiol 74:213247.[CrossRef][Web of Science][Medline]
Chapin JK (1986) Laminar differences in sizes, shapes, and response profiles of cutaneous receptive fields in the rat SI cortex. Exp Brain Res 62:549559.[Web of Science][Medline]
Chapin JK, Lin CS (1984) Mapping the body representation in the SI cortex of anesthetized and awake rats. J Comp Neurol 229:199213.[CrossRef][Web of Science][Medline]
Chapin JK, Sadeq M, Guise JL (1987) Corticocortical connections within the primary somatosensory cortex of the rat. J Comp Neurol 263:326346.[CrossRef][Web of Science][Medline]
Cusick CG (1996) Extensive cortical reorganization following sciatic nerve injury in adult rats versus restricted reorganization after neonatal injury: implications for spatial and temporal limits on somatosensory plasticity. Prog Brain Res 108:379390.[Web of Science][Medline]
Derdikman D, Hildesheim R, Ahissar E, Arieli A, Grinvald A (2003) Imaging spatiotemporal dynamics of surround inhibition in the barrels somatosensory cortex. J Neurosci 23:31003105.
DiCarlo JJ, Johnson KO (2002) Receptive field structure in cortical area 3b of the alert monkey. Behav Brain Res 135:167178.[CrossRef][Web of Science][Medline]
Dupont E, Canu MH, Langlet C, Falempin M (2001) Time course of recovery of the somatosensory map following hindpaw sensory deprivation in the rat. Neurosci Lett 309:121124.[CrossRef][Web of Science][Medline]
Durbin R, Mitchison G (1990) A dimensional dimension reduction framework for understanding cortical maps. Nature 343:644647.[CrossRef][Medline]
Ebner FF, Armstrong-James MA (1990) Intracortical processes regulating the integration of sensory information. Prog Brain Res 86:129141.[Medline]
Edin BB (1990) Finger joint movement sensitivity of non-cutaneous mechanoreceptor afferents in the human radial nerve. Exp Brain Res 82:417422.[Web of Science][Medline]
Edin BB (1992) Quantitative analysis of static strain sensitivity in human mechanoreceptors from hairy skin. J Neurophysiol 67:11051113.
Edin BB, Abbs JH (1991) Finger movement responses of cutaneous mechanoreceptors in the dorsal skin of the human hand. J Neurophysiol 65:657670.
Eggermont JJ, Aertsen AM, Hermes DJ, Johannesma PI (1981) Spectro-temporal characterization of auditory neurons: redundant or necessary. Hear Res 5:109121.[CrossRef][Web of Science][Medline]
Ego-Stengel V, Mello E Souza T, Jacob V, Shulz DE (2005) Spatiotemporal characteristics of neuronal sensory integration in the barrel cortex of the rat. J Neurophysiol 93:14501467.
Erzurumlu RS, Killackey HP (1980) Diencephalic projections of the subnucleus interpolaris of the brainstem trigeminal complex in the rat. Neuroscience 5:18911901.[CrossRef][Web of Science][Medline]
Essick GK, Edin BB (1995) Receptor encoding of moving tactile stimuli in humans. II. The mean response of individual low-threshold mechanoreceptors to motion across the receptive field. J Neurosci 15:848864.[Abstract]
Favorov OV, Diamond ME (1990) Demonstration of discrete place-defined columns segregates in the cat SI. J Comp Neurol 298:97112.[CrossRef][Web of Science][Medline]
Favorov O, Whitsel BL (1988a) Spatial organization of the peripheral input to area 1 cell columns. I. The detection of segregates. Brain Res 472:2542.[Medline]
Favorov O, Whitsel BL (1988b) Spatial organization of the peripheral input to area 1 cell columns. II. The forelimb representation achieved by a mosaic of segregates. Brain Res 472:4356.[Medline]
Favorov OV, Diamond ME, Whitsel BL (1987) Evidence for a mosaic representation of the body surface in area 3b of the somatic cortex of cat. Proc Natl Acad Sci USA 84:66066610.
Florence SL, Taub HB, Kaas JH (1998) Large-scale sprouting of cortical connections after peripheral injury in adult macaque monkeys. Science 282:11171121.
Foffani G, Moxon KA (2004) PSTH-based classification of sensory stimuli using ensembles of single neurons. J Neurosci Methods 135:107120.[CrossRef][Web of Science][Medline]
Foffani G, Tutunculer B, Moxon KA (2004) Role of spike timing in the forelimb somatosensory cortex of the rat. J Neurosci 24:72667271.
Fox K, Wright N, Wallace H, Glazewski S (2003) The origin of cortical surround receptive fields studied in the barrel cortex. J Neurosci 23:83808391.
Friedberg MH, Lee SM, Ebner FF (1999) Modulation of receptive field properties of thalamic somatosensory neurons by the depth of anesthesia. J Neurophysiol 81:22432252.
Ghazanfar AA, Nicolelis MA (1997) Nonlinear processing of tactile information in the thalamocortical loop. J Neurophysiol 78:506510.
Ghazanfar AA, Nicolelis MA (1999) Spatiotemporal properties of layer V neurons of the rat primary somatosensory cortex. Cereb Cortex 9:348361.
Ghazanfar AA, Nicolelis MA (2001) Feature article: the structure and function of dynamic cortical and thalamic receptive fields. Cereb Cortex 11:183193.
Grossberg S, Seitz A (2003) Laminar development of receptive fields, maps and columns in visual cortex: the coordinating role of the subplate. Cereb Cortex 13:852863.
Higley MJ, Contreras D (2003) Nonlinear integration of sensory responses in the rat barrel cortex: an intracellular study in vivo. J Neurosci 23:10190200.
Hubel DH, Wiesel TN (1962) Receptive fields, binocular interaction and functional architecture in the cat's visual cortex. J Physiol 160:106154.
Hubel DH, Wiesel TN (1968) Receptive fields and functional architecture of monkey striate cortex. J Physiol 195:215243.
Jain N, Diener PS, Coq JO, Kaas JH (2003) Patterned activity via spinal dorsal quadrant inputs is necessary for the formation of organized somatosensory maps. J Neurosci 23:1032110330.
Jones EG, Pons TP (1998) Thalamic and brainstem contributions to large-scale plasticity of primate somatosensory cortex. Science 282:11211125.
Kaas JH, Collins CE (2003) Anatomic and functional reorganization of somatosensory cortex in mature primates after peripheral nerve and spinal cord injury. Adv Neurol 93:8795.[Medline]
Kaas JH, Merzenich MM, Killackey HP (1983) The reorganization of somatosensory cortex following peripheral nerve damage in adult and developing mammals. Annu Rev Neurosci 6:325356.[CrossRef][Web of Science][Medline]
Kelahan AM, Ray RH, Carson LV, Massey CE, Doetsch GS (1981) Functional reorganization of adult raccoon somatosensory cerebral cortex following neonatal digit amputation. Brain Res 223:152159.[CrossRef][Medline]
Killackey HP, Belford GR (1979) The formation of afferent patterns in the somatosensory cortex of the neonatal rat. J Comp Neurol 183:285303.[CrossRef][Web of Science][Medline]
Killackey HP, Ebner FF (1972) Two different types of thalamocortical projections to a single cortical area in mammals. Brain Behav Evol 6:141169.[Web of Science][Medline]
Krupa DJ, Wiest MC, Shuler MG, Laubach M, Nicolelis MA (2004) Layer-specific somatosensory cortical activation during active tactile discrimination. Science 304:19891992.
Lane RD, Killackey HP, Rhoades RW (1997) Blockade of GABAergic inhibition reveals reordered cortical somatotopic maps in rats that sustained neonatal forelimb removal. J Neurophysiol 77:27232735.
Laskin SE, Spencer WA (1979) Cutaneous masking. II. Geometry of excitatory and inhibitory receptive fields of single units in somatosensory cortex of the cat. J Neurophysiol 42:10611082.
Linden JF, Schreiner CE (2003) Columnar transformations in auditory cortex? A comparison to visual and somatosensory cortices. Cereb Cortex 13:8389.
Lund JS, Angelucci A, Bressloff PC (2003) Anatomical substrates for functional columns in macaque monkey primary visual cortex. Cereb Cortex 13:1524.
Manns ID, Sakmann B, Brecht M (2004) Sub- and suprathreshold receptive field properties of pyramidal neurons in layers 5A and 5B of rat somatosensory barrel cortex. J Physiol 556:601622.
McKenna TM, Light AR, Whitsel BL (1984) Neurons with unusual response and receptive-field properties in upper laminae of cat SI cortex. J Neurophysiol 51:10551076.
McKinley PA, Jenkins WM, Smith JL, Merzenich MM (1987) Age-dependent capacity for somatosensory cortex reorganization in chronic spinal cats. Brain Res 428:136139.[Medline]
McLean J, Waterhouse BD (1994) Noradrenergic modulation of cat area 17 neuronal responses to moving visual stimuli. Brain Res 667:8397.[CrossRef][Web of Science][Medline]
Merzenich MM, Kaas JH, Wall JT, Sur M, Nelson RJ, Felleman DJ (1983a) Progression of change following median nerve section in the cortical representation of the hand in areas 3b and 1 in adult owl and squirrel monkeys. Neuroscience 10:639665.[CrossRef][Web of Science][Medline]
Merzenich MM, Kaas JH, Wall J, Nelson RJ, Sur M, Felleman D (1983b) Topographic reorganization of somatosensory cortical areas 3b and 1 in adult monkeys following restricted deafferentation. Neuroscience 8:3355.[CrossRef][Web of Science][Medline]
Mirabella G, Battiston S, Diamond ME (2001) Integration of multiple-whisker inputs in rat somatosensory cortex. Cereb Cortex 11:164170.
Montague PR (1996) The resource consumption principle: attention and memory in volumes of neural tissue. Proc Natl Acad Sci USA 93:36193623.
Moore CI, Nelson SB (1998) Spatio-temporal subthreshold receptive fields in the vibrissa representation of rat primary somatosensory cortex. J Neurophysiol 80:28822892.
Mountcastle VB (1957) Modality and topographic properties of single neurons of cat's somatic sensory cortex. J Neurophysiol 20:408434.
Mountcastle VB (1997) The columnar organization of the neocortex. Brain 120:701722.
Movshon JA, Thompson ID, Tolhurst DJ (1978) Spatial and temporal contrast sensitivity of neurones in areas 17 and 18 of the cat's visual cortex. J Physiol 283:101120.
Moxon KA, Foffani G (2003) Neurorobotics as a means for assessing brain function after spinal injury. In: Proceedings of the 1st International IEEE EMBS Conference on Neural Engineering, 2022 March 2003 (Wolf LJ, Strock JL, eds.), pp. 9598. Piscataway, NJ: IEEE.
Nicolelis MA, Chapin JK (1994) Spatiotemporal structure of somatosensory responses of many-neuron ensembles in the rat ventral posterior medial nucleus of the thalamus. J Neurosci 14:35113532.[Abstract]
Nicolelis MA, Ghazanfar AA, Stambaugh CR, Oliveira LM, Laubach M, Chapin JK, Nelson RJ, Kaas JH (1998) Simultaneous encoding of tactile information by three primate cortical areas. Nat Neurosci 1:621630.[CrossRef][Web of Science][Medline]
Olavarria J, van Sluyters RC, Killackey HP (1984) Evidence for the complementary organization of callosal and thalamic connections within rat somatosensory cortex. Brain Res 291:364368.[CrossRef][Web of Science][Medline]
Panzeri S, Petersen RS, Schultz SR, Lebedev M, Diamond ME (2001) The role of spike timing in the coding of stimulus location in rat somatosensory cortex. Neuron 29:769777.[CrossRef][Web of Science][Medline]
Panzeri S, Petroni F, Petersen RS, Diamond ME (2003) Decoding neuronal population activity in rat somatosensory cortex: role of columnar organization. Cereb Cortex 13:4552.
Petersen RS, Diamond ME (2000) Spatialtemporal distribution of whisker-evoked activity in rat somatosensory cortex and the coding of stimulus location. J Neurosci 20:61356143.
Pidoux B, Verley R (1979) Projections on the cortical somatic I barrel subfield from ipsilateral vibrissae in adult rodents. Electroencephalogr Clin Neurophysiol 46:715726.[CrossRef][Web of Science][Medline]
Plumbley MD (1999) Do cortical maps adapt to optimize information density? Network: Comput Neural Syst 10:4158.[CrossRef]
Raizada RD, Grossberg S (2003) Towards a theory of the laminar architecture of cerebral cortex: computational clues from the visual system. Cereb Cortex 13:100113.
Reid SN, Juraska JM (1991) The cytoarchitectonic boundaries of the monocular and binocular areas of the rat primary visual cortex. Brain Res 563:293296.[CrossRef][Web of Science][Medline]
Ringach DL, Hawken MJ, Shapley R (1997) Dynamics of orientation tuning in macaque primary visual cortex. Nature 387:281284.[CrossRef][Medline]
Schiller PH, Finlay BL, Volman SF (1976) Quantitative studies of single-cell properties in monkey striate cortex. I. Spatiotemporal organization of receptive fields. J Neurophysiol 39:12881319.
Shimegi S, Ichikawa T, Akasaki T, Sato H (1999) Temporal characteristics of response integration evoked by multiple whisker stimulations in the barrel cortex of rats. J Neurosci 19:1016410175.
Shimegi S, Akasaki T, Ichikawa T, Sato H (2000) Physiological and anatomical organization of multiwhisker response interactions in the barrel cortex of rats. J Neurosci 20:62416248.
Shuler MG, Krupa DJ, Nicolelis MA (2001) Bilateral integration of whisker information in the primary somatosensory cortex of rats. J Neurosci 21:52515261.
Shuler MG, Krupa DJ, Nicolelis MA (2002) Integration of bilateral whisker stimuli in rats: role of the whisker barrel cortices. Cereb Cortex 12:8697.
Simons DJ (1985) Temporal and spatial integration in the rat SI vibrissa cortex. J Neurophysiol 54:615635.
Simons DJ, Carvell GE (1989) Thalamocortical response transformation in the rat vibrissa/barrel system. J Neurophysiol 61:311330.
Sneath PHA, Sokal RR (1973) Numerical taxonomy the principles and practice of numerical classification. San Francisco, CA: W.H. Freeman.
Staiger JF, Kotter R, Zilles K, Luhmann HJ (2000) Laminar characteristics of functional connectivity in rat barrel cortex revealed by stimulation with caged-glutamate. Neurosci Res 37:4958.[CrossRef][Web of Science][Medline]
Stojic AS, Lane RD, Killackey HP, Qadri BA, Rhoades RW (1998) Thalamocortical and intracortical projections to the forelimb-stump SI representation of rats that sustained neonatal forelimb removal. J Comp Neurol 401:187204.[CrossRef][Web of Science][Medline]
Suga N, O'Neill WE, Kujirai K, Manabe T (1983) Specificity of combination-sensitive neurons for processing of complex biosonar signals in auditory cortex of the mustached bat. J Neurophysiol 49:15731626.
Temereanca S, Simons DJ (2004) Functional topography of corticothalamic feedback enhances thalamic spatial response tuning in the somatosensory whisker/barrel system. Neuron 41:639651.[CrossRef][Web of Science][Medline]
Thomson AM, Bannister AP (2003) Interlaminar connections in the neocortex. Cereb Cortex 13:514.
Van Essen DC (1997) A tension-based theory of morphogenesis and compact wiring in the central nervous system. Nature 385:313318.[CrossRef][Medline]
Wang L, Millecchia R, Brown PB (1997) Correlation of peripheral innervation density and dorsal horn map scale. J Neurophysiol 78:689702.
Waters RS, Li CX, McCandlish CA (1995) Relationship between the organization of the forepaw barrel subfield and the representation of the forepaw in layer IV of rat somatosensory cortex. Exp Brain Res 103:183197.[Web of Science][Medline]
Wheeler BC (1999) Automatic discrimination of single units. In: Methods for neural ensemble recordings (Nicolelis MAL, ed.), pp. 6177. New York: CRC Press.
Whishaw IQ, Gorny B (1994) Arpeggio and fractionated digit movements used in prehension by rats. Behav Brain Res 60:1524.[CrossRef][Web of Science][Medline]
Whishaw IQ, Pellis SM (1990) The structure of skilled forelimb reaching in the rat: a proximally driven movement with a single distal rotatory component. Behav Brain Res 41:4959.[CrossRef][Web of Science][Medline]
Whishaw IQ, O'Connor WT, Dunnett SB (1986) The contributions of motor cortex, nigrostriatal dopamine and caudate-putamen to skilled forelimb use in the rat. Brain 109:805843.
Wirth C, Luscher HR (2004) Spatiotemporal evolution of excitation and inhibition in the rat barrel cortex investigated with multielectrode arrays. J Neurophysiol 91:16351647.
Wise SP, Jones EG (1976) The organization and postnatal development of the commissural projection of the rat somatic sensory cortex. J Comp Neurol 168:313343.[CrossRef][Web of Science][Medline]
Xing J, Gerstein GL (1996) Networks with lateral connectivity. I. dynamic properties mediated by the balance of intrinsic excitation and inhibition. J Neurophysiol 75:184199.
Zarei M, Stephenson JD (1996) Ipsilateral and bilateral receptive fields in rat primary somatosensory cortex. Neuroreport 7:647651.[Medline]
Zhu JJ, Connors BW (1999) Intrinsic firing patterns and whisker-evoked synaptic responses of neurons in the rat barrel cortex. J Neurophysiol 81:11711183.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
T. Kao, J. S. Shumsky, M. Murray, and K. A. Moxon Exercise Induces Cortical Plasticity after Neonatal Spinal Cord Injury in the Rat J. Neurosci., June 10, 2009; 29(23): 7549 - 7557. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Foffani, M. L. Morales-Botello, and J. Aguilar Spike Timing, Spike Count, and Temporal Information for the Discrimination of Tactile Stimuli in the Rat Ventrobasal Complex J. Neurosci., May 6, 2009; 29(18): 5964 - 5973. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Foffani, J. K. Chapin, and K. A. Moxon Computational Role of Large Receptive Fields in the Primary Somatosensory Cortex J Neurophysiol, July 1, 2008; 100(1): 268 - 280. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. H. Fagg, N. G. Hatsopoulos, V. de Lafuente, K. A. Moxon, S. Nemati, J. M. Rebesco, R. Romo, S. A. Solla, J. Reimer, D. Tkach, et al. Biomimetic Brain Machine Interfaces for the Control of Movement J. Neurosci., October 31, 2007; 27(44): 11842 - 11846. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||













