On the sensitivity of the diffusion MRI signal to brain activity in response to a motor cortex paradigm

Diffusion functional MRI (dfMRI) is a promising technique to map functional activations by acquiring diffusion-weighed spin-echo images. In previous studies, dfMRI showed higher spatial accuracy at activation mapping compared to classic functional MRI approaches. However, it remains unclear whether dfMRI measures result from changes in the intra-/extracellular environment, perfusion and/or T2 values. We designed an acquisition/quantification scheme to disentangle such effects in the motor cortex during a finger tapping paradigm. dfMRI was acquired at specific diffusion weightings to selectively suppress perfusion and free-water diffusion, then times series of the apparent diffusion coefficient (ADC-fMRI) and of the perfusion signal fraction (IVIM-fMRI) were derived. ADC-fMRI provided ADC estimates sensitive to changes in perfusion and free-water volume, but not to T2/T2* values. With IVIM-fMRI we isolated the perfusion contribution to ADC, while suppressing T2 effects. Compared to conventional gradient-echo BOLD fMRI, activation maps obtained with dfMRI and ADC-fMRI had smaller clusters, and the spatial overlap between the three techniques was below 50%. Increases of perfusion fractions were observed during task in both dfMRI and ADC-fMRI activations. Perfusion effects were more prominent with ADC-fMRI than with dfMRI but were significant in less than 25% of activation ROIs. Taken together, our results suggest that the sensitivity to task of dfMRI derives from a decrease of hindered diffusion and an increase of the pseudo-diffusion signal fraction, leading to different, more confined spatial activation patterns compared to classic functional MRI.Comment: Submitted to peer-reviewed journa

Real-Time Decoding of Brain Responses to Visuospatial Attention Using 7T fMRI

Brain-Computer interface technologies mean to create new communication channels between our mind and our environment, independent of the motor system, by detecting and classifying self regulation of local brain activity. BCIs can provide patients with severe paralysis a means to communicate and to live more independent lives. There has been a growing interest in using invasive recordings for BCI to improve the signal quality. This also potentially gives access to new control strategies previously inaccessible by non-invasive methods. However, before surgery, the best implantation site needs to be determined. The blood-oxygen-level dependent signal changes measured with fMRI have been shown to agree well spatially with those found with invasive electrodes, and are the best option for pre-surgical localization. We show, using real-time fMRI at 7T, that eye movement-independent visuospatial attention can be used as a reliable control strategy for BCIs. At this field strength even subtle signal changes can be detected in single trials thanks to the high contrast-to-noise ratio. A group of healthy subjects were instructed to move their attention between three (two peripheral and one central) spatial target regions while keeping their gaze fixated at the center. The activated regions were first located and thereafter the subjects were given real-time feedback based on the activity in these regions. All subjects managed to regulate local brain areas without training, which suggests that visuospatial attention is a promising new target for intracranial BCI. ECoG data recorded from one epilepsy patient showed that local changes in gamma-power can be used to separate the three classes

Shape and volume changes of the superior lateral ventricle after electroconvulsive therapy measured with ultra-high field MRI

The subventricular zone (SVZ) of the lateral ventricles harbors neuronal stem cells in adult mammals. Rodent studies report neurogenic effects in the SVZ of electroconvulsive stimulation. We hypothesize that if this finding translates to depressed patients undergoing electroconvulsive therapy (ECT), this would be reflected in shape changes at the SVZ. Using T1-weighted MR images acquired at ultra-high field strength (7T), the shape and volume of the ventricles were compared from pre to post ECT after 10 ECT sessions (in patients twice weekly) or 5 weeks apart (controls) using linear mixed models with age and gender as covariates. Ventricle shape significantly changed and volume significantly decreased over time in patients for the left ventricle, but not in controls. The decrease in volume of the ventricles was associated to a decrease in depression scores, and an increase in the left dentate gyrus, However, the shape changes of the ventricles were not restricted to the neurogenic niche in the lateral walls of the ventricles, providing no clear evidence for neurogenesis as sole explanation of volume changes in the ventricles after ECT

CADASIL Affects Multiple Aspects of Cerebral Small Vessel Function on 7T-MRI

International audienceObjective: Cerebral small vessel diseases (cSVDs) are a major cause of stroke and dementia. We used cutting-edge 7T-MRI techniques in patients with Cerebral Autosomal Dominant Arteriopathy with Subcortical Infarcts and Leukoencephalopathy (CADASIL), to establish which aspects of cerebral small vessel function are affected by this monogenic form of cSVD. Methods: We recruited 23 CADASIL patients (age 51.1 AE 10.1 years, 52% women) and 13 age-and sex-matched controls (46.1 AE 12.6, 46% women). Small vessel function measures included: basal ganglia and centrum semiovale perforating artery blood flow velocity and pulsatility, vascular reactivity to a visual stimulus in the occipital cortex and reactivity to hypercapnia in the cortex, subcortical gray matter, white matter, and white matter hyperintensities. Results: Compared with controls, CADASIL patients showed lower blood flow velocity and higher pulsatility index within perforating arteries of the centrum semiovale (mean difference À 0.09 cm/s, p = 0.03 and 0.20, p = 0.009) and basal ganglia (mean difference À 0.98 cm/s, p = 0.003 and 0.17, p = 0.06). Small vessel reactivity to a short visual stimulus was decreased (blood-oxygen-level dependent [BOLD] mean difference À0.21%, p = 0.04) in patients, while reactivity to hypercapnia was preserved in the cortex, subcortical gray matter, and normal appearing white matter. Among patients, reactivity to hypercapnia was decreased in white matter hyperintensities compared to normal appearing white matter (BOLD mean difference À0.29%, p = 0.02). Interpretation: Multiple aspects of cerebral small vessel function on 7T-MRI were abnormal in CADASIL patients, indicative of increased arteriolar stiffness and regional abnormalities in reactivity, locally also in relation to white matter injury. These observations provide novel markers of cSVD for mechanistic and intervention studies

A resting-state fMRI pattern of spinocerebellar ataxia type 3 and comparison with F-18-FDG PET

Spinocerebellar ataxia type 3 (SCA3) is a rare genetic neurodegenerative disease. The neurobiological basis of SCA3 is still poorly understood, and up until now resting-state fMRI (rs-fMRI) has not been used to study this disease. In the current study we investigated (multi-echo) rs-fMRI data from patients with genetically confirmed SCA3 (n = 17) and matched healthy subjects (n = 16). Using independent component analysis (ICA) and subsequent regression with bootstrap resampling, we identified a pattern of differences between patients and healthy subjects, which we coined the fMRI SCA3 related pattern (fSCA3-RP) comprising cerebellum, anterior striatum and various cortical regions. Individual fSCA3-RP scores were highly correlated with a previously published F-18-FDG PET pattern found in the same sample (rho = 0.78, P = 0.0003). Also, a high correlation was found with the Scale for Assessment and Rating of Ataxia scores (r = 0.63, P = 0.007). No correlations were found with neuropsychological test scores, nor with levels of grey matter atrophy. Compared with the F-18-FDG PET pattern, the fSCA3-RP included a more extensive contribution of the mediofrontal cortex, putatively representing changes in default network activity. This rs-fMRI identification of additional regions is proposed to reflect a consequence of the nature of the BOLD technique, enabling measurement of dynamic network activity, compared to the more static F-18-FDG PET methodology. Altogether, our findings shed new light on the neural substrate of SCA3, and encourage further validation of the fSCA3-RP to assess its potential contribution as imaging biomarker for future research and clinical use

Zooming in on cerebral small vessel function in small vessel diseases with 7T MRI: Rationale and design of the “ZOOM@SVDs” study

Background: Cerebral small vessel diseases (SVDs) are a major cause of stroke and dementia. Yet, specific treatment strategies are lacking in part because of a limited understanding of the underlying disease processes. There is therefore an urgent need to study SVDs at their core, the small vessels themselves. Objective: This paper presents the rationale and design of the ZOOM@SVDs study, which aims to establish measures of cerebral small vessel dysfunction on 7T MRI as novel disease markers of SVDs. Methods: ZOOM@SVDs is a prospective observational cohort study with two years follow-up. ZOOM@SVDs recruits participants with Cerebral Autosomal Dominant Arteriopathy with Subcortical Infarcts and Leukoencephalopathy (CADASIL, N = 20), sporadic SVDs (N = 60), and healthy controls (N = 40). Participants undergo 7T brain MRI to assess different aspects of small vessel function including small vessel reactivity, cerebral perforating artery flow, and pulsatility. Extensive work-up at baseline and follow-up further includes clinical and neuropsychological assessment as well as 3T brain MRI to assess conventional SVD imaging markers. Measures of small vessel dysfunction are compared between patients and controls, and related to the severity of clinical and conventional MRI manifestations of SVDs. Discussion: ZOOM@SVDs will deliver novel markers of cerebral small vessel function in patients with monogenic and sporadic forms of SVDs, and establish their relation with disease burden and progression. These small vessel markers can support etiological studies in SVDs and may serve as surrogate outcome measures in future clinical trials to show target engagement of drugs directed at the small vessels

Neuronal Models for EEG–fMRI Integration

Human brain activity can be measured in many ways, providing different views into its functions. Common measurements of human brain function, such as functional MRI (fMRI) and electroencephalography (EEG), are often integrated to get the best spatial and temporal resolution, but these measurements frequently show differences. This chapter therefore considers how the pooling of signals across populations of neurons affects these measurements. We consider a modeling framework in which fMRI and field potential signals integrate across transmembrane currents in a population of neurons, because synaptic activity is thought to be the largest drive of the fMRI signal and transmembrane potentials are the biggest source of field potentials. We formulate computations that provide simplified abstractions that inform how levels of synaptic activity or synchrony across neurons contribute to fMRI or electrophysiological signals and how to interpret deviations or similarities between measurements. The modeling framework and computations highlight that the level of activity in a neuronal population influences each measurement, but synchrony only has a large effect on field potentials and not on the fMRI signal. An application to data from human visual cortex explains why certain signals correlate and others do not. Advancing the fundamental understanding of how different measurements integrate neuronal activity will be important to combine fMRI and EEG measurements to better understand human brain function

Laminar fMRI: What can the time domain tell us?

The rapid developments in functional MRI (fMRI) acquisition methods and hardware technologies in recent years, particularly at high field (≥7 T), have enabled unparalleled visualization of functional detail at a laminar or columnar level, bringing fMRI close to the intrinsic resolution of brain function. These advances highlight the potential of high resolution fMRI to be a valuable tool to study the fundamental processing performed in cortical micro-circuits, and their interactions such as feedforward and feedback processes. Notably, because fMRI measures neuronal activity via hemodynamics, the ultimate resolution it affords depends on the spatial specificity of hemodynamics to neuronal activity at a detailed spatial scale, and by the evolution of this specificity over time. Several laminar (≤1 mm spatial resolution) fMRI studies have examined spatial characteristics of the measured hemodynamic signals across cortical depth, in light of understanding or improving the spatial specificity of laminar fMRI. Few studies have examined temporal features of the hemodynamic response across cortical depth. Temporal features of the hemodynamic response offer an additional means to improve the specificity of fMRI, and could help target neuronal processes and neurovascular coupling relationships across laminae, for example by differences in the onset times of the response across cortical depth. In this review, we discuss factors that affect the timing of neuronal and hemodynamic responses across laminae, touching on the neuronal laminar organization, and focusing on the laminar vascular organization. We provide an overview of hemodynamics across the cortical vascular tree based on optical imaging studies, and review temporal aspects of hemodynamics that have been examined across cortical depth in high spatiotemporal resolution fMRI studies. Last, we discuss the limits and potential of high spatiotemporal resolution fMRI to study laminar neurovascular coupling and neuronal processes

Effect sizes of BOLD CVR, resting-state signal fluctuations and time delay measures for the assessment of hemodynamic impairment in carotid occlusion patients

Background and purpose: The BOLD signal amplitude as a response to a hypercapnia stimulus is commonly used to assess cerebrovascular reserve. Despite recent advances, the implementation remains cumbersome and alternative ways to assess hemodynamic impairment are desirable. Resting-state BOLD signal fluctuations (rsBOLD) have been proposed however data on its sensitivity and dependence on baseline venous cerebral blood volume (vCBV) is limited. The primary aim of this study was to compare the effect sizes of resting-state and hypercapnia induced BOLD signal changes in the detection of hemodynamic impairment. The second aim of the study was to assess the dependence of BOLD signal variability on vCBV. Materials and methods: Fifteen patients with internal carotid artery occlusive disease and 15 matched healthy controls were included in this study. The BOLD signal was derived from a dual-echo gradient-echo echo-planar sequence during hypercapnia (HC) and hyperoxia (HO) gas modulations. BOLD (fractional) amplitude of low frequency fluctuations ((f)ALFF) was compared to HC-BOLD, BOLD response delays derived from time delay analysis and ΔBOLD in response to progressively increasing HC. Effect sizes (i.e. the standard mean difference between patients and controls) were calculated. HO-BOLD was used to estimate vCBV, and its contribution to the variability in rsBOLD signal was evaluated. Results: The effect sizes of ALFF and fALFF (0.61 and 0.72) were lower than the effect sizes related to hypercapnia-based hemodynamic assessment analysis; 1.62, 1.56 and 0.90 for HC-BOLD, BOLD response delays and ΔBOLD in response to progressively increasing HC. A moderate relation was found between (f)ALFF and HC-BOLD in controls (R2 of 0.61 and 0.42), but this relation decreased in patients (R2 of 0.33 and 0.15). (f)ALFF did not differ between patients and controls whereas HC-BOLD did (p < 0.005). The ΔBOLD response to progressively increasing HC was significantly different in between patients and controls for ΔEtCO2 values ≥ 2 mmHg (at +2 mmHg F(1, 18) = 5.85, p = 0.026). Up to 31% and 53% of the variance in the ALFF and HC-BOLD spatial distribution could be explained by HO-BOLD. Conclusion: ALFF and fALFF demonstrated a moderate effect size to detect hemodynamic impairment whereas the effect size was large for methods employing a hypercapnia-based vascular stress stimulus. Based on our analysis of BOLD signal change as a response to a progressively increasing hypercapnia stimulus we can argue that a hypercapnia stimulus of at least 2 mmHg above baseline EtCO2 is necessary to evaluate hemodynamic impairment. We also demonstrated that a substantial amount of information imbedded in the rsBOLD and HC-BOLD was explained by HO-BOLD. HO-BOLD can serve as a proxy for vCBV and this thus indicates that one should be careful when adopting these techniques in disease cases with compromised CBV

Accelerating Brain Imaging Using a Silent Spatial Encoding Axis

PURPOSE: To characterize the acceleration capabilities of a silent head insert gradient axis that operates at the inaudible frequency of 20 kHz and a maximum gradient amplitude of 40 mT/m without inducing peripheral nerve stimulation. METHODS: The silent gradient axis' acquisitions feature an oscillating gradient in the phase-encoding direction that is played out on top of a cartesian readout, similarly as done in Wave-CAIPI. The additional spatial encoding fills k-space in readout lanes allowing for the acquisition of fewer phase-encoding steps without increasing aliasing artifacts. Fully sampled 2D gradient echo datasets were acquired both with and without the silent readout. All scans were retrospectively undersampled (acceleration factors R = 1 to 12) to compare conventional SENSE acceleration and acceleration using the silent gradient. The silent gradient amplitude and the readout bandwidth were varied to investigate the effect on artifacts and g-factor. RESULTS: The silent readout reduced the g-factor for all acceleration factors when compared to SENSE acceleration. Increasing the silent gradient amplitude from 31.5 mT/m to 40 mT/m at an acceleration factor of 10 yielded a reduction in the average g-factor (gavg ) from 1.3 ± 0.14 (gmax = 1.9) to 1.1 ± 0.09 (gmax = 1.6). Furthermore, reducing the number of cycles increased the readout bandwidth and the g-factor that reached gavg = 1.5 ± 0.16 for a readout bandwidth of 651 Hz/pixel and an acceleration factor of R = 8. CONCLUSION: A silent gradient axis enables high acceleration factors up to R = 10 while maintaining a g-factor close to unity (gavg = 1.1 and gmax = 1.6) and can be acquired with clinically relevant readout bandwidths