Nonlinear statistical properties of fMRI signals in human visual cortex during resting-state

The purpose of the current work is to study nonlinear dynamics in neuronal activity within human brain visual cortex based on blood-oxygen-level dependent (BOLD) contrast imaging. In particular, based on functional magnetic resonance imaging (fMRI) signals, measures of fractality, complexity, and state disorder are estimated from central and peripheral eccentricity bands across three visual areas. Statistical results from analysis of 48750 resting-state fMRI signals show evidence that nonlinear dynamics of neuronal activity in resting-state in central and peripheral eccentricity bands of human visual cortex are persistent. However, they exhibit heterogeneous variability across eccentricity bands and visual areas. Also, information content in first visual area is more ordered than in the second one, whilst information content in the third visual area is the least ordered. These interesting nonlinear statistical properties are a further step toward understanding neuronal activity and nonlinear dynamics in human brain visual cortex.     

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₁₅).     

Imaging the neural circuitry and chemical control of aggressive motivation

Background  

With the advent of functional magnetic resonance imaging (fMRI) in awake animals it is possible to resolve patterns of neuronal activity across the entire brain with high spatial and temporal resolution. Synchronized changes in neuronal activity across multiple brain areas can be viewed as functional neuroanatomical circuits coordinating the thoughts, memories and emotions for particular behaviors. To this end, fMRI in conscious rats combined with 3D computational analysis was used to identifying the putative distributed neural circuit involved in aggressive motivation and how this circuit is affected by drugs that block aggressive behavior.     

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.     

Preoperative localization of the sensorimotor area using independent component analysis of resting-state fMRI

Analysis of resting-state functional magnetic resonance imaging (fMRI) data is based on detecting low-frequency signal fluctuations in functionally connected brain areas. These synchronous fluctuations in resting-state networks have been observed in several studies with healthy subjects. In this study, we explored if independent component analysis (ICA) can be used to localize the sensorimotor area from resting-state fMRI data in patients with brain tumors. Finger-tapping activation task and resting-state blood-oxygenation-level-dependent fMRI data were acquired from 8 patients with brain tumors and 10 healthy volunteers. Sensorimotor task independent components (IC_task) were used to verify resting-state independent components (IC_rest) individually. In addition, sensorimotor IC_rests were compared between the groups and no significant differences were detected in volume, spatial correlation or temporal correlation. These results show that it is possible to localize a sensorimotor area from resting-state data using ICA in patients with brain tumors. This offers a complementary method for assessing the sensorimotor area in subjects with brain tumors who have difficulties in performing motor paradigms.     

Sensitivity limitations of high-resolution perfusion-based human fMRI at 7 Tesla

The study of the brain's functional organization at laminar and columnar level of the cortex with blood oxygenation-level dependent (BOLD) functional MRI (fMRI) is affected by the contribution of large veins downstream from the microvascular response to brain activity. Blood volume- and especially perfusion-based techniques may reduce this problem because of their reduced sensitivity to venous effects, but may not allow the same spatial resolution because of smaller signal changes associated with brain activity. Here we investigated the practical resolution limits of perfusion-weighted fMRI in human visual stimulation experiments. For this purpose, we used a highly sensitive, single-shot perfusion labeling (SSPL) technique at 7 T and compared sensitivity to detect visual activation at low (2 mm, n = 10) and high (1 mm, n = 8) nominal isotropic spatial, and 3 s temporal, resolution with BOLD in 5½-minute-long experiments. Despite the smaller absolute signal change with activation, 2 mm resolution SSPL yielded comparable sensitivity to BOLD. This was attributed to a superior suppression of physiological noise with SSPL. However, at 1 mm nominal resolution, SSPL sensitivity fell on average at least 42% below that of BOLD, and detection of visual activation was compromised. This is explained by the fact that at high resolution, with both techniques, typically thermal noise rather than physiological noise dominates sensitivity. The observed sensitivity loss implies that to perform 1-mm resolution, perfusion weighted fMRI with a robustness similar to BOLD, scan times that are almost 3 times longer than the comparable BOLD experiment are required. This is in line with or slightly better than previous comparisons between perfusion-weighted fMRI and BOLD. The lower sensitivity has to be weighed against the spatial fidelity advantages of high-resolution perfusion-weighted fMRI.     

Midazolam sedation increases fluctuation and synchrony of the resting brain BOLD signal

The blood oxygen level-dependent (BOLD) magnetic resonance signal of functional brain cortices is dominated by very low frequency (VLF) fluctuations in anesthetized child patients. The temporal synchrony of the BOLD signal is also higher in anesthetized children compared with awake adults. The origin of the synchronous fluctuations can be related to maturation, pathological status or the anesthesia used in the imaging. Two of the three confounding variables (maturation and pathology) were controlled in this study. The effect of midazolam (4+/-0.8 mg) sedation on the BOLD signal was assessed in 12 healthy adults (aged 24+/-1.5 years) at 1.5 T. The VLF fluctuation power and temporal synchrony of the BOLD signal increased significantly after the sedation in the auditory and visual cortices. The fast Fourier transformation power spectral baseline fit parameters of the BOLD signal were also found to change significantly after sedation. It is concluded that the VLF fluctuation and temporal synchrony of the BOLD signal become increased after sedation in functional brain regions.     

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.     

Regional homogeneity of resting-state fMRI contributes to both neurovascular and task activation variations

The task induced blood oxygenation level dependent signal changes observed using functional magnetic resonance imaging (fMRI) are critically dependent on the relationship between neuronal activity and hemodynamic response. Therefore, understanding the nature of neurovascular coupling is important when interpreting fMRI signal changes evoked via task. In this study, we used regional homogeneity (ReHo), a measure of local synchronization of the BOLD time series, to investigate whether the similarities of one voxel with the surrounding voxels are a property of neurovascular coupling. FMRI scans were obtained from fourteen subjects during bilateral finger tapping (FTAP), digit–symbol substitution (DSST) and periodic breath holding (BH) paradigm. A resting-state scan was also obtained for each of the subjects for 4 min using identical imaging parameters. Inter-voxel correlation analyses were conducted between the resting-state ReHo, resting-state amplitude of low frequency fluctuations (ALFF), BH responses and task activations within the masks related to task activations. There was a reliable mean voxel-wise spatial correlation between ReHo and other neurovascular variables (BH responses and ALFF). We observed a moderate correlation between ReHo and task activations (FTAP: r = 0.32; DSST: r = 0.22) within the task positive network and a small yet reliable correlation within the default mode network (DSST: r = − 0.08). Subsequently, a linear regression was used to estimate the contribution of ReHo, ALFF and BH responses to the task activated voxels. The unique contribution of ReHo was minimal. The results suggest that regional synchrony of the BOLD activity is a property that can explain the variance of neurovascular coupling and task activations; but its contribution to task activations can be accounted for by other neurovascular factors such as the ALFF.     

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

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

Functional phantom for fMRI: a feasibility study

Functional magnetic resonance imaging (fMRI) reveals changes in blood oxygen level-dependent (BOLD) signal after considerable processing. This paper describes the implementation and testing of an fMRI phantom where electric current applied to a thin wire within a proton-rich medium substituted BOLD distortion of the magnetic field; the scanner detects these two distortions as practically identical signal changes. The magnitude of the change depended on the current strength. The phantom has a number of possible applications. Signal changes across sessions, days, instruments and individuals could be monitored. Placing the phantom close to a subject during an fMRI experiment could allow differentiating sensitivity changes in the scanner due to instrumentation from changes in the subject's state and performance during the experiment. The spatial extent of brain activations and effects of various changes in the chain of image formation could be analyzed using current-induced "activations". Furthermore, the phantom could expedite fMRI sequence development by reducing the need to scan human subjects, who introduce uncertainty to the signal. Thus, this fMRI phantom could be useful for both cognitive fMRI studies and scanner calibration.     

Learning-based structurally-guided construction of resting-state functional correlation tensors

Functional magnetic resonance imaging (fMRI) measures changes in blood-oxygenation-level-dependent (BOLD) signals to detect brain activities. It has been recently reported that the spatial correlation patterns of resting-state BOLD signals in the white matter (WM) also give WM information often measured by diffusion tensor imaging (DTI). These correlation patterns can be captured using functional correlation tensor (FCT), which is analogous to the diffusion tensor (DT) obtained from DTI. In this paper, we propose a noise-robust FCT method aiming at further improving its quality, and making it eligible for further neuroscience study. The novel FCT estimation method consists of three major steps: First, we estimate the initial FCT using a patch-based approach for BOLD signal correlation to improve the noise robustness. Second, by utilizing the relationship between functional and diffusion data, we employ a regression forest model to learn the mapping between the initial FCTs and the corresponding DTs using the training data. The learned forest can then be applied to predict the DTI-like tensors given the initial FCTs from the testing fMRI data. Third, we re-estimate the enhanced FCT by utilizing the DTI-like tensors as a feedback guidance to further improve FCT computation. We have demonstrated the utility of our enhanced FCTs in Alzheimer's disease (AD) diagnosis by identifying mild cognitive impairment (MCI) patients from normal subjects.     

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.     

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.     

High spatial resolution BOLD fMRI using simultaneous multislice excitation with echo-shifting gradient echo at 7 Tesla

We introduce an accelerated gradient echo (GRE) sequence combining simultaneous multislice excitation (SMS) with echo-shifting technique for high spatial resolution blood oxygen level dependent (BOLD) functional MRI (fMRI). The simulation was conducted to optimize scan parameters. To validate the feasibility of the proposed technique, the visual and motor task experiments were performed at 7.0 Tesla (T). The single-shot EPI sequence was also applied in comparison with the proposed technique. The simulation results showed that an optimized flip angle of 9° provided maximal BOLD contrast for our scanning scheme, allowing low power deposition and SMS acceleration factor of 5. Additionally, parallel acquisition imaging with acceleration factor of 2 was utilized, which allowed a total acceleration factor of 10 in volunteer study. The experiment results showed that geometric distortion-free BOLD images with voxel size of 1.0 × 1.0 × 2.5 mm³ were obtained. Significant brain activation was identified in both visual and motor task experiments, which were in accordance with previous investigations. The proposed technique has potential for high spatial resolution fMRI at ultra-high field because of its sufficient BOLD sensitivity as well as improved acquisition speed over conventional GRE-based techniques.     

Activation mapping as a percentage of local excitation: fMRI stability within scans, between scans and across field strengths

Functional magnetic resonance imaging (fMRI) does not typically yield highly reproducible maps of brain activation. Maps can vary significantly even with constant scanning parameters and consistent task performance conditions (Liu et al., Magn. Reson. Med., 2004, 52:751-760). Reproducibility is even more of a problem when comparing fMRI signal magnitude and spatial extent of activation across scans involving different task performance levels, scan durations, pulse sequences or magnetic field strengths. In this report, the consistency of fMRI was reexamined by considering the relative spatial and temporal distribution of fMRI blood oxygen level dependent (BOLD) activation signals separately from the absolute magnitude of the activation signal in each brain area. Subjects repeatedly performed the same simple motor task but under a variety of imaging conditions, using both spiral and standard echo-planar pulse sequences and at 1.5- and 4.0-T magnetic field strengths. The results demonstrate that the absolute amplitude of BOLD statistical activation signals varied significantly across time and scanning conditions, but the relative spatial pattern of BOLD activation was highly reproducible across all conditions. Analysis of realistic simulated fMRI data sets indicates that stability of relative activation patterns could provide a useful tool for assessing the accuracy of fMRI maps.     

Frontoparietal activity with minimal decision and control in the awake macaque at 7 T

Previous imaging work has identified a frontoparietal network in the human brain involved in many different cognitive functions, as well as in simple updates of attended information. To determine whether a similar network is present in the monkey brain and direct future electrophysiological recordings, we examined the activation of frontoparietal areas during visual stimulation in the awake, fixating monkey. We measured activity with BOLD fMRI in three animals and analyzed the data individually for each animal and at group level. We found reliable activations in lateral prefrontal and parietal areas, even though task-related decision making was minimal, as a response to simple update of visual information. These activations were significant for each individual animal, as well as at group level. Similar to human imaging results the update of visual input was enough to activate an extensive network of frontoparietal cortex in the macaque brain, a network which is normally associated with complex cognitive control processes.     

Evaluation of preprocessing steps to compensate for magnetic field distortions due to body movements in BOLD fMRI

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is currently the dominant technique for non-invasive investigation of brain functions. One of the challenges with BOLD fMRI, particularly at high fields, is compensation for the effects of spatiotemporally varying magnetic field inhomogeneities (ΔB₀) caused by normal subject respiration and, in some studies, movement of the subject during the scan to perform tasks related to the functional paradigm. The presence of ΔB₀ during data acquisition distorts reconstructed images and introduces extraneous fluctuations in the fMRI time series that decrease the BOLD contrast-to-noise ratio. Optimization of the fMRI data-processing pipeline to compensate for geometric distortions is of paramount importance to ensure high quality of fMRI data. To investigate ΔB₀ caused by subject movement, echo-planar imaging scans were collected with and without concurrent motion of a phantom arm. The phantom arm was constructed and moved by the experimenter to emulate forearm motions while subjects remained still and observed a visual stimulation paradigm. These data were then subjected to eight different combinations of preprocessing steps. The best preprocessing pipeline included navigator correction, a complex phase regressor and spatial smoothing. The synergy between navigator correction and phase regression reduced geometric distortions better than either step in isolation and preconditioned the data to make them more amenable to the benefits of spatial smoothing. The combination of these steps provided a 10% increase in t-statistics compared to only navigator correction and spatial smoothing and reduced the noise and false activations in regions where no legitimate effects would occur.     

Depression recognition using resting-state and event-related fMRI signals

Purpose

This paper aimed to develop a method for depression detection using blood-oxygen-level-dependent (BOLD) response estimated from event-related signals and resting-state functional magnetic resonance imaging (fMRI) signals together.

Materials and Methods

Thirteen patients with unipolar depression and matched healthy subjects were recruited. Resting state data of each subject were collected. Thereafter, event-related paradigm was undertaken using sad facial stimuli. The resting-state fMRI signal was deemed as the baseline of each subject's activity. Coefficient marks were designed to sort and select temporal independent components of event-related signals. Thereafter, stimulus-evoked BOLD response components inside event-related signal were extracted and taken as features to discriminate depressive patients from healthy controls.

Results

Accuracy rate for depression recognition was 77.27% with P value of .017 for whole-brain analysis and 81.82% with P value of .009 for region-of-interest analysis. The effectiveness and the superiority of the proposed method for disease recognition were demonstrated via the performance comparison with three other typical methods.

Conclusions

The proposed model was effective in depression recognition.     

Fractional amplitude of low-frequency fluctuation changes in monkeys with spinal cord injury: A resting-state fMRI study

Purpose

Although functional magnetic resonance imaging (fMRI) has revealed that spinal cord injury (SCI) causes anomalous changes in task-induced brain activation, its effect during the resting state remains unclear. The aim of this study is to explore the changes of the brain resting-state function in non-human primates with unilateral SCI.

Materials and methods

Eleven adult female rhesus monkeys were subjected to resting-state fMRI: five with unilateral thoracic SCI and six healthy monkeys, to obtain the fractional amplitude of low-frequency fluctuations (fALFF) of the blood oxygenation level-dependent (BOLD) contrast signal to determine the influence of SCI on the cerebral resting-state function.

Results

The SCI-induced fALFF vary significantly in several encephalic regions, including the left cerebellum, the left thalamus, the right lateral geniculate nucleus, the right superior parietal lobule, and the posterior cingulate gyrus.

Conclusion

Analysis of the resting-state fMRI provides evidence of abnormal spontaneous brain activations in primates with SCI, which may help us understand the pathophysiologic mechanisms underlying the changes in neural plasticity in the central nervous system after SCI.