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

Detection and characterization of resting state functional networks in squirrel monkey brain

Resting-state fMRI based on analyzing BOLD signals is widely used to derive functional networks in the brain and how they alter during disease or injury conditions. Resting-state networks can also be used to study brain functional connectomes across species, which provides insights into brain evolution. The squirrel monkey (SM) is a non-human primate (NHP) that is widely used as a preclinical model for experimental manipulations to understand the organization and functioning of the brain. We derived resting-state networks from the whole brain of anesthetized SMs using Independent Component Analysis of BOLD acquisitions. We detected 15 anatomically constrained resting-state networks localized in the cortical and subcortical regions as well as in the white-matter. Networks encompassing visual, somatosensory, executive control, sensorimotor, salience and default mode regions, and subcortical networks including the Hippocampus-Amygdala, thalamus, basal-ganglia and brainstem region correspond well with previously detected networks in humans and NHPs. The connectivity pattern between the networks also agrees well with previously reported seed-based resting-state connectivity of SM brain. This study demonstrates that SMs share remarkable homologous network organization with humans and other NHPs, thereby providing strong support for their suitability as a translational animal model for research and additional insight into brain evolution across species.

Disrupting nociceptive information processing flow through transcranial focused ultrasound neuromodulation of thalamic nuclei

BACKGROUND

MRI-guided transcranial focused ultrasound (MRgFUS) as a next-generation neuromodulation tool can precisely target and stimulate deep brain regions with high spatial selectivity. Combined with MR-ARFI (acoustic radiation force imaging) and using fMRI BOLD signal as functional readouts, our previous studies have shown that low-intensity FUS can excite or suppress neural activity in the somatosensory cortex.

OBJECTIVE

To investigate whether low-intensity FUS can suppress nociceptive heat stimulation-induced responses in thalamic nuclei during hand stimulation, and to determine how this suppression influences the information processing flow within nociception networks.

FINDINGS

BOLD fMRI activations evoked by 47.5 °C heat stimulation of hand were detected in 24 cortical regions, which belong to sensory, affective, and cognitive nociceptive networks. Concurrent delivery of low-intensity FUS pulses (650 kHz, 550 kPa) to the predefined heat nociceptive stimulus-responsive thalamic centromedial_parafascicular (CM_para), mediodorsal (MD), ventral_lateral (VL_ and ventral_lateral_posteroventral (VLpv) nuclei suppressed their heat responses. Off-target cortical areas exhibited reduced, enhanced, or no significant fMRI signal changes, depending on the specific areas. Differentiable thalamocortical information flow during the processing of nociceptive heat input was observed, as indicated by the time to reach 10% or 30% of the heat-evoked BOLD signal peak. Suppression of thalamic heat responses significantly altered nociceptive processing flow and direction between the thalamus and cortical areas. Modulation of contralateral versus ipsilateral areas by unilateral thalamic activity differed. Signals detected in high-order cortical areas, such as dorsal frontal (DFC) and ventrolateral prefrontal (vlPFC) cortices, exhibited faster response latencies than sensory areas.

CONCLUSIONS

The concurrent delivery of FUS suppressed nociceptive heat response in thalamic nuclei and disrupted the nociceptive network. This study offers new insights into the causal functional connections within the thalamocortical networks and demonstrates the modulatory effects of low-intensity FUS on nociceptive information processing.

Infrared neural stimulation in human cerebral cortex

BACKGROUND

Modulation of brain circuits by electrical stimulation has led to exciting and powerful therapies for diseases such as Parkinson's. Because human brain organization is based in mesoscale (millimeter-scale) functional nodes, having a method that can selectively target such nodes could enable more precise, functionally specific stimulation therapies. Infrared Neural Stimulation (INS) is an emerging stimulation technology that stimulates neural tissue via delivery of tiny heat pulses. In nonhuman primates, this optical method provides focal intensity-dependent stimulation of the brain without tissue damage. However, whether INS application to the human central nervous system (CNS) is similarly effective is unknown.

OBJECTIVE

To examine the effectiveness of INS on human cerebral cortex in intraoperative setting and to evaluate INS damage threshholds.

METHODS

Five epileptic subjects undergoing standard lobectomy for epilepsy consented to this study. Cortical response to INS was assessed by intrinsic signal optical imaging (OI, a method that detects changes in tissue reflectance due to neuronal activity). A custom integrated INS and OI system was developed specifically for short-duration INS and OI acquisition during surgical procedures. Single pulse trains of INS with intensities from 0.2 to 0.8 J/cm² were delivered to the somatosensory cortex and responses were recorded via optical imaging. Following tissue resection, histological analysis was conducted to evaluate damage threshholds.

RESULTS

As assessed by OI, and similar to results in monkeys, INS induced responses in human cortex were highly focal (millimeter sized) and led to relative suppression of nearby cortical sites. Intensity dependence was observed at both stimulated and functionally connected sites. Histological analysis of INS-stimulated human cortical tissue provided damage threshold estimates.

CONCLUSION

This is the first study demonstrating application of INS to human CNS and shows feasibility for stimulating single cortical nodes and associated sites and provided INS damage threshold estimates for cortical tissue. Our results suggest that INS is a promising tool for stimulation of functionally selective mesoscale circuits in the human brain, and may lead to advances in the future of precision medicine.

Differential dose responses of transcranial focused ultrasound at brain regions indicate causal interactions

We have previously shown that focused ultrasound (FUS) pulses in low pressure range exerted bidirectional and brain state-dependent neuromodulation in the nonhuman primate somatosensory cortices by fMRI. Here we aim to gain insights about the proposed neuron selective modulation of FUS and probe feedforward versus feedback interactions by simultaneously quantifying the stimulus (FUS pressures: 925, 425, 250 kPa) and response (% BOLD fMRI changes) function at the targeted area 3a/3b and off-target cortical areas at 7T. In resting-state, lowered intensities of FUS resulted in decreased fMRI signal changes at the target area 3a/3b and off-target area 1/2, S2, MCC, insula and auditory cortex, and no signal difference in thalamic VPL and MD nuclei. In activated states, concurrent high-intensity FUS significantly enhanced touch-evoked signals in area 1/2. Medium- and low-intensity FUS significantly suppressed touch-evoked BOLD signals in all areas except in the auditory cortex, VPL and MD thalamic nuclei. Distinct state dependent and dose-response curves led us to hypothesize that FUS's neuromodulatory effects may be mediated through preferential activation of different populations of neurons. Area 3a/3b may have distinct causal feedforward and feedback interactions with Area 1/2, S2, MCC, insula, and VPL. FUS offers a noninvasive neural stimulation tool for dissecting brain circuits and probing causal functional connections.

Critical factors in achieving fine-scale functional MRI: Removing sources of inadvertent spatial smoothing

Ultra‐high Field (≥7T) functional magnetic resonance imaging (UHF‐fMRI) provides opportunities to resolve fine‐scale features of functional architecture such as cerebral cortical columns and layers, in vivo. While the nominal resolution of modern fMRI acquisitions may appear to be sufficient to resolve these features, several common data preprocessing steps can introduce unwanted spatial blurring, especially those that require interpolation of the data. These resolution losses can impede the detection of the fine‐scale features of interest. To examine quantitatively and systematically the sources of spatial resolution losses occurring during preprocessing, we used synthetic fMRI data and real fMRI data from the human visual cortex—the spatially interdigitated human V2 “thin” and “thick” stripes. The pattern of these cortical columns lies along the cortical surface and thus can be best appreciated using surface‐based fMRI analysis. We used this as a testbed for evaluating strategies that can reduce spatial blurring of fMRI data. Our results show that resolution losses can be mitigated at multiple points in preprocessing pathway. We show that unwanted blur is introduced at each step of volume transformation and surface projection, and can be ameliorated by replacing multi‐step transformations with equivalent single‐step transformations. Surprisingly, the simple approaches of volume upsampling and of cortical mesh refinement also helped to reduce resolution losses caused by interpolation. Volume upsampling also serves to improve motion estimation accuracy, which helps to reduce blur. Moreover, we demonstrate that the level of spatial blurring is nonuniform over the brain—knowledge which is critical for interpreting data in high‐resolution fMRI studies. Importantly, our study provides recommendations for reducing unwanted blurring during preprocessing as well as methods that enable quantitative comparisons between preprocessing strategies. These findings highlight several underappreciated sources of a spatial blur. Individually, the factors that contribute to spatial blur may appear to be minor, but in combination, the cumulative effects can hinder the interpretation of fine‐scale fMRI and the detectability of these fine‐scale features of functional architecture.

Photoacoustic imaging of squirrel monkey cortical and subcortical brain regions during peripheral electrical stimulation

The investigation of neuronal activity in non-human primate models is of critical importance due to their genetic similarity to human brains. In this study, we tested the feasibility of using photoacoustic imaging for the detection of cortical and subcortical responses due to peripheral electrical stimulation in a squirrel monkey model. Photoacoustic computed tomography and photoacoustic microscopy were applied on squirrel monkeys for real-time deep subcortical imaging and optical-resolution cortical imaging, respectively. The electrically evoked hemodynamic changes in primary somatosensory cortex, premotor cortices, primary motor cortex, and underlying subcortical areas were measured. Hemodynamic responses were observed in both cortical and subcortical brain areas at the cortices during external stimulation, demonstrating the feasibility of photoacoustic technique for functional imaging of non-human primate brain.

A thermal nociceptive patch in the S2 cortex of nonhuman primates: a combined functional magnetic resonance imaging and electrophysiology study

ABSTRACT

Human functional magnetic resonance imaging (fMRI) and behavioral studies have established the roles of cortical areas along the Sylvian fissure in sensing subjective pain. Yet, little is known about how sensory aspects of painful information are represented and processed by neurons in these regions and how their electrophysiological activities are related to fMRI signals. The current study aims to partially address this critical knowledge gap by performing fMRI-guided microelectrode mapping and recording studies in the homologous region of the parietal operculum in squirrel monkeys under light anesthesia. In each animal studied (n = 8), we detected mesoscale mini-networks for heat nociception in cortical regions around the lateral sulcus. Within the network, we discovered a ∼1.5 × 1.5-mm2-sized cortical patch that solely contained heat nociceptive neurons that aligned with the heat fMRI activation locus. These neurons responded slowly to thermal (heat and cold) nociceptive stimuli exclusively, continued firing for several seconds after the succession of stimulation, and exhibited multidigit receptive fields and high spontaneous firing rates. Similar to the fMRI responses, increasing temperatures in the nociceptive range led to a nonlinear increase in firing rates. The finding of a clustering of heat nociceptive neurons provides novel insights into the unique functional organization of thermal nociception in the S2 subregion of the primate brain. With fMRI, it supports the existence of a modality-preferred heat nociceptive patch that is spatially separated and intermingled with touch patches containing neurons with comparable receptive fields and the presence of functionally distinct mini-networks in primate opercular cortex.

Chemogenetic sensory fMRI reveals behaviorally relevant bidirectional changes in primate somatosensory network

Concurrent genetic neuromodulation and functional magnetic resonance imaging (fMRI) in primates has provided a valuable opportunity to assess the modified brain-wide operation in the resting state. However, its application to link the network operation with behavior still remains challenging. Here, we combined chemogenetic silencing of the primary somatosensory cortex (SI) with tactile fMRI and related behaviors in macaques. Focal chemogenetic silencing of functionally identified SI hand region impaired grasping behavior. The same silencing also attenuated hand stimulation-evoked fMRI signal at both the local silencing site and the anatomically and/or functionally connected downstream grasping network, suggesting altered network operation underlying the induced behavioral impairment. Furthermore, the hand region silencing unexpectedly disinhibited foot representation with accompanying behavioral hypersensitization. These results demonstrate that focal chemogenetic silencing with sensory fMRI in macaques unveils bidirectional network changes to generate multifaceted behavioral impairments, thereby opening a pivotal window toward elucidating the causal network operation underpinning higher brain functions in primates.

Common functional localizers to enhance NHP & cross-species neuroscience imaging research

Functional localizers are invaluable as they can help define regions of interest, provide cross-study comparisons, and most importantly, allow for the aggregation and meta-analyses of data across studies and laboratories. To achieve these goals within the non-human primate (NHP) imaging community, there is a pressing need for the use of standardized and validated localizers that can be readily implemented across different groups. The goal of this paper is to provide an overview of the value of localizer protocols to imaging research and we describe a number of commonly used or novel localizers within NHPs, and keys to implement them across studies. As has been shown with the aggregation of resting-state imaging data in the original PRIME-DE submissions, we believe that the field is ready to apply the same initiative for task-based functional localizers in NHP imaging. By coming together to collect large datasets across research group, implementing the same functional localizers, and sharing the localizers and data via PRIME-DE, it is now possible to fully test their robustness, selectivity and specificity. To do this, we reviewed a number of common localizers and we created a repository of well-established localizer that are easily accessible and implemented through the PRIME-RE platform.

Longitudinal fMRI measures of cortical reactivation and hand use with and without training after sensory loss in primates

In a series of previous studies, we demonstrated that damage to the dorsal column in the cervical spinal cord deactivates the contralateral somatosensory hand cortex and impairs hand use in a reach-to-grasp task in squirrel monkeys. Nevertheless, considerable cortical reactivation and behavioral recovery occurs over the following weeks to months after lesion. This timeframe may also be a window for targeted therapies to promote cortical reactivation and functional reorganization, aiding in the recovery process. Here we asked if and how task specific training of an impaired hand would improve behavioral recovery and cortical reorganization in predictable ways, and if recovery related cortical changes would be detectable using noninvasive functional magnetic resonance imaging (fMRI). We further asked if invasive neurophysiological mapping reflected fMRI results. A reach-to-grasp task was used to test impairment and recovery of hand use before and after dorsal column lesions (DC-lesion). The activation and organization of the affected primary somatosensory cortex (area 3b) was evaluated with two types of fMRI - either blood oxygenation level dependent (BOLD) or cerebral blood volume (CBV) with a contrast agent of monocrystalline iron oxide nanocolloid (MION) - before and after DC-lesion. At the end of the behavioral and fMRI studies, microelectrode recordings in the somatosensory areas 3a, 3b and 1 were used to characterize neuronal responses and verify the somatotopy of cortical reactivations. Our results indicate that even after nearly complete DC lesions, monkeys had both considerable post-lesion behavioral recovery, as well as cortical reactivation assessed with fMRI followed by extracellular recordings. Generalized linear regression analyses indicate that lesion extent is correlated with the behavioral outcome, as well as with the difference in the percent signal change from pre-lesion peak activation in fMRI. Monkeys showed behavioral recovery and nearly complete cortical reactivation by 9-12 weeks post-lesion (particularly when the DC-lesion was incomplete). Importantly, the specific training group revealed trends for earlier behavioral recovery and had higher magnitude of fMRI responses to digit stimulation by 5-8 weeks post-lesion. Specific kinematic measures of hand movements in the selected retrieval task predicted recovery time and related to lesion characteristics better than overall task performance success. For measures of cortical reactivation, we found that CBV scans provided stronger signals to vibrotactile digit stimulation as compared to BOLD scans, and thereby may be the preferred non-invasive way to study the cortical reactivation process after sensory deprivations from digits. When the reactivation of cortex for each of the digits was considered, the reactivation by digit 2 stimulation as measured with microelectrode maps and fMRI maps was best correlated with overall behavioral recovery.

Bidirectional and state-dependent modulation of brain activity by transcranial focused ultrasound in non-human primates

Transcranial focused ultrasound (FUS) stimulation under MRI guidance, coupled with functional MRI (fMRI) monitoring of effects, offers a precise, noninvasive technology to dissect functional brain circuits and to modulate altered brain functional networks in neurological and psychiatric disorders. Here we show that ultrasound at moderate intensities modulated neural activity bi-directionally. Concurrent sonication of somatosensory areas 3a/3b with 250 kHz FUS suppressed the fMRI signals produced there by peripheral tactile stimulation, while at the same time eliciting fMRI activation at inter-connected, off-target brain regions. Direct FUS stimulation of the cortex resulted in different degrees of BOLD signal changes across all five off-target regions, indicating that its modulatory effects on active and resting neurons differed. This is the first demonstration of the dual suppressive and excitative modulations of FUS on a specific functional circuit and of ability of concurrent FUS and MRI to evaluate causal interactions between functional circuits with neuron-class selectivity.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.

Longitudinal changes in DTI parameters of specific spinal white matter tracts correlate with behavior following spinal cord injury in monkeys

This study aims to evaluate how parameters derived from diffusion tensor imaging reflect axonal disruption and demyelination in specific white matter tracts within the spinal cord of squirrel monkeys following traumatic injuries, and their relationships to function and behavior. After a unilateral section of the dorsal white matter tract of the cervical spinal cord, we found that both lesioned dorsal and intact lateral tracts on the lesion side exhibited prominent disruptions in fiber orientation, integrity and myelination. The degrees of pathological changes were significantly more severe in segments below the lesion than above. The lateral tract on the opposite (non-injured) side was minimally affected by the injury. Over time, RD, FA, and AD values of the dorsal and lateral tracts on the injured side closely tracked measurements of the behavioral recovery. This unilateral section of the dorsal spinal tract provides a realistic model in which axonal disruption and demyelination occur together in the cord. Our data show that specific tract and segmental FA and RD values are sensitive to the effects of injury and reflect specific behavioral changes, indicating their potential as relevant indicators of recovery or for assessing treatment outcomes. These observations have translational value for guiding future studies of human subjects with spinal cord injuries.

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.

Investigating the Influence of Spatial Constraints on Ultimate Receive Coil Performance for Monkey Brain MRI at 7 T

The RF receive coil array has become increasingly vital in current MR imaging practice due to its extended spatial coverage, maintained high SNR, and improved capability of accelerating data acquisition. The performance of a coil array is intrinsically determined by the current patterns generated in coil elements as well as by the induced electromagnetic fields inside the object. Investigations of the ultimate performance constrained by a specific coil space, which defines all possible current patterns flowing within, offer the opportunity to evaluate coil-space parameters (i.e., coverage, coil-to-object distance, layer thickness, and coil element type) without the necessity of considering the realistic coil element geometry, coil elements layout, and number of receive channels in modeling. In this paper, to mimic 7-T monkey RF head coil design, seven hypothetical ultimate coil arrays with different coil-space configurations were mounted over a numerical macaque head model; by using Huygens's surface approximation method, the influences of coil-space design parameters were systematically investigated through evaluating the spatial constrained ultimate intrinsic SNR and ultimate g-factor. Moreover, simulations were also conducted by using four coil arrays with limited number of loop-only elements, in order to explore to what extent the ultimate coil performance can be achieved by using practical coil designs, and hence several guidelines in RF coil design for monkey brain imaging at 7 T have been tentatively concluded. It is believed that the present analysis will offer important implications in novel receive array design for monkey brain MR imaging at ultra-high field.

Neuromodulation of sensory networks in monkey brain by focused ultrasound with MRI guidance and detection

Focused ultrasound (FUS) has gained recognition as a technique for non-invasive neuromodulation with high spatial precision and the ability to both excite and inhibit neural activity. Here we demonstrate that MRI-guided FUS is capable of exciting precise targets within areas 3a/3b in the monkey brain, causing downstream activations in off-target somatosensory and associated brain regions which are simultaneously detected by functional MRI. The similarity between natural tactile stimulation-and FUS- evoked fMRI activation patterns suggests that FUS likely can excite populations of neurons and produce associated spiking activities that may be subsequently transmitted to other functionally related touch regions. The across-region differences in fMRI signal changes relative to area 3a/3b between tactile and FUS conditions also indicate that FUS modulated the tactile network differently. The significantly faster rising (>1 sec) fMRI signals elicited by direct FUS stimulation at the targeted cortical region suggest that a different neural hemodynamic coupling mechanism may be involved in generating fMRI signals. This is the first demonstration of imaging neural excitation effects of FUS with BOLD fMRI on a specific functional circuit in non-human primates.

Realistic models of apparent dynamic changes in resting-state connectivity in somatosensory cortex

Variations over time in resting-state correlations in blood oxygenation level-dependent (BOLD) signals from different cortical areas may indicate changes in brain functional connectivity. However, apparent variations over time may also arise from stationary signals when the sample duration is finite. Recently, a vector autoregressive (VAR) null model has been proposed to simulate real functional magnetic resonance imaging (fMRI) data, which provides a robust stationary model for identifying possible temporal dynamic changes in functional connectivity. In this work, we propose a simpler model that uses a filtered stationary dataset. The filtered stationary model generates statistically stationary time series from random data with a single prescribed correlation coefficient that is calculated as the average over the entire time series. In addition, we propose a dynamic model, which is better able to replicate real fMRI connectivity, estimated from monkey brain studies, than the two stationary models. We compare simulated results using these three models with the behavior of primary somatosensory cortex (S1) networks in anesthetized squirrel monkeys at high field (9.4 T), using a sliding window correlation analysis. We found that at short window sizes, both stationary models reproduced the distribution of correlations of real signals well, but at longer window sizes, a dynamic model reproduced the distribution of correlations of real signals better than the stationary models. While stationary models replicate several features of real data, a close representation of the behavior of resting-state data acquired from somatosensory cortex of non-human primates is obtained only when a dynamic correlation is introduced, suggesting dynamic variations in connectivity are real. Hum Brain Mapp 37:3897-3910, 2016. © 2016 Wiley Periodicals, Inc.

Cortical Merging in S1 as a Substrate for Tactile Input Grouping

Perception is a reconstruction process guided by rules based on knowledge about the world. Little is known about the neural implementation of the rules of object formation in the tactile sensory system. When two close tactile stimuli are delivered simultaneously on the skin, subjects feel a unique sensation, spatially centered between the two stimuli. Voltage-sensitive dye imaging (VSDi) and electrophysiological recordings [local field potentials (LFPs) and single units] were used to extract the cortical representation of two-point tactile stimuli in the primary somatosensory cortex of anesthetized Long-Evans rats. Although layer 4 LFP responses to brief costimulation of the distal region of two digits resembled the sum of individual responses, approximately one-third of single units demonstrated merging-compatible changes. In contrast to previous intrinsic optical imaging studies, VSD activations reflecting layer 2/3 activity were centered between the representations of the digits stimulated alone. This merging was found for every tested distance between the stimulated digits. We discuss this laminar difference as evidence that merging occurs through a buildup stream and depends on the superposition of inputs, which increases with successive stages of sensory processing. These findings show that layers 2/3 are involved in the grouping of sensory inputs. This process that could be inscribed in the cortical computing routine and network organization is likely to promote object formation and implement perception rules.

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.

High-resolution functional MRI identified distinct global intrinsic functional networks of nociceptive posterior insula and S2 regions in squirrel monkey brain

Numerous functional imaging and electrophysiological studies in humans and animals indicate that the two contiguous areas of secondary somatosensory cortex (S2) and posterior insula (pIns) are core regions in nociceptive processing and pain perception. In this study, we tested the hypothesis that the S2-pIns connection serves as a hub for connecting distinct sensory and affective nociceptive processing networks in the squirrel monkey brain. At 9.4T, we first mapped the brain regions that respond to nociceptive heat stimuli with high-resolution fMRI, and then used seed-based resting-state fMRI (rsfMRI) analysis to delineate and refine the global intrinsic functional connectivity circuits of the proximal S2 and pIns regions. In each subject, nociceptive (47.5°C) heat-evoked fMRI activations were detected in many brain regions, including primary somatosensory (S1), S2, pIns, area 7b, anterior cingulate cortex (ACC), primary motor cortex, prefrontal cortex, supplementary motor area, thalamus, and caudate. Using the heat-evoked fMRI activation foci in S2 and pIns as the seeds, voxel-wise whole-brain resting-state functional connectivity (rsFC) analysis revealed strong functional connections between contralateral S2 and pIns, as well as their corresponding regions in the ipsilateral hemisphere. Spatial similarity and overlap analysis identified each region as part of two distinct intrinsic functional networks with 7% overlap: sensory S2-S1-area 7b and affective pIns-ACC-PCC networks. Moreover, a high degree of overlap was observed between the combined rsFC maps of nociceptive S2 and pIns regions and the nociceptive heat-evoked activation map. In summary, our study provides evidence for the existence of two distinct intrinsic functional networks for S2 and pIns nociceptive regions, and these two networks are joined via the S2-pIns connection. Brain regions that are involved in processing nociceptive inputs are also highly interconnected at rest. The presence of robust and distinct S1-S2-area 7b and pIns-ACC-PCC rsFC networks under anesthesia underscores their fundamental roles in processing nociceptive information.

Biophysical and neural basis of resting state functional connectivity: Evidence from non-human primates

Functional MRI (fMRI) has evolved from simple observations of regional changes in MRI signals caused by cortical activity induced by a task or stimulus, to task-free acquisitions of images in a resting state. Such resting state signals contain low frequency fluctuations which may be correlated between voxels, and strongly correlated regions are deemed to reflect functional connectivity within synchronized circuits. Resting state functional connectivity (rsFC) measures have been widely adopted by the neuroscience community, and are being used and interpreted as indicators of intrinsic neural circuits and their functional states in a broad range of applications, both basic and clinical. However, there has been relatively little work reported that validates whether inter-regional correlations in resting state fluctuations of fMRI (rsfMRI) signals actually measure functional connectivity between brain regions, or to establish how MRI data correlate with other metrics of functional connectivity. In this mini-review, we summarize recent studies of rsFC within mesoscopic scale cortical networks (100μm-10mm) within a well defined functional region of primary somatosensory cortex (S1), as well as spinal cord and brain white matter in non-human primates, in which we have measured spatial patterns of resting state correlations and validated their interpretation with electrophysiological signals and anatomic connections. Moreover, we emphasize that low frequency correlations are a general feature of neural systems, as evidenced by their presence in the spinal cord as well as white matter. These studies demonstrate the valuable role of high field MRI and invasive measurements in an animal model to inform the interpretation of human imaging studies.

Positive and Negative Somatotopic BOLD Responses in Contralateral Versus Ipsilateral Penfield Homunculus

One of the basic properties of sensory cortices is their topographical organization. Most imaging studies explored this organization using the positive blood oxygenation level-dependent (BOLD) signal. Here, we studied the topographical organization of both positive and negative BOLD in contralateral and ipsilateral primary somatosensory cortex (S1). Using phase-locking mapping methods, we verified the topographical organization of contralateral S1, and further showed that different body segments elicit pronounced negative BOLD responses in both hemispheres. In the contralateral hemisphere, we found a sharpening mechanism in which stimulation of a given body segment triggered a gradient of activation with a significant deactivation in more remote areas. In the ipsilateral cortex, deactivation was not only located in the homolog area of the stimulated parts but rather was widespread across many parts of S1. Additionally, analysis of resting-state functional magnetic resonance imaging signal showed a gradient of connectivity to the neighboring contralateral body parts as well as to the ipsilateral homologous area for each body part. Taken together, our results indicate a complex pattern of baseline and activity-dependent responses in the contralateral and ipsilateral sides. Both primary sensory areas were characterized by unique negative BOLD responses, suggesting that they are an important component in topographic organization of sensory cortices.

Effects of isoflurane anesthesia on resting-state fMRI signals and functional connectivity within primary somatosensory cortex of monkeys

INTRODUCTION

Correlated low-frequency fluctuations of resting-state functional magnetic resonance imaging (rsfMRI) signals have been widely used for inferring intrinsic brain functional connectivity (FC). In animal studies, accurate estimate of anesthetic effects on rsfMRI signals is demanded for reliable interpretations of FC changes. We have previously shown that inter-regional FC can reliably delineate local millimeter-scale circuits within digit representations of primary somatosensory cortex (S1) subregions (areas 3a, 3b, and 1) in monkeys under isoflurane anesthesia. The goals of this study are to determine (1) the general effects of isoflurane on rsfMRI signals in the S1 circuit and (2) whether the effects are functional- and regional- dependent, by quantifying the relationships between isoflurane levels, power and inter-regional correlation coefficients in digit and face regions of distinct S1 subregions.

METHODS

Functional MRI data were collected from male adult squirrel monkeys at three different isoflurane levels (1.25%, 0.875%, and 0.5%). All scans were acquired on a 9.4T magnet with a 3-cm-diameter surface transmit-receive coil centered over the S1 cortex. Power and seed-based inter-regional functional connectivity analyses were subsequently performed.

RESULTS

As anesthesia level increased, we observed (1) diminishing amplitudes of signal fluctuations, (2) reduced power of fluctuations in the low-frequency band used for connectivity measurements, (3) decreased inter-voxel connectivity around seed regions, and (4) weakened inter-regional FC across all pairs of regions of interest (digit-to-digit). The low-frequency power measures derived from rsfMRI signals from control muscle regions, however, did not exhibit any isoflurane level-related changes. Within the isoflurane dosage range we tested, the inter-regional functional connectivity differences were still detectable, and the effects of isoflurane did not differ across region-of-interest (ROI) pairs.

CONCLUSION

Our data demonstrate that isoflurane induced similar dose-dependent suppressive effects on the power of rsfMRI signals and local fine-scale FC across functionally related but distinct S1 subregions.

Altered Spatiotemporal Dynamics of Cortical Activation to Tactile Stimuli in Somatosensory Area 3b and Area 1 of Monkeys after Spinal Cord Injury

Reactivation of deafferented cortex plays a key role in mediating the recovery of lost functions, although the precise mechanism is not fully understood. This study simultaneously characterized the dynamic spatiotemporal features of tactile responses in areas 3b and 1 before and 6–8 weeks after partial dorsal column lesion (DCL), and examined how the reactivation relates to the recovery of simple hand use in squirrel monkeys. A combination of high spatiotemporal resolution functional intrinsic optical imaging, microelectrode mapping, behavioral assessment, and tracer histology methods were used. Compared with the normal cortex, we found that the responses of deafferented areas 3b and 1 to 3 s of continuous 8 Hz tactile stimulation of a single digit were significantly weaker and more transient. This finding indicates a loss of response to sustained tactile stimuli. The activation area enlarged for areas 3b and 1 in both directions along digit representation (medial–lateral) and across areas (anterior–posterior). All subjects showed behavioral deficits in a food reaching-grasping-retrieving task within the first 5 weeks after DCL, but recovered at the time when optical images were acquired. Summarily, we showed that these populations of cortical neurons responded to peripheral tactile inputs, albeit in significantly altered manners in each area, several weeks after deafferentation. We propose that compromised ascending driven inputs, impaired lateral inhibition, and local integration of input signals may account for the altered spatiotemporal dynamics of the reactivated areas 3b and 1 cortices. Further investigation with large sample sizes is needed to fully characterize the effects of deafferentation on area 1 activation size.

Spatiotemporal trajectories of reactivation of somatosensory cortex by direct and secondary pathways after dorsal column lesions in squirrel monkeys

After lesions of the somatosensory dorsal column (DC) pathway, the cortical hand representation can become unresponsive to tactile stimuli, but considerable responsiveness returns over weeks of post-lesion recovery. The reactivation suggests that preserved subthreshold sensory inputs become potentiated and axon sprouting occurs over time to mediate recovery. Here, we studied the recovery process in 3 squirrel monkeys, using high-resolution cerebral blood volume-based functional magnetic resonance imaging (CBV-fMRI) mapping of contralateral somatosensory cortex responsiveness to stimulation of distal finger pads with low and high level electrocutaneous stimulation (ES) before and 2, 4, and 6weeks after a mid-cervical level contralateral DC lesion. Both low and high intensity ES of digits revealed the expected somatotopy of the area 3b hand representation in pre-lesion monkeys, while in areas 1 and 3a, high intensity stimulation was more effective in activating somatotopic patterns. Six weeks post-lesion, and irrespective of the severity of loss of direct DC inputs (98%, 79%, 40%), somatosensory cortical area 3b of all three animals showed near complete recovery in terms of somatotopy and responsiveness to low and high intensity ES. However there was significant variability in the patterns and amplitudes of reactivation of individual digit territories within and between animals, reflecting differences in the degree of permanent and/or transient silencing of primary DC and secondary inputs 2weeks post-lesion, and their spatio-temporal trajectories of recovery between 2 and 6weeks. Similar variations in the silencing and recovery of somatotopy and responsiveness to high intensity ES in areas 3a and 1 are consistent with individual differences in damage to and recovery of DC and spinocuneate pathways, and possibly the potentiation of spinothalamic pathways. Thus, cortical deactivation and subsequent reactivation depends not only on the degree of DC lesion, but also on the severity and duration of loss of secondary as well as primary inputs revealed by low and high intensity ES.

Effects of anesthesia on resting state BOLD signals in white matter of non-human primates

Resting state functional magnetic resonance imaging (rsfMRI) has been widely used to measure functional connectivity between cortical regions of the brain. However, there have been minimal reports of bold oxygenation level dependent (BOLD) signals in white matter, and even fewer attempts to detect resting state connectivity. Recently, there has been growing evidence that suggests that reliable detection of white matter BOLD signals may be possible. We have previously shown that nearest neighbor inter-voxel correlations of resting state BOLD signal fluctuations in white matter are anisotropic and can be represented by a functional correlation tensor, but the biophysical origins of these signal variations are not clear. We aimed to assess whether MRI signal fluctuations in white matter vary for different baseline levels of neural activity. We performed imaging studies on live squirrel monkeys under different levels of isoflurane anesthesia at 9.4T. We found 1) the fractional power (0.01-0.08Hz) in white matter was between 60 to 75% of the level in gray matter; 2) the power in both gray and white matter low frequencies decreased monotonically in similar manner with increasing levels of anesthesia; 3) the distribution of fractional anisotropy values of the functional tensors in white matter were significantly higher than those in gray matter; and 4) the functional tensor eigenvalues decreased with increasing level of anesthesia. Our results suggest that as anesthesia level changes baseline neural activity, white matter signal fluctuations behave similarly to those in gray matter, and functional tensors in white matter are affected in parallel.

fMRI at High Spatial Resolution: Implications for BOLD-Models

As high-resolution functional magnetic resonance imaging (fMRI) and fMRI of cortical layers become more widely used, the question how well high-resolution fMRI signals reflect the underlying neural processing, and how to interpret laminar fMRI data becomes more and more relevant. High-resolution fMRI has shown laminar differences in cerebral blood flow (CBF), volume (CBV), and neurovascular coupling. Features and processes that were previously lumped into a single voxel become spatially distinct at high resolution. These features can be vascular compartments such as veins, arteries, and capillaries, or cortical layers and columns, which can have differences in metabolism. Mesoscopic models of the blood oxygenation level dependent (BOLD) response therefore need to be expanded, for instance, to incorporate laminar differences in the coupling between neural activity, metabolism and the hemodynamic response. Here we discuss biological and methodological factors that affect the modeling and interpretation of high-resolution fMRI data. We also illustrate with examples from neuropharmacology and the negative BOLD response how combining BOLD with CBF- and CBV-based fMRI methods can provide additional information about neurovascular coupling, and can aid modeling and interpretation of high-resolution fMRI.

In Vivo Mapping of Cortical Columnar Networks in the Monkey with Focal Electrical and Optical Stimulation

There are currently largescale efforts to understand the brain as a connection machine. However, there has been little emphasis on understanding connection patterns between functionally specific cortical columns. Here, we review development and application of focal electrical and optical stimulation methods combined with optical imaging and fMRI mapping in the non-human primate. These new approaches, when applied systematically on a large scale, will elucidate functionally specific intra-areal and inter-areal network connection patterns. Such functionally specific network data can provide accurate views of brain network topology.

Differential fMRI Activation Patterns to Noxious Heat and Tactile Stimuli in the Primate Spinal Cord

Mesoscale local functional organizations of the primate spinal cord are largely unknown. Using high-resolution fMRI at 9.4 T, we identified distinct interhorn and intersegment fMRI activation patterns to tactile versus nociceptive heat stimulation of digits in lightly anesthetized monkeys. Within a spinal segment, 8 Hz vibrotactile stimuli elicited predominantly fMRI activations in the middle part of ipsilateral dorsal horn (iDH), along with significantly weaker activations in ipsilateral (iVH) and contralateral (cVH) ventral horns. In contrast, nociceptive heat stimuli evoked widespread strong activations in the superficial part of iDH, as well as in iVH and contralateral dorsal (cDH) horns. As controls, only weak signal fluctuations were detected in the white matter. The iDH responded most strongly to both tactile and heat stimuli, whereas the cVH and cDH responded selectively to tactile versus nociceptive heat, respectively. Across spinal segments, iDH activations were detected in three consecutive segments in both tactile and heat conditions. Heat responses, however, were more extensive along the cord, with strong activations in iVH and cDH in two consecutive segments. Subsequent subunit B of cholera toxin tracer histology confirmed that the spinal segments showing fMRI activations indeed received afferent inputs from the stimulated digits. Comparisons of the fMRI signal time courses in early somatosensory area 3b and iDH revealed very similar hemodynamic stimulus-response functions. In summary, we identified with fMRI distinct segmental networks for the processing of tactile and nociceptive heat stimuli in the cervical spinal cord of nonhuman primates. Significance statement: This is the first fMRI demonstration of distinct intrasegmental and intersegmental nociceptive heat and touch processing circuits in the spinal cord of nonhuman primates. This study provides novel insights into the local functional organizations of the primate spinal cord for pain and touch, information that will be valuable for designing and optimizing therapeutic interventions for chronic pain management.

The postcentral sulcal complex and the transverse postcentral sulcus and their relation to sensorimotor functional organization

It has been demonstrated that the postcentral sulcus, which forms the posterior boundary of the sensorimotor region, is a complex of distinct sulcal segments. Although the general somatotopic arrangement in the human sensorimotor cortex is relatively well known, we do not know whether the different segments of the postcentral sulcus relate in a systematic way to the sensorimotor functional representations. Participants were scanned with functional magnetic resonance imaging while they made movements of different body parts and the location of functional activity was examined on a subject-by-subject basis with respect to the morphological features of the postcentral sulcus. The findings demonstrate that the postcentral sulcus of each subject may be divided into five segments and there is a tight relationship between sensorimotor representations of different body parts and specific segments of the postcentral sulcus. The results also addressed the issue of the transverse postcentral sulcus, a short sulcus that is present within the ventral part of the postcentral gyrus in some brains. It was shown that, when present, this sulcus is functionally related to the oral (mouth and tongue) sensorimotor representation. When this sulcus is not present, the inferior postcentral sulcus which is also related to the oral representation is longer. Thus, the sulcal morphology provides an improved framework for functional assignments in individual subjects.

Injury alters intrinsic functional connectivity within the primate spinal cord

Recent demonstrations of correlated low-frequency MRI signal variations between subregions of the spinal cord at rest in humans, similar to those found in the brain, suggest that such resting-state functional connectivity constitutes a common feature of the intrinsic organization of the entire central nervous system. We report our detection of functional connectivity within the spinal cords of anesthetized squirrel monkeys at rest and show that the strength of connectivity within these networks is altered by the effects of injuries. By quantifying the low-frequency MRI signal correlations between different horns within spinal cord gray matter, we found distinct functional connectivity relationships between the different sensory and motor horns, a pattern that was similar to activation patterns evoked by nociceptive heat or tactile stimulation of digits. All horns within a single spinal segment were functionally connected, with the strongest connectivity occurring between ipsilateral dorsal and ventral horns. Each horn was strongly connected to the same horn on neighboring segments, but this connectivity reduced drastically along the spinal cord. Unilateral injury to the spinal cord significantly weakened the strength of the intrasegment horn-to-horn connectivity only on the injury side and in slices below the lesion. These findings suggest resting-state functional connectivity may be a useful biomarker of functional integrity in injured and recovering spinal cords.

Brain process for perception of the "out of the body" tactile illusion for virtual object interaction

"Out of the body" tactile illusion refers to the phenomenon in which one can perceive tactility as if emanating from a location external to the body without any stimulator present there. Taking advantage of such a tactile illusion is one way to provide and realize richer interaction feedback without employing and placing actuators directly at all stimulation target points. However, to further explore its potential, it is important to better understand the underlying physiological and neural mechanism. As such, we measured the brain wave patterns during such tactile illusion and mapped out the corresponding brain activation areas. Participants were given stimulations at different levels with the intention to create veridical (i.e., non-illusory) and phantom sensations at different locations along an external hand-held virtual ruler. The experimental data and analysis indicate that both veridical and illusory sensations involve, among others, the parietal lobe, one of the most important components in the tactile information pathway. In addition, we found that as for the illusory sensation, there is an additional processing resulting in the delay for the ERP (event-related potential) and involvement by the limbic lobe. These point to regarding illusion as a memory and recognition task as a possible explanation. The present study demonstrated some basic understanding; how humans process "virtual" objects and the way associated tactile illusion is generated will be valuable for HCI (Human-Computer Interaction).

Parallel functional reorganizations of somatosensory areas 3b and 1, and S2 following spinal cord injury in squirrel monkeys

Multiple somatosensory cortices of adult primates reorganize following spinal cord injury, but little is known about the temporal dynamics and inter-areal differences of the reorganization. Using longitudinal high-resolution fMRI in combination with microelectrode recordings and tracer histology, we previously illustrated a two-phase dynamic spatial reorganization of digit representations in area 3b within weeks after a unilateral lesion of the dorsal column in squirrel monkeys (Chen et al., 2012). Here we report that higher-order area 1 and secondary somatosensory cortex (S2) underwent similar spatial reorganizations, which were characterized by shifted and expanded digit activations at week 4 after lesion, which then shifted back and contracted by week 8. In addition, the responsiveness of areas 3b and 1, and S2, as measured by the magnitude of the BOLD signal change to tactile stimuli, was reduced markedly at 4 weeks and then recovered to ~50% of the prelesion level at 8 weeks, a time when behavioral recovery was complete, as assessed by successful food retrieval rates. Across animals, the extents of spatial reorganizations and changes in cortical responsiveness and activation sizes in all three areas were correlated with the degree of afferent disruption. In summary, our data show that more severe afferent disruption was associated with greater cortical plasticity and behavioral impairment. Reorganization that occurred in area 3b, area 1, and S2 were similar across most measures.

Robust data driven model order estimation for independent component analysis of FMRI data with low contrast to noise

Independent component analysis (ICA) has been successfully utilized for analysis of functional MRI (fMRI) data for task related as well as resting state studies. Although it holds the promise of becoming an unbiased data-driven analysis technique, a few choices have to be made prior to performing ICA, selection of a method for determining the number of independent components (nIC) being one of them. Choice of nIC has been shown to influence the ICA maps, and various approaches (mostly relying on information theoretic criteria) have been proposed and implemented in commonly used ICA analysis packages, such as MELODIC and GIFT. However, there has been no consensus on the optimal method for nIC selection, and many studies utilize arbitrarily chosen values for nIC. Accurate and reliable determination of true nIC is especially important in the setting where the signals of interest contribute only a small fraction of the total variance, i.e. very low contrast-to-noise ratio (CNR), and/or very focal response. In this study, we evaluate the performance of different model order selection criteria and demonstrate that the model order selected based upon bootstrap stability of principal components yields more reliable and accurate estimates of model order. We then demonstrate the utility of this fully data-driven approach to detect weak and focal stimulus-driven responses in real data. Finally, we compare the performance of different multi-run ICA approaches using pseudo-real data.

Distinct fine-scale fMRI activation patterns of contra- and ipsilateral somatosensory areas 3b and 1 in humans

Inter-areal and ipsilateral cortical responses to tactile stimulation have not been well described in human S1 cortex. By taking advantage of the high signal-to-noise ratio at 7 T, we quantified blood oxygenation level dependent (BOLD) response patterns and time courses to tactile stimuli on individual distal finger pads at a fine spatial scale, and examined whether there are inter-areal (area 3b versus area 1) and interhemispheric response differences to unilateral tactile stimulation in healthy human subjects. We found that 2-Hz tactile stimulation of individual fingertips evoked detectable BOLD signal changes in both contralateral and ipsilateral area 3b and area 1. Contralateral digit activations were organized in an orderly somatotopic manner, and BOLD responses in area 3b were more digit selective than those in area 1. However, the area of cortex that was responsive to stimulation of a single digit (stimulus-response field) was similar across areas. In the ipsilateral hemisphere, response magnitudes in both areas 3b and 1 were significantly weaker than those of the contralateral hemisphere. Digit activations exhibited no clear somatotopic organizational pattern in either area 3b or area 1, yet digit selectivity was retained in area 1 but not in area 3b. The observation of distinct digit-selective responses of contralateral area 3b versus area 1 supports a higher order function of contralateral area 1 in spatial integration. In contrast, ipsilateral cortices may play a less discriminative role in the perception of unilateral tactile sensation in humans.

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.

The relationship of anatomical and functional connectivity to resting-state connectivity in primate somatosensory cortex

Studies of resting-state activity in the brain have provoked critical questions about the brain's functional organization, but the biological basis of this activity is not clear. Specifically, the relationships between interregional correlations in resting-state measures of activity, neuronal functional connectivity and anatomical connectivity are much debated. To investigate these relationships, we have examined both anatomical and steady-state functional connectivity within the hand representation of primary somatosensory cortex (areas 3b and 1) in anesthetized squirrel monkeys. The comparison of three data sets (fMRI, electrophysiological, and anatomical) indicate two primary axes of information flow within the SI: prominent interdigit interactions within area 3b and predominantly homotopic interactions between area 3b and area 1. These data support a strikingly close relationship between baseline functional connectivity and anatomical connections. This study extends findings derived from large-scale cortical networks to the realm of local millimeter-scale networks.

Large-scale reorganization of the somatosensory cortex of adult macaque monkeys revealed by fMRI

Somatosensory cortex of adult primates undergoes topographic reorganization following spinal cord or peripheral nerve injuries. Electrophysiological studies in monkeys show that after chronic lesions of dorsal columns of the spinal cord at cervical levels, there is an expansion of face representation into the deafferented hand region of area 3b of cortex. However, these techniques can sample only a limited portion of the brain. In order to help understand mechanisms of brain reorganization use of noninvasive tools in non-human primate experimental model is important. Use of blood oxygen level dependent-functional magnetic resonance imaging (BOLD-fMRI) to study brain reorganization in non-human primates has been extremely limited. Here, we show that in monkeys with long-term unilateral lesions of the dorsal columns at cervical levels, tactile stimulation of the chin showed BOLD activation in the deafferented hand region of contralesional area 3b in the post-central gyrus. In a monkey with a partial lesion of the dorsal columns, stimulations of both hand and chin activated the partially deafferented hand region. We also show that the somatotopic organization in the non-deafferented ipsilesional somatosensory cortex remained normal.

Illusory and veridical mapping of tactile objects in the primary somatosensory and posterior parietal cortex

While several behavioral and neuroscience studies have explored visual, auditory, and cross-modal illusions, information about the phenomenology and neural correlates of somatosensory illusions is meager. By combining psychophysics and somatosensory evoked potentials, we explored in healthy humans the neural correlates of 2 compelling tactuo-proprioceptive illusions, namely Aristotle (1 object touching the contact area between 2 crossed fingers is perceived as 2 lateral objects) and Reverse illusions (2 lateral objects are perceived as 1 between crossed-fingers object). These illusions likely occur because of the tactuo-proprioceptive conflict induced by fingers being crossed in a non-natural posture. We found that different regions in the somatosensory stream exhibit different proneness to the illusions. Early electroencephalographic somatosensory activity (at 20 ms) originating in the primary somatosensory cortex (S1) reflects the phenomenal rather than the physical properties of the stimuli. Notably, later activity (around 200 ms) originating in the posterior parietal cortex is higher when subjects resist the illusions. Thus, while S1 activity is related to illusory perception, PPC acts as a conflict resolver that recodes tactile events from somatotopic to spatiotopic frames of reference and ultimately enables veridical perception.

Aristotle's illusion reveals interdigit functional somatosensory alterations in focal hand dystonia

In focal hand dystonia, the cortical somatosensory representation of the fingers is abnormal, with overlapping receptive fields and reduced interdigit separation. These abnormalities are associated with deficits in sensory perception, as previously demonstrated by applying tactile stimuli to one finger at a time. What is still unknown is whether the sensory deficits can be observed when tactile perception involves more than one finger. To address this issue, we applied 'Aristotle's illusion' to 15 patients with focal hand dystonia, 15 patients with dystonia not affecting the hand (blepharospasm and cervical dystonia) and 15 healthy control subjects. In this illusion, one object touching the contact point of two crossed fingertips is perceived as two objects by a blindfolded subject. The same object placed between two parallel fingertips is correctly perceived as one. The illusory doubling sensation is because of the fact that the contact point between the crossed fingers consists of non-adjacent and functionally unrelated skin regions, which usually send sensory signals to separate spots in the somatosensory cortex. In our study, participants were touched by one sphere between the second-third digits, the second-fourth digits and the fourth-fifth digits of both hands, either in crossed or in parallel position, and had to refer whether they felt one or two stimuli. The percentage of 'two stimuli' responses was an index of the illusory doubling. Both healthy control subjects and dystonic patients presented Aristotle's illusion when the fingers were crossed. However, patients with focal hand dystonia presented a significant reduction of the illusion when the sphere was placed between the crossed fourth and fifth digits of the affected hand. This reduction correlated with the severity of motor disease at the fingers. Similar findings were not observed in non-hand dystonia and control groups. The reduction of Aristotle's illusion in non-affected fingers and its preservation in affected fingers suggests dissociation between the abnormal processing of sensory signals and the motor impairment. Based on previous evidence showing that the sensory signals coming from the fourth digit determine lower activation in the somatosensory cortex than those coming from the fifth digit, we suggest that in the crossed position, the tactile information conveyed by the fifth digit prevailed over the fourth digit, thus resulting in the perception of one stimulus. The reduction of the illusory doubling perception, therefore, may represent the functional correlate of the different level of activation between the fourth and the fifth digit in the somatosensory cortex.

Optical imaging of visually guided reaching in macaque posterior parietal cortex

Sensorimotor transformation for reaching movements in primates requires a large network of visual, parietal, and frontal cortical areas. We performed intrinsic optical imaging over posterior parietal cortex including areas 7a and the dorsal perilunate in macaque monkeys during visually guided hand movements. Reaching was performed while foveating one of nine static reach targets; thus eye-position-varied concurrently with reach position. The hemodynamic reflectance signal was analyzed during specific phases of the task including pre-reach, reach, and touch epochs. The eye position maps changed substantially as the task progressed: First, direction of spatial tuning shifted from a weak preference close to the center to the lower eye positions in both cortical areas. Overall tuning strength was greater in area 7a. Second, strength of spatial tuning increased from the early pre-reach to the later touch epoch. These consistent temporal changes suggest that dynamic properties of the reflectance signal were modulated by task parameters. The peak amplitude and peak delay of the reflectance signal showed considerable differences between eye position but were similar between areas. Compared with a detection task using a lever response, the reach task yielded higher amplitudes and longer delays. These findings demonstrate a spatially tuned topographical representation for reaching in both areas and suggest a strong synergistic combination of various feedback signals that result in a spatially tuned amplification of the hemodynamic response in posterior parietal cortex.

The functional architecture of S1 during touch observation described with 7 T fMRI

Recent studies indicate that the primary somatosensory cortex (S1) is active not only when touch is physically perceived but also when it is merely observed to be experienced by another person. This social responsivity of S1 has important implications for our understanding of S1 functioning. However, S1 activity during touch observation has not been characterized in great detail to date. We focused on two features of the S1 functional architecture during touch observation, namely the topographical arrangement of index and middle finger receptive fields (RFs), and their dynamic shrinkage during concurrent activation. Both features have important implications for human behavior. We conducted two fMRI studies at 7 T, one where touch was physically perceived, and one where touch was observed. In the two experiments, participants either had their index finger and/or middle finger stimulated using paintbrushes, or just observed similar touch events on video. Our data show that observing and physically experiencing touch elicits overlapping activity changes in S1. In addition, observing touch to the index finger or the middle finger alone evoked topographically arranged activation foci in S1. Importantly, when co-activated, the index and middle finger RFs not only shrank during physical touch perception, but also during touch observation. Our data, therefore, indicate a similarity between the functional architecture of S1 during touch observation and physical touch perception with respect to single-digit topography and RF shrinkage. These results may allow the tentative conclusion that even primary somatosensory experiences, such as physical touch perception, can be shared amongst individuals.

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.

Columnar distribution of activity dependent gabaergic depolarization in sensorimotor cortical neurons

Background

GABA, the major inhibitory neurotransmitter in CNS, has been demonstrated to paradoxically produce excitation even in mature brain. However activity-dependent form of GABA excitation in cortical neurons has not been observed. Here we report that after an intense electrical stimulation adult cortical neurons displayed a transient GABA excitation that lasted for about 30s.

Results

Whole-cell patch recordings were performed to evaluate the effects of briefly applied GABA on pyramidal neurons in adult rodent sensorimotor cortical slice before and after 1 s, 20 Hz suprathreshold electrical stimulation of the junction between layer 6 and the underlying white matter (L6/WM stimulation). Immediately after L6/WM stimulation, GABA puffs produced neuronal depolarization in the center of the column-shaped region. However, both prior to or 30s after stimulation GABA puffs produced hyperpolarization of neurons. 2-photon imaging in neurons infected with adenovirus carrying a chloride sensor Clomeleon revealed that GABA induced depolarization is due to an increase in [Cl^-]_i after stimulation. To reveal the spatial extent of excitatory action of GABA, isoguvacine, a GABA_A receptors agonist, was applied right after stimulation while monitoring the intracellular Ca²⁺ concentration in pyramidal neurons. Isoguvacine induced an increase in [Ca²⁺]_i in pyramidal neurons especially in the center of the column but not in the peripheral regions of the column. The global pattern of the Ca²⁺ signal showed a column-shaped distribution along the stimulation site.

Conclusion

These results demonstrate that the well-known inhibitory transmitter GABA rapidly switches from hyperpolarization to depolarization upon synaptic activity in adult somatosensory cortical neurons.

Spatiotemporal characteristics and vascular sources of neural-specific and -nonspecific fMRI signals at submillimeter columnar resolution

The neural specificity of hemodynamic-based functional magnetic resonance imaging (fMRI) signals is dependent on both the vascular regulation and the sensitivity of the applied fMRI technique to different types and sizes of blood vessels. In order to examine the specificity of MRI-detectable hemodynamic responses, submillimeter blood oxygenation level-dependent (BOLD) and cerebral blood volume (CBV) fMRI studies were performed in a well-established cat orientation column model at 9.4 T. Neural-nonspecific and -specific signals were separated by comparing the fMRI responses of orthogonal orientation stimuli. The BOLD response was dominantly neural-nonspecific, mostly originating from pial and intracortical emerging veins, and thus was highly correlated with baseline blood volume. Uneven baseline CBV may displace or distort small functional domains in high-resolution BOLD maps. The CBV response in the parenchyma exhibited dual spatiotemporal characteristics, a fast and early neural-nonspecific response (with 4.3-s time constant) and a slightly slower and delayed neural-specific response (with 9.4-s time constant). The nonspecific CBV signal originates from early-responding arteries and arterioles, while the specific CBV response, which is not correlated with baseline blood volume, arises from late-responding microvessels including small pre-capillary arterioles and capillaries. Our data indicate that although the neural specificity of CBV fMRI signals is dependent on stimulation duration, high-resolution functional maps can be obtained from steady-state CBV studies.

Layer-specific BOLD activation in awake monkey V1 revealed by ultra-high spatial resolution functional magnetic resonance imaging

The laminar structure of the cortex has previously been explored both in non-human primates and human subjects using high-resolution functional magnetic resonance imaging (fMRI). However, whether the spatial specificity of the blood-oxygenation-level-dependent (BOLD) fMRI is sufficiently high to reveal lamina specific organization in the cortex reliably is still unclear. In this study we demonstrate for the first time the detection of such layer-specific activation in awake monkeys at the spatial resolution of 200 × 200 × 1000 μm(3) in a vertical 4.7 T scanner. Results collected in trained monkeys are high in contrast-to-noise ratio and low in motion artifacts. Isolation of laminar activation was aided by choosing the optimal slice orientation and thickness using a novel pial vein pattern analysis derived from optical imaging. We found that the percent change of GE-BOLD signal is the highest at a depth corresponding to layer IV. Changes in the middle layers (layer IV) were 30% greater than changes in the top layers (layers I-III), and 32% greater than the bottom layers (layers V/VI). The laminar distribution of BOLD signal correlates well with neural activity reported in the literature. Our results suggest that the high intrinsic spatial resolution of GE-BOLD signal is sufficient for mapping sub-millimeter functional structures in awake monkeys. This degree of spatial specificity will be useful for mapping both laminar activations and columnar structures in the cerebral cortex.

Columnar specificity of microvascular oxygenation and blood flow response in primary visual cortex: evaluation by local field potential and spiking activity

The relation of cortical microcirculation, oxygen metabolism, and underlying neuronal network activity remains poorly understood. Anatomical distribution of cortical microvasculature and its relationship to cortical functional domains suggests that functional organizations may be revealed by mapping cerebral blood flow responses. However, there is little direct experimental evidence and a lack of electrophysiological evaluation. In this study, we mapped ocular-dominance columns in primary visual cortex (V1) of anesthetized macaques with capillary flow-based laser speckle contrast imaging and deoxyhemoglobin-based intrinsic optical imaging. In parallel, the local field potentials (LFPs) and spikes were recorded from a linear array of eight microelectrodes, carefully positioned into left and right eye columns in V1. We found differential activation maps of blood flow, after masking large superficial draining vessels, exhibited a column-like pattern similar as the oximetric maps. Both the activated spikes and γ-band LFP demonstrated corresponding eye preference, consistent with the imaging maps. Our results present direct support in favor of previous proposals that the regulation of microcirculation can be as fine as the submillimeter scale, suggesting that cortical vasculature is functionally organized at the columnar level in a manner appropriate for supplying energy demands of functionally specific neuronal populations.

A novel device for the study of somatosensory information processing

Current methods for applying multi-site vibratory stimuli to the skin typically involve the use of multiple, individual vibrotactile stimulators. Limitations of such an arrangement include difficulty with both positioning the stimuli as well as ensuring that stimuli are delivered in a synchronized and deliberate manner. Previously, we reported a two-site tactile stimulator that was developed in order to solve these problems (Tannan et al., 2007a). Due to both the success of that novel stimulator and the limitations that were inherent in that device, we designed and fabricated a four-site stimulator that provides a number of advantages over the previous version. First, the device can stimulate four independent skin sites and is primarily designed for stimulating the digit tips. Second, the positioning of the probe tips has been re-designed to provide better ergonomic hand placement. Third, the device is much more portable than the previously reported stimulator. Fourth, the stimulator head has a much smaller footprint on the table or surface where it resides. To demonstrate the capacity of the device for delivering tactile stimulation at four independent sites, a finger agnosia protocol, in the presence and absence of conditioning stimuli, was conducted on seventeen healthy control subjects. The study demonstrated that with increasing amplitudes of vibrotactile conditioning stimuli concurrent with the agnosia test, inaccuracies of digit identification increased, particularly at digits D3 and D4. The results are consistent with prior studies that implicated synchronization of adjacent and near-adjacent cortical ensembles with conditioning stimuli in impacting TOJ performance (Tommerdahl et al., 2007a,b).

The postcentral gyrus shows sustained fMRI activation during the tactile motion aftereffect

The tactile motion aftereffect (tMAE) is a perceptual illusion in which a stationary stimulus feels as though it is moving when presented following adaptation to a unidirectionally moving tactile stimulus. Using functional magnetic resonance imaging (fMRI), we localized the brain areas responsive to tactile motion and then investigated whether these areas underlie the tMAE. Tactile stimulation was delivered to the glabrous surface of the right hand by means of a plastic cylinder with a square-wave patterned surface. In the tactile motion localizer, we contrasted periods in which the cylinder rotated at 15 rpm with periods of rest (stationary contact). Activation was observed in the contralateral (left) thalamus, postcentral gyrus, and parietal operculum. In the tMAE experiment, the cylinder rotated at 15 or 60 rpm for 2 min. The 60-rpm speed induced reliable tMAEs, whereas the 15-rpm speed did not. Of the areas activated by the tactile motion localizer, only the postcentral gyrus showed a sustained fMRI response following the offset of 60-rpm (but not 15-rpm) stimulation, presumably reflecting the illusory perception of motion.