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

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

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

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

Fast CSF MRI for brain segmentation; Cross-validation by comparison with 3D T1-based brain segmentation methods

OBJECTIVE: In previous work we have developed a fast sequence that focusses on cerebrospinal fluid (CSF) based on the long T2 of CSF. By processing the data obtained with this CSF MRI sequence, brain parenchymal volume (BPV) and intracranial volume (ICV) can be automatically obtained. The aim of this study was to assess the precision of the BPV and ICV measurements of the CSF MRI sequence and to validate the CSF MRI sequence by comparison with 3D T1-based brain segmentation methods. MATERIALS AND METHODS: Ten healthy volunteers (2 females; median age 28 years) were scanned (3T MRI) twice with repositioning in between. The scan protocol consisted of a low resolution (LR) CSF sequence (0:57min), a high resolution (HR) CSF sequence (3:21min) and a 3D T1-weighted sequence (6:47min). Data of the HR 3D-T1-weighted images were downsampled to obtain LR T1-weighted images (reconstructed imaging time: 1:59 min). Data of the CSF MRI sequences was automatically segmented using in-house software. The 3D T1-weighted images were segmented using FSL (5.0), SPM12 and FreeSurfer (5.3.0). RESULTS: The mean absolute differences for BPV and ICV between the first and second scan for CSF LR (BPV/ICV: 12±9/7±4cc) and CSF HR (5±5/4±2cc) were comparable to FSL HR (9±11/19±23cc), FSL LR (7±4, 6±5cc), FreeSurfer HR (5±3/14±8cc), FreeSurfer LR (9±8, 12±10cc), and SPM HR (5±3/4±7cc), and SPM LR (5±4, 5±3cc). The correlation between the measured volumes of the CSF sequences and that measured by FSL, FreeSurfer and SPM HR and LR was very good (all Pearson's correlation coefficients >0.83, R2 .67-.97). The results from the downsampled data and the high-resolution data were similar. CONCLUSION: Both CSF MRI sequences have a precision comparable to, and a very good correlation with established 3D T1-based automated segmentations methods for the segmentation of BPV and ICV. However, the short imaging time of the fast CSF MRI sequence is superior to the 3D T1 sequence on which segmentation with established methods is performed

Phase contrast MRI measurements of net cerebrospinal fluid flow through the cerebral aqueduct are confounded by respiration

BACKGROUND: Net cerebrospinal fluid (CSF) flow through the cerebral aqueduct may serve as a marker of CSF production in the lateral ventricles, and changes that occur with aging and in disease. PURPOSE: To investigate the confounding influence of the respiratory cycle on net CSF flow and stroke volume measurements. STUDY TYPE: Cross-sectional study. SUBJECTS: Twelve young, healthy subjects (seven male, age range 19-39 years, average age 28.3 years). FIELD STRENGTH/SEQUENCE: Phase contrast MRI (PC-MRI) measurements were performed at 7T, with and without respiratory gating on expiration and on inspiration. All measurements were repeated. ASSESSMENT: Net CSF flow and stroke volume in the aqueduct, over the cardiac cycle, was determined. STATISTICAL TESTS: Repeatability was determined using the intraclass correlation coefficient (ICC) and linear regression analysis between the repeated measurements. Repeated measures analysis of variance (ANOVA) was performed to compare the measurements during inspiration/expiration/no gating. Linear regression analysis was performed between the net CSF flow difference (inspiration minus expiration) and stroke volume. RESULTS: Net CSF flow (average ± standard deviation) was 0.64 ± 0.32 mL/min (caudal) during expiration, 0.12 ± 0.49 mL/min (cranial) during inspiration, and 0.31 ± 0.18 mL/min (caudal) without respiratory gating. Respiratory gating did not affect stroke volume measurements (41 ± 18, 42 ± 19, 42 ± 19 μL/cycle for expiration, no respiratory gating and inspiration, respectively). Repeatability was best during inspiration (ICC = 0.88/0.56/-0.31 for gating on inspiration/expiration/no gating). A positive association was found between average stroke volume and net flow difference between inspiration and expiration (R = 0.678/0.605, P = 0.015/0.037 for the first/second repeated measurement). DATA CONCLUSION: Measured net CSF flow is confounded by respiration effects. Therefore, net CSF flow measurements with PC-MRI cannot in isolation be directly linked to CSF production. LEVEL OF EVIDENCE: 1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018

Improved Selectivity in 7 T Digit Mapping Using VASO-CBV

Functional magnetic resonance imaging (fMRI) at Ultra-high field (UHF, ≥ 7 T) benefits from significant gains in the BOLD contrast-to-noise ratio (CNR) and temporal signal-to-noise ratio (tSNR) compared to conventional field strengths (3 T). Although these improvements enabled researchers to study the human brain to unprecedented spatial resolution, the blood pooling effect reduces the spatial specificity of the widely-used gradient-echo BOLD acquisitions. In this context, vascular space occupancy (VASO-CBV) imaging may be advantageous since it is proposed to have a higher spatial specificity than BOLD. We hypothesized that the assumed higher specificity of VASO-CBV imaging would translate to reduced overlap in fine-scale digit representation maps compared to BOLD-based digit maps. We used sub-millimeter resolution VASO fMRI at 7 T to map VASO-CBV and BOLD responses simultaneously in the motor and somatosensory cortices during individual finger movement tasks. We assessed the cortical overlap in different ways, first by calculating similarity coefficient metrics (DICE and Jaccard) and second by calculating selectivity measures. In addition, we demonstrate a consistent topographical organization of the targeted digit representations (thumb-index-little finger) in the motor areas. We show that the VASO-CBV responses yielded less overlap between the digit clusters than BOLD, and other selectivity measures were higher for VASO-CBV too. In summary, these results were consistent across metrics and participants, confirming the higher spatial specificity of VASO-CBV compared to BOLD

BOLD specificity and dynamics evaluated in humans at 7 T: comparing gradient-echo and spin-echo hemodynamic responses.

High-field gradient-echo (GE) BOLD fMRI enables very high resolution imaging, and has great potential for detailed investigations of brain function. However, as spatial resolution increases, confounds due to signal from non-capillary vessels increasingly impact the fidelity of GE BOLD fMRI signals. Here we report on an assessment of the microvascular weighting of the GE BOLD response across the cortical depth in human cortex using spin-echo fMRI which is thought to be dominated by microvasculature (albeit not completely). BOLD responses were measured with a hemodynamic impulse response (HRF) obtained from the spin-echo (SE) and gradient-echo (GE) BOLD contrast using very short stimuli (0.25 s) and a fast event-related functional paradigm. We show that the onset (≈ 1.25 s) and the rising slope of the GE and SE HRFs are strikingly similar for voxels in deep gray matter presumably containing the most metabolically demanding neurons (layers III-IV). This finding provides a strong indication that the onset of the GE HRF in deep gray matter is predominantly associated with microvasculature