人脑初级视皮层内事件相关型视觉刺激fMRI BOLD响应特性初步研究

基于脑血氧水平依赖(BOLD)对比的功能磁共振成像（fMRI）方法是研究脑功能活动的一种重要的无损伤探测手段，BOLD的机制及它与脑神经活动关系一直是国际上十分活跃的研究领域，尤其是在实验研究方面. 视觉刺激所引起的脑的初级视皮层(V1) 的BOLD响应的时间特性已有较多的研究，但是在这些研究中，有的结论认为BOLD对视觉刺激的响应是线性的，有的结论却是相反的. 我们采用事件关联型核磁共振功能成像（ER-fMRI）方法，研究不同的短暂视觉刺激持续时间下的BOLD响应，得到了视觉刺激下的脑激活图. 同时找出了视觉刺激时间分别是1、2、3、4、5和6 s的V1区的BOLD响应曲线，并初步显示出BOLD响应与刺激持续时间的线性关系.     

Dynamics and nonlinearities of the BOLD response at very short stimulus durations

In designing a functional imaging experiment or analyzing data, it is typically assumed that task duration and hemodynamic response are linearly related to each other. However, numerous human and animal studies have previously reported a deviation from linearity for short stimulus durations (<4 s). Here, we investigated nonlinearities of blood-oxygenation-level-dependent (BOLD) signals following visual stimulation of 5 to 1000 ms duration at two different luminance levels in human subjects. It was found that (a) a BOLD response to stimulus durations as short as 5 ms can be reliably detected; this stimulus duration is shorter than employed in any previous study investigating BOLD signal time courses; (b) the responses are more nonlinear than in any other previous study: the BOLD response to 1000 ms stimulation is only twice as large as the BOLD response to 5 ms stimulation although 200 times more photons were projected onto the retina; (c) the degree of nonlinearity depends on stimulus intensity; that is, nonlinearities have to be characterized not only by stimulus duration but also by stimulus features like luminance. These findings are especially of most practical importance in rapid event-related functional magnetic resonance imaging (fMRI) experimental designs. In addition, an 'initial dip' response--thought to be generated by a rapid increase in cerebral metabolic rate of oxygen metabolism (CMRO2) relative to cerebral blood flow--was observed and shown to colocalize well with the positive BOLD response. Highly intense stimulation, better tolerated by human subjects for short stimulus durations, causes early CMRO2 increase, and thus, the experimental design utilized in this study is better for detecting the initial dip than standard fMRI designs. These results and those from other groups suggest that short stimulation combined with appropriate experimental designs allows neuronal events and interactions to be examined by BOLD signal analysis, despite its slow evolution.     

Coupling of simultaneously acquired electrophysiological and haemodynamic responses during visual stimulation

We investigate the relationship between the temporal variation in the magnitude of occipital visual evoked potentials (VEPs) and of haemodynamic measures of brain activity obtained using both blood oxygenation level dependent (BOLD) and perfusion sensitive (ASL) functional magnetic resonance imaging (fMRI). Volunteers underwent a continuous BOLD fMRI scan and/or a continuous perfusion-sensitive (gradient and spin echo readout) ASL scan, during which 30 second blocks of contrast reversing visual stimuli (at 4 Hz) were interleaved with 30 second blocks of rest (visual fixation). Electroencephalography (EEG) and fMRI were simultaneously recorded and following EEG artefact cleaning, VEPs were averaged across the whole stimulation block (120 reversals, VEP₁₂₀) and at a finer timescale (15 reversals, VEP₁₅). Both BOLD and ASL time-series were linearly modelled to establish: (1) the mean response to visual stimulation, (2) transient responses at the start and end of each stimulation block, (3) the linear decrease between blocks, (4) the nonlinear between-block variation (covariation with VEP₁₂₀), (5) the linear decrease within block and (6) the nonlinear variation within block (covariation with VEP₁₅).     

Origins of intersubject variability of blood oxygenation level dependent and arterial spin labeling fMRI: implications for quantification of brain activity

Accurate localization of brain activity using blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) has been challenged because of the large BOLD signal within distal veins. Arterial spin labeling (ASL) techniques offer greater sensitivity to the microvasculature but possess low temporal resolution and limited brain coverage. In this study, we show that the physiological origins of BOLD and ASL depend on whether percent change or statistical significance is being considered. For BOLD and ASL fMRI data collected during a simple unilateral hand movement task, we found that in the area of the contralateral motor cortex the centre of gravity (CoG) of the intersubject coefficient of variation (CV) of BOLD fMRI was near the brain surface for percent change in signal, whereas the CoG of the intersubject CV for Z-score was in close proximity of sites of brain activity for both BOLD and ASL. These findings suggest that intersubject variability of BOLD percent change is vascular in origin, whereas the origin of inter-subject variability of Z-score is neuronal for both BOLD and ASL. For longer duration tasks (12 s or greater), however, there was a significant correlation between BOLD and ASL percent change, which was not evident for short duration tasks (6 s). These findings suggest that analyses directly comparing percent change in BOLD signal between pre-defined regions of interest using short duration stimuli, as for example in event-related designs, may be heavily weighted by large-vessel responses rather than neuronal responses.     

Measurement and characterization of the human spinal cord SEEP response using event-related spinal fMRI

Although event-related fMRI is able to reliably detect brief changes in brain activity and is now widely used throughout systems and cognitive neuroscience, there have been no previous reports of event-related spinal cord fMRI. This is likely attributable to the various technical challenges associated with spinal fMRI (e.g., imaging a suitable length of the cord, reducing image artifacts from the vertebrae and intervertebral discs, and dealing with physiological noise from spinal cord motion). However, with many of these issues now resolved, the largest remaining impediment for event-related spinal fMRI is a deprived understanding of the spinal cord fMRI signal time course. Therefore, in this study, we used a proton density-weighted HASTE sequence, with functional contrast based on signal enhancement by extravascular water protons (SEEP), and a motion-compensating GLM analysis to (i) characterize the SEEP response function in the human cervical spinal cord and (ii) demonstrate the feasibility of event-related spinal fMRI. This was achieved by applying very brief (1 s) epochs of 22°C thermal stimulation to the palm of the hand and measuring the impulse response function. Our results suggest that the spinal cord SEEP response (time to peak ≈8 s; FWHM ≈4 s; and probably lacking pre- and poststimulus undershoots) is slower than previous estimates of SEEP or BOLD responses in the brain, but faster than previously reported spinal cord BOLD responses. Finally, by detecting and mapping consistent signal-intensity changes within and across subjects, and validating these regions with a block-designed experiment, this study represents the first successful demonstration of event-related spinal fMRI.     

A nonlinear BOLD model accounting for refractory effect by applying the longitudinal relaxation in NMR to the linear BOLD model

A mathematical model to regress the nonlinear blood oxygen level-dependent (BOLD) fMRI signal has been developed by incorporating the refractory effect into the linear BOLD model of the biphasic gamma variate function. The refractory effect was modeled as a relaxation of two separate BOLD capacities corresponding to the biphasic components of the BOLD signal in analogy with longitudinal relaxation of magnetization in NMR. When tested with the published fMRI data of finger tapping, the nonlinear BOLD model with the refractory effect reproduced the nonlinear BOLD effects such as reduced poststimulus undershoot and saddle pattern in a prolonged stimulation as well as the reduced BOLD signal for repetitive stimulation.     

Integration of EEG source imaging and fMRI during continuous viewing of natural movies

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are noninvasive neuroimaging tools which can be used to measure brain activity with excellent temporal and spatial resolution, respectively. By combining the neural and hemodynamic recordings from these modalities, we can gain better insight into how and where the brain processes complex stimuli, which may be especially useful in patients with different neural diseases. However, due to their vastly different spatial and temporal resolutions, the integration of EEG and fMRI recordings is not always straightforward. One fundamental obstacle has been that paradigms used for EEG experiments usually rely on event-related paradigms, while fMRI is not limited in this regard. Therefore, here we ask whether one can reliably localize stimulus-driven EEG activity using the continuously varying feature intensities occurring in natural movie stimuli presented over relatively long periods of time. Specifically, we asked whether stimulus-driven aspects in the EEG signal would be co-localized with the corresponding stimulus-driven BOLD signal during free viewing of a movie. Secondly, we wanted to integrate the EEG signal directly with the BOLD signal, by estimating the underlying impulse response function (IRF) that relates the BOLD signal to the underlying current density in the primary visual area (V1). We made sequential fMRI and 64-channel EEG recordings in seven subjects who passively watched 2-min-long segments of a James Bond movie. To analyze EEG data in this natural setting, we developed a method based on independent component analysis (ICA) to reject EEG artifacts due to blinks, subject movement, etc., in a way unbiased by human judgment. We then calculated the EEG source strength of this artifact-free data at each time point of the movie within the entire brain volume using low-resolution electromagnetic tomography (LORETA). This provided for every voxel in the brain (i.e., in 3D space) an estimate of the current density at every time point. We then carried out a correlation between the time series of visual contrast changes in the movie with that of EEG voxels. We found the most significant correlations in visual area V1, just as seen in previous fMRI studies (Bartels A, Zeki, S, Logothetis NK. Natural vision reveals regional specialization to local motion and to contrast-invariant, global flow in the human brain. Cereb Cortex 2008;18(3):705–717), but on the time scale of milliseconds rather than of seconds. To obtain an estimate of how the EEG signal relates to the BOLD signal, we calculated the IRF between the BOLD signal and the estimated current density in area V1. We found that this IRF was very similar to that observed using combined intracortical recordings and fMRI experiments in nonhuman primates. Taken together, these findings open a new approach to noninvasive mapping of the brain. It allows, firstly, the localization of feature-selective brain areas during natural viewing conditions with the temporal resolution of EEG. Secondly, it provides a tool to assess EEG/BOLD transfer functions during processing of more natural stimuli. This is especially useful in combined EEG/fMRI experiments, where one can now potentially study neural-hemodynamic relationships across the whole brain volume in a noninvasive manner.     

Functional MR imaging assessment of a non-responsive brain injured patient.

Functional magnetic resonance imaging (fMRI) was requested to assist in the evaluation of a comatose 38-year-old woman who had sustained multiple cerebral contusions from a motor vehicle accident. Previous electrophysiologic studies suggested absence of thalamocortical processing in response to median nerve stimulation. Whole-brain fMRI was performed utilizing visual, somatosensory, and auditory stimulation paradigms. Results demonstrated intact task-correlated sensory and cognitive blood oxygen level dependent (BOLD) hemodynamic response to stimuli. Electrodiagnostic studies were repeated and evoked potentials indicated supratentorial recovery in the cerebrum. At 3-months post trauma the patient had recovered many cognitive & sensorimotor functions, accurately reflecting the prognostic fMRI evaluation. These results indicate that fMRI examinations may provide a useful evaluation for brain function in non-responsive brain trauma patients.     

Effects of mild hypoxic hypoxia on poststimulus undershoot of blood-oxygenation-level-dependent fMRI signal in the human visual cortex

Characteristics of the blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal poststimulus undershoot in the visual cortex were studied at varying levels of arterial blood oxygen saturation (Ysat). Undershoot with an amplitude of -0.6+/-0.2% appeared after positive BOLD response (+1.7+/-0.5%) under control conditions. Cerebral blood volume (CBV), as determined with vascular-space-occupancy-dependent fMRI, increased by 26-43% during the positive BOLD peak, but the CBV proceeded at baseline level during the BOLD poststimulus undershoot. Mild hypoxic hypoxia (Ysat ranging from 0.82 to 0.89) had no effect on the amplitude or duration of poststimulus undershoot in activated BOLD pixels. Hypoxia did not influence CBV during the BOLD poststimulus undershoot. In contrast, the positive BOLD signal at the level of all activated pixels was smaller in hypoxia than in normoxia. The present results show that the BOLD poststimulus undershoot is not influenced by curtailed oxygen availability and that, during the undershoot, CBV is not different from resting state.     

Steady-state activation in somatosensory cortex after changes in stimulus rate during median nerve stimulation

Passive electrical stimulation activates various human somatosensory cortical systems including the contralateral primary somatosensory area (SI), bilateral secondary somatosensory area (SII) and bilateral insula. The effect of stimulation frequency on blood oxygenation level-dependent (BOLD) activity remains unclear. We acquired 3-T functional magnetic resonance imaging (fMRI) in eight healthy volunteers during electrical median nerve stimulation at frequencies of 1, 3 and 10 Hz. During stimulation BOLD signal changes showed activation in the contralateral SI, bilateral SII and bilateral insula. Results of fMRI analysis showed that these areas were progressively active with the increase of rate of stimulation. As a major finding, the contralateral SI showed an increase of peak of BOLD activation from 1 to 3 Hz but reached a plateau during 10-Hz stimulation. Our finding is of interest for basic research and for clinical applications in subjects unable to perform cognitive tasks in the fMRI scanner.     

Impact of intravenous nicotine on BOLD signal response to photic stimulation

Functional magnetic resonance imaging (fMRI) is increasingly being applied in the study of brain effects of nicotine. In addition, because tobacco smoking is common, many subjects studied with fMRI for other reasons may have appreciable levels of nicotine in plasma and brain during scanning. However, there is concern that the vascular effects of nicotine may alter the coupling between blood oxygen level dependent (BOLD) signal and neuronal activity. The objective of this study was to test for evidence of alteration of BOLD signal response of occipital cortex, a region with a relatively low concentration of neuronal nicotine receptors, to photic stimulation during intravenous infusion of nicotine. Nine nicotine dependent healthy smokers were withdrawn from nicotine under controlled conditions and then scanned while receiving photic stimulation and successive intravenous infusions of saline and nicotine. No evidence for an effect of nicotine on BOLD signal response to photic stimulation was detected at the doses studied. This observation suggests that nicotine does not alter the coupling between BOLD signal and neuronal activity in the visual cortex.     

Identification of activated regions during a language task

Functional magnetic resonance imaging (fMRI) techniques are based on the assumption that changes in neural activity are accompanied by modulation in the blood-oxygenation-level-dependent (BOLD) signal. In addition to conventional increases in BOLD signals, sustained negative BOLD signal changes are occasionally observed in many fMRI experiments, which show regions of cortex that seem to respond in antiphase with primary stimulus. The existence of this so-called negative BOLD response (NBR) has been observed and investigated in many functional studies. Several theoretical mechanisms have been proposed to account for it, but its origin has never been fully explained. In this study, the variability of fMRI activation, including the sources of the negative BOLD signal, during phonological and semantic language tasks, was investigated in six right-handed healthy subjects. We found significant activations in the brain regions, mainly in the left hemisphere, involved in the language stimuli [prominent in the inferior frontal gyrus, approximately Brodmann Areas (BA)7, BA44, BA45 and BA47, and in the precuneus]. Moreover, we observed activations in motor regions [precentral gyrus and supplementary motor area (SMA)], a result that suggests a specific role of these areas (particularly the SMA) in language processing. Functional analysis have also shown that certain brain regions, including the posterior cingulate cortex and the anterior cingulate cortex, have consistently greater activity during resting states compared to states of performing cognitive tasks. In our study, we observed diffuse NBR at the cortical level and a stronger negative response in correspondence to the main sinuses. These phenomena seem to be unrelated to a specific neural activity, appearing to be expressions of a mechanical variation in hemodynamics. We discussed about the importance of these responses that are anticorrelated with the stimulus. Our data suggest that particular care must be considered in the interpretation of fMRI findings, especially in the case of presurgical studies.     

Effects of hypoglycemia on human brain activation measured with fMRI

Functional magnetic resonance imaging (fMRI) was used to measure the effects of acute hypoglycemia caused by passive sensory stimulation on brain activation. Visual stimulation was used to generate blood-oxygen-level-dependent (BOLD) contrast, which was monitored during hyperinsulinemic hypoglycemic and euglycemic clamp studies. Hypoglycemia (50 +/- 1 mg glucose/dl) decreased the fMRI signal relative to euglycemia in 10 healthy human subjects: the fractional signal change was reduced by 28 +/- 12% (P < .05). These changes were reversed when euglycemia was restored. These data provide a basis of comparison for studies that quantify hypoglycemia-related changes in fMRI activity during cognitive tasks based on visual stimuli and demonstrate that variations in blood glucose levels may modulate BOLD signals in the healthy brain.     

Neonatal auditory activation detected by functional magnetic resonance imaging

The objective of this study was to detect auditory cortical activation in non-sedated neonates employing functional magnetic resonance imaging (fMRI). Using echo-planar functional brain imaging, subjects were presented with a frequency-modulated pure tone; the BOLD signal response was mapped in 5 mm-thick slices running parallel to the superior temporal gyrus. Twenty healthy neonates (13 term, 7 preterm) at term and 4 adult control subjects. Blood oxygen level-dependent (BOLD) signal in response to auditory stimulus was detected in all 4 adults and in 14 of the 20 neonates. FMRI studies of adult subjects demonstrated increased signal in the superior temporal regions during auditory stimulation. In contrast, signal decreases were detected during auditory stimulation in 9 of 14 newborns with BOLD response. fMRI can be used to detect brain activation with auditory stimulation in human infants.     

Neuronal event detection in fMRI time series using iterative deconvolution techniques

An iterative estimation algorithm for deconvolution of neuronal activity from Blood Oxygen Level Dependent (BOLD) time series data is presented. The algorithm requires knowledge of the hemodynamic impulse response function but does not require knowledge of the stimulation function. The method uses majorization-minimization of a cost function to find an optimal solution to the inverse problem. The cost function includes penalties for the l₁ norm, total variation and negativity. The algorithm is able to identify the occurrence of neuronal activity bursts from BOLD time series accurately. The accuracy of the algorithm was tested in simulations and experimental fMRI data using blocked and event-related designs. The simulations revealed that the algorithm is most sensitive to contrast-to-noise ratio levels and to errors in the assumed hemodynamic model and least sensitive to autocorrelation in the noise. Within normal fMRI conditions, the method is effective for event detection.     

Percept-related activity in the human somatosensory system: functional magnetic resonance imaging studies

In this paper, we review blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) studies addressing the neural correlates of touch, thermosensation, pain and the mechanisms of their cognitive modulation in healthy human subjects. There is evidence that fMRI signal changes can be elicited in the parietal cortex by stimulation of single mechanoceptive afferent fibers at suprathreshold intensities for conscious perception. Positive linear relationships between the amplitude or the spatial extents of BOLD fMRI signal changes, stimulus intensity and the perceived touch or pain intensity have been described in different brain areas. Some recent fMRI studies addressed the role of cortical areas in somatosensory perception by comparing the time course of cortical activity evoked by different kinds of stimuli with the temporal features of touch, heat or pain perception. Moreover, parametric single-trial functional MRI designs have been adopted in order to disentangle subprocesses within the nociceptive system.

Available evidence suggest that studies that combine fMRI with psychophysical methods may provide a valuable approach for understanding complex perceptual mechanisms and top-down modulation of the somatosensory system by cognitive factors specifically related to selective attention and to anticipation. The brain networks underlying somatosensory perception are complex and highly distributed. A deeper understanding of perceptual-related brain mechanisms therefore requires new approaches suited to investigate the spatial and temporal dynamics of activation in different brain regions and their functional interaction.     

Neural activity-induced modulation of BOLD poststimulus undershoot independent of the positive signal

Despite intense research on the blood oxygenation level-dependent (BOLD) signal underlying functional magnetic resonance imaging, our understanding of its physiological basis is far from complete. In this study, it was investigated whether the so-called poststimulus BOLD signal undershoot is solely a passive vascular effect or actively induced by neural responses. Prolonged static and flickering black-white checkerboard stimulation with isoluminant grey screen as baseline condition were employed on eight human subjects. Within the same region of interest, the positive BOLD time courses for static and flickering stimuli were identical over the entire stimulus duration. In contrast, the static stimuli exhibited no poststimulus BOLD signal undershoot, whereas the flickering stimuli caused a strong BOLD poststimulus undershoot. To ease the interpretation, we performed an additional study measuring both BOLD signal and cerebral blood flow (CBF) using arterial spin labeling. Also for CBF, a difference in the poststimulus period was found for the two stimuli. Thus, a passive blood volume effect as the only contributor to the poststimulus undershoot comes short in explaining the BOLD poststimulus undershoot phenomenon for this particular experiment. Rather, an additional active neuronal activation or deactivation can strongly modulate the BOLD poststimulus behavior. In summary, the poststimulus time course of BOLD signal could potentially be used to differentiate neuronal activity patterns that are otherwise indistinguishable using the positive evoked response.     

Comparing BOLD fMRI signal changes in the awake and anesthetized rat during electrical forepaw stimulation

The difference between awake curarized and alpha-chloralose anesthetized animals was studied with respect to the BOLD signal response in an fMRI experiment. By studying the activation of the cortex upon electrical forepaw stimulation in the same rat, but following consecutively applied curarization and alpha-chloralose anesthesia protocols, it was possible to compare quantitatively the effect of both immobilization protocols on the fMRI data. The largest BOLD signal change as a result of forepaw stimulation was found in the awake condition, however the activated areas are less specific than those in the anesthetized state leaving it more difficult to interpret.     

Direct measurement of oxygen extraction with fMRI using 6% CO2 inhalation

The blood-oxygenation-level-dependent (BOLD) signal is an indirect hemodynamic signal that is sensitive to cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen. Therefore, the BOLD signal amplitude and dynamics cannot be interpreted unambiguously without additional physiological measurements, and thus, there remains a need for a functional magnetic resonance imaging (fMRI) signal, which is more closely related to the underlying neuronal activity. In this study, we measured CBF with continuous arterial spin labeling, CBV with an exogenous contrast agent and BOLD combined with intracortical electrophysiological recording in the primary visual cortex of the anesthetized monkey. During inhalation of 6% CO2, it was observed that CBF and CBV are not further increased by a visual stimulus, although baseline CBF for 6% CO2 is below the maximal value of CBF. In contrast, the electrophysiological response to the stimulation was found to be preserved during hypercapnia. As a consequence, the simultaneously measured BOLD signal responds negatively to a visual stimulation for 6% CO2 inhalation in the same voxels responding positively during normocapnia. These observations suggest that the fMRI response to a sensory stimulus for 6% CO2 inhalation occurs in the absence of a hemodynamic response, and it therefore directly reflects oxygen extraction into the tissue.