Real-time cardiac synchronization with fixed volume frame rate for reducing physiological instabilities in 3D FMRI

A characterization of cardiac-induced noise in R2 * maps of the brain

PURPOSE

Cardiac pulsation increases the noise level in brain maps of the transverse relaxation rate R2 *. Cardiac-induced noise is challenging to mitigate during the acquisition of R2 * mapping data because its characteristics are unknown. In this work, we aim to characterize cardiac-induced noise in brain maps of the MRI parameter R2 *.

METHODS

We designed a sampling strategy to acquire multi-echo 3D data in 12 intervals of the cardiac cycle, monitored with a fingertip pulse-oximeter. We measured the amplitude of cardiac-induced noise in this data and assessed the effect of cardiac pulsation on R2 * maps computed across echoes. The area of k-space that contains most of the cardiac-induced noise in R2 * maps was then identified. Based on these characteristics, we introduced a tentative sampling strategy that aims to mitigate cardiac-induced noise in R2 * maps of the brain.

RESULTS

In inferior brain regions, cardiac pulsation accounts for R2 * variations of up to 3 s-1 across the cardiac cycle (i.e., ∼35% of the overall variability). Cardiac-induced fluctuations occur throughout the cardiac cycle, with a reduced intensity during the first quarter of the cycle. A total of 50% to 60% of the overall cardiac-induced noise is localized near the k-space center (k < 0.074 mm-1 ). The tentative cardiac noise mitigation strategy reduced the variability of R2 * maps across repetitions by 11% in the brainstem and 6% across the whole brain.

CONCLUSION

We provide a characterization of cardiac-induced noise in brain R2 * maps that can be used as a basis for the design of mitigation strategies during data acquisition.

Improving robustness of 3D multi-shot EPI by structured low-rank reconstruction of segmented CAIPI sampling for fMRI at 7T

Three-dimensional (3D) encoding methods are increasingly being explored as alternatives to two-dimensional (2D) multi-slice acquisitions in fMRI, particularly in cases where high isotropic resolution is needed. 3D multi-shot EPI acquisition, as the workhorse of 3D fMRI imaging, is susceptible to physiological fluctuations which can induce inter-shot phase variations, and thus reducing the achievable tSNR, negating some of the benefit of 3D encoding. This issue can be particularly problematic at ultra-high fields like 7T, which have more severe off-resonance effects. In this work, we aim to improve the temporal stability of 3D multi-shot EPI at 7T by improving its robustness to inter-shot phase variations. We presented a 3D segmented CAIPI sampling trajectory ("seg-CAIPI") and an improved reconstruction method based on Hankel structured low-rank matrix recovery. Simulation and in-vivo results demonstrate that the combination of the seg-CAIPI sampling scheme and the proposed structured low-rank reconstruction is a promising way to effectively reduce the unwanted temporal variance induced by inter-shot physiological fluctuations, and thus improve the robustness of 3D multi-shot EPI for fMRI.

Cerebrovascular blood oxygenation level dependent pulsatility at baseline and following acute exercise among healthy adolescents

Arterial stiffness is linked to cerebral small vessel damage and neurodegeneration, but barriers to accessing deep cerebrovascular anatomy limit our ability to assess the brain. This study describes an adaptation of a cardiac-related scrubbing method as a means of generating blood oxygenation level-dependent pulsatility maps based on the cardiac cycle. We examine BOLD pulsatility at rest, based on the non-parametric deviation from null metric, as well as changes following acute physiological stress from 20 min of moderate-intensity cycling in 45 healthy adolescents. We evaluate the influence of repetition time (TR) and echo time (TE) using simulated and multi-echo empirical data, respectively. There were tissue-specific and voxel-wise BOLD pulsatility decreases 20 min following exercise cessation. BOLD pulsatility detection was comparable over a range of TR and TE values when scan volumes were kept constant; however, short TRs (≤500 ms) and TEs (∼14 ms) acquisitions would yield the most efficient detection. Results suggest cardiac-related BOLD pulsatility may represent a robust and easily adopted method of mapping cerebrovascular pulsatility with voxel-wise resolution.

Methods for cleaning the BOLD fMRI signal

Blood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.

Optimization of 4D vessel-selective arterial spin labeling angiography using balanced steady-state free precession and vessel-encoding

Vessel-selective dynamic angiograms provide a wealth of useful information about the anatomical and functional status of arteries, including information about collateral flow and blood supply to lesions. Conventional x-ray techniques are invasive and carry some risks to the patient, so non-invasive alternatives are desirable. Previously, non-contrast dynamic MRI angiograms based on arterial spin labeling (ASL) have been demonstrated using both spoiled gradient echo (SPGR) and balanced steady-state free precession (bSSFP) readout modules, but no direct comparison has been made, and bSSFP optimization over a long readout period has not been fully explored. In this study bSSFP and SPGR are theoretically and experimentally compared for dynamic ASL angiography. Unlike SPGR, bSSFP was found to have a very low ASL signal attenuation rate, even when a relatively large flip angle and short repetition time were used, leading to a threefold improvement in the measured signal-to-noise ratio (SNR) efficiency compared with SPGR. For vessel-selective applications, SNR efficiency can be further improved over single-artery labeling methods by using a vessel-encoded pseudo-continuous ASL (VEPCASL) approach. The combination of a VEPCASL preparation with a time-resolved bSSFP readout allowed the generation of four-dimensional (4D; time-resolved three-dimensional, 3D) vessel-selective cerebral angiograms in healthy volunteers with 59 ms temporal resolution. Good quality 4D angiograms were obtained in all subjects, providing comparable structural information to 3D time-of-flight images, as well as dynamic information and vessel selectivity, which was shown to be high. A rapid 1.5 min dynamic two-dimensional version of the sequence yielded similar image features and would be suitable for a busy clinical protocol. Preliminary experiments with bSSFP that included the extracranial vessels showed signal loss in regions of poor magnetic field homogeneity. However, for intracranial vessel-selective angiography, the proposed bSSFP VEPCASL sequence is highly SNR efficient and could provide useful information in a range of cerebrovascular diseases. © 2016 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

Evaluation of 2D multiband EPI imaging for high-resolution, whole-brain, task-based fMRI studies at 3T: Sensitivity and slice leakage artifacts

Functional magnetic resonance imaging (fMRI) studies that require high-resolution whole-brain coverage have long scan times that are primarily driven by the large number of thin slices acquired. Two-dimensional multiband echo-planar imaging (EPI) sequences accelerate the data acquisition along the slice direction and therefore represent an attractive approach to such studies by improving the temporal resolution without sacrificing spatial resolution. In this work, a 2D multiband EPI sequence was optimized for 1.5 mm isotropic whole-brain acquisitions at 3 T with 10 healthy volunteers imaged while performing simultaneous visual and motor tasks. The performance of the sequence was evaluated in terms of BOLD sensitivity and false-positive activation at multiband (MB) factors of 1, 2, 4, and 6, combined with in-plane GRAPPA acceleration of 2 × (GRAPPA 2), and the two reconstruction approaches of Slice-GRAPPA and Split Slice-GRAPPA. Sensitivity results demonstrate significant gains in temporal signal-to-noise ratio (tSNR) and t-score statistics for MB 2, 4, and 6 compared to MB 1. The MB factor for optimal sensitivity varied depending on anatomical location and reconstruction method. When using Slice-GRAPPA reconstruction, evidence of false-positive activation due to signal leakage between simultaneously excited slices was seen in one instance, 35 instances, and 70 instances over the ten volunteers for the respective accelerations of MB 2 × GRAPPA 2, MB 4 × GRAPPA 2, and MB 6 × GRAPPA 2. The use of Split Slice-GRAPPA reconstruction suppressed the prevalence of false positives significantly, to 1 instance, 5 instances, and 5 instances for the same respective acceleration factors. Imaging protocols using an acceleration factor of MB 2 × GRAPPA 2 can be confidently used for high-resolution whole-brain imaging to improve BOLD sensitivity with very low probability for false-positive activation due to slice leakage. Imaging protocols using higher acceleration factors (MB 3 or MB 4 × GRAPPA 2) can likely provide even greater gains in sensitivity but should be carefully optimized to minimize the possibility of false activations.

•

MB factors 1, 2, 4, and 6 and two reconstructions were evaluated for fMRI performance.

•

MB accelerations 2, 4, and 6 improved BOLD sensitivity metrics over MB 1.

•

False-positive activation due to signal leakage was seen at high accelerations.

•

Use of Split Slice-GRAPPA reconstruction significantly reduces false positives.

Prospective motion correction of 3D echo-planar imaging data for functional MRI using optical tracking

We evaluated the performance of an optical camera based prospective motion correction (PMC) system in improving the quality of 3D echo-planar imaging functional MRI data. An optical camera and external marker were used to dynamically track the head movement of subjects during fMRI scanning. PMC was performed by using the motion information to dynamically update the sequence's RF excitation and gradient waveforms such that the field-of-view was realigned to match the subject's head movement. Task-free fMRI experiments on five healthy volunteers followed a 2 × 2 × 3 factorial design with the following factors: PMC on or off; 3.0 mm or 1.5 mm isotropic resolution; and no, slow, or fast head movements. Visual and motor fMRI experiments were additionally performed on one of the volunteers at 1.5 mm resolution comparing PMC on vs PMC off for no and slow head movements. Metrics were developed to quantify the amount of motion as it occurred relative to k-space data acquisition. The motion quantification metric collapsed the very rich camera tracking data into one scalar value for each image volume that was strongly predictive of motion-induced artifacts. The PMC system did not introduce extraneous artifacts for the no motion conditions and improved the time series temporal signal-to-noise by 30% to 40% for all combinations of low/high resolution and slow/fast head movement relative to the standard acquisition with no prospective correction. The numbers of activated voxels (p < 0.001, uncorrected) in both task-based experiments were comparable for the no motion cases and increased by 78% and 330%, respectively, for PMC on versus PMC off in the slow motion cases. The PMC system is a robust solution to decrease the motion sensitivity of multi-shot 3D EPI sequences and thereby overcome one of the main roadblocks to their widespread use in fMRI studies.

•

Data quality was assessed in terms of tSNR, RMSE, and task fMRI t-statistics.

•

Correction improved tSNR by 30% to 40%.

•

Percent reduction in RMSE as a function of motion level was characterized.

•

Number of active voxels for motor and visual tasks increased significantly.

Optimizing RetroICor and RetroKCor corrections for multi-shot 3D FMRI acquisitions

Physiological noise, if unaccounted for, can drastically reduce the statistical significance of detected activation in FMRI. In this paper, we systematically optimize physiological noise regressions for multi-shot 3D FMRI data. First, we investigate whether 3D FMRI data are best corrected in image space (RetroICor) or k-space (RetroKCor), in which each k-space segment can be assigned its unique physiological phase. In addition, the optimal regressor set is determined using the Bayesian Information Criterion (BIC) for a variety of 3D acquisitions corresponding to different image contrasts and k-space readouts. Our simulations and experiments indicate that: (a) k-space corrections are more robust when performed on real/imaginary than magnitude/phase data; (b) k-space corrections do not outperform image-space corrections, despite the ability to synchronize physiological phase to acquisition time more accurately; and (c) the optimal model varied considerably between the various acquisition techniques. These results suggest the use of a tailored set of volume-wide regressors, determined by BIC or other selection criteria, that achieves optimal balance between variance reduction and potential over-fitting.

FMRI using balanced steady-state free precession (SSFP)

Steady-state free precession (SSFP) is a highly-efficient MRI pulse sequence that has been a fairly recent arrival in the functional MRI realm. Several methods for using balanced SSFP to detect the BOLD signal have been proposed to date and will be discussed in this review. After a brief introduction to the general properties of SSFP, this review describes the quite different approaches of transition-band and pass-band SSFP in terms of functional contrast mechanism. It then discusses the potential advantages of these techniques, followed by their challenges and shortcomings. Finally, it gives an overview of some applications considered to date and the author's perspective on where these techniques are headed. In the spirit of this special issue, the author also includes some of the personal history underlying her own explorations in this area.

High-resolution functional MRI at 3 T: 3D/2D echo-planar imaging with optimized physiological noise correction

High-resolution functional MRI (fMRI) offers unique possibilities for studying human functional neuroanatomy. Although high-resolution fMRI has proven its potential at 7 T, most fMRI studies are still performed at rather low spatial resolution at 3 T. We optimized and compared single-shot two-dimensional echo-planar imaging (EPI) and multishot three-dimensional EPI high-resolution fMRI protocols. We extended image-based physiological noise correction from two-dimensional EPI to multishot three-dimensional EPI. The functional sensitivity of both acquisition schemes was assessed in a visual fMRI experiment. The physiological noise correction increased the sensitivity significantly, can be easily applied, and requires simple recordings of pulse and respiration only. The combination of three-dimensional EPI with physiological noise correction provides exceptional sensitivity for 1.5 mm high-resolution fMRI at 3 T, increasing the temporal signal-to-noise ratio by more than 25% compared to two-dimensional EPI. Magn Reson Med, 2013. © 2012 Wiley Periodicals, Inc.

Advances in High-Field BOLD fMRI

This review article examines the current state of BOLD fMRI at a high magnetic field strength of 7 Tesla. The following aspects are covered: a short description of the BOLD contrast, spatial and temporal resolution, BOLD sensitivity, localization and spatial specificity, technical challenges as well as an outlook on future developments are given. It is shown that the main technical challenges of performing BOLD fMRI at high magnetic field strengths-namely development of array coils, imaging sequences and parallel imaging reconstruction-have been solved successfully. The combination of these developments has lead to the availability of high-resolution BOLD fMRI protocols that are able to cover the whole brain with a repetition time (TR) shorter than 3 s. The structural information available from these high-resolution fMRI images itself is already very detailed, which helps to co-localize structure and function. Potential future applications include whole-brain connectivity analysis on a laminar resolution and single subject examinations.