Deconvolution filtering: Temporal smoothing revisited

Inferences made from analysis of BOLD data regarding neural processes are potentially confounded by multiple competing sources: cardiac and respiratory signals, thermal effects, scanner drift, and motion-induced signal intensity changes. To address this problem, we propose deconvolution filtering, a process of systematically deconvolving and reconvolving the BOLD signal via the hemodynamic response function such that the resultant signal is composed of maximally likely neural and neurovascular signals. To test the validity of this approach, we compared the accuracy of BOLD signal variants (i.e., unfiltered, deconvolution filtered, band-pass filtered, and optimized band-pass filtered BOLD signals) in identifying useful properties of highly confounded, simulated BOLD data: (1) reconstructing the true, unconfounded BOLD signal, (2) correlation with the true, unconfounded BOLD signal, and (3) reconstructing the true functional connectivity of a three-node neural system. We also tested this approach by detecting task activation in BOLD data recorded from healthy adolescent girls (control) during an emotion processing task.     

Fine-scale functional connectivity in somatosensory cortex revealed by high-resolution fMRI

High-resolution functional magnetic resonance imaging (fMRI) at high field (9.4 T) has been used to measure functional connectivity between subregions within the primary somatosensory (SI) cortex of the squirrel monkey brain. The hand-face region within the SI cortex of the squirrel monkey has been previously well mapped with functional imaging and electrophysiological and anatomical methods, and the orderly topographic map of the hand region is characterized by a lateral to medial representation of individual digits in four subregions of areas 3a, 3b, 1 and 2. With submillimeter resolution, we are able to detect not only the separate islands of activation corresponding to vibrotactile stimulations of single digits but also, in subsequent acquisitions, the degree of correlation between voxels within the SI cortex in the resting state. The results suggest that connectivity patterns are very similar to stimulus-driven distributions of activity and that connectivity varies on the scale of millimeters within the same primary region. Connectivity strength is not a reflection of global larger-scale changes in blood flow and is not directly dependent on distance between regions. Preliminary electrophysiological recordings agree well with the fMRI data. In human studies at 7 T, high-resolution fMRI may also be used to identify the same subregions and assess responses to sensory as well as painful stimuli, and to measure connectivity dynamically before and after such stimulations.     

Liver functional magnetic resonance imaging analysis using a latent variables approach

The liver is a highly vascular organ with a dual blood supply, and it performs a remarkable number of vital functions. Here, we show, through measurement of blood oxygen level‐dependent (BOLD) signal, that liver arterial and hepatic portal blood supplies can be modulated through hyperoxia exposure and by consumption of a standardized meal, respectively. As such, we suggest that hyperoxia modulates the hepatic arterial BOLD signal, whereas a controlled meal changes predominantly the hepatic portal BOLD signal. The hemodynamics of the dual liver blood supplies in response to the aforementioned challenges are complex and variable across subjects, making a general linear model‐based analysis difficult. Therefore, we present the application of two local (at each voxel) hemodynamic response‐independent techniques—principal component analysis and partial least squares—to observe the hypothesized reduction in BOLD contrast during cycles of hyperoxic breathing, when comparing preprandial versus postprandial states in a normally functioning liver. We illustrate the ability of our techniques to differentiate between healthy and diseased livers with an analysis of 17 subjects—11 with normal livers and 6 with liver disease (hepatitis or cirrhosis). Our local analysis can correctly classify all of the subjects. Copyright © 2012 John Wiley & Sons, Ltd.     

Responses of dopaminergic, serotonergic and noradrenergic networks to acute levo-tetrahydropalmatine administration in naïve rats detected at 9.4 T

Aim

The aim of this study was to understand the neuropharmacological characteristics of levo-tetrahydropalmatine (l-THP), a recently found potential treatment for drug addiction, and discover its neural correlates and sites of action.

Methods

High-field pharmacological magnetic resonance imaging (phMRI) was used to detect activation induced by acute l-THP administration in the naïve rat brain at dose levels of 5, 10, 20 and 40 mg/kg.

Results

Interestingly, the pharmacological profile of l-THP selectively binds to the receptors of the dopaminergic, serotonergic and noradrenergic systems. Using the phMRI method, it was demonstrated that l-THP selectively activated the key brain regions of the dopaminergic, serotonergic and noradrenergic systems in a dose-dependent manner.

Conclusion

Numerous studies suggest a critical role of monoamines in the behavioral, pharmacological and addictive properties of psychostimulants. It is suggested that l-THP holds great potential to be a therapeutic medication for drug addiction.     

SNR and functional sensitivity of BOLD and perfusion-based fMRI using arterial spin labeling with spiral SENSE at 3 T

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) studies using parallel imaging to reduce the readout window have reported a loss in temporal signal-to-noise ratio (SNR) that is less than would be expected given a purely thermal noise model. In this study, the impact of parallel imaging on the noise components and functional sensitivity of both BOLD and perfusion-based fMRI data was investigated. Dual-echo arterial spin labeling data were acquired on five subjects using sensitivity encoding (SENSE), at reduction factors (R) of 1, 2 and 3. Direct recording of cardiac and respiratory activity during data acquisition enabled the retrospective removal of physiological noise. The temporal SNR of the perfusion time series closely followed the thermal noise prediction of a √R loss in SNR as the readout window was shortened, with temporal SNR values (relative to the R=1 data) of 0.72 and 0.56 for the R=2 and R=3 data, respectively, after accounting for physiological noise. However, the BOLD temporal SNR decreased more slowly than predicted even after accounting for physiological noise, with relative temporal SNR values of 0.80 and 0.63 for the R=2 and R=3 data, respectively. Spectral analysis revealed that the BOLD trends were dominated by low-frequency fluctuations, which were not dominant in the perfusion data due to signal processing differences. The functional sensitivity, assessed using mean F values over activated regions of interest (ROIs), followed the temporal SNR trends for the BOLD data. However, results for the perfusion data were more dependent on the threshold used for ROI selection, most likely due to the inherently low SNR of functional perfusion data.     

Post-stimulus response in hemodynamics observed by functional magnetic resonance imaging--difference between the primary sensorimotor area and the supplementary motor area

The blood oxygen level dependency (BOLD) contrast is a useful tool for functional neuroimaging based on the hemodynamic response to neuronal activation. We observed different hemodynamic responses in the BOLD signal between the primary sensorimotor area (SM1) and the supplementary motor area (SMA) in the sequential finger movement task. In the SMA, a stronger initial overshoot and a post-stimulus overshoot were observed. It was hypothesized from the time course analysis that the stronger initial overshoot reflected the activation of the SMA for motor control programming in the initial phase. Although the post-stimulus overshoot may be partially explained by cerebral blood flow (CBF) cerebral blood volume (CBV) uncoupling, its mechanism remained unknown. In the SM1, only the initial overshoot was observed and the level of BOLD signal was almost constant after the initial overshoot during the task period. These observations suggested that the BOLD signal is characterized by both CBF-CBV uncoupling and the neuronal activation characteristics in each region.     

The effects of caffeine ingestion on cortical areas: Functional imaging study

The effect of caffeine as a cognitive enhancer is well known; however, caffeine-induced changes in the cortical regions are still not very clear. Therefore, in this study, we conducted an investigation of the activation and deactivation with blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) and of metabolic activity change with positron emission tomography (PET) in the human brain.     

Modelling of calcium dynamics in brain energy metabolism and Alzheimer's disease

Functional imaging techniques play a major role in the study of brain activation by monitoring the changes in blood flow and energy metabolism. In order to interpret functional neuroimaging data better, the existing mathematical models describing the links that may exist between electrical activity, energy metabolism and hemodynamics in literature are thoroughly analyzed for their advantages and disadvantages in terms of their prediction of available experimental data. Then, these models are combined within a single model that includes membrane ionic currents, glycolysis, mitochondrial activity, exchanges through the blood-brain barrier, as well as brain hemodynamics. Particular attention is paid to the transport and storage of calcium ions in neurons since calcium is not only an important molecule for signalling in neurons, but it is also essential for memory storage. Multiple efforts have underlined the importance of calcium dependent cellular processes in the biochemical characterization of Alzheimer's disease (AD), suggesting that abnormalities in calcium homeostasis might be involved in the pathophysiology of the disease. The ultimate goal of this study is to investigate the hypotheses about the physiological or biochemical changes in health and disease and to correlate them to measurable physiological parameters obtained from functional neuroimaging data as in the time course of blood oxygenation level dependent (BOLD) signal. When calcium dynamics are included in the model, both BOLD signal and metabolite concentration profiles are shown to exhibit temporal behaviour consistent with the experimental data found in literature. In the case of Alzheimer's disease, the effect of halved cerebral blood flow increase results in a negative BOLD signal implying suppressed neural activity.     

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.     

Testing methodologies for the nonlinear analysis of causal relationships in neurovascular coupling

We investigated the use and implementation of a nonlinear methodology for establishing which changes in neurophysiological signals cause changes in the blood oxygenation level-dependent (BOLD) contrast measured in functional magnetic resonance imaging. Unlike previous analytical approaches, which used linear correlation to establish covariations between neural activity and BOLD, we propose a directed information-theoretic measure, the transfer entropy, which can elucidate even highly nonlinear causal relationships between neural activity and BOLD signal. In this study we investigated the practicality of such an analysis given the limited data samples that can be collected experimentally due to the low temporal resolution of BOLD signals. We implemented several algorithms for the estimation of transfer entropy and we tested their effectiveness using simulated local field potentials (LFPs) and BOLD data constructed to match the main statistical properties of real LFP and BOLD signals measured simultaneously in monkey primary visual cortex. We found that using the advanced methods of entropy estimation implemented and described here, a transfer entropy analysis of neurovascular coupling based on experimentally attainable data sets is feasible.