Transient and sustained components of the sensorimotor BOLD response in fMRI

Differential Effects of Sustained Manual Pressure Stimulation According to Site of Action

Sustained pressure stimulation of the body surface has been used in several physiotherapeutic techniques, such as reflex locomotion therapy. Clinical observations of global motor responses and subsequent motor behavioral changes after stimulation in certain sites suggest modulation of central sensorimotor control, however, the neuroanatomical correlates remain undescribed. We hypothesized that different body sites would specifically influence the sensorimotor system during the stimulation. We tested the hypothesis using functional magnetic resonance imaging (fMRI) in thirty healthy volunteers (mean age 24.2) scanned twice during intermittent manual pressure stimulation, once at the right lateral heel according to reflex locomotion therapy, and once at the right lateral ankle (control site). A flexible modeling approach with finite impulse response basis functions was employed since non-canonical hemodynamic response was expected. Subsequently, a clustering algorithm was used to separate areas with differential timecourses. Stimulation at both sites induced responses throughout the sensorimotor system that could be mostly separated into two anti-correlated subsystems with transient positive or negative signal change and rapid adaptation, although in heel stimulation, insulo-opercular cortices and pons showed sustained activation. In direct voxel-wise comparison, heel stimulation was associated with significantly higher activation levels in the contralateral primary motor cortex and decreased activation in the posterior parietal cortex. Thus, we demonstrate that the manual pressure stimulation affects multiple brain structures involved in motor control and the choice of stimulation site impacts the shape and amplitude of the blood oxygenation level-dependent response. We further discuss the relationship between the affected structures and behavioral changes after reflex locomotion therapy.

Comparison of fMRI analysis methods for heterogeneous BOLD responses in block design studies

A large number of fMRI studies have shown that the temporal dynamics of evoked BOLD responses can be highly heterogeneous. Failing to model heterogeneous responses in statistical analysis can lead to significant errors in signal detection and characterization and alter the neurobiological interpretation. However, to date it is not clear that, out of a large number of options, which methods are robust against variability in the temporal dynamics of BOLD responses in block-design studies. Here, we used rodent optogenetic fMRI data with heterogeneous BOLD responses and simulations guided by experimental data as a means to investigate different analysis methods' performance against heterogeneous BOLD responses. Evaluations are carried out within the general linear model (GLM) framework and consist of standard basis sets as well as independent component analysis (ICA). Analyses show that, in the presence of heterogeneous BOLD responses, conventionally used GLM with a canonical basis set leads to considerable errors in the detection and characterization of BOLD responses. Our results suggest that the 3rd and 4th order gamma basis sets, the 7th to 9th order finite impulse response (FIR) basis sets, the 5th to 9th order B-spline basis sets, and the 2nd to 5th order Fourier basis sets are optimal for good balance between detection and characterization, while the 1st order Fourier basis set (coherence analysis) used in our earlier studies show good detection capability. ICA has mostly good detection and characterization capabilities, but detects a large volume of spurious activation with the control fMRI data.

Functional Transcranial Doppler Ultrasound for Measurement of Hemispheric Lateralization During Visual Memory and Visual Search Cognitive Tasks

Functional transcranial Doppler ultrasound (fTCD) is a noninvasive sensing modality that measures cerebral blood flow velocity (CBFV) with high temporal resolution. CBFV change is correlated to changes in cerebral oxygen uptake, enabling fTCD to measure brain activity and lateralization with high accuracy. However, few studies have examined the relationship of CBFV change during visual search and visual memory tasks. Here a protocol to compare lateralization between these two similar cognitive tasks using fTCD is demonstrated. Ten healthy volunteers (age 21±2 years) were shown visual scenes on a computer and performed visual search and visual memory tasks while CBFV in the bilateral middle cerebral arteries was monitored with fTCD. Each subject completed 40 trials, consisting of baseline (25 s), calibration (variable), instruction (2.5 s), and task (20 s) epochs. Lateralization was computed for each task by calculating the bilateral CBFV envelope percent change from baseline and subtracting the right side from the left side. The results showed significant lateralization ( ) of the visual memory and visual search tasks, with memory reaching lateralization of 1.6% and search reaching lateralization of 0.5%, suggesting that search is more right lateralized (and therefore may be related to "holistic" or global perception) and memory is more left lateralized (and therefore may be related to local perception). This method could be used to compare cerebral activity for any related cognitive tasks as long as the same stimulus is used in all tasks. The protocol is straightforward and the equipment is inexpensive, introducing a low-cost high temporal resolution technique to further study lateralization of the brain.

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.

Task Dependence, Tissue Specificity, and Spatial Distribution of Widespread Activations in Large Single-Subject Functional MRI Datasets at 7T

It was recently shown that when large amounts of task-based blood oxygen level-dependent (BOLD) data are combined to increase contrast- and temporal signal-to-noise ratios, the majority of the brain shows significant hemodynamic responses time-locked with the experimental paradigm. Here, we investigate the biological significance of such widespread activations. First, the relationship between activation extent and task demands was investigated by varying cognitive load across participants. Second, the tissue specificity of responses was probed using the better BOLD signal localization capabilities of a 7T scanner. Finally, the spatial distribution of 3 primary response types--namely positively sustained (pSUS), negatively sustained (nSUS), and transient--was evaluated using a newly defined voxel-wise waveshape index that permits separation of responses based on their temporal signature. About 86% of gray matter (GM) became significantly active when all data entered the analysis for the most complex task. Activation extent scaled with task load and largely followed the GM contour. The most common response type was nSUS BOLD, irrespective of the task. Our results suggest that widespread activations associated with extremely large single-subject functional magnetic resonance imaging datasets can provide valuable information about the functional organization of the brain that goes undetected in smaller sample sizes.

Identifying effective connectivity parameters in simulated fMRI: a direct comparison of switching linear dynamic system, stochastic dynamic causal, and multivariate autoregressive models

The number and variety of connectivity estimation methods is likely to continue to grow over the coming decade. Comparisons between methods are necessary to prune this growth to only the most accurate and robust methods. However, the nature of connectivity is elusive with different methods potentially attempting to identify different aspects of connectivity. Commonalities of connectivity definitions across methods upon which base direct comparisons can be difficult to derive. Here, we explicitly define "effective connectivity" using a common set of observation and state equations that are appropriate for three connectivity methods: dynamic causal modeling (DCM), multivariate autoregressive modeling (MAR), and switching linear dynamic systems for fMRI (sLDSf). In addition while deriving this set, we show how many other popular functional and effective connectivity methods are actually simplifications of these equations. We discuss implications of these connections for the practice of using one method to simulate data for another method. After mathematically connecting the three effective connectivity methods, simulated fMRI data with varying numbers of regions and task conditions is generated from the common equation. This simulated data explicitly contains the type of the connectivity that the three models were intended to identify. Each method is applied to the simulated data sets and the accuracy of parameter identification is analyzed. All methods perform above chance levels at identifying correct connectivity parameters. The sLDSf method was superior in parameter estimation accuracy to both DCM and MAR for all types of comparisons.

Effects of age on negative BOLD signal changes in the primary somatosensory cortex

In addition to a contralateral activation of the primary and secondary somatosensory cortices, peripheral sensory stimulation has been shown to elicit responses in the ipsilateral primary somatosensory cortex (SI). In particular, evidence is accumulating that processes of interhemispheric inhibition as depicted by negative blood oxygenation level dependent (BOLD) signal changes are part of somatosensory processes. The aim of the study was to analyze age-related differences in patterns of cerebral activation in the somatosensory system in general and processes of interhemispheric inhibition in particular. For this, a functional magnetic resonance imaging (fMRI) study was performed including 14 younger (mean age 23.3±0.9years) and 13 healthy older participants (mean age 73.2±8.3years). All subjects were scanned during peripheral electrical median nerve stimulation (40Hz) to obtain BOLD responses in the somatosensory system. Moreover, the individual current perception threshold (CPT) as a quantitative measure of sensory function was determined in a separate psychophysical testing. Significant increases in BOLD signal across the entire group could be measured within the contralateral SI, in the bilateral secondary somatosensory cortex (SII), the contralateral supplementary motor area and the insula. Negative BOLD signal changes were delineated in ipsilateral SI/MI as well as in the ipsilateral thalamus and basal ganglia. After comparing the two groups, only the cortical deactivation in ipsilateral SI in the early stimulation phase as well as the activation in contralateral SI and SII in the late stimulation block remained as statistically significant differences between the two groups. The psychophysical experiments yielded a significant age-dependent effect of CPT change with less difference in the older group which is in line with the significantly smaller alterations in maximal BOLD signal change in the contra- and ipsilateral SI found between the two groups. Healthy aging seems to be associated with a decrease in intracerebral inhibition as reflected by smaller negative BOLD signal changes during fMRI tasks. This finding could constitute an important link between age-related neurophysiological changes and behavioral alterations in humans.