Journal Papers

F

Akhondi-Asl, Alireza, Onur Afacan, Mukund Balasubramanian, Robert Mulkern, and Simon Warfield. 2016. “Fast Myelin Water Fraction Estimation Using 2D Multislice CPMG”. Magn Reson Med 76 (4): 1301-13. https://doi.org/10.1002/mrm.26034.
PURPOSE: T2 relaxometry based on multiexponential fitting to a single slice multiecho sequence has been the most common MRI technique for myelin water fraction mapping, where the short T2 is associated with myelin water. However, very long acquisition times and physically unrealistic models for T2 distribution are limitations of this approach. We present a novel framework for myelin imaging which substantially increases the imaging speed and myelin water fraction estimation accuracy. METHOD: We used the 2D multislice Carr-Purcell-Meiboom-Gill sequence to increase the volume coverage. To compensate for nonideal slice profiles, we numerically solved the Bloch equations for a range of T2 and B1 inhomogeneity scales to construct the bases for the estimation of the T2 distribution. We used a finite mixture of continuous parametric distributions to describe the complete T2 spectrum and used the constrained variable projection optimization algorithm to estimate myelin water fraction. To validate our model, synthetic, phantom, and in vivo brain experiments were conducted. RESULTS: Using the Bloch equations, we can model the slice profile and construct the forward model of the T2 curve. Our method estimated myelin water fraction with smaller error than the nonnegative least squares algorithm. CONCLUSIONS: The proposed framework can be used for reliable whole brain myelin imaging with a resolution of 2×2×4  mm3 in ≈17  min. Magn Reson Med 76:1301-1313, 2016. © 2015 Wiley Periodicals, Inc.
Sui, Yao, Onur Afacan, Ali Gholipour, and Simon Warfield. (2021) 2021. “Fast and High-Resolution Neonatal Brain MRI Through Super-Resolution Reconstruction From Acquisitions With Variable Slice Selection Direction”. Front Neurosci 15: 636268. https://doi.org/10.3389/fnins.2021.636268.
The brain of neonates is small in comparison to adults. Imaging at typical resolutions such as one cubic mm incurs more partial voluming artifacts in a neonate than in an adult. The interpretation and analysis of MRI of the neonatal brain benefit from a reduction in partial volume averaging that can be achieved with high spatial resolution. Unfortunately, direct acquisition of high spatial resolution MRI is slow, which increases the potential for motion artifact, and suffers from reduced signal-to-noise ratio. The purpose of this study is thus that using super-resolution reconstruction in conjunction with fast imaging protocols to construct neonatal brain MRI images at a suitable signal-to-noise ratio and with higher spatial resolution than can be practically obtained by direct Fourier encoding. We achieved high quality brain MRI at a spatial resolution of isotropic 0.4 mm with 6 min of imaging time, using super-resolution reconstruction from three short duration scans with variable directions of slice selection. Motion compensation was achieved by aligning the three short duration scans together. We applied this technique to 20 newborns and assessed the quality of the images we reconstructed. Experiments show that our approach to super-resolution reconstruction achieved considerable improvement in spatial resolution and signal-to-noise ratio, while, in parallel, substantially reduced scan times, as compared to direct high-resolution acquisitions. The experimental results demonstrate that our approach allowed for fast and high-quality neonatal brain MRI for both scientific research and clinical studies.

E

Machado-Rivas, Fedel, Camilo Jaimes, Benoit Scherrer, Leslie Benson, Mark Gorman, Simon Warfield, and Onur Afacan. 2022. “Evaluation of White Matter Microstructure in Pediatric Onset Multiple Sclerosis With Diffusion Compartment Imaging”. J Neuroimaging. https://doi.org/10.1111/jon.13038.
BACKGROUND AND PURPOSE: Pediatric-onset multiple sclerosis (POMS) shows earlier axonal involvement and greater axonal loss than in adults. We aim to characterize the white matter (WM) microstructural changes in POMS using a diffusion compartment imaging (DCI) model and compare it to standard diffusion tensor imaging (DTI). METHODS: Eleven patients (2 males, mean age 18.8 ± 3.9 years) with a diagnosis of relapsing and remitting POMS (mean age at disease onset 13.8 ± 2.9 years, mean duration 5.1 ± 1.9 years) and healthy controls (8 males, mean age 26.4 ± 6.5 years) were recruited and imaged at 3 T. A 90-gradient set Cube and Sphere acquisition and a novel DCI model known as DIstribution of Anisotropic MicrOstructural eNvironments with Diffusion-weighted imaging (DIAMOND) were used to calculate a single anisotropic compartment, an isotropic compartment, and a free diffusion compartment. Lesions and contralateral normal-appearing white matter (NAWM) in patients and whole brain WM for controls were labeled. RESULTS: Eleven patients and 11 controls were recruited. When comparing lesions and contralateral NAWM in patients using DCI, compartmental axial diffusivity, radial diffusivity (cRD), and mean diffusivity (cMD) were higher in lesions. Conversely, compartmental fractional anisotropy (cFA) and heterogeneity index were lower in lesions. An analysis of DTI equivalents showed the same trends. In whole-brain NAWM of patients compared to controls, cRD and cMD were higher and cFA was lower in patients. CONCLUSION: Lesions in POMS can be accurately characterized by a DCI model. Incipient changes in NAWM seen in DCI may not be readily observable by DTI.
Afacan, Onur, Burak Erem, Diona Roby, Noam Roth, Amir Roth, Sanjay Prabhu, and Simon Warfield. (2016) 2016. “Evaluation of Motion and Its Effect on Brain Magnetic Resonance Image Quality in Children”. Pediatr Radiol 46 (12): 1728-35. https://doi.org/10.1007/s00247-016-3677-9.
BACKGROUND: Motion artifacts pose significant problems for the acquisition of MR images in pediatric populations. OBJECTIVE: To evaluate temporal motion metrics in MRI scanners and their effect on image quality in pediatric populations in neuroimaging studies. MATERIALS AND METHODS: We report results from a large pediatric brain imaging study that shows the effect of motion on MRI quality. We measured motion metrics in 82 pediatric patients, mean age 13.4 years, in a T1-weighted brain MRI scan. As a result of technical difficulties, 5 scans were not included in the subsequent analyses. A radiologist graded the images using a 4-point scale ranging from clinically non-diagnostic because of motion artifacts to no motion artifacts. We used these grades to correlate motion parameters such as maximum motion, mean displacement from a reference point, and motion-free time with image quality. RESULTS: Our results show that both motion-free time (as a ratio of total scan time) and average displacement from a position at a fixed time (when the center of k-space was acquired) were highly correlated with image quality, whereas maximum displacement was not as good a predictor. Among the 77 patients whose motion was measured successfully, 17 had average displacements of greater than 0.5 mm, and 11 of those (14.3%) resulted in non-diagnostic images. Similarly, 14 patients (18.2%) had less than 90% motion-free time, which also resulted in non-diagnostic images. CONCLUSION: We report results from a large pediatric study to show how children and young adults move in the MRI scanner and the effect that this motion has on image quality. The results will help the motion-correction community in better understanding motion patterns in pediatric populations and how these patterns affect MR image quality.
Ortinau, Cynthia, Caitlin Rollins, Ali Gholipour, Hyuk Jin Yun, Mackenzie Marshall, Borjan Gagoski, Onur Afacan, et al. 2019. “Early-Emerging Sulcal Patterns Are Atypical in Fetuses With Congenital Heart Disease”. Cereb Cortex 29 (8): 3605-16. https://doi.org/10.1093/cercor/bhy235.
Fetuses with congenital heart disease (CHD) have third trimester alterations in cortical development on brain magnetic resonance imaging (MRI). However, the intersulcal relationships contributing to global sulcal pattern remain unknown. This study applied a novel method for examining the geometric and topological relationships between sulci to fetal brain MRIs from 21-30 gestational weeks in CHD fetuses (n = 19) and typically developing (TD) fetuses (n = 17). Sulcal pattern similarity index (SI) to template fetal brain MRIs was determined for the position, area, and depth for corresponding sulcal basins and intersulcal relationships for each subject. CHD fetuses demonstrated altered global sulcal patterns in the left hemisphere compared with TD fetuses (TD [SI, mean ± SD]: 0.822 ± 0.023, CHD: 0.795 ± 0.030, P = 0.002). These differences were present in the earliest emerging sulci and were driven by differences in the position of corresponding sulcal basins (TD: 0.897 ± 0.024, CHD: 0.878 ± 0.019, P = 0.006) and intersulcal relationships (TD: 0.876 ± 0.031, CHD: 0.857 ± 0.018, P = 0.033). No differences in cortical gyrification index, mean curvature, or surface area were present. These data suggest our methods may be more sensitive than traditional measures for evaluating cortical developmental alterations early in gestation.

D

Wallace, Tess, Jonathan Polimeni, Jason Stockmann, Scott Hoge, Tobias Kober, Simon Warfield, and Onur Afacan. 2021. “Dynamic Distortion Correction for Functional MRI Using FID Navigators”. Magn Reson Med 85 (3): 1294-1307. https://doi.org/10.1002/mrm.28505.
PURPOSE: To develop a method for slice-wise dynamic distortion correction for EPI using rapid spatiotemporal B0 field measurements from FID navigators (FIDnavs) and to evaluate the efficacy of this new approach relative to an established data-driven technique. METHODS: A low-resolution reference image was used to create a forward model of FIDnav signal changes to enable estimation of spatiotemporal B0 inhomogeneity variations up to second order from measured FIDnavs. Five volunteers were scanned at 3 T using a 64-channel coil with FID-navigated EPI. The accuracy of voxel shift measurements and geometric distortion correction was assessed for experimentally induced magnetic field perturbations. The temporal SNR was evaluated in EPI time-series acquired at rest and with a continuous nose-touching action, before and after image realignment. RESULTS: Field inhomogeneity coefficients and voxel shift maps measured using FIDnavs were in excellent agreement with multi-echo EPI measurements. The FID-navigated distortion correction accurately corrected image geometry in the presence of induced magnetic field perturbations, outperforming the data-driven approach in regions with large field offsets. In functional MRI scans with nose touching, FIDnav-based correction yielded temporal SNR gains of 30% in gray matter. Following image realignment, which accounted for global image shifts, temporal SNR gains of 3% were achieved. CONCLUSIONS: Our proposed application of FIDnavs enables slice-wise dynamic distortion correction with high temporal efficiency. We achieved improved signal stability by leveraging the encoding information from multichannel coils. This approach can be easily adapted to other EPI-based sequences to improve temporal SNR for a variety of clinical and research applications.

C

Kurugol, Sila, Onur Afacan, Richard Lee, Catherine Seager, Michael Ferguson, Deborah Stein, Reid Nichols, et al. (2020) 2020. “Correction To: Prospective Pediatric Study Comparing Glomerular Filtration Rate Estimates Based on Motion-Robust Dynamic Contrast-Enhanced Magnetic Resonance Imaging and Serum Creatinine (eGFR) to Tc DTPA”. Pediatr Radiol 50 (5): 755-56. https://doi.org/10.1007/s00247-020-04654-9.
The originally published version of this article contained a typographical error. In the text under the subheading "Dynamic contrast-enhanced MRI method, post-processing, and MR-GFR calculation" and in Table 1 the intravenous injection rate of gadobutrol was incorrectly listed as 0.2 mL/s.

B

Coll-Font, Jaume, Onur Afacan, Jeanne Chow, Richard Lee, Alto Stemmer, Simon Warfield, and Sila Kurugol. 2020. “Bulk Motion-Compensated DCE-MRI for Functional Imaging of Kidneys in Newborns”. J Magn Reson Imaging 52 (1): 207-16. https://doi.org/10.1002/jmri.27021.
BACKGROUND: Evaluation of kidney function in newborns with hydronephrosis is important for clinical decisions. Dynamic contrast-enhanced (DCE) MRI can provide the necessary anatomical and functional information. Golden angle dynamic radial acquisition and compressed sensing reconstruction provides sufficient spatiotemporal resolution to achieve accurate parameter estimation for functional imaging of kidneys. However, bulk motion during imaging (rigid or nonrigid movement of the subject resulting in signal dropout) remains an unresolved challenge. PURPOSE: To evaluate a motion-compensated (MoCo) DCE-MRI technique for robust evaluation of kidney function in newborns. Our method includes: 1) motion detection, 2) motion-robust image reconstruction, 3) joint realignment of the volumes, and 4) tracer-kinetic (TK) model fitting to evaluate kidney function parameters. STUDY TYPE: Retrospective. SUBJECTS: Eleven newborn patients (ages <6 months, 6 female). FIELD STRENGTH/SEQUENCE: 3T; dynamic "stack-of-stars" 3D fast low-angle shot (FLASH) sequence using a multichannel body-matrix coil. ASSESSMENT: We evaluated the proposed technique in terms of the signal-to-noise ratio (SNR) of the reconstructed images, the presence of discontinuities in the contrast agent concentration time curves due to motion with a total variation (TV) metric and the goodness of fit of the TK model, and the standard variation of its parameters. STATISTICAL TESTS: We used a paired t-test to compare the MoCo and no-MoCo results. RESULTS: The proposed MoCo method successfully detected motion and improved the SNR by 3.3 (P = 0.012) and decreased TV by 0.374 (P = 0.017) across all subjects. Moreover, it decreased nRMSE of the TK model fit for the subjects with less than five isolated bulk motion events in 6 minutes (mean 1.53, P = 0.043), but not for the subjects with more frequent events or no motion (P = 0.745 and P = 0.683). DATA CONCLUSION: Our results indicate that the proposed MoCo technique improves the image quality and accuracy of the TK model fit for subjects who present isolated bulk motion events. LEVEL OF EVIDENCE: 3 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;52:207-216.
Hesper, Tobias, Evgeny Bulat, Sarah Bixby, Alireza Akhondi-Asl, Onur Afacan, Patricia Miller, Garrett Bowen, Simon Warfield, and Young-Jo Kim. (2017) 2017. “Both 3-T DGEMRIC and Acetabular-Femoral T2 Difference May Detect Cartilage Damage at the Chondrolabral Junction”. Clin Orthop Relat Res 475 (4): 1058-65. https://doi.org/10.1007/s11999-016-5136-1.
BACKGROUND: In addition to case reports of gadolinium-related toxicities, there are increasing theoretical concerns about the use of gadolinium for MR imaging. As a result, there is increasing interest in noncontrast imaging techniques for biochemical cartilage assessment. Among them, T2 mapping holds promise because of its simplicity, but its biophysical interpretation has been controversial. QUESTIONS/PURPOSES: We sought to determine whether (1) 3-T delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) and T2 mapping are both capable of detecting cartilage damage at the chondrolabral junction in patients with femoroacetabular impingement (FAI); and (2) whether there is a correlation between these two techniques for acetabular and femoral head cartilage assessment. METHODS: Thirty-one patients with hip-related symptoms resulting from FAI underwent a preoperative 3-T MRI of their hip that included dGEMRIC and T2 mapping (symptomatic group, 16 women, 15 men; mean age, 27 ± 8 years). Ten volunteers with no symptoms according to the WOMAC served as a control (asymptomatic group, seven women, three men; mean age, 28 ± 3 years). After morphologic cartilage assessment, acetabular and femoral head cartilages were graded according to the modified Outerbridge grading criteria. In the midsagittal plane, single-observer analyses of precontrast T1 values (volunteers), the dGEMRIC index (T1Gd, patients), and T2 mapping values (everyone) were compared in acetabular and corresponding femoral head cartilage at the chondrolabral junction of each hip by region-of-interest analysis. RESULTS: In the symptomatic group, T1Gd and T2 values were lower in the acetabular cartilage compared with corresponding femoral head cartilage (T1Gd: 515 ± 165 ms versus 650 ± 191 ms, p < 0.001; T2: 39 ± 8 ms versus 46 ± 10 ms, p < 0.001). In contrast, the asymptomatic group demonstrated no differences in T1 and T2 values for the acetabular and femoral cartilages with the numbers available (T1: 861 ± 130 ms versus 860 ± 182 ms, p = 0.98; T2: 43 ± 7 ms versus 42 ± 6 ms, p = 0.73). No correlation with the numbers available was noted between the modified Outerbridge grade and T1, T1Gd, or T2 as well as between T2 and either T1 or T1Gd. CONCLUSIONS: Without the need for contrast media application, T2 mapping may be a viable alternative to dGEMRIC when assessing hip cartilage at the chondrolabral junction. However, acquisition-related phenomena as well as regional variations in the microstructure of hip cartilage necessitate an internal femoral head cartilage control when interpreting these results. LEVEL OF EVIDENCE: Level IV, diagnostic study.
Mitsouras, Dimitris, Robert Mulkern, Onur Afacan, Dana Brooks, and Frank Rybicki. (2007) 2007. “Basis Function Cross-Correlations for Robust K-Space Sample Density Compensation, With Application to the Design of Radiofrequency Excitations”. Magn Reson Med 57 (2): 338-52. https://doi.org/10.1002/mrm.21125.
The problem of k-space sample density compensation is restated as the normalization of the independent information that can be expressed by the ensemble of Fourier basis functions corresponding to the trajectory. Specifically, multiple samples (complex exponential functions) may be contributing to each independent information element (independent basis function). Normalization can be accomplished by solving a linear system based on the cross-correlation matrix of the underlying Fourier basis functions. The solution to this system is straightforward and can be obtained without resorting to discretization since the cross-correlations of Fourier basis functions are analytically known. Furthermore, no restrictions are placed on the k-space trajectory and its point-spread function. Additionally, the linear system can be used to elucidate key trade-offs involved in k-space trajectory design. The approach can be used to compensate samples acquired for image reconstruction or designed for low flip angle radiofrequency (RF) excitation. Here it is demonstrated for the latter application, using reversed spiral trajectories. In this case the linear system approach enables one to easily incorporate additional constraints such as smoothness to the solution. For typical RF excitation durations (<20 ms) it is shown that density compensation can even be achieved without numerical iteration.