Journal Papers

S

Taimouri, Vahid, Onur Afacan, Jeannette Perez-Rossello, Michael Callahan, Robert Mulkern, Simon Warfield, and Moti Freiman. (2015) 2015. “Spatially constrained incoherent motion method improves diffusion-weighted MRI signal decay analysis in the liver and spleen”. Med Phys 42 (4): 1895-903. https://doi.org/10.1118/1.4915495.
PURPOSE: To evaluate the effect of the spatially constrained incoherent motion (SCIM) method on improving the precision and robustness of fast and slow diffusion parameter estimates from diffusion-weighted MRI in liver and spleen in comparison to the independent voxel-wise intravoxel incoherent motion (IVIM) model. METHODS: We collected diffusion-weighted MRI (DW-MRI) data of 29 subjects (5 healthy subjects and 24 patients with Crohn's disease in the ileum). We evaluated parameters estimates' robustness against different combinations of b-values (i.e., 4 b-values and 7 b-values) by comparing the variance of the estimates obtained with the SCIM and the independent voxel-wise IVIM model. We also evaluated the improvement in the precision of parameter estimates by comparing the coefficient of variation (CV) of the SCIM parameter estimates to that of the IVIM. RESULTS: The SCIM method was more robust compared to IVIM (up to 70% in liver and spleen) for different combinations of b-values. Also, the CV values of the parameter estimations using the SCIM method were significantly lower compared to repeated acquisition and signal averaging estimated using IVIM, especially for the fast diffusion parameter in liver (CVIV IM = 46.61 ± 11.22, CVSCIM = 16.85 ± 2.160, p < 0.001) and spleen (CVIV IM = 95.15 ± 19.82, CVSCIM = 52.55 ± 1.91, p < 0.001). CONCLUSIONS: The SCIM method characterizes fast and slow diffusion more precisely compared to the independent voxel-wise IVIM model fitting in the liver and spleen.
Sui, Yao, Onur Afacan, Ali Gholipour, and Simon Warfield. 2020. “SLIMM: Slice localization integrated MRI monitoring”. Neuroimage 223: 117280. https://doi.org/10.1016/j.neuroimage.2020.117280.
Functional MRI (fMRI) is extremely challenging to perform in subjects who move because subject motion disrupts blood oxygenation level dependent (BOLD) signal measurement. It has become common to use retrospective framewise motion detection and censoring in fMRI studies to eliminate artifacts arising from motion. Data censoring results in significant loss of data and statistical power unless the data acquisition is extended to acquire more data not corrupted by motion. Acquiring more data than is necessary leads to longer than necessary scan duration, which is more expensive and may lead to additional subject non-compliance. Therefore, it is well established that real-time prospective motion monitoring is crucial to ensure data quality and reduce imaging costs. In addition, real-time monitoring of motion allows for feedback to the operator and the subject during the acquisition, to enable intervention to reduce the subject motion. The most widely used form of motion monitoring for fMRI is based on volume-to-volume registration (VVR), which quantifies motion as the misalignment between subsequent volumes. However, motion is not constrained to occur only at the boundaries of volume acquisition, but instead may occur at any time. Consequently, each slice of an fMRI acquisition may be displaced by motion, and assessment of whole volume to volume motion may be insensitive to both intra-volume and inter-volume motion that is revealed by displacement of the slices. We developed the first slice-by-slice self-navigated motion monitoring system for fMRI by developing a real-time slice-to-volume registration (SVR) algorithm. Our real-time SVR algorithm, which is the core of the system, uses a local image patch-based matching criterion along with a Levenberg-Marquardt optimizer, all accelerated via symmetric multi-processing, with interleaved and simultaneous multi-slice acquisition schemes. Extensive experimental results on real motion data demonstrated that our fast motion monitoring system, named Slice Localization Integrated MRI Monitoring (SLIMM), provides more accurate motion measurements than a VVR based approach. Therefore, SLIMM offers improved online motion monitoring which is particularly important in fMRI for challenging patient populations. Real-time motion monitoring is crucial for online data quality control and assurance, for enabling feedback to the subject and the operator to act to mitigate motion, and in adaptive acquisition strategies that aim to ensure enough data of sufficient quality is acquired without acquiring excess data.
Afacan, Onur, Scott Hoge, Tess Wallace, Ali Gholipour, Sila Kurugol, and Simon Warfield. 2020. “Simultaneous Motion and Distortion Correction Using Dual-Echo Diffusion-Weighted MRI”. J Neuroimaging 30 (3): 276-85. https://doi.org/10.1111/jon.12708.
BACKGROUND AND PURPOSE: Geometric distortions resulting from large pose changes reduce the accuracy of motion measurements and interfere with the ability to generate artifact-free information. Our goal is to develop an algorithm and pulse sequence to enable motion-compensated, geometric distortion compensated diffusion-weighted MRI, and to evaluate its efficacy in correcting for the field inhomogeneity and position changes, induced by large and frequent head motions. METHODS: Dual echo planar imaging (EPI) with a blip-reversed phase encoding distortion correction technique was evaluated in five volunteers in two separate experiments and compared with static field map distortion correction. In the first experiment, dual-echo EPI images were acquired in two head positions designed to induce a large field inhomogeneity change. A field map and a distortion-free structural image were acquired at each position to assess the ability of dual-echo EPI to generate reliable field maps and enable geometric distortion correction in both positions. In the second experiment, volunteers were asked to move to multiple random positions during a diffusion scan. Images were reconstructed using the dual-echo correction and a slice-to-volume registration (SVR) registration algorithm. The accuracy of SVR motion estimates was compared to externally measured ground truth motion parameters. RESULTS: Our results show that dual-echo EPI can produce slice-level field maps with comparable quality to field maps generated by the reference gold standard method. We also show that slice-level distortion correction improves the accuracy of SVR algorithms as slices acquired at different orientations have different levels of distortion, which can create errors in the registration process. CONCLUSIONS: Dual-echo acquisitions with blip-reversed phase encoding can be used to generate slice-level distortion-free images, which is critical for motion-robust slice to volume registration. The distortion corrected images not only result in better motion estimates, but they also enable a more accurate final diffusion image reconstruction.
Vasylechko, Serge Didenko, Simon Warfield, Onur Afacan, and Sila Kurugol. 2022. “Self-supervised IVIM DWI parameter estimation with a physics based forward model”. Magn Reson Med 87 (2): 904-14. https://doi.org/10.1002/mrm.28989.
PURPOSE: To assess the robustness and repeatability of intravoxel incoherent motion model (IVIM) parameter estimation for the diffusion-weighted MRI in the abdominal organs under the constraints of noisy diffusion signal using a novel neural network method. METHODS: Clinically acquired abdominal scans of Crohn's disease patients were retrospectively analyzed with regions segmented in the kidney cortex, spleen, liver, and bowel. A novel IVIM parameter fitting method based on the principle of a physics guided self-supervised convolutional neural network that does not require reference parameter estimates for training was compared to a conventional non-linear least squares (NNLS) algorithm, and a voxelwise trained artificial neural network (ANN). RESULTS: Results showed substantial increase in parameter robustness to the noise corrupted signal. In an intra-session repeatability experiment, the proposed method showed reduced coefficient of variation (CoV) over multiple acquisitions in comparison to conventional NLLS method and comparable performance to ANN. The use of D and f estimates from the proposed method led to the smallest misclassification error in linear discriminant analysis for characterization between normal and abnormal Crohn's disease bowel tissue. The fitting of D ∗ parameter remains to be challenging. CONCLUSION: The proposed method yields robust estimates of D and f IVIM parameters under the constraints of noisy diffusion signal. This indicates a potential for the use of the proposed method in conjunction with accelerated DW-MRI acquisition strategies, which would typically result in lower signal to noise ratio.
Sui, Yao, Onur Afacan, Camilo Jaimes, Ali Gholipour, and Simon Warfield. 2022. “Scan-Specific Generative Neural Network for MRI Super-Resolution Reconstruction”. IEEE Trans Med Imaging 41 (6): 1383-99. https://doi.org/10.1109/TMI.2022.3142610.
The interpretation and analysis of Magnetic resonance imaging (MRI) benefit from high spatial resolution. Unfortunately, direct acquisition of high spatial resolution MRI is time-consuming and costly, which increases the potential for motion artifact, and suffers from reduced signal-to-noise ratio (SNR). Super-resolution reconstruction (SRR) is one of the most widely used methods in MRI since it allows for the trade-off between high spatial resolution, high SNR, and reduced scan times. Deep learning has emerged for improved SRR as compared to conventional methods. However, current deep learning-based SRR methods require large-scale training datasets of high-resolution images, which are practically difficult to obtain at a suitable SNR. We sought to develop a methodology that allows for dataset-free deep learning-based SRR, through which to construct images with higher spatial resolution and of higher SNR than can be practically obtained by direct Fourier encoding. We developed a dataset-free learning method that leverages a generative neural network trained for each specific scan or set of scans, which in turn, allows for SRR tailored to the individual patient. With the SRR from three short duration scans, we achieved high quality brain MRI at an isotropic spatial resolution of 0.125 cubic mm with six minutes of imaging time for T2 contrast and an average increase of 7.2 dB (34.2%) in SNR to these short duration scans. Motion compensation was achieved by aligning the three short duration scans together. We assessed our technique on simulated MRI data and clinical data acquired from 15 subjects. Extensive experimental results demonstrate that our approach achieved superior results to state-of-the-art methods, while in parallel, performed at reduced cost as scans delivered with direct high-resolution acquisition.

R

Coll-Font, Jaume, Onur Afacan, Scott Hoge, Harsha Garg, Kumar Shashi, Bahram Marami, Ali Gholipour, Jeanne Chow, Simon Warfield, and Sila Kurugol. 2021. “Retrospective Distortion and Motion Correction for Free-Breathing DW-MRI of the Kidneys Using Dual-Echo EPI and Slice-to-Volume Registration”. J Magn Reson Imaging 53 (5): 1432-43. https://doi.org/10.1002/jmri.27473.
BACKGROUND: Diffusion-weighted MRI (DW-MRI) of the kidneys is a technique that provides information about the microstructure of renal tissue without requiring exogenous contrasts such as gadolinium, and it can be used for diagnosis in cases of renal disease and assessing response-to-therapy. However, physiological motion and large geometric distortions due to main B0 field inhomogeneities degrade the image quality, reduce the accuracy of quantitative imaging markers, and impede their subsequent clinical applicability. PURPOSE: To retrospectively correct for geometric distortion for free-breathing DW-MRI of the kidneys at 3T, in the presence of a nonstatic distortion field due to breathing and bulk motion. STUDY TYPE: Prospective. SUBJECTS: Ten healthy volunteers (ages 29-38, four females). FIELD STRENGTH/SEQUENCE: 3T; DW-MR dual-echo echo-planar imaging (EPI) sequence (10 b-values and 17 directions) and a T2 volume. ASSESSMENT: The distortion correction was evaluated subjectively (Likert scale 0-5) and numerically with cross-correlation between the DW images at b = 0 s/mm2 and a T2 volume. The intravoxel incoherent motion (IVIM) and diffusion tensor (DTI) model-fitting performance was evaluated using the root-mean-squared error (nRMSE) and the coefficient of variation (CV%) of their parameters. STATISTICAL TESTS: Statistical comparisons were done using Wilcoxon tests. RESULTS: The proposed method improved the Likert scores by 1.1 ± 0.8 (P < 0.05), the cross-correlation with the T2 reference image by 0.13 ± 0.05 (P < 0.05), and reduced the nRMSE by 0.13 ± 0.03 (P < 0.05) and 0.23 ± 0.06 (P < 0.05) for IVIM and DTI, respectively. The CV% of the IVIM parameters (slow and fast diffusion, and diffusion fraction for IVIM and mean diffusivity, and fractional anisotropy for DTI) was reduced by 2.26 ± 3.98% (P = 6.971 × 10-2 ), 11.24 ± 26.26% (P = 6.971 × 10-2 ), 4.12 ± 12.91% (P = 0.101), 3.22 ± 0.55% (P < 0.05), and 2.42 ± 1.15% (P < 0.05). DATA CONCLUSION: The results indicate that the proposed Di + MoCo method can effectively correct for time-varying geometric distortions and for misalignments due to breathing motion. Consequently, the image quality and precision of the DW-MRI model parameters improved. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY STAGE: 1.
Afacan, Onur, Tess Wallace, and Simon Warfield. 2020. “Retrospective correction of head motion using measurements from an electromagnetic tracker”. Magn Reson Med 83 (2): 427-37. https://doi.org/10.1002/mrm.27934.
PURPOSE: To investigate the feasibility of using an electromagnetic (EM) tracker to estimate rigid body head motion parameters, and using these measurements to retrospectively reduce motion artifacts. THEORY AND METHODS: A clinically used MPRAGE sequence was modified to measure motion using the EM tracking system once per repetition time. A retrospective k-space based motion correction algorithm that corrects for phase ramps (translation in image domain) and rotation of 3D k-space (rotation in image domain) was developed, using the parameters recorded using an EM tracker. The accuracy of the EM tracker for the purpose of motion measurement and correction was tested in phantoms, volunteers, and pediatric patients. RESULTS: Position localization was accurate to the order of 200 microns compared with registration localization in a phantom study. The quality of reconstructed images was assessed by computing the root mean square error, the structural similarity metric and average edge strength. Image quality improved consistently when motion correction was applied in both volunteer scans with deliberate head motion and in pediatric patient scans. In patients, the average edge strength improved significantly with retrospective motion correction, compared with images with no correction applied. CONCLUSIONS: EM tracking was effective in measuring head motion in the MRI scanner with high accuracy, and enabled retrospective reconstruction to improve image quality by reducing motion artifacts.
Rollins, Caitlin, Cynthia Ortinau, Christian Stopp, Kevin Friedman, Wayne Tworetzky, Borjan Gagoski, Clemente Velasco-Annis, et al. 2021. “Regional Brain Growth Trajectories in Fetuses with Congenital Heart Disease”. Ann Neurol 89 (1): 143-57. https://doi.org/10.1002/ana.25940.
OBJECTIVE: Congenital heart disease (CHD) is associated with abnormal brain development in utero. We applied innovative fetal magnetic resonance imaging (MRI) techniques to determine whether reduced fetal cerebral substrate delivery impacts the brain globally, or in a region-specific pattern. Our novel design included two control groups, one with and the other without a family history of CHD, to explore the contribution of shared genes and/or fetal environment to brain development. METHODS: From 2014 to 2018, we enrolled 179 pregnant women into 4 groups: "HLHS/TGA" fetuses with hypoplastic left heart syndrome (HLHS) or transposition of the great arteries (TGA), diagnoses with lowest fetal cerebral substrate delivery; "CHD-other," with other CHD diagnoses; "CHD-related," healthy with a CHD family history; and "optimal control," healthy without a family history. Two MRIs were obtained between 18 and 40 weeks gestation. Random effect regression models assessed group differences in brain volumes and relationships to hemodynamic variables. RESULTS: HLHS/TGA (n = 24), CHD-other (50), and CHD-related (34) groups each had generally smaller brain volumes than the optimal controls (71). Compared with CHD-related, the HLHS/TGA group had smaller subplate (-13.3% [standard error = 4.3%], p < 0.01) and intermediate (-13.7% [4.3%], p < 0.01) zones, with a similar trend in ventricular zone (-7.1% [1.9%], p = 0.07). These volumetric reductions were associated with lower cerebral substrate delivery. INTERPRETATION: Fetuses with CHD, especially those with lowest cerebral substrate delivery, show a region-specific pattern of small brain volumes and impaired brain growth before 32 weeks gestation. The brains of fetuses with CHD were more similar to those of CHD-related than optimal controls, suggesting genetic or environmental factors also contribute. ANN NEUROL 2021;89:143-157.
Balachandrasekaran, Arvind, Alexander Cohen, Onur Afacan, Simon Warfield, and Ali Gholipour. 2022. “Reducing the Effects of Motion Artifacts in fMRI: A Structured Matrix Completion Approach”. IEEE Trans Med Imaging 41 (1): 172-85. https://doi.org/10.1109/TMI.2021.3107829.
Functional MRI (fMRI) is widely used to study the functional organization of normal and pathological brains. However, the fMRI signal may be contaminated by subject motion artifacts that are only partially mitigated by motion correction strategies. These artifacts lead to distance-dependent biases in the inferred signal correlations. To mitigate these spurious effects, motion-corrupted volumes are censored from fMRI time series. Censoring can result in discontinuities in the fMRI signal, which may lead to substantial alterations in functional connectivity analysis. We propose a new approach to recover the missing entries from censoring based on structured low rank matrix completion. We formulated the artifact-reduction problem as the recovery of a super-resolved matrix from unprocessed fMRI measurements. We enforced a low rank prior on a large structured matrix, formed from the samples of the time series, to recover the missing entries. The recovered time series, in addition to being motion compensated, are also slice-time corrected at a fine temporal resolution. To achieve a fast and memory-efficient solution for our proposed optimization problem, we employed a variable splitting strategy. We validated the algorithm with simulations, data acquired under different motion conditions, and datasets from the ABCD study. Functional connectivity analysis showed that the proposed reconstruction resulted in connectivity matrices with lower errors in pair-wise correlation than non-censored and censored time series based on a standard processing pipeline. In addition, seed-based correlation analyses showed improved delineation of the default mode network. These demonstrate that the method can effectively reduce the adverse effects of motion in fMRI analysis.
Wallace, Tess, Tobias Kober, Jason Stockmann, Jonathan Polimeni, Simon Warfield, and Onur Afacan. 2022. “Real-time shimming with FID navigators”. Magn Reson Med. https://doi.org/10.1002/mrm.29421.
PURPOSE: To implement a method for real-time field control using rapid FID navigator (FIDnav) measurements and evaluate the efficacy of the proposed approach for mitigating dynamic field perturbations and improving T 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality. METHODS: FIDnavs were embedded in a gradient echo sequence and a subject-specific linear calibration model was generated on the scanner to facilitate rapid shim updates in response to measured FIDnav signals. To confirm the accuracy of FID-navigated field updates, phantom and volunteer scans were performed with online updates of the scanner B0 shim settings. To evaluate improvement in T 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality with real-time shimming, 10 volunteers were scanned at 3T while performing deep-breathing and nose-touching tasks designed to modulate the B0 field. Quantitative image quality metrics were compared with and without FID-navigated field control. An additional volunteer was scanned at 7T to evaluate performance at ultra-high field. RESULTS: Applying measured FIDnav shim updates successfully compensated for applied global and linear field offsets in phantoms and across all volunteers. FID-navigated real-time shimming led to a substantial reduction in field fluctuations and a consequent improvement in T 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality in volunteers performing deep-breathing and nose-touching tasks, with 7.57% ± 6.01% and 8.21% ± 10.90% improvement in peak SNR and structural similarity, respectively. CONCLUSION: FIDnavs facilitate rapid measurement and application of field coefficients for slice-wise B0 shimming. The proposed approach can successfully counteract spatiotemporal field perturbations and substantially improves T 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality, which is important for a variety of clinical and research applications, particularly at ultra-high field.