Motion robust respiratory-resolved 3D radial flow MRI and its application in neonatal congenital heart disease

Advanced imaging of fetal cardiac function

Over recent decades, a variety of advanced imaging techniques for assessing cardiovascular physiology and cardiac function in adults and children have been applied in the fetus. In many cases, technical development has been required to allow feasibility in the fetus, while an appreciation of the unique physiology of the fetal circulation is required for proper interpretation of the findings. This review will focus on recent advances in fetal echocardiography and cardiovascular magnetic resonance (CMR), providing examples of their application in research and clinical settings. We will also consider future directions for these technologies, including their ongoing technical development and potential clinical value.

Gaussian Processes for real-time 3D motion and uncertainty estimation during MR-guided radiotherapy

Respiratory motion during radiotherapy causes uncertainty in the tumor's location, which is typically addressed by an increased radiation area and a decreased dose. As a result, the treatments' efficacy is reduced. The recently proposed hybrid MR-linac scanner holds the promise to efficiently deal with such respiratory motion through real-time adaptive MR-guided radiotherapy (MRgRT). For MRgRT, motion-fields should be estimated from MR-data and the radiotherapy plan should be adapted in real-time according to the estimated motion-fields. All of this should be performed with a total latency of maximally 200 ms, including data acquisition and reconstruction. A measure of confidence in such estimated motion-fields is highly desirable, for instance to ensure the patient's safety in case of unexpected and undesirable motion. In this work, we propose a framework based on Gaussian Processes to infer 3D motion-fields and uncertainty maps in real-time from only three readouts of MR-data. We demonstrated an inference frame rate up to 69 Hz including data acquisition and reconstruction, thereby exploiting the limited amount of required MR-data. Additionally, we designed a rejection criterion based on the motion-field uncertainty maps to demonstrate the framework's potential for quality assurance. The framework was validated in silico and in vivo on healthy volunteer data (n=5) acquired using an MR-linac, thereby taking into account different breathing patterns and controlled bulk motion. Results indicate end-point-errors with a 75th percentile below 1 mm in silico, and a correct detection of erroneous motion estimates with the rejection criterion. Altogether, the results show the potential of the framework for application in real-time MR-guided radiotherapy with an MR-linac.

Volumetric Fetal Flow Imaging With Magnetic Resonance Imaging

Fetal development relies on a complex circulatory network. Accurate assessment of flow distribution is important for understanding pathologies and potential therapies. In this paper, we demonstrate a method for volumetric imaging of fetal flow with magnetic resonance imaging (MRI). Fetal MRI faces challenges: small vascular structures, unpredictable motion, and inadequate traditional cardiac gating methods. Here, orthogonal multislice stacks are acquired with accelerated multidimensional radial phase contrast (PC) MRI. Slices are reconstructed into flow sensitive time-series images with motion correction and image-based cardiac gating. They are then combined into a dynamic volume using slice-to-volume reconstruction (SVR) while resolving interslice spatiotemporal coregistration. Compared to prior methods, this approach achieves higher spatiotemporal resolution ( 1×1×1 mm³, ~30 ms) with reduced scan time - important features for the quantification of flow through small fetal structures. Validation is demonstrated in adults by comparing SVR with 4D radial PCMRI (flow bias and limits of agreement: -1.1 ml/s and [-11.8 9.6] ml/s). Feasibility is demonstrated in late gestation fetuses by comparing SVR with 2D Cartesian PCMRI (flow bias and limits of agreement: -0.9 ml/min/kg and [-39.7 37.8] ml/min/kg). With SVR, we demonstrate complex flow pathways (such as parallel flow streams in the proximal inferior vena cava, preferential shunting of blood from the ductus venosus into the left atrium, and blood from the brain leaving the heart through the main pulmonary artery) for the first time in human fetal circulation. This method allows for comprehensive evaluation of the fetal circulation and enables future studies of fetal physiology.

Clinical Feasibility of Structural and Functional MRI in Free-Breathing Neonates and Infants

BACKGROUND

Evaluation of structural lung abnormalities with magnetic resonance imaging (MRI) has previously been shown to be predictive of clinical neonatal outcomes in preterm birth. MRI during free-breathing with phase-resolved functional lung (PREFUL) may allow for complimentary functional information without exogenous contrast.

PURPOSE

To investigate the feasibility of structural and functional pulmonary MRI in a cohort of neonates and infants with no cardiorespiratory disease. Macrovascular pulmonary blood flows were also evaluated.

STUDY TYPE

Prospective.

POPULATION

Ten term infants with no clinically defined cardiorespiratory disease were imaged. Infants recruited from the general population and neonatal intensive care unit (NICU) were studied.

FIELD STRENGTH/SEQUENCE

T₁ -weighted VIBE, T₂ -weighted BLADE uncorrected for motion. Ultrashort echo time (UTE) and 3D-flow data were acquired during free-breathing with self-navigation and retrospective reconstruction. Single slice 2D-gradient echo (GRE) images were acquired during free-breathing for PREFUL analysis. Imaging was performed at 3 T.

ASSESSMENT

T₁ , T₂ , and UTE images were scored according to the modified Ochiai scheme by three pediatric body radiologists. Ventilation/perfusion-weighted maps were extracted from free-breathing GRE images using PREFUL analysis. Ventilation and perfusion defect percent (VDP, QDP) were calculated from the segmented ventilation and perfusion-weighted maps. Time-averaged cardiac blood velocities from three-dimensional-flow were evaluated in major pulmonary arteries and veins.

STATISTICAL TEST

Intraclass correlation coefficient (ICC).

RESULTS

The ICC of replicate structural scores was 0.81 (95% CI: 0.45-0.95) across three observers. Elevated Ochiai scores, VDP, and QDP were observed in two NICU participants. Excluding these participants, mean ± standard deviation structural scores were 1.2 ± 0.8, while VDP and QDP were 1.0% ± 1.1% and 0.4% ± 0.5%, respectively. Main pulmonary arterial blood flows normalized to body surface area were 3.15 ± 0.78 L/min/m² .

DATA CONCLUSION

Structural and functional pulmonary imaging is feasible using standard clinical MRI hardware (commercial whole-body 3 T scanner, table spine array, and flexible thoracic array) in free-breathing infants.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Open or closed: Changes in ductus arteriosus flow patterns at birth using 4D flow MRI in newborn piglets

The ductus arteriosus (DA) functionally closes during the transition from fetal to postnatal life because of lung aeration and corresponding cardiovascular changes. The thorough and explicit measurement and visualization of the right and left heart output during this transition has not been previously accomplished. We combined 4D flow MRI (dynamic volumetric blood flow measurements) and T2 relaxometry (oxygen delivery quantification) in surgically instrumented newborn piglets to assess the DA. This was performed in Large White-Landrace-Duroc piglets (n = 34) spanning four age groups: term-9 days, term-3, term+1, and term+5. Subject's DA status was classified using 4D flow: closed DA, forward DA flow, reversed DA flow, and bidirectional DA flow. Over all states, vessel diameters and flows normalized to body weight increased with age (for example in the ascending aorta flow-term-9: 126.6 ± 45.4; term+5: 260.2 ± 80.0 ml/min per kg; p = 0.0005; ascending aorta diameter-term-9: 5.2 ± 0.8; term+5: 7.7 ± 0.4 mm; p = 0.0004). In subjects with reversed DA blood flow there was lower common carotid artery blood flow (forward: 37.5 ± 9.0; reversed: 20.0 ± 7.4 ml/min per kg; p = 0.032). Linear regression revealed that as net DA flow decreases, common carotid artery flow decreases (R2  = 0.32, p = 0.004), and left (R2  = 0.33, p = 0.003) and right (R2  = 0.34, p = 0.003) pulmonary artery flow increases. Bidirectional DA blood flow changed oxygen saturation as determined by MRI between the ascending and descending aorta (ascending aorta: 90.1% ± 8.4%; descending aorta: 75.6% ± 14.2%; p < 0.05). Expanded use of these techniques in preterm animal models will aid in providing new understandings of normal versus abnormal DA transition, as well as to test the effectiveness of related clinical interventions.

Using 5D flow MRI to decode the effects of rhythm on left atrial 3D flow dynamics in patients with atrial fibrillation

PURPOSE

This study used a 5D flow framework to explore the influence of arrhythmia on thrombogenic hemodynamic parameters in patients with atrial fibrillation (AF).

METHODS

A fully self-gated, 3D radial, highly accelerated free-running 5D flow sequence with interleaved four-point velocity-encoding was acquired using an in vitro arrhythmic flow phantom and in 25 patients with a history of AF (68 ± 8 y, 6 female). Self-gating signals were used to calculate AF burden, bin data, and tag each k-space line with its RRLength . Data were binned as an RR-resolved dataset with four RR-interval bins (RR1-RR4, short-to-long) for compressed sensing reconstruction. AF burden was calculated as interquartile range of all intrascan RR-intervals divided by median RR-interval, and left atrial (LA) stasis as the percent of the cardiac cycle where the velocity was <0.1 m/s.

RESULTS

In vitro results demonstrated successful recovery of RR-binned flow curves using RR-resolved 5D flow compared to a real-time PC reference standard. In vivo, 5D flow was acquired in 8:48 minutes. AF burden was significantly correlated with 5D flow-derived peak (PV) and mean (MV) velocity and stasis (|ρ| = 0.54-0.75, P < .001). Sensitivity analyses determined a threshold for low versus high AF burden at 9.7%. High burden patients had increased LA mean stasis (up to +42%, P < .01), and lower MV and PV (-30%, -40.6%, respectively, P < .01). RR4 deviated furthest from respiratory-resolved reconstruction (end-expiration) with increased mean stasis (7.6% ± 14.0%, P = .10) and decreased PV (-12.7 ± 14.2%, P = .09).

CONCLUSIONS

RR-resolved 5D flow can capture temporal and RR-resolved 3D hemodynamics in <10 minutes and offers a novel approach to investigate arrhythmias.

Additional value and new insights by four-dimensional flow magnetic resonance imaging in congenital heart disease: application in neonates and young children

Cardiovascular MRI has become an essential imaging modality in children with congenital heart disease (CHD) in the last 15-20 years. With use of appropriate sequences, it provides important information on cardiovascular anatomy, blood flow and function for initial diagnosis and post-surgical or -interventional monitoring in children. Although considered as more sophisticated and challenging than CT, in particular in neonates and infants, MRI is able to provide information on intra- and extracardiac haemodynamics, in contrast to CT. In recent years, four-dimensional (4-D) flow MRI has emerged as an additional MR technique for retrospective assessment and visualisation of blood flow within the heart and any vessel of interest within the acquired three-dimensional (3-D) volume. Its application in young children requires special adaptations for the smaller vessel size and faster heart rate compared to adolescents or adults. In this article, we provide an overview of 4-D flow MRI in various types of complex CHD in neonates and infants to demonstrate its potential indications and beneficial application for optimised individual cardiovascular assessment. We focus on its application in clinical routine cardiovascular workup and, in addition, show some examples with pathologies other than CHD to highlight that 4-D flow MRI yields new insights in disease understanding and therapy planning. We shortly review the essentials of 4-D flow data acquisition, pre- and post-processing techniques in neonates, infants and young children. Finally, we conclude with some details on accuracy, limitations and pitfalls of the technique.

5D Flow MRI: A Fully Self-gated, Free-running Framework for Cardiac and Respiratory Motion-resolved 3D Hemodynamics

PURPOSE

To implement, validate, and apply a self-gated free-running whole-heart five-dimensional (5D) flow MRI framework to evaluate respiration-driven effects on three-dimensional (3D) hemodynamics in a clinical setting.

MATERIALS AND METHODS

In this prospective study, a free-running five-dimensional (5D) flow sequence was implemented with 3D radial sampling, self-gating, and a compressed-sensing reconstruction. The 5D flow was evaluated in a pulsatile phantom and adult participants with aortic and/or valvular disease who were enrolled between May and August 2019. Conventional twofold-accelerated four-dimensional (4D) flow of the thoracic aorta with navigator gating was performed as a reference comparison. Continuous parameters were evaluated for parameter normality and were compared between conventional 4D flow and 5D flow using a signed-rank or two-tailed paired t test. Differences between respiratory states were evaluated using a repeated-measure analysis of variance or a nonparametric Friedman test.

RESULTS

A total of 20 adult participants (mean age, 49 years ± 17 [standard deviation]; 18 men and two women) were included. In vitro 5D flow results showed excellent agreement with conventional 4D flow-derived values (peak and net flow, <7% difference over all quantified planes). Whole-heart 5D flow data were collected in all participants in 7.65 minutes ± 0.35 (acceleration rate = 36.0-76.9) versus 9.88 minutes ± 3.17 for conventional aortic 4D flow. In vivo, 5D flow demonstrated moderate agreement with conventional 4D flow but demonstrated overestimation in net flow and peak velocity (up to 26% and 12%, respectively) in the ascending aorta and underestimation (<12%) in the arch and descending aorta. Respiratory-resolved analyses of caval veins showed significantly increased net and peak flow in the inferior vena cava in end inspiration compared with end expiration, and the opposite trend was shown in the superior vena cava.

CONCLUSION

A free-running 5D flow MRI framework consistently captured cardiac and respiratory motion-resolved 3D hemodynamics in less than 8 minutes. Supplemental material is available for this article. © RSNA, 2020.