4D flow cardiovascular magnetic resonance consensus statement

Mitral Annular Disease at Cardiac MRI: What to Know and Look For

Accurate evaluation of the mitral valve (MV) apparatus is essential for understanding the mechanisms of MV disease across various clinical scenarios. The mitral annulus (MA) is a complex and crucial structure that supports MV function; however, conventional imaging techniques have limitations in fully capturing the entirety of the MA. Moreover, recognizing annular changes might aid in identifying patients who may benefit from advanced cardiac imaging and interventions. Multimodality cardiovascular imaging plays a major role in the diagnosis, prognosis, and management of MV disease. Transthoracic echocardiography is the first-line modality for evaluation of the MA, but it has limitations. Cardiac MRI (CMR) has emerged as a robust imaging modality for assessing annular changes, with distinct advantages over other imaging techniques, including accurate flow and volumetric quantification and assessment of variations in the measurements and shape of the MA during the cardiac cycle. Mitral annular disjunction (MAD) is defined as atrial displacement of the hinge point of the MV annulus away from the ventricular myocardium, a condition that is now more frequently diagnosed and studied owing to recent technical advances in cardiac imaging. However, several unresolved issues regarding MAD, such as the functional significance of pathologic disjunction and how this disjunction advances in the clinical course, require further investigation. The authors review the role of CMR in the assessment of MA disease, with a focus on MAD and its functional implications in MV prolapse and mitral regurgitation. ©RSNA, 2024 Supplemental material is available for this article. See the invited commentary by Stojanovska and Fujikura in this issue.

Impact of 4D-Flow CMR Parameters on Functional Evaluation of Fontan Circulation

We sought to evaluate the potential clinical role of 4D-flow cardiac magnetic resonance (CMR)-derived energetics and flow parameters in a cohort of patients' post-Fontan palliation. In patients with Fontan circulation who underwent 4D-Flow CMR, streamlines distribution was evaluated, as well a 4D-flow CMR-derived energetics parameters as kinetic energy (KE) and energy loss (EL) normalized by volume. EL/KE index as a marker of flow efficiency was also calculated. Cardiopulmonary exercise test (CPET) was also performed in a subgroup of patients. The population study included 55 patients (mean age 22 ± 11 years). The analysis of the streamlines revealed a preferential distribution of the right superior vena cava flow for the right pulmonary artery (62.5 ± 35.4%) and a mild preferential flow for the left pulmonary artery (52.3 ± 40.6%) of the inferior vena cave-pulmonary arteries (IVC-PA) conduit. Patients with heart failure (HF) presented lower IVC/PA-conduit flow (0.75 ± 0.5 vs 1.3 ± 0.5 l/min/m2, p = 0.004) and a higher mean flow-jet angle of the IVC-PA conduit (39.2 ± 22.8 vs 15.2 ± 8.9, p < 0.001) than the remaining patients. EL/KE index correlates inversely with VO2/kg/min: R: - 0.45, p = 0.01 peak, minute ventilation (VE) R: - 0.466, p < 0.01, maximal voluntary ventilation: R:0.44, p = 0.001 and positively with the physiological dead space to the tidal volume ratio (VD/VT) peak: R: 0.58, p < 0.01. From our data, lower blood flow in IVC/PA conduit and eccentric flow was associated with HF whereas higher EL/KE index was associated with reduced functional capacity and impaired lung function. Larger studies are needed to confirm our results and to further improve the prognostic role of the 4D-Flow CMR in this challenging population.

Tricuspid valve flow measurement using a deep learning framework for automated valve-tracking 2D phase contrast

PURPOSE

Tricuspid valve flow velocities are challenging to measure with cardiovascular MR, as the rapidly moving valvular plane prohibits direct flow evaluation, but they are vitally important to diastolic function evaluation. We developed an automated valve-tracking 2D method for measuring flow through the dynamic tricuspid valve.

METHODS

Nine healthy subjects and 2 patients were imaged. The approach uses a previously trained deep learning network, TVnet, to automatically track the tricuspid valve plane from long-axis cine images. Subsequently, the tracking information is used to acquire 2D phase contrast (PC) with a dynamic (moving) acquisition plane that tracks the valve. Direct diastolic net flows evaluated from the dynamic PC sequence were compared with flows from 2D-PC scans acquired in a static slice localized at the end-systolic valve position, and also ventricular stroke volumes (SVs) using both planimetry and 2D PC of the great vessels.

RESULTS

The mean tricuspid valve systolic excursion was 17.8 ± 2.5 mm. The 2D valve-tracking PC net diastolic flow showed excellent correlation with SV by right-ventricle planimetry (bias ± 1.96 SD = -0.2 ± 10.4 mL, intraclass correlation coefficient [ICC] = 0.92) and aortic PC (-1.0 ± 13.8 mL, ICC = 0.87). In comparison, static tricuspid valve 2D PC also showed a strong correlation but had greater bias (p = 0.01) versus the right-ventricle SV (10.6 ± 16.1 mL, ICC = 0.61). In most (8 of 9) healthy subjects, trace regurgitation was measured at begin-systole. In one patient, valve-tracking PC displayed a high-velocity jet (380 cm/s) with maximal velocity agreeing with echocardiography.

CONCLUSION

Automated valve-tracking 2D PC is a feasible route toward evaluation of tricuspid regurgitant velocities, potentially solving a major clinical challenge.

Pulmonary 4D-flow MRI imaging in landrace pigs under rest and stress

4D-flow MRI is a promising technique for assessing vessel hemodynamics. However, its utilization is currently limited by the lack of reference values, particularly for pulmonary vessels. In this work, we have analysed flow and velocity in the pulmonary trunk (PT), left and right pulmonary arteries (LPA and RPA, respectively) in Landrace pigs at both rest and stress through the software MEVISFlow. Nine healthy Landrace pigs were acutely instrumented closed-chest and transported to the CMR facility for evaluation. After rest measurements, dobutamine was administered to achieve a 25% increase in heart rate compared to rest. 4D-flow MRI images have been analysed through MEVISFlow by two independent observers. Inter- and intra-observer reproducibility was quantified using intraclass correlation coefficient. A significant difference between rest and stress regarding flow and velocity in all the pulmonary vessels was observed. Mean flow increased 55% in PT, 75% in LPA and 40% in RPA. Mean peak velocity increased 55% in PT, 75% in LPA and 66% in RPA. A good-to-excellent reproducibility was observed in rest and stress for flow measurements in all three arteries. An excellent reproducibility for velocity was found in PT at rest and stress, a good one for LPA and RPA at rest, while poor reproducibility was found at stress. The current study showed that pulmonary flow and velocity assessed through 4D-flow MRI follow the physiological alterations during cardiac cycle and after stress induced by dobutamine. A clinical translation to assess pulmonary diseases with 4D-flow MRI under stress conditions needs investigation.

Patient-Specific Numerical Simulations of Coronary Artery Hemodynamics and Biomechanics: A Pathway to Clinical Use

PURPOSE

Numerical models that simulate the behaviors of the coronary arteries have been greatly improved by the addition of fluid-structure interaction (FSI) methods. Although computationally demanding, FSI models account for the movement of the arterial wall and more adequately describe the biomechanical conditions at and within the arterial wall. This offers greater physiological relevance over Computational Fluid Dynamics (CFD) models, which assume the walls do not move or deform. Numerical simulations of patient-specific cases have been greatly bolstered by the use of imaging modalities such as Computed Tomography Angiography (CTA), Magnetic Resonance Imaging (MRI), Optical Coherence Tomography (OCT), and Intravascular Ultrasound (IVUS) to reconstruct accurate 2D and 3D representations of artery geometries. The goal of this study was to conduct a comprehensive review on CFD and FSI models on coronary arteries, and evaluate their translational potential.

METHODS

This paper reviewed recent work on patient-specific numerical simulations of coronary arteries that describe the biomechanical conditions associated with atherosclerosis using CFD and FSI models. Imaging modality for geometry collection and clinical applications were also discussed.

RESULTS

Numerical models using CFD and FSI approaches are commonly used to study biomechanics within the vasculature. At high temporal and spatial resolution (compared to most cardiac imaging modalities), these numerical models can generate large amount of biomechanics data.

CONCLUSIONS

Physiologically relevant FSI models can more accurately describe atherosclerosis pathogenesis, and help to translate biomechanical assessment to clinical evaluation.

Repaired coarctation of the aorta does not affect four-dimensional flow metrics in bicuspid aortic valve disease

OBJECTIVES

The objective of this study was primarily to compare four-dimensional flow magnetic resonance imaging metrics in the ascending aorta (AA) of patients with right-left fusion type bicuspid aortic valve (RL-BAV) and repaired coarctation of the aorta (CoA) to RL-BAV without CoA. Metrics of patients with RL-BAV were also compared to the matched group of patients with common tricuspid aortic valve (TAV).

METHODS

Eleven patients with RL-BAV and CoA, 11 patients with RL-BAV without CoA and 22 controls with TAV were investigated. Peak velocity (cm/s), peak flow (ml/s) and flow displacement (%) were analysed at 5 pre-defined AA levels. In addition, regional wall shear stress (WSS, mN/m2), circumferential WSS (WSSc) and axial WSS (WSSa) at all levels were quantified in 6 sectors of the aortic circle. Averaged WSS values on each level (WSSavg, WSSc, avg and WSSa, avg) were calculated as well.

RESULTS

Peak velocity at the proximal tubular AA was significantly lower in BAV and CoA group (P = 0.047) compared to BAV without CoA. In addition, the WSSa, avg was found to be higher for the BAV and CoA group at proximal AA respectively (P = 0.040). No other significant differences were found between these groups. BAV group's peak velocity was higher at every level (P < 0.001-0.004) compared to TAV group. Flow displacement was significantly higher for the BAV group at every level (P < 0.001) besides at the most distal level. All averaged WSS values were significantly higher in BAV patients in distal AA (P < 0.001-0.018).

CONCLUSIONS

Repaired CoA does not relevantly alter four-dimensional flow metrics in the AA of patients with RL-BAV. However, RL-BAV majorly alters flow dynamics in the AA when compared to patients with TAV.

CLINICAL TRIAL REGISTRATION NUMBER

https://www.clinicaltrials.gov/study/NCT05065996, Unique Protocol ID 5063566.

Toward accurate and fast velocity quantification with 3D ultrashort TE phase-contrast imaging

PURPOSE

Traditional phase-contrast MRI is affected by displacement artifacts caused by non-synchronized spatial- and velocity-encoding time points. The resulting inaccurate velocity maps can affect the accuracy of derived hemodynamic parameters. This study proposes and characterizes a 3D radial phase-contrast UTE (PC-UTE) sequence to reduce displacement artifacts. Furthermore, it investigates the displacement of a standard Cartesian flow sequence by utilizing a displacement-free synchronized-single-point-imaging MR sequence (SYNC-SPI) that requires clinically prohibitively long acquisition times.

METHODS

3D flow data was acquired at 3T at three different constant flow rates and varying spatial resolutions in a stenotic aorta phantom using the proposed PC-UTE, a Cartesian flow sequence, and a SYNC-SPI sequence as reference. Expected displacement artifacts were calculated from gradient timing waveforms and compared to displacement values measured in the in vitro flow experiments.

RESULTS

The PC-UTE sequence reduces displacement and intravoxel dephasing, leading to decreased geometric distortions and signal cancellations in magnitude images, and more spatially accurate velocity quantification compared to the Cartesian flow acquisitions; errors increase with velocity and higher spatial resolution.

CONCLUSION

PC-UTE MRI can measure velocity vector fields with greater accuracy than Cartesian acquisitions (although pulsatile fields were not studied) and shorter scan times than SYNC-SPI. As such, this approach is superior to traditional Cartesian 3D and 4D flow MRI when spatial misrepresentations cannot be tolerated, for example, when computational fluid dynamics simulations are compared to or combined with in vitro or in vivo measurements, or regional parameters such as wall shear stress are of interest.

Tetralogy of Fallot regurgitation energetics and kinetics: an intracardiac flow analysis of the right ventricle using computational fluid dynamics

Repaired Tetralogy of Fallot (rTOF) patients suffer from pulmonary regurgitation and may require pulmonary valve replacement (PVR). Cardiac magnetic resonance imaging (cMRI) guides therapy, but conventional measurements do not quantify the intracardiac flow effects from pulmonary regurgitation or PVR. This study investigates intracardiac flow parameters of the right ventricle (RV) of rTOF by computational fluid dynamics (CFD). cMRI of rTOF patients and controls were retrospectively included. Feature-tracking captured RV endocardial contours from long-axis/short-axis cine. Ventricular motion was reconstructed via diffeomorphic mapping, serving as domain boundary for CFD simulations. Vorticity (1/s), viscous energy loss (ELoss, mJ/L) and turbulent kinetic energy (TKE, mJ/L) were quantified in RV outflow tract (RVOT) and RV inflow. These parameters were normalized against total RV kinetic energy (KE) and RV inflow vorticity to derive dimensionless metrics. Vorticity contours by Q-criterion were qualitatively compared. rTOF patients (n = 15) had mean regurgitant fraction 38 ± 12% and RV size 162 ± 35 mL/m2. Compared to controls (n = 12), rTOF had increased RVOT vorticity (142.6 ± 75.6/s vs. 40.4 ± 11.8/s, p < 0.0001), Eloss (55.6 ± 42.5 vs. 5.2 ± 4.4 mJ/L, p = 0.0004), and TKE (5.7 ± 5.9 vs. 0.84 ± 0.46 mJ/L, p = 0.0003). After PVR, there was decrease in normalized RVOT Eloss/TKE (p = 0.0009, p = 0.029) and increase in normalized tricuspid inflow vorticity/KE (p = 0.0136, p = 0.043), corresponding to reorganization of the "donut"-shaped tricuspid ring-vortex. The intracardiac flow in rTOF patients can be simulated to determine the impact of PVR and improve the clinical indications guided by cardiac imaging.

Fully Automated Valve Segmentation for Blood Flow Assessment From 4D Flow MRI Including Automated Cardiac Valve Tracking and Transvalvular Velocity Mapping

BACKGROUND

Automated 4D flow MRI valvular flow quantification without time-consuming manual segmentation might improve workflow.

PURPOSE

Compare automated valve segmentation (AS) to manual (MS), and manually corrected automated segmentation (AMS), in corrected atrioventricular septum defect (c-AVSD) patients and healthy volunteers, for assessing net forward volume (NFV) and regurgitation fraction (RF).

STUDY TYPE

Retrospective.

POPULATION

27 c-AVSD patients (median, 23 years; interquartile range, 16-31 years) and 24 healthy volunteers (25 years; 12.5-36.5 years).

FIELD STRENGTH/SEQUENCE

Whole-heart 4D flow MRI and cine steady-state free precession at 3T.

ASSESSMENT

After automatic valve tracking, valve annuli were segmented on time-resolved reformatted trans-valvular velocity images by AS, MS, and AMS. NFV was calculated for all valves, and RF for right and left atrioventricular valves (RAVV and LAVV). NFV variation (standard deviation divided by mean NFV) and NFV differences (NFV difference of a valve vs. mean NFV of other valves) expressed internal NFV consistency.

STATISTICAL TESTS

Comparisons between methods were assessed by Wilcoxon signed-rank tests, and intra/interobserver variability by intraclass correlation coefficients (ICCs). P < 0.05 was considered statistically significant, with multiple testing correction.

RESULTS

AMS mean analysis time was significantly shorter compared with MS (5.3 ± 1.6 minutes vs. 9.1 ± 2.5 minutes). MS NFV variation (6.0%) was significantly smaller compared with AMS (6.3%), and AS (8.2%). Median NFV difference of RAVV, LAVV, PV, and AoV between segmentation methods ranged from -0.7-1.0 mL, -0.5-2.8 mL, -1.1-3.6 mL, and - 3.1--2.1 mL, respectively. Median RAVV and LAVV RF, between 7.1%-7.5% and 3.8%-4.3%, respectively, were not significantly different between methods. Intraobserver/interobserver agreement for AMS and MS was strong-to-excellent for NFV and RF (ICC ≥0.88).

DATA CONCLUSION

MS demonstrates strongest internal consistency, followed closely by AMS, and AS. Automated segmentation, with or without manual correction, can be considered for 4D flow MRI valvular flow quantification.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 3.

Association of Impaired Left Ventricular Mitral Filling from 4D Flow Cardiac MRI and Prognosis of Hypertrophic Cardiomyopathy

Purpose To investigate whether the peak early filling rate normalized to the filling volume (PEFR/FV) estimated from four-dimensional (4D) flow cardiac MRI may be used to assess impaired left ventricular (LV) filling and predict clinical outcomes in individuals with hypertrophic cardiomyopathy (HCM). Materials and Methods Cardiac MRI with a 4D flow sequence and late gadolinium enhancement (LGE), as well as echocardiography, was performed in 88 individuals: 44 participants with HCM from a French prospective registry (ClinicalTrials.gov; NCT01091480) and 44 healthy volunteers matched for age and sex. In participants with HCM, a composite primary end point was assessed at follow-up, including unexplained syncope, new-onset atrial fibrillation, hospitalization for congestive heart failure, ischemic stroke, sustained ventricular arrhythmia, septal reduction therapy, and cardiac death. A Cox proportional hazard model was used to analyze associations with the primary end point. Results PEFR/FV was significantly lower in the HCM group (mean age, 51.8 years ± 18.5 [SD]; 29 male participants) compared with healthy volunteers (mean, 3.35 sec-1 ± 0.99 [0.90-5.20] vs 4.42 sec-1 ± 1.68 [2.74-11.86]; P < .001) and correlated with both B-type natriuretic peptide (BNP) level (r = -0.31; P < .001) and the ratio of pulsed Doppler early transmitral inflow to Doppler tissue imaging annulus velocities (E/E'; r = -0.54; P < .001). At a median follow-up of 2.3 years (IQR, 1.7-3.3 years), the primary end point occurred in 14 (32%) participants. A PEFR/FV of 2.61 sec-1 or less was significantly associated with occurrence of the primary end point (hazard ratio, 9.46 [95% CI: 2.61, 45.17; P < .001] to 15.21 [95% CI: 3.51, 80.22; P < .001]), independently of age, BNP level, E/E', LGE extent, and LV and left atrial strain according to successive bivariate models. Conclusion In HCM, LV filling evaluated with 4D flow cardiac MRI correlated with Doppler and biologic indexes of diastolic dysfunction and predicted clinical outcomes. Keywords: Diastolic Function, Left Ventricular Filling, Hypertrophic Cardiomyopathy, Cardiac MRI, 4D Flow Sequence Clinical trial registration no. NCT01091480 Supplemental material is available for this article. © RSNA, 2024.

Ascending Aorta 4D Time to Peak Distention Sexual Dimorphism and Association with Coronary Plaque Burden Severity in Women

Coronary artery disease (CAD) risk and plaque scores are often subjective and biased, particularly in mid-age asymptomatic women, whose CAD risk assessment has been historically underestimated. In this study, a new automatic ascending aorta time-to-peak-distention (TPD) analysis was developed for fast screening and as an independent surrogate for subclinical atherosclerosis in asymptomatic women. CCTA was obtained in 50 asymptomatic adults. Plaque burden segment involvement score (SIS) and automatic TPD were obtained from all subjects. Logistic regression analyses were performed to investigate the association between CAD risk scores and TPD with severe coronary plaque burden (SIS>5). TPD, individually, was found to be a significant predictor of SIS>5. Additionally, sex was a significant effect modifier of TPD, with a stronger statistically significant association with women. Four-dimensional aortic time-to-peak distention could supplement conventional CCTA analysis and offer a quick objective screening tool for plaque burden severity and CAD risk stratification, especially in women.

Diastolic function assessment with four-dimensional flow cardiovascular magnetic resonance using automatic deep learning E/A ratio analysis

BACKGROUND

Diastolic left ventricular (LV) dysfunction is a powerful contributor to the symptoms and prognosis of patients with heart failure. In patients with depressed LV systolic function, the E/A ratio, the ratio between the peak early (E) and the peak late (A) transmitral flow velocity, is the first step to defining the grade of diastolic dysfunction. Doppler echocardiography (echo) is the preferred imaging technique for diastolic function assessment, while cardiovascular magnetic resonance (CMR) is less established as a method. Previous four-dimensional (4D) Flow-based studies have looked at the E/A ratio proximal to the mitral valve, requiring manual interaction. In this study, we compare an automated, deep learning-based and two semi-automated approaches for 4D Flow CMR-based E/A ratio assessment to conventional, gold-standard echo-based methods.

METHODS

Ninety-seven subjects with chronic ischemic heart disease underwent a cardiac echo followed by CMR investigation. 4D Flow-based E/A ratio values were computed using three different approaches; two semi-automated, assessing the E/A ratio by measuring the inflow velocity (MVvel) and the inflow volume (MVflow) at the mitral valve plane, and one fully automated, creating a full LV segmentation using a deep learning-based method with which the E/A ratio could be assessed without constraint to the mitral plane (LVvel).

RESULTS

MVvel, MVflow, and LVvel E/A ratios were strongly associated with echocardiographically derived E/A ratio (R2 = 0.60, 0.58, 0.72). LVvel peak E and A showed moderate association to Echo peak E and A, while MVvel values were weakly associated. MVvel and MVflow EA ratios were very strongly associated with LVvel (R2 = 0.84, 0.86). MVvel peak E was moderately associated with LVvel, while peak A showed a strong association (R2 = 0.26, 0.57).

CONCLUSION

Peak E, peak A, and E/A ratio are integral to the assessment of diastolic dysfunction and may expand the utility of CMR studies in patients with cardiovascular disease. While underestimation of absolute peak E and A velocities was noted, the E/A ratio measured with all three 4D Flow methods was strongly associated with the gold standard Doppler echocardiography. The automatic, deep learning-based method performed best, with the most favorable runtime of ∼40 seconds. As both semi-automatic methods associated very strongly to LVvel, they could be employed as an alternative for estimation of E/A ratio.

Comparison of Echo-Planar Imaging and Compressed Sensing in the Estimation of Flow Metrics from Aortic 4D Flow MR Imaging: A Healthy Volunteer Study

PURPOSE

Prolonged scanning of time-resolved 3D phase-contrast MRI (4D flow MRI) limits its routine use in clinical practice. An echo-planar imaging (EPI)-based sequence and compressed sensing can reduce the scan duration. We aimed to determine the impact of EPI for 4D flow MRI on the scan duration, image quality, and quantitative flow metrics.

METHODS

This was a prospective study of 15 healthy volunteers (all male, mean age 33 ± 5 years). Conventional sensitivity encoding (SENSE), EPI with SENSE (EPI), and compressed SENSE (CS) (reduction factors: 6 and 12, respectively) were scanned.Scan duration, qualitative indexes of image quality, and quantitative flow parameters of net flow volume, maximum flow velocity, wall shear stress (WSS), and energy loss (EL) in the ascending aorta were assessed. Two-dimensional phase-contrast cine MRI (2D-PC) was considered the gold standard of net flow volume and maximum flow velocity.

RESULTS

Compared to SENSE, EPI and CS12 shortened scan durations by 71% and 73% (EPI, 4 min 39 sec; CS6, 7 min 29 sec; CS12, 4 min 14 sec; and SENSE, 15 min 51 sec). Visual image quality was significantly better for EPI than for SENSE and CS (P < 0.001). The net flow volumes obtained with SENSE, EPI, and CS12 and those obtained with 2D-PC were correlated well (r = 0.950, 0.871, and 0.850, respectively). However, the maximum velocity obtained with EPI was significantly underestimated (P < 0.010). The average WSS was significantly higher with EPI than with SENSE, CS6, and CS12 (P < 0.001, P = 0.040, and P = 0.012, respectively). The EL was significantly lower with EPI than with CS6 and CS12 (P = 0.002 and P = 0.007, respectively).

CONCLUSION

EPI reduced the scan duration, improved visual image quality, and was associated with more accurate net flow volume than CS. However, the flow velocity, WSS, and EL values obtained with EPI and other sequences may not be directly comparable.

Evaluating a Phase-Specific Approach to Aortic Flow: A 4D Flow MRI Study

BACKGROUND

Aortic flow parameters can be quantified using 4D flow MRI. However, data are sparse on how different methods of analysis influence these parameters and how these parameters evolve during systole.

PURPOSE

To assess multiphase segmentations and multiphase quantification of flow-related parameters in aortic 4D flow MRI.

STUDY TYPE

Prospective.

POPULATION

40 healthy volunteers (50% male, 28.9 ± 5.0 years) and 10 patients with thoracic aortic aneurysm (80% male, 54 ± 8 years).

FIELD STRENGTH/SEQUENCE

4D flow MRI with a velocity encoded turbo field echo sequence at 3 T.

ASSESSMENT

Phase-specific segmentations were obtained for the aortic root and the ascending aorta. The whole aorta was segmented in peak systole. In all aortic segments, time to peak (TTP; for flow velocity, vorticity, helicity, kinetic energy, and viscous energy loss) and peak and time-averaged values (for velocity and vorticity) were calculated.

STATISTICAL TESTS

Static vs. phase-specific models were assessed using Bland-Altman plots. Other analyses were performed using phase-specific segmentations for aortic root and ascending aorta. TTP for all parameters was compared to TTP of flow rate using paired t-tests. Time-averaged and peak values were assessed using Pearson correlation coefficient. P < 0.05 was considered statistically significant.

RESULTS

In the combined group, velocity in static vs. phase-specific segmentations differed by 0.8 cm/sec for the aortic root, and 0.1 cm/sec (P = 0.214) for the ascending aorta. Vorticity differed by 167 sec-1  mL-1 (P = 0.468) for the aortic root, and by 59 sec-1  mL-1 (P = 0.481) for the ascending aorta. Vorticity, helicity, and energy loss in the ascending aorta, aortic arch, and descending aorta peaked significantly later than flow rate. Time-averaged velocity and vorticity values correlated significantly in all segments.

DATA CONCLUSION

Static 4D flow MRI segmentation yields comparable results as multiphase segmentation for flow-related parameters, eliminating the need for time-consuming multiple segmentations. However, multiphase quantification is necessary for assessing peak values of aortic flow-related parameters.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 3.

Quantitative 4D flow MRI-derived thoracic aortic normal values of 2D flow MRI parameters in healthy volunteers

PURPOSE

 To utilize 4 D flow MRI to acquire normal values of "conventional 2 D flow MRI parameters" in healthy volunteers in order to replace multiple single 2 D flow measurements with a single 4 D flow acquisition.

MATERIALS AND METHODS

 A kt-GRAPPA accelerated 4 D flow sequence was used. Flow volumes were assessed by forward (FFV), backward (BFV), and net flow volumes (NFV) [ml/heartbeat] and flow velocities by axial (VAX) and absolute velocity (VABS) [m/s] in 116 volunteers (58 females, 43 ± 13 years). The aortic regurgitant fraction (RF) was calculated.

RESULTS

 The sex-neutral mean FFV, BFV, NFV, and RF in the ascending aorta were 93.5 ± 14.8, 3.6 ± 2.8, 89.9 ± 0.6 ml/heartbeat, and 3.9 ± 2.9 %, respectively. Significantly higher values were seen in males regarding FFV, BFV, NFV and RF, but there was no sex dependency regarding VAX and VABS. The mean maximum VAX was lower (1.01 ± 0.31 m/s) than VABS (1.23 ± 0.35 m/s). We were able to determine normal ranges for all intended parameters.

CONCLUSION

 This study provides quantitative 4 D flow-derived thoracic aortic normal values of 2 D flow parameters in healthy volunteers. FFV, BFV, NFV, and VAX did not differ significantly from single 2 D flow acquisitions and could therefore replace time-consuming multiple single 2 D flow acquisitions. VABS should not be used interchangeably.

KEY POINTS

  · 4 D flow MRI can be used to replace 2 D flow MRI measurements.. · The parameter absolute velocities can be assessed by 4 D flow MRI.. · There are sex-dependent differences regarding forward, backward, net aortic blood flow and the aortic valve regurgitant fraction..

Four-dimensional flow cardiovascular magnetic resonance aortic cross-sectional pressure changes and their associations with flow patterns in health and ascending thoracic aortic aneurysm

BACKGROUND

Ascending thoracic aortic aneurysm (ATAA) is a silent and threatening dilation of the ascending aorta (AscAo). Maximal aortic diameter which is currently used for ATAA patients management and surgery planning has been shown to inadequately characterize risk of dissection in a large proportion of patients. Our aim was to propose a comprehensive quantitative evaluation of aortic morphology and pressure-flow-wall associations from four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) data in healthy aging and in patients with ATAA.

METHODS

We studied 17 ATAA patients (64.7 ± 14.3 years, 5 females) along with 17 age- and sex-matched healthy controls (59.7 ± 13.3 years, 5 females) and 13 younger healthy subjects (33.5 ± 11.1 years, 4 females). All subjects underwent a CMR exam, including 4D flow and three-dimensional anatomical images of the aorta. This latter dataset was used for aortic morphology measurements, including AscAo maximal diameter (iDMAX) and volume, indexed to body surface area. 4D flow MRI data were used to estimate 1) cross-sectional local AscAo spatial (∆PS) and temporal (∆PT) pressure changes as well as the distance (∆DPS) and time duration (∆TPT) between local pressure peaks, 2) AscAo maximal wall shear stress (WSSMAX) at peak systole, and 3) AscAo flow vorticity amplitude (VMAX), duration (VFWHM), and eccentricity (VECC).

RESULTS

Consistency of flow and pressure indices was demonstrated through their significant associations with AscAo iDMAX (WSSMAX:r = -0.49, p < 0.001; VECC:r = -0.29, p = 0.045; VFWHM:r = 0.48, p < 0.001; ∆DPS:r = 0.37, p = 0.010; ∆TPT:r = -0.52, p < 0.001) and indexed volume (WSSMAX:r = -0.63, VECC:r = -0.51, VFWHM:r = 0.53, ∆DPS:r = 0.54, ∆TPT:r = -0.63, p < 0.001 for all). Intra-AscAo cross-sectional pressure difference, ∆PS, was significantly and positively associated with both VMAX (r = 0.55, p = 0.002) and WSSMAX (r = 0.59, p < 0.001) in the 30 healthy subjects (48.3 ± 18.0 years). Associations remained significant after adjustment for iDMAX, age, and systolic blood pressure. Superimposition of ATAA patients to normal aging trends between ∆PS and WSSMAX as well as VMAX allowed identifying patients with substantially high pressure differences concomitant with AscAo dilation.

CONCLUSION

Local variations in pressures within ascending aortic cross-sections derived from 4D flow MRI were associated with flow changes, as quantified by vorticity, and with stress exerted by blood on the aortic wall, as quantified by wall shear stress. Such flow-wall and pressure interactions might help for the identification of at-risk patients.

Evaluation of left ventricular and left atrial volumetric function from native MR multislice 4D flow magnitude data

OBJECTIVES

To assess the feasibility, precision, and accuracy of left ventricular (LV) and left atrial (LA) volumetric function evaluation from native magnetic resonance (MR) multislice 4D flow magnitude images.

MATERIALS & METHODS

In this prospective study, 60 subjects without signs or symptoms of heart failure underwent 3T native cardiac MR multislice 4D flow and bSSFP-cine realtime imaging. LV and LA volumetric function parameters were evaluated from 4D flow magnitude (4D flow-cine) and bSSFP-cine data using standard software to obtain end-diastolic volume (EDV), end-systolic volume (ESV), ejection-fraction (EF), stroke-volume (SV), LV muscle mass (LVM), LA maximum volume, LA minimum volume, and LA total ejection fraction (LATEF). Stroke volumes derived from both imaging methods were further compared to 4D pulmonary artery flow-derived net forward volumes (NFV). Methods were compared by correlation and Bland-Altman analysis.

RESULTS

Volumetric function parameters from 4D flow-cine and bSSFP-cine showed high to very high correlations (r = 0.83-0.98). SV, LA volumes and LATEF did not differ between methods. LV end-diastolic and end-systolic volumes were slightly underestimated (EDV: -2.9 ± 5.8 mL; ESV: -2.3 ± 3.8 mL), EF was slightly overestimated (EF: 0.9 ± 2.6%), and LV mass was considerably overestimated (LVM: 39.0 ± 11.4 g) by 4D flow-cine imaging. SVs from both methods correlated very highly with NFV (r = 0.91 in both cases) and did not differ from NFV.

CONCLUSION

Native multislice 4D flow magnitude data allows precise evaluation of LV and LA volumetric parameters; however, apart from SV, LV volumetric parameters demonstrate bias and need to be referred to their respective normal values.

CLINICAL RELEVANCE STATEMENT

Volumetric function assessment from native multislice 4D flow magnitude images can be performed with routinely used clinical software, facilitating the application of 4D flow as a one-stop-shop functional cardiac MR exam, providing consistent, simultaneously acquired, volume and flow data.

KEY POINTS

• Native multislice 4D flow imaging allows evaluation of volumetric left ventricular and atrial function parameters. • Left ventricular and left atrial function parameters derived from native multislice 4D flow data correlate highly with corresponding standard cine-derived parameters. • Multislice 4D flow-derived volumetric stroke volume and net forward volume do not differ.

Characterization of pulmonary arterial stiffness using cardiac MRI

Pulmonary arterial stiffness (PAS) is a pathologic hallmark of all types of pulmonary hypertension (PH). Cardiac MRI (CMR), a gold-standard imaging modality for the evaluation of pulmonary flow, biventricular morphology and function has been historically reserved for the longitudinal clinical follow-up, PH phenotyping purposes, right ventricular evaluation, and research purposes. Over the last two decades, numerous indices combining invasive catheterization and non-invasive CMR have been utilized to phenotype the character and severity of PAS in different types of PH and to assess its clinically prognostic potential with encouraging results. Many recent studies have demonstrated a strong role of CMR derived PAS markers in predicting long-term clinical outcomes and improving currently gold standard risk assessment provided by the REVEAL calculator. With the utilization of a machine learning strategies, strong diagnostic and prognostic performance of CMR reported in multicenter studies, and ability to detect PH at early stages, the non-invasive assessment of PAS is on verge of routine clinical utilization. In this review, we focus on appraising important CMR studies interrogating PAS over the last 20 years, describing the benefits and limitations of different PAS indices, and their pathophysiologic relevance to pulmonary vascular remodeling. We also discuss the role of CMR and PAS in clinical surveillance and phenotyping of PH, and the long-term future goal to utilize PAS as a biomarker to aid with more targeted therapeutic management.

Revolutionizing vascular imaging: trends and future directions of 4D flow MRI based on a 20-year bibliometric analysis

Background

Four-dimensional flow magnetic resonance imaging (4D flow MRI) is a promising new technology with potential clinical value in hemodynamic quantification. Although an increasing number of articles on 4D flow MRI have been published over the past decades, few studies have statistically analyzed these published articles. In this study, we aimed to perform a systematic and comprehensive bibliometric analysis of 4D flow MRI to explore the current hotspots and potential future directions.

Methods

The Web of Science Core Collection searched for literature on 4D flow MRI between 2003 and 2022. CiteSpace was utilized to analyze the literature data, including co-citation, cooperative network, cluster, and burst keyword analysis.

Results

A total of 1,069 articles were extracted for this study. The main research hotspots included the following: quantification and visualization of blood flow in different clinical settings, with keywords such as "cerebral aneurysm", "heart", "great vessel", "tetralogy of Fallot", "portal hypertension", and "stiffness"; optimization of image acquisition schemes, such as "resolution" and "reconstruction"; measurement and analysis of flow components and patterns, as indicated by keywords "pattern", "KE", "WSS", and "fluid dynamics". In addition, international consensus for metrics derived from 4D flow MRI and multimodality imaging may also be the future research direction.

Conclusions

The global domain of 4D flow MRI has grown over the last 2 decades. In the future, 4D flow MRI will evolve towards becoming a relatively short scan duration with adequate spatiotemporal resolution, expansion into the diagnosis and treatment of vascular disease in other related organs, and a shift in focus from vascular structure to function. In addition, artificial intelligence (AI) will assist in the clinical promotion and application of 4D flow MRI.

Robust and fast stochastic 4D flow vector-field signature technique for quantifying composite flow dynamics from 4D flow MRI: Application to left atrial flow in atrial fibrillation

4D flow MRI is an emerging imaging modality that maps voxel-wise blood flow information as velocity vector fields that is acquired in 7-dimensional image volumes (3 spatial dimensions + 3 velocity directions + time). Blood flow in the cardiovascular system is often complex and composite involving multiple flow dynamics and patterns (e.g., vortex flow, jets, stagnating flow) that occur and interact simultaneously. The spectrum of such complex flow dynamics is embedded in the velocity vector field dynamics derived from 4D Flow MRI. However, current flow metrics cannot fully measure high-dimensional vector-field data and embedded complex composite flow data. Instead, these methods need to break down the vector-field data into secondary scalar fields of individual flow components using fluid dynamics operators. These methods are gradient-based and sensitive to data uncertainties, and only focus on individual flow components of the overall composite flow, therefore potentially underestimating the severity of overall flow changes associated with cardiovascular diseases. To address these limitations, in MICCAI 2021, we introduced a novel comprehensive stochastic 4D Flow vector-field signature technique that works directly on the entire spatiotemporal velocity vector field. This technique uses efficient stochastic gradient-free interrogation of multi-million flow vector-pairs per patient to derive the patient's unique flow profile of the complex composite flow alterations and in real-time processing. The signature technique's probabilistic gradient-free formulation should allow for highly robust quantification despite inherent errors in 4D flow MRI acquisitions. Here, we extend the application of the 4D flow vector-field signature technique to the left atrium to analyze complex composite flow changes in patients with atrial fibrillation. In 128 subjects, we performed extensive sensitivity testing and determined that the vector-field signature technique is highly robust to typical sources of data uncertainties in 4D flow MRI: degradation in spatiotemporal resolution, added Gaussian noise, and segmentation errors. We demonstrate the excellent generalizability of the stochastic convergence from the aorta to the left atrium and between different 4D Flow MRI acquisition protocols. We compare the robustness of our technique to existing advanced flow quantification metrics of kinetic energy, vorticity, and energy loss demonstrating a superior performance of up-to 14-fold. Our results show the potential diagnostic and clinical utility of our signature technique in identifying distinctly altered composite flow signatures in atrial fibrillation patients independent of existing flow metrics.

Impact of the evaluation method on 4D flow-derived diastolic transmitral and myocardial peak velocities: Comparison with echocardiography

PURPOSE

To compare agreement of different evaluation methods of magnetic resonance (MR) 4D flow-derived diastolic transmitral and myocardial peak velocities as well as their ratios, using echocardiography as reference.

METHODS

In this prospective study, 60 subjects without symptoms of cardiovascular disease underwent echocardiography and non-contrast 3 T MR 4D flow imaging of the heart. Early- (E) and late-diastolic (A) transmitral peak filling velocities were evaluated from 4D flow data using three different strategies: 1) at the mitral valve tips in short-axis orientation (SA-method), 2) between the mitral valve tips in 4-chamber orientation (4-chamber-method), and 3) as maximal velocities in the transmitral inflow volume (max-velocity-method). Septal, lateral and average early-diastolic myocardial peak velocities (e') were derived from the myocardial tissue in the vicinity of the mitral valve. 4D flow parameters were compared with echocardiography by correlation and Bland-Altman analysis.

RESULTS

All 4D flow-derived E, A and E/A values correlated with echocardiography (r = 0.65-0.73, 0.75-0.83 and 0.74-0.86, respectively). While the SA- and 4-chamber-methods substantially underestimated E and A compared to echocardiography (p < 0.001), the max-velocity-method provided E (p = 0.13) and E/A (p = 0.07) without significant bias. Septal, lateral and average e' from 4D flow as well as the max-velocity-method-derived E/e' correlated with echocardiographic measurements (r = 0.64-0.81) and showed no significant bias (p = 0.26-0.54).

CONCLUSION

MR 4D flow imaging allows precise and accurate evaluation of transmitral and myocardial peak velocities for characterization of LV diastolic function without significant bias to echocardiography, when transmitral velocities are assessed from the transmitral inflow volume. This enables the use of validated echocardiography threshold values.

Using Gaussian process for velocity reconstruction after coronary stenosis applicable in positron emission particle tracking: An in-silico study

Accurate velocity reconstruction is essential for assessing coronary artery disease. We propose a Gaussian process method to reconstruct the velocity profile using the sparse data of the positron emission particle tracking (PEPT) in a biological environment, which allows the measurement of tracer particle velocity to infer fluid velocity fields. We investigated the influence of tracer particle quantity and detection time interval on flow reconstruction accuracy. Three models were used to represent different levels of stenosis and anatomical complexity: a narrowed straight tube, an idealized coronary bifurcation with stenosis, and patient-specific coronary arteries with a stenotic left circumflex artery. Computational fluid dynamics (CFD), particle tracking, and the Gaussian process of kriging were employed to simulate and reconstruct the pulsatile flow field. The study examined the error and uncertainty in velocity profile reconstruction after stenosis by comparing particle-derived flow velocity with the CFD solution. Using 600 particles (15 batches of 40 particles) released in the main coronary artery, the time-averaged error in velocity reconstruction ranged from 13.4% (no occlusion) to 161% (70% occlusion) in patient-specific anatomy. The error in maximum cross-sectional velocity at peak flow was consistently below 10% in all cases. PEPT and kriging tended to overestimate area-averaged velocity in higher occlusion cases but accurately predicted maximum cross-sectional velocity, particularly at peak flow. Kriging was shown to be useful to estimate the maximum velocity after the stenosis in the absence of negative near-wall velocity.

Comparison of aortic blood flow rotational direction in healthy volunteers and patients with bicuspid aortic valves using volumetric velocity-sensitive cardiovascular magnetic resonance imaging

Background

The rotational direction (RD) of helical blood flow can be classified as either a clockwise (RD+) or counter-clockwise (RD-) flow. We hypothesized that this simple classification might not be sufficient for analysis in vivo and a simultaneous existence of RD+/- may occur. We utilized volumetric velocity-sensitive cardiovascular magnetic resonance imaging (4D flow MRI) to analyze rotational blood flow in the thoracic aorta.

Methods

Forty volunteers (22 females; mean age, 41±16 years) and seventeen patients with bicuspid aortic valves (BAVs) (9 females; mean age, 42±14 years) were prospectively included. The RDs and the calculation of the rotating blood volumes (RBVs) in the thoracic aorta were performed using a pathline-projection strategy.

Results

We could confirm a mainly clockwise RD in the ascending, descending aorta and in the aortic arch. Furthermore, we found a simultaneous existence of RD+/RD-. The RD+/--volume in the ascending aorta was significantly higher in BAV patients, the mean RD+/RD- percentage was approximately 80%/20% vs. 60%/40% in volunteers (P<0.01). The maximum RBV always occurred during systole. There was significantly more clockwise than counter-clockwise rotational flow in the ascending aorta (P<0.01) and the aortic arch (P<0.01), but no significant differences in the descending aorta (P=0.48).

Conclusions

A simultaneous occurrence of RD+/RD- indicates that a simple categorization in either of both is insufficient to describe blood flow in vivo. Rotational flow in the ascending aorta and in the aortic arch differs significantly from flow in the descending aorta. BAV patients show significantly more clockwise rotating volume in the ascending aorta compared to healthy volunteers.

Contrast Agent-Free 3D Renal Ultrafast Doppler Imaging Reveals Vascular Dysfunction in Acute and Diabetic Kidney Diseases

To combat the irreversible decline in renal function associated with kidney disease, it is essential to establish non-invasive biomarkers for assessing renal microcirculation. However, the limited resolution and/or vascular sensitivity of existing diagnostic imaging techniques hinders the visualization of complex cortical vessels. Here, a 3D renal ultrafast Doppler (UFD) imaging system that uses a high ultrasound frequency (18 MHz) and ultrahigh frame rate (1 KHz per slice) to scan the entire volume of a rat's kidney in vivo is demonstrated. The system, which can visualize the full 3D renal vascular branching pyramid at a resolution of 167 µm without any contrast agent, is used to chronically and noninvasively monitor kidneys with acute kidney injury (AKI, 3 days) and diabetic kidney disease (DKD, 8 weeks). Multiparametric UFD analyses (e.g., vessel volume occupancy (VVO), fractional moving blood volume (FMBV), vessel number density (VND), and vessel tortuosity (VT)) describe rapid vascular rarefaction from AKI and long-term vascular degeneration from DKD, while the renal pathogeneses are validated by in vitro blood serum testing and stained histopathology. This work demonstrates the potential of 3D renal UFD to offer valuable insights into assessing kidney perfusion levels for future research in diabetes and kidney transplantation.

Hemi-Fontan and bidirectional Glenn operations result in flow-mediated viscous energy loss at the time of stage II palliation

Background

Superior cavopulmonary connection (SCPC) for stage II palliation of hypoplastic left heart syndrome (HLHS) is achieved most frequently by either a bidirectional Glenn (BDG) or hemi-Fontan (HF) operation. The comparison of flow hemodynamic efficiency at the region of surgical reconstruction and in proximal pulmonary arteries has been evaluated primarily using computational modeling techniques with conflicting reports. The purpose of this descriptive study was to compare flow hemodynamics following stage II (BDG vs HF) using 4-dimensional flow magnetic resonance imaging (4D-Flow MRI) with particular focus on flow-mediated viscous energy loss (EL') under matched hemodynamic conditions.

Methods

Patients with hypoplastic left heart syndrome (HLHS) who underwent either HF or BDG as part of stage II palliation underwent pre-Fontan 4D-Flow MRI. Patients were matched by the pulmonary vascular resistance index, net superior vena cava (SVC) flow, right pulmonary artery (RPA) and left pulmonary artery (LPA) size, and age. Maximum EL' throughout the cardiac cycle was calculated along the SVC-RPA and SVC-LPA tracts.

Results

Eight patients who underwent HF as part of their stage II single ventricle palliation were matched with 8 patients who underwent BDG. There were no differences between the 2 groups in median volumetric indices, including end-diastolic volume (P = .278) and end-systolic volume (P = .213). Moreover, no differences were observed in ejection fraction (P = .091) and cardiac index (P = .324). There also were no differences in peak EL' measured along the SVC-RPA tract (median, 0.05 mW for HF vs 0.04 mW for BDG; P = .365) or along the SVC-LPA tract (median, 0.05 mW vs 0.04 mW; P = .741).

Conclusions

The second stage of surgical palliation of HLHS using either HF or BDG results in similar flow-mediated viscous energy loss throughout the SCPC junction. 4D-Flow MRI and computational methods should be applied together to investigate flow hemodynamic patterns throughout the Fontan palliation and overall efficiency of the Fontan circuit.

Advanced Imaging Technologies for Assessing Tetralogy of Fallot: Insights Into Flow Dynamics

Tetralogy of Fallot is the most common cyanotic congenital heart defect requiring surgical repair. Although surgical interventions have significantly reduced mortality, postrepair complications, such as pulmonary valve regurgitation and stenosis, may lead to adverse outcomes, including right ventricular dysfunction and increased risks of morbidity and mortality. This review explores the potential of advanced imaging technologies, including 4-dimensional-flow magnetic resonance imaging and high-frame-rate echocardiography, in providing valuable insights into blood flow dynamics and energy parameters. Quantitative measures, such as energy loss and vorticity, along with qualitative flow analysis, can provide additional insights into adverse haemodynamics at a potentially earlier and more reversible stage. Furthermore, personalized patient-specific information from these imaging modalities aids in guiding treatment decisions and monitoring postoperative interventions effectively. By characterizing flow patterns, these advanced imaging techniques hold great promise in improving the assessment and management of tetralogy of Fallot, providing tailored insights. However, further research and longitudinal studies are required to fully establish their clinical utility and potential impact on patient care.

Highly accelerated compressed sensing 4D flow MRI in congenital and acquired heart disease: comparison of aorta and main pulmonary artery flow parameters with conventional 4D flow MRI in children and young adults

BACKGROUND

Four-dimensional flow (4D flow) MRI has become a clinically utilized cardiovascular flow assessment tool. However, scans can be lengthy and may require anesthesia in younger children. Adding compressed sensing can decrease scan time, but its impact on hemodynamic data accuracy needs additional assessment.

OBJECTIVE

To compare 4D flow hemodynamics acquired with and without compressed sensing.

MATERIALS AND METHODS

Twenty-seven patients (median age: 13 [IQR: 9.5] years) underwent conventional and compressed sensing cardiovascular 4D flow following informed consent. Conventional 4D flow was performed using parallel imaging and an acceleration factor of 2. Compressed sensing 4D flow was performed with an acceleration factor of 7.7. Regions of interest were placed to compare flow parameters in the ascending aorta and main pulmonary artery. Paired Student's t-tests, Wilcoxon signed-rank tests, Bland-Altman plots, and intraclass correlation coefficients were conducted. A P-value of < 0.05 was considered statistically significant.

RESULTS

Mean scan acquisition time was reduced by 59% using compressed sensing (3.4 vs. 8.2 min, P < 0.001). Flow quantification was similar for compressed sensing and conventional 4D flow for the ascending aorta net flow: 47 vs. 49 ml/beat (P = 0.28); forward flow: 49 vs. 50 ml/beat (P = 0.07), and main pulmonary artery net flow: 49 vs. 51 ml/beat (P = 0.18); forward flow: 50 vs. 55 ml/beat (P = 0.07). Peak systolic velocity was significantly underestimated by compressed sensing 4D flow in the ascending aorta: 114 vs. 128 cm/s (P < 0.001) and main pulmonary artery: 106 vs. 112 cm/s (P = 0.02).

CONCLUSION

For both the aorta and main pulmonary artery, compressed sensing 4D flow provided equivalent net and forward flow values compared to conventional 4D flow but underestimated peak systolic velocity. By reducing scan time, compressed sensing 4D flow may decrease the need for anesthesia and increase scanner output without significantly compromising data integrity.

Generalized super-resolution 4D Flow MRI-using ensemble learning to extend across the cardiovascular system

4D Flow Magnetic Resonance Imaging (4D Flow MRI) is a non-invasive measurement technique capable of quantifying blood flow across the cardiovascular system. While practical use is limited by spatial resolution and image noise, incorporation of trained super-resolution (SR) networks has potential to enhance image quality post-scan. However, these efforts have predominantly been restricted to narrowly defined cardiovascular domains, with limited exploration of how SR performance extends across the cardiovascular system; a task aggravated by contrasting hemodynamic conditions apparent across the cardiovasculature. The aim of our study was to explore the generalizability of SR 4D Flow MRI using a combination of heterogeneous training sets and dedicated ensemble learning. With synthetic training data generated across three disparate domains (cardiac, aortic, cerebrovascular), varying convolutional base and ensemble learners were evaluated as a function of domain and architecture, quantifying performance on both in-silico and acquired in-vivo data from the same three domains. Results show that both bagging and stacking ensembling enhance SR performance across domains, accurately predicting high-resolution velocities from low-resolution input data in-silico. Likewise, optimized networks successfully recover native resolution velocities from downsampled in-vivo data, as well as show qualitative potential in generating denoised SR-images from clinicallevel input data. In conclusion, our work presents a viable approach for generalized SR 4D Flow MRI, with ensemble learning extending utility across various clinical areas of interest.

A time-consistent stabilized finite element method for fluids with applications to hemodynamics

Several finite element methods for simulating incompressible flows rely on the streamline upwind Petrov-Galerkin stabilization (SUPG) term, which is weighted by [Formula: see text]. The conventional formulation of [Formula: see text] includes a constant that depends on the time step size, producing an overall method that becomes exceedingly less accurate as the time step size approaches zero. In practice, such method inconsistency introduces significant error in the solution, especially in cardiovascular simulations, where small time step sizes may be required to resolve multiple scales of the blood flow. To overcome this issue, we propose a consistent method that is based on a new definition of [Formula: see text]. This method, which can be easily implemented on top of an existing streamline upwind Petrov-Galerkin and pressure stabilizing Petrov-Galerkin method, involves the replacement of the time step size in [Formula: see text] with a physical time scale. This time scale is calculated in a simple operation once every time step for the entire computational domain from the ratio of the L2-norm of the acceleration and the velocity. The proposed method is compared against the conventional method using four cases: a steady pipe flow, a blood flow through vascular anatomy, an external flow over a square obstacle, and a fluid-structure interaction case involving an oscillatory flexible beam. These numerical experiments, which are performed using linear interpolation functions, show that the proposed formulation eliminates the inconsistency issue associated with the conventional formulation in all cases. While the proposed method is slightly more costly than the conventional method, it significantly reduces the error, particularly at small time step sizes. For the pipe flow where an exact solution is available, we show the conventional method can over-predict the pressure drop by a factor of three. This large error is almost completely eliminated by the proposed formulation, dropping to approximately 1% for all time step sizes and Reynolds numbers considered.

Intracardiac shunt assessment using CT coronary angiography

BACKGROUND

Detection of intracardiac shunts using CT Coronary Angiography (CTCA) is currently based on anatomical demonstration of defects. We assessed a novel technique using a standard CTCA test bolus in detecting shunts independent of anatomical assessment and to provide an estimate of Qp/Qs.

METHODS

We retrospectively reviewed 51 CTCAs: twenty-one from patients with known simple left to right intracardiac shunts with contemporaneous functional assessment (using CMR) within 6 months, 20 controls with structurally normal hearts, and 10 patients with shunt repairs. From the dynamic acquisition of a test bolus, we measured mean Hounsfield Units (HU) in various anatomical structures. We created time/density curves from the test bolus data, and calculated disappearance time (DT) from the ascending aorta (deriving a Qp/Qs), peak ascending aortic HU, and mean coefficient of variation of the arterial curves, and compared these with the Qp/Qs from the respective CMR.

RESULTS

Patients with intracardiac shunts had significantly higher test bolus derived Qp/Qs compared with both the controls, and the repaired shunt comparator group. There was a very strong agreement between the test bolus derived Qp/Qs, and Qp/Qs as measured by CMR (Intraclass correlation 0.89). Mean bias was 0.032 ​± ​0.341 (95% limits of agreement -0.64 to 0.70). Interobserver, and intraobserver agreement of the disappearance time was excellent (0.99, 0.99 (reader 1) and 1.00 (reader 2) respectively).

CONCLUSION

In this proof-of-concept study, we demonstrate a novel technique to detect, and to estimate severity of left to right intracardiac shunts on routine Cardiac CT.

4Dflow-VP-Net: A deep convolutional neural network for noninvasive estimation of relative pressures in stenotic flows from 4D flow MRI

PURPOSE

To estimate relative transvalvular pressure gradient (TVPG) noninvasively from 4D flow MRI.

METHODS

A novel deep learning-based approach is proposed to estimate pressure gradient across stenosis from four-dimensional flow MRI (4D flow MRI) velocities. A deep neural network 4D flow Velocity-to-Presure Network (4Dflow-VP-Net) was trained to learn the spatiotemporal relationship between velocities and pressure in stenotic vessels. Training data were simulated by computational fluid dynamics (CFD) for different pulsatile flow conditions under an aortic flow waveform. The network was tested to predict pressure from CFD-simulated velocity data, in vitro 4D flow MRI data, and in vivo 4D flow MRI data of patients with both moderate and severe aortic stenosis. TVPG derived from 4Dflow-VP-Net was compared to catheter-based pressure measurements for available flow rates, in vitro and Doppler echocardiography-based pressure measurement, in vivo.

RESULTS

Relative pressures calculated by 4Dflow-VP-Net and in vitro pressure catheterization revealed strong correlation (r2  = 0.91). Correlations analysis of TVPG from reference CFD and 4Dflow-VP-Net for 450 simulated flow conditions showed strong correlation (r2  = 0.99). TVPG from in vitro MRI had a correlation coefficient of r2  = 0.98 with reference CFD. 4Dflow-VP-Net, applied to 4D flow MRI in 16 patients, showed comparable TVPG measurement with Doppler echocardiography (r2  = 0.85). Bland-Altman analysis of TVPG measurements showed mean bias and limits of agreement of -0.20 ± 2.07 mmHg and 0.19 ± 0.45 mmHg for CFD-simulated velocities and in vitro 4D flow velocities. In patients, overestimation of Doppler echocardiography relative to TVPG from 4Dflow-VP-Net (10.99 ± 6.77 mmHg) was observed.

CONCLUSION

The proposed approach can predict relative pressure in both in vitro and in vivo 4D flow MRI of aortic stenotic patients with high fidelity.

Obesity dominates early effects on cardiac structure and arterial stiffness in people with type 2 diabetes

OBJECTIVE

Type 2 diabetes (T2D) and obesity are global epidemics leading to excess cardiovascular disease (CVD). This study investigates standard and novel cardiac MRI parameters to detect subclinical cardiac and central vascular dysfunction in inactive people with and without T2D.

METHODS

Physically inactive age and BMI-similar premenopausal women and men with ( n  = 22) and without [ n  = 34, controls with overweight/obesity (CWO)] uncomplicated T2D were compared to an age-similar and sex-similar reference control cohort ( n  = 20). Left ventricular (LV) structure, function, and aortic stiffness were assessed by MRI. Global arterial pulse wave velocity (PWV) was assessed using carotid-to-femoral applanation tonometry. Regional PWV was measured via 2D phase-contrast MRI and 4D flow MRI.

RESULTS

Global arterial PWV did not differ between CWO and T2D. 2D PC-MRI PWV in the ascending aorta was higher in people with T2D compared with CWOs ( P  < 0.01). 4D flow PWV in the thoracic aorta was higher in CWO ( P  < 0.01), and T2D ( P  < 0.001) compared with RC. End-diastolic volume, end-systolic volume, stroke volume, and cardiac output were lower in CWO and T2D groups compared with reference control.

CONCLUSION

Subclinical changes in arterial stiffening and cardiac remodeling in inactive CWO and T2D compared with reference control support obesity and/or physical inactivity as determinants of incipient CVD complications in uncomplicated T2D. Future studies should determine the mechanistic causes of the CVD complications in greater detail in order to create therapeutic targets.

CLINICAL TRIAL REGISTRATION

Cardiovascular Mechanisms of Exercise Intolerance in Diabetes and the Role of Sex (NCT03419195).

Assessing atrial myopathy with cardiac magnetic resonance imaging in embolic stroke of undetermined source

BACKGROUND

Left atrial myopathy has been implicated in atrial fibrillation (AF)-related stroke and embolic stroke of undetermined source (ESUS).

OBJECTIVE

To use advanced cardiac magnetic resonance (CMR) imaging techniques, including left atrial (LA) strain and 4D flow CMR, to identify atrial myopathy in patients with ESUS.

METHODS

20 patients with ESUS and no AF or other cause for stroke, and 20 age and sex-matched controls underwent CMR with 4D flow analysis. Markers of LA myopathy were assessed including LA size, volume, ejection fraction, and strain. 4D flow CMR was performed to measure novel markers of LA stasis such as LA velocities and the LA residence time distribution time constant (RTDtc). These markers of LA myopathy were compared between the two groups.

RESULTS

There was no significant difference in: CMR-calculated LA velocities or LA total, passive or active ejection fractions between the groups. There was no significant difference in CMR-derived reservoir, conduit or contractile average longitudinal strain between the ESUS and control groups (22.9 vs 22.6%, p=0.379, 11.2 ± 3.5 vs 12.4 ± 2.6% p=0.224, 10.8 ± 3.2 vs 10.4 ± 2.3%, p=0.625 respectively). Similarly, RTDtc was not significantly longer in ESUS patients compared to controls (1.3 ± 0.2 vs 1.2 ± 0.2, p=0.1).

CONCLUSIONS

There were no significant differences in any CMR marker of atrial myopathy in ESUS patients compared to healthy controls, likely reflecting the multiple possible aetiologies of ESUS suggesting that the role LA myopathy plays in ESUS is smaller than previously thought.

Left atrium 4D-flow segmentation with high-resolution contrast-enhanced magnetic resonance angiography

Background

Atrial fibrillation (AF) leads to intracardiac thrombus and an associated risk of stroke. Phase-contrast cardiovascular magnetic resonance (CMR) with flow-encoding in all three spatial directions (4D-flow) provides a time-resolved 3D volume image with 3D blood velocity, which brings individual hemodynamic information affecting thrombus formation. As the resolution and contrast of 4D-flow are limited, we proposed a semi-automated 4D-flow segmentation method for the left atrium (LA) using a standard-of-care contrast-enhanced magnetic resonance angiography (CE-MRA) and registration technique.

Methods

LA of 54 patients with AF were segmented from 4D-flow taken in sinus rhythm using two segmentation methods. (1) Phase-contrast magnetic resonance angiography (PC-MRA) was calculated from 4D-flow, and LA was segmented slice-by-slice manually. (2) LA and other structures were segmented from CE-MRA and transformed into 4D-flow coordinates by registration with the mutual information method. Overlap of volume was tested by the Dice similarity coefficient (DSC) and the average symmetric surface distance (ASSD). Mean velocity and stasis were calculated to compare the functional property of LA from two segmentation methods.

Results

LA volumes from segmentation on CE-MRA were strongly correlated with PC-MRA volume, although mean CE-MRA volumes were about 10% larger. The proposed registration scheme resulted in visually successful registration in 76% of cases after two rounds of registration. The mean of DSC of the registered cases was 0.770 ± 0.045, and the mean of ASSD was 2.704 mm ± 0.668 mm. Mean velocity had no significant difference between the two segmentation methods, and mean stasis had a 3.3% difference.

Conclusion

The proposed CE-MRA segmentation and registration method can generate segmentation for 4D-flow images. This method will facilitate 4D-flow analysis for AF patients by making segmentation easier and overcoming the limit of resolution.

4D-Flow MRI intracardiac flow analysis considering different subtypes of pulmonary hypertension

Intracardiac flow hemodynamic patterns have been considered to be an early sign of diastolic dysfunction. In this study we investigated right ventricular (RV) diastolic dysfunction between patients with pulmonary arterial hypertension (PAH) and pulmonary hypertension with chronic lung disease (PH-CLD) via 4D-Flow cardiac MRI (CMR). Patients underwent prospective, comprehensive CMR for function and size including 4D-Flow CMR protocol for intracardiac flow visualization and analysis. RV early filling phase and peak atrial phase vorticity (E-vorticity and A-vorticity) values were calculated in all patients. Patients further underwent comprehensive Doppler and tissue Doppler evaluation for the RV diastolic dysfunction. In total 13 patients with PAH, 15 patients with PH-CLD, and 10 control subjects underwent the 4D-Flow CMR and echocardiography evaluation for RV diastolic dysfunction. Reduced E-vorticity differentiated PAH and PH-CLD from healthy controls (both p < 0.01) despite the same Doppler E values. E-vorticity was further decreased in PAH patients when compared to PH-CLD group (p < 0.05) with similar Doppler and tissue Doppler markers of diastolic dysfunction. A-vorticity was decreased in both PAH and PH-CLD groups compared to controls but with no difference between the disease groups. E-vorticity correlated with ejection fraction (R = 0.60, p < 0.001), end-systolic volume (R = 0.50, p = 0.001), stroke volume (R = 0.42, p = 0.007), and cardiac output (R = 0.30, p = 0.027). Intracardiac flow analysis using 4D-Flow CMR derived vorticity is a sensitive method to differentiate diastolic dysfunction in patients with different PH etiology and similar Doppler echocardiography profile.

Mitral valve regurgitation assessed by intraventricular CMR 4D-flow: a systematic review on the technological aspects and potential clinical applications

Cardiac magnetic resonance (CMR) four-dimensional (4D) flow is a novel method for flow quantification potentially helpful in management of mitral valve regurgitation (MVR). In this systematic review, we aimed to depict the clinical role of intraventricular 4D-flow in MVR. The reproducibility, technical aspects, and comparison against conventional techniques were evaluated. Published studies on SCOPUS, MEDLINE, and EMBASE were included using search terms on 4D-flow CMR in MVR. Out of 420 screened articles, 18 studies fulfilled our inclusion criteria. All studies (n = 18, 100%) assessed MVR using 4D-flow intraventricular annular inflow (4D-flowAIM) method, which calculates the regurgitation by subtracting the aortic forward flow from the mitral forward flow. Thereof, 4D-flow jet quantification (4D-flowjet) was assessed in 5 (28%), standard 2D phase-contrast (2D-PC) flow imaging in 8 (44%) and the volumetric method (the deviation of left ventricle stroke volume and right ventricular stroke volume) in 2 (11%) studies. Inter-method correlations among the 4 MVR quantification methods were heterogeneous across studies, ranging from moderate to excellent correlations. Two studies compared 4D-flowAIM to echocardiography with moderate correlation. In 12 (63%) studies the reproducibility of 4D-flow techniques in quantifying MVR was studied. Thereof, 9 (75%) studies investigated the reproducibility of the 4D-flowAIM method and the majority (n = 7, 78%) reported good to excellent intra- and inter-reader reproducibility. Intraventricular 4D-flowAIM provides high reproducibility with heterogeneous correlations to conventional quantification methods. Due to the absence of a gold standard and unknown accuracies, future longitudinal outcome studies are needed to assess the clinical value of 4D-flow in the clinical setting of MVR.

Comprehensive assessment of molecular function, tissue characterization, and hemodynamic performance by non-invasive hybrid imaging: Potential role of cardiac PETMR

Noninvasive cardiovascular imaging plays a key role in diagnosis and patient management including monitoring treatment efficacy. The usefulness of noninvasive cardiovascular imaging has been extensively studied and shown to have high diagnostic reliability and prognostic significance, while the nondiagnostic results frequently encountered with single imaging modality require complementary or alternative imaging techniques. Hybrid cardiac imaging was initially introduced to integrate anatomical and functional information to enhance the diagnostic performance, and lately employed as a strategy for comprehensive assessment of the underlying pathophysiology of diseases. More recently, the utility of computed tomography has grown in diversity, and emerged from being an exploratory technique allowing functional measurement such as stress dynamic perfusion. Cardiac magnetic resonance imaging (CMR) is widely accepted as a robust tool for evaluation of cardiac function, fibrosis, and edema, yielding high spatial resolution and soft-tissue contrast. However, the use of intravenous contrast materials is typically required for accurate diagnosis with these imaging modalities, despite the associated risk of renal toxicity. Nuclear cardiology, established as a molecular imaging technique, has advantages in visualization of the disease-specific biological process at cellular level using numerous probes without requiring contrast materials. Various imaging modalities should be appropriately used sequentially to assess concomitant disease and the progression over time. Therefore, simultaneous evaluation combining high spatial resolution and disease-specific imaging probe is a useful approach to identify the regional activity and the stage of the disease. Given the recent advance and potential of multiparametric CMR and novel nuclide tracers, hybrid positron emission tomography MR is becoming an ideal tool for disease-specific imaging.

Flow evaluation software for four-dimensional flow MRI: a reliability and validation study

PURPOSE

Four-dimensional time-resolved phase-contrast cardiovascular magnetic resonance imaging (4D flow MRI) enables blood flow quantification in multiple vessels, which is crucial for patients with congenital heart disease (CHD). We investigated net flow volumes in the ascending aorta and pulmonary arteries by four different postprocessing software packages for 4D flow MRI in comparison with 2D cine phase-contrast measurements (2D PC).

MATERIAL AND METHODS

4D flow and 2D PC datasets of 47 patients with biventricular CHD (median age 16, range 0.6-52 years) were acquired at 1.5 T. Net flow volumes in the ascending aorta, the main, right, and left pulmonary arteries were measured using four different postprocessing software applications and compared to offset-corrected 2D PC data. Reliability of 4D flow postprocessing software was assessed by Bland-Altman analysis and intraclass correlation coefficient (ICC). Linear regression of internal flow controls was calculated. Interobserver reproducibility was evaluated in 25 patients.

RESULTS

Correlation and agreement of flow volumes were very good for all software compared to 2D PC (ICC ≥ 0.94; bias ≤ 5%). Internal controls were excellent for 2D PC (r ≥ 0.95, p < 0.001) and 4D flow (r ≥ 0.94, p < 0.001) without significant difference of correlation coefficients between methods. Interobserver reliability was good for all vendors (ICC ≥ 0.94, agreement bias < 8%).

CONCLUSION

Haemodynamic information from 4D flow in the large thoracic arteries assessed by four commercially available postprocessing applications matches routinely performed 2D PC values. Therefore, we consider 4D flow MRI-derived data ready for clinical use in patients with CHD.

Fluid-structure interactions of peripheral arteries using a coupled in silico and in vitro approach

Vascular compliance is considered both a cause and a consequence of cardiovascular disease and a significant factor in the mid- and long-term patency of vascular grafts. However, the biomechanical effects of localised changes in compliance cannot be satisfactorily studied with the available medical imaging technologies or surgical simulation materials. To address this unmet need, we developed a coupled silico-vitro platform which allows for the validation of numerical fluid-structure interaction results as a numerical model and physical prototype. This numerical one-way and two-way fluid-structure interaction study is based on a three-dimensional computer model of an idealised femoral artery which is validated against patient measurements derived from the literature. The numerical results are then compared with experimental values collected from compliant arterial phantoms via direct pressurisation and ring tensile testing. Phantoms within a compliance range of 1.4-68.0%/100 mmHg were fabricated via additive manufacturing and silicone casting, then mechanically characterised via ring tensile testing and optical analysis under direct pressurisation with moderately statistically significant differences in measured compliance ranging between 10 and 20% for the two methods. One-way fluid-structure interaction coupling underestimated arterial wall compliance by up to 14.7% compared with two-way coupled models. Overall, Solaris™ (Smooth-On) matched the compliance range of the numerical and in vivo patient models most closely out of the tested silicone materials. Our approach is promising for vascular applications where mechanical compliance is especially important, such as the study of diseases which commonly affect arterial wall stiffness, such as atherosclerosis, and the model-based design, surgical training, and optimisation of vascular prostheses.

Intraventricular 4D flow cardiovascular magnetic resonance for assessing patients with heart failure with preserved ejection fraction: a pilot study

Diagnosing heart failure with preserved ejection fraction (HFpEF) remains challenging. Intraventricular four-dimensional flow (4D flow) phase-contrast cardiovascular magnetic resonance (CMR) can assess different components of left ventricular (LV) flow including direct flow, delayed ejection, retained inflow and residual volume. This could be utilised to identify HFpEF. This study investigated if intraventricular 4D flow CMR could differentiate HFpEF patients from non-HFpEF and asymptomatic controls. Suspected HFpEF patients and asymptomatic controls were recruited prospectively. HFpEF patients were confirmed using European Society of Cardiology (ESC) 2021 expert recommendations. Non-HFpEF patients were diagnosed if suspected HFpEF patients did not fulfil ESC 2021 criteria. LV direct flow, delayed ejection, retained inflow and residual volume were obtained from 4D flow CMR images. Receiver operating characteristic (ROC) curves were plotted. 63 subjects (25 HFpEF patients, 22 non-HFpEF patients and 16 asymptomatic controls) were included in this study. 46% were male, mean age 69.8 ± 9.1 years. CMR 4D flow derived LV direct flow and residual volume could differentiate HFpEF vs combined group of non-HFpEF and asymptomatic controls (p < 0.001 for both) as well as HFpEF vs non-HFpEF patients (p = 0.021 and p = 0.005, respectively). Among the 4 parameters, direct flow had the largest area under curve (AUC) of 0.781 when comparing HFpEF vs combined group of non-HFpEF and asymptomatic controls, while residual volume had the largest AUC of 0.740 when comparing HFpEF and non-HFpEF patients. CMR 4D flow derived LV direct flow and residual volume show promise in differentiating HFpEF patients from non-HFpEF patients.

The Inner Complexities of Multimodal Medical Data: Bitmap-Based 3D Printing for Surgical Planning Using Dynamic Physiology

Motivated by the need to develop more informative and data-rich patient-specific presurgical planning models, we present a high-resolution method that enables the tangible replication of multimodal medical data. By leveraging voxel-level control of multimaterial three-dimensional (3D) printing, our method allows for the digital integration of disparate medical data types, such as functional magnetic resonance imaging, tractography, and four-dimensional flow, overlaid upon traditional magnetic resonance imaging and computed tomography data. While permitting the explicit translation of multimodal medical data into physical objects, this approach also bypasses the need to process data into mesh-based boundary representations, alleviating the potential loss and remodeling of information. After evaluating the optical characteristics of test specimens generated with our correlative data-driven method, we culminate with multimodal real-world 3D-printed examples, thus highlighting current and potential applications for improved surgical planning, communication, and clinical decision-making through this approach.

Age- and sex-specific reference values of biventricular flow components and kinetic energy by 4D flow cardiovascular magnetic resonance in healthy subjects

BACKGROUND

Advances in four-dimensional flow cardiovascular magnetic resonance (4D flow CMR) have allowed quantification of left ventricular (LV) and right ventricular (RV) blood flow. We aimed to (1) investigate age and sex differences of 4D flow CMR-derived LV and RV relative flow components and kinetic energy (KE) parameters indexed to end-diastolic volume (KEiEDV) in healthy subjects; and (2) assess the effects of age and sex on these parameters.

METHODS

We performed 4D flow analysis in 163 healthy participants (42% female; mean age 43 ± 13 years) of a prospective registry study (NCT03217240) who were free of cardiovascular diseases. Relative flow components (direct flow, retained inflow, delayed ejection flow, residual volume) and multiple phasic KEiEDV (global, peak systolic, average systolic, average diastolic, peak E-wave, peak A-wave) for both LV and RV were analysed.

RESULTS

Compared with men, women had lower median LV and RV residual volume, and LV peak and average systolic KEiEDV, and higher median values of RV direct flow, RV global KEiEDV, RV average diastolic KEiEDV, and RV peak E-wave KEiEDV. ANOVA analysis found there were no differences in flow components, peak and average systolic, average diastolic and global KEiEDV for both LV and RV across age groups. Peak A-wave KEiEDV increased significantly (r = 0.458 for LV and 0.341 for RV), whereas peak E-wave KEiEDV (r = - 0.355 for LV and - 0.318 for RV), and KEiEDV E/A ratio (r = - 0.475 for LV and - 0.504 for RV) decreased significantly, with age.

CONCLUSION

These data using state-of-the-art 4D flow CMR show that biventricular flow components and kinetic energy parameters vary significantly by age and sex. Age and sex trends should be considered in the interpretation of quantitative measures of biventricular flow. Clinical trial registration  https://www.

CLINICALTRIALS

gov . Unique identifier: NCT03217240.

Haemodynamic effects of hypertension and type 2 diabetes: Insights from a 4D flow MRI-based personalized cardiovascular mathematical model

Type 2 diabetes (T2D) and hypertension increase the risk of cardiovascular diseases mediated by whole-body changes to metabolism, cardiovascular structure and haemodynamics. The haemodynamic changes related to hypertension and T2D are complex and subject-specific, however, and not fully understood. We aimed to investigate the haemodynamic mechanisms in T2D and hypertension by comparing the haemodynamics between healthy controls and subjects with T2D, hypertension, or both. For all subjects, we combined 4D flow magnetic resonance imaging data, brachial blood pressure and a cardiovascular mathematical model to create a comprehensive subject-specific analysis of central haemodynamics. When comparing the subject-specific haemodynamic parameters between the four groups, the predominant haemodynamic difference is impaired left ventricular relaxation in subjects with both T2D and hypertension compared to subjects with only T2D, only hypertension and controls. The impaired relaxation indicates that, in this cohort, the long-term changes in haemodynamic load of co-existing T2D and hypertension cause diastolic dysfunction demonstrable at rest, whereas either disease on its own does not. However, through subject-specific predictions of impaired relaxation, we show that altered relaxation alone is not enough to explain the subject-specific and group-related differences; instead, a combination of parameters is affected in T2D and hypertension. These results confirm previous studies that reported more adverse effects from the combination of T2D and hypertension compared to either disease on its own. Furthermore, this shows the potential of personalized cardiovascular models in providing haemodynamic mechanistic insights and subject-specific predictions that could aid in the understanding and treatment planning of patients with T2D and hypertension. KEY POINTS: The combination of 4D flow magnetic resonance imaging data and a cardiovascular mathematical model allows for a comprehensive analysis of subject-specific haemodynamic parameters that otherwise cannot be derived non-invasively. Using this combination, we show that diastolic dysfunction in subjects with both type 2 diabetes (T2D) and hypertension is the main group-level difference between controls, subjects with T2D, subjects with hypertension, and subjects with both T2D and hypertension. These results suggest that, in this relatively healthy population, the additional load of both hypertension and T2D affects the haemodynamic function of the left ventricle, whereas each disease on its own is not enough to cause significant effects under resting conditions. Finally, using the subject-specific model, we show that the haemodynamic effects of diastolic dysfunction alone are not sufficient to explain all the observed haemodynamic differences. Instead, additional subject-specific variations in cardiac and vascular function combine to explain the complex haemodynamics of subjects affected by hypertension and/or T2D.

A Vector Fitting Approach for the Automated Estimation of Lumped Boundary Conditions of 1D Circulation Models

PURPOSE

The choice of appropriate boundary conditions is a crucial step in the development of cardiovascular models for blood flow simulations. The three-element Windkessel model is usually employed as a lumped boundary condition, providing a reduced order representation of the peripheral circulation. However, the systematic estimation of the Windkessel parameters remains an open problem. Moreover, the Windkessel model is not always adequate to model blood flow dynamics, which often require more elaborate boundary conditions. In this study, we propose a method for the estimation of the parameters of high order boundary conditions, including the Windkessel model, from pressure and flow rate waveforms at the truncation point. Moreover, we investigate the effect of adopting higher order boundary conditions, corresponding to equivalent circuits with more than one storage element, on the accuracy of the model.

METHOD

The proposed technique is based on Time-Domain Vector Fitting, a modeling algorithm that, given samples of the input and output of a system, such as pressure and flow waveforms, can derive a differential equation approximating their relation.

RESULTS

The capabilities of the proposed method are tested on a 1D circulation model consisting of the 55 largest human systemic arteries, to demonstrate its accuracy and its usefulness to estimate boundary conditions with order higher than the traditional Windkessel models. The proposed method is compared to other common estimation techniques, and its robustness in parameter estimation is verified in presence of noisy data and of physiological changes of aortic flow rate induced by mental stress.

CONCLUSION

Results suggest that the proposed method is able to accurately estimate boundary conditions of arbitrary order. Higher order boundary conditions can improve the accuracy of cardiovascular simulations, and Time-Domain Vector Fitting can automatically estimate them.

Automatic 4D Flow MRI Segmentation Using the Standardized Difference of Means Velocity

We present a method to automatically segment 4D flow magnetic resonance imaging (MRI) by identifying net flow effects using the standardized difference of means (SDM) velocity. The SDM velocity quantifies the ratio between the net flow and observed flow pulsatility in each voxel. Vessel segmentation is performed using an F-test, identifying voxels with significantly higher SDM velocity values than background voxels. We compare the SDM segmentation algorithm against pseudo-complex difference (PCD) intensity segmentation of 4D flow measurements in in vitro cerebral aneurysm models and 10 in vitro Circle of Willis (CoW) datasets. We also compared the SDM algorithm to convolutional neural network (CNN) segmentation in 5 thoracic vasculature datasets. The in vitro flow phantom geometry is known, while the ground truth geometries for the CoW and thoracic aortas are derived from high-resolution time-of-flight (TOF) magnetic resonance angiography and manual segmentation, respectively. The SDM algorithm demonstrates greater robustness than PCD and CNN approaches and can be applied to 4D flow data from other vascular territories. The SDM to PCD comparison demonstrated an approximate 48% increase in sensitivity in vitro and 70% increase in the CoW, respectively; the SDM and CNN sensitivities were similar. The vessel surface derived from the SDM method was 46% closer to the in vitro surfaces and 72% closer to the in vitro TOF surfaces than the PCD approach. The SDM and CNN approaches both accurately identify vessel surfaces. The SDM algorithm is a repeatable segmentation method, enabling reliable computation of hemodynamic metrics associated with cardiovascular disease.

Four-Dimensional Flow MR Imaging: Technique and Advances

4D Flow MRI is an advanced imaging technique for comprehensive non-invasive assessment of the cardiovascular system. The capture of the blood velocity vector field throughout the cardiac cycle enables measures of flow, pulse wave velocity, kinetic energy, wall shear stress, and more. Advances in hardware, MRI data acquisition and reconstruction methodology allow for clinically feasible scan times. The availability of 4D Flow analysis packages allows for more widespread use in research and the clinic and will facilitate much needed multi-center, multi-vendor studies in order to establish consistency across scanner platforms and to enable larger scale studies to demonstrate clinical value.

Clinical Applications of Four-Dimensional Flow MRI

Four-dimensional flow MRI is a powerful phase contrast technique used for assessing three-dimensional (3D) blood flow dynamics. By acquiring a time-resolved velocity field, it enables flexible retrospective analysis of blood flow that can include qualitative 3D visualization of complex flow patterns, comprehensive assessment of multiple vessels, reliable placement of analysis planes, and calculation of advanced hemodynamic parameters. This technique provides several advantages over routine two-dimensional flow imaging techniques, allowing it to become part of clinical practice at major academic medical centers. In this review, we present the current state-of-the-art cardiovascular, neurovascular, and abdominal applications.

Mechanisms of Aortic Dilation in Patients With Bicuspid Aortic Valve: JACC State-of-the-Art Review

Bicuspid aortic valve is the most common congenital heart disease and exposes patients to an increased risk of aortic dilation and dissection. Aortic dilation is a slow, silent process, leading to a greater risk of aortic dissection. The prevention of adverse events together with optimization of the frequency of the required lifelong imaging surveillance are important for both clinicians and patients and motivated extensive research to shed light on the physiopathologic processes involved in bicuspid aortic valve aortopathy. Two main research hypotheses have been consolidated in the last decade: one supports a genetic basis for the increased prevalence of dilation, in particular for the aortic root, and the second supports the damaging impact on the aortic wall of altered flow dynamics associated with these structurally abnormal valves, particularly significant in the ascending aorta. Current opinion tends to rule out mutually excluding causative mechanisms, recognizing both as important and potentially clinically relevant.