Fine-scale organization of SI (area 3b) in the squirrel monkey revealed with intrinsic optical imaging

Motor actions are spatially organized in motor and dorsal premotor cortex

Frontal motor areas are central to controlling voluntary movements. In non-human primates, the motor areas contain independent, somatotopic, representations of the forelimb (i.e., motor maps). But are the neural codes for actions spatially organized within those forelimb representations? Addressing this question would provide insight into the poorly understood structure-function relationships of the cortical motor system. Here, we tackle the problem using high-resolution optical imaging and motor mapping in motor (M1) and dorsal premotor (PMd) cortex. Two macaque monkeys performed an instructed reach-to-grasp task while cortical activity was recorded with intrinsic signal optical imaging (ISOI). The spatial extent of activity in M1 and PMd was then quantified in relation to the forelimb motor maps, which we obtained from the same hemisphere with intracortical microstimulation. ISOI showed that task-related activity was concentrated in patches that collectively overlapped <40% of the M1 and PMd forelimb representations. The spatial organization of the patches was consistent across task conditions despite small variations in forelimb use. Nevertheless, the largest condition differences in forelimb use were reflected in the magnitude of cortical activity. Distinct time course profiles from patches in arm zones and patches in hand zones suggest functional differences within the forelimb representations. The results collectively support an organizational framework wherein the forelimb representations contain subzones enriched with neurons tuned for specific actions. Thus, the often-overlooked spatial dimension of neural activity appears to be an important organizing feature of the neural code in frontal motor areas.

Similar somatotopy for active and passive digit representation in primary somatosensory cortex

Scientists traditionally use passive stimulation to examine the organisation of primary somatosensory cortex (SI). However, given the close, bidirectional relationship between the somatosensory and motor systems, active paradigms involving free movement may uncover alternative SI representational motifs. Here, we used 7 Tesla functional magnetic resonance imaging to compare hallmark features of SI digit representation between active and passive tasks which were unmatched on task or stimulus properties. The spatial location of digit maps, somatotopic organisation, and inter-digit representational structure were largely consistent between tasks, indicating representational consistency. We also observed some task differences. The active task produced higher univariate activity and multivariate representational information content (inter-digit distances). The passive task showed a trend towards greater selectivity for digits versus their neighbours. Our findings highlight that, while the gross features of SI functional organisation are task invariant, it is important to also consider motor contributions to digit representation.

Perception of microstimulation frequency in human somatosensory cortex

Microstimulation in the somatosensory cortex can evoke artificial tactile percepts and can be incorporated into bidirectional brain–computer interfaces (BCIs) to restore function after injury or disease. However, little is known about how stimulation parameters themselves affect perception. Here, we stimulated through microelectrode arrays implanted in the somatosensory cortex of two human participants with cervical spinal cord injury and varied the stimulus amplitude, frequency, and train duration. Increasing the amplitude and train duration increased the perceived intensity on all tested electrodes. Surprisingly, we found that increasing the frequency evoked more intense percepts on some electrodes but evoked less-intense percepts on other electrodes. These different frequency–intensity relationships were divided into three groups, which also evoked distinct percept qualities at different stimulus frequencies. Neighboring electrode sites were more likely to belong to the same group. These results support the idea that stimulation frequency directly controls tactile perception and that these different percepts may be related to the organization of somatosensory cortex, which will facilitate principled development of stimulation strategies for bidirectional BCIs.

A touch of hierarchy: population receptive fields reveal fingertip integration in Brodmann areas in human primary somatosensory cortex

Several neuroimaging studies have shown the somatotopy of body part representations in primary somatosensory cortex (S1), but the functional hierarchy of distinct subregions in human S1 has not been adequately addressed. The current study investigates the functional hierarchy of cyto-architectonically distinct regions, Brodmann areas BA3, BA1, and BA2, in human S1. During functional MRI experiments, we presented participants with vibrotactile stimulation of the fingertips at three different vibration frequencies. Using population Receptive Field (pRF) modeling of the fMRI BOLD activity, we identified the hand region in S1 and the somatotopy of the fingertips. For each voxel, the pRF center indicates the finger that most effectively drives the BOLD signal, and the pRF size measures the spatial somatic pooling of fingertips. We find a systematic relationship of pRF sizes from lower-order areas to higher-order areas. Specifically, we found that pRF sizes are smallest in BA3, increase slightly towards BA1, and are largest in BA2, paralleling the increase in visual receptive field size as one ascends the visual hierarchy. Additionally, we find that the time-to-peak of the hemodynamic response in BA3 is roughly 0.5 s earlier compared to BA1 and BA2, further supporting the notion of a functional hierarchy of subregions in S1. These results were obtained during stimulation of different mechanoreceptors, suggesting that different afferent fibers leading up to S1 feed into the same cortical hierarchy.

High-density mapping of primate digit representations with a 1152-channelµECoG array

Objective.Advances in brain-machine interfaces (BMIs) are expected to support patients with movement disorders. Electrocorticogram (ECoG) measures electrophysiological activities over a large area using a low-invasive flexible sheet placed on the cortex. ECoG has been considered as a feasible signal source of the clinical BMI device. To capture neural activities more precisely, the feasibility of higher-density arrays has been investigated. However, currently, the number of electrodes is limited to approximately 300 due to wiring difficulties, device size, and system costs.Approach.We developed a high-density recording system with a large coverage (14 × 7 mm²) and using 1152 electrodes by directly integrating dedicated flexible arrays with the neural-recording application-specific integrated circuits and their interposers.Main results.Comparative experiments with a 128-channel array demonstrated that the proposed device could delineate the entire digit representation of a nonhuman primate. Subsampling analysis revealed that higher-amplitude signals can be measured using higher-density arrays.Significance.We expect that the proposed system that simultaneously establishes large-scale sampling, high temporal-precision of electrophysiology, and high spatial resolution comparable to optical imaging will be suitable for next-generation brain-sensing technology.

Optical imaging reveals functional domains in primate sensorimotor cortex

Motor cortex (M1) and somatosensory cortex (S1) are central to arm and hand control. Efforts to understand encoding in M1 and S1 have focused on temporal relationships between neural activity and movement features. However, it remains unclear how the neural activity is spatially organized within M1 and S1. Optical imaging methods are well-suited for revealing the spatio-temporal organization of cortical activity, but their application is sparse in monkey sensorimotor cortex. Here, we investigate the effectiveness of intrinsic signal optical imaging (ISOI) for measuring cortical activity that supports arm and hand control in a macaque monkey. ISOI revealed spatial domains that were active in M1 and S1 in response to instructed reaching and grasping. The lateral M1 domains overlapped the hand representation and contained a population of neurons with peak firing during grasping. In contrast, the medial M1 domain overlapped the arm representation and a population of neurons with peak firing during reaching. The S1 domain overlapped the hand representations of areas 1 and 2 and a population of neurons with peak firing upon hand contact with the target. Our single unit recordings indicate that ISOI domains report the locations of spatial clusters of functionally related neurons. ISOI is therefore an effective tool for surveilling the neocortex for “hot zones” of activity that supports movement. Combining the strengths of ISOI with other imaging modalities (e.g., fMRI, 2-photon) and with electrophysiological methods can open new frontiers in understanding the spatio-temporal organization of cortical signals involved in movement control.

Mapping mesoscale cortical connectivity in monkey sensorimotor cortex with optical imaging and microstimulation

To map in vivo cortical circuitry at the mesoscale, we applied a novel approach to map interareal functional connectivity. Electrical intracortical microstimulation (ICMS) in conjunction with optical imaging of intrinsic signals (OIS) was used map functional connections in somatosensory cortical areas in anesthetized squirrel monkeys. ICMS produced activations that were focal and that displayed responses which were stimulation intensity dependent. ICMS in supragranular layers of Brodmann Areas 3b, 1, 2, 3a, and M1 evoked interareal activation patterns that were topographically appropriate and appeared consistent with known anatomical connectivity. Specifically, ICMS revealed Area 3b connections with Area 1; Area 1 connections with Areas 2 and 3a; Area 2 connections with Areas 1, 3a, and M1; Area 3a connections with Areas M1, 1, and 2; and M1 connections with Areas 3a, 1, and 2. These somatosensory connectivity patterns were reminiscent of feedforward patterns observed anatomically, although feedback contributions are also likely present. Further consistent with anatomical connectivity, intra-areal and intra-areal patterns of activation were patchy with patch sizes of 200-300 μm. In summary, ICMS with OIS is a novel approach for mapping interareal and intra-areal connections in vivo. Comparisons with feedforward and feedback anatomical connectivity are discussed.

Focal infrared neural stimulation with high-field functional MRI: A rapid way to map mesoscale brain connectomes

We have developed a way to map brain-wide networks using focal pulsed infrared neural stimulation in ultrahigh-field magnetic resonance imaging (MRI). The patterns of connections revealed are similar to those of connections previously mapped with anatomical tract tracing methods. These include connections between cortex and subcortical locations and long-range cortico-cortical connections. Studies of local cortical connections reveal columnar-sized laminar activation, consistent with feed-forward and feedback projection signatures. This method is broadly applicable and can be applied to multiple areas of the brain in different species and across different MRI platforms. Systematic point-by-point application of this method may lead to fundamental advances in our understanding of brain connectomes.

Neural Basis of Touch and Proprioception in Primate Cortex

The sense of proprioception allows us to keep track of our limb posture and movements and the sense of touch provides us with information about objects with which we come into contact. In both senses, mechanoreceptors convert the deformation of tissues-skin, muscles, tendons, ligaments, or joints-into neural signals. Tactile and proprioceptive signals are then relayed by the peripheral nerves to the central nervous system, where they are processed to give rise to percepts of objects and of the state of our body. In this review, we first examine briefly the receptors that mediate touch and proprioception, their associated nerve fibers, and pathways they follow to the cerebral cortex. We then provide an overview of the different cortical areas that process tactile and proprioceptive information. Next, we discuss how various features of objects-their shape, motion, and texture, for example-are encoded in the various cortical fields, and the susceptibility of these neural codes to attention and other forms of higher-order modulation. Finally, we summarize recent efforts to restore the senses of touch and proprioception by electrically stimulating somatosensory cortex. © 2018 American Physiological Society. Compr Physiol 8:1575-1602, 2018.

Connectivity of neuronal populations within and between areas of primate somatosensory cortex

Functions of the cerebral cortex emerge via interactions of horizontally distributed neuronal populations within and across areas. However, the connectional underpinning of these interactions is not well understood. The present study explores the circuitry of column-size cortical domains within the hierarchically organized somatosensory cortical areas 3b and 1 using tract tracing and optical intrinsic signal imaging (OIS). The anatomical findings reveal that feedforward connections exhibit high topographic specificity, while intrinsic and feedback connections have a more widespread distribution. Both intrinsic and inter-areal connections are topographically oriented across the finger representations. Compared to area 3b, the low clustering of connections and small cortical magnification factor supports that the circuitry of area 1 scaffolds a sparse functional representation that integrates peripheral information from a large area that is fed back to area 3b. Fast information exchange between areas is ensured by thick axons forming a topographically organized, reciprocal pathway. Moreover, the highest density of projecting neurons and groups of axon arborization patches corresponds well with the size and locations of the functional population response reported by OIS. The findings establish connectional motifs at the mesoscopic level that underpin the functional organization of the cerebral cortex.

Discrete Modules and Mesoscale Functional Circuits for Thermal Nociception within Primate S1 Cortex

This study addresses one long-standing question of whether functional separations are preserved for somatosensory modalities of touch, heat, and cold nociception within primate primary somatosensory (S1) cortex. This information is critical for understanding how the nature of pain is represented in the primate brain. Using a combination of submillimeter-resolution fMRI and microelectrode local field potential (LFP) and spike recordings, we identified spatially segregated cortical zones for processing touch and nociceptive heat and cold stimuli in somatotopically appropriate areas 3a, 3b, 1, and 2 of S1 in male monkeys. The distances between zones were comparable (∼3.4 mm) across stimulus modalities (heat, cold, and tactile), indicating the existence of uniform, modality-specific modules. Stimulus-evoked LFP maps validated the fMRI maps in areas 3b and 1. Isolation of heat and cold nociceptive neurons from the fMRI zones confirmed the validity of using fMRI to probe nociceptive regions and circuits. Resting-state fMRI analysis revealed distinct intrinsic functional circuits among functionally related zones. We discovered distinct modular structures and networks for thermal nociception within S1 cortex, a finding that has significant implications for studying chronic pain syndromes and guiding the selection of neuromodulation targets for chronic pain management.SIGNIFICANCE STATEMENT Primate S1 subregions contain discrete heat and cold nociceptive modules. Modules with the same properties exhibit strong functional connection. Nociceptive fMRI response coincides with LFP and spike activities of nociceptive neurons. Functional separation of heat and cold pain is retained within primate S1 cortex.

Somatosensory maps

Somatosensory areas containing topographic maps of the body surface are a major feature of parietal cortex. In primates, parietal cortex contains four somatosensory areas, each with its own map, with the primary cutaneous map in area 3b. Rodents have at least three parietal somatosensory areas. Maps are not isomorphic to the body surface, but magnify behaviorally important skin regions, which include the hands and face in primates, and the whiskers in rodents. Within each map, intracortical circuits process tactile information, mediate spatial integration, and support active sensation. Maps may also contain fine-scale representations of touch submodalities, or direction of tactile motion. Functional representations are more overlapping than suggested by textbook depictions of map topography. The whisker map in rodent somatosensory cortex is a canonic system for studying cortical microcircuits, sensory coding, and map plasticity. Somatosensory maps are plastic throughout life in response to altered use or injury. This chapter reviews basic principles and recent findings in primate, human, and rodent somatosensory maps.

Correlated Disruption of Resting-State fMRI, LFP, and Spike Connectivity between Area 3b and S2 following Spinal Cord Injury in Monkeys

This study aims to understand how functional connectivity (FC) between areas 3b and S2 alters following input deprivation and the neuronal basis of disrupted FC of resting-state fMRI signals. We combined submillimeter fMRI with microelectrode recordings to localize the deafferented digit regions in areas 3b and S2 by mapping tactile stimulus-evoked fMRI activations before and after cervical dorsal column lesion in each male monkey. An average afferent disruption of 97% significantly reduced fMRI, local field potential (LFP), and spike responses to stimuli in both areas. Analysis of resting-state fMRI signal correlation, LFP coherence, and spike cross-correlation revealed significantly reduced functional connectivity between deafferented areas 3b and S2. The degrees of reductions in stimulus responsiveness and FC after deafferentation differed across fMRI, LFP, and spiking signals. The reduction of FC was much weaker than that of stimulus-evoked responses. Whereas the largest stimulus-evoked signal drop (∼80%) was observed in LFP signals, the greatest FC reduction was detected in the spiking activity (∼30%). fMRI signals showed mild reductions in stimulus responsiveness (∼25%) and FC (∼20%). The overall deafferentation-induced changes were quite similar in areas 3b and S2 across signals. Here we demonstrated that FC strength between areas 3b and S2 was much weakened by dorsal column lesion, and stimulus response reduction and FC disruption in fMRI covary with those of LFP and spiking signals in deafferented areas 3b and S2. These findings have important implications for fMRI studies aiming to probe FC alterations in pathological conditions involving deafferentation in humans.SIGNIFICANCE STATEMENT By directly comparing fMRI, local field potential, and spike signals in both tactile stimulation and resting states before and after severe disruption of dorsal column afferent, we demonstrated that reduction in fMRI responses to stimuli is accompanied by weakened resting-state fMRI functional connectivity (FC) in input-deprived and reorganized digit regions in area 3b of the S1 and S2. Concurrent reductions in local field potential and spike FC validated the use of resting-state fMRI signals for probing neural intrinsic FC alterations in pathological deafferented cortex, and indicated that disrupted FC between mesoscale functionally highly related regions may contribute to the behavioral impairments.

Integration of vibrotactile frequency information beyond the mechanoreceptor channel and somatotopy

A wide variety of tactile sensations arise from the activation of several types of mechanoreceptor-afferent channels scattered all over the body, and their projections create a somatotopic map in the somatosensory cortex. Recent findings challenge the traditional view that tactile signals from different mechanoreceptor-channels/locations are independently processed in the brain, though the contribution of signal integration to perception remains obscure. Here we show that vibrotactile frequency perception is functionally enriched by signal integration across different mechanoreceptor channels and separate skin locations. When participants touched two sinusoidal vibrations of far-different frequency, which dominantly activated separate channels with the neighboring fingers or the different hand and judged the frequency of one vibration, the perceived frequency shifted toward the other (assimilation effect). Furthermore, when the participants judged the frequency of the pair as a whole, they consistently reported an intensity-based interpolation of the two vibrations (averaging effect). Both effects were similar in magnitude between the same and different hand conditions and significantly diminished by asynchronous presentation of the vibration pair. These findings indicate that human tactile processing is global and flexible in that it can estimate the ensemble property of a large-scale tactile event sensed by various receptors distributed over the body.

Cortical Representation of Pain and Touch: Evidence from Combined Functional Neuroimaging and Electrophysiology in Non-human Primates

Human functional MRI studies in acute and various chronic pain conditions have revolutionized how we view pain, and have led to a new theory that complex multi-dimensional pain experience (sensory-discriminative, affective/motivational, and cognitive) is represented by concurrent activity in widely-distributed brain regions (termed a network or pain matrix). Despite these breakthrough discoveries, the specific functions proposed for these regions remain elusive, because detailed electrophysiological characterizations of these regions in the primate brain are lacking. To fill in this knowledge gap, we have studied the cortical areas around the central and lateral sulci of the non-human primate brain with combined submillimeter resolution functional imaging (optical imaging and fMRI) and intracranial electrophysiological recording. In this mini-review, I summarize and present data showing that the cortical circuitry engaged in nociceptive processing is much more complex than previously recognized. Electrophysiological evidence supports the engagement of a distinct nociceptive-processing network within SI (i.e., areas 3a, 3b, 1 and 2), SII, and other areas along the lateral sulcus. Deafferentation caused by spinal cord injury profoundly alters the relationships between fMRI and electrophysiological signals. This finding has significant implications for using fMRI to study chronic pain conditions involving deafferentation in humans.

High Spatiotemporal Resolution ECoG Recording of Somatosensory Evoked Potentials with Flexible Micro-Electrode Arrays

Electrocorticogram (ECoG) has great potential as a source signal, especially for clinical BMI. Until recently, ECoG electrodes were commonly used for identifying epileptogenic foci in clinical situations, and such electrodes were low-density and large. Increasing the number and density of recording channels could enable the collection of richer motor/sensory information, and may enhance the precision of decoding and increase opportunities for controlling external devices. Several reports have aimed to increase the number and density of channels. However, few studies have discussed the actual validity of high-density ECoG arrays. In this study, we developed novel high-density flexible ECoG arrays and conducted decoding analyses with monkey somatosensory evoked potentials (SEPs). Using MEMS technology, we made 96-channel Parylene electrode arrays with an inter-electrode distance of 700 μm and recording site area of 350 μm². The arrays were mainly placed onto the finger representation area in the somatosensory cortex of the macaque, and partially inserted into the central sulcus. With electrical finger stimulation, we successfully recorded and visualized finger SEPs with a high spatiotemporal resolution. We conducted offline analyses in which the stimulated fingers and intensity were predicted from recorded SEPs using a support vector machine. We obtained the following results: (1) Very high accuracy (~98%) was achieved with just a short segment of data (~15 ms from stimulus onset). (2) High accuracy (~96%) was achieved even when only a single channel was used. This result indicated placement optimality for decoding. (3) Higher channel counts generally improved prediction accuracy, but the efficacy was small for predictions with feature vectors that included time-series information. These results suggest that ECoG signals with high spatiotemporal resolution could enable greater decoding precision or external device control.

Neural Categorization of Vibrotactile Frequency in Flutter and Vibration Stimulations: An fMRI Study

As the use of wearable haptic devices with vibrating alert features is commonplace, an understanding of the perceptual categorization of vibrotactile frequencies has become important. This understanding can be substantially enhanced by unveiling how neural activity represents vibrotactile frequency information. Using functional magnetic resonance imaging (fMRI), this study investigated categorical clustering patterns of the frequency-dependent neural activity evoked by vibrotactile stimuli with gradually changing frequencies from 20 to 200 Hz. First, a searchlight multi-voxel pattern analysis (MVPA) was used to find brain regions exhibiting neural activities associated with frequency information. We found that the contralateral postcentral gyrus (S1) and the supramarginal gyrus (SMG) carried frequency-dependent information. Next, we applied multidimensional scaling (MDS) to find low-dimensional neural representations of different frequencies obtained from the multi-voxel activity patterns within these regions. The clustering analysis on the MDS results showed that neural activity patterns of 20-100 Hz and 120-200 Hz were divided into two distinct groups. Interestingly, this neural grouping conformed to the perceptual frequency categories found in the previous behavioral studies. Our findings therefore suggest that neural activity patterns in the somatosensory cortical regions may provide a neural basis for the perceptual categorization of vibrotactile frequency.

Rate and timing of cortical responses driven by separate sensory channels

The sense of touch comprises multiple sensory channels that each conveys characteristic signals during interactions with objects. These neural signals must then be integrated in such a way that behaviorally relevant information about the objects is preserved. To understand the process of integration, we implement a simple computational model that describes how the responses of neurons in somatosensory cortex—recorded from awake, behaving monkeys—are shaped by the peripheral input, reconstructed using simulations of neuronal populations that reproduce natural spiking responses in the nerve with millisecond precision. First, we find that the strength of cortical responses is driven by one population of nerve fibers (rapidly adapting) whereas the timing of cortical responses is shaped by the other (Pacinian). Second, we show that input from these sensory channels is integrated in an optimal fashion that exploits the disparate response behaviors of different fiber types.

DOI:http://dx.doi.org/10.7554/eLife.10450.001

Our sense of touch depends upon receptors in our skin that send signals to the brain about the objects with which we interact. Different types of touch receptors respond in different ways when we grasp and manipulate objects; for example, by altering the strength of their response or its timing.

Saal et al. have now investigated how neurons in a part of the brain called the somatosensory cortex, which processes touch signals from the hand, respond to the information from the different receptor types. First, recordings were made of the electrical activity of the touch receptors and the neurons in the brain of monkeys. Using this data, Saal et al. built computer models that allow the response of neurons in the brain to be predicted from the responses of the touch receptors.

The models showed that signals from different types of touch receptors shape the response of neurons in the brain in different ways. One receptor type controls how strong a neuron’s response will be, while another one controls the precise timing of the response. Further investigation revealed that this way of combining the signals from the different receptors preserves as much information as possible about objects and thereby helps the brain to process information acquired by touch quickly and efficiently.

Future experiments will examine how touch is represented in two structures that lie between the receptors and the somatosensory cortex: one in the brainstem, the other in a brain region called the thalamus.

DOI:http://dx.doi.org/10.7554/eLife.10450.002

Fish do not feel pain and its implications for understanding phenomenal consciousness

Phenomenal consciousness or the subjective experience of feeling sensory stimuli is fundamental to human existence. Because of the ubiquity of their subjective experiences, humans seem to readily accept the anthropomorphic extension of these mental states to other animals. Humans will typically extrapolate feelings of pain to animals if they respond physiologically and behaviourally to noxious stimuli. The alternative view that fish instead respond to noxious stimuli reflexly and with a limited behavioural repertoire is defended within the context of our current understanding of the neuroanatomy and neurophysiology of mental states. Consequently, a set of fundamental properties of neural tissue necessary for feeling pain or experiencing affective states in vertebrates is proposed. While mammals and birds possess the prerequisite neural architecture for phenomenal consciousness, it is concluded that fish lack these essential characteristics and hence do not feel pain.

Connectivity of somatosensory cortical area 1 forms an anatomical substrate for the emergence of multifinger receptive fields and complex feature selectivity in the squirrel monkey (Saimiri sciureus)

Converging evidence shows that interaction of digit-specific input, which is required to form global tactile percepts, begins as early as area 3b in the primary somatosensory cortex with the involvement of intrinsic lateral connections. How tactile processing is further elaborated in area 1, the next stage of the somatosensory cortical hierarchy, is less understood. This question was investigated by studying the tangential distribution of intrinsic and interareal connections of finger representations of area 1. Retrogradely labeled cell densities and anterogradely labeled fibers and terminal patches were plotted and quantified with respect to the hand representation by combining tract tracing with electrophysiological mapping and intrinsic signal optical imaging in somatosensory areas. Intrinsic connections of distal finger pad representations of area 1 spanned the representation of multiple digits indicating strong cross-digit connectivity. Area 1 distal finger pad regions also established high-density connections with homotopic regions of areas 3b and 2. Although similar to area 3b, connections of area 1 distributed more widely and covered a larger somatotopic representation including more proximal parts of the finger representations. The lateral connectivity pattern of area 1 is a suitable anatomical substrate of the emergence of multifinger receptive fields, complex feature selectivity, and invariant stimulus properties of the neurons.

Orientation domain diversity in macaque area V2

Single orientation domains in primary (V1) and second (V2) visual cortical areas are known to encode the orientation of visual contours. However, the visual world contains multiple and complex contour types. How do these domains handle such complexity? Using optical imaging methods, we have examined orientation response to two types of contours: real (luminance-defined) and illusory (inferred). We find that, unlike area V1, there are multiple types of orientation domain in V2. These include “real only” domains, “higher-order” domains (which respond to an orientation whether real or illusory), and other domains with nonmatching real/illusory orientation preference. We suggest that this plurality of orientation domains in V2 enables the complexities of local and global contour extraction.

Convergence across tactile afferent types in primary and secondary somatosensory cortices

Integration of information by convergence of inputs onto sensory cortical neurons is a requisite for processing higher-order stimulus features. Convergence across defined peripheral input classes has generally been thought to occur at levels beyond the primary sensory cortex, however recent work has shown that this does not hold for the convergence of slowly-adapting and rapidly-adapting inputs in primary somatosensory cortex. We have used a new analysis method for multi-unit recordings, to show convergence of inputs deriving from the rapidly-adapting and Pacinian channels in a proportion of neurons in both primary and secondary somatosensory cortex in the anaesthetised cat. We have validated this method using single-unit recordings. The secondary somatosensory cortex has a greater proportion of sites that show convergence of this type than primary somatosensory cortex. These findings support the hypothesis that the more complex features processed in higher cortical areas require a greater degree of convergence across input classes, but also shows that this convergence is apparent in the primary somatosensory cortex.

Intrinsic horizontal connections process global tactile features in the primary somatosensory cortex: neuroanatomical evidence

To understand manual tactile functions in primates, it is essential to explore the interactions between the finger pad representations in somatosensory cortex. To this end, we used optical imaging and electrophysiological mapping to guide neuroanatomical tracer injections into distal digit tip representations of Brodmann area 3b in the squirrel monkey. Retrogradely labeled cell densities and anterogradely labeled fibers and terminal patches in somatosensory areas were plotted and quantified with respect to tangential distribution. Within area 3b, reciprocal patchy distribution of anterograde and retrograde labeling spanned the representation of the distal pad of multiple digits, indicating strong cross-digit connectivity. Inter-areal connections revealed bundles of long-range fibers projecting anteroposteriorly, connecting area 3b with clusters of labeled neurons and terminal axon arborizations in area 1. Inter-areal linkage appeared to be largely confined to the representation of the injected finger. These findings provide the neuroanatomical basis for the interaction between distal finger pad representations observed by recent electrophysiological studies. We propose that intra-areal connectivity may be heavily involved in interdigit integration such as shape discrimination, whereas long-range inter-areal connections may subserve active touch in a digit-specific manner.

Temporal jitter of the BOLD signal reveals a reliable initial dip and improved spatial resolution

fMRI, one of the most important noninvasive brain imaging methods, relies on the blood oxygen level-dependent (BOLD) signal, whose precise underpinnings are still not fully understood. It is a widespread assumption that the components of the hemodynamic response function (HRF) are fixed relative to each other in time, leading most studies as well as analysis tools to focus on trial-averaged responses, thus using or estimating a condition- or location-specific "canonical HRF". In the current study, we examined the nature of the variability of the BOLD response and asked in particular whether the positive BOLD peak is subject to trial-to-trial temporal jitter. Our results show that the positive peak of the stimulus-evoked BOLD signal exhibits a trial-to-trial temporal jitter on the order of seconds. Moreover, the trial-to-trial variability can be exploited to uncover the initial dip in the majority of voxels by pooling trial responses with large peak latencies. Initial dips exposed by this procedure possess higher spatial resolution compared to the positive BOLD signal in the human visual cortex. These findings allow for the reliable observation of fMRI signals that are physiologically closer to neural activity, leading to improvements in both temporal and spatial resolution.

The organization of submodality-specific touch afferent inputs in the vibrissa column

The rodent tactile vibrissae are innervated by several different types of touch sensory neurons. The central afferents of all touch neurons from one vibrissa collectively project to a columnar structure called a barrelette in the brainstem. Delineating how distinct types of sensors connect to second-order neurons within each barrelette is critical for understanding tactile information coding and processing. Using genetic and viral techniques, we labeled slowly adapting (SA) mechanosensory neurons, rapidly adapting (RA) mechanosensory neurons, afferent synapses, and second-order projection neurons with four different fluorescent markers to examine their connectivity. We discovered that within each vibrissa column, individual sensory neurons project collaterals to multiply distributed locations, inputs from SA and RA afferents are spatially intermixed without any discernible stereotypy or topography, and second-order projection neurons receive convergent SA and RA inputs. Our findings reveal a "one-to-many and many-to-one" connectivity scheme and the circuit architecture for tactile information processing at the first-order synapses.

Optical imaging of cortical networks via intracortical microstimulation

Understanding cortical organization is key to understanding brain function. Distinct neural networks underlie the functional organization of the cerebral cortex; however, little is known about how different nodes in the cortical network interact during perceptual processing and motor behavior. To study cortical network function we examined whether the optical imaging of intrinsic signals (OIS) reveals the functional patterns of activity evoked by electrical cortical microstimulation. We examined the effects of current amplitude, train duration, and depth of cortical stimulation on the hemodynamic response to electrical microstimulation (250-Hz train, 0.4-ms pulse duration) in anesthetized New World monkey somatosensory cortex. Electrical stimulation elicited a restricted cortical response that varied according to stimulation parameters and electrode depth. Higher currents of stimulation recruited more areas of cortex than smaller currents. The largest cortical responses were seen when stimulation was delivered around cortical layer 4. Distinct local patches of activation, highly suggestive of local projections, around the site of stimulation were observed at different depths of stimulation. Thus we find that specific electrical stimulation parameters can elicit activation of single cortical columns and their associated columnar networks, reminiscent of anatomically labeled networks. This novel functional tract tracing method will open new avenues for investigating relationships of local cortical organization.

Frequency-dependent patterns of somatosensory cortical responses to vibrotactile stimulation in humans: a fMRI study

In the human mechanosensation system, rapidly adapting afferents project sensory signals of flutter (5-50Hz) to the contralateral primary somatosensory cortex (S1) and bilateral secondary somatosensory cortex (S2) whereas Pacinian afferents project sensory signals of vibration (50-400Hz) to bilateral S2. However, it remains largely unknown how somatosensory cortical activity changes as a function of vibrotactile frequency. This functional magnetic resonance imaging (fMRI) study investigated frequency dependency of somatosensory cortical activity in humans by applying vibrotactile stimulation with various frequencies (20-200Hz) to the index finger. We found more frequency-dependent voxels in the upper bank of the lateral sulcus (LS) of S2 than in S1 and the posterior parietal cortex of S2. Our statistical spatial clustering analysis showed that two groups of positively or negatively frequency-dependent voxels formed distinct clusters, most clearly in the LS. Using a cortical separability index, we reaffirmed that somatosensory cortical activity was most separable at 50Hz, previously known to demarcate flutter and vibration. Our results suggest that the LS (S2) may play an important role in processing vibrotactile frequency information and that the somatosensory cortex may include spatially localized neural assemblies specialized to higher or lower vibrotactile frequency.

Crossmodal interactions of haptic and visual texture information in early sensory cortex

Both visual and haptic information add to the perception of surface texture. While prior studies have reported crossmodal interactions of both sensory modalities at the behavioral level, neuroimaging studies primarily investigated texture perception in separate visual and haptic paradigms. These experimental designs, however, only allowed to identify overlap in both sensory processing streams but no interaction of visual and haptic texture processing. By varying texture characteristics in a bimodal task, the current study investigated how these crossmodal interactions are reflected at the cortical level. We used fMRI to compare cortical activation in response to matching versus non-matching visual-haptic texture information. We expected that passive simultaneous presentation of matching visual-haptic input would be sufficient to induce BOLD responses graded with varying texture characteristics. Since no cognitive evaluation of the stimuli was required, we expected to find changes primarily at a rather early processing stage. Our results confirmed our assumptions by showing crossmodal interactions of visual-haptic texture information in early somatosensory and visual cortex. However, the nature of the crossmodal effects was slightly different in both sensory cortices. In early visual cortex, matching visual-haptic information increased the average activation level and induced parametric BOLD signal variations with varying texture characteristics. In early somatosensory cortex only the latter was true. These results challenge the notion that visual and haptic texture information is processed independently and indicate a crossmodal interaction of sensory information already at an early cortical processing stage.

Dynamic causal modeling suggests serial processing of tactile vibratory stimuli in the human somatosensory cortex--an fMRI study

Sensitivity to location and frequency of tactile stimuli is a characterizing feature of human primary (S1), and secondary (S2) somatosensory cortices. Recent evidence suggests that S1 is predominantly receptive to stimulus location, while S2 is attuned to stimulus frequency. Although it is well established in humans that tactile frequency information is relayed serially from S1 to S2, a recent study, using functional magnetic resonance imaging (fMRI) in combination with dynamic causal modeling (DCM), suggested that somatosensory inputs are processed in parallel in S1 and S2. In the present fMRI/DCM study, we revisited this controversy and investigated the specialization of the human somatosensory cortical areas with regard to tactile stimulus representations, as well as their effective connectivity. During brain imaging, 14 participants performed a somatosensory discrimination task on vibrotactile stimuli. Importantly, the model space for DCM was chosen to allow for direct inference on the question of interest by systematically varying the information transmission from pure parallel to pure serial implementations. Bayesian model comparison on the level of model families strongly favors a serial, instead of a parallel processing route for tactile stimulus information along the somatosensory pathway. Our fMRI/DCM data thus support previous suggestions of a sequential information transmission from S1 to S2 in humans.

Dynamic reorganization of digit representations in somatosensory cortex of nonhuman primates after spinal cord injury

Somatosensory cortices of adult primates reactivate over time after sensory loss. The time course and the neural mechanisms underlying the cortical reactivation are not well understood. Here we report that longitudinal high-resolution functional magnetic resonance imaging (fMRI) studies on anesthetized squirrel monkeys revealed dynamic reorganizations of digit activations in area 3b, within 2 months after severely disrupting afferent inputs by dorsal column section. We found that digit regions in which inputs were severely disrupted exhibited fMRI tactile responses. Reorganization was characterized by an early moving away phase and a late returning phase, as indicated by spatial shifts of individual digit activation centers in relation to the pre-lesion activation sites. Subsequent optical imaging studies confirmed fMRI activations, and dense microelectrode penetrations identified weak neuronal activity at the reactivated sites. Activation zones detected by fMRI and optical imaging were significantly larger in input-deprived than normal input single-digit regions and were larger than regions defined by neuronal spiking activity. This study captures the dynamic reorganization of digit representations after dorsal column lesions and reveals differences between functional imaging and microelectrode recording maps. Our observations suggest that subthreshold activity plays an important role in the reactivation of deafferented cortex and could promote behavioral recovery.

Anesthesia and the quantitative evaluation of neurovascular coupling

Anesthesia has broad actions that include changing neuronal excitability, vascular reactivity, and other baseline physiologies and eventually modifies the neurovascular coupling relationship. Here, we review the effects of anesthesia on the spatial propagation, temporal dynamics, and quantitative relationship between the neural and vascular responses to cortical stimulation. Previous studies have shown that the onset latency of evoked cerebral blood flow (CBF) changes is relatively consistent across anesthesia conditions compared with variations in the time-to-peak. This finding indicates that the mechanism of vasodilation onset is less dependent on anesthesia interference, while vasodilation dynamics are subject to this interference. The quantitative coupling relationship is largely influenced by the type and dosage of anesthesia, including the actions on neural processing, vasoactive signal transmission, and vascular reactivity. The effects of anesthesia on the spatial gap between the neural and vascular response regions are not fully understood and require further attention to elucidate the mechanism of vascular control of CBF supply to the underlying focal and surrounding neural activity. The in-depth understanding of the anesthesia actions on neurovascular elements allows for better decision-making regarding the anesthetics used in specific models for neurovascular experiments and may also help elucidate the signal source issues in hemodynamic-based neuroimaging techniques.

Second-order receptive fields reveal multidigit interactions in area 3b of the macaque monkey

Linear receptive field (RF) models of area 3b neurons reveal a three-component structure: a central excitatory region flanked by two inhibitory regions that are spatially and temporally nonoverlapping with the excitation. Previous studies also report that there is an "infield" inhibitory region throughout the neuronal RF, which is a nonlinear interactive (second order) effect whereby stimuli lagging an input to the excitatory region are suppressed. Thus linear models may be inaccurate approximations of the neurons' true RFs. In this study, we characterize the RFs of area 3b neurons, using a second-order quadratic model. Data were collected from 80 neurons of two awake, behaving macaque monkeys while a random dot pattern was scanned simultaneously across the distal pads of digits D2, 3, and 4. We used an iterative method derived from matching pursuit to identify a set of linear and nonlinear terms with significant effects on the neuronal response. For most neurons (65/80), the linear component of the quadratic RF was characterized by a single excitatory region on the dominant digit. Interactions within the dominant digit were characterized by two quadratic filters that capture the spatial aspects of the interactive infield inhibition. Interactions between the dominant (most responsive) digit and its adjacent digit(s) formed the largest class of cross-digit interactions. The results demonstrate that a significant part of area 3b responses is due to nonlinear mechanisms, and furthermore, the data support the notion that area 3b neurons have "nonclassical RF"-like input from adjacent fingers, indicating that area 3b plays a role in integrating shape inputs across digits.

Differential fMRI activation to noxious heat and tactile stimuli in parasylvian areas of new world monkeys

Emerging evidence supports an important role of posterior parasylvian areas in both pain and touch processing. Whether there are separate or shared networks for these sensations remains controversial. The present study compared spatial patterns of brain activation in response to unilateral nociceptive heat (47.5°C) or innocuous tactile stimulation (8-Hz vibration) to digits through high-resolution functional magnetic resonance imaging (fMRI) in squirrel monkeys. In addition, the temporal profile of heat-stimulus-evoked fMRI Blood Oxygenation Level Dependent (BOLD) signal changes was characterized. By examining high-resolution fMRI and histological measures at both the individual and the group levels, we found that both nociceptive heat and tactile stimuli elicited activation in bilateral secondary somatosensory and ventral parietal areas (S2/PV) and in ipsilateral ventral somatosensory areas (VS) and retroinsula (Ri). Bilateral posterior insular cortex (pIns) and area 7b responded preferentially to nociceptive heat stimulation. Single voxels within each activation cluster showed robust BOLD signal changes during each block of nociceptive stimulation. Across animals (n=11), nociceptive response magnitudes of contralateral VS and pIns and ipsilateral Ri were significantly greater than corresponding areas in the opposite hemisphere. In sum, both distinct and shared areas in regions surrounding the posterior sylvian fissure were activated in response to nociceptive and tactile inputs in nonhuman primates.

Optical imaging in galagos reveals parietal-frontal circuits underlying motor behavior

The posterior parietal cortex (PPC) of monkeys and prosimian galagos contains a number of subregions where complex, behaviorally meaningful movements, such as reaching, grasping, and body defense, can be evoked by electrical stimulation with long trains of electrical pulses through microelectrodes. Shorter trains of pulses evoke no or simple movements. One possibility for the difference in effectiveness of intracortical microstimulation is that long trains activate much larger regions of the brain. Here, we show that long-train stimulation of PPC does not activate widespread regions of frontal motor and premotor cortex but instead, produces focal, somatotopically appropriate activations of frontal motor and premotor cortex. Shorter stimulation trains activate the same frontal foci but less strongly, showing that longer stimulus trains do not produce less specification. Because the activated sites in frontal cortex correspond to the locations of direct parietal-frontal anatomical connections from the stimulated PPC subregions, the results show the usefulness of optical imaging in conjunction with electrical stimulation in showing functional pathways between nodes in behavior-specific cortical networks. Thus, long-train stimulation is effective in evoking ethologically relevant sequences of movements by activating nodes in a cortical network for a behaviorally relevant period rather than spreading activation in a nonspecific manner.

Intra- and inter-subject variability of high field fMRI digit maps in somatosensory area 3b of new world monkeys

This study evaluates the intra- and inter-subject variability of digit maps in area 3b of anesthetized squirrel monkeys. Maps were collected using high field blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI). BOLD responses to individual digit stimulations were mapped and their response properties (location, area of activation, % signal change, time to peak response) were compared within and across imaging sessions separated by up to 20 months. During single digit stimulation using a block design, the spatiotemporal response of the BOLD signal for individual runs within and across sessions and animals was well conserved, with a time to peak BOLD response of 20+/-4 s. The variability in the center of BOLD activation in area 3b was 0.41+/-0.24 mm (mean+/-SD) across individual 5-7 min runs within a scanning session and 0.55+/-0.15 mm across sessions. The average signal change across all animals, runs and sessions was 0.62+/-0.38%, and varied 32% within and 40% across sessions. In a comparison of the stability and reproducibility of the area of single digit activation obtained using three approaches, use of a fixed statistical threshold (P<10(-5)) yielded an average area of 4.8+/-3.5 mm(2) (mean+/-SD), adaptive statistical thresholding 1.32+/-1.259 mm(2) (mean+/-SD), and combined fixed statistical and adaptive BOLD signal amplitude 4.4+/-2.5 mm(2) (mean+/-SD) across image runs and sessions. The somatotopic organization was stable within animals across sessions, while across animals, there was some variation in overall activation pattern and inter-digit distances. These results confirm that BOLD activation maps of single digits in area 3b as characterized by activation center, signal amplitudes, and temporal profile are very stable. The activation sizes determined by various criteria are the most variable measure in this preparation, but adaptive statistical thresholding appears to yield the most stable and reproducible maps. This study serves as a baseline assessment of the limits imposed on the detection of plastic changes by experimental variations of the digit BOLD fMRI activation maps in normal animals, and as an indicator of the likely performance limits in human studies.

Optical imaging of digit topography in individual awake and anesthetized squirrel monkeys

Topographic maps and columnar structures are fundamental to cortical sensory information processing. Most of the knowledge about detailed topographic maps and columnar structure comes mainly from experiments conducted on anesthetized animals. Towards the goal of evaluating whether topographic maps change with respect to behavioral demands, we used intrinsic signal optical imaging in alert monkeys to examine the spatial specificity of cortical topographic representation. Specifically, the somatotopies of neighboring distal finger pad representation in areas 3b and 1 were examined in the same awake and anesthetized squirrel monkey. In comparison to the anesthetized animal, we found larger cortical activation sizes in the alert animal in area 3b, where activation widths were found to overlap with even non-adjacent digits. This may suggest that in the alert animal, there is less inhibition across the somatotopic map within area 3b.

Convergence of submodality-specific input onto neurons in primary somatosensory cortex

At the somatosensory periphery, slowly adapting type 1 (SA1) and rapidly adapting (RA) afferents respond very differently to step indentations: SA1 afferents respond throughout the entire stimulus interval (sustained response), whereas RA afferents respond only at stimulus onset (on response) and offset (off response). We recorded the responses of cortical neurons to step indentations and found many neurons in areas 3b and 1 to exhibit properties that are intermediate between these two extremes: These neurons responded during the sustained portion of the stimulus and also at the offset of the stimulus. Several lines of evidence indicate that these neurons, which exist in large proportions even at these early stages of somatosensory cortical processing, receive input from both populations of afferents. First, we show that many cortical neurons have both a significant sustained response and a significant off response. Second, the strength of the off response is uncorrelated with that of the sustained response, which is to be expected if sustained and off responses stem from different populations of afferent fibers. Third, the bulk of the variance in cortical responses to step indentations can be accounted for using a linear combination of both SA1 and RA responses. Finally, we show that the off response in cortical neurons does not reflect rebound from inhibition. We conclude that the convergence of modality specific input onto individual neurons is common in primary somatosensory cortex and discuss how this conclusion might be reconciled with previous findings.

Area-specific representation of mechanical nociceptive stimuli within SI cortex of squirrel monkeys

While functional imaging studies in humans have consistently reported activation of primary somatosensory cortex (SI) with painful stimuli, the specific roles of subdivisions of areas 3a, 3b, and 1 within SI during pain perception are largely unknown, particularly in the representation of mechanical evoked pain. In this study, we investigated how modality, location, and intensity of nociceptive stimuli are represented within SI by using high-spatial resolution optical imaging of intrinsic signals in Pentothal-anesthetized squirrel monkeys. Perceptually comparable mechanical nociceptive and innocuous tactile stimuli were delivered by indenting the glabrous skin of the distal finger pads with 0.2 and 2mm diameter probes, respectively. Within each of areas 3a, 3b, and 1, activations to mechanical nociceptive stimulation of individual distal finger pads were spatially distinct and somatotopically organized. We observed differential cortical activation patterns. Areas 3a, 3b, and 1 were all activated during mechanical nociceptive stimulation and were modulated by nociceptive stimulus intensity. However, with innocuous tactile stimulation, mainly areas 3b and 1 exhibited response modulation with different levels of stimulation. In summary, mechanical nociceptive inputs are area-specific and topographically represented within SI. We propose that all areas of SI are implicated in encoding the features of mechanical nociception, where areas 3a and 3b are distinctively involved in coding nociceptive and pressure sensation components of stimulation.

Responses of areas 3b and 1 in anesthetized squirrel monkeys to single- and dual-site stimulation of the digits

Stimulation of the skin evokes topographically organized activation in somatosensory cortex. This representation is context dependent, however, since a different cortical topography is observed in area 3b when stimulated with complex tactile stimuli that evoke the von Békésy funneling illusion. Here we report on the population responses, as observed with intrinsic optical imaging, of area 1 and area 3b in the anesthetized squirrel monkey to pressure indentation of distal finger pads. Individual finger pad stimulation revealed that area 1 exhibited a smaller magnification factor than 3b, as evidenced by a smaller area of activation elicited by distal finger pad stimulation. Effects of paired finger pad stimulation produced largely similar effects in area 1 and area 3b. Paired finger pad stimulation produced reductions in the area of digit activation in area 1, suggesting the presence of lateral inhibition and funneling of information in area 1. Suppressive effects were stronger for paired stimulations at adjacent than at nonadjacent sites. Single-unit recordings revealed a mixture of either a summation or a suppression of the response to paired finger stimulation, compared with single finger pad stimulation of the primary digit. However, the average population response showed that paired finger pad stimulation resulted in response suppression. Based on this study and previous studies, we suggest the presence of at least three distinct ranges of lateral inhibition in areas 3b and 1.

Whole-hand sensorimotor area: cortical stimulation localization and correlation with functional magnetic resonance imaging

OBJECT

The pli de passage moyen (PPM) is an omega-shaped cortical landmark bulging into the central sulcus. There has been considerable interest in the PPM given that hand motor and sensory tasks have been found on functional magnetic resonance (fMR) imaging to activate the structure. Note, however, that the cortical function subserved by the PPM is not completely understood. Finger and thumb function are somatotopically organized over the central area and encompass a larger cortical surface than the anatomical PPM. Therefore, a sensory or motor hand area within the PPM would be redundant with the somatotopically organized digit function in the primary sensorimotor cortex. In this study the authors aimed to clarify the function subserved by the PPM and further evaluate hand area function in the primary sensorimotor cortex.

METHODS

To further elucidate the function subserved by the PPM, patients underwent cortical stimulation in the region of the PPM as well as fMR imaging-demonstrated activation of the hand area. Two separate analytical methods were used to correlate hand area functional imaging with whole-hand sensory and motor responses induced by cortical stimulation.

RESULTS

A relationship of the anatomical PPM with cortical stimulation responses as well as hand fMR imaging activation was observed.

CONCLUSIONS

A strong relationship was identified between the PPM, whole-hand sensory and motor stimulation responses, and fMR imaging hand activation. Whole-hand motor and whole-hand sensory cortical regions were identified in the primary sensorimotor cortex. It was localized to the PPM and exists in addition to the somatotopically organized finger and thumb sensory and motor areas.

High-resolution fMRI maps of cortical activation in nonhuman primates: correlation with intrinsic signal optical images

One of the most widely used functional brain mapping tools is blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI). This method has contributed to new understandings of the functional roles of different areas in the human brain. However, its ability to map cerebral cortex at high spatial (submillimeter) resolution is still unknown. Other methods such as single- and multiunit electrophysiology and intrinsic signal optical imaging have revealed submillimeter resolution of sensory topography and cortical columnar activations. However, they are limited either by spatial scale (electrophysiology characterizes only local groups of neurons) or by the inability to monitor deep structures in the brain (i.e., cortical regions buried in sulci or subcortical structures). A method that could monitor all regions of the brain at high spatial resolution would be ideal. This capacity would open the doors to investigating, for example, how networks of cerebral cortical columns relate to or produce behavior. In this article we demonstrate that, without benefit of contrast agents, at a magnetic field strength of 9.4 tesla, BOLD fMRI can reveal millimeter-sized topographic maps of digit representation in the somatosensory cortex of the anesthetized squirrel monkey. Furthermore, by mapping the "funneling illusion," it is possible to detect even submillimeter shifts in activation in the cortex. Our data suggest that at high magnetic field strength, the positive BOLD signal can be used to reveal high spatial resolution maps of brain activity, a finding that weakens previous notions about the ultimate spatial specificity of the positive BOLD signal.

High-resolution maps of real and illusory tactile activation in primary somatosensory cortex in individual monkeys with functional magnetic resonance imaging and optical imaging

Although blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) has been widely used to explore human brain function, questions remain regarding the ultimate spatial resolution of positive BOLD fMRI, and indeed the extent to which functional maps revealed by positive BOLD correlate spatially with maps obtained with other high-spatial-resolution mapping techniques commonly used in animals, such as optical imaging of intrinsic signal (OIS) and single-unit electrophysiology. Here, we demonstrate that the positive BOLD signal at 9.4T can reveal the fine topography of individual fingerpads in single-condition activation maps in nonhuman primates. These digit maps are similar to maps obtained from the same animal using intrinsic optical imaging. Furthermore, BOLD fMRI reliably resolved submillimeter spatial shifts in activation in area 3b previously identified with OIS (Chen et al., 2003) as neural correlates of the "funneling illusion." These data demonstrate that at high field, high-spatial-resolution topographic maps can be achieved using the positive BOLD signal, weakening previous notions regarding the spatial specificity of the positive BOLD signal.

Dependence of BOLD signal change on tactile stimulus intensity in SI of primates

Recently, we have demonstrated that the fine-digit topography (millimeter sized) previously identified in the primary somatosensory cortex (SI), using electrophysiology and intrinsic signal optical imaging, can also be mapped with submillimeter resolution using blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging at high field. In the present study, we have examined the dependence of BOLD signal response on stimulus intensity in two subregions of SI, Areas 3b and 1. In a region(s)-of-interest (ROI) analysis of Area 3b, BOLD signal amplitude increased linearly with increasing amplitude of an 8-Hz vibrotactile stimulus, and BOLD signal was sustained throughout the stimulation period. In contrast, in Area 1, a significant BOLD signal response was only observed with more intense stimuli, and ROI analysis of the dependence of BOLD response showed no significant dependence on stimulus intensity. In addition, activation was not sustained throughout the period of stimulation. Differing responses of Areas 3b and 1 suggest potentially divergent roles for subregions of SI cortices in vibrotactile intensity encoding. Moreover, this study underscores the importance of imaging at small spatial scales. In this case, such high-resolution imaging allows differentiation between area-specific roles in intensity encoding and identifies anatomic targets for detailed electrophysiological studies of somatosensory neuronal populations with different coding properties. These experiments illustrate the value of nonhuman primates for characterizing the dependence of the BOLD signal response on stimulus parameters and on underlying neural response properties.

Long-term optical imaging of intrinsic signals in anesthetized and awake monkeys

Some exciting new efforts to use intrinsic signal optical imaging methods for long-term studies in anesthetized and awake monkeys are reviewed. The development of such methodologies opens the door for studying behavioral states such as attention, motivation, memory, emotion, and other higher-order cognitive functions. Long-term imaging is also ideal for studying changes in the brain that accompany development, plasticity, and learning. Although intrinsic imaging lacks the temporal resolution offered by dyes, it is a high spatial resolution imaging method that does not require application of any external agents to the brain. The bulk of procedures described here have been developed in the monkey but can be applied to the study of surface structures in any in vivo preparation.

Optical imaging of contrast response in Macaque monkey V1 and V2

Our studies on brightness information processing in Macaque monkey visual cortex suggest that the thin stripes in the secondary visual area (V2) are preferentially activated by brightness stimuli (such as full field luminance modulation and illusory edge-induced brightness modulation). To further examine this possibility, we used intrinsic signal optical imaging to examine contrast response of different functional domains in primary and secondary visual areas (V1 and V2). Color and orientation stimuli were used to map functional domains in V1 (color domains, orientation domains) and V2 (thin stripes, thick/pale stripes). To examine contrast response, sinusoidal gratings at different contrasts and spatial frequencies were presented. We find that, consistent with previous studies, the optical signal increased systematically with contrast level. Unlike single-unit responses, optical signals for both color domains and orientation domains in V1 exhibit linear contrast response functions, thereby providing a large dynamic range for V1 contrast response. In contrast to domains in V1, domains in V2 exhibit nonlinear responses, characterized by high gain at low contrasts, saturating at a mid-high contrast levels. At high contrasts, thin stripes exhibit increasing response, whereas thick/pale stripes saturate, consistent with a strong parvocellular input to thin stripes. These findings suggest that, with respect to contrast encoding, thin stripes have a larger dynamic range than thick/pale stripes and further support a role for thin stripes in processing of brightness information.

Optical imaging of SI topography in anesthetized and awake squirrel monkeys

Orderly topographic maps in the primary somatosensory cortex (SI) serve as an anchor for our understanding of somatosensory cortical organization. However, this view is mostly based on data collected in the anesthetized animal. Less is known about these topographies in the awake primate. Even less is known about the relative activations of different subdivisions of SI (areas 3a, 3b, 1, and 2). Toward the goal of understanding the functional activation of SI, we conducted intrinsic signal optical imaging of areas 3b and 1 in awake squirrel monkeys. Monkeys were imaged repeatedly for a period of >2 years in awake and anesthetized states in response to vibrotactile and electrocutaneous stimuli presented to individual fingerpads. During this period, we found stable somatotopic maps in both the anesthetized and awake states, consistent with electrophysiologically recorded maps in areas 3b and 1 in the anesthetized state. In the awake animal, signal sizes were larger, but variability was greater, leading to decreased signal-to-noise ratios. Topographic activations were larger (in both area and amplitude) in the awake animal, suggesting either a less precise topography and/or more complex integration. This brings into question the role of a precise topographic map during behavior. In addition, whereas in the anesthetized animal strongest imaging signals were obtained from area 3b, in the awake animal, area 1 activation dominated over that in area 3b. Differences in relative dominance of area 3b versus area 1 suggest that inter-areal interactions in the alert animal differ substantially from that in the anesthetized animal.

Stimulus-dependent spatial patterns of response in SI cortex

BACKGROUND

Recently we reported that vibrotactile flutter stimulation of a skin locus at different amplitudes evokes an optical response confined to the same local region of the primary somatosensory cortex (SI), where its overall magnitude varies proportionally to the flutter amplitude. In this report, we characterize the impact of the flutter amplitude on the spatial patterns of activity evoked within the responding SI region.

RESULTS

In order to characterize the spatial pattern of activity within the responding SI region, images of the flutter-evoked SI optical response were segmented and analyzed with spatial frequency analysis. The analysis revealed that: (1) dominant spatial frequencies in the optical intrinsic signal emerge within the responding SI region within 3-5 sec of stimulus onset; (2) the stimulus-evoked activity is spatially organized in a form of several roughly parallel, anterior-posteriorly extended waves, spaced 0.4-0.5 mm apart; (3) the waves themselves exhibit spatial periodicities along their long axis; and (4) depending on the flutter stimulus amplitude, these periodicities can range from fine 0.15 mm "ripples" at 50 microm amplitude to well-developed 0.5 mm fluctuations at the amplitude of 400 microm.

CONCLUSION

The observed spatiointensive fractionation on a sub-macrocolumnar scale of the SI response to skin stimulation might be the product of local competitive interactions within the stimulus-activated SI region and may be a feature that could yield novel insights into the functional interactions that take place in SI cortex.