Is optical imaging spectroscopy a viable measurement technique for the investigation of the negative BOLD phenomenon? A concurrent optical imaging spectroscopy and fMRI study at high field (7T)

Traditionally functional magnetic resonance imaging (fMRI) has been used to map activity in the human brain by measuring increases in the Blood Oxygenation Level Dependent (BOLD) signal. Often accompanying positive BOLD fMRI signal changes are sustained negative signal changes. Previous studies investigating the neurovascular coupling mechanisms of the negative BOLD phenomenon have used concurrent 2D-optical imaging spectroscopy (2D-OIS) and electrophysiology (Boorman et al., 2010). These experiments suggested that the negative BOLD signal in response to whisker stimulation was a result of an increase in deoxy-haemoglobin and reduced multi-unit activity in the deep cortical layers. However, Boorman et al. (2010) did not measure the BOLD and haemodynamic response concurrently and so could not quantitatively compare either the spatial maps or the 2D-OIS and fMRI time series directly. Furthermore their study utilised a homogeneous tissue model in which is predominantly sensitive to haemodynamic changes in more superficial layers. Here we test whether the 2D-OIS technique is appropriate for studies of negative BOLD. We used concurrent fMRI with 2D-OIS techniques for the investigation of the haemodynamics underlying the negative BOLD at 7 Tesla. We investigated whether optical methods could be used to accurately map and measure the negative BOLD phenomenon by using 2D-OIS haemodynamic data to derive predictions from a biophysical model of BOLD signal changes. We showed that despite the deep cortical origin of the negative BOLD response, if an appropriate heterogeneous tissue model is used in the spectroscopic analysis then 2D-OIS can be used to investigate the negative BOLD phenomenon.

BMI not WHR modulates BOLD fMRI responses in a sub-cortical reward network when participants judge the attractiveness of human female bodies

In perceptual terms, the human body is a complex 3d shape which has to be interpreted by the observer to judge its attractiveness. Both body mass and shape have been suggested as strong predictors of female attractiveness. Normally body mass and shape co-vary, and it is difficult to differentiate their separate effects. A recent study suggested that altering body mass does not modulate activity in the reward mechanisms of the brain, but shape does. However, using computer generated female body-shaped greyscale images, based on a Principal Component Analysis of female bodies, we were able to construct images which covary with real female body mass (indexed with BMI) and not with body shape (indexed with WHR), and vice versa. Twelve observers (6 male and 6 female) rated these images for attractiveness during an fMRI study. The attractiveness ratings were correlated with changes in BMI and not WHR. Our primary fMRI results demonstrated that in addition to activation in higher visual areas (such as the extrastriate body area), changing BMI also modulated activity in the caudate nucleus, and other parts of the brain reward system. This shows that BMI, not WHR, modulates reward mechanisms in the brain and we infer that this may have important implications for judgements of ideal body size in eating disordered individuals.

Assessment of the Spatial and Temporal Distribution of Functional Connectivity in Resting-State BOLD fMRI

Recent research into resting-state functional magnetic resonance imaging (fMRI) has shown that the brain is very active during rest. This thesis work utilizes blood oxygenation level dependent (BOLD) signals to investigate the spatial and temporal functional network information found within resting-state data, and aims to investigate the feasibility of extracting functional connectivity networks using different methods as well as the dynamic variability within some of the methods. Furthermore, this work looks into producing valid networks using a sparsely-sampled sub-set of the original data.

In this work we utilize four main methods: independent component analysis (ICA), principal component analysis (PCA), correlation, and a point-processing technique. Each method comes with unique assumptions, as well as strengths and limitations into exploring how the resting state components interact in space and time.

Correlation is perhaps the simplest technique. Using this technique, resting-state patterns can be identified based on how similar the time profile is to a seed region’s time profile. However, this method requires a seed region and can only identify one resting state network at a time. This simple correlation technique is able to reproduce the resting state network using subject data from one subject’s scan session as well as with 16 subjects.

Independent component analysis, the second technique, has established software programs that can be used to implement this technique. ICA can extract multiple components from a data set in a single analysis. The disadvantage is that the resting state networks it produces are all independent of each other, making the assumption that the spatial pattern of functional connectivity is the same across all the time points. ICA is successfully able to reproduce resting state connectivity patterns for both one subject and a 16 subject concatenated data set.

Using principal component analysis, the dimensionality of the data is compressed to find the directions in which the variance of the data is most significant. This method utilizes the same basic matrix math as ICA with a few important differences that will be outlined later in this text. Using this method, sometimes different functional connectivity patterns are identifiable but with a large amount of noise and variability.

To begin to investigate the dynamics of the functional connectivity, the correlation technique is used to compare the first and second halves of a scan session. Minor differences are discernable between the correlation results of the scan session halves. Further, a sliding window technique is implemented to study the correlation coefficients through different sizes of correlation windows throughout time. From this technique it is apparent that the correlation level with the seed region is not static throughout the scan length.

The last method introduced, a point processing method, is one of the more novel techniques because it does not require analysis of the continuous time points. Here, network information is extracted based on brief occurrences of high or low amplitude signals within a seed region. Because point processing utilizes less time points from the data, the statistical power of the results is lower. There are also larger variations in DMN patterns between subjects. In addition to boosted computational efficiency, the benefit of using a point-process method is that the patterns produced for different seed regions do not have to be independent of one another.

This work compares four unique methods of identifying functional connectivity patterns. ICA is a technique that is currently used by many scientists studying functional connectivity patterns. The PCA technique is not optimal for the level of noise and the distribution of the data sets. The correlation technique is simple and obtains good results, however a seed region is needed and the method assumes that the DMN regions is correlated throughout the entire scan. Looking at the more dynamic aspects of correlation changing patterns of correlation were evident. The last point-processing method produces a promising results of identifying functional connectivity networks using only low and high amplitude BOLD signals.

The Stability of the Blood Oxygenation Level-Dependent Functional MRI Response to Motor Tasks Is Altered in Patients With Chronic Ischemic Stroke

Background and Purpose-Functional MRI is a powerful tool to investigate recovery of brain function in patients with stroke. An inherent assumption in functional MRI data analysis is that the blood oxygenation level-dependent (BOLD) signal is stable over the course of the examination. In this study, we evaluated the validity of such assumption in patients with chronic stroke. Methods-Fifteen patients performed a simple motor task with repeated epochs using the paretic and the unaffected hand in separate runs. The corresponding BOLD signal time courses were extracted from the primary and supplementary motor areas of both hemispheres. Statistical maps were obtained by the conventional General Linear Model and by a parametric General Linear Model. Results-Stable BOLD amplitude was observed when the task was executed with the unaffected hand. Conversely, the BOLD signal amplitude in both primary and supplementary motor areas was progressively attenuated in every patient when the task was executed with the paretic hand. The conventional General Linear Model analysis failed to detect brain activation during movement of the paretic hand. However, the proposed parametric General Linear Model corrected the misdetection problem and showed robust activation in both primary and supplementary motor areas. Conclusions-The use of data analysis tools that are built on the premise of a stable BOLD signal may lead to misdetection of functional regions and underestimation of brain activity in patients with stroke. The present data urge the use of caution when relying on the BOLD response as a marker of brain reorganization in patients with stroke. (Stroke. 2010; 41:1921-1926.)

The resting-state neurovascular coupling relationship: rapid changes in spontaneous neural activity in the somatosensory cortex are associated with haemodynamic fluctuations that resemble stimulus-evoked haemodynamics

Although promise exists for patterns of resting-state blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) brain connectivity to be used as biomarkers of early brain pathology, a full understanding of the nature of the relationship between neural activity and spontaneous fMRI BOLD fluctuations is required before such data can be correctly interpreted. To investigate this issue, we combined electrophysiological recordings of rapid changes in multi-laminar local field potentials from the somatosensory cortex of anaesthetized rats with concurrent two-dimensional optical imaging spectroscopy measurements of resting-state haemodynamics that underlie fluctuations in the BOLD fMRI signal. After neural ‘events’ were identified, their time points served to indicate the start of an epoch in the accompanying haemodynamic fluctuations. Multiple epochs for both neural ‘events’ and the accompanying haemodynamic fluctuations were averaged. We found that the averaged epochs of resting-state haemodynamic fluctuations taken after neural ‘events’ closely resembled the temporal profile of stimulus-evoked cortical haemodynamics. Furthermore, we were able to demonstrate that averaged epochs of resting-state haemodynamic fluctuations resembling the temporal profile of stimulus-evoked haemodynamics could also be found after peaks in neural activity filtered into specific electroencephalographic frequency bands (theta, alpha, beta, and gamma). This technique allows investigation of resting-state neurovascular coupling using methodologies that are directly comparable to that developed for investigating stimulus-evoked neurovascular responses.

From EEG to BOLD: Brain mapping and estimating transfer functions in simultaneous EEG-fMRI acquisitions

Simultaneous acquisition of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) aims to disentangle the description of brain processes by exploiting the advantages of each technique. Most studies in this field focus on exploring the relationships between fMRI signals and the power spectrum at some specific frequency bands (alpha, beta, etc.). On the other hand, brain mapping of EEG signals (e.g., interictal spikes in epileptic patients) usually assumes an haemodynamic response function for a parametric analysis applying the GLM, as a rough approximation. The integration of the information provided by the high spatial resolution of MR images and the high temporal resolution of EEG may be improved by referencing them by transfer functions, which allows the identification of neural driven areas without strong assumptions about haemodynamic response shapes or brain haemodynamic`s homogeneity. The difference on sampling rate is the first obstacle for a full integration of EEG and fMRI information. Moreover, a parametric specification of a function representing the commonalities of both signals is not established. In this study, we introduce a new data-driven method for estimating the transfer function from EEG signal to fMRI signal at EEG sampling rate. This approach avoids EEG subsampling to fMRI time resolution and naturally provides a test for EEG predictive power over BOLD signal fluctuations, in a well-established statistical framework. We illustrate this concept in resting state (eyes closed) and visual simultaneous fMRI-EEG experiments. The results point out that it is possible to predict the BOLD fluctuations in occipital cortex by using EEG measurements. (C) 2010 Elsevier Inc. All rights reserved.

Comparison of spatial and temporal pattern for fMRI obtained with BOLD and arterial spin labeling

Functional magnetic resonance imaging (fMRI) is presently either performed using blood oxygenation level-dependent (BOLD) contrast or using cerebral blood flow (CBF), measured with arterial spin labeling (ASL) technique. The present fMRI study aimed to provide practical hints to favour one method over the other. It involved three different acquisition methods during visual checkerboard stimulation on nine healthy subjects: 1) CBF contrast obtained from ASL, 2) BOLD contrast extracted from ASL and 3) BOLD contrast from Echo planar imaging. Previous findings were replicated; i) no differences between the three measurements were found in the location of the activated region; ii) differences were found in the temporal characteristics of the signals and iii) BOLD has significantly higher sensitivity than ASL perfusion. ASL fMRI was favoured when the investigation demands for perfusion and task related signal changes. BOLD fMRI is more suitable in conjunction with fast event related design.

An event-related fMRI study of cognitive deficits in Huntington's disease

The BOLD contrast signal history determined by lagged Unear correlation has a significant contribution to functional connectivity in activation data sets. It has been demonstrated that in resting state fMRI data, the major contribution to synchronous correlation between functionally connected areas arises from low frequency contributions (

BOLD responses to trigeminal nerve stimulation.

The current study investigates a new model of barrel cortex activation using stimulation of the infraorbital branch of the trigeminal nerve. A robust and reproducible activation of the rat barrel cortex was obtained following trigeminal nerve stimulation. Blood oxygen level-dependent (BOLD) effects were obtained in the primary somatosensory barrel cortex (S1BF), the secondary somatosensory cortex (S2) and the motor cortex. These cortical areas were reached from afferent pathways from the trigeminal ganglion, the trigeminal nuclei and thalamic nuclei from which neurons project their axons upon whisker stimulation. The maximum BOLD responses were obtained for a stimulus frequency of 1 Hz, a stimulus pulse width of 100 μs and for current intensities between 1.5 and 3 mA. The BOLD response was nonlinear as a function of frequency and current intensity. Additionally, modeling BOLD responses in the rat barrel cortex from separate cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)) measurements showed good agreement with the shape and amplitude of measured BOLD responses as a function of stimulus frequency and will potentially allow to identify the sources of BOLD nonlinearities. Activation of the rat barrel cortex using trigeminal nerve stimulation will contribute to the interpretation of the BOLD signals from functional magnetic resonance imaging studies.

The history and role of long duration stimulation in fMRI.

During the past 20 years, BOLD fMRI has developed towards a central and fundamental tool in neuroscience. It has been shown that the BOLD response provides an indicator of neuronal activity in the brain. Consequently, for an accurate interpretation of findings in BOLD MRI experiments and to draw meaningful conclusions about the temporal evolution of neural events, a deep understanding of the nature of the BOLD contrast has become of essential importance. Since the dynamics of the major direct determinants of the BOLD signal (CBF, CBV and CMRO(2)) range between seconds and minutes, long duration stimulation was an early key strategy needed to study and understand the BOLD characteristics. This paper summarizes and discusses the thoughts and rationales of the long duration stimulation studies.

Effect of repetitive arm cycling following botulinum toxin injection for poststroke spasticity: evidence from FMRI.

BACKGROUND AND OBJECTIVE: Investigations were performed to establish if repetitive arm cycling training enhances the antispastic effect of intramuscular botulinum toxin (BTX) injections in postischemic spastic hemiparesis. Effects on cerebral activation were evaluated by functional magnetic resonance imaging (fMRI). METHODS: Eight chronic spastic hemisyndrome patients (49 ± 10 years) after middle cerebral artery infarction (5.5 ± 2.7 years) were investigated. BTX was injected into the affected arm twice, 6 months apart. Spasticity was assessed using the Ashworth Scale and range of motion before and 3 months after BTX injections. Images were analyzed using Brain Voyager QX 1.8, and fMRI signal changes were corrected for multiple comparisons. RESULTS: During passive movements of affected and nonaffected hands, fMRI activity was increased bilaterally in the sensorimotor cortex (MISI), secondary somatosensory areas (SII), and supplementary motor area predominantly in the contralesional hemisphere, compared with the rest. Following repetitive arm cycling, fMRI activity increased further in MISI of the lesioned hemisphere and SII of the contralesional hemisphere. For patients with residual motor activity, treatment-related fMRI activity increases were associated with reduced spasticity; in completely plegic patients, there was no fMRI activity change in SII but increased spasticity after training. CONCLUSION: Increased activity in SII of the contralesional hemisphere and in MISI of the lesioned hemisphere reflect a treatment-induced effect in the paretic arm. It is hypothesized that the increased BOLD activity results from increased afferent information related to the antispastic BTX effect reinforced by training.

Modulation of visually evoked cortical fMRI responses by phase of ongoing occipital alpha oscillations

Using simultaneous electroencephalography as a measure of ongoing activity and functional magnetic resonance imaging (fMRI) as a measure of the stimulus-driven neural response, we examined whether the amplitude and phase of occipital alpha oscillations at the onset of a brief visual stimulus affects the amplitude of the visually evoked fMRI response. When accounting for intrinsic coupling of alpha amplitude and occipital fMRI signal by modeling and subtracting pseudo-trials, no significant effect of prestimulus alpha amplitude on the evoked fMRI response could be demonstrated. Regarding the effect of alpha phase, we found that stimuli arriving at the peak of the alpha cycle yielded a lower blood oxygenation level-dependent (BOLD) fMRI response in early visual cortex (V1/V2) than stimuli presented at the trough of the cycle. Our results therefore show that phase of occipital alpha oscillations impacts the overall strength of a visually evoked response, as indexed by the BOLD signal. This observation complements existing evidence that alpha oscillations reflect periodic variations in cortical excitability and suggests that the phase of oscillations in postsynaptic potentials can serve as a mechanism of gain control for incoming neural activity. Finally, our findings provide a putative neural basis for observations of alpha phase dependence of visual perceptual performance.

Ressonância magnética funcional: mapeamento do córtex motor através da técnica BOLD

A ressonância magnética funcional (RMf) é hoje uma ferramenta fundamental na investigação funcional do cérebro humano, quer em indivíduos saudáveis quer em pacientes com patologias diversas. É uma técnica complexa que necessita de uma aplicação cuidada e rigorosa, e uma compreensão dos mecanismos biofísicos a ela subjacentes, de modo a serem obtidos resultados fiáveis e com melhor aceitação clínica. O efeito BOLD (Blood Oxygenation Level Dependent) é o método mais utilizado para medir e estudar a actividade cerebral e baseia-se nas alterações das propriedades magnéticas da molécula hemoglobina. Com este Projecto propomo-nos optimizar um protocolo de RMf realizada com o efeito BOLD, em voluntários saudáveis, de modo a que este possa ser aplicado em futuros estudos de pacientes com patologias. ABSTRACT - Nowadays functional magnetic resonance imaging (fMRI) is a fundamental tool for the research of human brain function of healthy subjects or patients with several pathologies. It is a complex technique that requires a careful and rigorous application, and an understanding of its biophysical mechanisms, so that reliable results can be obtained with better clinical acceptance. The BOLD effect (Blood Oxygenation Level Dependent) is the most widely used method to measure and study the brain activity and its based on changes in magnetic properties of the hemoglobin molecule. The aim of this project was to optimize a BOLD fMRI protocol on healthy subjects, so it can be applied in future studies of patients with pathologies.

Ressonância magnética funcional: mapeamento do córtex motor através do efeito BOLD

Introdução – A ressonância magnética funcional (RMf) é hoje uma ferramenta fundamental na investigação funcional do cérebro humano, quer em indivíduos saudáveis quer em doentes com patologias diversas. É uma técnica complexa que necessita de uma aplicação cuidada e rigorosa e uma compreensão dos mecanismos biofísicos, de modo a serem obtidos resultados fiáveis e com melhor aceitação clínica. O efeito BOLD (Blood Oxygenation Level Dependent), que se baseia nas propriedades magnéticas da hemoglobina, é o método mais utilizado para medir a atividade cerebral por RMf. Objetivos – Otimizar um protocolo de RMf por efeito BOLD em voluntários saudáveis para mapeamento do córtex motor, de modo a que possa ser aplicado no futuro em doentes com patologias diversas. Metodologia – Foram estudados 34 voluntários saudáveis divididos em 2 grupos de estudo: BOLD 1 e BOLD 2. Com vista à otimização, foram testados no subgrupo BOLD 1 diferentes paradigmas e no subgrupo BOLD 2 foi estudada a influência do tempo de eco (TE). Para as várias condições foram comparados os volumes da região ativada e os níveis de ativação obtidos. Resultados/Discussão – O córtex motor foi identificado em todos os voluntários estudados. Não foram detetadas diferenças estatisticamente significativas quando comparados os resultados obtidos com os diferentes parâmetros de aquisição. Conclusão – O protocolo foi otimizado tendo em conta o nível de conforto reportado pelos voluntários. Uma vez que se pretende aplicar este mesmo protocolo no estudo de doentes, este fator torna-se particularmente relevante.

Regional reproducibility of calibrated BOLD functional MRI: implications for the study of cognition and plasticity.

Calibrated BOLD fMRI is a promising alternative to the classic BOLD contrast due to its reduced venous sensitivity and greater physiological specificity. The delayed adoption of this technique for cognitive studies may stem partly from a lack of information on the reproducibility of these measures in the context of cognitive tasks. In this study we have explored the applicability and reproducibility of a state-of-the-art calibrated BOLD technique using a complex functional task at 7 tesla. Reproducibility measures of BOLD, CBF, CMRO2 flow-metabolism coupling n and the calibration parameter M were compared and interpreted for three ROIs. We found an averaged intra-subject variation of CMRO2 of 8% across runs and 33% across days. BOLD (46% across runs, 36% across days), CBF (33% across runs, 46% across days) and M (41% across days) showed significantly higher intra-subject variability. Inter-subject variability was found to be high for all quantities, though CMRO2 was the most consistent across brain regions. The results of this study provide evidence that calibrated BOLD may be a viable alternative for longitudinal and cognitive MRI studies.