4D flow MRI velocity uncertainty quantification.

PURPOSE: An automatic method is presented for estimating 4D flow MRI velocity measurement uncertainty in each voxel. The velocity distance (VD) metric, a statistical distance between the measured velocity and local error distribution, is introduced as a novel measure of 4D flow MRI velocity measurement quality. METHODS: The method uses mass conservation to assess the local velocity error variance and the standardized difference of means (SDM) velocity to estimate the velocity error correlations. VD is evaluated as the Mahalanobis distance between the local velocity measurement and the local error distribution. The uncertainty model is validated synthetically and tested in vitro under different flow resolutions and noise levels. The VD's application is demonstrated on two in vivo thoracic vasculature 4D flow datasets. RESULTS: Synthetic results show the proposed uncertainty quantification method is sensitive to aliased regions across various velocity-to-noise ratios and assesses velocity error correlations in four- and six-point acquisitions with correlation errors at or under 3.2%. In vitro results demonstrate the method's sensitivity to spatial resolution, venc settings, partial volume effects, and phase wrapping error sources. Applying VD to assess in vivo 4D flow MRI in the aorta demonstrates the expected increase in measured velocity quality with contrast administration and systolic flow. CONCLUSION: The proposed 4D flow MRI uncertainty quantification method assesses velocity measurement error owing to sources including noise, intravoxel phase dispersion, and velocity aliasing. This method enables rigorous comparison of 4D flow MRI datasets obtained in longitudinal studies, across patient populations, and with different MRI systems.

GRASP reconstruction amplified with view-sharing and KWIC filtering reduces underestimation of peak velocity in highly-accelerated real-time phase-contrast MRI: A preliminary evaluation in pediatric patients with congenital heart disease.

PURPOSE: To develop a highly-accelerated, real-time phase contrast (rtPC) MRI pulse sequence with 40 fps frame rate (25 ms effective temporal resolution). METHODS: Highly-accelerated golden-angle radial sparse parallel (GRASP) with over regularization may result in temporal blurring, which in turn causes underestimation of peak velocity. Thus, we amplified GRASP performance by synergistically combining view-sharing (VS) and k-space weighted image contrast (KWIC) filtering. In 17 pediatric patients with congenital heart disease (CHD), the conventional GRASP and the proposed GRASP amplified by VS and KWIC (or GRASP + VS + KWIC) reconstruction for rtPC MRI were compared with respect to clinical standard PC MRI in measuring hemodynamic parameters (peak velocity, forward volume, backward volume, regurgitant fraction) at four locations (aortic valve, pulmonary valve, left and right pulmonary arteries). RESULTS: The proposed reconstruction method (GRASP + VS + KWIC) achieved better effective spatial resolution (i.e., image sharpness) compared with conventional GRASP, ultimately reducing the underestimation of peak velocity from 17.4% to 6.4%. The hemodynamic metrics (peak velocity, volumes) were not significantly (p > 0.99) different between GRASP + VS + KWIC and clinical PC, whereas peak velocity was significantly (p < 0.007) lower for conventional GRASP. RtPC with GRASP + VS + KWIC also showed the ability to assess beat-to-beat variation and detect the highest peak among peaks. CONCLUSION: The synergistic combination of GRASP, VS, and KWIC achieves 25 ms effective temporal resolution (40 fps frame rate), while minimizing the underestimation of peak velocity compared with conventional GRASP.

Multiyear Interval Changes in Aortic Wall Shear Stress in Patients with Bicuspid Aortic Valve Assessed by 4D Flow MRI.

BACKGROUND: In patients with bicuspid aortic valve (BAV), 4D flow MRI can quantify regions exposed to abnormal aortic hemodynamics, including high wall shear stress (WSS), a known stimulus for arterial wall dysfunction. However, the long-term multiscan reproducibility of 4D flow MRI-derived hemodynamic parameters is unknown. PURPOSE: To investigate the long-term stability of 4D flow MRI-derived peak velocity, WSS, and WSS-derived heatmaps in patients with BAV undergoing multiyear surveillance imaging. STUDY TYPE: Retrospective. POPULATION: 20 BAV patients (mean age 48.4 ± 13.9 years; 14 males) with five 4D flow MRI scans, with intervals of at least 6 months between scans, and 125 controls (mean age: 50.7 ± 15.8 years; 67 males). FIELD STRENGTH/SEQUENCE: 1.5 and 3.0T, prospectively ECG and respiratory navigator-gated aortic 4D flow MRI. ASSESSMENT: Automated AI-based 4D flow analysis pipelines were used for data preprocessing, aorta 3D segmentation, and quantification of ascending aorta (AAo) peak velocity, peak systolic WSS, and heatmap-derived relative area of elevated WSS compared to WSS ranges in age and sex-matched normative control populations. Growth rate was derived from the maximum AAo diameters measured on the first and fifth MRI scans. STATISTICAL TESTS: One-way repeated measures analysis of variance. P < 0.05 indicated significance. RESULTS: One hundred 4D flow MRI exams (five per patient) were analyzed. The mean total follow-up duration was 5.5 ± 1.1 years, and the average growth rate was 0.3 ± 0.2 mm/year. Peak velocity, peak systolic WSS, and relative area of elevated WSS did not change significantly over the follow-up period (P = 0.64, P = 0.69, and P = 0.35, respectively). The patterns and areas of elevated WSS demonstrated good reproducibility on semiquantitative assessment. CONCLUSION: 4D flow MRI-derived peak velocity, WSS, and WSS-derived heatmaps showed good multiyear and multiscan stability in BAV patients with low aortic growth rates. These findings underscore the reliability of these metrics in monitoring BAV patients for potential risk of dilation. LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY: Stage 1.

Environmental Sustainability and MRI: Challenges, Opportunities, and a Call for Action.

The environmental impact of magnetic resonance imaging (MRI) has recently come into focus. This includes its enormous demand for electricity compared to other imaging modalities and contamination of water bodies with anthropogenic gadolinium related to contrast administration. Given the pressing threat of climate change, addressing these challenges to improve the environmental sustainability of MRI is imperative. The purpose of this review is to discuss the challenges, opportunities, and the need for action to reduce the environmental impact of MRI and prepare for the effects of climate change. The approaches outlined are categorized as strategies to reduce greenhouse gas (GHG) emissions from MRI during production and use phases, approaches to reduce the environmental impact of MRI including the preservation of finite resources, and development of adaption plans to prepare for the impact of climate change. Co-benefits of these strategies are emphasized including lower GHG emission and reduced cost along with improved heath and patient satisfaction. Although MRI is energy-intensive, there are many steps that can be taken now to improve the environmental sustainability of MRI and prepare for the effects of climate change. On-going research, technical development, and collaboration with industry partners are needed to achieve further reductions in MRI-related GHG emissions and to decrease the reliance on finite resources. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 6.

MRI Investigation of the Association of Left Atrial and Left Atrial Appendage Hemodynamics with Silent Brain Infarction.

BACKGROUND: Left atrial (LA) myopathy is thought to be associated with silent brain infarctions (SBI) through changes in blood flow hemodynamics leading to thrombogenesis. 4D-flow MRI enables in-vivo hemodynamic quantification in the left atrium (LA) and LA appendage (LAA). PURPOSE: To determine whether LA and LAA hemodynamic and volumetric parameters are associated with SBI. STUDY TYPE: Prospective observational study. POPULATION: A single-site cohort of 125 Participants of the multiethnic study of atherosclerosis (MESA), mean age: 72.3 ± 7.2 years, 56 men. FIELD STRENGTH/SEQUENCE: 1.5T. Cardiac MRI: Cine balanced steady state free precession (bSSFP) and 4D-flow sequences. Brain MRI: T1- and T2-weighted SE and FLAIR. ASSESSMENT: Presence of SBI was determined from brain MRI by neuroradiologists according to routine diagnostic criteria in all participants without a history of stroke based on the MESA database. Minimum and maximum LA volumes and ejection fraction were calculated from bSSFP data. Blood stasis (% of voxels <10 cm/sec) and peak velocity (cm/sec) in the LA and LAA were assessed by a radiologist using an established 4D-flow workflow. STATISTICAL TESTS: Student's t test, Mann-Whitney U test, one-way ANOVA, chi-square test. Multivariable stepwise logistic regression with automatic forward and backward selection. Significance level P < 0.05. RESULTS: 26 (20.8%) had at least one SBI. After Bonferroni correction, participants with SBI were significantly older and had significantly lower peak velocities in the LAA. In multivariable analyses, age (per 10-years) (odds ratio (OR) = 1.99 (95% confidence interval (CI): 1.30-3.04)) and LAA peak velocity (per cm/sec) (OR = 0.87 (95% CI: 0.81-0.93)) were significantly associated with SBI. CONCLUSION: Older age and lower LAA peak velocity were associated with SBI in multivariable analyses whereas volumetric-based measures from cardiac MRI or cardiovascular risk factors were not. Cardiac 4D-flow MRI showed potential to serve as a novel imaging marker for SBI. LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY: Stage 2.

Interval changes in four-dimensional flow-derived in vivo hemodynamics stratify aortic growth in type B aortic dissection patients.

BACKGROUND: Aortic diameter growth in type B aortic dissection (TBAD) is associated with progressive aortic dilation, resulting in increased mortality in patients with both de novo TBAD (dnTBAD) and residual dissection after type A dissection repair (rTAAD). Preemptive thoracic endovascular aortic repair may improve mortality in patients with TBAD, although it is unclear which patients may benefit most from early intervention. In vivo hemodynamic assessment using four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has been used to characterize TBAD patients with growing aortas. In this longitudinal study, we investigated whether changes over time in 4D flow-derived true and false lumen (TL and FL) hemodynamic parameters correlate with aortic growth rate, which is a marker of increased risk. METHODS: We retrospectively identified TBAD patients with baseline and follow-up 4D flow CMR at least 120 days apart. Patients with TBAD intervention before baseline or between scans were excluded. 4D flow CMR data analysis included segmentation of the TL and FL, followed by voxel-wise calculation of TL and FL total kinetic energy (KE), maximum velocity (MV), mean forward flow (FF), and mean reverse flow (RF). Changes over time (Δ) were calculated for all hemodynamic parameters. Maximal diameter in the descending aorta was measured from magnetic resonance angiogram images acquired at the time of 4D flow. Aortic growth rate was defined as the change in diameter divided by baseline diameter and standardized to scan interval. RESULTS: Thirty-two patients met inclusion criteria (age: 56.9 ± 14.1 years, female: 13, n = 19 rTAAD, n = 13 dnTBAD). Mean follow-up time was 538 days (range: 135-1689). Baseline aortic diameter did not correlate with growth rate. In the entire cohort, Δ FL MV (Spearman's rho [rho] = 0.37, p = 0.04) and Δ FL RF (rho = 0.45, p = 0.01) correlated with growth rate. In rTAAD only, Δ FL MV (rho = 0.48, p = 0.04) and Δ FL RF (rho = 0.51, p = 0.03) correlated with growth rate, while in dnTBAD only, Δ TL KE (rho = 0.63, p = 0.02) and Δ TL MV (rho = 0.69, p = 0.01) correlated with growth rate. CONCLUSION: 4D flow-derived longitudinal hemodynamic changes correlate with aortic growth rate in TBAD and may provide additional prognostic value for risk stratification. 4D flow MRI could be integrated into existing imaging protocols to allow for the identification of TBAD patients who would benefit from preemptive surgical or endovascular intervention.

Respiratory-resolved five-dimensional flow cardiovascular magnetic resonance : In-vivo validation and respiratory-dependent flow changes in healthy volunteers and patients with congenital heart disease.

BACKGROUND: This study aimed to validate respiratory-resolved five-dimensional (5D) flow cardiovascular magnetic resonance (CMR) against real-time two-dimensional (2D) phase-contrast MRI, assess the impact of number of respiratory states, and measure the impact of respiration on hemodynamics in congenital heart disease (CHD) patients. METHODS: Respiratory-resolved 5D flow MRI-derived net and peak flow measurements were compared to real-time 2D phase-contrast MRI-derived measurements in 10 healthy volunteers. Pulmonary-to-systemic flow ratios (Qp:Qs) were measured in 19 CHD patients and aortopulmonary collateral burden was measured in 5 Fontan patients. Additionally, the impact of number of respiratory states on measured respiratory-driven net flow changes was investigated in 10 healthy volunteers and 19 CHD patients (shunt physiology, n = 11, single ventricle disease [SVD], n = 8). RESULTS: There was good agreement between 5D flow MRI and real-time 2D phase-contrast-derived net and peak flow. Respiratory-driven changes had a good correlation (rho = 0.64, p < 0.001). In healthy volunteers, fewer than four respiratory states reduced measured respiratory-driven flow changes in veins (5.2 mL/cycle, p < 0.001) and arteries (1.7 mL/cycle, p = 0.05). Respiration drove substantial venous net flow changes in SVD (64% change) and shunt patients (57% change). Respiration had significantly greater impact in SVD patients compared to shunt patients in the right and left pulmonary arteries (46% vs 15%, p = 0.003 and 59% vs 20%, p = 0.002). Qp:Qs varied by 37 ± 24% over respiration in SVD patients and 12 ± 20% in shunt patients. Aortopulmonary collateral burden varied by 118 ± 84% over respiration in Fontan patients. The smallest collateral burden was measured during active inspiration in all patients and the greatest burden was during active expiration in four of five patients. Reduced respiratory resolution blunted measured flow changes in the caval veins of shunt and SVD patients (p < 0.005). CONCLUSIONS: Respiratory-resolved 5D flow MRI measurements agree with real-time 2D phase contrast. Venous measurements are sensitive to number of respiratory states, whereas arterial measurements are more robust. Respiration has a substantial impact on caval vein flow, Qp:Qs, and collateral burden in CHD patients.

Combined free-running four-dimensional anatomical and flow magnetic resonance imaging with native contrast using Synchronization of Neighboring Acquisitions by Physiological Signals.

BACKGROUND: Four-dimensional (4D) flow magnetic resonance imaging (MRI) often relies on the injection of gadolinium- or iron-oxide-based contrast agents to improve vessel delineation. In this work, a novel technique is developed to acquire and reconstruct 4D flow data with excellent dynamic visualization of blood vessels but without the need for contrast injection. Synchronization of Neighboring Acquisitions by Physiological Signals (SyNAPS) uses pilot tone (PT) navigation to retrospectively synchronize the reconstruction of two free-running three-dimensional radial acquisitions, to create co-registered anatomy and flow images. METHODS: Thirteen volunteers and two Marfan syndrome patients were scanned without contrast agent using one free-running fast interrupted steady-state (FISS) sequence and one free-running phase-contrast MRI (PC-MRI) sequence. PT signals spanning the two sequences were recorded for retrospective respiratory motion correction and cardiac binning. The magnitude and phase images reconstructed, respectively, from FISS and PC-MRI, were synchronized to create SyNAPS 4D flow datasets. Conventional two-dimensional (2D) flow data were acquired for reference in ascending (AAo) and descending aorta (DAo). The blood-to-myocardium contrast ratio, dynamic vessel area, net volume, and peak flow were used to compare SyNAPS 4D flow with Native 4D flow (without FISS information) and 2D flow. A score of 0-4 was given to each dataset by two blinded experts regarding the feasibility of performing vessel delineation. RESULTS: Blood-to-myocardium contrast ratio for SyNAPS 4D flow magnitude images (1.5 ± 0.3) was significantly higher than for Native 4D flow (0.7 ± 0.1, p < 0.01) and was comparable to 2D flow (2.3 ± 0.9, p = 0.02). Image quality scores of SyNAPS 4D flow from the experts (M.P.: 1.9 ± 0.3, E.T.: 2.5 ± 0.5) were overall significantly higher than the scores from Native 4D flow (M.P.: 1.6 ± 0.6, p = 0.03, E.T.: 0.8 ± 0.4, p < 0.01) but still significantly lower than the scores from the reference 2D flow datasets (M.P.: 2.8 ± 0.4, p < 0.01, E.T.: 3.5 ± 0.7, p < 0.01). The Pearson correlation coefficient between the dynamic vessel area measured on SyNAPS 4D flow and that from 2D flow was 0.69 ± 0.24 for the AAo and 0.83 ± 0.10 for the DAo, whereas the Pearson correlation between Native 4D flow and 2D flow measurements was 0.12 ± 0.48 for the AAo and 0.08 ± 0.39 for the DAo. Linear correlations between SyNAPS 4D flow and 2D flow measurements of net volume (r2 = 0.83) and peak flow (r2 = 0.87) were larger than the correlations between Native 4D flow and 2D flow measurements of net volume (r2 = 0.79) and peak flow (r2 = 0.76). CONCLUSION: The feasibility and utility of SyNAPS were demonstrated for joint whole-heart anatomical and flow MRI without requiring electrocardiography gating, respiratory navigators, or contrast agents. Using SyNAPS, a high-contrast anatomical imaging sequence can be used to improve 4D flow measurements that often suffer from poor delineation of vessel boundaries in the absence of contrast agents.

Cardiac magnetic resonance imaging in detection of progressive graft dysfunction in pediatric heart transplantation.

BACKGROUND: Chronic graft failure (CGF) in pediatric heart transplant (PHT) is multifactorial and may present with findings of fibrosis and microvessel disease (MVD) on endomyocardial biopsy (EMB). There is no optimal CGF surveillance method. We evaluated associations between cardiac magnetic resonance imaging (CMR) and historical/EMB correlates of CGF to assess CMR's utility as a surveillance method. METHODS: Retrospective analysis of PHT undergoing comprehensive CMR between September 2015 and January 2022 was performed. EMB within 6 months was graded for fibrosis (scale 0-5) and MVD (number of capillaries with stenotic wall thickening per field of view). Correlation analysis and logistic regression were performed. RESULTS: Forty-seven PHT with median age at CMR of 15.7 years (11.6, 19.3) and time from transplant of 6.4 years (4.1, 11.0) were studied. Cardiac allograft vasculopathy (CAV) was present in 11/44 (22.0%) and historical rejection in 14/41 (34.2%). CAV was associated with higher global T2 (49.0 vs. 47.0 ms; p = 0.038) and peak T2 (57.0 vs. 53.0 ms; p = 0.013) on CMR. Historical rejection was associated with higher global T2 (49.0 vs. 47.0 ms; p = 0.007) and peak T2 (57.0 vs. 53.0 ms; p = 0.03) as well as global extracellular volume (31.0 vs. 26.3%; p = 0.03). Higher fibrosis score on EMB correlated with smaller indexed left ventricular mass (rho = -0.34; p = 0.019) and greater degree of MVD with lower indexed left ventricular end-diastolic volume (rho = -0.35; p = 0.017). CONCLUSION: Adverse ventricular remodeling and abnormal myocardial characteristics on CMR are present in PHT with CAV, historical rejection, as well as greater fibrosis and MVD on EMB. CMR has the potential use for screening of CGF.

Focused navigation for respiratory-motion-corrected free-running radial 4D flow MRI.

PURPOSE: To validate a respiratory motion correction method called focused navigation (fNAV) for free-running radial whole-heart 4D flow MRI. METHODS: Using fNAV, respiratory signals derived from radial readouts are converted into three orthogonal displacements, which are then used to correct respiratory motion in 4D flow datasets. Hundred 4D flow acquisitions were simulated with non-rigid respiratory motion and used for validation. The difference between generated and fNAV displacement coefficients was calculated. Vessel area and flow measurements from 4D flow reconstructions with (fNAV) and without (uncorrected) motion correction were compared to the motion-free ground-truth. In 25 patients, the same measurements were compared between fNAV 4D flow, 2D flow, navigator-gated Cartesian 4D flow, and uncorrected 4D flow datasets. RESULTS: For simulated data, the average difference between generated and fNAV displacement coefficients was 0.04 ± $$ \pm $$ 0.32 mm and 0.31 ± $$ \pm $$ 0.35 mm in the x and y directions, respectively. In the z direction, this difference was region-dependent (0.02 ± $$ \pm $$ 0.51 mm up to 5.85 ± $$ \pm $$ 3.41 mm). For all measurements (vessel area, net volume, and peak flow), the average difference from ground truth was higher for uncorrected 4D flow datasets (0.32 ± $$ \pm $$ 0.11 cm2 , 11.1 ± $$ \pm $$ 3.5 mL, and 22.3 ± $$ \pm $$ 6.0 mL/s) than for fNAV 4D flow datasets (0.10 ± $$ \pm $$ 0.03 cm2 , 2.6 ± $$ \pm $$ 0.7 mL, and 5.1 ± 0 $$ \pm 0 $$ .9 mL/s, p < 0.05). In vivo, average vessel area measurements were 4.92 ± $$ \pm $$ 2.95 cm2 , 5.06 ± $$ \pm $$ 2.64 cm2 , 4.87 ± $$ \pm $$ 2.57 cm2 , 4.87 ± $$ \pm $$ 2.69 cm2 , for 2D flow and fNAV, navigator-gated and uncorrected 4D flow datasets, respectively. In the ascending aorta, all 4D flow datasets except for the fNAV reconstruction had significantly different vessel area measurements from 2D flow. Overall, 2D flow datasets demonstrated the strongest correlation to fNAV 4D flow for both net volume (r2  = 0.92) and peak flow (r2  = 0.94), followed by navigator-gated 4D flow (r2  = 0.83 and r2  = 0.86, respectively), and uncorrected 4D flow (r2  = 0.69 and r2  = 0.86, respectively). CONCLUSION: fNAV corrected respiratory motion in vitro and in vivo, resulting in fNAV 4D flow measurements that are comparable to those derived from 2D flow and navigator-gated Cartesian 4D flow datasets, with improvements over those from uncorrected 4D flow.

Cine MRI-Derived Radiomics Features of the Cardiac Blood Pool: Periodicity, Specificity, and Reproducibility.

BACKGROUND: Although radiomics features of the left ventricular wall have been used to assess cardiac diseases, radiomics features of the cardiac blood pool have been relatively ignored. PURPOSE: To test the hypothesis that cine MRI-derived radiomics features of the cardiac blood pool are associated with cardiac function and motion. STUDY TYPE: Retrospective. POPULATION: A total of 26 healthy volunteers (51.2 ± 15.6 years; 17 males). FIELD STRENGTH/SEQUENCE: A 1.5 T/balanced steady-state free precession (bSSFP). ASSESSMENT: The radiomics features (107 features in seven classes) of the blood pool of the left/right ventricle/atrium (LV/RV/LA/RA) were extracted on four-chamber cine images (25 phases). Conventional cardiac function parameters (volumes, ejection fraction [EF] and longitudinal strain) were assessed in each cardiac chamber. Intraobserver- and interobserver agreements of radiomics features of all chambers acquired at all phases were assessed, as well as scan-rescan agreement in a subset of 13 volunteers. STATISTICAL TESTS: Pearson correlation coefficients (r) were used to assess the associations between peak values of radiomics features and end-diastolic (or maximal) volume, end-systolic (or minimal) volume, EF, and longitudinal strain of corresponding chambers. Good intraobserver, interobserver, and scan-rescan agreements for radiomics features acquired were defined as intraclass correlation coefficient (ICC) > 0.7 or coefficient of variation (CoV) < 20%. RESULTS: Most radiomics features of the blood pool varied periodically throughout the cardiac cycle. Peak values of chamber-specific blood pool radiomics features were correlated with traditional cardiac function and motion indices of corresponding chambers (r: 0.4-0.87). Ninety-three (87%), 86 (80%), and 73 (68%) radiomics features demonstrated good intraobserver, interobserver, and scan-rescan reproducibility, respectively. CONCLUSION: Cine MRI-derived radiomics features within LV/RV/LA/RA are associated with traditional cardiac function and motion indices of corresponding chambers and may have the potential to become novel quantitative imaging biomarkers in cardiovascular medicine. EVIDENCE LEVEL: 3. TECHNICAL EFFICACY: 1.

Right Ventricular Remodeling Assessed by MRI in Duchenne Muscular Dystrophy.

BACKGROUND: In Duchenne muscular dystrophy (DMD), the right ventricle (RV) tends to be relatively well preserved, but characterization remains difficult due to its complex architecture. Tissue phase mapping (TPM) is a phase contrast cine MRI technique that allows for multidirectional assessment of myocardial velocities. PURPOSE: To use TPM to elucidate relationships between myocardial structure, function, and clinical variables in DMD. STUDY TYPE: Retrospective. SUBJECTS: A total of 20 patients with muscular dystrophy (median age: 16 years); 18 age-matched normal controls (median age: 15 years). FIELD STRENGTH/SEQUENCE: Three-directional velocity encoded cine gradient echo sequence (TPM) at 1.5 T, balanced steady-state free procession (bSSFP), T1 mapping with extracellular volume (ECV), and late gadolinium enhancement (LGE). ASSESSMENT: TPM in basal, mid, and apical short-axis planes was performed as part of a standard MRI study with collection of clinical data. Radial, circumferential, and longitudinal velocities (Vr, Vφ, and Vz, respectively) and corresponding time to peak (TTP) velocities were quantified from TPM and used to calculate RV twist as well as intraventricular and interventricular dyssynchrony. The correlations between TPM velocities, myocardial structure/function, and clinical variables were assessed. STATISTICAL TEST: Unpaired t-test, Wilcoxon rank-sum test, Bland-Altman analyses were used for comparisons between DMD patients and controls and between DMD subgroups. Pearson's test was used for correlations (r). Significance level: P < 0.05. RESULTS: Compared to controls, DMD patients had preserved RV ejection fraction (RVEF 53% ± 8%) but significantly increased interventricular dyssynchrony (Vφ: 0.49 ± 0.21 vs. 0.72 ± 0.17). Within the DMD cohort, RV dyssynchrony significantly increased with lower LV ejection fraction (intraventricular Vr and Vz: r = -0.49; interventricular Vz: r = 0.48). In addition, RV intraventricular dyssynchrony significantly increased with older age (Vz: r = 0.67). DATA CONCLUSION: RV remodeling in DMD occurs in the context of preserved RVEF. Within DMD, this abnormal RV deformation is associated with older age and decreased LVEF. EVIDENCE LEVEL: 4. TECHNICAL EFFICACY: Stage 2.

Comparison of Improved Unidirectional Dual Velocity-Encoding MRI Methods.

BACKGROUND: In phase-contrast (PC) MRI, several dual velocity encoding methods have been proposed to robustly increase velocity-to-noise ratio (VNR), including a standard dual-VENC (SDV), an optimal dual-VENC (ODV), and bi- and triconditional methods. PURPOSE: To develop a correction method for the ODV approach and to perform a comparison between methods. STUDY TYPE: Case-control study. POPULATION: Twenty-six volunteers. FIELD STRENGTH/SEQUENCE: 1.5 T phase-contrast MRI with VENCs of 50, 75, and 150 cm/second. ASSESSMENT: Since we acquired single-VENC protocols, we used the background phase from high-VENC (VENCH ) to reconstruct the low-VENC (VENCL ) phase. We implemented and compared the unwrapping methods for different noise levels and also developed a correction of the ODV method. STATISTICAL TESTS: Shapiro-Wilk's normality test, two-way analysis of variance with homogeneity of variances was performed using Levene's test, and the significance level was adjusted by Tukey's multiple post hoc analysis with Bonferroni (P < 0.05). RESULTS: Statistical analysis revealed no extreme outliers, normally distributed residuals, and homogeneous variances. We found statistically significant interaction between noise levels and the unwrapping methods. This implies that the number of non-unwrapped pixels increased with the noise level. We found that for ß = VENCL /VENCH  = 1/2, unwrapping methods were more robust to noise. The post hoc test showed a significant difference between the ODV corrected and the other methods, offering the best results regarding the number of unwrapped pixels. DATA CONCLUSIONS: All methods performed similarly without noise, but the ODV corrected method was more robust to noise at the price of a higher computational time. LEVEL OF EVIDENCE: 4 TECHNICAL EFFICACY STAGE: 1.

Assessment of Beat-To-Beat Variability in Left Atrial Hemodynamics Using Real Time Phase Contrast MRI in Patients With Atrial Fibrillation.

BACKGROUND: Hemodynamic assessment of left atrial (LA) flow using phase contrast MRI provides insight into thromboembolic risk in atrial fibrillation (AF). However, conventional flow imaging techniques are averaged over many heartbeats. PURPOSE: To evaluate beat-to-beat variability and LA hemodynamics in patients with AF using real time phase contrast (RTPC) MRI. STUDY TYPE: Prospective. SUBJECTS: Thirty-five patients with history of AF (68 ± 10 years, nine female), 10 healthy controls (57 ± 19 years, four female). FIELD STRENGTH/SEQUENCE: 5T, 2D RTPC with through-plane velocity-encoded gradient echo sequence and 4D flow MRI with three-directional velocity-encoded gradient echo sequence. ASSESSMENT: RTPC was continuously acquired for a mid-LA slice in all subjects. 4D flow data were interpolated at the RTPC location and normally projected for comparison with RTPC. RR intervals extracted from RTPC were used to calculate heart rate variability (HRV = interquartile range over median × 100%). Patients were classified into low (<9.7%) and high (>9.7%) HRV groups. LA peak/mean velocity and stasis (%velocities < 5.8 cm/sec) were calculated from segmented 2D images. Variability in RTPC flow metrics was quantified by coefficient of variation (CV) over all cycles. STATISTICAL TESTS: Pearson's correlation coefficient (r), Bland-Altman analysis, Kruskal-Wallis test. A P value < 0.05 was considered statistically significant. RESULTS: RTPC and 4D flow measurements were strongly/significantly correlated for all hemodynamic parameters (R2  = 0.75-0.83) in controls. Twenty-four patients had low HRV (mean = 4 ± 2%) and 11 patients had high HRV (27 ± 9%). In patients, increased HRV was significantly correlated with CV of peak velocity (r = 0.67), mean velocity (r = 0.51), and stasis (r = 0.41). A stepwise decrease in peak/mean velocity and increase in stasis was observed when comparing controls vs. low HRV vs. high HRV groups. Mean velocity and stasis differences were significant for control vs. high HRV groups. CONCLUSIONS: RTPC may be suitable for assessing the impact of HRV on hemodynamics and provide insight for AF management in highly arrhythmic patients. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 2.

Quantitative Assessment of Regional Pulmonary Transit Times in Pulmonary Hypertension.

BACKGROUND: Pulmonary hypertension (PH) contributes to restricted flow through the pulmonary circulation characterized by elevated mean pulmonary artery pressure acquired from invasive right heart catheterization (RHC). MRI may provide a noninvasive alternative for diagnosis and characterization of PH. PURPOSE: To characterize PH via quantification of regional pulmonary transit times (rPTT). STUDY TYPE: Retrospective. POPULATION: A total of 43 patients (58% female); 24 controls (33% female). RHC-confirmed patients classified as World Health Organization (WHO) subgroups 1-4. FIELD STRENGTH/SEQUENCE: A 1.5 T/time-resolved contrast-enhanced MR Angiography (CE-MRA). ASSESSMENT: CE-MRA data volumes were combined into a 4D matrix (3D resolution + time). Contrast agent arrival time was defined as the peak in the signal-intensity curve generated for each voxel. Average arrival times within a vessel region of interest (ROI) were normalized to the main pulmonary artery ROI (t0 ) for eight regions to define rPTT for all subjects. Subgroup analysis included grouping the four arterial and four venous regions. Intraclass correlation analysis completed for reproducibility. STATISTICAL TESTS: Analysis of covariance with age as covariate. A priori Student's t-tests or Wilcoxon rank-sum test; α = 0.05. Results compared to controls unless noted. Significant without listing P value. ICC ran as two-way absolute agreement model with two observers. RESULTS: PH patients demonstrated elevated rPTT in all vascular regions; average rPTT increase in arterial and venous branches was 0.85 ± 0.15 seconds (47.7%) and 1.0 ± 0.18 seconds (16.9%), respectively. Arterial rPTT was increased for all WHO subgroups; venous regions were elevated for subgroups 2 and 4 (group 1, P = 0.86; group 3, P = 0.32). No significant rPTT differences were found between subgroups (P = 0.094-0.94). Individual vessel ICC values ranged from 0.58 to 0.97. DATA CONCLUSION: Noninvasive assessment of PH using standard-of-care time-resolved CE-MRA can detect increased rPTT in PH patients of varying phenotypes compared to controls. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY: Stage 3.

Highly Accelerated Compressed-Sensing 4D Flow for Intracardiac Flow Assessment.

BACKGROUND: Four-dimensional (4D) flow MRI allows for the quantification of complex flow patterns; however, its clinical use is limited by its inherently long acquisition time. Compressed sensing (CS) is an acceleration technique that provides substantial reduction in acquisition time. PURPOSE: To compare intracardiac flow measurements between conventional and CS-based highly accelerated 4D flow acquisitions. STUDY TYPE: Prospective. SUBJECTS: Fifty healthy volunteers (28.0 ± 7.1 years, 24 males). FIELD STRENGTH/SEQUENCE: Whole heart time-resolved 3D gradient echo with three-directional velocity encoding (4D flow) with conventional parallel imaging (factor 3) as well as CS (factor 7.7) acceleration at 3 T. ASSESSMENT: 4D flow MRI data were postprocessed by applying a valve tracking algorithm. Acquisition times, flow volumes (mL/cycle) and diastolic function parameters (ratio of early to late diastolic left ventricular peak velocities [E/A] and ratio of early mitral inflow velocity to mitral annular early diastolic velocity [E/e']) were quantified by two readers. STATISTICAL TESTS: Paired-samples t-test and Wilcoxon rank sum test to compare measurements. Pearson correlation coefficient (r), Bland-Altman-analysis (BA) and intraclass correlation coefficient (ICC) to evaluate agreement between techniques and readers. A P value < 0.05 was considered statistically significant. RESULTS: A significant improvement in acquisition time was observed using CS vs. conventional accelerated acquisition (6.7 ± 1.3 vs. 12.0 ± 1.3 min). Net forward flow measurements for all valves showed good correlation (r > 0.81) and agreement (ICCs > 0.89) between conventional and CS acceleration, with 3.3%-8.3% underestimation by the CS technique. Evaluation of diastolic function showed 3.2%-17.6% error: E/A 2.2 [1.9-2.4] (conventional) vs. 2.3 [2.0-2.6] (CS), BA bias 0.08 [-0.81-0.96], ICC 0.82; and E/e' 4.6 [3.9-5.4] (conventional) vs. 3.8 [3.4-4.3] (CS), BA bias -0.90 [-2.31-0.50], ICC 0.89. DATA CONCLUSION: Analysis of intracardiac flow patterns and evaluation of diastolic function using a highly accelerated 4D flow sequence prototype is feasible, but it shows underestimation of flow measurements by approximately 10%. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

Global Aortic Pulse Wave Velocity is Unchanged in Bicuspid Aortopathy With Normal Valve Function but Elevated in Patients With Aortic Valve Stenosis: Insights From a 4D Flow MRI Study of 597 Subjects.

BACKGROUND: Aortopathy is common with bicuspid aortic valve (BAV), and underlying intrinsic tissue abnormalities are believed causative. Valve-mediated hemodynamics are altered in BAV and may contribute to aortopathy and its progression. The contribution of intrinsic tissue defects versus altered hemodynamics to aortopathy progression is not known. PURPOSE: To investigate relative contributions of tissue-innate versus hemodynamics in progression of BAV aortopathy. STUDY TYPE: Retrospective. SUBJECTS: Four hundred seventy-three patients with aortic dilatation (diameter ≥40 mm; comprised of 281 BAV with varied AS severity, 192 tricuspid aortic valve [TAV] without AS) and 124 healthy controls. Subjects were 19-91 years (141/24% female). FIELD STRENGTH/SEQUENCE: 1.5T, 3T; time-resolved gradient-echo 3D phase-contrast (4D flow) MRI. ASSESSMENT: A surrogate measure for global aortic wall stiffness, pulse wave velocity (PWV), was quantified from MRI by standardized, automated technique based on through-plane flow cross-correlation maximization. Comparisons were made between BAV patients with aortic dilatation and varying aortic valve stenosis (AS) severity and healthy subjects and aortopathy patients with normal TAV. STATISTICAL TESTS: Multivariable regression, analysis of covariance (ANCOVA), Tukey's, student's (t), Mann-Whitney (U) tests, were used with significance levels P < 0.05 or P < 0.01 for post-hoc Bonferroni-corrected t/U tests. Bland-Altman and ICC calculations were performed. RESULTS: Multivariable regression showed age with the most significant association for increased PWV in all groups (increase 0.073-0.156 m/sec/year, R2  = 0.30-48). No significant differences in aortic PWV were observed between groups without AS (P = 0.20-0.99), nor were associations between PWV and regurgitation or Sievers type observed (P = 0.60, 0.31 respectively). In contrast, BAV AS patients demonstrated elevated PWV and a significant relationship for AS severity with increased PWV (covariate: age, R2  = 0.48). BAV and TAV patients showed no association between aortic diameter and PWV (P = 0.73). DATA CONCLUSION: No significant PWV differences were observed between BAV patients with normal valve function and control groups. However, AS severity and age in BAV patients were directly associated with PWV increases. EVIDENCE LEVEL: 3 TECHNICAL EFFICACY: Stage 3.

Hemodynamic Evaluation of Type B Aortic Dissection Using Compressed Sensing Accelerated 4D Flow MRI.

BACKGROUND: 4D Flow MRI is a quantitative imaging technique to evaluate blood flow patterns; however, it is unclear how compressed sensing (CS) acceleration would impact aortic hemodynamic quantification in type B aortic dissection (TBAD). PURPOSE: To investigate CS-accelerated 4D Flow MRI performance compared to GRAPP-accelerated 4D Flow MRI (GRAPPA) to evaluate aortic hemodynamics in TBAD. STUDY TYPE: Prospective. POPULATION: Twelve TBAD patients, two volunteers. FIELD STRENGTH/SEQUENCE: 1.5T, 3D time-resolved cine phase-contrast gradient echo sequence. ASSESSMENT: GRAPPA (acceleration factor [R] = 2) and two CS-accelerated (R = 7.7 [CS7.7] and 10.2 [CS10.2]) 4D Flow MRI scans were acquired twice for interscan reproducibility assessment. Voxelwise kinetic energy (KE), peak velocity (PV), forward flow (FF), reverse flow (RF), and stasis were calculated. Plane-based mid-lumen flows were quantified. Imaging times were recorded. TESTS: Repeated measures analysis of variance, Pearson correlation coefficients (r), intraclass correlation coefficients (ICC). P < 0.05 indicated statistical significance. RESULTS: The KE and FF in true lumen (TL) and PV in false lumen (FL) did not show difference among three acquisition types (P = 0.818, 0.065, 0.284 respectively). The PV and stasis in TL were higher, KE, FF, and RF in FL were lower, and stasis was higher in GRAPPA compared to CS7.7 and CS10.2. The RF was lower in GRAPPA compared to CS10.2. The correlation coefficients were strong in TL (r = [0.781-0.986]), and low to strong in FL (r = [0.347-0.948]). The ICC levels demonstrated moderate to excellent interscan reproducibility (0.732-0.989). The FF and net flow in mid-descending aorta TL were significantly different between CS7.7 and CS10.2. CONCLUSION: CS-accelerated 4D Flow MRI has potential for clinical utilization with shorter scan times in TBAD. Our results suggest similar hemodynamic trends between acceleration types, but CS-acceleration impacts KE, FF, RF, and stasis more in FL. EVIDENCE LEVEL: 1 Technical Efficacy: Stage 2.

Fully-automated deep learning-based flow quantification of 2D CINE phase contrast MRI.

OBJECTIVES: Time-resolved, 2D-phase-contrast MRI (2D-CINE-PC-MRI) enables in vivo blood flow analysis. However, accurate vessel contour delineation (VCD) is required to achieve reliable results. We sought to evaluate manual analysis (MA) compared to the performance of a deep learning (DL) application for fully-automated VCD and flow quantification and corrected semi-automated analysis (corSAA). METHODS: We included 97 consecutive patients (age = 52.9 ± 16 years, 41 female) with 2D-CINE-PC-MRI imaging on 1.5T MRI systems at sinotubular junction (STJ), and 28/97 also received 2D-CINE-PC at main pulmonary artery (PA). A cardiovascular radiologist performed MA (reference) and corSAA (built-in tool) in commercial software for all cardiac time frames (median: 20, total contours per analysis: 2358 STJ, 680 PA). DL-analysis automatically performed VCD, followed by net flow (NF) and peak velocity (PV) quantification. Contours were compared using Dice similarity coefficients (DSC). Discrepant cases (> ± 10 mL or > ± 10 cm/s) were reviewed in detail. RESULTS: DL was successfully applied to 97% (121/125) of the 2D-CINE-PC-MRI series (STJ: 95/97, 98%, PA: 26/28, 93%). Compared to MA, mean DSC were 0.91 ± 0.02 (DL), 0.94 ± 0.02 (corSAA) at STJ, and 0.85 ± 0.08 (DL), 0.93 ± 0.02 (corSAA) at PA; this indicated good to excellent DL-performance. Flow quantification revealed similar NF at STJ (p = 0.48) and PA (p > 0.05) between methods while PV assessment was significantly different (STJ: p < 0.001, PA: p = 0.04). A detailed review showed noisy voxels in MA and corSAA impacted PV results. Overall, DL analysis compared to human assessments was accurate in 113/121 (93.4%) cases. CONCLUSIONS: Fully-automated DL-analysis of 2D-CINE-PC-MRI provided flow quantification at STJ and PA at expert level in > 93% of cases with results being available instantaneously. KEY POINTS: • Deep learning performed flow quantification on clinical 2D-CINE-PC series at the sinotubular junction and pulmonary artery at the expert level in > 93% of cases. • Location detection and contouring of the vessel boundaries were performed fully-automatic with results being available instantaneously compared to human assessments which approximately takes three minutes per location. • The evaluated tool indicates usability in daily practice.

Late-gadolinium enhancement is common in older pediatric heart transplant recipients and is associated with lower ejection fraction.

BACKGROUND: Chronic graft failure and cumulative rejection history in pediatric heart transplant recipients (PHTR) are associated with myocardial fibrosis on endomyocardial biopsy (EMB). Cardiovascular magnetic resonance imaging (CMR) is a validated, non-invasive method to detect myocardial fibrosis via the presence of late gadolinium enhancement (LGE). In adult heart transplant recipients, LGE is associated with increased risk of future adverse clinical events including hospitalization and death. We describe the prevalence, pattern, and extent of LGE on CMR in a cohort of PHTR and its associations with recipient and graft characteristics. METHODS: This was a retrospective study of consecutive PHTR who underwent CMR over a 6-year period at a single center. Two independent reviewers assessed the presence and distribution of left ventricular (LV) LGE using the American Heart Association (AHA) 17-segment model. LGE quantification was performed on studies with visible fibrosis (LGE+). Patient demographics, clinical history, and CMR-derived volumetry and ejection fractions were obtained. RESULTS: Eighty-one CMR studies were performed on 59 unique PHTR. Mean age at CMR was 14.8 ± 6.2 years; mean time since transplant was 7.3 ± 5.0 years. The CMR indication was routine surveillance (without a clinical concern based on laboratory parameters, echocardiography, or cardiac catheterization) in 63% (51/81) of studies. LGE was present in 36% (29/81) of PHTR. In these LGE + studies, patterns included inferoseptal in 76% of LGE + studies (22/29), lateral wall in 41% (12/29), and diffuse, involving > 4 AHA segments, in 21% (6/29). The mean LV LGE burden as a percentage of myocardial mass was 18.0 ± 9.0%. When reviewing only the initial CMR per PHTR (n = 59), LGE + patients were older (16.7 ± 2.9 vs. 12.8 ± 4.6 years, p = 0.001), with greater time since transplant (8.3 ± 5.4 vs. 5.7 ± 3.9 years, p = 0.041). These patients demonstrated higher LV end-systolic volume index (LVESVI) (34.7 ± 11.7 vs. 28.7 ± 6.1 ml/m2, p = 0.011) and decreased LV ejection fraction (LVEF) (56.2 ± 8.1 vs. 60.6 ± 5.3%, p = 0.015). There were no significant differences in history of moderate/severe rejection (p = 0.196) or cardiac allograft vasculopathy (CAV) (p = 0.709). CONCLUSIONS: LV LGE was present in approximately one third of PHTR, more commonly in older patients with longer time since transplantation. Grafts with LGE have lower LVEF. CMR-derived LGE may aid in surveillance of chronic graft failure in PHTR.