Application of k-space energy spectrum analysis to susceptibility field mapping and distortion correction in gradient-echo EPI

Unsupervised deep learning model for correcting Nyquist ghosts of single-shot spatiotemporal encoding

PURPOSE

To design an unsupervised deep learning (DL) model for correcting Nyquist ghosts of single-shot spatiotemporal encoding (SPEN) and evaluate the model for real MRI applications.

METHODS

The proposed method consists of three main components: (1) an unsupervised network that combines Residual Encoder and Restricted Subspace Mapping (RERSM-net) and is trained to generate a phase-difference map based on the even and odd SPEN images; (2) a spin physical forward model to obtain the corrected image with the learned phase difference map; and (3) cycle-consistency loss that is explored for training the RERSM-net.

RESULTS

The proposed RERSM-net could effectively generate smooth phase difference maps and correct Nyquist ghosts of single-shot SPEN. Both simulation and real in vivo MRI experiments demonstrated that our method outperforms the state-of-the-art SPEN Nyquist ghost correction method. Furthermore, the ablation experiments of generating phase-difference maps show the advantages of the proposed unsupervised model.

CONCLUSION

The proposed method can effectively correct Nyquist ghosts for the single-shot SPEN sequence.

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.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO 2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Comparison and evaluation of distortion correction techniques on an MR-guided radiotherapy system

PURPOSE

To evaluate two distortion correction techniques for diffusion-weighted single-shot echo-planar imaging (DW-ssEPI) on a 0.35 T magnetic resonance-guided radiotherapy (MRgRT) system.

METHODS

The effects of sequence optimization through enabling parallel imaging (PI) and selecting appropriate bandwidth on spatial distortion were first evaluated on the 0.35 T MRgRT system using a spatial integrity phantom. Field map (FM) and reversed gradient (RG) corrections were then performed on the optimized protocol to further reduce distortion. An open-source toolbox was used to quantify the spatial displacement before and after distortion correction. To evaluate ADC accuracy and repeatability of the optimized protocol, as well as impacts of distortion correction on ADC values, the optimized protocol was scanned twice on a diffusion phantom. The calculated ADC values were compared with reference ADCs using paired t-test. Intraclass correlation coefficient (ICC) between the two repetitions, as well as between before and after FM/RG correction was calculated to evaluate ADC repeatability and effects of distortion correction. Six patients were recruited to assess the in-vivo performance. The optimal distortion correction technique was identified by visual assessment. To quantify distortion reduction, tumor and critical structures were contoured on the turbo spin echo (TSE) image (reference image), the DW-ssEPI image, and the distortion corrected images independently by two radiation oncologists. Mean distance to agreement (MDA) and DICE coefficient between contours on the reference images and the diffusion images were calculated. Tumor apparent diffusion coefficient (ADC) values from the original DW-ssEPI images and the distortion corrected images were compared using Bland-Altman analysis.

RESULTS

Sequence optimization played a vital role in improving the spatial integrity, and spatial distortion was proportional to the total readout time. Before the correction, distortion of the optimized protocol (PI and high bandwidth) was 1.50 ± 0.89 mm in a 100 mm radius and 2.21 ± 1.39 mm in a 175 mm radius for the central plane. FM corrections reduced the distortions to 0.42 ± 0.27 mm and 0.67 ± 0.49 mm respectively, and RG reduced distortion to 0.40 ± 0.22 mm and 0.64 ± 0.47 mm, respectively. The optimized protocol provided accurate and repeatable ADC quantification on the diffusion phantom. The calculated ADC values were consistent before and after FM/RG correction. For the patient study, the FM correction was unable to reduce chemical shift artifacts whereas the RG method successfully mitigated the chemical shift. MDA reduced from 2.52 ± 1.29 mm to 1.11 ± 0.72 mm after the RG correction. The DICE coefficient increased from 0.80 ± 0.13 to 0.91 ± 0.06. A Bland-Altman plot showed that there was a good agreement between ADC measurements before and after application of the RG correction.

CONCLUSION

Two distortion correction techniques were evaluated on a commercial low-field MRgRT system. Overall, the RG correction was able to drastically improve spatial distortion and preserve ADC accuracy.

Deconvolution-based distortion correction of EPI using analytic single-voxel point-spread functions

PURPOSE

To develop a postprocessing algorithm that corrects geometric distortions due to spatial variations of the static magnetic field amplitude, B₀ , and effects from relaxation during signal acquisition in EPI.

THEORY AND METHODS

An analytic, complex point-spread function is deduced for k-space trajectories of EPI variants and applied to corresponding acquisitions in a resolution phantom and in human volunteers at 3 T. With the analytic point-spread function and experimental maps of B₀ (and, optionally, the effective transverse relaxation time, T 2 * ) as input, a point-spread function matrix operator is devised for distortion correction by a Thikonov-regularized deconvolution in image space. The point-spread function operator provides additional information for an appropriate correction of the signal intensity distribution. A previous image combination algorithm for acquisitions with opposite phase blip polarities is adapted to the proposed method to recover destructively interfering signal contributions.

RESULTS

Applications of the proposed deconvolution-based distortion correction ("DecoDisCo") algorithm demonstrate excellent distortion corrections and superior performance regarding the recovery of an undistorted intensity distribution in comparison to a multifrequency reconstruction. Examples include full and partial Fourier standard EPI scans as well as double-shot center-out trajectories. Compared with other distortion-correction approaches, DecoDisCo permits additional deblurring to obtain sharper images in cases of significant T 2 * effects.

CONCLUSION

Robust distortion corrections in EPI acquisitions are feasible with high quality by regularized deconvolution with an analytic point-spread function. The general algorithm, which is publicly released on GitHub, can be straightforwardly adapted for specific EPI variants or other acquisition schemes.

Improved PRF-based MR thermometry using k-space energy spectrum analysis

PURPOSE

Proton resonance frequency (PRF) thermometry encodes information in the phase of MRI signals. A multiplicative factor converts phase changes into temperature changes, and this factor includes the TE. However, phase variations caused by B₀ and/or B₁ inhomogeneities can effectively change TE in ways that vary from pixel to pixel. This work presents how spatial phase variations affect temperature maps and how to correct for corresponding errors.

METHODS

A method called "k-space energy spectrum analysis" was used to map regions in the object domain to regions in the k-space domain. Focused ultrasound heating experiments were performed in tissue-mimicking gel phantoms under two scenarios: with and without proper shimming. The second scenario, with deliberately de-adjusted shimming, was meant to emulate B₀ inhomogeneities in a controlled manner. The TE errors were mapped and compensated for using k-space energy spectrum analysis, and corrected results were compared with reference results. Furthermore, a volunteer was recruited to help evaluate the magnitude of the errors being corrected.

RESULTS

The in vivo abdominal results showed that the TE and heating errors being corrected can readily exceed 10%. In phantom results, a linear regression between reference and corrected temperatures results provided a slope of 0.971 and R² of 0.9964. Analysis based on the Bland-Altman method provided a bias of -0.0977°C and 95% limits of agreement that were 0.75°C apart.

CONCLUSION

Spatially varying TE errors, such as caused by B₀ and/or B₁ inhomogeneities, can be detected and corrected using the k-space energy spectrum analysis method, for increased accuracy in proton resonance frequency thermometry.

The use of Fourier-domain analyses for unwrapping phase images of low SNR

PURPOSE

We report a new postprocessing procedure that uses Fourier-domain data analyses to improve the accuracy and reliability of phase unwrapping for MRI data of low SNR.

METHODS

The developed method first identifies the Fourier-domain energy peak locations corresponding to different image-domain areas from which a robust measurement of image-domain phase gradients can be obtained even for MRI data of low SNR. The phase-gradient information measured from critical brain regions using the above-mentioned Fourier-domain analysis is then combined with the conventional temporal-domain or spatial-domain phase-unwrapping procedure to remove phase wraps. The developed method was tested with MRI data obtained from 30 healthy adult volunteers, and its performance was quantitatively evaluated.

RESULTS

The developed Fourier-domain analysis could robustly quantify image-domain phase gradients even for MRI data with low SNR (e.g., SNR ≃ 2). Experimental results show that the Fourier-domain analyses could further reduce phase wrap artifact in data produced by the conventional temporal-domain or spatial-domain phase-unwrapping procedures.

CONCLUSION

Our results demonstrate that the developed phase-unwrapping method can reduce residual phase wraps resulting from conventional procedures in critical brain regions (e.g., near the air-tissue interfaces) and should prove valuable for studies that require accurate measurements of MRI phase values, such as QSM, B₀ field mapping, and temperature mapping.

The Impact of Echo Time Shifts and Temporal Signal Fluctuations on BOLD Sensitivity in Presurgical Planning at 7 T

OBJECTIVES

Gradients in the static magnetic field caused by tissues with differing magnetic susceptibilities lead to regional variations in the effective echo time, which modifies both image signal and BOLD sensitivity. Local echo time changes are not considered in the most commonly used metric for BOLD sensitivity, temporal signal-to-noise ratio (tSNR), but may be significant, particularly at ultrahigh field close to air cavities (such as the sinuses and ear canals) and near gross brain pathologies and postoperative sites.

MATERIALS AND METHODS

We have studied the effect of local variations in echo time and tSNR on BOLD sensitivity in 3 healthy volunteers and 11 patients with tumors, postoperative cavities, and venous malformations at 7 T. Temporal signal-to-noise ratio was estimated from a 5-minute run of resting state echo planar imaging with a nominal echo time of 22 milliseconds. Maps of local echo time were derived from the phase of a multiecho GE scan. One healthy volunteer performed 10 runs of a breath-hold task. The t-map from this experiment served as a criterion standard BOLD sensitivity measure. Two runs of a less demanding breath-hold paradigm were used for patients.

RESULTS

In all subjects, a strong reduction in the echo time (from 22 milliseconds to around 11 milliseconds) was found close to the ear canals and sinuses. These regions were characterized by high tSNR but low t-values in breath-hold t-maps. In some patients, regions of particular interest in presurgical planning were affected by reductions in the echo time to approximately 13-15 milliseconds. These included the primary motor cortex, Broca's area, and auditory cortex. These regions were characterized by high tSNR values (70 and above). Breath-hold results were corrupted by strong motion artifacts in all patients.

CONCLUSIONS

Criterion standard BOLD sensitivity estimation using hypercapnic experiments is challenging, especially in patient populations. Taking into consideration the tSNR, commonly used for BOLD sensitivity estimation, but ignoring local reductions in the echo time (eg, from 22 to 11 milliseconds), would erroneously suggest functional sensitivity sufficient to map BOLD signal changes. It is therefore important to consider both local variations in the echo time and temporal variations in signal, using the product metric of these two indices for instance. This should ensure a reliable estimation of BOLD sensitivity and to facilitate the identification of potential false-negative results. This is particularly true at high fields, such as 7 T and in patients with large pathologies and postoperative cavities.

Correction for Eddy Current-Induced Echo-Shifting Effect in Partial-Fourier Diffusion Tensor Imaging

In most diffusion tensor imaging (DTI) studies, images are acquired with either a partial-Fourier or a parallel partial-Fourier echo-planar imaging (EPI) sequence, in order to shorten the echo time and increase the signal-to-noise ratio (SNR). However, eddy currents induced by the diffusion-sensitizing gradients can often lead to a shift of the echo in k-space, resulting in three distinct types of artifacts in partial-Fourier DTI. Here, we present an improved DTI acquisition and reconstruction scheme, capable of generating high-quality and high-SNR DTI data without eddy current-induced artifacts. This new scheme consists of three components, respectively, addressing the three distinct types of artifacts. First, a k-space energy-anchored DTI sequence is designed to recover eddy current-induced signal loss (i.e., Type 1 artifact). Second, a multischeme partial-Fourier reconstruction is used to eliminate artificial signal elevation (i.e., Type 2 artifact) associated with the conventional partial-Fourier reconstruction. Third, a signal intensity correction is applied to remove artificial signal modulations due to eddy current-induced erroneous T₂^∗-weighting (i.e., Type 3 artifact). These systematic improvements will greatly increase the consistency and accuracy of DTI measurements, expanding the utility of DTI in translational applications where quantitative robustness is much needed.

Reference-free unwarping of single-shot spatiotemporally encoded MRI using asymmetric self-refocused echoes acquisition

This paper presents a phase evolution rewinding algorithm for correcting the geometric and intensity distortions in single-shot spatiotemporally encoded (SPEN) MRI with acquisition of asymmetric self-refocused echo trains. Using the field map calculated from the phase distribution of the source image, the off-resonance induced phase errors are successfully rewound through deconvolution. The alias-free partial Fourier transform reconstruction helps improve the signal-to-noise ratio of the field maps and the output images. The effectiveness of the proposed algorithm was validated through 7 T MRI experiments on a lemon, a water phantom, and in vivo rat head. SPEN imaging was evaluated using rapid acquisition by sequential excitation and refocusing (RASER) which produces uniform T2 weighting. The results indicate that the new technique can more robustly deal with the cases in which the images obtained with conventional single-shot spin-echo EPI are difficult to be restored due to serious field variations.

Reduced field-of-view imaging for single-shot MRI with an amplitude-modulated chirp pulse excitation and Fourier transform reconstruction

PURPOSE

We employ an amplitude-modulated chirp pulse to selectively excite spins in one or more regions of interest (ROIs) to realize reduced field-of-view (rFOV) imaging based on single-shot spatiotemporally encoded (SPEN) sequence and Fourier transform reconstruction.

MATERIALS AND METHODS

The proposed rFOV imaging method was theoretically analyzed and illustrated with numerical simulation and tested with phantom experiments and in vivo rat experiments. In addition, point spread function was applied to demonstrate the feasibility of the proposed method. To evaluate the proposed method, the rFOV results were compared with those obtained using the EPI method with orthogonal RF excitation.

RESULTS

The simulation and experimental results show that the proposed method can image one or two separated ROIs along the SPEN dimension in a single shot with higher spatial resolution, less sensitive to field inhomogeneity, and practically no aliasing artifacts. In addition, the proposed method may produce rFOV images with comparable signal-to-noise ratio to the rFOV EPI images.

CONCLUSION

The proposed method is promising for the applications under severe susceptibility heterogeneities and for imaging separate ROIs simultaneously.

Accurate field mapping in the presence of B0 inhomogeneities, applied to MR thermometry

PURPOSE

To describe how B0 inhomogeneities can cause errors in proton resonance frequency (PRF) shift thermometry, and to correct for these errors.

METHODS

With PRF thermometry, measured phase shifts are converted into temperature measurements through the use of a scaling factor proportional to the echo time, TE. However, B0 inhomogeneities can deform, spread, and translate MR echoes, potentially making the "true" echo time vary spatially within the imaged object and take on values that differ from the prescribed TE value. Acquisition and reconstruction methods able to avoid or correct for such errors are presented.

RESULTS

Tests were performed in a gel phantom during sonication, and temperature measurements were made with proper shimming as well as with intentionally introduced B0 inhomogeneities. Errors caused by B0 inhomogeneities were observed, described, and corrected by the proposed methods. No statistical difference was found between the corrected results and the reference results obtained with proper shimming, while errors by more than 10% in temperature elevation were corrected for. The approach was also applied to an abdominal in vivo dataset.

CONCLUSION

Field variations induce errors in measured field values, which can be detected and corrected. The approach was validated for a PRF thermometry application.

Alterations in brain transition metals in Huntington disease: an evolving and intricate story

BACKGROUND

Aberrant accumulation of transition metals in the brain may have an early and important role in the pathogenesis of several neurodegenerative disorders, including Huntington disease (HD).

OBJECTIVE

To comprehensively evaluate and validate the distribution of metal deposition in the brain using advanced magnetic resonance imaging methods from the premanifest through symptomatic stages of HD.

DESIGN

Observational study.

SETTING

University imaging center.

PARTICIPANTS

Twenty-eight HD expanded gene carriers, 34 patients with symptomatic HD, and 56 age- and sex-matched healthy control subjects were included in the study.

INTERVENTIONS

Participants underwent magnetic resonance imaging for the quantification of the phase evolution of susceptibility-weighted images.

MAIN OUTCOME MEASURES

To verify the identity of the metals responsible for the changes in the phase evolution of the susceptibility signal in the brain and to assess correlations with systemic levels. Inductively coupled plasma mass spectrometry was used to measure transition metal concentrations in postmortem brains.

RESULTS

In the basal ganglia, progressive increases in the phase evolution were found in HD, beginning in premanifest individuals who were far from expected onset and increasing with proximity to expected onset and thereafter. Increases in the cerebral cortex were regionally selective and present only in symptomatic HD. Increases were verified by excessive deposition of brain iron, but a complex alteration in other transition metals was found.

CONCLUSION

An important and early role of altered metal homeostasis is suggested in the pathogenesis of HD.

A correction method for streak artifacts in gradient-echo EPI using spin-echo EPI reference data

OBJECTIVE

To analyze the streak artifacts in a gradient-echo echo planar imaging (GE-EPI) sequence and to propose a correction method for the Nyquist ghost artifacts that does not cause streak artifacts in the GE-EPI imaging.

MATERIALS AND METHODS

Several GE-EPI imaging experiments with various reference scans, using both GE-EPI and SE-EPI scan data, were performed to analyze the streak artifacts and to investigate the spin dephasing phenomena of the GE-EPI reference scan. In addition, the analysis based on the spin dephasing was undertaken in order to demonstrate that the SE-EPI reference data can be used for the correction of the GE-EPI main scan data.

RESULTS

The experimental results confirmed that the improvement of the reference data using either signal averaging or a large flip angle cannot guarantee perfect correction of the streak artifact if the noise is not completely removed. Due to the main field inhomogeneity, the spins of the GE-EPI reference data were dephased in multiple echo signals. The proposed correction method, which uses a SE-EPI reference scan for the GE-EPI images, eliminates the N/2 ghost artifacts without producing streak artifacts.

CONCLUSION

It is believed that the proposed phase error correction scheme can improve the EPI performance in high field MRIs with higher magnetic field inhomogeneities.

Positive contrast technique for the detection and quantification of superparamagnetic iron oxide nanoparticles in MRI

In vivo detection and quantification of cells labeled with superparamagnetic iron oxide (SPIO) nanoparticles has been attracting increasing attention. In particular, positive contrast methods, such as susceptibility gradient mapping (SGM) and phase gradient mapping (PGM), have been proposed for the improved detection of SPIO nanoparticles. In this study, a different implementation of the PGM method is introduced; it calculates the phase gradient in the image space using a fast Fourier transform without the need for phase unwrapping. We first compared positive contrast generation between the PGM and SGM methods, which estimates the susceptibility gradient in k space through echo shift measurements. Next, PGM was applied to quantify SPIO concentrations by fitting the resulting phase gradient maps to those of a theoretical model. MR experiments were conducted using a 3-T magnet scanner to acquire two datasets: the first was acquired from a gelatin phantom with three SPIO-doped vials of different concentrations, and the second was obtained in vivo from a nude rat with SPIO-labeled C6 glioma cells implanted in the flanks. The sensitivity of the PGM and SGM methods was compared using various factors, including different SPIO concentrations, TEs and signal-to-noise ratios. Based on the theoretical model of an infinite cylinder, the results demonstrated that, without loss of spatial resolution, the PGM method presents positive contrast maps with a higher sensitivity than SGM at medium and low SPIO concentrations, whereas SGM is more sensitive than PGM at longer TEs. The quantification of SPIO concentrations using the phantom dataset was also reported. On the basis of the same infinite cylinder model, it was shown that the PGM method provides an accurate estimation of SPIO concentration.

Mapping cortical representations of the rodent forepaw and hindpaw with BOLD fMRI reveals two spatial boundaries

Electrical stimulation of the rat forepaw and hindpaw was employed to study the spatial distribution of BOLD fMRI. Averaging of multiple fMRI sessions significantly improved the spatial stability of the BOLD signal and enabled quantitative determination of the boundaries of the BOLD fMRI maps. The averaged BOLD fMRI signal was distributed unevenly over the extent of the map and the data at the boundaries could be modeled with major and minor spatial components. Comparison of three-dimensional echo-planar imaging (EPI) fMRI at isotropic 300 μm resolution demonstrated that the border locations of the major spatial component of BOLD signal did not overlap between the forepaw and hindpaw maps. Interestingly, the border positions of the minor BOLD fMRI spatial components extended significantly into neighboring representations. Similar results were found for cerebral blood volume (CBV) weighted fMRI obtained using iron oxide particles, suggesting that the minor spatial components may not be due to vascular mislocalization typically associated with BOLD fMRI. Comparison of the BOLD fMRI maps of the forepaw and hindpaw to histological determination of these representations using cytochrome oxidase (CO) staining demonstrated that the major spatial component of the BOLD fMRI activation maps accurately localizes the borders. Finally, 2-3 weeks following peripheral nerve denervation, cortical reorganization/plasticity at the boundaries of somatosensory limb representations in adult rat brain was studied. Denervation of the hindpaw caused a growth in the major component of forepaw representation into the adjacent border of hindpaw representation, such that fitting to two components no longer led to a better fit as compared to using one major component. The border of the representation after plasticity was the same as the border of its minor component in the absence of any plasticity. It is possible that the minor components represent either vascular effects that extend from the real neuronal representations or the neuronal communication between neighboring regions. Either way the results will be useful for studying mechanisms of plasticity that cause alterations in the boundaries of neuronal representations.

Noninvasive assessment of atherosclerotic plaque progression in ApoE-/- mice using susceptibility gradient mapping

BACKGROUND

Macrophages have been identified as a major contributor to plaque development and destabilization in atherosclerosis. The aim of this study was to noninvasively assess uptake of citrate coated very small iron oxide particles at different stages of plaque development in the brachiocephalic artery of apoE(-/-) mice. Susceptibility gradient mapping (SGM) was applied to generate positive contrast images and to quantify iron oxide uptake.

METHODS AND RESULTS

ApoE(-/-) mice were fed a high-fat diet for 4, 8, or 12 weeks; 300 μmol Fe/kg was injected 24 and 48 hours before final MRI. Increasing very small iron oxide particle uptake was observed over the course of atherosclerotic plaque development. Simultaneous administration of pravastatin led to a significant decrease in very small iron oxide particle uptake, assessed by mass spectroscopy and histology. SGM-MRI allowed the generation of positive contrast images, and magnitudes (mT/m) of contrast enhancement in SG parameter maps significantly correlated with the absolute iron oxide content (R(2)=0.70, P<0.05) and the macrophage density (R(2)=0.71, P<0.05).

CONCLUSIONS

This study shows an increase in iron oxide uptake (measured by in vivo SGM-MRI, histology, and mass spectroscopy) with the progression of plaque development in an apoE(-/-) mouse model of accelerated atherosclerosis. Positive contrast provided by SGM-MRI allowed for a clear visualization of intraplaque iron oxide depositions, and magnitudes (mT/m) of contrast enhancement in SG parameter maps allowed for the quantification of intraplaque iron oxide particles.

Application of k-space energy spectrum analysis for inherent and dynamic B0 mapping and deblurring in spiral imaging

Spiral imaging is vulnerable to spatial and temporal variations of the amplitude of the static magnetic field (B(0)) caused by susceptibility effects, eddy currents, chemical shifts, subject motion, physiological noise, and system instabilities, resulting in image blurring. Here, a novel off-resonance correction method is proposed to address these issues. A k-space energy spectrum analysis algorithm is first applied to inherently and dynamically generate a B(0) map from the k-space data at each time point, without requiring any additional data acquisition, pulse sequence modification, or phase unwrapping. A simulated phase evolution rewinding algorithm and an automatic residual deblurring algorithm are then used to correct for the blurring caused by both spatial and temporal B(0) variations, resulting in a high spatial and temporal fidelity. This method is validated against conventional B(0) mapping and deblurring methods, and its advantages for dynamic MRI applications are demonstrated in functional MRI studies.

A joint acquisition-reconstruction paradigm for correcting inhomogeneity artifacts in MR echo planar imaging

One of the main sources of signal degradation in rapid MR acquisitions, such as Echo Planar Imaging (EPI), is magnetic field variations caused by field inhomogeneities and susceptibility gradients. If unaccounted for during the reconstruction process, this spatially-varying field can cause severe image artifacts. In this paper, we show that correcting for the resulting degradations can be formulated as a blind image deconvolution problem. We propose a novel joint acquisition-processing paradigm to solve this problem. We describe a practical implementation of this paradigm using a multi-image acquisition strategy and a corresponding joint estimation-reconstruction algorithm. The estimation step computes the spatial distribution of the field maps, while the reconstruction step yields a Minimum Mean Squared Error (MMSE) estimate of the imaged slice. Our simulations show that this proposed joint acquisition-reconstruction method is robust and efficient, offering factors of improvement in the quality of the reconstructed image as compared to other traditional methods.

Direct magnetic field estimation based on echo planar raw data

Gradient recalled echo echo planar imaging is widely used in functional magnetic resonance imaging. The fast data acquisition is, however, very sensitive to field inhomogeneities which manifest themselves as artifacts in the images. Typically used correction methods have the common deficit that the data for the correction are acquired only once at the beginning of the experiment, assuming the field inhomogeneity distribution B(0) does not change over the course of the experiment. In this paper, methods to extract the magnetic field distribution from the acquired k-space data or from the reconstructed phase image of a gradient echo planar sequence are compared and extended. A common derivation for the presented approaches provides a solid theoretical basis, enables a fair comparison and demonstrates the equivalence of the k-space and the image phase based approaches. The image phase analysis is extended here to calculate the local gradient in the readout direction and improvements are introduced to the echo shift analysis, referred to here as "k-space filtering analysis." The described methods are compared to experimentally acquired B(0) maps in phantoms and in vivo. The k-space filtering analysis presented in this work demonstrated to be the most sensitive method to detect field inhomogeneities.

Stimulated echo induced misestimates on diffusion tensor indices and its remedy

PURPOSE

To report possible erroneous estimates of diffusion parameters in the twice-refocused spin-echo (TRSE) technique, proposed to eliminate eddy-current-induced geometric distortions in diffusion-weighted echo-planar imaging, when stimulated echo signals are inappropriately included.

MATERIALS AND METHODS

Eleven subjects were included for imaging experiments on two 1.5 Tesla systems using the TRSE sequence. Three versions, two with unbalanced crusher gradients inserted to dephase the stimulated echo from the b = 0 images and one with balanced crusher gradients, were implemented. The apparent diffusion coefficients (ADC) and fractional anisotropy (FA) were derived and compared.

RESULTS

The ADCs obtained with unbalanced crusher gradients were closer to values reported in the literature. Stimulated echo led to ADC over-estimations by 34.2%, 50.4%, 54.0%, 51.5%, 24.0%, and 41.9% in the genu of corpus callosum, splenium of corpus callosum, bilateral corona radiata, internal capsule, mediofrontal gyrus, and the cuneus, respectively (P < 0.01), with concomitant reduction in FA in highly anisotropic regions. Over-estimations of diffusion coefficients were found to be roughly equal along all directions.

CONCLUSION

Formation of stimulated echo in the TRSE technique can lead to erroneous estimations of the diffusion parameters, even if no prominent morphological artifacts are seen.

Recovery of signal loss due to an in-plane susceptibility gradient in the gradient echo EPI through acquisition of extended phase-encoding lines

In gradient echo imaging the in-plane susceptibility gradient causes an echo shift which results in signal loss. The loss of signal becomes more severe in gradient echo EPI, due to the low amplitude of the gradient which is applied in the phase-encoding direction during a long echo train. As the readout gradient amplitude is set to be very high in gradient echo EPI, the echo shift in the readout direction is negligible compared to that in the phase-encoding direction. Traditionally, a z-shimming technique has been applied to the phase-encoding direction of gradient echo EPI to restore the lost signal. This technique, however, requires a significant increase of scan time, as is also the case with the through-plane z-shimming technique. A new approach that allows one to restore the lost signal is to acquire additional phase-encoding lines beyond the regular phase-encoding range. The extension of the phase-encoding lines prior to the regular phase-encoding range exploits the delay time for optimum echo time of the BOLD sensitivity. Therefore, scan time is increased only for the extended phase-encoding lines posterior to the regular phase-encoding range. This technique has been confirmed experimentally by imaging human subject's heads at 3T.

Susceptibility gradient quantization by MRI signal response mapping (SIRMA) to dephaser

PURPOSE

Susceptibility effects are a very efficient source of contrast in magnetic resonance imaging. However, detection is hampered by the fact the induced contrast is negative. In this work, the SIgnal Response MApping (SIRMA) to dephaser method is proposed to map susceptibility gradient to improve visualization.

METHODS

In conventional gradient echo acquisitions, the echo formation of susceptibility affected spins is shifted in k-space, the shift being proportional to the susceptibility gradient. Susceptibility gradients map can be produced by measuring this induced shifts. The SIRMA method measures these shifts from a series of dephased images collected with additional incremental dephasers. These additional dephasers correspond either to a slice refocusing gradient offset or to a reconstruction window off-centering. The signal intensity profile as a function of the additional dephaser was determined on a pixel-by-pixel basis from the ensemble of dephased images. Susceptibility affected voxels presented a signal response profile maximum shifted compared to nonaffected voxels ones. Shift magnitude and sign were measured for each pixel to determine susceptibility gradients and produce a susceptibility gradient map.

RESULTS

In vitro experiments demonstrated the ability of the method to map gradient inhomogeneities induced by a cylinder. Quantization accuracy was evaluated comparing SIRMA images and simulations performed on the well-characterized air filled cylinder model. Performances of the SIRMA method, evaluated in vitro on cylinders filled with various superparamagnetic iron oxide SPIO concentrations, showed limited influence of acquisition parameters. Robustness of the method was then assessed in vivo after an infusion of SPIO-loaded nanocapsules into the rat brain using a convection-enhanced drug delivery approach. The region of massive susceptibility gradient induced by the SPIO-loaded nanocapsules was clearly delineated on SIRMA maps and images were compared to T2* weighted images, Susceptibility Gradient Map (SGM), and histological Perl's staining slice. The potential for quantitative evaluation of SPIO distribution volume was demonstrated.

CONCLUSIONS

The proposed method is a promising technique for a wide range of applications especially in molecular or cellular imaging with respect to its quantitative nature and its computational simplicity.

Fieldmap-free retrospective registration and distortion correction for EPI-based functional imaging

We describe a method for correcting the distortions present in echo planar images (EPI) and registering the EPI to structural MRI. A fieldmap is predicted from an air / tissue segmentation of the MRI using a perturbation method and subsequently used to unwarp the EPI data. Shim and other missing parameters are estimated by registration. We obtain results that are similar to those obtained using fieldmaps, however neither fieldmaps, nor knowledge of shim coefficients is required.

Event-related single-shot volumetric functional magnetic resonance inverse imaging of visual processing

Developments in multi-channel radio-frequency (RF) coil array technology have enabled functional magnetic resonance imaging (fMRI) with higher degrees of spatial and temporal resolution. While modest improvement in temporal acceleration has been achieved by increasing the number of RF coils, the maximum attainable acceleration in parallel MRI acquisition is intrinsically limited only by the amount of independent spatial information in the combined array channels. Since the geometric configuration of a large-n MRI head coil array is similar to that used in EEG electrode or MEG SQUID sensor arrays, the source localization algorithms used in MEG or EEG source imaging can be extended to also process MRI coil array data, resulting in greatly improved temporal resolution by minimizing k-space traversal during signal acquisition. Using a novel approach, we acquire multi-channel MRI head coil array data and then apply inverse reconstruction methods to obtain volumetric fMRI estimates of blood oxygenation level dependent (BOLD) contrast at unprecedented whole-brain acquisition rates of 100 ms. We call this combination of techniques magnetic resonance Inverse Imaging (InI), a method that provides estimates of dynamic spatially-resolved signal change that can be used to construct statistical maps of task-related brain activity. We demonstrate the sensitivity and inter-subject reliability of volumetric InI using an event-related design to probe the hemodynamic signal modulations in primary visual cortex. Robust results from both single subject and group analyses demonstrate the sensitivity and feasibility of using volumetric InI in high temporal resolution investigations of human brain function.

Linear constraint minimum variance beamformer functional magnetic resonance inverse imaging

Accurate estimation of the timing of neural activity is required to fully model the information flow among functionally specialized regions whose joint activity underlies perception, cognition and action. Attempts to detect the fine temporal structure of task-related activity would benefit from functional imaging methods allowing higher sampling rates. Spatial filtering techniques have been used in magnetoencephalography source imaging applications. In this work, we use the linear constraint minimal variance (LCMV) beamformer localization method to reconstruct single-shot volumetric functional magnetic resonance imaging (fMRI) data using signals acquired simultaneously from all channels of a high density radio-frequency (RF) coil array. The LCMV beamformer method generalizes the existing volumetric magnetic resonance inverse imaging (InI) technique, achieving higher detection sensitivity while maintaining whole-brain spatial coverage and 100 ms temporal resolution. In this paper, we begin by introducing the LCMV reconstruction formulation and then quantitatively assess its performance using both simulated and empirical data. To demonstrate the sensitivity and inter-subject reliability of volumetric LCMV InI, we employ an event-related design to probe the spatial and temporal properties of task-related hemodynamic signal modulations in primary visual cortex. Compared to minimum-norm estimate (MNE) reconstructions, LCMV offers better localization accuracy and superior detection sensitivity. Robust results from both single subject and group analyses demonstrate the excellent sensitivity and specificity of volumetric InI in detecting the spatial and temporal structure of task-related brain activity.

Improved image reconstruction for partial Fourier gradient-echo echo-planar imaging (EPI)

The partial Fourier gradient-echo echo planar imaging (EPI) technique makes it possible to acquire high-resolution functional MRI (fMRI) data at an optimal echo time. This technique is especially important for fMRI studies at high magnetic fields, where the optimal echo time is short and may not be achieved with a full Fourier acquisition scheme. In addition, it has been shown that partial Fourier EPI provides better anatomic resolvability than full Fourier EPI. However, the partial Fourier gradient-echo EPI may be degraded by artifacts that are not usually seen in other types of imaging. Those unique artifacts in partial Fourier gradient-echo EPI, to our knowledge, have not yet been systematically evaluated. Here we use the k-space energy spectrum analysis method to understand and characterize two types of partial Fourier EPI artifacts. Our studies show that Type 1 artifact, originating from k-space energy loss, cannot be corrected with pure postprocessing, and Type 2 artifact can be eliminated with an improved reconstruction method. We propose a novel algorithm, that combines images obtained from two or more reconstruction schemes guided by k-space energy spectrum analysis, to generate partial Fourier EPI with greatly reduced Type 2 artifact. Quality control procedures for avoiding Type 1 artifact in partial Fourier EPI are also discussed.

Correction for geometric distortion and N/2 ghosting in EPI by phase labeling for additional coordinate encoding (PLACE)

Echo-planar imaging (EPI) is vulnerable to geometric distortion and N/2 ghosting. These artifacts can be analyzed with an intuitive k-t space tool, and here we propose a simple method for their correction. In a slightly modified additional EPI acquisition, we sample the k-t space with a shift in k(y) by adding a small area to the phase-encoding (PE) gradient. Physically, the added gradient area creates a relative phase ramp across the object and directly encodes the undistorted original y-coordinate of each voxel into a phase difference between two distorted complex images, in a method called "phase labeling for additional coordinate encoding" (PLACE). The phase information is then used to map the mismapped signals back to their original locations for geometric and intensity correction. Smoothing of expanded complex data matrix effectively reduces noise in the differential phase map and allows subpixel warping. The two acquired images can also be averaged to effectively suppress the N/2 ghost. Efficient correction for both artifacts can be achieved with three acquisitions. These acquisitions can also serve as reference scans to correct for geometric distortion and/or N/2 ghost artifacts on all images in a time series. The technique was successfully demonstrated in phantom and animal studies.