Search

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., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Classic galactosemia (OMIM# 230400) is an autosomal recessive disorder due to galactose-1-phosphate uridyltransferase deficiency. Newborn screening and prompt treatment with a galactose-free diet prevent the severe consequences of galactosemia, but clinical outcomes remain suboptimal. Five men and five women with classic galactosemia (mean age = 27.2 ± 5.47 years) received comprehensive neurological and neuropsychological evaluations, electroencephalogram (EEG) and magnetic resonance imaging (MRI). MRI data from nine healthy controls (mean age = 30.22 ± 3.52 years) were used for comparison measures. Galactosemia subjects experienced impaired memory, language processing, visual-motor skills, and increased anxiety. Neurological examinations revealed tremor and dysarthria in six subjects. In addition, there was ataxia in three subjects and six subjects had abnormal gait. Mean full scale IQ was 80.4 ± 17.3. EEG evaluations revealed right-sided abnormalities in five subjects and bilateral abnormalities in one subject. Compared to age- and gender-matched controls, subjects with galactosemia had reduced volume in left cerebellum white matter, bilateral putamen, and left superior temporal sulcus. Galactosemia patients also had lower fractional anisotropy and higher radial diffusivity values in the dorsal and ventral language networks compared to the controls. Furthermore, there were significant correlations between neuropsychological test results and the T1 volume and diffusivity scalars. Our findings help to identify anatomic correlates to motor control, learning and memory, and language in subjects with galactosemia. The results from this preliminary assessment may provide insights into the pathophysiology of this inborn error of metabolism., (© 2020 SSIEM.)

In MRI practice, it is inevitable to appropriately balance between image resolution, signal-to-noise ratio (SNR), and scan time. It has been shown that super-resolution reconstruction (SRR) is effective to achieve such a balance, and has obtained better results than direct high-resolution (HR) acquisition, for certain contrasts and sequences. The focus of this work was on constructing images with spatial resolution higher than can be practically obtained by direct Fourier encoding. A novel learning approach was developed, which was able to provide an estimate of the spatial gradient prior from the low-resolution (LR) inputs for the HR reconstruction. By incorporating the anisotropic acquisition schemes, the learning model was trained over the LR images themselves only. The learned gradients were integrated as prior knowledge into a gradient-guided SRR model. A closed-form solution to the SRR model was developed to obtain the HR reconstruction. Our approach was assessed on the simulated data as well as the data acquired on a Siemens 3T MRI scanner containing 45 MRI scans from 15 subjects. The experimental results demonstrated that our approach led to superior SRR over state-of-the-art methods, and obtained better images at lower or the same cost in scan time than direct HR acquisition.

The third trimester of pregnancy is a period of rapid development of fiber bundles in the fetal white matter. Using a recently developed motion-tracked slice-to-volume registration (MT-SVR) method, we aimed to quantify tract-specific developmental changes in apparent diffusion coefficient (ADC), fractional anisotropy (FA), and volume in third trimester healthy fetuses. To this end, we reconstructed diffusion tensor images from motion corrected fetal diffusion magnetic resonance imaging data. With an approved protocol, fetal MRI exams were performed on healthy pregnant women at 3 Tesla and included multiple (2-8) diffusion scans of the fetal head (1-2 b = 0 s/mm 2 images and 12 diffusion-sensitized images at b = 500 s/mm 2 ). Diffusion data from 32 fetuses (13 females) with median gestational age (GA) of 33 weeks 4 days were processed with MT-SVR and deterministic tractography seeded by regions of interest corresponding to 12 major fiber tracts. Multivariable regression analysis was used to evaluate the association of GA with volume, FA, and ADC for each tract. For all tracts, the volume and FA increased, and the ADC decreased with GA. Associations reached statistical significance for: FA and ADC of the forceps major; volume and ADC for the forceps minor; FA, ADC, and volume for the cingulum; ADC, FA, and volume for the uncinate fasciculi; ADC of the inferior fronto-occipital fasciculi, ADC of the inferior longitudinal fasciculi; and FA and ADC for the corticospinal tracts. These quantitative results demonstrate the complex pattern and rates of tract-specific, GA-related microstructural changes of the developing white matter in human fetal brain., (© 2020 The Authors. Human Brain Mapping published by Wiley Periodicals, Inc.)

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., (© 2019 International Society for Magnetic Resonance in Medicine.)

Background: The right ventricle (RV) often fails when functioning as the systemic ventricle, but the cause is not understood. We tested the hypothesis that myofiber organization is abnormal in the failing systemic right ventricle., Methods: We used diffusion-weighted cardiovascular magnetic resonance imaging to examine 3 failing hearts explanted from young patients with a systemic RV and one structurally normal heart with postnatally acquired RV hypertrophy for comparison. Diffusion compartment imaging was computed to separate the free diffusive component representing free water from an anisotropic component characterizing the orientation and diffusion characteristics of myofibers. The orientation of each anisotropic compartment was displayed in glyph format and used for qualitative description of myofibers and for construction of tractograms. The helix angle was calculated across the ventricular walls in 5 locations and displayed graphically. Scalar parameters (fractional anisotropy and mean diffusivity) were compared among specimens., Results: The hypertrophied systemic RV has an inner layer, comprising about 2/3 of the wall, composed of hypertrophied trabeculae and an epicardial layer of circumferential myofibers. Myofibers within smaller trabeculae are aligned and organized with parallel fibers while larger, composite bundles show marked disarray, largely between component trabeculae. We observed a narrow range of helix angles in the outer, compact part of the wall consistent with aligned, approximately circumferential fibers. However, there was marked variation of helix angle in the inner, trabecular part of the wall consistent with marked variation in fiber orientation. The apical whorl was disrupted or incomplete and we observed myocardial whorls or vortices at other locations. Fractional anisotropy was lower in abnormal hearts while mean diffusivity was more variable, being higher in 2 but lower in 1 heart, compared to the structurally normal heart., Conclusions: Myofiber organization is abnormal in the failing systemic RV and might be an important substrate for heart failure and arrhythmia. It is unclear if myofiber disorganization is due to hemodynamic factors, developmental problems, or both.

Background: Current methods to estimate glomerular filtration rate (GFR) have shortcomings. Estimates based on serum creatinine are known to be inaccurate in the chronically ill and during acute changes in renal function. Gold standard methods such as inulin and 99m Tc diethylenetriamine pentaacetic acid (DTPA) require blood or urine sampling and thus can be difficult to perform in children. Motion-robust radial volumetric interpolated breath-hold examination (VIBE) dynamic contrast-enhanced MRI represents a novel tool for estimating GFR that has not been validated in children., Objective: The purpose of our study was to determine the feasibility and accuracy of GFR measured by motion-robust radial VIBE dynamic contrast-enhanced MRI compared to estimates by serum creatinine (eGFR) and 99m Tc DTPA in children., Materials and Methods: We enrolled children, 0-18 years of age, who were undergoing both a contrast-enhanced MRI and nuclear medicine 99m Tc DTPA glomerular filtration rate (NM-GFR) within 2 weeks of each other. Enrolled children consented to an additional 6-min dynamic contrast-enhanced MRI scan using the motion-robust high spatiotemporal resolution prototype dynamic radial VIBE sequence (Siemens, Erlangen, Germany) at 3 tesla (T). The images were reconstructed offline with high temporal resolution (~3 s/volume) using compressed sensing image reconstruction including regularization in temporal dimension to improve image quality and reduce streaking artifacts. Images were then automatically post-processed using in-house-developed software. Post-processing steps included automatic segmentation of kidney parenchyma and aorta using convolutional neural network techniques and tracer kinetic model fitting using the Sourbron two-compartment model to calculate the MR-based GFR (MR-GFR). The NM-GFR was compared to MR-GFR and estimated GFR based on serum creatinine (eGFR) using Pearson correlation coefficient and Bland-Altman analysis., Results: Twenty-one children (7 female, 14 male) were enrolled between February 2017 and May 2018. Data from six of these children were not further analyzed because of deviations from the MRI protocol. Fifteen patients were analyzed (5 female, 10 male; average age 5.9 years); the method was technically feasible in all children. The results showed that the MR-GFR correlated with NM-GFR with a Pearson correlation coefficient (r-value) of 0.98. Bland-Altman analysis (i.e. difference of MR-GFR and NM-GFR versus mean of NM-GFR and MR-GFR) showed a mean difference of -0.32 and reproducibility coefficient of 18 with 95% confidence interval, and the coefficient of variation of 6.7% with values between -19 (-1.96 standard deviation) and 18 (+1.96 standard deviation). In contrast, serum creatinine compared with NM-GFR yielded an r-value of 0.73. Bland-Altman analysis (i.e. difference of eGFR and NM-GFR versus mean of NM-GFR and eGFR) showed a mean difference of 2.9 and reproducibility coefficient of 70 with 95% confidence interval, and the coefficient of variation of 25% with values between -67 (-1.96 standard deviation) and 73 (+1.96 standard deviation)., Conclusion: MR-GFR is a technically feasible and reliable method of measuring GFR when compared to the reference standard, NM-GFR by serum 99m Tc DTPA, and MR-GFR is more reliable than estimates based on serum creatinine.

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., (© 2020 by the American Society of Neuroimaging.)

Objective: To describe a technique for performing magnetic resonance urogram (MRU) in infants without sedation or anesthesia., Methods: Eighteen infants underwent MRU in the absence of sedating medications using a 'feed and wrap' technique (FW-MRU). Dynamic contrast enhanced images were obtained. Dynamic radial VIBE and compressed sensing image reconstruction were used to correct for motion artifact., Results: Seventeen of the 18 patients had successful FW-MRU. Feed and wrap' magnetic resonance urogram provided high-quality anatomic and functional renal data., Conclusion: Initial experience with FW-MRU demonstrates it to be a promising anesthesia-free modality for obtaining anatomic and functional imaging of the urinary tract in infants., (Copyright © 2019 Journal of Pediatric Urology Company. All rights reserved.)

Purpose: To measure spatiotemporal B 0 field changes in real time using FID navigators (FIDnavs) and to demonstrate the efficacy of retrospectively correcting high-resolution T 2 * -weighted images using a novel FIDnav framework., Methods: A forward model of the complex FIDnav signals was generated by simulating the effect of changes in the underlying B 0 inhomogeneity coefficients, with spatial encoding provided by a multi-channel reference image. Experiments were performed at 3T to assess the accuracy of B 0 field estimates from FIDnavs acquired from a 64-channel head coil under different shim settings and in 5 volunteers performing deep-breathing and nose-touching tasks designed to modulate the B 0 field. Second-order, in-plane spherical harmonic (SH) inhomogeneity coefficients estimated from FIDnavs were incorporated into an iterative reconstruction to retrospectively correct 2D gradient-echo images acquired in both axial and sagittal planes., Results: Spatiotemporal B 0 field changes measured from rapidly acquired FIDnavs were in good agreement with the results of second-order SH fitting to the measured field maps. FIDnav field estimates accounted for a significant proportion of the ΔB 0 variance induced by deep breathing (64 ± 21%) and nose touching (67 ± 34%) across all volunteers. Ghosting, blurring, and intensity modulation artifacts in T 2 * -weighted images, induced by spatiotemporal field changes, were visibly reduced following retrospective correction with FIDnav inhomogeneity coefficients., Conclusions: Spatially resolved B 0 inhomogeneity changes up to second order can be characterized in real time using the proposed approach. Retrospective FIDnav correction substantially improves T 2 * -weighted image quality in the presence of strong B 0 field modulations, with potential for real-time shimming., (© 2019 International Society for Magnetic Resonance in Medicine.)

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., (© 2019 International Society for Magnetic Resonance in Medicine.)

In this work, we proposed a novel image-based MRI super-resolution reconstruction (SRR) approach based on anisotropic acquisition schemes. We achieved superior reconstruction to state-of-the-art work by introducing a new multi-scale gradient field prior that guides the reconstruction of the high-resolution (HR) image. The prior improves both spatial smoothness and edge preservation. The inverse of the forward model of image formation is used to propagate the gradient guidance from the low-resolution (LR) images to the HR image space. The gradient fields over multiple scales were exploited for more accurate edge localization in the reconstruction. The proposed SRR allows inter-volume motion during the MRI scans and can incorporate with the LR images with arbitrary orientations and displacements in the frequency space, such as orthogonal and origin-shifted scans. The proposed approach was evaluated on the synthetic data as well as the data acquired on a Siemens 3T MRI scanner containing 45 MRI scans from 14 subjects. The evaluation results demonstrate that our proposed prior leads to improved SRR as compared to state-of-the-art priors, and that the proposed SRR obtains better results at lower or the same cost in scan time than direct HR acquisition. In particular, the anatomical structures of hippocampus can be clearly shown in our reconstructed images. This is a significant improvement for the in vivo studies of the hippocampus.

Early identification of kidney function deterioration is essential to determine which newborn patients with dilation of the renal pelvis (hydronephrosis) should undergo surgery. Kidney function can be measured by fitting a tracer kinetic (TK) model onto a series of Dynamic Contrast Enhanced (DCE) MR images and deriving the glomerular filtration rate (GFR) from the TK model. Unfortunately, heavy breathing and large bulk motion events create outliers and misalignments that introduce large errors in the TK estimates. Moreover, aligning the series of DCE images is not trivial due to the contrast differences between them and the undersampling artifacts due to fast imaging. We present a bulk motion detection and a linear time invariant (LTI) model-based motion correction approach for DCE-MRI alignment that leverages the temporal dynamics of the DCE data at each voxel. We evaluate our approach on 10 newborn patients that underwent DCE imaging without sedation. For each patient, we reconstructed the sequence of DCE images, detected and removed the volumes corrupted by motion using a self navigation approach, aligned the sequence using our approach and fitted the TK model to compute GFR. The results show that our approach correctly aligned all volumes and improved the TK model fit and, on average, reducing the normalized root-mean-squared error by 0.17.

Fetal magnetic resonance imaging (MRI) has been gaining increasing interest in both clinical radiology and research. Echoplanar imaging (EPI) offers a unique potential, as it can be used to acquire images very fast. It can be used to freeze motion, or to get multiple images with various contrast mechanisms that allow studying the microstructure and function of the fetal brain and body organs. In this article, we discuss the current clinical and research applications of fetal EPI. This includes T2*-weighted imaging to better identify blood products and vessels, using diffusion-weighted MRI to investigate connections of the developing brain and using functional MRI (fMRI) to identify the functional networks of the developing brain. EPI can also be used as an alternative structural sequence when banding or standing wave artifacts adversely affect the mainstream sequences used routinely in structural fetal MRI. We also discuss the challenges with EPI acquisitions, and potential solutions. As EPI acquisitions are inherently sensitive to susceptibility artifacts, geometric distortions limit the use of high-resolution EPI acquisitions. Also, interslice motion and transmit and receive field inhomogeneities may create significant artifacts in fetal EPI. We conclude by discussing promising research directions to overcome these challenges to improve the use of EPI in clinical and research applications.

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., (© The Author(s) 2018. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Objective: Diffusion tensor imaging (DTI) of the white matter is a biomarker for neurological disease burden in tuberous sclerosis complex (TSC). To clarify the basis of abnormal diffusion in TSC, we correlated ex vivo high-resolution diffusion imaging with histopathology in four tissue types: cortex, tuber, perituber, and white matter., Methods: Surgical specimens of three children with TSC were scanned in a 3T or 7T MRI with a structural image isotropic resolution of 137-300 micron, and diffusion image isotropic resolution of 270-1,000 micron. We stained for myelin (luxol fast blue, LFB), gliosis (glial fibrillary acidic protein, GFAP), and neurons (NeuN) and registered the digitized histopathology slides (0.686 micron resolution) to MRI for visual comparison. We then performed colocalization analysis in four tissue types in each specimen. Finally, we applied a linear mixed model (LMM) for pooled analysis across the three specimens., Results: In white matter and perituber regions, LFB optical density measures correlated with fractional anisotropy (FA) and inversely with mean diffusivity (MD). In white matter only, GFAP correlated with MD, and inversely with FA. In tubers and in the cortex, there was little variation in mean LFB and GFAP signal intensity, and no correlation with MRI metrics. Neuronal density correlated with MD. In the analysis of the combined specimens, the most robust correlation was between white matter MD and LFB metrics., Interpretation: In TSC, diffusion imaging abnormalities in microscopic tissue types correspond to specific histopathological markers. Across all specimens, white matter diffusivity correlates with myelination., (© 2019 The Authors. Annals of Clinical and Translational Neurology published by Wiley Periodicals, Inc on behalf of American Neurological Association.)

Purpose: To achieve motion-robust diffusion compartment imaging (DCI) in near continuously moving subjects based on simultaneous multi-slice, diffusion-weighted brain MRI., Methods: Simultaneous multi-slice (SMS) acquisition enables fast and dense sampling of k- and q-space. We propose to achieve motion-robust DCI via slice-level motion correction by exploiting the rigid coupling between simultaneously acquired slices. This coupling provides 3D coverage of the anatomy that substantially constraints the slice-to-volume alignment problem. This is incorporated into an explicit model of motion dynamics that handles continuous and large subject motion in robust DCI reconstruction., Results: We applied the proposed technique, called Motion Tracking based on Simultanous Multislice Registration (MT-SMR) to multi b-value SMS diffusion-weighted brain MRI of healthy volunteers and motion-corrupted scans of 20 pediatric subjects. Quantitative and qualitative evaluation based on fractional anisotropy in unidirectional fiber regions, and DCI in crossing-fiber regions show robust reconstruction in the presence of motion., Conclusion: The proposed approach has the potential to extend routine use of SMS DCI in very challenging populations, such as young children, newborns, and non-cooperative patients., (© 2018 International Society for Magnetic Resonance in Medicine.)

This study used a swine model of mildly hypothermic prolonged circulatory arrest and found that the addition of 2.4% inhaled hydrogen gas to inspiratory gases during and after the ischemic insult significantly decreased neurologic and renal injury compared with controls. With proper precautions, inhalational hydrogen may be administered safely through conventional ventilators and may represent a complementary therapy that can be easily incorporated into current workflows. In the future, inhaled hydrogen may diminish the sequelae of ischemia that occurs in congenital heart surgery, cardiac arrest, extracorporeal life-support events, acute myocardial infarction, stroke, and organ transplantation.

Altered structural fetal brain development has been linked to neuro-developmental disorders. These structural alterations can be potentially detected in utero using diffusion tensor imaging (DTI). However, acquisition and reconstruction of in utero fetal brain DTI remains challenging. Until now, motion-robust DTI methods have been employed for reconstruction of in utero fetal DTIs. However, due to the unconstrained fetal motion and permissible in utero acquisition times, these methods yielded limited success and have typically resulted in noisy DTIs. Consequently, atlases and methods that could enable groupwise studies, multi-modality imaging, and computer-aided diagnosis from in utero DTIs have not yet been developed. This paper presents the first DTI atlas of the fetal brain computed from in utero diffusion-weighted images. For this purpose an algorithm for computing an unbiased spatiotemporal DTI atlas, which integrates kernel-regression in age with a diffeomorphic tensor-to-tensor registration of motion-corrected and reconstructed individual fetal brain DTIs, was developed. Our new algorithm was applied to a set of 67 fetal DTI scans acquired from healthy fetuses each scanned at a gestational age between 21 and 39 weeks. The neurodevelopmental trends in the fetal brain, characterized by the atlas, were qualitatively and quantitatively compared with the observations reported in prior ex vivo and in utero studies, and with results from imaging gestational-age equivalent preterm infants. Our major findings revealed early presence of limbic fiber bundles, followed by the appearance and maturation of projection pathways (characterized by an age related increase in FA) during late 2nd and early 3rd trimesters. During the 3rd trimester association fiber bundles become evident. In parallel with the appearance and maturation of fiber bundles, from 21 to 39 gestational weeks gradual disappearance of the radial coherence of the telencephalic wall was qualitatively identified. These results and analyses show that our DTI atlas of the fetal brain is useful for reliable detection of major neuronal fiber bundle pathways and for characterization of the fetal brain reorganization that occurs in utero. The atlas can also serve as a useful resource for detection of normal and abnormal fetal brain development in utero., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Purpose: To develop a novel framework for rapid, intrinsic head motion measurement in MRI using FID navigators (FIDnavs) from a multichannel head coil array., Methods: FIDnavs encode substantial rigid-body motion information; however, current implementations require patient-specific training with external tracking data to extract quantitative positional changes. In this work, a forward model of FIDnav signals was calibrated using simulated movement of a reference image within a model of the spatial coil sensitivities. A FIDnav module was inserted into a nonselective 3D FLASH sequence, and rigid-body motion parameters were retrospectively estimated every readout time using nonlinear optimization to solve the inverse problem posed by the measured FIDnavs. This approach was tested in simulated data and in 7 volunteers, scanned at 3T with a 32-channel head coil array, performing a series of directed motion paradigms., Results: FIDnav motion estimates achieved mean absolute errors of 0.34 ± 0.49 mm and 0.52 ± 0.61° across all subjects and scans, relative to ground-truth motion measurements provided by an electromagnetic tracking system. Retrospective correction with FIDnav motion estimates resulted in substantial improvements in quantitative image quality metrics across all scans with intentional head motion., Conclusions: Quantitative rigid-body motion information can be effectively estimated using the proposed FIDnav-based approach, which represents a practical method for retrospective motion compensation in less cooperative patient populations., (© 2018 International Society for Magnetic Resonance in Medicine.)

Diffusion tensor imaging (DTI) based group analysis has helped uncover the impact of white matter injuries in a wide range of studies involving subjects from preterm neonates to adults. The application of these methods to fetal cohorts, however, has been hampered by the challenging nature of in utero fetal DTI caused by unconstrained fetal motion, limited scan times, and limited signal-to-noise ratio. We present a framework that addresses these issues to systematically evaluate group differences in fetal cohorts. A motion-robust DTI computation approach with a new unbiased DTI template construction method is unified with kernel-regression in age and tensor-specific registration to normalize DTI volumes in an unbiased space. A robust statistical approach is used to map region-specific group differences to the medial representation of the tracts of interest. The proposed approach was applied and showed, for the first time, differences in local white matter fractional anisotropy based on in utero DTI of fetuses with congenital heart disease and age-matched healthy controls. This paper suggests the need for fetal-specific pipelines to be used for DTI-based group analysis involving fetal cohorts.

The development of interventions to prevent congenital Zika syndrome (CZS) has been limited by the lack of an established nonhuman primate model. Here we show that infection of female rhesus monkeys early in pregnancy with Zika virus (ZIKV) recapitulates many features of CZS in humans. We infected 9 pregnant monkeys with ZIKV, 6 early in pregnancy (weeks 6-7 of gestation) and 3 later in pregnancy (weeks 12-14 of gestation), and compared findings with uninfected controls. 100% (6 of 6) of monkeys infected early in pregnancy exhibited prolonged maternal viremia and fetal neuropathology, including fetal loss, smaller brain size, and histopathologic brain lesions, including microcalcifications, hemorrhage, necrosis, vasculitis, gliosis, and apoptosis of neuroprogenitor cells. High-resolution MRI demonstrated concordant lesions indicative of deep gray matter injury. We also observed spinal, ocular, and neuromuscular pathology. Our data show that vascular compromise and neuroprogenitor cell dysfunction are hallmarks of CZS pathogenesis, suggesting novel strategies to prevent and to treat this disease., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Diffusion weighted magnetic resonance imaging, or DWI, is one of the most promising tools for the analysis of neural microstructure and the structural connectome of the human brain. The application of DWI to map early development of the human connectome in-utero, however, is challenged by intermittent fetal and maternal motion that disrupts the spatial correspondence of data acquired in the relatively long DWI acquisitions. Fetuses move continuously during DWI scans. Reliable and accurate analysis of the fetal brain structural connectome requires careful compensation of motion effects and robust reconstruction to avoid introducing bias based on the degree of fetal motion. In this paper we introduce a novel robust algorithm to reconstruct in-vivo diffusion-tensor MRI (DTI) of the moving fetal brain and show its effect on structural connectivity analysis. The proposed algorithm involves multiple steps of image registration incorporating a dynamic registration-based motion tracking algorithm to restore the spatial correspondence of DWI data at the slice level and reconstruct DTI of the fetal brain in the standard (atlas) coordinate space. A weighted linear least squares approach is adapted to remove the effect of intra-slice motion and reconstruct DTI from motion-corrected data. The proposed algorithm was tested on data obtained from 21 healthy fetuses scanned in-utero at 22-38 weeks gestation. Significantly higher fractional anisotropy values in fiber-rich regions, and the analysis of whole-brain tractography and group structural connectivity, showed the efficacy of the proposed method compared to the analyses based on original data and previously proposed methods. The results of this study show that slice-level motion correction and robust reconstruction is necessary for reliable in-vivo structural connectivity analysis of the fetal brain. Connectivity analysis based on graph theoretic measures show high degree of modularity and clustering, and short average characteristic path lengths indicative of small-worldness property of the fetal brain network. These findings comply with previous findings in newborns and a recent study on fetuses. The proposed algorithm can provide valuable information from DWI of the fetal brain not available in the assessment of the original 2D slices and may be used to more reliably study the developing fetal brain connectome., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Quantitative body DW-MRI can detect abdominal abnormalities as well as monitor response-to-therapy for applications including cancer and inflammatory bowel disease with increased accuracy. Parameter estimates are obtained by fitting a forward model of DW-MRI signal decay to the observed data acquired with several b-values. The DW-MRI signal decay models typically used do not account for respiratory, cardiac and peristaltic motion, however, which may deteriorate the accuracy and robustness of parameter estimates. In this work, we introduce a new model of DW-MRI signal decay that explicitly accounts for motion. Specifically, we estimated motion-compensated model parameters by simultaneously solving image registration and model estimation (SIR-ME) problems utilizing the interdependence of acquired volumes along the diffusion-weighting dimension. To accomplish this, we applied the SIR-ME model to the in-vivo DW-MRI data sets of 26 Crohn's disease (CD) patients and achieved improved precision of the estimated parameters by reducing the coefficient of variation by 8%, 24% and 8% for slow diffusion (D), fast diffusion (D*) and fast diffusion fraction (f) parameters respectively, compared to parameters estimated with independent registration in normal-appearing bowel regions. Moreover, the parameters estimated with the SIR-ME model reduced the error rate in classifying normal and abnormal bowel loops to 12% for D and 10% for f parameter with a reduction in error rate by 13% and 11% for D and f parameters, respectively, compared to the error rate in classifying parameter estimates obtained with independent registration. The experiments in DW-MRI of liver in 20 subjects also showed that the SIR-ME model improved the precision of parameter estimation by reducing the coefficient of variation to 7% for D, 23% for D*, and 8% for the f parameter. Using the SIR-ME model, the coefficient of variation was reduced by 4%, 14% and 6% for D, D* and f parameters, respectively, compared to parameters estimated with independent registration. These results demonstrate that the proposed SIR-ME model improves the accuracy and robustness of quantitative body DW-MRI in characterizing tissue microstructure., (Copyright © 2017 Elsevier B.V. All rights reserved.)

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 (T1 Gd , 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, T1 Gd and T2 values were lower in the acetabular cartilage compared with corresponding femoral head cartilage (T1 Gd : 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, T1 Gd , or T2 as well as between T2 and either T1 or T1 Gd ., 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.

Longitudinal characterization of early brain growth in-utero has been limited by a number of challenges in fetal imaging, the rapid change in size, shape and volume of the developing brain, and the consequent lack of suitable algorithms for fetal brain image analysis. There is a need for an improved digital brain atlas of the spatiotemporal maturation of the fetal brain extending over the key developmental periods. We have developed an algorithm for construction of an unbiased four-dimensional atlas of the developing fetal brain by integrating symmetric diffeomorphic deformable registration in space with kernel regression in age. We applied this new algorithm to construct a spatiotemporal atlas from MRI of 81 normal fetuses scanned between 19 and 39 weeks of gestation and labeled the structures of the developing brain. We evaluated the use of this atlas and additional individual fetal brain MRI atlases for completely automatic multi-atlas segmentation of fetal brain MRI. The atlas is available online as a reference for anatomy and for registration and segmentation, to aid in connectivity analysis, and for groupwise and longitudinal analysis of early brain growth.

In this work, we introduce a novel motion-robust spatially constrained parameter estimation (MOSCOPE) technique for kidney diffusion-weighted MRI. The proposed motion compensation technique does not require a navigator, trigger, or breath-hold but only uses the intrinsic features of the acquired data to track and compensate for motion to reconstruct precise models of the renal diffusion signal. We have developed a technique for physiological motion tracking based on robust state estimation and sequential registration of diffusion sensitized slices acquired within 200ms. This allows a sampling rate of 5Hz for state estimation in motion tracking that is sufficiently faster than both respiratory and cardiac motion rates in children and adults, which range between 0.8 to 0.2Hz, and 2.5 to 1Hz, respectively. We then apply the estimated motion parameters to data from each slice and use motion-compensated data for 1) robust intra-voxel incoherent motion (IVIM) model estimation in the kidney using a spatially constrained model fitting approach, and 2) robust weighted least squares estimation of the diffusion tensor model. Experimental results, including precision of IVIM model parameters using bootstrap-sampling and in-vivo whole kidney tractography, showed significant improvement in precision and accuracy of these models using the proposed method compared to models based on the original data and volumetric registration.

Purpose: To evaluate the feasibility of using diffusion-weighted magnetic resonance imaging (DW-MRI) to assess the fetal lung apparent diffusion coefficient (ADC) at 3 Tesla (T)., Materials and Methods: Seventy-one pregnant women (32 second trimester, 39 third trimester) were scanned with a twice-refocused Echo-planar diffusion-weighted imaging sequence with 6 different b-values in 3 orthogonal diffusion orientations at 3T. After each scan, a region-of-interest (ROI) mask was drawn to select a region in the fetal lung and an automated robust maximum likelihood estimation algorithm was used to compute the ADC parameter. The amount of motion in each scan was visually rated., Results: When scans with unacceptable levels of motion were eliminated, the lung ADC values showed a strong association with gestational age (P < 0.01), increasing dramatically between 16 and 27 weeks and then achieving a plateau around 27 weeks., Conclusion: We show that to get reliable estimates of ADC values of fetal lungs, a multiple b-value acquisition, where motion is either corrected or considered, can be performed. J. Magn. Reson. Imaging 2016;44:1650-1655., (© 2016 International Society for Magnetic Resonance in Medicine.)

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., Competing Interests: Compliance with ethical standards: Conflicts of interest None

Simultaneous multi-slice (SMS) echo-planar imaging has had a huge impact on the acceleration and routine use of diffusion-weighted MRI (DWI) in neuroimaging studies in particular the human connectome project; but also holds the potential to facilitate DWI of moving subjects, as proposed by the new technique developed in this paper. We present a novel registration-based motion tracking technique that takes advantage of the multi-plane coverage of the anatomy by simultaneously acquired slices to enable robust reconstruction of neural microstructure from SMS DWI of moving subjects. Our technique constitutes three main components: 1) motion tracking and estimation using SMS registration, 2) detection and rejection of intra-slice motion, and 3) robust reconstruction. Quantitative results from 14 volunteer subject experiments and the analysis of motion-corrupted SMS DWI of 6 children indicate robust reconstruction in the presence of continuous motion and the potential to extend the use of SMS DWI in very challenging populations.

This work proposes a novel approach for motion-robust diffusion-weighted (DW) brain MRI reconstruction through tracking temporal head motion using slice-to-volume registration. The slice-level motion is estimated through a filtering approach that allows tracking the head motion during the scan and correcting for out-of-plane inconsistency in the acquired images. Diffusion-sensitized image slices are registered to a base volume sequentially over time in the acquisition order where an outlier-robust Kalman filter, coupled with slice-to-volume registration, estimates head motion parameters. Diffusion gradient directions are corrected for the aligned DWI slices based on the computed rotation parameters and the diffusion tensors are directly estimated from the corrected data at each voxel using weighted linear least squares. The method was evaluated in DWI scans of adult volunteers who deliberately moved during scans as well as clinical DWI of 28 neonates and children with different types of motion. Experimental results showed marked improvements in DWI reconstruction using the proposed method compared to the state-of-the-art DWI analysis based on volume-to-volume registration. This approach can be readily used to retrieve information from motion-corrupted DW imaging data.

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., (© 2015 Wiley Periodicals, Inc.)

Purpose: Current diagnosis of fetal posterior fossa anomalies by sonography and conventional MRI is limited by fetal position, motion, and by two-dimensional (2D), rather than three-dimensional (3D), representation. In this study, we aimed to validate the use of a novel magnetic resonance imaging (MRI) technique, 3D super-resolution motion-corrected MRI, to image the fetal posterior fossa., Methods: From a database of pregnant women who received fetal MRIs at our institution, images of 49 normal fetal brains were reconstructed. Six measurements of the cerebellum, vermis, and pons were obtained for all cases on 2D conventional and 3D reconstructed MRI, and the agreement between the two methods was determined using concordance correlation coefficients. Concordance of axial and coronal measurements of the transcerebellar diameter was also assessed within each method., Results: Between the two methods, the concordance of measurements was high for all six structures (P < .001), and was highest for larger structures such as the transcerebellar diameter. Within each method, agreement of axial and coronal measurements of the transcerebellar diameter was superior in 3D reconstructed MRI compared to 2D conventional MRI (P < .001)., Conclusions: This comparison study validates the use of 3D super-resolution motion-corrected MRI for imaging the fetal posterior fossa, as this technique results in linear measurements that have high concordance with 2D conventional MRI measurements. Lengths of the transcerebellar diameter measured within a 3D reconstruction are more concordant between imaging planes, as they correct for fetal motion and orthogonal slice acquisition. This technique will facilitate further study of fetal abnormalities of the posterior fossa., (Copyright © 2016 by the American Society of Neuroimaging.)

Quantitative diffusion-weighted MR imaging (DW-MRI) of the body enables characterization of the tissue microenvironment by measuring variations in the mobility of water molecules. The diffusion signal decay model parameters are increasingly used to evaluate various diseases of abdominal organs such as the liver and spleen. However, previous signal decay models (i.e., mono-exponential, bi-exponential intra-voxel incoherent motion (IVIM) and stretched exponential models) only provide insight into the average of the distribution of the signal decay rather than explicitly describe the entire range of diffusion scales. In this work, we propose a probability distribution model of incoherent motion that uses a mixture of Gamma distributions to fully characterize the multi-scale nature of diffusion within a voxel. Further, we improve the robustness of the distribution parameter estimates by integrating spatial homogeneity prior into the probability distribution model of incoherent motion (SPIM) and by using the fusion bootstrap solver (FBM) to estimate the model parameters. We evaluated the improvement in quantitative DW-MRI analysis achieved with the SPIM model in terms of accuracy, precision and reproducibility of parameter estimation in both simulated data and in 68 abdominal in-vivo DW-MRIs. Our results show that the SPIM model not only substantially reduced parameter estimation errors by up to 26%; it also significantly improved the robustness of the parameter estimates (paired Student's t-test, p < 0.0001) by reducing the coefficient of variation (CV) of estimated parameters compared to those produced by previous models. In addition, the SPIM model improves the parameter estimates reproducibility for both intra- (up to 47%) and inter-session (up to 30%) estimates compared to those generated by previous models. Thus, the SPIM model has the potential to improve accuracy, precision and robustness of quantitative abdominal DW-MRI analysis for clinical applications., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Purpose: To compare and evaluate the use of super-resolution reconstruction (SRR), in frequency, image, and wavelet domains, to reduce through-plane partial voluming effects in magnetic resonance imaging., Methods: The reconstruction of an isotropic high-resolution image from multiple thick-slice scans has been investigated through techniques in frequency, image, and wavelet domains. Experiments were carried out with thick-slice T2-weighted fast spin echo sequence on the Academic College of Radiology MRI phantom, where the reconstructed images were compared to a reference high-resolution scan using peak signal-to-noise ratio (PSNR), structural similarity image metric (SSIM), mutual information (MI), and the mean absolute error (MAE) of image intensity profiles. The application of super-resolution reconstruction was then examined in retrospective processing of clinical neuroimages of ten pediatric patients with tuberous sclerosis complex (TSC) to reduce through-plane partial voluming for improved 3D delineation and visualization of thin radial bands of white matter abnormalities., Results: Quantitative evaluation results show improvements in all evaluation metrics through super-resolution reconstruction in the frequency, image, and wavelet domains, with the highest values obtained from SRR in the image domain. The metric values for image-domain SRR versus the original axial, coronal, and sagittal images were PSNR = 32.26 vs 32.22, 32.16, 30.65; SSIM = 0.931 vs 0.922, 0.924, 0.918; MI = 0.871 vs 0.842, 0.844, 0.831; and MAE = 5.38 vs 7.34, 7.06, 6.19. All similarity metrics showed high correlations with expert ranking of image resolution with MI showing the highest correlation at 0.943. Qualitative assessment of the neuroimages of ten TSC patients through in-plane and out-of-plane visualization of structures showed the extent of partial voluming effect in a real clinical scenario and its reduction using SRR. Blinded expert evaluation of image resolution in resampled out-of-plane views consistently showed the superiority of SRR compared to original axial and coronal image acquisitions., Conclusions: Thick-slice 2D T2-weighted MRI scans are part of many routine clinical protocols due to their high signal-to-noise ratio, but are often severely affected by through-plane partial voluming effects. This study shows that while radiologic assessment is performed in 2D on thick-slice scans, super-resolution MRI reconstruction techniques can be used to fuse those scans to generate a high-resolution image with reduced partial voluming for improved postacquisition processing. Qualitative and quantitative evaluation showed the efficacy of all SRR techniques with the best results obtained from SRR in the image domain. The limitations of SRR techniques are uncertainties in modeling the slice profile, density compensation, quantization in resampling, and uncompensated motion between scans.

Magnetic resonance (MR) imaging provides a unique in-vivo capability of visualizing tissue in the human brain non-invasively, which has tremendously improved patient care over the past decades. However, there are still prominent artifacts, such as intensity inhomogeneities due to the use of an array of receiving coils (RC) to measure the MR signal or noise amplification due to accelerated imaging strategies. It is critical to mitigate these artifacts for both visual inspection and quantitative analysis. The cornerstone to address this issue pertains to the knowledge of coil sensitivity profiles (CSP) of the RCs, which describe how the measured complex signal decays with the distance to the RC. Existing methods for CSP estimation share a number of limitations: (i) they primarily focus on CSP magnitude, while it is known that the solution to the MR image reconstruction problem involves complex CSPs and (ii) they only provide point estimates of the CSPs, which makes the task of optimizing the parameters and acquisition protocol for their estimation difficult. In this paper, we propose a novel statistical framework for estimating complex-valued CSPs. We define a CSP estimator that uses spatial smoothing and additional body coil data for phase normalization. The main contribution is to provide detailed information on the statistical distribution of the CSP estimator, which yields automatic determination of the optimal degree of smoothing for ensuring minimal bias and provides guidelines to the optimal acquisition strategy.

Objective: To determine normative ranges for fetal ocular biometrics between 19 and 38 weeks gestational age (GA) using volumetric MRI reconstruction., Method: The 3D images of 114 healthy fetuses between 19 and 38 weeks GA were created using super-resolution volume reconstructions from MRI slice acquisitions. These 3D images were semi-automatically segmented to measure fetal orbit volume, binocular distance (BOD), interocular distance (IOD), and ocular diameter (OD)., Results: All biometry correlated with GA (Volume, Pearson's correlation coefficient (CC) = 0.9680; BOD, CC = 0.9552; OD, CC = 0.9445; and IOD, CC = 0.8429), and growth curves were plotted against linear and quadratic growth models. Regression analysis showed quadratic models to best fit BOD, IOD, and OD and a linear model to best fit volume., Conclusion: Orbital volume had the greatest correlation with GA, although BOD and OD also showed strong correlation. The normative data found in this study may be helpful for the detection of congenital fetal anomalies with more consistent measurements than are currently available. © 2015 John Wiley & Sons, Ltd., (© 2015 John Wiley & Sons, Ltd.)

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.

Non-invasive characterization of water molecule's mobility variations by quantitative analysis of diffusion-weighted MRI (DW-MRI) signal decay in the abdomen has the potential to serve as a biomarker in gastrointestinal and oncological applications. Accurate and reproducible estimation of the signal decay model parameters is challenging due to the presence of respiratory, cardiac, and peristalsis motion. Independent registration of each b-value image to the b-value=0 s/mm(2) image prior to parameter estimation might be sub-optimal because of the low SNR and contrast difference between images of varying b-value. In this work, we introduce a motion-compensated parameter estimation framework that simultaneously solves image registration and model estimation (SIR-ME) problems by utilizing the interdependence of acquired volumes along the diffusion weighting dimension. We evaluated the improvement in model parameters estimation accuracy using 16 in-vivo DW-MRI data sets of Crohn's disease patients by comparing parameter estimates obtained using the SIR-ME model to the parameter estimates obtained by fitting the signal decay model to the acquired DW-MRI images. The proposed SIR-ME model reduced the average root-mean-square error between the observed signal and the fitted model by more than 50%. Moreover, the SIR-ME model estimates discriminate between normal and abnormal bowel loops better than the standard parameter estimates.

Acquisition of a series of anisotropically oversampled acquisitions (so-called anisotropic "snapshots") and reconstruction in the image space has recently been proposed to increase the spatial resolution in diffusion weighted imaging (DWI), providing a theoretical 8x acceleration at equal signal-to-noise ratio (SNR) compared to conventional dense k-space sampling. However, in most works, each DW image is reconstructed separately and the fact that the DW images constitute different views of the same anatomy is ignored. In addition, current approaches are limited by their inability to reconstruct a high resolution (HR) acquisition from snapshots with different subsets of diffusion gradients: an isotropic HR gradient image cannot be reconstructed if one .of its anisotropic snapshots is missing, for example due to intra-scan motion, even if other snapshots for this gradient were successfully acquired. In this work, we propose a novel multi-snapshot DWI reconstruction technique that simultaneously achieves HR reconstruction and local tissue model estimation while enabling reconstruction from snapshots containing different subsets of diffusion gradients, providing increased robustness to patient motion and potential for acceleration. Our approach is formalized as a joint probabilistic model with missing observations, from which interactions between missing snapshots, HR reconstruction and a generic tissue model naturally emerge. We evaluate our approach with synthetic simulations, simulated multi-snapshot scenario and in vivo multi-snapshot imaging. We show that (1) our combined approach ultimately provides both better HR reconstruction and better tissue model estimation and (2) the error in the case of missing snapshots can be quantified. Our novel multi-snapshot technique will enable improved high spatial characterization of the brain connectivity and microstructure in vivo.

Fetal magnetic resonance imaging (MRI) examinations have become well-established procedures at many institutions and can serve as useful adjuncts to ultrasound (US) exams when diagnostic doubts remain after US. Due to fetal motion, however, fetal MRI exams are challenging and require the MR scanner to be used in a somewhat different mode than that employed for more routine clinical studies. Herein we review the techniques most commonly used, and those that are available, for fetal MRI with an emphasis on the physics of the techniques and how to deploy them to improve success rates for fetal MRI exams. By far the most common technique employed is single-shot T2-weighted imaging due to its excellent tissue contrast and relative immunity to fetal motion. Despite the significant challenges involved, however, many of the other techniques commonly employed in conventional neuro- and body MRI such as T1 and T2*-weighted imaging, diffusion and perfusion weighted imaging, as well as spectroscopic methods remain of interest for fetal MR applications. An effort to understand the strengths and limitations of these basic methods within the context of fetal MRI is made in order to optimize their use and facilitate implementation of technical improvements for the further development of fetal MR imaging, both in acquisition and post-processing strategies.

Quantitative assessment of myelin density in the white matter is an emerging tool for neurodegenerative disease related studies such as multiple sclerosis and Schizophrenia. For the last two decades, T2 relaxometry based on multi-exponential fitting to a single slice multi-echo sequence has been the most common MRI technique for myelin water fraction (MWF) mapping, where the short T2 is associated with myelin water. However, modeling the spectrum of the relaxations as the sum of large number of impulse functions with unknown amplitudes makes the accuracy and robustness of the estimated MWF's questionable. In this paper, we introduce a novel model with small number of parameters to simultaneously characterize transverse relaxation rate spectrum and B1 inhomogeneity at each voxel. We use mixture of three Wald distributions with unknown mixture weights, mean and shape parameters to represent the distribution of the relative amount of water in between myelin sheets, tissue water, and cerebrospinal fluid. The parameters of the model are estimated using the variable projection method and are used to extract the MWF at each voxel. In addition, we use Extended Phase Graph (EPG) method to compensate for the stimulated echoes caused by B1 inhomogeneity. To validate our model, synthetic and real brain experiments were conducted where we have compared our novel algorithm with the non-negative least squares (NNLS) as the state-of-the-art technique in the literature. Our results indicate that we can estimate MWF map with substantially higher accuracy as compared to the NNLS method.

Diffusion-weighted MRI images acquired with multiple b-values have the potential to improve diagnostic accuracy by increasing the conspicuity of lesions and inflammatory activity with background suppression. Unfortunately, the inherently low signal-to-noise ratio (SNR) of DW-MRI reduces enthusiasm for using these images for diagnostic purposes. Moreover, lengthy acquisition times limit our ability to improve the quality of multi b-value DW-MRI images by multiple excitations acquisition and signal averaging at each b-value. To offset these limitations, we propose the Simultaneous Model Estimation and Image Reconstruction (SMEIR) for DW-MRI, which substantially improves the quality of multi b-value DW-MRI images without increasing acquisition times. Our model introduces the physiological signal decay model of DW-MRI as a constraint in the reconstruction of the DW-MRI images. An in-vivo experiment using 6 low-quality DW-MRI datasets of a healthy subject showed that SMEIR reconstruction of low-quality data improved SNR by 55% in the liver and by 41% in the kidney without increasing acquisition times. We also demonstrated the clinical impact of our SMEIR reconstruction by increasing the conspicuity of inflamed bowel regions in DW-MRI of 12 patients with Crohn's disease. The contrast-to-noise ratio (CNR) of the inflamed regions in the SMEIR images was higher by 12.6% relative to CNR in the original DW-MRI images.

The desire to understand complex mental processes using functional MRI drives development of imaging techniques that scan the whole human brain at a high spatial and temporal resolution. In this work, an accelerated multishot three-dimensional echo-planar imaging sequence is proposed to increase the temporal resolution of these studies. A combination of two modern acceleration techniques, UNFOLD and GRAPPA is used in the secondary phase encoding direction to reduce the scan time effectively. The sequence (repetition time of 1.02 s) was compared with standard two-dimensional echo-planar imaging (3 s) and multishot three-dimensional echo-planar imaging (3 s) sequences with both block design and event-related functional MRI paradigms. With the same experimental setup and imaging time, the temporal resolution improvement with our sequence yields similar activation regions in the block design functional MRI paradigm with slightly increased t-scores. Moreover, additional information on the timing of rapid dynamic changes was extracted from accelerated images for the case of the event related complex mental paradigm., (Copyright © 2011 Wiley Periodicals, Inc.)

Awake animal imaging is becoming an important tool in behavioral neuroscience and preclinical drug discovery. Non-invasive ultra-high-field, functional magnetic resonance imaging (fMRI) provides a window to the mind, making it possible to image changes in brain activity across distributed, integrated neural circuits with high temporal and spatial resolution. In theory, changes in brain function, anatomy, and chemistry can be recorded in the same animal from early life into old age under stable or changing environmental conditions. This prospective capability of animal imaging to follow changes in brain neurobiology after genetic or environmental insult has great value to the fields of psychiatry and neurology and probably stands as the key advantage of MRI over other methods in the neuroscience toolbox. In addition, awake animal imaging offers the ability to record signal changes across the entire brain in seconds. When combined with the use of 3D segmented, annotated, brain atlases, and computational analysis, it is possible to reconstruct distributed, integrated neural circuits or 'fingerprints' of brain activity. These fingerprints can be used to characterize the activity and function of new psychotherapeutics in preclinical development and to study the neurobiology of integrated neural circuits controlling cognition and emotion. In this review, we describe the methods used to image awake animals and the recent advances in the radiofrequency electronics, pulse sequences, and the development of 3D segmented atlases and software for image analysis. Results from pharmacological MRI studies and from studies using provocation paradigms to elicit emotional responses are provided as a small sample of the number of different applications possible with awake animal imaging.

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., (Copyright (c) 2007 Wiley-Liss, Inc.)