Amplitude-dependency of response of SI cortex to flutter stimulation
© Simons et al; licensee BioMed Central Ltd. 2005
Received: 07 February 2005
Accepted: 21 June 2005
Published: 21 June 2005
It is established that increasing the amplitude of a flutter stimulus increases its perceived intensity. Although many studies have examined this phenomenon with regard to the responding afferent population, the way in which the intensity of a stimulus is coded in primary somatosensory cortex (SI) remains unclear.
Optical intrinsic signal (OIS) imaging was used to study the evoked responses in SI of anesthetized squirrel monkeys by 25 Hz sinusoidal vertical skin displacement stimulation. Stimuli were 10 sec duration with a 50 sec inter-stimulus interval. Stimulus amplitude ranged from 50 to 400 microns and different amplitudes were interleaved. Control levels of activity were measured in the absence of stimulation, and used to compare with activation levels evoked by the different stimulus amplitudes. Stimulation of a discrete skin site on the forelimb evoked a prominent increase in absorbance within the forelimb representational region in cytoarchitectonic areas 3b and 1 of the contralateral hemisphere. An increase in stimulus amplitude led to a proportional increase in the magnitude of the absorbance increase in this region of areas 3b and 1 while surrounding cortex underwent a decrease in absorbance. Correlation maps revealed that as stimulus amplitude is increased, the spatial extent of the activated region in SI remains relatively constant, and the activity within this region increases progressively. Additionally, as stimulus amplitude is increased to suprathreshold levels, activity in the surround of the activated SI territory decreases, suggesting an increase in inhibition of neuronal activity within these regions.
Increasing the amplitude of a flutter stimulus leads to a proportional increase in absorbance within the forelimb representational region of SI. This most likely reflects an increase in the firing rate of neurons in this region of SI. The relatively constant spatial extent of this stimulus-evoked increase in absorbance suggests that an increase in the amplitude of a 25 Hz skin stimulus does not evoke a larger area of SI neuronal activation due to an amplitude-dependent lateral inhibitory effect that spatially funnels the responding SI neuronal population.
The way in which stimulus intensity is represented in primary somatosensory (SI) cortex has remained an intriguing question in the study of the cortical correlates of perception. Although there have been numerous studies of the SI response to changes in stimulus intensity, few have focused on the response at the population level of analysis. Thus, a dearth of information about the global SI response to changes in stimulus intensity exists – current knowledge of the subject depends almost exclusively on reconstruction of predictions of the SI neuronal population response from afferent recordings [1–5] and single unit cortical recordings [6, 7].
A number of studies have examined the global SI response using imaging techniques such as fMRI (functional Magnetic Resonance Imaging) [8–10] and MEG (MagnetoEncephaloGraphy)[11, 12]. In general, results from these studies indicate that increases in stimulus intensity are accompanied by increases in the intensity of the evoked signal as well as increases in the activated volume of cortex. As a result these studies predicted that amplitude might be coded not only by the average firing rates of individual SI neurons, but also by the total aggregate of responding neurons. However, although each of these imaging techniques provide measures which are indirectly related to neuronal activity, their resolution is limited in two important ways. First, it is difficult to determine the nature of the neuronal activity being imaged (whether or not it is excitatory or inhibitory); and second, both fMRI and MEG studies have limited spatial resolution, which is typically on the order of ~1 mm2 [8, 11].
Recently, Chen and colleagues used the optical intrinsic signal (OIS) to demonstrate that a proportionally greater (larger magnitude) response is evoked in SI of squirrel monkeys as the amplitude (as measured by force) of a skin stimulus is increased. However, the primary focus of their report was that the response to simultaneous stimulation of multiple adjacent sites on the skin produced a smaller, more intense region of SI activation than would be normally predicted by summation of the two, and their findings did not detail the effects of amplitude on the dynamics of the response to a flutter stimulus at a single site on the skin. In this report, we extend the aforementioned work by imaging the OIS evoked in SI cortex of squirrel monkeys by a range of amplitudes of skin flutter stimulation. The results suggest that increasing the amplitude of a skin flutter stimulus evokes a proportionally larger absorbance increase in SI that remains confined to the same SI territory. In addition, it was found that increasing the amplitude of flutter evokes a large decrease in absorbance in the territory that borders the activated region of SI. Neurons in the SI region that demonstrate decreased absorbance in response to flutter stimulation are proposed to undergo stimulus-evoked inhibition and to contribute importantly to the SI processing of high-amplitude skin flutter stimuli.
This study evaluated the global SI response to different amplitudes of flutter stimulation by imaging the optical intrinsic signal (OIS). The OIS indirectly reflects both cortical neuronal spike discharge activity and the local, subthreshold changes in neuronal membrane potential evoked by sensory stimulation [14–16]. As a result, the observed tendency for absorbance in the same localized region of area 3b to increase with increasing stimulus amplitude (Figures 3 &4) most likely is due to the amplitude-dependence of the average firing rate of neurons in the same region . An increase in SI absorbance in response to an increase in stimulus intensity has been reported previously by Chen et al . The observed increase we observed in the stimulus-evoked SI OIS that accompanies increases in stimulus amplitude is well described by a linear function.
One important distinction between previous work done using the OIS and the present study is our use of near-infrared wavelengths for illumination during acquisition of the OIS. The OIS obtained using infrared light has been shown to be highly correlated with light scattering effects that accompany astrocyte swelling subsequent to the clearance of extracellular K+ and neurotransmitter [14, 15, 18] and local increases in blood volume [19, 20]. Although the OIS at this wavelength may be influenced by changes in hemoglobin concentration and oxygenation, it is likely that contributions from these factors are small when compared with light scattering effects. Additionally, OIS imaging using near-infrared illumination not only minimizes the contributions of artifacts introduced by changes in the vasculature (which can dominate the OIS at lower wavelengths)  but the time course of the OIS detected at shorter wavelengths (600 nm) is markedly different (shorter) than the protracted OIS observed in this study [13, 15, 19].
Previous studies conducted by this laboratory have reported that the SI optical response evoked by an extended period (>1 sec) of flutter stimulation not only consists of an increase in absorbance in the region that receives its input from the skin site that was stimulated, but also decreases in absorbance (frequently to levels well below-background) that occur in the surrounding cortex . The present study demonstrates that the below-background component of the SI optical response to flutter stimulation is particularly evident at large stimulus amplitudes (figures 5, 7 &8). However, unlike the increase in absorbance evoked by flutter, the relationship between the magnitude of the stimulus-evoked decrease in absorbance and stimulus amplitude is not satisfactorily described by a linear function. Indeed, the results shown in figure 6 suggest that this component of the optical response to skin flutter is either absent or extremely small at small stimulus amplitudes, and remains maximal or near-maximal across a wide range of suprathreshold stimulus amplitudes. Interestingly, this stimulus-evoked decrease in absorbance (in the surround) appears later in time than the increases in absorbance (within the ROI) and as amplitude is increased it tends to develop earlier after stimulus onset.
The correlation maps shown in figure 7 provide a comprehensive overview of the time course of absorbance at each location in the image. The optical signal at each pixel is cross-correlated with a known input signal. In this case, it is the average absorbance measured within the ROI. The assigned coefficient of determination indicates the degree of similarity between a pixel's time course of absorbance and the input signal. Therefore pixels with a large positive correlation undergo increases in absorbance very similar to the input signal, while pixels with a large negative correlation undergo a decrease in absorbance which more closely resembles the opposite (negative) of the input signal. Figure 8 suggests that at high amplitudes of stimulation the ROI in SI becomes more homogeneously activated with longer stimulus duration. Some evidence for this is indicated by the large discrepancy (between low- and high-amplitude surface plots) in the standard deviations measured within the ROI. Further studies are required to investigate absorbance distribution and patterning within the ROI.
Examination of spatial histograms (figure 2) and the maps in figures 5 &7 also reveal that the size of the SI region that undergoes an increase in absorbance does not increase with increasing stimulus amplitude, but rather remains relatively constant. Regardless of stimulus amplitude, the activated cortical region appears circular in shape and occupies an area approximately 2 mm in diameter (figures 4 &6). Within the ROI average absorbance increases progressively with increasing stimulus duration. The dimensions of the SI region activated by flutter stimulation observed in this study contrast sharply with results presented in previous studies which demonstrated activation within a 1 mm region in diameter [13, 22]. One possible explanation for this discrepancy is the level and type of general anesthesia used in the different studies. Previous studies have reported that anesthetics (e.g. ketamine) which block NMDA receptors or enhance GABAA receptor mediated inhibition (barbiturates), significantly reduce the dimensions of the receptive field of individual SI neurons; actions that would reduce the size of the responding SI neuronal population . Chen et al. previously reported similar (~2 mm) sized regions of activation in response to flutter stimulation of the digit tips in squirrel monkey anesthetized with isofluorane, as well as showing in the same report that use of pentothal anesthetic confined the response to a much smaller region (~1 mm) .
It has been suggested that the amplitude of skin flutter stimulation is coded by both the number of activated SI neurons as well as by their level of spike discharge activity . This suggestion is based largely in part on the fact that larger-amplitude stimuli, through transduction of the laterally-transmitted mechanical wave produced by sinusoidal skin displacement, recruit larger numbers of RA afferents and therefore lead to a spatially more widely distributed pattern of afferent input to SI cortex. Combined metabolic tracer and neurophysiological studies have shown that the initial response to a repetitive tactile stimulus occupies an extremely large cortical territory. As the repetitive mechanical stimulation is continued, however, the response is quickly sculpted by cortical inhibitory mechanisms, leading to an activity pattern that becomes confined to a relatively restricted region in SI [25–27]. The results obtained in the present study and results previously reported by other researchers, lead us to suggest that stimulus amplitude contributes importantly to the shaping (via lateral inhibitory mechanisms) of the SI response to protracted skin flutter.
This study investigated the SI response to flutter stimulation of the skin using the OIS. An increase of the amplitude of the flutter stimulus was associated with an increase in absorbance within the region of SI cortex that receives its input from the stimulated skin field. The relationship between the maximal change in absorbance and stimulus amplitude is well characterized by a linear function within the range of amplitudes studied. Measurement of the spatial extent of the activated SI region showed that higher amplitudes of stimulation do not produce a more extensive region of SI activation. Instead, as amplitude is increased, while average peak absorbance within the same ~2 mm diameter SI region increases with amplitude of stimulation, the region of surrounding cortex undergoes a prominent decrease (frequently to levels well below background) in absorbance. Further studies are required to establish the relationship between the effect of different amplitudes of skin flutter stimulation on SI absorbance and SI neuroelectrical activity.
Subjects & preparation
All methods and procedures are consistent with USPHS policies and guidelines on animal care and welfare in biomedical research. They were reviewed and approved by an institutional committee prior to initiation of the experiments. Experiments were conducted in 10 squirrel monkeys. Following induction of anesthesia with 4% halothane in a 50/50 mix of nitrous oxide (N2O) and oxygen, the trachea was intubated. A veterinary anesthesia machine (Forreger Compac-75) provided an anesthetic gas mix whose composition could be adjusted (typically 1.5–3.0% halothane in 50/50 N2O/oxygen) to maintain a stable level of surgical anesthesia. Methylprednisolone sodium succinate (20 mg/kg) and gentamicin sulfate (2.5 mg/kg) were injected intramuscularly to lessen the probability of halothane-induced cerebral edema and prevent bacterial septicemia, respectively. Placement of a valved catheter into a superficial hindlimb vein enabled administration of 5%glucose, 0.9%NaCl, and drugs.
A 1.5 cm opening was trephined in the skull overlying SI cortex. A recording chamber (25 mm i.d.) was placed over the opening and cemented to the skull with dental acrylic. Wound margins were infiltrated with local anesthetic, closed with sutures and bandaged, and the dura overlying SI was resected. After the completion of all surgical procedures subjects were immobilized with norcuron (loading dose 0.25–0.5 mg/kg,i.v.; maintenance dose 0.025–0.05 mg/kg/hr). From this point on, the animal was ventilated with a 50/50 mix of N2O and oxygen and the concentration of halothane was adjusted (typically between 0.5 and 1.0%) to maintain heart rate, blood pressure, and the EEG at values consistent with general anesthesia. Rate and depth of ventilation were adjusted to maintain end-tidal CO2 between 3.0 and 4.5%. Under these experimental/anesthetic conditions both SI neuron spontaneous and stimulus-evoked spike discharge activity patterns are highly reproducible over even prolonged (>1 hr) time periods.
OIS imaging and stimulus protocol
After obtaining a photograph of the exposed cortical surface, the recording chamber was filled with artificial cerebrospinal fluid, and hydraulically sealed using a clear glass plate. In each of five experiments, the OIS evoked in SI by cutaneous flutter stimuli on the thenar region of the forelimb was recorded. The flutter stimulus was delivered at 50, 100, 200 and 400 μm and was interleaved in order to prevent conditioning of the response. The imaging system consisted of a computer-interfaced CCD camera (Quantix 540 from Roper Scentific), light source, guide and filters required for near-infrared (833 nm) illumination of the cortical surface, a focusing device, and a recording chamber capped by an optical window (for additional methodological details see Tommerdahl et al, 1999a,b). Images of the exposed anterior parietal and surrounding cortical surface were acquired 200 ms before stimulus onset ("reference images") and continuously thereafter for 22 s after stimulus onset ("post-stimulus images") at a rate of one image every 0.9–1.4 s. Exposure time was 200 ms. Difference images were generated by subtracting each pre-stimulus image from its corresponding post-stimulus image. Averaged difference images typically show regions of both increased light absorption (decreased reflectance) and decreased light absorption (increased reflectance) which are believed widely (e.g., Grinvald 1985; Grinvald et al. 1991a,b) to be accompanied by increases and decreases in neuronal activation, respectively.
Although onset of the OIS is delayed at longer wavelengths [19, 20], use of near-infrared illumination minimizes the contributions to OIS images of the changes in blood flow and flow/volume that normally accompany cortical neuronal activation, and thus may be more useful in the description of spatial characteristics of the signal . All images were examined prior to their inclusion for analysis. Images containing random high amplitude noise were excluded, and the remaining trials (typically ~20) were averaged to improve the signal to noise ratio. OIS images were analyzed using custom routines written in Matlab.
Correlation maps were constructed for comparison of spatial characteristics of the OIS response. This method of analysis has been previously described in detail . Briefly, maps were constructed by choosing a reference region within the imaged field and computing the intensity correlation rij between the reflectance value of each pixel (i, j) and the average reflectance value within the reference region over the time from stimulus onset to stimulus offset. The region selected as the reference was defined by a boxel (π mm2 area) centered on the region of interest (ROI). Each pixel (i, j) on the correlation map is represented by a correlation coefficient r ij (-1 <r < 1; - 1 indicates negative correlation; + 1 indicates positive. The statistical significance of each of the correlations was tested with the standard t-test.
Histological procedures/identification of cytoarchitectural boundaries
At the conclusion of the experiment, the imaged cortical region was removed immediately following intracardial perfusion with saline and fixative. The region then was blocked, postfixed, cryoprotected, frozen, sectioned serially at 30 μm, and the sections stained with cresyl fast violet. The boundaries between adjacent cytoarchitectonic areas were identified by scanning individual sagittal sections separated by no more than 300 mm and were plotted at high resolution using a microscope with a drawing tube attachment. The resulting plots then were used to reconstruct a two-dimensional surface map of the cytoarchitectonic boundaries within the region studied with optical and neurophysiological recording methods. The locations of microelectrode tracks and electrolytic lesions evident in the histological sections were projected radially to the pial surface and transferred to the map of cytoarchitectonic boundaries reconstructed from the same sections. As the final step, the cytoarchitectonic boundaries (along with the locations of microelectrode tracks and lesions whenever present) identified in each brain were mapped onto the images of the stimulus-evoked intrinsic signal obtained from the same subject, using fiducial points (made by postmortem applications of india ink or needle stabs) as well as morphological landmarks (e.g., blood vessels and sulci evident both in the optical images and in histological sections). Locations of cytoarchitectonic boundaries were identified using established criteria [29–31].
This work was supported, in part, by US Army Research Office grant P43077-LS (M. Tommerdahl, P.I.), NIH NS050587 (M. Tommerdahl, P.I.) and NIH NS35222 (B. Whitsel, P.I.).
- Johnson K: Reconstruction of population response to a vibratory stimulus in quickly adapting mechanoreceptive afferent fiber population innervating glabrous skin of the monkey. J Neurophysiol. 1974, 37 (1): 48-72.PubMedGoogle Scholar
- Connor C, Hsiao S, Phillips J, Johnson K: Tactile roughness: neural codes that account for psychophysical magnitude estimates. J Neurosci. 1990, 10 (12): 3823-3836.PubMedGoogle Scholar
- Guclu B, Bolanowski S: Modelling population responses of rapidly-adapting mechanoreceptive fibers. J Comput Neurosci. 2002, 12 (3): 201-218. 10.1023/A:1016535413000.View ArticlePubMedGoogle Scholar
- Whitsel B, Kelly E, Delemos K, Xu M, Quibrera P: Stability of rapidly adapting afferent entrainment vs responsivity. Somatosens Mot Res. 2000, 17 (1): 13-31. 10.1080/08990220070265.View ArticlePubMedGoogle Scholar
- Werner G, Mountcastle V: Neural activity in mechanoreceptive cutaneous afferents: stimulus-response relations, weber functions, and information transmission. J Neurophysiol. 1965, 28: 359-397.PubMedGoogle Scholar
- Mountcastle V, Poggio G, Werner G: The relation of thalamic cell response to peripheral stimuli varied over an intensive continuum. J Neurophysiol. 1963, 26: 807-834.PubMedGoogle Scholar
- Simons D: Response properties of vibrissa units in rat SI somatosensory neocortex. J Neurophysiol. 1978, 41 (3): 798-820.PubMedGoogle Scholar
- Nelson A, Staines W, Graham S, McIlroy W: Activation in SI and SII; the influence of vibrotactile amplitude during passive and task-relevant stimulation. Cogn Brain Res. 2004, 19: 174-184. 10.1016/j.cogbrainres.2003.11.013.View ArticleGoogle Scholar
- Arthurs O, Williams E, Carpenter T, Pickard J, Boniface S: Linear coupling between functional magnetic resonance imaging and evoked potential amplitude in human somatosensory cortex. Neuroscience. 2000, 101 (4): 803-806. 10.1016/S0306-4522(00)00511-X.View ArticlePubMedGoogle Scholar
- Backes W, Mess W, van Kranen-Mastenbroek V, Reulen J: Somatosensory cortex responses to median nerve stimulation: fMRI effects of current amplitude and selective attention. Clin Neurophysiol. 2000, 111: 1738-1744. 10.1016/S1388-2457(00)00420-X.View ArticlePubMedGoogle Scholar
- Torquati K, Pizzella V, Penna S, Franciotti R, Babiloni C, Rossini P, Romani G: Comparison between SI and SII responses as a function of stimulus intensity. NeuroReport. 2002, 13: 813-819. 10.1097/00001756-200205070-00016.View ArticlePubMedGoogle Scholar
- Iguchi Y, Hoshi Y, Tanosaki M, Taira M, Hashimoto I: Selective attention regulates spatial and intensity information processing in the human primary somatosensory cortex. NeuroReport. 2002, 13: 2335-2339. 10.1097/00001756-200212030-00033.View ArticlePubMedGoogle Scholar
- Chen LM, Friedman RM, Roe Aw: Optical imaging of a tactile illusion in area 3b of the primary somatosensory cortex. Science. 2003, 302: 881-885. 10.1126/science.1087846.View ArticlePubMedGoogle Scholar
- Holthoff K, Witte O: Intrinsic optical signals in rat neocortical slices measured with near-infrared dark-field microscopy reveal changes in extracellular space. J Neurosci. 16 (8): 2740-2749. 1996, Apr 15Google Scholar
- Kohn A, Metz C, Quibrera M, Tommerdahl M, Whitsel B: Functional neocortical microcircuitry demonstrated with intrinsic signal optical imaging in vitro. Neuroscience. 2000, 95 (1): 51-62. 10.1016/S0306-4522(99)00385-1.View ArticlePubMedGoogle Scholar
- MacVicar BA, Hochman D: Imaging of synaptically evoked intrinsic optical signals in hippocampal slices. J Neurosci. 1991, 11: 1458-1469.PubMedGoogle Scholar
- Whitsel BL, Kelly EF, Quibrera M, Tommerdahl M, Li Y, Favorov OV, Xu M, Metz CB: Time-dependence of SI RA neuron response to cutaneous flutter stimulation. Somatosens Mot Res. 2003, 20 (1): 45-69. 10.1080/0899022031000083834.View ArticlePubMedGoogle Scholar
- Lee J, Tommerdahl M, Favorov OV, Whitsel BL: The optically recorded response of the superficial dorsal horn: dissociation from neuronal activity, sensitivity to formalin-evoked skin nociceptor activation. J Neurophysiol. 2005, Mar 2Google Scholar
- Frostig RD, Lieke EE, Ts'o DY, Grinvald A: Cortical functional architecture and local coupling between neuronal activity and the microcirculation revealed by in vivo high-resolution optical imaging of intrinsic signals. Proc Natl Acad Sci USA. 1990, 87 (16): 6082-6086.PubMed CentralView ArticlePubMedGoogle Scholar
- Ba A, Guiou M, Pouratian N, Muthialu A, Rex D, Cannestra A, Chen J, Toga AW: Multiwavelength optical intrinsic signal imaging of cortical spreading depression. J Neurophysiol. 2002, 88 (5): 2726-2735.View ArticlePubMedGoogle Scholar
- Tommerdahl M, Favorov OV, Whitsel BL: Optical imaging of intrinsic signals in somatosensory cortex. Behav Brain Res. 2002, 135: 83-91. 10.1016/S0166-4328(02)00159-6.View ArticlePubMedGoogle Scholar
- Shoham D, Grinvald A: The cortical representation of the hand in macaque and human area S-I: high resolution optical imaging. J Neurosci. 2001, 21 (17): 6820-6825.PubMedGoogle Scholar
- McKenna T, Whitsel BL, Dreyer D: Anterior parietal cortical topographic organization in macaque monkey: a reevaluation. J Neurophysiol. 1982, 48 (2): 289-317.PubMedGoogle Scholar
- Chen LM, Friedman RM, Ramsden BM, LaMotte RH, Roe AW: Fine-scale organization of SI (Area 3b) in the squirrel monkey revealed with intrinsic optical imaging. J Neurophysiol. 2001, 86: 3011-3029.PubMedGoogle Scholar
- Juliano SL, Whitsel BL: A combined 2-deoxyglucose and neurophysiological study of primate somatosensory cortex. J Comp Neurol. 1987, 263: 514-525. 10.1002/cne.902630405.View ArticlePubMedGoogle Scholar
- Whitsel BL, Favorov OV, Tommerdahl M, Diamond M, Juliano S, Kelly D: Dynamic processes govern the somatosensory cortical response to natural stimulation. Sensory Processing in the Mammalian Brain. Edited by: Lund JS. 1989, New York: Oxford Univ. Press, 79-107.Google Scholar
- Tommerdahl M, Whitsel B: Optical imaging of intrinsic signals in somatosensory cortex. Somesthesis and the Neurobiology of Somatosensory Cortex. Edited by: Franzen O, Johansson R, Terenius L. 1996, Basel: Birkhauser Verlag AB, 369-384.View ArticleGoogle Scholar
- Tommerdahl M, Whitsel B, Favorov O, Metz C, BL OQ: Responses of contralateral SI and SII in cat to same site cutaneous flutter versus vibration. J Neurophysiol. 1999, 82: 1982-1992.PubMedGoogle Scholar
- Jones E, Porter R: What is area 3a?. Brain Res. 1980, 203 (1): 1-43.View ArticlePubMedGoogle Scholar
- Powell T, Mountcastle V: The cytoarchitecture of the postcentral gyrus of the monkey Macaca mulatta. Bull Johns Hopkins Hosp. 1959, 105: 108-131.PubMedGoogle Scholar
- Sur M, Nelson R, Kaas J: Representations of the body surface in cortical areas 3b and 1 of squirrel monkeys: comparisons with other primates. J Comp Neurol. 211 (2): 177-192. 10.1002/cne.902110207. 1982, Oct 20Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.