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

The current status and trend of the functional magnetic resonance combined with stimulation in animals

As a non-radiative, non-invasive imaging technique, functional magnetic resonance imaging (fMRI) has excellent effects on studying the activation of blood oxygen levels and functional connectivity of the brain in human and animal models. Compared with resting-state fMRI, fMRI combined with stimulation could be used to assess the activation of specific brain regions and the connectivity of specific pathways and achieve better signal capture with a clear purpose and more significant results. Various fMRI methods and specific stimulation paradigms have been proposed to investigate brain activation in a specific state, such as electrical, mechanical, visual, olfactory, and direct brain stimulation. In this review, the studies on animal brain activation using fMRI combined with different stimulation methods were retrieved. The instruments, experimental parameters, anesthesia, and animal models in different stimulation conditions were summarized. The findings would provide a reference for studies on estimating specific brain activation using fMRI combined with stimulation.

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.

Optimization of a quadrature birdcage coil for functional imaging of squirrel monkey brain at 9.4T

A quadrature transmit/receive birdcage coil was optimized for squirrel monkey functional imaging at the high field of 9.4 T. The coil length was chosen to gain maximum coil efficiency/signal-to-noise ratio (SNR) and meanwhile provide enough homogenous RF field in the whole brain area. Based on the numerical simulation results, a 16-rung high-pass birdcage coil with the optimal length of 9 cm was constructed and evaluated on phantom and in vivo experiments. Compared to a general-purpose non-optimized coil, it exhibits approximately 25% in vivo SNR improvement. In addition to the volume coil, details about how to design and construct the associated animal preparation system were provided.

Variability and Reproducibility of Directed and Undirected Functional MRI Connectomes in the Human Brain

A growing number of studies are focusing on methods to estimate and analyze the functional connectome of the human brain. Graph theoretical measures are commonly employed to interpret and synthesize complex network-related information. While resting state functional MRI (rsfMRI) is often employed in this context, it is known to exhibit poor reproducibility, a key factor which is commonly neglected in typical cohort studies using connectomics-related measures as biomarkers. We aimed to fill this gap by analyzing and comparing the inter- and intra-subject variability of connectivity matrices, as well as graph-theoretical measures, in a large (n = 1003) database of young healthy subjects which underwent four consecutive rsfMRI sessions. We analyzed both directed (Granger Causality and Transfer Entropy) and undirected (Pearson Correlation and Partial Correlation) time-series association measures and related global and local graph-theoretical measures. While matrix weights exhibit a higher reproducibility in undirected, as opposed to directed, methods, this difference disappears when looking at global graph metrics and, in turn, exhibits strong regional dependence in local graphs metrics. Our results warrant caution in the interpretation of connectivity studies, and serve as a benchmark for future investigations by providing quantitative estimates for the inter- and intra-subject variabilities in both directed and undirected connectomic measures.

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.

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.

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.

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.

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.

Fine-scale functional connectivity in somatosensory cortex revealed by high-resolution fMRI

High-resolution functional magnetic resonance imaging (fMRI) at high field (9.4 T) has been used to measure functional connectivity between subregions within the primary somatosensory (SI) cortex of the squirrel monkey brain. The hand-face region within the SI cortex of the squirrel monkey has been previously well mapped with functional imaging and electrophysiological and anatomical methods, and the orderly topographic map of the hand region is characterized by a lateral to medial representation of individual digits in four subregions of areas 3a, 3b, 1 and 2. With submillimeter resolution, we are able to detect not only the separate islands of activation corresponding to vibrotactile stimulations of single digits but also, in subsequent acquisitions, the degree of correlation between voxels within the SI cortex in the resting state. The results suggest that connectivity patterns are very similar to stimulus-driven distributions of activity and that connectivity varies on the scale of millimeters within the same primary region. Connectivity strength is not a reflection of global larger-scale changes in blood flow and is not directly dependent on distance between regions. Preliminary electrophysiological recordings agree well with the fMRI data. In human studies at 7 T, high-resolution fMRI may also be used to identify the same subregions and assess responses to sensory as well as painful stimuli, and to measure connectivity dynamically before and after such stimulations.

Reproducibility of functional MRI localization within the human somatosensory cortex

PURPOSE

To evaluate the effect of fMRI localization approaches and region size on the reproducibility of digit localization in the human somatosensory cortex.

MATERIALS AND METHODS

Vibrotactile stimulation was applied to digits 2 and 4 producing cortical activation sites relating to each digit. Thirteen subjects were scanned twice on separate occasions in a 3 Tesla scanner using a voxel size of 2 mm. Regions of activity were thresholded to different sizes varying from 50 to 1000 voxels. Three measures of position were acquired from these regions: center of gravity (COG), center co-ordinate and peak voxel. To account for registration errors, Euclidean distance between the two digits was computed. Reproducibility was determined in terms of the 95% confidence interval for individual position in X, Y, and Z and also the distance between the two digit locations.

RESULTS

Region size of 200 most significant voxels was shown to have the best reproducibility. Center co-ordinate proved to be the most precise form of localizing activity with a 95% CI of 2.1 mm, 2.6 mm, and 3.1 mm in the X, Y, and Z axes. Euclidean distance between the center co-ordinates of the two digit activation sites was shown to be a reliable means of overcoming registration errors with a 95% CI of 1.7 mm.

CONCLUSION

This study shows a high level of reproducibility for fMRI localization in the somatosensory system.

Mapping plasticity in the forepaw digit barrel subfield of rat brains using functional MRI

The topographic organization of the forepaw barrel subfield in layer IV of rat primary somatosensory cortex (S1) is a good model for studying neural function and plasticity. The goal of this study was to test the feasibility of functional MRI (fMRI) to map the forepaw digit representations in the S1 of the rat and its plasticity after digit amputation. Three dimensional echo-planar imaging with 300 micron isotropic resolution at 11.7 T was used to achieve high signal-to-noise ratios and laminar layer resolution. By alternating electrical stimulation of the 2nd (D2) and 4th (D4) digits, functional activation in layer IV of the barrel subfields could be distinguished using a differential analysis. Furthermore, 2 and a half months after the amputation of the 3rd digit in baby rats, the overlapping area between D2 and D4 representations was increased. This indicates that the forepaw barrel subfield previously associated with the ablated digit is now associated with the representation of nearby digits, which is consistent with studies using electrophysiology and cytochrome oxidase staining.

High-resolution functional magnetic resonance imaging mapping of noxious heat and tactile activations along the central sulcus in New World monkeys

This study mapped the fine-scale functional representation of tactile and noxious heat stimuli in cortical areas around the central sulcus of anesthetized squirrel monkeys by using high-resolution blood oxygen level-dependent (BOLD) fMRI at 9.4T. Noxious heat (47.5°C) stimulation of digits evoked multiple spatially distinct and focal BOLD activations. Consistent activations were observed in areas 3a, 3b, 1, and 2, whereas less frequent activation was present in M1. Compared with tactile activations, thermal nociceptive activations covered more area and formed multiple foci within each functional area. In general, noxious heat activations in area 3b did not colocalize with tactile responses. The spatial relationships of heat and tactile activations in areas 3a and 1/2 varied across animals. Subsequent electrophysiological mapping confirmed that the evoked heat and tactile BOLD signals were somatotopically appropriate. The magnitude and temporal profiles of the BOLD signals to noxious heat stimuli differed across cortical areas. Comparatively late-peaking but stronger signals were observed in areas 3b and 2, whereas earlier-peaking but weaker signals were observed in areas 3a, 1, and M1. In sum, this study not only confirmed the involvement of somatosensory areas of 3a, 3b, and 1, but also identified the engagements of area 2 and M1 in the processing of heat nociceptive inputs. Differential BOLD response profiles of the individual cortical areas along the central sulcus suggest that these areas play different roles in the encoding of nociceptive inputs. Thermal nociceptive and tactile inputs may be processed by different clusters of neurons in different areas. To critically bridge animal and human pain studies, human fMRI was related to primate fMRI and electrophysiology of nociceptive processing, examining the functional role of the primary somatosensory cortex in heat nociception and demonstrating that subregion areas 3a, 3b, 1, 2, and M1 are responsive to noxious heat stimuli.