A cluster mass permutation test with contextual enhancement for fMRI activation detection

Gaussian random field (GRF)-based methods are commonly used for statistical inference and to control the family-wise error rate (FWE) in neuroimaging. They require that the error fields are reasonable lattice approximations to an underlying continuous multivariate Gaussian random field and have differentiable and invertible spatial autocorrelation function. Permutation test estimates the distribution of the test statistic from the data and adjusts automatically for the FWE. Here we present a new analysis procedure, the cluster mass permutation test with contextual enhancement (CMPCE), and compare it to GRF. In CMPCE, the data are first pre-whitened to remove temporal autocorrelations. The FWE rates, the cluster detection probability and delineation accuracy of CMPCE and GRF were compared using measured null data and null data containing simulated activations. We also applied both methods to an fMRI experiment where tactile somatosensory stimulation into the right hand was used. When analyzing the FWE using null data, both CMPCE and GRF gave significantly higher FWEs (CMPCE up to 0.12, GRF up to 0.18) than the nominal significance level 0.05, indicating that the pre-whitening, motion correction or high-pass filtering partially failed. In the simulated activation data, CMPCE gave less falsely classified voxels for the same cluster detection probability level than GRF. The maximal cluster detection probability was on the other hand higher in the GRF-based method. Both methods gave qualitatively similar results in the tactile fMRI data. CMPCE seems to be a promising fMRI analysis method, especially if high delineation accuracy is required.

Contextual clustering for analysis of functional MRI data

We present a contextual clustering procedure for statistical parametric maps (SPM) calculated from time varying three-dimensional images. The algorithm can be used for the detection of neural activations from functional magnetic resonance images (fMRI). An important characteristic of SPM is that the intensity distribution of background (nonactive area) is known whereas the distributions of activation areas are not. The developed contextual clustering algorithm divides an SPM into background and activation areas so that the probability of detecting false activations by chance is controlled, i.e., hypothesis testing is performed. Unlike the much used voxel-by-voxel testing, neighborhood information is utilized, an important difference. This is achieved by using a Markov random field prior and iterated conditional modes (ICM) algorithm. However, unlike in the conventional use of ICM algorithm, the classification is based only on the distribution of background. The results from our simulations and human fMRI experiments using visual stimulation demonstrate that a better sensitivity is achieved with a given specificity in comparison to the voxel-by-voxel thresholding technique. The algorithm is computationally efficient and can be used to detect and delineate objects from a noisy background in other applications.

Reproducibility of fMRI: effect of the use of contextual information

We studied the effect of use of contextual information on the reproducibility of the results in analysis of fMRI data. We used data from a repeated simple motor fMRI experiment. In the first approach, statistical parametric maps were computed from a spatially unsmoothed data and thresholded using a Bonferroni corrected threshold. In the second approach, the maps were computed from a spatially unsmoothed data but were segmented into nonactive and active regions using a spatial contextual clustering method. In the third approach, the statistical parametric maps were computed from spatially smoothed data and thresholded, using, optionally, a spatial extent threshold. The variation in the classification was largest in the Bonferroni thresholded statistical parametric maps. There were no significant differences in variation between statistical parametric maps generated with all the other methods. In addition to reproducibility, the detection rates of weak simulated activations in the presence of measured scanner and physiological noise were investigated. Contextual clustering method was the most sensitive method, while the least sensitive method was the Bonferroni corrected thresholding. Using simulated data, we demonstrated that the contextual clustering method preserves the shapes of activation regions better than the method using spatial smoothing of the data.

Working memory of auditory localization

To investigate brain mechanisms of sound location memory, we studied the distribution of brain activation with functional magnetic resonance imaging (fMRI) in subjects performing an audiospatial n-back task with three memory load levels. Working memory processing of audiospatial information activated areas in the superior, middle and inferior frontal gyri, and in the posterior parietal and middle temporal cortices. In a control experiment, fMRI during audio- and visuospatial 2-back task performances revealed only few differentially activated subregions between the two tasks. These results demonstrate that working memory processing of auditory locations involves a distributed network of brain areas and suggest that mnemonic processing of audio- and visuospatial information is directed along a common neural pathway in the posterior parietal and prefrontal cortices.

Spatiotemporal activity of a cortical network for processing visual motion revealed by MEG and fMRI

A sudden change in the direction of motion is a particularly salient and relevant feature of visual information. Extensive research has identified cortical areas responsive to visual motion and characterized their sensitivity to different features of motion, such as directional specificity. However, relatively little is known about responses to sudden changes in direction. Electrophysiological data from animals and functional imaging data from humans suggest a number of brain areas responsive to motion, presumably working as a network. Temporal patterns of activity allow the same network to process information in different ways. The present study in humans sought to determine which motion-sensitive areas are involved in processing changes in the direction of motion and to characterize the temporal patterns of processing within this network of brain regions. To accomplish this, we used both magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI). The fMRI data were used as supplementary information in the localization of MEG sources. The change in the direction of visual motion was found to activate a number of areas, each displaying a different temporal behavior. The fMRI revealed motion-related activity in areas MT+ (the human homologue of monkey middle temporal area and possibly also other motion sensitive areas next to MT), a region near the posterior end of the superior temporal sulcus (pSTS), V3A, and V1/V2. The MEG data suggested additional frontal sources. An equivalent dipole model for the generators of MEG signals indicated activity in MT+, starting at 130 ms and peaking at 170 ms after the reversal of the direction of motion, and then again at approximately 260 ms. Frontal activity began 0-20 ms later than in MT+, and peaked approximately 180 ms. Both pSTS and FEF+ showed long-duration activity continuing over the latency range of 200-400 ms. MEG responses in the region of V3A and V1/V2 were relatively small, and peaked at longer latencies than the initial peak in MT+. These data revealed characteristic patterns of activity in this cortical network for processing sudden changes in the direction of visual motion.

Activation of multiple cortical areas in response to somatosensory stimulation: combined magnetoencephalographic and functional magnetic resonance imaging

We combined information from functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) to assess which cortical areas and in which temporal order show macroscopic activation after right median nerve stimulation. Five healthy subjects were studied with the two imaging modalities, which both revealed significant activation in the contra- and ipsilateral primary somatosensory cortex (SI), the contra- and ipsilateral opercular areas, the walls of the contralateral postcentral sulcus (PoCS), and the contralateral supplementary motor area (SMA). In fMRI, two separate foci of activation in the opercular cortex were discerned, one posteriorly in the parietal operculum (PO), and one anteriorly near the insula or frontal operculum (anterior operculum, AO). The activation sites from fMRI were used to constrain the solution of the inverse problem of MEG, which allowed us to construct a model of the temporal sequence of activation of the different sites. According to this model, the mean onset latency for significant activation at the contralateral SI was 20 msec (range, 17-22 msec), followed by activation of PoCS at 23 msec (range, 21-25 msec). The contralateral PO was activated at 26 msec (range, 19-32 msec) and AO at 33 msec (range, 22-51 msec). The contralateral SMA became active at 36 msec (range, 24-48 msec). The ipsilateral SI, PO, and AO became activated at 54-67 msec. We conclude that fMRI provides a useful means to constrain the inverse problem of MEG, allowing the construction of spatiotemporal models of cortical activation, which may have significant implications for the understanding of cortical network functioning.

Activation of multiple cortical areas in response to somatosensory stimulation: combined magnetoencephalographic and functional magnetic resonance imaging

We combined information from functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) to assess which cortical areas and in which temporal order show macroscopic activation after right median nerve stimulation. Five healthy subjects were studied with the two imaging modalities, which both revealed significant activation in the contra- and ipsilateral primary somatosensory cortex (SI), the contra- and ipsilateral opercular areas, the walls of the contralateral postcentral sulcus (PoCS), and the contralateral supplementary motor area (SMA). In fMRI, two separate foci of activation in the opercular cortex were discerned, one posteriorly in the parietal operculum (PO), and one anteriorly near the insula or frontal operculum (anterior operculum, AO). The activation sites from fMRI were used to constrain the solution of the inverse problem of MEG, which allowed us to construct a model of the temporal sequence of activation of the different sites. According to this model, the mean onset latency for significant activation at the contralateral SI was 20 msec (range, 17-22 msec), followed by activation of PoCS at 23 msec (range, 21-25 msec). The contralateral PO was activated at 26 msec (range, 19-32 msec) and AO at 33 msec (range, 22-51 msec). The contralateral SMA became active at 36 msec (range, 24-48 msec). The ipsilateral SI, PO, and AO became activated at 54-67 msec. We conclude that fMRI provides a useful means to constrain the inverse problem of MEG, allowing the construction of spatiotemporal models of cortical activation, which may have significant implications for the understanding of cortical network functioning.

Distribution of cortical activation during visuospatial n-back tasks as revealed by functional magnetic resonance imaging

Human neuroimaging studies conducted during visuospatial working memory tasks have inconsistently detected activation in the prefrontal cortical areas depending presumably on the type of memory and control tasks employed. We used functional magnetic resonance imaging to study brain activation related to the performance of a visuospatial n-back task with different memory loads (0-back, 1-back and 2-back tasks). Comparison of the 2-back versus 0-back tasks revealed consistent, bilateral activation in the medial frontal gyrus (MFG), superior frontal sulcus and adjacent cortical tissue (SFS/SFG) in all subjects and in six out of seven subjects in the intraparietal sulcus (IPS). Activation was also detected in the inferior frontal gyrus, medially in the superior frontal gyrus, precentral gyrus, superior and inferior parietal lobuli, occipital visual association areas, anterior and posterior cingulate areas and in the insula. Comparison between the 1-back versus 0-back tasks revealed activation only in a few brain areas. Activation in the MFG, SFS/SFG and IPS appeared dependent on memory load. The results suggest that the performance of a visuospatial working memory task engages a network of distributed brain areas and that areas in the dorsal visual pathway are engaged in mnemonic processing of visuospatial information.

Somatosensory evoked magnetic fields from primary sensorimotor cortex in juvenile neuronal ceroid lipofuscinosis

The present study evaluated neurophysiologic function of the primary sensorimotor cortex in juvenile neuronal ceroid lipofuscinosis. A 122-channel magnetometer, which allowed studies of the somatosensory system in millimeter and millisecond precision, was used to record somatosensory evoked magnetic fields to median nerve stimulation from 10 patients and their matched control subjects. In both patients and controls, the somatosensory evoked magnetic fields from primary sensorimotor area typically consisted of N20m, P35m, and P60m deflections. In the patients, N20m was significantly delayed, whereas P35m peaked earlier than in the control subjects. The source strengths for N20m and P35m were greater in the patients than in the controls. Both deflections showed a significant positive correlation with the disease duration: the sources were stronger in the older patients than in the younger ones. P60m deflections were normal or reduced in the patients. The results indicated increased thalamocortical excitability in the sensorimotor cortex in juvenile neuronal ceroid lipofuscinosis.

Effects of interstimulus interval on somatosensory evoked magnetic fields (SEFs): a hypothesis concerning SEF generation at the primary sensorimotor cortex

Cerebral responses evoked by peripheral stimuli are known to depend critically on the interstimulus interval (ISI). Here we report on the effects of ISI on somatosensory evoked magnetic fields (SEFs) to right median nerve stimulation, obtained in 9 healthy adults with ISIs of 0.15 0.3, 1,3 and 5 s. At the contralateral (left) primary sensorimotor cortex (SMI), the first cortical response, N20m, was stable between the ISIs 0.3 and 5 s, but slightly attenuated at the shortest ISI of 0.15 s. In contrast, the P35m and P60m deflections were very sensitive to changes of the ISI, declining steadily with shortening of the ISI throughout the entire range. These deflections were frequently undetectable at the shortest ISI of 0.15 s. Concomitant with the reductions of P35m and P60m, an N45m deflection was enhanced toward the short ISIs. Responses from second somatosensory cortex (SII) and posterior parietal cortex (PPC) were seen only with ISIs of 1 s or greater, being strongest at the 5 s ISI. Based on known effects of the ISI on intracellular evoked potentials, we present the following tentative model for the generation mechanism of the SMI response: N20m represents early excitatory postsynaptic potentials (EPSPs), P35m early inhibitory postsynaptic potentials (IPSPs), N45m secondary EPSPs and P60m late IPSPs in pyramidal neurones of area 3b. For practical purposes, SEFs from SMI can be obtained with short ISIs, while responses from SII and PPC require an ISI of at least 1 s.

Significance of the second somatosensory cortex in sensorimotor integration: enhancement of sensory responses during finger movements

The functional significance of the second somatosensory cortex (SII) is poorly understood. However, lesion and cortical stimulation studies indicate that SII may be involved in sensory aspects of tactile learning and in movement control. In the present study, we explored a possible role of SII in sensorimotor integration in humans using a multichannel magnetometer. Somatosensory evoked fields (SEFs) from SII to electrical stimulation of left and right median nerves were recorded in six healthy volunteers during rest and in different test conditions. Continuous cutaneous stimulation of the right hand or face reduced the SEFs to both left and right median nerve stimulation. Right-sided finger movements increased the SEFs to right, but not left, median nerve stimulation. The responses were equally enhanced by simple finger flexion movement and by a complex finger sequence. The suppression of SEFs by competing cutaneous inputs from different areas of the body indicates that the neurones underlying the responses receive inputs from large, bilateral receptive fields. The enhancement of sensory reactions to signals from the actively moving limb but not to those from the opposite limb indicates a spatial tuning of the SII neurones to behaviourally relevant input channels, also suggesting that SII is important for the integration of sensory information to motor programmes.

Activation of ipsilateral primary sensorimotor cortex by median nerve stimulation

We report evidence for activation of ipsilateral primary sensorimotor cortex (SMI) after median nerve stimulation recorded with magnetoencephalography (MEG). We measured somatosensory evoked magnetic fields (SEFs) to median nerve stimulation with a 122-channel helmet-shaped magnetometer in 10 healthy subjects. In five, the magnetic field patterns suggested long-latency activation of the ipsilateral SMI. Source locations found by current dipole fitting corresponded to the SMI hand area, as determined by contralateral stimulation. Further evidence for the origin of the ipsilateral responses in SMI was provided by the suppression of these responses during movement of the contralateral fingers. Sensory input to ipsilateral SMI could play a role in sensorimotor integration of bilateral movements.