Diffusion magnetic resonance imaging (dMRI) is pivotal for probing the microstructure of the rapidly-developing fetal brain. However, fetal motion during scans and its interaction with magnetic field inhomogeneities result in artifacts and data scattering across spatial and angular domains. The effects of those artifacts are more pronounced in high-angular resolution fetal dMRI, where signal-to-noise ratio is very low. Those effects lead to biased estimates and compromise the consistency and reliability of dMRI analysis. This work presents High Angular resolution diffusion Imaging reconsTruction and Correction approacH (HAITCH), the first and the only publicly available tool to correct and reconstruct multi-shell high-angular resolution fetal dMRI data. HAITCH offers several technical advances that include a blip-reversed dual-echo acquisition for dynamic distortion correction, advanced motion correction formodel-freeand robust reconstruction, optimized multi-shell design for enhanced information capture and increased tolerance to motion, and outlier detection for improved reconstruction fidelity. The framework is open-source, flexible, and can be used to process any type of fetal dMRI data, including single-echo or single-shell acquisitions, but is most effective when used with multi-shell multi-echo fetal dMRI data that cannot be processed with any of the existing tools. Validation experiments on real fetal dMRI scans demonstrate significant improvements and accurate correction across diverse fetal ages and motion levels. HAITCH successfully removes artifacts and reconstructs high-fidelity fetal dMRI data suitable for advanced diffusion modeling, including fiber orientation distribution function estimation. These advancements pave the way for more reliable analysis of the fetal brain microstructure and tractography under challenging imaging conditions.
Journal Papers
2025
Diffusion-weighted imaging (DWI) during MR enterography helps identify bowel inflammation in Crohn's disease (CD). However, image quality is compromised by geometric distortions from B0 field variations and physiological motion, making it challenging for radiologists to correlate findings between DWI and structural images. Traditional correction methods using reversed polarity scans are ineffective due to motion between acquisitions, which limits accurate estimation of intravoxel incoherent motion (IVIM) parameters. We propose a dual-echo echo-planar imaging (EPI) method that retrospectively corrects both geometric distortions and motion in 3T bowel DWI by accounting for field changes during peristalsis and breathing. We added a 5- to 7-min dual-echo EPI DW sequence (eight b-values, six directions) to the clinical MR enterography protocol of 21 patients with suspected CD at 3T MRI. Distortion correction was applied based on dynamically estimated fields from dual-echo DWI, followed by intra-volume registration between odd-even slices and inter-volume registration for motion correction. Two experienced board-certified radiologists evaluated the severity of the disease using simplified magnetic resonance index of activity (MaRIA) scores. Based on their consensus scores, patients were categorized into three groups: no active disease (MaRIA score = 0), active disease (MaRIA score = 1-2), and severe disease (MaRIA score = 3-5). The proposed DWI correction pipeline improved DWI/T2-weighted image Dice similarity from 0.73 to 0.89, enabling better correlation of findings between structural and DW-MR images and enhancing DWI's clinical value. Corrected IVIM parameters showed stronger correlations with MaRIA scores (D: ρ = -0.93; f: ρ = -0.94, p < 0.001) compared to uncorrected parameters (D: ρ = -0.68, p = 0.001; f: ρ = -0.35, p = 0.118). Diagnostic sensitivity increased from 0.44 to 0.89, while parameter uncertainty decreased from 35.58% to 19.31% for D and 63.48% to 40.40% for f (p < 0.001). These improvements strengthen quantitative IVIM imaging for CD assessment, potentially reducing reliance on contrast imaging while offering enhanced tissue perfusion and diffusion insights.
2024
BACKGROUND: A lack of in utero imaging data hampers our understanding of the connections in the human fetal brain. Generalizing observations from postmortem subjects and premature newborns is inaccurate due to technical and biological differences.
PURPOSE: To evaluate changes in fetal brain structural connectivity between 23 and 35 weeks postconceptional age using a spatiotemporal atlas of diffusion tensor imaging (DTI).
STUDY TYPE: Retrospective.
POPULATION: Publicly available diffusion atlases, based on 60 healthy women (age 18-45 years) with normal prenatal care, from 23 and 35 weeks of gestation.
FIELD STRENGTH/SEQUENCE: 3.0 Tesla/DTI acquired with diffusion-weighted echo planar imaging (EPI).
ASSESSMENT: We performed whole-brain fiber tractography from DTI images. The cortical plate of each diffusion atlas was segmented and parcellated into 78 regions derived from the Edinburgh Neonatal Atlas (ENA33). Connectivity matrices were computed, representing normalized fiber connections between nodes. We examined the relationship between global efficiency (GE), local efficiency (LE), small-worldness (SW), nodal efficiency (NE), and betweenness centrality (BC) with gestational age (GA) and with laterality.
STATISTICAL TESTS: Linear regression was used to analyze changes in GE, LE, NE, and BC throughout gestation, and to assess changes in laterality. The t-tests were used to assess SW. P-values were corrected using Holm-Bonferroni method. A corrected P-value <0.05 was considered statistically significant.
RESULTS: Network analysis revealed a significant weekly increase in GE (5.83%/week, 95% CI 4.32-7.37), LE (5.43%/week, 95% CI 3.63-7.25), and presence of SW across GA. No significant hemisphere differences were found in GE (P = 0.971) or LE (P = 0.458). Increasing GA was significantly associated with increasing NE in 41 nodes, increasing BC in 3 nodes, and decreasing BC in 2 nodes.
DATA CONCLUSION: Extensive network development and refinement occur in the second and third trimesters, marked by a rapid increase in global integration and local segregation.
LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY: Stage 2.
We utilized motion-corrected diffusion tensor imaging (DTI) to evaluate microstructural changes in healthy fetal brains during the late second and third trimesters. Data were derived from fetal magnetic resonance imaging scans conducted as part of a prospective study spanning from 2013 March to 2019 May. The study included 44 fetuses between the gestational ages (GAs) of 23 and 36 weeks. We reconstructed fetal brain DTI using a motion-tracked slice-to-volume registration framework. Images were segmented into 14 regions of interest (ROIs) through label propagation using a fetal DTI atlas, with expert refinement. Statistical analysis involved assessing changes in fractional anisotropy (FA) and mean diffusivity (MD) throughout gestation using mixed-effects models, and identifying points of change in trajectory for ROIs with nonlinear trends. Results showed significant GA-related changes in FA and MD in all ROIs except in the thalamus' FA and corpus callosum's MD. Hemispheric asymmetries were found in the FA of the periventricular white matter (pvWM), intermediate zone, and subplate and in the MD of the ganglionic eminence and pvWM. This study provides valuable insight into the normal patterns of development of MD and FA in the fetal brain. These changes are closely linked with cytoarchitectonic changes and display indications of early functional specialization.
PURPOSE: To develop a self-navigated motion compensation strategy for 3D radial MRI that can compensate for continuous head motion by measuring rigid body motion parameters with high temporal resolution from the central k-space acquisition point (self-encoded FID navigator) in each radial spoke.
METHODS: A forward model was created from low-resolution calibration data to simulate the effect of relative motion between the coil sensitivity profiles and the underlying object on the self-encoded FID navigator signal. Trajectory deviations were included in the model as low spatial-order field variations. Three volunteers were imaged at 3 T using a modified 3D gradient-echo sequence acquired with a Kooshball trajectory while performing abrupt and continuous head motion. Rigid body-motion parameters were estimated from the central k-space signal of each spoke using a least-squares fitting algorithm. The accuracy of self-navigated motion parameters was assessed relative to an established external tracking system. Quantitative image quality metrics were computed for images with and without retrospective correction using external and self-navigated motion measurements.
RESULTS: Self-encoded FID navigators achieved mean absolute errors of 0.69 ± 0.82 mm and 0.73 ± 0.87° relative to external tracking for maximum motion amplitudes of 12 mm and 10°. Retrospective correction of the 3D radial data resulted in substantially improved image quality for both abrupt and continuous motion paradigms, comparable to external tracking results.
CONCLUSIONS: Accurate rigid body motion parameters can be rapidly obtained from self-encoded FID navigator signals in 3D radial MRI to continuously correct for head movements. This approach is suitable for robust neuroanatomical imaging in subjects that exhibit patterns of large and frequent motion.
We introduce a generative model for synthesis of large scale 3D datasets for quantitative parameter mapping of myelin water fraction (MWF). Our model combines a MR physics signal decay model with an accurate probabilistic multi-component parametric T2 model. We synthetically generate a wide variety of high quality signals and corresponding parameters from a wide range of naturally occurring prior parameter values. To capture spatial variation, the generative signal decay model is combined with a generative spatial model conditioned on generic tissue segmentations. Synthesized 3D datasets can be used to train any convolutional neural network (CNN) based architecture for MWF estimation. Our source code is available at: https://github.com/quin-med-harvard-edu/synthmap Reduction of acquisition time at the expense of lower SNR, as well as accuracy and repeatability of MWF estimation techniques, are key factors that affect the adoption of MWF mapping in clinical practice. We demonstrate that the synthetically trained CNN provides superior accuracy over the competing methods under the constraints of naturally occurring noise levels as well as on the synthetically generated images at low SNR levels. Normalized root mean squared error (nRMSE) is less than 7% on synthetic data, which is significantly lower than competing methods. Additionally, the proposed method yields a coefficient of variation (CoV) that is at least 4x better than the competing method on intra-session test-retest reference dataset.
Succinic semialdehyde dehydrogenase deficiency (SSADHD) is an inherited metabolic disorder of γ-aminobutyrate (GABA) catabolism. Cerebral waste clearance along glymphatic perivascular spaces depends on aquaporin 4 (AQP4) water channels, the function of which was shown to be influenced by GABA. Sleep disturbances are associated independently with SSADHD and glymphatic dysfunction. This study aimed to determine whether indices of the hyperGABAergic state characteristic of SSADHD coincide with glymphatic dysfunction and sleep disturbances and to explicate the modulatory effect that GABA may have on the glymphatic system. The study included 42 individuals (21 with SSADHD; 21 healthy controls) who underwent brain MRIs and magnetic resonance spectroscopy (MRS) for assessment of glymphatic dysfunction and cortical GABA, plasma GABA measurements, and circadian clock gene expression. The SSADHD subjects responded to an additional Children's Sleep Habits Questionnaire (CSHQ). Compared with the control group, SSADHD subjects did not differ in sex and age but had a higher severity of enlarged perivascular spaces in the centrum semiovale (p < 0.001), basal ganglia (p = 0.01), and midbrain (p = 0.001), as well as a higher MRS-derived GABA/NAA peak (p < 0.001). Within the SSADHD group, the severity of glymphatic dysfunction was specific for a lower MRS-derived GABA/NAA (p = 0.04) and lower plasma GABA (p = 0.004). Additionally, the degree of their glymphatic dysfunction correlated with the CSHQ-estimated sleep disturbances scores (R = 5.18, p = 0.03). In the control group, EPVS burden did not correlate with age or cerebral and plasma GABA values. The modulatory effect that GABA may exert on the glymphatic system has therapeutic implications for sleep-related disorders and neurodegenerative conditions associated with glymphatic dysfunction.
Accurately measuring renal function is crucial for pediatric patients with kidney conditions. Traditional methods have limitations, but dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides a safe and efficient approach for detailed anatomical evaluation and renal function assessment. However, motion artifacts during DCE-MRI can degrade image quality and introduce misalignments, leading to unreliable results. This study introduces a motion-compensated reconstruction technique for DCE-MRI data acquired using golden-angle radial sampling. Our proposed method achieves three key objectives: (1) identifying and removing corrupted data (outliers) using a Gaussian process model fitting with a k -space center navigator, (2) efficiently clustering the data into motion phases and performing interphase registration, and (3) utilizing a novel formulation of motion-compensated radial reconstruction. We applied the proposed motion correction (MoCo) method to DCE-MRI data affected by varying degrees of motion, including both respiratory and bulk motion. We compared the outcomes with those obtained from the conventional radial reconstruction. Our evaluation encompassed assessing the quality of images, concentration curves, and tracer kinetic model fitting, and estimating renal function. The proposed MoCo reconstruction improved the temporal signal-to-noise ratio for all subjects, with a 21.8% increase on average, while total variation values of the aorta, right, and left kidney concentration were improved for each subject, with 32.5%, 41.3%, and 42.9% increases on average, respectively. Furthermore, evaluation of tracer kinetic model fitting indicated that the median standard deviation of the estimated filtration rate ( σ F T ), mean normalized root-mean-squared error (nRMSE), and chi-square goodness-of-fit of tracer kinetic model fit were decreased from 0.10 to 0.04, 0.27 to 0.24, and, 0.43 to 0.27, respectively. The proposed MoCo technique enabled more reliable renal function assessment and improved image quality for detailed anatomical evaluation in the case of bulk and respiratory motion during the acquisition of DCE-MRI.
OBJECTIVES: The T1-weighted GRE (gradient recalled echo) sequence with the Dixon technique for water/fat separation is an essential component of abdominal MRI (magnetic resonance imaging), useful in detecting tumors and characterizing hemorrhage/fat content. Unfortunately, the current implementation of this sequence suffers from several problems: (1) low resolution to maintain high pixel bandwidth and minimize chemical shift; (2) image blurring due to respiratory motion; (3) water/fat swapping due to the natural ambiguity between fat and water peaks; and (4) off-resonance fat blurring due to the multipeak nature of the fat spectrum. The goal of this study was to evaluate the image quality of water/fat separation using a high-resolution 3-point Dixon golden angle radial acquisition with retrospective motion compensation and multipeak fat modeling in children undergoing abdominal MRI.
MATERIALS AND METHODS: Twenty-two pediatric patients (4.2 ± 2.3 years) underwent abdominal MRI on a 3 T scanner with routine abdominal protocol and with a 3-point Dixon radial-VIBE (volumetric interpolated breath-hold examination) sequence. Field maps were calculated using 3D graph-cut optimization followed by fat and water calculation from k-space data by iteratively solving an optimization problem. A 6-peak fat model was used to model chemical shifts in k-space. Residual respiratory motion was corrected through soft-gating by weighting each projection based on the estimated respiratory motion from the center of the k-space. Reconstructed images were reviewed by 3 pediatric radiologists on a PACS (picture archiving and communication systems) workstation. Subjective image quality and water/fat swapping artifact were scored by each pediatric radiologist using a 5-point Likert scale. The VoL (variance of Laplacian) of the reconstructed images was used to objectively quantify image sharpness.
RESULTS: Based on the overall Likert scores, the images generated using the described method were significantly superior to those reconstructed by the conventional 2-point Dixon technique ( P < 0.05). Water/fat swapping artifact was observed in 14 of 22 patients using 2-point Dixon, and this artifact was not present when using the proposed method. Image sharpness was significantly improved using the proposed framework.
CONCLUSIONS: In smaller patients, a high-quality water/fat separation with sharp visualization of fine details is critical for diagnostic accuracy. High-resolution golden angle radial-VIBE 3-point Dixon acquisition with 6-peak fat model and soft-gated motion correction offers improved image quality at the expense of an additional 1-minute acquisition time. Thus, this technique offers the potential to replace the conventional 2-point Dixon technique.
2023
Complex engineered systems are often equipped with suites of sensors and ancillary devices that monitor their performance and maintenance needs. MRI scanners are no different in this regard. Some of the ancillary devices available to support MRI equipment, the ones of particular interest here, have the distinction of actually participating in the image acquisition process itself. Most commonly, such devices are used to monitor physiological motion or variations in the scanner's imaging fields, allowing the imaging and/or reconstruction process to adapt as imaging conditions change. "Classic" examples include electrocardiography (ECG) leads and respiratory bellows to monitor cardiac and respiratory motion, which have been standard equipment in scan rooms since the early days of MRI. Since then, many additional sensors and devices have been proposed to support MRI acquisitions. The main physical properties that they measure may be primarily "mechanical" (eg acceleration, speed, and torque), "acoustic" (sound and ultrasound), "optical" (light and infrared), or "electromagnetic" in nature. A review of these ancillary devices, as currently available in clinical and research settings, is presented here. In our opinion, these devices are not in competition with each other: as long as they provide useful and unique information, do not interfere with each other and are not prohibitively cumbersome to use, they might find their proper place in future suites of sensors. In time, MRI acquisitions will likely include a plurality of complementary signals. A little like the microbiome that provides genetic diversity to organisms, these devices can provide signal diversity to MRI acquisitions and enrich measurements. Machine-learning (ML) algorithms are well suited at combining diverse input signals toward coherent outputs, and they could make use of all such information toward improved MRI capabilities. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.