Signatures of Electrical Stimulation Driven Network Interactions in the Human Limbic System.

Stimulation-evoked signals are starting to be used as biomarkers to indicate the state and health of brain networks. The human limbic network, often targeted for brain stimulation therapy, is involved in emotion and memory processing. Previous anatomic, neurophysiological, and functional studies suggest distinct subsystems within the limbic network (Rolls, 2015). Studies using intracranial electrical stimulation, however, have emphasized the similarities of the evoked waveforms across the limbic network. We test whether these subsystems have distinct stimulation-driven signatures. In eight patients (four male, four female) with drug-resistant epilepsy, we stimulated the limbic system with single-pulse electrical stimulation. Reliable corticocortical evoked potentials (CCEPs) were measured between hippocampus and the posterior cingulate cortex (PCC) and between the amygdala and the anterior cingulate cortex (ACC). However, the CCEP waveform in the PCC after hippocampal stimulation showed a unique and reliable morphology, which we term the "limbic Hippocampus-Anterior nucleus of the thalamus-Posterior cingulate, HAP-wave." This limbic HAP-wave was visually distinct and separately decoded from the CCEP waveform in ACC after amygdala stimulation. Diffusion MRI data show that the measured end points in the PCC overlap with the end points of the parolfactory cingulum bundle rather than the parahippocampal cingulum, suggesting that the limbic HAP-wave may travel through fornix, mammillary bodies, and the anterior nucleus of the thalamus (ANT). This was further confirmed by stimulating the ANT, which evoked the same limbic HAP-wave but with an earlier latency. Limbic subsystems have unique stimulation-evoked signatures that may be used in the future to help network pathology diagnosis.SIGNIFICANCE STATEMENT The limbic system is often compromised in diverse clinical conditions, such as epilepsy or Alzheimer's disease, and characterizing its typical circuit responses may provide diagnostic insight. Stimulation-evoked waveforms have been used in the motor system to diagnose circuit pathology. We translate this framework to limbic subsystems using human intracranial stereo EEG (sEEG) recordings that measure deeper brain areas. Our sEEG recordings describe a stimulation-evoked waveform characteristic to the memory and spatial subsystem of the limbic network that we term the "limbic HAP-wave." The limbic HAP-wave follows anatomic white matter pathways from hippocampus to thalamus to the posterior cingulum and shows promise as a distinct biomarker of signaling in the human brain memory and spatial limbic network.

Calibration of concomitant field offsets using phase contrast MRI for asymmetric gradient coils.

PURPOSE: Asymmetric gradient coils introduce zeroth- and first-order concomitant field terms, in addition to higher-order terms common to both asymmetric and symmetric gradients. Salient to compensation strategies is the accurate calibration of the concomitant field spatial offset parameters for asymmetric coils. A method that allows for one-time calibration of the offset parameters is described. THEORY AND METHODS: A modified phase contrast pulse sequence with single-sided bipolar flow encoding is proposed to calibrate the offsets for asymmetric, transverse gradient coils. By fitting the measured phase offsets to different gradient amplitudes, the spatial offsets were calculated by fitting the phase variation. This was used for calibrating real-time pre-emphasis compensation of the zeroth- and first-order concomitant fields. RESULTS: Image quality improvement with the proposed corrections was demonstrated in phantom and healthy volunteers with non-Cartesian and Cartesian trajectory acquisitions. Concomitant field compensation using the calibrated offsets resulted in a residual phase error <3% at the highest gradient amplitude and demonstrated substantial reduction of image blur and slice position/selection artifacts. CONCLUSIONS: The proposed implementation provides an accurate method for calibrating spatial offsets that can be used for real-time concomitant field compensation of zeroth and first-order terms, substantially reducing artifacts without retrospective correction or sequence specific waveform modifications.

Improved Resting-State Functional MRI Using Multi-Echo Echo-Planar Imaging on a Compact 3T MRI Scanner with High-Performance Gradients.

In blood-oxygen-level-dependent (BOLD)-based resting-state functional (RS-fMRI) studies, usage of multi-echo echo-planar-imaging (ME-EPI) is limited due to unacceptable late echo times when high spatial resolution is used. Equipped with high-performance gradients, the compact 3T MRI system (C3T) enables a three-echo whole-brain ME-EPI protocol with smaller than 2.5 mm isotropic voxel and shorter than 1 s repetition time, as required in landmark fMRI studies. The performance of the ME-EPI was comprehensively evaluated with signal variance reduction and region-of-interest-, seed- and independent-component-analysis-based functional connectivity analyses and compared with a counterpart of single-echo EPI with the shortest TR possible. Through the multi-echo combination, the thermal noise level is reduced. Functional connectivity, as well as signal intensity, are recovered in the medial orbital sulcus and anterior transverse collateral sulcus in ME-EPI. It is demonstrated that ME-EPI provides superior sensitivity and accuracy for detecting functional connectivity and/or brain networks in comparison with single-echo EPI. In conclusion, the high-performance gradient enabled high-spatial-temporal resolution ME-EPI would be the method of choice for RS-fMRI study on the C3T.

Rapid Geometry-Corrected Echo-Planar Diffusion Imaging at Ultrahigh Field: Fusing View Angle Tilting and Point-Spread Function Mapping.

PURPOSE: Severe geometric distortions induced by tissue susceptibility, water-fat chemical shift, and eddy currents pose a substantial obstacle in single-shot EPI, especially for high-resolution imaging at ultrahigh field. View angle tilting (VAT)-EPI can mitigate in-plane distortion. However, the accompanied strong image blurring prevented its widespread applications. On the other hand, point-spread function mapping (PSF)-EPI can correct distortion and blurring accurately but requires prolonged scan time. We present fused VAT-PSF-EPI and possibilities for acceleration. METHODS: MR signal equations were explicitly derived to quantify image blurring in VAT-EPI and the maximum acceleration capacity in VAT-PSF-EPI. To validate the theoretical prediction, phantom measurements with varying in-plane parallel imaging factors, slice thicknesses, and RF pulses were conducted at 7 Tesla. In addition, in vivo human brain scans were acquired with T2 and diffusion weighting to assess distortion and blurring correction. RESULTS: VAT can effectively suppress distortion, and the introduced image blurring is corrected through PSF encoding. Up to fourfold acceleration (only 5 shots) in VAT-PSF-EPI was achieved compared with standard PSF-EPI without VAT. VAT-induced signal loss was mitigated by adjusting the sequence parameters and EPI resolution. In vivo T2 -weighted EPI data with 1.4 mm3 resolution demonstrate immunity to water-fat chemical shift-induced distortion. Very high-spatial resolution diffusion-weighted EPI (0.7 × 0.7 × 2.8 mm3 and 1.2 mm3 ) demonstrates the immunity to eddy current-induced distortion. CONCLUSION: VAT-PSF-EPI is a novel spin-echo EPI-based sequence for fast high-resolution diffusion imaging at ultrahigh field.

Impact of prospective motion correction, distortion correction methods and large vein bias on the spatial accuracy of cortical laminar fMRI at 9.4 Tesla.

Functional imaging with sub-millimeter spatial resolution is a basic requirement for assessing functional MRI (fMRI) responses across different cortical depths and is used extensively in the emerging field of laminar fMRI. Such studies seek to investigate the detailed functional organization of the brain and may develop to a new powerful tool for human neuroscience. However, several studies have shown that measurement of laminar fMRI responses can be biased by the image acquisition and data processing strategies. In this work, measurements with three different gradient-echo EPI BOLD fMRI protocols with a voxel size down to 650 â€‹µm isotropic were performed at 9.4 â€‹T. We estimated how prospective motion correction can help to improve spatial accuracy by reducing the number of spatial resampling steps in postprocessing. In addition, we demonstrate key requirements for accurate geometric distortion correction to ensure that distortion correction maps are properly aligned to the functional data and that strong variations of distortions near large veins can lead to signal overlays which cannot be corrected for during postprocessing. Furthermore, this study illustrates the spatial extent of bias induced by pial and other larger veins in laminar BOLD experiments. Since these issues under investigation affect studies performed with more conventional spatial resolutions, the methods applied in this work may also help to improve the understanding of the BOLD signal more broadly.

The effect of spiral trajectory correction on pseudo-continuous arterial spin labeling with high-performance gradients on a compact 3T scanner.

PURPOSE: To demonstrate the feasibility of pseudo-continuous arterial-spin-labeled (pCASL) imaging with 3D fast-spin-echo stack-of-spirals on a compact 3T scanner (C3T), to perform trajectory correction for eddy-current-induced deviations in the spiral readout of pCASL imaging, and to assess the correction effect on perfusion-related images with high-performance gradients (80 mT/m, 700T/m/s) of the C3T. METHODS: To track eddy-current-induced artifacts with Archimedean spiral readout, the spiral readout in pCASL imaging was performed with 5 different peak gradient slew rate (Smax ) values ranging from 70 to 500 T/m/s. The trajectory for each Smax was measured using a dynamic field camera and applied in a density-compensated gridding image reconstruction in addition to the nominal trajectory. The effect of the trajectory correction was assessed with perfusion-weighted (ΔM) images and proton-density-weighted images as well as cerebral blood flow (CBF) maps, obtained from 10 healthy volunteers. RESULTS: Blurring artifact on ΔM images was mitigated by the trajectory correction. CBF values on the left and right calcarine cortices showed no significant difference after correction. Also, the signal-to-noise ratio of ΔM images improved, on average, by 7.6% after correction (P < .001). The greatest improvement of 12.1% on ΔM images was achieved with a spiral readout using Smax of 300~400 T/m/s. CONCLUSION: Eddy currents can cause spiral trajectory deviation, which leads to deformation of the CBF map even in cases of low value Smax . The trajectory correction for spiral-readout-based pCASL produces more reliable results for perfusion imaging. These results suggest that pCASL is feasible on C3T with high-performance gradients.

Distortion-free imaging: A double encoding method (DIADEM) combined with multiband imaging for rapid distortion-free high-resolution diffusion imaging on a compact 3T with high-performance gradients.

BACKGROUND: Distortion-free, high-resolution diffusion imaging using DIADEM (Distortion-free Imaging: A Double Encoding Method), proposed recently, has great potential for clinical applications. However, it can suffer from prolonged scan times and its reliability for quantitative diffusion imaging has not been evaluated. PURPOSE: To investigate the clinical feasibility of DIADEM-based high-resolution diffusion imaging on a novel compact 3T (C3T) by evaluating the reliability of quantitative diffusion measurements and utilizing both the high-performance gradients (80 mT/m, 700 T/m/s) and the sequence optimization with the navigator acquisition window reduction and simultaneous multislice (multiband) imaging. STUDY TYPE: Prospective feasibility study. PHANTOM/SUBJECTS: Diffusion quality control phantom scans to evaluate the reliability of quantitative diffusion measurements; 36 normal control scans for B0 -field mapping; six healthy and two patient subject scans with a brain tumor for comparisons of diffusion and anatomical imaging. FIELD STRENGTH/SEQUENCE: 3T; the standard single-shot echo-planar-imaging (EPI), multishot DIADEM diffusion, and anatomical (2D-FSE [fast-spin-echo], 2D-FLAIR [fluid-attenuated-inversion-recovery], and 3D-MPRAGE [magnetization prepared rapid acquisition gradient echo]) imaging. ASSESSMENT: The scan time reduction, the reliability of quantitative diffusion measurements, and the clinical efficacy for high-resolution diffusion imaging in healthy control and brain tumor volunteers. STATISTICAL TEST: Bland-Altman analysis. RESULTS: The scan time for high in-plane (0.86 mm2 ) resolution, distortion-free, and whole brain diffusion imaging were reduced from 10 to 5 minutes with the sequence optimizations. All of the mean apparent diffusion coefficient (ADC) values in phantom were within the 95% confidence interval in the Bland-Altman plot. The proposed acquisition with a total off-resonance coverage of 597.2 Hz wider than the expected bandwidth of 500 Hz in human brain could yield a distortion-free image without foldover artifacts. Compared with EPI, therefore, this approach allowed direct image matching with the anatomical images and enabled improved delineation of the tumor boundaries. DATA CONCLUSION: The proposed high-resolution diffusion imaging approach is clinically feasible on C3T due to a combination of hardware and sequence improvements. LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY: Stage 1 J. Magn. Reson. Imaging 2020;51:296-310.

Dynamic 2D self-phase-map Nyquist ghost correction for simultaneous multi-slice echo planar imaging.

PURPOSE: To develop a reconstruction pipeline that intrinsically accounts for both simultaneous multislice echo planar imaging (SMS-EPI) reconstruction and dynamic slice-specific Nyquist ghosting correction in time-series data. METHODS: After 1D slice-group average phase correction, the separate polarity (i.e., even and odd echoes) SMS-EPI data were unaliased by slice GeneRalized Autocalibrating Partial Parallel Acquisition. Both the slice-unaliased even and odd echoes were jointly reconstructed using a model-based framework, extended for SMS-EPI reconstruction that estimates a 2D self-phase map, corrects dynamic slice-specific phase errors, and combines data from all coils and echoes to obtain the final images. RESULTS: The percentage ghost-to-signal ratios (%GSRs) and its temporal variations for MB3Ry 2 with a field of view/4 shift in a human brain obtained by the proposed dynamic 2D and standard 1D phase corrections were 1.37 ± 0.11 and 2.66 ± 0.16, respectively. Even with a large regularization parameter λ applied in the proposed reconstruction, the smoothing effect in fMRI activation maps was comparable to a very small Gaussian kernel size 1 × 1 × 1 mm3 . CONCLUSION: The proposed reconstruction pipeline reduced slice-specific phase errors in SMS-EPI, resulting in reduction of GSR. It is applicable for functional MRI studies because the smoothing effect caused by the regularization parameter selection can be minimal in a blood-oxygen-level-dependent activation map.

Reduced acoustic noise in diffusion tensor imaging on a compact MRI system.

PURPOSE: To investigate the feasibility of substantially reducing acoustic noise while performing diffusion tensor imaging (DTI) on a compact 3T (C3T) MRI scanner equipped with a 42-cm inner-diameter asymmetric gradient. METHODS: A-weighted acoustic measurements were made using 10 mT/m-amplitude sinusoidal waveforms, corresponding to echo-planar imaging (EPI) echo spacing of 0.25 to 5.0 ms, on a conventional, whole-body 3T MRI and on the C3T. Acoustic measurements of DTI with trapezoidal EPI waveforms were then made at peak gradient performance on the C3T (80 mT/m amplitude, 700 T/m/s slew rate) and at derated performance (33 mT/m, 10 to 50 T/m/s) for acoustic noise reduction. DTI was acquired in two different phantoms and in seven human subjects, with and without gradient-derating corresponding to multi- and single-shot acquisitions, respectively. RESULTS: Sinusoidal waveforms on the C3T were quieter by 8.5 to 15.6 A-weighted decibels (dBA) on average as compared to the whole-body MRI. The derated multishot DTI acquisition noise level was only 8.7 dBA (at 13 T/m/s slew rate) above ambient, and was quieter than non-derated, single-shot DTI by 22.3 dBA; however, the scan time was almost quadrupled. Although derating resulted in negligible diffusivity differences in the phantoms, small biases in diffusivity measurements were observed in human subjects (apparent diffusion coefficient = +9.3 ± 8.8%, fractional anisotropy = +3.2 ± 11.2%, radial diffusivity = +9.4 ± 16.8%, parallel diffusivity = +10.3 ± 8.4%). CONCLUSION: The feasibility of achieving reduced acoustic noise levels with whole-brain DTI on the C3T MRI was demonstrated. Magn Reson Med 79:2902-2911, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

High-resolution distortion-free diffusion imaging using hybrid spin-warp and echo-planar PSF-encoding approach.

High-resolution diffusion-weighted imaging (DWI) has great potential to provide unique information about tissue microstructure in-vivo. Although single-shot echo-planar imaging (EPI) is a most popular tool for DWI, its application for high-resolution DWI is limited due to T2* blurring and susceptibility- and eddy-current-induced geometric distortions, especially at ultra-high field (UHF) such as 7T. In this study, we adapt a hybrid spin-warp and echo-planar encoding strategy inspired by point spread function (PSF) mapping and optimize it for high-resolution and distortion-free diffusion imaging applications. More specifically, a 2D navigator echo is added into the original sequence for shot-to-shot motion-induced phase error estimation and correction. The spatial encoding is shared between the PSF and the EPI phase encoding dimension allowing short echo trains to preserve the diffusion and navigator signals efficiently at UHF, where T2 decay is relatively fast. In addition, variable k-space spacing was applied in the PSF dimension and combined with parallel imaging in the EPI-PE dimension to further accelerate the PSF acquisition. The results demonstrate that this method can yield isotropic submillimeter resolution without T2* blurring and geometric distortions at 7T and enables a clear and detailed delineation of human brain structures in-vivo with the diffusion contrasts. In addition, results of the proposed approach for high-resolution diffusion imaging at 3T are presented.

Correction of metal-induced susceptibility artifacts for functional MRI during deep brain stimulation.

Functional magnetic resonance imaging (fMRI) is an emerging tool for investigating brain activation associated with, or modulated by, deep brain stimulation (DBS). However, DBS-fMRI generally suffers from severe susceptibility to artifacts in regions near the metallic stimulation electrodes, as well as near tissue/air boundaries of the brain. These result in strong intensity and geometric distortions along the phase-encoding (PE) (i.e., blipped) direction in gradient-echo echo-planar imaging (GE-EPI). Distortion presents a major challenge to conducting reliable data analysis and in interpreting the findings. A recent study showed that the point spread function (PSF) mapping-based reverse gradient approach has a potential to correct for distortions not only in spin-echo EPI, but also in GE-EPI acquired in both the forward and reverse PE directions. In this study, we adapted that approach in order to minimize severe metal-induced susceptibility artifacts for DBS-fMRI, and to evaluate the performance of the approach in a phantom study and a large animal DBS-fMRI study. The method combines the distortion-corrected GE-EPI pair with geometrically different intensity distortions due to the opposing encoding directions. The results demonstrate that the approach can minimize susceptibility artifacts that appear around the metallic electrodes, as well as in the regions near the tissue/air boundaries in the brain. We also demonstrated that an accurate geometric correction is important in improving BOLD contrast in the group dataset, especially in regions where strong susceptibility artifacts appear.

Model-based iterative reconstruction for single-shot EPI at 7T.

PURPOSE: To describe a model-based reconstruction strategy for single-shot echo planar imaging (EPI) that intrinsically accounts for k-space nonuniformity, Nyquist ghosting, and geometric distortions during rather than before or after image reconstruction. METHODS: Ramp sampling and inhomogeneous B0 field-induced distortion cause the EPI samples to lie on a non-Cartesian grid, thus requiring the nonuniform fast Fourier transform. Additionally, a 2D Nyquist ghost phase correction without the need for extra navigator acquisition is included in the proposed reconstruction. Coil compression is also incorporated to reduce the computational load. The proposed method is applied to phantom and human brain MRI data. RESULTS: The results demonstrate that Nyquist ghosting and geometric distortions are reduced by the proposed reconstruction. The proposed 2D phase correction is superior to a conventional 1D correction. The reductions of both artifacts lead to improved temporal signal-to-noise ratio (tSNR). The virtual coil results suggest that the processing time can be reduced by up to 75%, with a mean tSNR loss of only 3.2% when using 8-virtual instead of 32-physical coils for twofold undersampled data. CONCLUSION: The proposed reconstruction improves the quality (ghosting, geometry, and tSNR) of EPI without requiring calibration data for Nyquist ghost correction. Magn Reson Med 78:2250-2264, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

High-resolution diffusion MRI at 7T using a three-dimensional multi-slab acquisition.

High-resolution diffusion MRI can provide the ability to resolve small brain structures, enabling investigations of detailed white matter architecture. A major challenge for in vivo high-resolution diffusion MRI is the low signal-to-noise ratio. In this work, we combine two highly compatible methods, ultra-high field and three-dimensional multi-slab acquisition to improve the SNR of high-resolution diffusion MRI. As each kz plane is encoded using a single-shot echo planar readout, scan speeds of the proposed technique are similar to the commonly used two-dimensional diffusion MRI. In-plane parallel acceleration is applied to reduce image distortions. To reduce the sensitivity of auto-calibration signal data to subject motion and respiration, several new adaptions of the fast low angle excitation echo-planar technique (FLEET) that are suitable for 3D multi-slab echo planar imaging are proposed and evaluated. A modified reconstruction scheme is proposed for auto-calibration with the most robust method, Slice-FLEET acquisition, to make it compatible with navigator correction of motion induced phase errors. Slab boundary artefacts are corrected using the nonlinear slab profile encoding method recently proposed by our group. In vivo results demonstrate that using 7T and three-dimensional multi-slab acquisition with improved auto-calibration signal acquisition and nonlinear slab boundary artefacts correction, high-quality diffusion MRI data with ~1mm isotropic resolution can be achieved.

PSF mapping-based correction of eddy-current-induced distortions in diffusion-weighted echo-planar imaging.

PURPOSE: To accurately correct diffusion-encoding direction-dependent eddy-current-induced geometric distortions in diffusion-weighted echo-planar imaging (DW-EPI) and to minimize the calibration time at 7 Tesla (T). METHODS: A point spread function (PSF) mapping based eddy-current calibration method is newly presented to determine eddy-current-induced geometric distortions even including nonlinear eddy-current effects within the readout acquisition window. To evaluate the temporal stability of eddy-current maps, calibration was performed four times within 3 months. Furthermore, spatial variations of measured eddy-current maps versus their linear superposition were investigated to enable correction in DW-EPIs with arbitrary diffusion directions without direct calibration. For comparison, an image-based eddy-current correction method was additionally applied. Finally, this method was combined with a PSF-based susceptibility-induced distortion correction approach proposed previously to correct both susceptibility and eddy-current-induced distortions in DW-EPIs. RESULTS: Very fast eddy-current calibration in a three-dimensional volume is possible with the proposed method. The measured eddy-current maps are very stable over time and very similar maps can be obtained by linear superposition of principal-axes eddy-current maps. High resolution in vivo brain results demonstrate that the proposed method allows more efficient eddy-current correction than the image-based method. CONCLUSION: The combination of both PSF-based approaches allows distortion-free images, which permit reliable analysis in diffusion tensor imaging applications at 7T.

Thermometry mapping during CT-guided thermal ablations: proof of feasibility and internal validation using spectral CT.

Objective. Conventional computed tomography (CT) imaging does not provide quantitative information on local thermal changes during percutaneous ablative therapy of cancerous and benign tumors, aside from few qualitative, visual cues. In this study, we have investigated changes in CT signal across a wide range of temperatures and two physical phases for two different tissue mimicking materials, each.Approach. A series of experiments were conducted using an anthropomorphic phantom filled with water-based gel and olive oil, respectively. Multiple, clinically used ablation devices were applied to locally cool or heat the phantom material and were arranged in a configuration that produced thermal changes in regions with inconsequential amounts of metal artifact. Eight fiber optic thermal sensors were positioned in the region absent of metal artifact and were used to record local temperatures throughout the experiments. A spectral CT scanner was used to periodically acquire and generate electron density weighted images. Average electron density weighted values in 1 mm3volumes of interest near the temperature sensors were computed and these data were then used to calculate thermal volumetric expansion coefficients for each material and phase.Main results. The experimentally determined expansion coefficients well-matched existing published values and variations with temperature-maximally differing by 5% of the known value. As a proof of concept, a CT-generated temperature map was produced during a heating time point of the water-based gel phantom, demonstrating the capability to map changes in electron density weighted signal to temperature.Significance. This study has demonstrated that spectral CT can be used to estimate local temperature changes for different materials and phases across temperature ranges produced by thermal ablations.

Model-based iterative reconstruction for direct imaging with point spread function encoded echo planar MRI.

BACKGROUND: Echo planar imaging (EPI) is a fast measurement technique commonly used in magnetic resonance imaging (MRI), but is highly sensitive to measurement non-idealities in reconstruction. Point spread function (PSF)-encoded EPI is a multi-shot strategy which alleviates distortion, but acquisition of encodings suitable for direct distortion-free imaging prolongs scan time. In this work, a model-based iterative reconstruction (MBIR) framework is introduced for direct imaging with PSF-EPI to improve image quality and acceleration potential. METHODS: An MBIR platform was developed for accelerated PSF-EPI. The reconstruction utilizes a subspace representation, is regularized to promote local low-rankedness (LLR), and uses variable splitting for efficient iteration. Comparisons were made against standard reconstructions from prospectively accelerated PSF-EPI data and with retrospective subsampling. Exploring aggressive partial Fourier acceleration of the PSF-encoding dimension, additional comparisons were made against an extension of Homodyne to direct PSF-EPI in numerical experiments. A neuroradiologists' assessment was completed comparing images reconstructed with MBIR from retrospectively truncated data directly against images obtained with standard reconstructions from non-truncated datasets. RESULTS: Image quality results were consistently superior for MBIR relative to standard and Homodyne reconstructions. As the MBIR signal model and reconstruction allow for arbitrary sampling of the PSF space, random sampling of the PSF-encoding dimension was also demonstrated, with quantitative assessments indicating best performance achieved through nonuniform PSF sampling combined with partial Fourier. With retrospective subsampling, MBIR reconstructs high-quality images from sub-minute scan datasets. MBIR was shown to be superior in a neuroradiologists' assessment with respect to three of five performance criteria, with equivalence for the remaining two. CONCLUSIONS: A novel image reconstruction framework is introduced for direct imaging with PSF-EPI, enabling arbitrary PSF space sampling and reconstruction of diagnostic-quality images from highly accelerated PSF-encoded EPI data.

Minimizing susceptibility-induced BOLD sensitivity loss in multi-band accelerated fMRI using point spread function mapping and gradient reversal.

Objective. Interleaved reverse-gradient fMRI (RG-fMRI) with a point-spread-function (PSF) mapping-based distortion correction scheme has the potential to minimize signal loss in echo-planar-imaging (EPI). In this work, the RG-fMRI is further improved by imaging protocol optimization and application of reverse Fourier acquisition.Approach. Multi-band imaging was adapted for RG-fMRI to improve the temporal and spatial resolution. To better understand signal dropouts in forward and reverse EPIs, a simple theoretical relationship between echo shift and geometric distortion was derived and validated by the reliable measurements using PSF mapping method. After examining practical imaging protocols for RG-fMRI in three subjects on both a conventional whole-body and a high-performance compact 3 T, the results were compared and the feasibility to further improve the RG-fMRI scheme were explored. High-resolution breath-holding RG-fMRI was conducted with nine subjects on the compact 3 T and the fMRI reliability improvement in high susceptibility brain regions was demonstrated. Finally, reverse Fourier acquisition was applied to RG-fMRI, and its benefit was assessed by a simulation study based on the breath-holding RG-fMRI data.Main results. The temporal and spatial resolution of the multi-band RG-fMRI became feasible for whole-brain fMRI. Echo shift measurements from PSF mapping well estimated signal dropout effects in the EPI pair and were useful to further improve the RG-fMRI scheme. Breath-holding RG-fMRI demonstrated improved fMRI reliability in high susceptibility brain regions. Reverse partial Fourier acquisition omitting the late echoes could further improve the temporal or spatial resolution for RG-fMRI without noticeable signal degradation and spatial resolution loss.Significance. With the improved imaging scheme, RG-fMRI could reliably investigate the functional mechanisms of the human brain in the temporal and frontal areas suffering from susceptibility-induced functional sensitivity loss.

Locally Low-Rank Denoising of Multi-Echo Functional MRI Data With Application in Resting-State Analysis.

OBJECTIVES: Locally low-rank (LLR) denoising of functional magnetic resonance imaging (fMRI) time series image data is extended to multi-echo (ME) data. The proposed method extends the capabilities of non-physiologic noise suppression beyond single-echo applications with a dedicated ME algorithm. MATERIALS AND METHODS: Following an institutional review board (IRB) approved protocol, resting-state fMRI data were acquired in 7 healthy subjects. A compact 3T scanner enabled whole-brain acquisition of multiband ME fMRI data at high spatial resolution (1.4 × 1.4 × 2.8 mm 3 ) with a 1810 ms repetition time (TR). Image data were denoised with ME-LLR preceding functional processing. The results of connectivity maps generated from denoised data were compared with maps generated with equivalent processing of non-denoised images. To assess ME-LLR as a method to reduce scan time, comparisons were made between maps computed from image data with full and retrospectively truncated durations. Assessments were completed with seed-based connectivity analyses using echo-combined image data. In a feasibility assessment, nondenoised and denoised full-duration echo-combined data were equivalently processed with independent component analysis (ICA) and compared. RESULTS: ME-LLR denoising yielded strengthened resting-state network connectivity maps after nuisance regression and seed-based connectivity analysis. In assessing ME-LLR as a scan reduction mechanism, maps generated from denoised data at half scan time showed comparable quality with maps generated from full-duration, non-denoised data, at both single subject and group levels. ME-LLR substantially increased temporal signal-to-noise ratio (tSNR) for image data respective to each individual echo and for image data after nuisance regression. Among echo-specific image volumes, increases in tSNR yielded by ME-LLR were most pronounced for image data with the longest echo time and thereby lowest SNR. ICA showed resting-state networks consistently identified between non-denoised and denoised data, with clearer demarcation of networks for ME-LLR. CONCLUSIONS: ME-LLR is demonstrated to suppress non-physiologic noise, enhance functional connectivity map quality, and could potentially facilitate scan time reduction in ME-fMRI.

Enhanced clinical task-based fMRI metrics through locally low-rank denoising of complex-valued data.

OBJECTIVE: This study investigates a locally low-rank (LLR) denoising algorithm applied to source images from a clinical task-based functional MRI (fMRI) exam before post-processing for improving statistical confidence of task-based activation maps. METHODS: Task-based motor and language fMRI was obtained in eleven healthy volunteers under an IRB approved protocol. LLR denoising was then applied to raw complex-valued image data before fMRI processing. Activation maps generated from conventional non-denoised (control) data were compared with maps derived from LLR-denoised image data. Four board-certified neuroradiologists completed consensus assessment of activation maps; region-specific and aggregate motor and language consensus thresholds were then compared with nonparametric statistical tests. Additional evaluation included retrospective truncation of exam data without and with LLR denoising; a ROI-based analysis tracked t-statistics and temporal SNR (tSNR) as scan durations decreased. A test-retest assessment was performed; retest data were matched with initial test data and compared for one subject. RESULTS: fMRI activation maps generated from LLR-denoised data predominantly exhibited statistically significant (p = 4.88×10-4 to p = 0.042; one p = 0.062) increases in consensus t-statistic thresholds for motor and language activation maps. Following data truncation, LLR data showed task-specific increases in t-statistics and tSNR respectively exceeding 20 and 50% compared to control. LLR denoising enabled truncation of exam durations while preserving cluster volumes at fixed thresholds. Test-retest showed variable activation with LLR data thresholded higher in matching initial test data. CONCLUSION: LLR denoising affords robust increases in t-statistics on fMRI activation maps compared to routine processing, and offers potential for reduced scan duration while preserving map quality.

Evaluation of a compact 3 T MRI scanner for patients with implanted devices.

Access to high-quality MR exams is severely limited for patients with some implanted devices due to labeled MR safety conditions, but small-bore systems can overcome this limitation. For example, a compact 3 T MR scanner (C3T) with high-performance gradients can acquire exams of the head, extremities, and infants. Because of its reduced bore size and the patient being advanced only partially into the bore, the associated electromagnetic (EM) fields drop off rapidly caudal to the head, compared to whole-body systems. Therefore, some patients with MR conditional implanted devices can safely receive 3 T brain exams on the C3T using its strong gradients and a multiple-channel receive coil, while a corresponding exam on whole-body MR is precluded. The purpose of this study is to evaluate the performance of a small-bore scanner for subjects with MR conditional spinal or sacral nerve stimulators, or abandoned cardiac implantable electronic device (CIED) leads. The spatial dependence of specific absorption rate (SAR) on the C3T was compared to whole-body scanners. A device assessment tool was developed and applied to evaluate MR safety individually on the C3T for 12 subjects with implanted devices or abandoned CIED leads. Once MR safety was established, the subjects received a C3T brain exam along with their clinical, 1.5 T exam. The resulting images were graded by three board-certified neuroradiologists. The C3T exams were well-tolerated with no adverse events, and significantly outperformed the whole-body 1.5 T exams in terms of overall image quality.