Placental MRI: Development of an MRI compatible ex vivo system for whole placenta dual perfusion

PURPOSE

Placental dysfunction plays a key role in diseases that affect the fetus in utero and after birth. Aiming to develop a platform for validating in vivo placental MRI and investigations into placental physiology, we designed and built a prototype MRI-compatible perfusion chamber with an integrated MRI receive coil for high SNR ex vivo placental imaging.

PRINCIPAL RESULTS

After optimizing placenta vascular clearing and perfusion protocols, we performed contrast enhanced MR angiography and MR relaxometry on eight carefully selected placentas while they were perfused via the umbilical arteries (UAs). Additionally, two of these placentas underwent maternal perfusion via the intervillous space (IVS). Despite striving for homogenous perfusion across the whole placenta, imaging results were highly heterogeneous for both UA and IVS perfused placentas. By histology, we observed blood congestion in the villi in regions that showed low UA perfusion during MRI. In two placentas prominent chorionic arteries followed by adjacent veins underwent contrast enhancement in the absence of villous capillary blush. The single placenta from a pregnancy affected by IUGR had the most homogeneous villous capillary perfusion.

MAJOR CONCLUSIONS

A dual perfusion system for ex vivo placentas compatible with MRI permitted assessment of UA and IVS placental perfusion. We observed spatial UA perfusion heterogeneity and evidence for arteriovenous shunting in placentas from normal pregnancies and deliveries, but relative villous capillary perfusion homogeneity in a single IUGR placenta. Future work will focus on system optimization, followed by physiological manipulation and validation of in vivo placental MRI.

Individualized SAR calculations using computer vision-based MR segmentation and a fast electromagnetic solver

PURPOSE

We propose a fast, patient-specific workflow for on-line specific absorption rate (SAR) supervision. An individualized electromagnetic model is created while the subject is on the table, followed by rapid SAR estimates for that individual. Our goal is an improved correspondence between the patient and model, reducing reliance on general anatomical body models.

METHODS

A 3D fat-water 3T acquisition (~2 minutes) is automatically segmented using a computer vision algorithm (~1 minute) into what we found to be the most important electromagnetic tissue classes: air, bone, fat, and soft tissues. We then compute the individual's EM field exposure and global and local SAR matrices using a fast electromagnetic integral equation solver. We assess the approach in 10 volunteers and compare to the SAR seen in a standard generic body model (Duke).

RESULTS

The on-the-table workflow averaged 7'44″. Simulation of the simplified Duke models confirmed that only air, bone, fat, and soft tissue classes are needed to estimate global and local SAR with an error of 6.7% and 2.7%, respectively, compared to the full model. In contrast, our volunteers showed a 16.0% and 20.3% population variability in global and local SAR, respectively, which was mostly underestimated by the Duke model.

CONCLUSION

Timely construction and deployment of a patient-specific model is computationally feasible. The benefit of resolving the population heterogeneity compared favorably to the modest modeling error incurred. This suggests that individualized SAR estimates can improve electromagnetic safety in MRI and possibly reduce conservative safety margins that account for patient-model mismatch, especially in non-standard patients.

Placental MRI: Effect of maternal position and uterine contractions on placental BOLD MRI measurements

INTRODUCTION

Before using blood-oxygen-level-dependent magnetic resonance imaging (BOLD MRI) during maternal hyperoxia as a method to detect individual placental dysfunction, it is necessary to understand spatiotemporal variations that represent normal placental function. We investigated the effect of maternal position and Braxton-Hicks contractions on estimates obtained from BOLD MRI of the placenta during maternal hyperoxia.

METHODS

For 24 uncomplicated singleton pregnancies (gestational age 27-36 weeks), two separate BOLD MRI datasets were acquired, one in the supine and one in the left lateral maternal position. The maternal oxygenation was adjusted as 5 min of room air (21% O₂), followed by 5 min of 100% FiO₂. After datasets were corrected for signal non-uniformities and motion, global and regional BOLD signal changes in R₂* and voxel-wise Time-To-Plateau (TTP) in the placenta were measured. The overall placental and uterine volume changes were determined across time to detect contractions.

RESULTS

In mothers without contractions, increases in global placental R₂* in the supine position were larger compared to the left lateral position with maternal hyperoxia. Maternal position did not alter global TTP but did result in regional changes in TTP. 57% of the subjects had Braxton-Hicks contractions and 58% of these had global placental R₂* decreases during the contraction.

CONCLUSION

Both maternal position and Braxton-Hicks contractions significantly affect global and regional changes in placental R₂* and regional TTP. This suggests that both factors must be taken into account in analyses when comparing placental BOLD signals over time within and between individuals.

Design and implementation of a low-cost, tabletop MRI scanner for education and research prototyping

While access to a laboratory MRI system is ideal for teaching MR physics as well as many aspects of signal processing, providing multiple MRI scanners can be prohibitively expensive for educational settings. To address this need, we developed a small, low-cost, open-interface tabletop MRI scanner for academic use. We constructed and tested 20 of these scanners for parallel use by teams of 2-3 students in a teaching laboratory. With simplification and down-scaling to a 1 cm FOV, fully-functional scanners were achieved within a budget of $10,000 USD each. The design was successful for teaching MR principles and basic signal processing skills and serves as an accessible testbed for more advanced MR research projects. Customizable GUIs, pulse sequences, and reconstruction code accessible to the students facilitated tailoring the scanner to the needs of laboratory exercise. The scanners have been used by >800 students in 6 different courses and all designs, schematics, sequences, GUIs, and reconstruction code is open-source.

Comparison of CBF Measured with Combined Velocity-Selective Arterial Spin-Labeling and Pulsed Arterial Spin-Labeling to Blood Flow Patterns Assessed by Conventional Angiography in Pediatric Moyamoya

BACKGROUND AND PURPOSE

Imaging CBF is important for managing pediatric moyamoya. Traditional arterial spin-labeling MR imaging detects delayed transit thorough diseased arteries but is inaccurate for measuring perfusion because of these delays. Velocity-selective arterial spin-labeling is insensitive to transit delay and well-suited for imaging Moyamoya perfusion. This study assesses the accuracy of a combined velocity-selective arterial spin-labeling and traditional pulsed arterial spin-labeling CBF approach in pediatric moyamoya, with comparison to blood flow patterns on conventional angiography.

MATERIALS AND METHODS

Twenty-two neurologically stable pediatric patients with moyamoya and 5 asymptomatic siblings without frank moyamoya were imaged with velocity-selective arterial spin-labeling, pulsed arterial spin-labeling, and DSA (patients). Qualitative comparison was performed, followed by a systematic comparison using ASPECTS-based scoring. Quantitative pulsed arterial spin-labeling CBF and velocity-selective arterial spin-labeling CBF for the middle cerebral artery, anterior cerebral artery, and posterior cerebral artery territories were also compared.

RESULTS

Qualitatively, velocity-selective arterial spin-labeling perfusion maps reflect the DSA parenchymal phase, regardless of postinjection timing. Conversely, pulsed arterial spin-labeling maps reflect the DSA appearance at postinjection times closer to the arterial spin-labeling postlabeling delay, regardless of vascular phase. ASPECTS comparison showed excellent agreement (88%, κ = 0.77, P < .001) between arterial spin-labeling and DSA, suggesting velocity-selective arterial spin-labeling and pulsed arterial spin-labeling capture key perfusion and transit delay information, respectively. CBF coefficient of variation, a marker of perfusion variability, was similar for velocity-selective arterial spin-labeling in patient regions of delayed-but-preserved perfusion compared to healthy asymptomatic sibling regions (coefficient of variation = 0.30 versus 0.26, respectively, Δcoefficient of variation = 0.04), but it was significantly different for pulsed arterial spin-labeling (coefficient of variation = 0.64 versus 0.34, Δcoefficient of variation = 0.30, P < .001).

CONCLUSIONS

Velocity-selective arterial spin-labeling offers a powerful approach to image perfusion in pediatric moyamoya due to transit delay insensitivity. Coupled with pulsed arterial spin-labeling for transit delay information, a volumetric MR imaging approach capturing key DSA information is introduced.

Individual variation in simulated fetal SAR assessed in multiple body models

PURPOSE

We generate 12 models from 4 pregnant individuals to evaluate individual differences in local specific absorption rate (SAR) for differing body habitus and fetal and maternal positions.

METHODS

Structural MR images from 4 pregnant subjects (including supine and left-lateral maternal positions) were manually segmented to create 12 body models by rotating the fetus, modifying the fat content, and altering the maternal arm position in 1 of the subjects. Electromagnetic simulations modeled at 3 Tesla determined the average and peak local SAR in the maternal trunk, fetus, fetal brain, and amniotic fluid.

RESULTS

We observed a significant range of fetal and maternal peak local SAR across the models (maternal trunk: 19.14-44.03 watts/kg, fetus: 9.93-18.79 watts/kg, fetal brain 3.36-10.3 watts/kg). We found that maternal body habitus changes introduced a significant variation in the maternal peak local SAR but not the fetal local SAR. However, the maternal position (either rotating the mother to left-lateral position or altering the arm position) introduced changes in fetal peak local SAR (range: 11.9-17.9 watts/kg). Rotating the fetus also introduced variation in the fetal and fetal brain peak local SAR.

CONCLUSION

The observed variation in SAR emphasizes the need for more anatomical models to enable better safety management of individuals during fetal MRI, including a wider range of gestational ages.

Placental MRI: Developing Accurate Quantitative Measures of Oxygenation

The Human Placenta Project has focused attention on the need for noninvasive magnetic resonance imaging (MRI)-based techniques to diagnose and monitor placental function throughout pregnancy. The hope is that the management of placenta-related pathologies would be improved if physicians had more direct, real-time measures of placental health to guide clinical decision making. As oxygen alters signal intensity on MRI and oxygen transport is a key function of the placenta, many of the MRI methods under development are focused on quantifying oxygen transport or oxygen content of the placenta. For example, measurements from blood oxygen level-dependent imaging of the placenta during maternal hyperoxia correspond to outcomes in twin pregnancies, suggesting that some aspects of placental oxygen transport can be monitored by MRI. Additional methods are being developed to accurately quantify baseline placental oxygenation by MRI relaxometry. However, direct validation of placental MRI methods is challenging and therefore animal studies and ex vivo studies of human placentas are needed. Here we provide an overview of the current state of the art of oxygen transport and quantification with MRI. We suggest that as these techniques are being developed, increased focus be placed on ensuring they are robust and reliable across individuals and standardized to enable predictive diagnostic models to be generated from the data. The field is still several years away from establishing the clinical benefit of monitoring placental function in real time with MRI, but the promise of individual personalized diagnosis and monitoring of placental disease in real time continues to motivate this effort.

Computer-Vision Techniques for Water-Fat Separation in Ultra High-Field MRI Local Specific Absorption Rate Estimation

OBJECTIVE

The purpose of this paper is to prove that computer-vision techniques allow synthesizing water-fat separation maps for local specific absorption rate (SAR) estimation, when patient-specific water-fat images are not available.

METHODS

We obtained ground truth head models by using patient-specific water-fat images. We obtained two different label-fusion water-fat models generating a water-fat multiatlas and applying the STAPLE and local-MAP-STAPLE label-fusion methods. We also obtained patch-based water-fat models applying a local group-wise weighted combination of the multiatlas. Electromagnetic (EM) simulations were performed, and B1+ magnitude and 10 g averaged SAR maps were generated.

RESULTS

We found local approaches provide a high DICE overlap (72.6 ± 10.2% fat and 91.6 ± 1.5% water in local-MAP-STAPLE, and 68.8 ± 8.2% fat and 91.1 ± 1.0% water in patch-based), low Hausdorff distances (18.6 ± 7.7 mm fat and 7.4 ± 11.2 mm water in local-MAP-STAPLE, and 16.4 ± 8.5 mm fat and 7.2 ± 11.8 mm water in patch-based) and a low error in volume estimation (15.6 ± 34.4% fat and 5.6 ± 4.1% water in the local-MAP-STAPLE, and 14.0 ± 17.7% fat and 4.7 ± 2.8% water in patch-based). The positions of the peak 10 g-averaged local SAR hotspots were the same for every model.

CONCLUSION

We have created patient-specific head models using three different computer-vision-based water-fat separation approaches and compared the predictions of B1+ field and SAR distributions generated by simulating these models. Our results prove that a computer-vision approach can be used for patient-specific water-fat separation, and utilized for local SAR estimation in high-field MRI.

SIGNIFICANCE

Computer-vision approaches can be used for patient-specific water-fat separation and for patient specific local SAR estimation, when water-fat images of the patient are not available.

Simultaneous multislice magnetic resonance fingerprinting with low-rank and subspace modeling

Magnetic resonance fingerprinting (MRF) is a new quantitative imaging paradigm that enables simultaneous acquisition of multiple magnetic resonance tissue parameters (e.g., T₁, T₂, and spin density). Recently, MRF has been integrated with simultaneous multislice (SMS) acquisitions to enable volumetric imaging with faster scan time. In this paper, we present a new image reconstruction method based on low-rank and subspace modeling for improved SMS-MRF. Here the low-rank model exploits strong spatiotemporal correlation among contrast-weighted images, while the subspace model captures the temporal evolution of magnetization dynamics. With the proposed model, the image reconstruction problem is formulated as a convex optimization problem, for which we develop an algorithm based on variable splitting and the alternating direction method of multipliers. The performance of the proposed method has been evaluated by numerical experiments, and the results demonstrate that the proposed method leads to improved accuracy over the conventional approach. Practically, the proposed method has a potential to allow for a 3× speedup with minimal reconstruction error, resulting in less than 5 sec imaging time per slice.

Motion-robust sub-millimeter isotropic diffusion imaging through motion corrected generalized slice dithered enhanced resolution (MC-gSlider) acquisition

PURPOSE

To develop an efficient MR technique for ultra-high resolution diffusion MRI (dMRI) in the presence of motion.

METHODS

gSlider is an SNR-efficient high-resolution dMRI acquisition technique. However, subject motion is inevitable during a prolonged scan for high spatial resolution, leading to potential image artifacts and blurring. In this study, an integrated technique termed Motion Corrected gSlider (MC-gSlider) is proposed to obtain high-quality, high-resolution dMRI in the presence of large in-plane and through-plane motion. A motion-aware reconstruction with spatially adaptive regularization is developed to optimize the conditioning of the image reconstruction under difficult through-plane motion cases. In addition, an approach for intra-volume motion estimation and correction is proposed to achieve motion correction at high temporal resolution.

RESULTS

Theoretical SNR and resolution analysis validated the efficiency of MC-gSlider with regularization, and aided in selection of reconstruction parameters. Simulations and in vivo experiments further demonstrated the ability of MC-gSlider to mitigate motion artifacts and recover detailed brain structures for dMRI at 860 μm isotropic resolution in the presence of motion with various ranges.

CONCLUSION

MC-gSlider provides motion-robust, high-resolution dMRI with a temporal motion correction sensitivity of 2 s, allowing for the recovery of fine detailed brain structures in the presence of large subject movements.

Computation of ultimate SAR amplification factors for radiofrequency hyperthermia in non-uniform body models: impact of frequency and tumour location

PURPOSE

We introduce a method for calculation of the ultimate specific absorption rate (SAR) amplification factors (uSAF) in non-uniform body models. The uSAF is the greatest possible SAF achievable by any hyperthermia (HT) phased array for a given frequency, body model and target heating volume.

METHODS

First, we generate a basis-set of solutions to Maxwell's equations inside the body model. We place a large number of electric and magnetic dipoles around the body model and excite them with random amplitudes and phases. We then compute the electric fields created in the body model by these excitations using an ultra-fast volume integral solver called MARIE. We express the field pattern that maximises the SAF in the target tumour as a linear combination of these basis fields and optimise the combination weights so as to maximise SAF (concave problem). We compute the uSAFs in the Duke body models at 10 frequencies in the 20-900 MHz range and for twelve 3 cm-diameter tumours located at various depths in the head and neck.

RESULTS

For both shallow and deep tumours, the frequency yielding the greatest uSAF was ∼900 MHz. Since this is the greatest frequency that we simulated, we hypothesise that the globally optimal frequency is actually greater.

CONCLUSIONS

The uSAFs computed in this work are very large (40-100 for shallow tumours and 4-17 for deep tumours), indicating that there is a large room for improvement of the current state-of-the-art head and neck HT devices.

Assessing the effects of subject motion on T2 relaxation under spin tagging (TRUST) cerebral oxygenation measurements using volume navigators

PURPOSE

Subject motion may cause errors in estimates of blood T2 when using the T2 -relaxation under spin tagging (TRUST) technique on noncompliant subjects like neonates. By incorporating 3D volume navigators (vNavs) into the TRUST pulse sequence, independent measurements of motion during scanning permit evaluation of these errors.

METHODS

The effects of integrated vNavs on TRUST-based T2 estimates were evaluated using simulations and in vivo subject data. Two subjects were scanned with the TRUST+vNav sequence during prescribed movements. Mean motion scores were derived from vNavs and TRUST images, along with a metric of exponential fit quality. Regression analysis was performed between T2 estimates and mean motion scores. Also, motion scores were determined from independent neonatal scans.

RESULTS

vNavs negligibly affected venous blood T2 estimates and better detected subject motion than fit quality metrics. Regression analysis showed that T2 is biased upward by 4.1 ms per 1 mm of mean motion score. During neonatal scans, mean motion scores of 0.6 to 2.0 mm were detected.

CONCLUSION

Motion during TRUST causes an overestimate of T2 , which suggests a cautious approach when comparing TRUST-based cerebral oxygenation measurements of noncompliant subjects. Magn Reson Med 78:2283-2289, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Functional oxygen extraction fraction (OEF) imaging with turbo gradient spin echo QUIXOTIC (Turbo QUIXOTIC)

PURPOSE

QUantitative Imaging of eXtraction of Oxygen and TIssue Consumption (QUIXOTIC) is a recent technique that measures voxel-wise oxygen extraction fraction (OEF) but suffers from long scan times, limiting its application. We implemented multiecho QUIXOTIC dubbed turbo QUIXOTIC (tQUIXOTIC) that reduces scan time eightfold and then applied it in functional MRI.

METHODS

tQUIXOTIC utilizes a novel turbo gradient spin echo readout enabling measurement of venular blood transverse relaxation rate in a single tag-control acquisition. Using tQUIXOTIC, we estimated cortical gray matter (GM) OEF, created voxel-by-voxel GM OEF maps, and quantified changes in visual cortex OEF during a blocked design flashing checkerboard visual stimulus. Contamination from cerebrospinal fluid partial volume averaging was estimated and corrected.

RESULTS

The average cortical GM OEF was estimated as 0.38 ± 0.06 (n = 8) using a 3.4-min acquisition. The average OEF in the visual cortex was estimated as 0.43 ± 0.04 at baseline and 0.35 ± 0.05 during activation, with an average %ΔOEF of -20%. These values are consistent with those of past studies.

CONCLUSION

tQUIXOTIC successfully estimated cortical GM OEF in clinical scan times and detected changes in OEF during blocked design visual stimulation. tQUIXOTIC will be useful to monitor regional OEF clinically and in blocked design or event-related functional MRI experiments. Magn Reson Med 79:2713-2723, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

In Vivo Quantification of Placental Insufficiency by BOLD MRI: A Human Study

Fetal health is critically dependent on placental function, especially placental transport of oxygen from mother to fetus. When fetal growth is compromised, placental insufficiency must be distinguished from modest genetic growth potential. If placental insufficiency is present, the physician must trade off the risk of prolonged fetal exposure to placental insufficiency against the risks of preterm delivery. Current ultrasound methods to evaluate the placenta are indirect and insensitive. We propose to use Blood-Oxygenation-Level-Dependent (BOLD) MRI with maternal hyperoxia to quantitatively assess mismatch in placental function in seven monozygotic twin pairs naturally matched for genetic growth potential. In-utero BOLD MRI time series were acquired at 29 to 34 weeks gestational age. Maps of oxygen Time-To-Plateau (TTP) were obtained in the placentas by voxel-wise fitting of the time series. Fetal brain and liver volumes were measured based on structural MR images. After delivery, birth weights were obtained and placental pathological evaluations were performed. Mean placental TTP negatively correlated with fetal liver and brain volumes at the time of MRI as well as with birth weights. Mean placental TTP positively correlated with placental pathology. This study demonstrates the potential of BOLD MRI with maternal hyperoxia to quantify regional placental function in vivo.

Improved magnetic resonance fingerprinting reconstruction with low-rank and subspace modeling

PURPOSE

This article introduces a constrained imaging method based on low-rank and subspace modeling to improve the accuracy and speed of MR fingerprinting (MRF).

THEORY AND METHODS

A new model-based imaging method is developed for MRF to reconstruct high-quality time-series images and accurate tissue parameter maps (e.g., T₁ , T₂ , and spin density maps). Specifically, the proposed method exploits low-rank approximations of MRF time-series images, and further enforces temporal subspace constraints to capture magnetization dynamics. This allows the time-series image reconstruction problem to be formulated as a simple linear least-squares problem, which enables efficient computation. After image reconstruction, tissue parameter maps are estimated via dictionary-based pattern matching, as in the conventional approach.

RESULTS

The effectiveness of the proposed method was evaluated with in vivo experiments. Compared with the conventional MRF reconstruction, the proposed method reconstructs time-series images with significantly reduced aliasing artifacts and noise contamination. Although the conventional approach exhibits some robustness to these corruptions, the improved time-series image reconstruction in turn provides more accurate tissue parameter maps. The improvement is pronounced especially when the acquisition time becomes short.

CONCLUSIONS

The proposed method significantly improves the accuracy of MRF, and also reduces data acquisition time. Magn Reson Med 79:933-942, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Mitigation of partial volume effects in susceptibility-based oxygenation measurements by joint utilization of magnitude and phase (JUMP)

PURPOSE

Susceptibility-based blood oxygenation measurements in small vessels of the brain derive from gradient echo (GRE) phase and can provide localized assessment of brain function and pathology. However, when vessel diameter becomes smaller than the acquisition voxel size, partial volume effects compromise these measurements. The purpose of this study was to develop a technique to improve the reliability of vessel oxygenation estimates in the presence of partial volume effects.

METHODS

Intravoxel susceptibility variations are present when a vessel and parenchyma experience partial volume effects, modifying the voxel's GRE phase signal and attenuating the GRE magnitude signal. Using joint utilization of magnitude and phase (JUMP), both vessel susceptibility and voxel partial volume fraction can be estimated, providing measurements of venous oxygen saturation ( Yv) in straight, nearly vertical vessels that have improved robustness to partial volume effects.

RESULTS

JUMP was demonstrated by estimating vessel Yv in numerical and in vivo experiments. Deviations from ground truth of Yv measurements in vessels tilted up to 30° from B0 were reduced by over 50% when using JUMP compared with phase-only techniques.

CONCLUSION

JUMP exploits both magnitude and phase data in GRE imaging to mitigate partial volume effects in estimation of vessel oxygenation. Magn Reson Med 77:1713-1727, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Single-step quantitative susceptibility mapping with variational penalties

Quantitative susceptibility mapping (QSM) estimates the underlying tissue magnetic susceptibility from the gradient echo (GRE) phase signal through background phase removal and dipole inversion steps. Each of these steps typically requires the solution of an ill-posed inverse problem and thus necessitates additional regularization. Recently developed single-step QSM algorithms directly relate the unprocessed GRE phase to the unknown susceptibility distribution, thereby requiring the solution of a single inverse problem. In this work, we show that such a holistic approach provides susceptibility estimation with artifact mitigation and develop efficient algorithms that involve simple analytical solutions for all of the optimization steps. Our methods employ total variation (TV) and total generalized variation (TGV) to jointly perform the background removal and dipole inversion in a single step. Using multiple spherical mean value (SMV) kernels of varying radii permits high-fidelity background removal whilst retaining the phase information in the cortex. Using numerical simulations, we demonstrate that the proposed single-step methods reduce the reconstruction error by up to 66% relative to the multi-step methods that involve SMV background filtering with the same number of SMV kernels, followed by TV- or TGV-regularized dipole inversion. In vivo single-step experiments demonstrate a dramatic reduction in dipole streaking artifacts and improved homogeneity of image contrast. These acquisitions employ the rapid three-dimensional echo planar imaging (3D EPI) and Wave-CAIPI (controlled aliasing in parallel imaging) trajectories for signal-to-noise ratio-efficient whole-brain imaging. Herein, we also demonstrate the multi-echo capability of the Wave-CAIPI sequence for the first time, and introduce an automated, phase-sensitive coil sensitivity estimation scheme based on a 4-s calibration acquisition. Copyright © 2016 John Wiley & Sons, Ltd.

Spatiotemporal alignment of in utero BOLD-MRI series

PURPOSE

To present a method for spatiotemporal alignment of in-utero magnetic resonance imaging (MRI) time series acquired during maternal hyperoxia for enabling improved quantitative tracking of blood oxygen level-dependent (BOLD) signal changes that characterize oxygen transport through the placenta to fetal organs.

MATERIALS AND METHODS

The proposed pipeline for spatiotemporal alignment of images acquired with a single-shot gradient echo echo-planar imaging includes 1) signal nonuniformity correction, 2) intravolume motion correction based on nonrigid registration, 3) correction of motion and nonrigid deformations across volumes, and 4) detection of the outlier volumes to be discarded from subsequent analysis. BOLD MRI time series collected from 10 pregnant women during 3T scans were analyzed using this pipeline. To assess pipeline performance, signal fluctuations between consecutive timepoints were examined. In addition, volume overlap and distance between manual region of interest (ROI) delineations in a subset of frames and the delineations obtained through propagation of the ROIs from the reference frame were used to quantify alignment accuracy. A previously demonstrated rigid registration approach was used for comparison.

RESULTS

The proposed pipeline improved anatomical alignment of placenta and fetal organs over the state-of-the-art rigid motion correction methods. In particular, unexpected temporal signal fluctuations during the first normoxia period were significantly decreased (P < 0.01) and volume overlap and distance between region boundaries measures were significantly improved (P < 0.01).

CONCLUSION

The proposed approach to align MRI time series enables more accurate quantitative studies of placental function by improving spatiotemporal alignment across placenta and fetal organs.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:403-412.

Construction and modeling of a reconfigurable MRI coil for lowering SAR in patients with deep brain stimulation implants

Post-operative MRI of patients with deep brain simulation (DBS) implants is useful to assess complications and diagnose comorbidities, however more than one third of medical centers do not perform MRIs on this patient population due to stringent safety restrictions and liability risks. A new system of reconfigurable magnetic resonance imaging head coil composed of a rotatable linearly-polarized birdcage transmitter and a close-fitting 32-channel receive array is presented for low-SAR imaging of patients with DBS implants. The novel system works by generating a region with low electric field magnitude and steering it to coincide with the DBS lead trajectory. We demonstrate that the new coil system substantially reduces the SAR amplification around DBS electrodes compared to commercially available circularly polarized coils in a cohort of 9 patient-derived realistic DBS lead trajectories. We also show that the optimal coil configuration can be reliably identified from the image artifact on B₁⁺ field maps. Our preliminary results suggest that such a system may provide a viable solution for high-resolution imaging of DBS patients in the future. More data is needed to quantify safety limits and recommend imaging protocols before the novel coil system can be used on patients with DBS implants.

The ultimate signal-to-noise ratio in realistic body models

PURPOSE

We compute the ultimate signal-to-noise ratio (uSNR) and G-factor (uGF) in a realistic head model from 0.5 to 21 Tesla.

METHODS

We excite the head model and a uniform sphere with a large number of electric and magnetic dipoles placed at 3 cm from the object. The resulting electromagnetic fields are computed using an ultrafast volume integral solver, which are used as basis functions for the uSNR and uGF computations.

RESULTS

Our generalized uSNR calculation shows good convergence in the sphere and the head and is in close agreement with the dyadic Green's function approach in the uniform sphere. In both models, the uSNR versus B₀ trend was linear at shallow depths and supralinear at deeper locations. At equivalent positions, the rate of increase of the uSNR with B₀ was greater in the sphere than in the head model. The uGFs were lower in the realistic head than in the sphere for acceleration in the anterior-posterior direction, but similar for the left-right direction.

CONCLUSION

The uSNR and uGFs are computable in nonuniform body models and provide fundamental performance limits for human imaging with close-fitting MRI array coils. Magn Reson Med 78:1969-1980, 2017. © 2016 International Society for Magnetic Resonance in Medicine.