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.

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.

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.

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.

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.

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.

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.

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.

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.

Symptom-specific differential motor network modulation by deep brain stimulation in Parkinson's disease.

OBJECTIVE: Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is an established neurosurgical treatment for the motor symptoms of Parkinson's disease (PD). While often highly effective, DBS does not always yield optimal therapeutic outcomes, and stimulation-induced adverse effects, including paresthesia, muscle contractions, and nausea/lightheadedness, commonly occur and can limit the efficacy of stimulation. Currently, objective metrics do not exist for monitoring neural changes associated with stimulation-induced therapeutic and adverse effects. METHODS: In the present study, the authors combined intraoperative functional MRI (fMRI) with STN DBS in 20 patients with PD to test the hypothesis that stimulation-induced blood oxygen level-dependent signals contained predictive information concerning the therapeutic and adverse effects of stimulation. RESULTS: As expected, DBS resulted in blood oxygen level-dependent activation in myriad motor regions, including the primary motor cortex, caudate, putamen, thalamus, midbrain, and cerebellum. Across the patients, DBS-induced improvements in contralateral Unified Parkinson's Disease Rating Scale tremor subscores correlated with activation of thalamic, brainstem, and cerebellar regions. In addition, improvements in rigidity and bradykinesia subscores correlated with activation of the primary motor cortex. Finally, activation of specific sensorimotor-related subregions correlated with the presence of DBS-induced adverse effects, including paresthesia and nausea (cerebellar cortex, sensorimotor cortex) and unwanted muscle contractions (caudate and putamen). CONCLUSIONS: These results suggest that DBS-induced activation patterns revealed by fMRI contain predictive information with respect to the therapeutic and adverse effects of DBS. The use of fMRI in combination with DBS therefore may hold translational potential to guide and improve clinical stimulator optimization in patients.

The benefit of high-performance gradients on echo planar imaging for BOLD-based resting-state functional MRI.

Improved gradient performance in an MRI system reduces distortion in echo planar imaging (EPI), which has been a key imaging method for functional studies. A lightweight, low-cryogen compact 3T MRI scanner (C3T) is capable of achieving 80 mT m-1 gradient amplitude with 700 T m-1 s-1 slew rate, in comparison with a conventional whole-body 3T MRI scanner (WB3T, 50 mT m-1 with 200 T m-1 s-1). We investigated benefits of the high-performance gradients in a high-spatial-resolution (1.5 mm isotropic) functional MRI study. Reduced echo spacing in the EPI pulse sequence inherently leads to less severe geometric distortion, which provided higher accuracy than with WB3T for registration between EPI and anatomical images. The cortical coverage of C3T datasets was improved by more accurate signal depiction (i.e. less dropout or pile-up). Resting-state functional analysis results showed that greater magnitude and extent in functional connectivity (FC) for the C3T than the WB3T when the selected seed region is susceptible to distortions, while the FC matrix for well-known brain networks showed little difference between the two scanners. This shows that the improved quality in EPI is particularly valuable for studying certain brain regions typically obscured by severe distortion.

Reducing PNS with minimal performance penalties via simple pulse sequence modifications on a high-performance compact 3T scanner.

One of the major concerns associated with high-performance gradients is peripheral nerve stimulation (PNS) of the subject during MRI exams. Since the installation, more than 680 volunteer subjects (patients and controls) have been scanned on a compact 3 T MRI system with high-performance gradients, capable of 80 mT m-1 gradient amplitude and 700 T m-1 s-1 slew rate simultaneously. Despite PNS concerns associated with the high-performance gradients, due to the smaller physical dimensions of the gradient coils, minimal or no PNS sensation was reported with most pulse sequences. The exception was PNS reported by only five of 252 subjects (about 2%) scanned with a specific 3D fast spin echo pulse sequence (3DFLAIR). Rather than derating the entire system performance across all pulse sequences and all gradient lobes, we addressed reported PNS effect with a simple and specific modification to the targeted lobes of the problematic pulse sequence. in addition, the PNS convolutional model was adapted to predict sequence-specific PNS threshold level and its reduction after derating. The effectiveness of the targeted pulse sequence modification was demonstrated by successfully re-scanning four of the subjects who previously reported PNS sensations without further reported PNS. The pulse sequence modification did not result in noticeable degradation of image quality or substantial increase in scan time. The results demonstrated that PNS was rarely reported on the compact 3 T, and when it was, utilizing a specific modification of the gradient waveform causing PNS was an effective strategy, rather than derating the performance of the entire gradient system.

Multivariate pattern classification on BOLD activation pattern induced by deep brain stimulation in motor, associative, and limbic brain networks.

Deep brain stimulation (DBS) has been shown to be an effective treatment for movement disorders and it is now being extended to the treatment of psychiatric disorders. Functional magnetic resonance imaging (fMRI) studies indicate that DBS stimulation targets dependent brain network effects, in networks that respond to stimulation. Characterizing these patterns is crucial for linking DBS-induced therapeutic and adverse effects. Conventional DBS-fMRI, however, lacks the sensitivity needed for decoding multidimensional information such as spatially diffuse patterns. We report here on the use of a multivariate pattern analysis (MVPA) to demonstrate that stimulation of three DBS targets (STN, subthalamic nucleus; GPi, globus pallidus internus; NAc, nucleus accumbens) evoked a sufficiently distinctive blood-oxygen-level-dependent (BOLD) activation in swine brain. The findings indicate that STN and GPi evoke a similar motor network pattern, while NAc shows a districted associative and limbic pattern. The findings show that MVPA could be effectively applied to overlapping or sparse BOLD patterns which are often found in DBS. Future applications are expected employ MVPA fMRI to identify the proper stimulation target dependent brain circuitry for a DBS outcome.

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.

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.

Resting-state functional connectivity modulates the BOLD activation induced by nucleus accumbens stimulation in the swine brain.

INTRODUCTION: While the clinical efficacy of deep brain stimulation (DBS) the treatment of motor-related symptoms is well established, the mechanism of action of the resulting cognitive and behavioral effects has been elusive. METHODS: By combining functional magnetic resonance imaging (fMRI) and DBS, we investigated the pattern of blood-oxygenation-level-dependent (BOLD) signal changes induced by stimulating the nucleus accumbens in a large animal model. RESULTS: We found that diffused BOLD activation across multiple functional networks, including the prefrontal, limbic, and thalamic regions during the stimulation, resulted in a significant change in inter-regional functional connectivity. More importantly, the magnitude of the modulation was closely related to the strength of the inter-regional resting-state functional connectivity. CONCLUSIONS: Nucleus accumbens stimulation affects the functional activity in networks that underlie cognition and behavior. Our study provides an insight into the nature of the functional connectivity, which mediates activation effect via brain networks.

Mapping BOLD Activation by Pharmacologically Evoked Tremor in Swine.

Harmaline-induced tremor is one of the most commonly utilized disease models for essential tremor (ET). However, the underlying neural networks involved in harmaline-induced tremor and the degree to which these are a representative model of the pathophysiologic mechanism of ET are incompletely understood. In this study, we evaluated the functional brain network effects induced by systemic injection of harmaline using pharmacological functional magnetic resonance imaging (ph-fMRI) in the swine model. With harmaline administration, we observed significant activation changes in cerebellum, thalamus, and inferior olivary nucleus (ION). In addition, inter-regional correlations in activity between cerebellum and deep cerebellar nuclei and between cerebellum and thalamus were significantly enhanced. These harmaline-induced effects gradually decreased with repeated administration of drug, replicating the previously demonstrated 'tolerance' effect. This study demonstrates that harmaline-induced tremor is associated with activity changes in brain regions previously implicated in humans with ET. Thus, harmaline-induction of tremor in the swine may be a useful model to explore the neurological effects of novel therapeutic agents and/or neuromodulation techniques for ET.

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.