4D flow imaging with MRI

Validation of ultrasound velocimetry and computational fluid dynamics for flow assessment in femoral artery stenotic disease

Purpose

To investigate the accuracy of high-framerate echo particle image velocimetry (ePIV) and computational fluid dynamics (CFD) for determining velocity vectors in femoral bifurcation models through comparison with optical particle image velocimetry (oPIV).

Approach

Separate femoral bifurcation models were built for oPIV and ePIV measurements of a non-stenosed (control) and a 75%-area stenosed common femoral artery. A flow loop was used to create triphasic pulsatile flow. In-plane velocity vectors were measured with oPIV and ePIV. Flow was simulated with CFD using boundary conditions from ePIV and additional duplex-ultrasound (DUS) measurements. Mean differences and 95%-limits of agreement (1.96*SD) of the velocity magnitudes in space and time were compared, and the similarity of vector complexity (VC) and time-averaged wall shear stress (TAWSS) was assessed.

Results

Similar flow features were observed between modalities with velocities up to 110 and 330    cm / s in the control and the stenosed model, respectively. Relative to oPIV, ePIV and CFD-ePIV showed negligible mean differences in velocity (< 3    cm / s ), with limits of agreement of ± 25    cm / s (control) and ± 34    cm / s (stenosed). CFD-DUS overestimated velocities with limits of agreements of 13 ± 40 and 16.1 ± 55    cm / s for the control and stenosed model, respectively. VC showed good agreement, whereas TAWSS showed similar trends but with higher values for ePIV, CFD-DUS, and CFD-ePIV compared to oPIV.

Conclusions

EPIV and CFD-ePIV can accurately measure complex flow features in the femoral bifurcation and around a stenosis. CFD-DUS showed larger deviations in velocities making it a less robust technique for hemodynamical assessment. The applied ePIV and CFD techniques enable two- and three-dimensional assessment of local hemodynamics with high spatiotemporal resolution and thereby overcome key limitations of current clinical modalities making them an attractive and cost-effective alternative for hemodynamical assessment in clinical practice.

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..

Comparison of computational fluid dynamics with transcranial Doppler ultrasound in response to physiological stimuli

Cerebrovascular haemodynamics are sensitive to multiple physiological stimuli that require synergistic response to maintain adequate perfusion. Understanding haemodynamic changes within cerebral arteries is important to inform how the brain regulates perfusion; however, methods for direct measurement of cerebral haemodynamics in these environments are challenging. The aim of this study was to assess velocity waveform metrics obtained using transcranial Doppler (TCD) with flow-conserving subject-specific three-dimensional (3D) simulations using computational fluid dynamics (CFD). Twelve healthy participants underwent head and neck imaging with 3 T magnetic resonance angiography. Velocity waveforms in the middle cerebral artery were measured with TCD ultrasound, while diameter and velocity were measured using duplex ultrasound in the internal carotid and vertebral arteries to calculate incoming cerebral flow at rest, during hypercapnia and exercise. CFD simulations were developed for each condition, with velocity waveform metrics extracted in the same insonation region as TCD. Exposure to stimuli induced significant changes in cardiorespiratory measures across all participants. Measured absolute TCD velocities were significantly higher than those calculated from CFD (P range < 0.001-0.004), and these data were not correlated across conditions (r range 0.030-0.377, P range 0.227-0.925). However, relative changes in systolic and time-averaged velocity from resting levels exhibited significant positive correlations when the distinct techniques were compared (r range 0.577-0.770, P range 0.003-0.049). Our data indicate that while absolute measures of cerebral velocity differ between TCD and 3D CFD simulation, physiological changes from resting levels in systolic and time-averaged velocity are significantly correlated between techniques.

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.

A Deep Learning Approach to Using Wearable Seismocardiography (SCG) for Diagnosing Aortic Valve Stenosis and Predicting Aortic Hemodynamics Obtained by 4D Flow MRI

In this paper, we explored the use of deep learning for the prediction of aortic flow metrics obtained using 4-dimensional (4D) flow magnetic resonance imaging (MRI) using wearable seismocardiography (SCG) devices. 4D flow MRI provides a comprehensive assessment of cardiovascular hemodynamics, but it is costly and time-consuming. We hypothesized that deep learning could be used to identify pathological changes in blood flow, such as elevated peak systolic velocity ([Formula: see text]) in patients with heart valve diseases, from SCG signals. We also investigated the ability of this deep learning technique to differentiate between patients diagnosed with aortic valve stenosis (AS), non-AS patients with a bicuspid aortic valve (BAV), non-AS patients with a mechanical aortic valve (MAV), and healthy subjects with a normal tricuspid aortic valve (TAV). In a study of 77 subjects who underwent same-day 4D flow MRI and SCG, we found that the [Formula: see text] values obtained using deep learning and SCGs were in good agreement with those obtained by 4D flow MRI. Additionally, subjects with non-AS TAV, non-AS BAV, non-AS MAV, and AS could be classified with ROC-AUC (area under the receiver operating characteristic curves) values of 92%, 95%, 81%, and 83%, respectively. This suggests that SCG obtained using low-cost wearable electronics may be used as a supplement to 4D flow MRI exams or as a screening tool for aortic valve disease.

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.

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.

Reconstruction and Validation of Arterial Geometries for Computational Fluid Dynamics Using Multiple Temporal Frames of 4D Flow-MRI Magnitude Images

PURPOSE

Segmentation and reconstruction of arterial blood vessels is a fundamental step in the translation of computational fluid dynamics (CFD) to the clinical practice. Four-dimensional flow magnetic resonance imaging (4D Flow-MRI) can provide detailed information of blood flow but processing this information to elucidate the underlying anatomical structures is challenging. In this study, we present a novel approach to create high-contrast anatomical images from retrospective 4D Flow-MRI data.

METHODS

For healthy and clinical cases, the 3D instantaneous velocities at multiple cardiac time steps were superimposed directly onto the 4D Flow-MRI magnitude images and combined into a single composite frame. This new Composite Phase-Contrast Magnetic Resonance Angiogram (CPC-MRA) resulted in enhanced and uniform contrast within the lumen. These images were subsequently segmented and reconstructed to generate 3D arterial models for CFD. Using the time-dependent, 3D incompressible Reynolds-averaged Navier-Stokes equations, the transient aortic haemodynamics was computed within a rigid wall model of patient geometries.

RESULTS

Validation of these models against the gold standard CT-based approach showed no statistically significant inter-modality difference regarding vessel radius or curvature (p > 0.05), and a similar Dice Similarity Coefficient and Hausdorff Distance. CFD-derived near-wall hemodynamics indicated a significant inter-modality difference (p > 0.05), though these absolute errors were small. When compared to the in vivo data, CFD-derived velocities were qualitatively similar.

CONCLUSION

This proof-of-concept study demonstrated that functional 4D Flow-MRI information can be utilized to retrospectively generate anatomical information for CFD models in the absence of standard imaging datasets and intravenous contrast.

Methods of ex vivo analysis of tissue status in vascularized composite allografts

Vascularized composite allotransplantation can improve quality of life and restore functionality. However, the complex tissue composition of vascularized composite allografts (VCAs) presents unique clinical challenges that increase the likelihood of transplant rejection. Under prolonged static cold storage, highly damage-susceptible tissues such as muscle and nerve undergo irreversible degradation that may render allografts non-functional. Skin-containing VCA elicits an immunogenic response that increases the risk of recipient allograft rejection. The development of quantitative metrics to evaluate VCAs prior to and following transplantation are key to mitigating allograft rejection. Correspondingly, a broad range of bioanalytical methods have emerged to assess the progression of VCA rejection and characterize transplantation outcomes. To consolidate the current range of relevant technologies and expand on potential for development, methods to evaluate ex vivo VCA status are herein reviewed and comparatively assessed. The use of implantable physiological status monitoring biochips, non-invasive bioimpedance monitoring to assess edema, and deep learning algorithms to fuse disparate inputs to stratify VCAs are identified.

4D Flow cardiovascular magnetic resonance consensus statement: 2023 update

Hemodynamic assessment is an integral part of the diagnosis and management of cardiovascular disease. Four-dimensional cardiovascular magnetic resonance flow imaging (4D Flow CMR) allows comprehensive and accurate assessment of flow in a single acquisition. This consensus paper is an update from the 2015 '4D Flow CMR Consensus Statement'. We elaborate on 4D Flow CMR sequence options and imaging considerations. The document aims to assist centers starting out with 4D Flow CMR of the heart and great vessels with advice on acquisition parameters, post-processing workflows and integration into clinical practice. Furthermore, we define minimum quality assurance and validation standards for clinical centers. We also address the challenges faced in quality assurance and validation in the research setting. We also include a checklist for recommended publication standards, specifically for 4D Flow CMR. Finally, we discuss the current limitations and the future of 4D Flow CMR. This updated consensus paper will further facilitate widespread adoption of 4D Flow CMR in the clinical workflow across the globe and aid consistently high-quality publication standards.

Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D Flow-MRI

Introduction: Patient-specific computational fluid dynamics (CFD) models permit analysis of complex intra-aortic hemodynamics in patients with aortic dissection (AD), where vessel morphology and disease severity are highly individualized. The simulated blood flow regime within these models is sensitive to the prescribed boundary conditions (BCs), so accurate BC selection is fundamental to achieve clinically relevant results. Methods: This study presents a novel reduced-order computational framework for the iterative flow-based calibration of 3-Element Windkessel Model (3EWM) parameters to generate patient-specific BCs. These parameters were calibrated using time-resolved flow information derived from retrospective four-dimensional flow magnetic resonance imaging (4D Flow-MRI). For a healthy and dissected case, blood flow was then investigated numerically in a fully coupled zero dimensional-three dimensional (0D-3D) numerical framework, where the vessel geometries were reconstructed from medical images. Calibration of the 3EWM parameters was automated and required ~3.5 min per branch. Results: With prescription of the calibrated BCs, the computed near-wall hemodynamics (time-averaged wall shear stress, oscillatory shear index) and perfusion distribution were consistent with clinical measurements and previous literature, yielding physiologically relevant results. BC calibration was particularly important in the AD case, where the complex flow regime was captured only after BC calibration. Discussion: This calibration methodology can therefore be applied in clinical cases where branch flow rates are known, for example, via 4D Flow-MRI or ultrasound, to generate patient-specific BCs for CFD models. It is then possible to elucidate, on a case-by-case basis, the highly individualized hemodynamics which occur due to geometric variations in aortic pathology high spatiotemporal resolution through CFD.

Imaging the vasculature of a beating heart by dynamic speckle: the challenge of a quasiperiodic motion

The spatial and temporal evolution of the field backscattered by a beating heart while illuminated with a coherent light reveals its macro- and microvascularization in real time. To perform these vascularization images, we use a recently published method of laser speckle imaging, based on the selective detection of spatially depolarized speckle field that is mainly generated by multiple scattering. We consider the calculation of the speckle contrast, by a spatial or temporal estimation. We show that the signal-to-noise ratio of the observed vascular structure can be noticeably increased by a postprocessing method implying the calculation of a motion field that allows the selection of similar frames extracted from different heartbeat periods. This later optimization reveals vascular microstructures with a spatial resolution of the order of 100    μ m .

Deep Learning in MRI-guided Radiation Therapy: A Systematic Review

MRI-guided radiation therapy (MRgRT) offers a precise and adaptive approach to treatment planning. Deep learning applications which augment the capabilities of MRgRT are systematically reviewed. MRI-guided radiation therapy offers a precise, adaptive approach to treatment planning. Deep learning applications which augment the capabilities of MRgRT are systematically reviewed with emphasis placed on underlying methods. Studies are further categorized into the areas of segmentation, synthesis, radiomics, and real time MRI. Finally, clinical implications, current challenges, and future directions are discussed.

Aortic wall shear stress in bicuspid aortic valve disease-10-year follow-up

Background

Bicuspid aortic valve (BAV) disease leads to deviant helical flow patterns especially in the mid-ascending aorta (AAo), potentially causing wall alterations such as aortic dilation and dissection. Among others, wall shear stress (WSS) could contribute to the prediction of long-term outcome of patients with BAV. 4D flow in cardiovascular magnetic resonance (CMR) has been established as a valid method for flow visualization and WSS estimation. The aim of this study is to reevaluate flow patterns and WSS in patients with BAV 10 years after the initial evaluation.

Methods

Fifteen patients (median age 34.0 years) with BAV were re-evaluated 10 years after the initial study from 2008/2009 using 4D flow by CMR. Our particular patient cohort met the same inclusion criteria as in 2008/2009, all without enlargement of the aorta or valvular impairment at that time. Flow patterns, aortic diameters, WSS and distensibility were calculated in different aortic regions of interest (ROI) with dedicated software tools.

Results

Indexed aortic diameters in the descending aorta (DAo), but especially in the AAo did not change in the 10-year period. Median difference 0.05 cm/m² (95% CI: 0.01 to 0.22; P=0.06) for AAo and median difference -0.08 cm/m² (95% CI: -0.12 to 0.01; P=0.07) for DAo. WSS values were lower in 2018/2019 at all measured levels. Aortic distensibility decreased by median 25.6% in the AAo, while stiffness increased concordantly (median +23.6%).

Conclusions

After a ten years' follow-up of patients with isolated BAV disease, indexed aortic diameters did not change in this patient cohort. WSS was lower compared to values generated 10 years earlier. Possibly a drop of WSS in BAV could serve as a marker for a benign long-term course and implementation of more conservative treatment strategies.

Multivariable Technique for the Evaluation of the Trans-stenotic Pressure Gradient

PURPOSE

This study establishes a reliable image-based multivariable technique for measuring the trans-stenotic pressure gradient.

METHODS

A self-made in vitro steady flow model based on adjustable velocities and stenotic properties were used as the experimental subject. The pre-stenotic flow velocity, severity, and length of the stenosis were used as the input variables. Based on equations used to fit the plots of the physically measured pressure gradient values versus each input variable, a multivariable formula for the pressure gradient measurement could then be derived. The flow model was scanned using velocity-encoded phase-contrast magnetic resonance imaging (PC-MRI) to validate the derived formula while simultaneously measuring the trans-stenotic pressure gradient. The correlation between the physically-measured pressure gradient values and the pressure gradient values calculated using the new formula were subsequently analyzed.

RESULTS

The results of linear regression analysis using the physically measured pressure gradient values for the new method were compared to values obtained using the simplified Bernoulli equation (R², 0.991, and 0.975, respectively). In a paired t-test, no statistically significant difference was found between the new method and the physical measurements.

CONCLUSIONS

The derived multivariable technique was found to reliably measure the trans-stenotic pressure gradient, with better performance than a traditional procedure based on the simplified Bernoulli equation.

Accelerated 4D-flow MRI with 3-point encoding enabled by machine learning

PURPOSE

To investigate the acceleration of 4D-flow MRI using a convolutional neural network (CNN) that produces three directional velocities from three flow encodings, without requiring a fourth reference scan measuring background phase.

METHODS

A fully 3D CNN using a U-net architecture was trained in a block-wise fashion to take complex images from three flow encodings and to produce three real-valued images for each velocity component. Using neurovascular 4D-flow scans (n = 144), the CNN was trained to predict velocities computed from four flow encodings by standard reconstruction including correction for residual background phase offsets. Methods to optimize loss functions were investigated, including magnitude, complex difference, and uniform velocity weightings. Subsequently, 3-point encoding was evaluated using cross validation of pixelwise correlation, flow measurements in major arteries, and in experiments with data at differing acceleration rates than the training data.

RESULTS

The CNN-produced 3-point velocities showed excellent agreements with the 4-point velocities, both qualitatively in velocity images, in flow rate measures, and quantitatively in regression analysis (slope = 0.96, R²  = 0.992). Optimizing the training to focus on vessel velocities rather than the global velocity error and improved the correlation of velocity within vessels themselves. The lowest error was observed when the loss function used uniform velocity weighting, in which the magnitude-weighted inverse of the velocity frequency uniformly distributed weighting across all velocity ranges. When applied to highly accelerated data, the 3-point network maintained a high correlation with ground truth data and demonstrated a denoising effect.

CONCLUSION

The 4D-flow MRI can be accelerated using machine learning requiring only three flow encodings to produce three-directional velocity maps with small errors.

Recent developments in modeling, imaging, and monitoring of cardiovascular diseases using machine learning

Cardiovascular diseases are the leading cause of mortality, morbidity, and hospitalization around the world. Recent technological advances have facilitated analyzing, visualizing, and monitoring cardiovascular diseases using emerging computational fluid dynamics, blood flow imaging, and wearable sensing technologies. Yet, computational cost, limited spatiotemporal resolution, and obstacles for thorough data analysis have hindered the utility of such techniques to curb cardiovascular diseases. We herein discuss how leveraging machine learning techniques, and in particular deep learning methods, could overcome these limitations and offer promise for translation. We discuss the remarkable capacity of recently developed machine learning techniques to accelerate flow modeling, enhance the resolution while reduce the noise and scanning time of current blood flow imaging techniques, and accurate detection of cardiovascular diseases using a plethora of data collected by wearable sensors.

A hybrid hierarchical strategy for registration of 7T TOF-MRI to 7T PC-MRI intracranial vessel data

PURPOSE

7T time-of-flight (TOF) MRI provides high resolution for the evaluation of cerebrovascular vessels and pathologies. In combination with 4D flow fields acquired with phase-contrast (PC) MRI, hemodynamic information can be extracted to enhance the analysis by providing direct measurements in the larger arteries or patient-specific boundary conditions. Hence, a registration between both modalities is required.

METHODS

To combine TOF and PC-MRI data, we developed a hybrid registration approach. Vessels and their centerlines are segmented from the TOF data. The centerline is fit to the intensity ridges of the lower resolved PC-MRI data, which provides temporal information. We used a metric that utilizes a scaled sum of weighted intensities and gradients on the normal plane. The registration is then guided by decoupled local affine transformations. It is applied hierarchically following the branching order of the vessel tree.

RESULTS

A landmark validation over Monte Carlo simulations yielded an average mean squared error of 184.73 mm and an average Hausdorff distance of 15.20 mm. The hierarchical traversal that transforms child vessels with their parents registers even small vessels not detectable in the PC-MRI.

CONCLUSION

The presented work combines high-resolution tomographic information from 7T TOF-MRI and measured flow data from 4D 7T PC-MRI scan for the arteries of the brain. This enables usage of patient-specific flow parameters for realistic simulations, thus supporting research in areas such as cerebral small vessel disease. Automatization and free deformations can help address the limiting error measures in the future.

Fetal 4D flow MRI of the great thoracic vessels at 3 Tesla using Doppler-ultrasound gating: a feasibility study

OBJECTIVES

To evaluate the feasibility of Doppler-ultrasound (DUS)-gated 4D flow MRI of the fetal great thoracic vessels at 3T in a clinical setting.

METHODS

Sixteen consecutive fetuses (range 30⁺⁴-38⁺⁵ weeks) with (n = 11) and without (n = 5) cardiovascular anomalies underwent 4D flow MRI of the great thoracic vessels at 3T. Direct fetal cardiac gating was obtained using a MR-compatible DUS device. 4D flow MRI-based visualisation and quantification of four target regions (ascending aorta (AAo), descending aorta (DAo), main pulmonary artery (MPA), and ductus arteriosus (DA)) were performed using dedicated software.

RESULTS

Fetal 4D flow MRI of the great thoracic vessels was successful in 12/16 fetuses (75%) by adopting clinical 4D flow MR protocols in combination with direct fetal cardiac DUS-gating. Four datasets were excluded due to artefacts by fetal movement or maternal breathing. 4D flow MRI-derived time-velocity curves revealed typical arterial blood flow patterns in the aorta. 4D flow quantification was achieved for the pre-defined target regions. Average velocity and flow volume were 21.1 ± 5.2 cm/s and 6.0 ± 3.1 mL/s in the AAo, 24.3 ± 6.7 cm/s and 8.4 ± 3.7 mL/s in the DAo, 21.9 ± 6.4 cm/s and 7.8 ± 4.2 mL/s in the MPA, and 23.4 ± 4.7 cm/s and 5.9 ± 3.6 mL/s in the DA, respectively.

CONCLUSIONS

Combination of DUS-gating of the fetal heart and 4D flow MRI allows comprehensive visualisation and quantification of haemodynamics in the fetal great thoracic vessels. DUS-gated fetal 4D flow MRI may provide a new diagnostic approach for prenatal assessment of blood flow haemodynamics.

KEY POINTS

• Fetal cardiac Doppler-ultrasound (DUS) gating and 4D flow MRI can be successfully combined. • DUS-gated fetal 4D flow MRI allowed visualisation and evaluation of streamline directionality, illustration of blood flow variations, and pulsatile arterial waveforms in the target vessels. • 4D flow MRI-based visualisation and quantification of the fetal great thoracic vessels were successful and flow metrics agreed with echocardiographic reference values.

Distribution and regional variation of wall shear stress in the curved middle cerebral artery using four-dimensional flow magnetic resonance imaging

BACKGROUND

To investigate the distribution and regional variation of wall shear stress (WSS) in the curved middle cerebral artery (MCA) in healthy individuals using four-dimensional (4D) flow magnetic resonance imaging (MRI).

METHODS

A total of 44 healthy participants (18 males; mean ages: 27.16±5.69 years) were included in this cross-sectional study. The WSS parameters of mean, minimum, and maximum values, the coefficient of variation of time-averaged WSS (TAWSSCV), and the maximum values of the oscillatory shear index (OSI) were calculated and compared in the curved proximal (M1) segments. Three cross-sectional planes were selected: the location perpendicular to the beginning of the long axis of the curved M1 segment of the MCA (proximal section), the most curved M1 location (curved M1 section), and the location before the insular (M2) segment bifurcation (distal section). The WSS and OSI parameters of the proximal, curved, and distal sections of the curved M1 segment were compared, including the inner and outer curvatures of the curved M1 section.

RESULTS

Of the curved M1 segments, the curved M1 section had significantly lower minimum TAWSS values than the proximal (P=0.031) and distal sections (P=0.002), and the curved M1 section had significantly higher maximum OSI values than the distal section (P=0.001). The TAWSSCV values at the curved M1 section were significantly higher than the proximal (P=0.001) and distal sections (P<0.001). At the curved M1 section, the inner curvature showed a significantly lower minimum TAWSS (P=0.013) and higher maximum OSI values (P=0.002) than the outer curvature.

CONCLUSIONS

There are distribution variation of WSS and OSI parameters at the curved M1 section of the curved MCA, and the inner curvature of the curved M1 section has the lowest WSS and highest OSI distribution. The local hemodynamic features of the curved MCA may be related to the predilection for atherosclerotic plaque development.

Wall Shear Stress Estimation for 4D Flow MRI Using Navier-Stokes Equation Correction

This study introduces a novel wall shear stress (WSS) estimation method for 4D flow MRI. The method improves the WSS accuracy by using the reconstructed pressure gradient and the flow-physics constraints to correct velocity gradient estimation. The method was tested on synthetic 4D flow data of analytical Womersley flow and flow in cerebral aneurysms and applied to in vivo 4D flow data acquired in cerebral aneurysms and aortas. The proposed method's performance was compared to the state-of-the-art method based on smooth-spline fitting of velocity profile and the WSS calculated from uncorrected velocity gradient. The proposed method improved the WSS accuracy by as much as 100% for the Womersley flow and reduced the underestimation of mean WSS by 39 to 50% for the synthetic aneurysmal flow. The predicted mean WSS from the in vivo aneurysmal data using the proposed method was 31 to 50% higher than the other methods. The predicted aortic WSS using the proposed method was 3 to 6 times higher than the other methods and was consistent with previous CFD studies and the results from recently developed methods that take into account the limited spatial resolution of 4D flow MRI. The proposed method improves the accuracy of WSS estimation from 4D flow MRI, which can help predict blood vessel remodeling and progression of cardiovascular diseases.

Cardiovascular deconditioning and impact of artificial gravity during 60-day head-down bed rest—Insights from 4D flow cardiac MRI

Microgravity has deleterious effects on the cardiovascular system. We evaluated some parameters of blood flow and vascular stiffness during 60 days of simulated microgravity in head-down tilt (HDT) bed rest. We also tested the hypothesis that daily exposure to 30 min of artificial gravity (1 g) would mitigate these adaptations. 24 healthy subjects (8 women) were evenly distributed in three groups: continuous artificial gravity, intermittent artificial gravity, or control. 4D flow cardiac MRI was acquired in horizontal position before (−9 days), during (5, 21, and 56 days), and after (+4 days) the HDT period. The false discovery rate was set at 0.05. The results are presented as median (first quartile; third quartile). No group or group × time differences were observed so the groups were combined. At the end of the HDT phase, we reported a decrease in the stroke volume allocated to the lower body (−30% [−35%; −22%]) and the upper body (−20% [−30%; +11%]), but in different proportions, reflected by an increased share of blood flow towards the upper body. The aortic pulse wave velocity increased (+16% [+9%; +25%]), and so did other markers of arterial stiffness (CAVI; CAVI0). In males, the time-averaged wall shear stress decreased (−13% [−17%; −5%]) and the relative residence time increased (+14% [+5%; +21%]), while these changes were not observed among females. Most of these parameters tended to or returned to baseline after 4 days of recovery. The effects of the artificial gravity countermeasure were not visible. We recommend increasing the load factor, the time of exposure, or combining it with physical exercise. The changes in blood flow confirmed the different adaptations occurring in the upper and lower body, with a larger share of blood volume dedicated to the upper body during (simulated) microgravity. The aorta appeared stiffer during the HDT phase, however all the changes remained subclinical and probably the sole consequence of reversible functional changes caused by reduced blood flow. Interestingly, some wall shear stress markers were more stable in females than in males. No permanent cardiovascular adaptations following 60 days of HDT bed rest were observed.

Pediatric magnetic resonance angiography: to contrast or not to contrast

Magnetic resonance (MR) angiography and MR venography imaging with contrast and non-contrast techniques are widely used for pediatric vascular imaging. However, as with any MRI examination, imaging the pediatric population can be challenging because of patient motion, which sometimes requires sedation. There are multiple benefits of non-contrast MR angiographic techniques, including the ability to repeat sequences if motion is present, the decreased need for sedation, and avoidance of potential risks associated with gadolinium administration and radiation exposure. Thus, MR angiography is an attractive alternative to CT or conventional catheter-based angiography in pediatric populations. Contrast-enhanced MR angiographic techniques have the advantage of increased signal to noise. Blood pool imaging allows long imaging times that result in high-spatial-resolution imaging, and thus high-quality diagnostic images. This article outlines the technique details, indications, benefits and downsides of non-contrast-enhanced and contrast-enhanced MR angiographic techniques to assist in protocol decision-making.

Compact pediatric cardiac magnetic resonance imaging protocols

Cardiac MRI is in many respects an ideal modality for pediatric cardiovascular imaging, enabling a complete noninvasive assessment of anatomy, morphology, function and flow in one radiation-free and potentially non-contrast exam. Nonetheless, traditionally lengthy and complex imaging acquisition strategies have often limited its broader use beyond specialized centers. In this review, the author presents practical cardiac MRI imaging protocols to facilitate the performance of succinct yet successful exams that provide the most salient clinical data for the majority of congenital and acquired pediatric cardiac disease. In addition, the author reviews newer and evolving techniques that permit more rapid but similarly diagnostic MRI, including compressed sensing and artificial intelligence/machine learning reconstruction, four-dimensional flow acquisition and blood pool contrast agents. With the modern armamentarium of cardiac MRI methods, the goal of compact yet comprehensive exams in children can now be realized.

Baseline 4D Flow-Derived in vivo Hemodynamic Parameters Stratify Descending Aortic Dissection Patients With Enlarging Aortas

Purpose

The purpose of our study was to assess the value of true lumen and false lumen hemodynamics compared to aortic morphological measurements for predicting adverse-aorta related outcomes (AARO) and aortic growth in patients with type B aortic dissection (TBAD).

Materials and Methods

Using an IRB approved protocol, we retrospectively identified patients with descending aorta (DAo) dissection at a large tertiary center. Inclusion criteria includes known TBAD with ≥ 6 months of clinical follow-up after initial presentation for TBAD or after ascending aorta intervention for patients with repaired type A dissection with residual type B aortic dissection (rTAAD). Patients with prior descending aorta intervention were excluded. The FL and TL of each patient were manually segmented from 4D flow MRI data, and 3D parametric maps of aortic hemodynamics were generated. Groups were divided based on (1) presence vs. absence of AARO and (2) growth rate ≥ vs. < 3 mm/year. True and false lumen kinetic energy (KE), stasis, peak velocity (PV), reverse/forward flow (RF/FF), FL to TL KE ratio, as well as index aortic diameter were compared between groups using the Mann–Whitney U or independent t-test.

Results

A total of n = 51 patients (age: 58.4 ± 15.0 years, M/F: 31/20) were included for analysis of AARO. This group contained n = 26 patients with TBAD and n = 25 patients with rTAAD. In the overall cohort, AARO patients had larger baseline diameters, lower FL-RF, FL stasis, TL-KE, TL-FF and TL-PV. Among patients with de novo TBAD, those with AAROs had larger baseline diameter, lower FL stasis and TL-PV. In both the overall cohort and in the subgroup of de novo TBAD, subjects with aortic growth ≥ 3mm/year, patients had a higher KE ratio.

Conclusion

Our study suggests that 4D flow MRI is a promising tool for TBAD evaluation that can provide information beyond traditional MRA or CTA. 4D flow has the potential to become an integral aspect of TBAD work-up, as hemodynamic assessment may allow earlier identification of at-risk patients who could benefit from earlier intervention.

Impact of extra-anatomical bypass on coarctation fluid dynamics using patient-specific lumped parameter and Lattice Boltzmann modeling

Accurate hemodynamic analysis is not only crucial for successful diagnosis of coarctation of the aorta (COA), but intervention decisions also rely on the hemodynamics assessment in both pre and post intervention states to minimize patient risks. Despite ongoing advances in surgical techniques for COA treatments, the impacts of extra-anatomic bypass grafting, a surgical technique to treat COA, on the aorta are not always benign. Our objective was to investigate the impact of bypass grafting on aortic hemodynamics. We investigated the impact of bypass grafting on aortic hemodynamics using a patient-specific computational-mechanics framework in three patients with COA who underwent bypass grafting. Our results describe that bypass grafting improved some hemodynamic metrics while worsened the others: (1) Doppler pressure gradient improved (decreased) in all patients; (2) Bypass graft did not reduce the flow rate substantially through the COA; (3) Systemic arterial compliance increased in patients #1 and 3 and didn't change (improve) in patient 3; (4) Hypertension got worse in all patients; (5) The flow velocity magnitude improved (reduced) in patient 2 and 3 but did not improve significantly in patient 1; (6) There were elevated velocity magnitude, persistence of vortical flow structure, elevated turbulence characteristics, and elevated wall shear stress at the bypass graft junctions in all patients. We concluded that bypass graft may lead to pseudoaneurysm formation and potential aortic rupture as well as intimal hyperplasia due to the persistent abnormal and irregular aortic hemodynamics in some patients. Moreover, post-intervention, exposures of endothelial cells to high shear stress may lead to arterial remodeling, aneurysm, and rupture.

Quantitative normal values of helical flow, flow jets and wall shear stress of healthy volunteers in the ascending aorta

OBJECTIVES

4D flow MRI enables quantitative assessment of helical flow. We sought to generate normal values and elucidate changes of helical flow (duration, volume, length, velocities and rotational direction) and flow jet (displacement, flow angle) as well as wall shear stress (WSS).

METHODS

We assessed the temporal helical existence (TH_EX), maximum helical volume (HV_max), accumulated helical volume (HV_acc), accumulated helical volume length (HVL_acc), maximum forward velocity (maxV_for), maximum circumferential velocity (maxV_circ), rotational direction (RD) and maximum wall shear stress (WSS) as reported elsewhere using the software tool Bloodline in 86 healthy volunteers (46 females, mean age 41 ± 13 years).

RESULTS

WSS decreased by 42.1% and maxV_for by 55.7% across age. There was no link between age and gender regarding the other parameters.

CONCLUSION

This study provides age-dependent normal values regarding WSS and maxV_for and age- and gender-independent normal values regarding TH_EX, HV_max, HV_acc, HVL_acc, RD and _maxV_circ.

KEY POINTS

• 4D flow provides numerous new parameters; therefore, normal values are mandatory. • Wall shear stress decreases over age. • Maximum helical forward velocity decreases over age.

Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications

Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.

Clinical Application of 4D Flow MR Imaging to Pulmonary Hypertension

Pulmonary hypertension (PH) is characterized by elevated pulmonary arterial pressure (PAP). Although right-heart catheterization is the gold standard method for the diagnosis of PH by definition, various less-invasive imaging tests have been used for screening, detection of underlying diseases-causing PH, and monitoring of diseases. Among them, 4D flow MRI is an emerging and unique imaging test that allows for comprehensive visualization of blood flow in the right heart and proximal pulmonary arteries. The characteristic blood flow pattern observed in patients with PH is vortical flow formation in the main pulmonary artery. Recent studies have proposed the use of these findings to determine not only the presence of PH but also estimate the mean PAP. Other applications of 4D flow MRI for PH include measurement of wall shear stress, helicity, and 3D flow balance in the pulmonary arteries. It is worth noting that 4D flow has also the potential for longitudinal follow-ups. In this review, the clinical definition of PH, summary of conventional imaging tests, characteristics of pulmonary arterial flow as shown by 4D flow MRI, and clinical application of 4D flow MRI in the management of patients with PH will be discussed.

Interstitial Fluid Behavior and Diseases

Living things comprise a typical hierarchical and porous medium, and their most fundamental logical architectures are interstitial structures encapsulating parenchymal structures. The recent discovery of the efficient transport mechanisms of interstitial streams has provided a new understanding of these complex activities. The substance transport of interstitial streams follows mesoscopic fluid behavior dynamics, which is intimately associated with material transfer in nanoconfined spaces and a unique signal transmission. Accordingly, the evaluation of interstitial stream transport behavior at the mesoscopic scale is essential. In this review, recent advances in physical and chemical properties, the substance transport model, and the characterization methods of interstitial streams at the mesoscopic scale, as well as the relationships between interstitial streams and disease are summarized. Interstitial stream transport can be used as a basis to fully mine hierarchal behavior in images to expand imaging behavior into an omics field. By starting from the perspective of soft matter, a new understanding can be gained of health and disease and quantitative physical markers for research, clinical diagnosis, and treatment can be provided, as well as prognosis evaluation in complex diseases such as cancer and Alzheimer's disease. This will provide a foundation for the development of medicine of soft matter.

Right Ventricular Flow Vorticity Relationships With Biventricular Shape in Adult Tetralogy of Fallot

Remodeling in adults with repaired tetralogy of Fallot (rToF) may occur due to chronic pulmonary regurgitation, but may also be related to altered flow patterns, including vortices. We aimed to correlate and quantify relationships between vorticity and ventricular shape derived from atlas-based analysis of biventricular shape. Adult rToF (n = 12) patients underwent 4D flow and cine MRI imaging. Vorticity in the RV was computed after noise reduction using a neural network. A biventricular shape atlas built from 95 rToF patients was used to derive principal component modes, which were associated with vorticity and pulmonary regurgitant volume (PRV) using univariate and multivariate linear regression. Univariate analysis showed that indexed PRV correlated with 3 modes (r = −0.55,−0.50, and 0.6, all p &lt; 0.05) associated with RV dilatation and an increase in basal bulging, apical bulging and tricuspid annulus tilting with more severe regurgitation, as well as a smaller LV and paradoxical movement of the septum. RV outflow and inflow vorticity were also correlated with these modes. However, total vorticity over the whole RV was correlated with two different modes (r = −0.62,−0.69, both p &lt; 0.05). Higher vorticity was associated with both RV and LV shape changes including longer ventricular length, a larger bulge beside the tricuspid valve, and distinct tricuspid tilting. RV flow vorticity was associated with changes in biventricular geometry, distinct from associations with PRV. Flow vorticity may provide additional mechanistic information in rToF remodeling. Both LV and RV shapes are important in rToF RV flow patterns.

In-vitro and In-Vivo Assessment of 4D Flow MRI Reynolds Stress Mapping for Pulsatile Blood Flow

Imaging hemodynamics play an important role in the diagnosis of abnormal blood flow due to vascular and valvular diseases as well as in monitoring the recovery of normal blood flow after surgical or interventional treatment. Recently, characterization of turbulent blood flow using 4D flow magnetic resonance imaging (MRI) has been demonstrated by utilizing the changes in signal magnitude depending on intravoxel spin distribution. The imaging sequence was extended with a six-directional icosahedral (ICOSA6) flow-encoding to characterize all elements of the Reynolds stress tensor (RST) in turbulent blood flow. In the present study, we aimed to demonstrate the feasibility of full RST analysis using ICOSA6 4D flow MRI under physiological conditions. First, the turbulence analysis was performed through in vitro experiments with a physiological pulsatile flow condition. Second, a total of 12 normal subjects and one patient with severe aortic stenosis were analyzed using the same sequence. The in-vitro study showed that total turbulent kinetic energy (TKE) was less affected by the signal-to-noise ratio (SNR), however, maximum principal turbulence shear stress (MPTSS) and total turbulence production (TP) had a noise-induced bias. Smaller degree of the bias was observed for TP compared to MPTSS. In-vivo study showed that the subject-variability on turbulence quantification was relatively low for the consistent scan protocol. The in vivo demonstration of the stenosis patient showed that the turbulence analysis could clearly distinguish the difference in all turbulence parameters as they were at least an order of magnitude larger than those from the normal subjects.

Divergence-Free Constrained Phase Unwrapping and Denoising for 4D Flow MRI Using Weighted Least-Squares

A novel divergence-free constrained phase unwrapping method was proposed and evaluated for 4D flow MRI. The unwrapped phase field was obtained by integrating the phase variations estimated from the wrapped phase data using weighted least-squares. The divergence-free constraint for incompressible blood flow was incorporated to regulate and denoise the resulting phase field. The proposed method was tested on synthetic phase data of left ventricular flow and in vitro 4D flow measurement of Poiseuille flow. The method was additionally applied to in vivo 4D flow measurements in the thoracic aorta from 30 human subjects. The performance of the proposed method was compared to the state-of-the-art 4D single-step Laplacian algorithm. The synthetic phase data were completely unwrapped by the proposed method for all the cases with velocity encoding (venc) as low as 20% of the maximum velocity and signal-to-noise ratio as low as 5. The in vitro Poiseuille flow data were completely unwrapped with a 60% increase in the velocity-to-noise ratio. For the in-vivo aortic datasets with venc ratio less than 0.4, the proposed method significantly improved the success rate by as much as 40% and reduced the velocity error levels by a factor of 10 compared to the state-of-the-art method. The divergence-free constrained method exhibits reliability and robustness on phase unwrapping and shows improved accuracy of velocity and hemodynamic quantities by unwrapping the low-venc 4D flow MRI data.

An optimal control approach to determine resistance-type boundary conditions from in-vivo data for cardiovascular simulations

The choice of appropriate boundary conditions is a fundamental step in computational fluid dynamics (CFD) simulations of the cardiovascular system. Boundary conditions, in fact, highly affect the computed pressure and flow rates, and consequently haemodynamic indicators such as wall shear stress (WSS), which are of clinical interest. Devising automated procedures for the selection of boundary conditions is vital to achieve repeatable simulations. However, the most common techniques do not automatically assimilate patient-specific data, relying instead on expensive and time-consuming manual tuning procedures. In this work, we propose a technique for the automated estimation of outlet boundary conditions based on optimal control. The values of resistive boundary conditions are set as control variables and optimized to match available patient-specific data. Experimental results on four aortic arches demonstrate that the proposed framework can assimilate 4D-Flow MRI data more accurately than two other common techniques based on Murray's law and Ohm's law.

Alterations in Intracardiac Flow Patterns Affect Mitral Leaflets Dynamics in a Model of Ischemic Mitral Regurgitation

PURPOSE

This study was to evaluate the effects of ischemic mitral regurgitation (IMR) on vortex formation and leaflet dynamics using an established porcine infarct model of IMR.

METHODS

Using direct coronary ligation, five animals were subjected to a posterolateral myocardial infarction (MI) followed by an MRI at 12-weeks post MI. MR imaging consisted of 4D time-resolved left ventricular (LV) flow, full coverage 2D LV cine, and high resolution 2D cine of mitral valve dynamics. Five additional naïve animals underwent identical imaging protocols to serve as controls. Image analysis was performed to obtain mitral transvalvular flows as well as LV volumes throughout the cardiac cycle. In addition, anterior to posterior mid-leaflet tip distances were measured throughout the cardiac cycle for determination of temporal leaflet dynamics.

RESULTS

It was found IMR caused asymmetric vortex ring formation with the anterior vortex having a lower vorticity relative to its posterior counterpart. In contrast, normal ventricles create symmetric and tightly curled vortices in the basal chamber just underneath the mitral leaflets which conserve kinetic energy and aid in effective ejection. IMR animals were also evaluated for leaflet separation and were found to have a greater leaflet opening and achieved peak vorticity and peak leaflet opening later than control animals.

CONCLUSION

In conclusion, this study shows the effects that altered vortex formation, due to IMR, can have on ventricular filling and leaflet dynamics. These findings have important implications for understanding blood flow through the dilated heart and how ring annuloplasty and volume reduction interventions may influence mitral valve dynamics.

Patient-Specific 3D Bioprinted Models of Developing Human Heart

The heart is the first organ to develop in the human embryo through a series of complex chronological processes, many of which critically rely on the interplay between cells and the dynamic microenvironment. Tight spatiotemporal regulation of these interactions is key in heart development and diseases. Due to suboptimal experimental models, however, little is known about the role of microenvironmental cues in the heart development. This study investigates the use of 3D bioprinting and perfusion bioreactor technologies to create bioartificial constructs that can serve as high-fidelity models of the developing human heart. Bioprinted hydrogel-based, anatomically accurate models of the human embryonic heart tube (e-HT, day 22) and fetal left ventricle (f-LV, week 33) are perfused and analyzed both computationally and experimentally using ultrasound and magnetic resonance imaging. Results demonstrate comparable flow hemodynamic patterns within the 3D space. We demonstrate endothelial cell growth and function within the bioprinted e-HT and f-LV constructs, which varied significantly in varying cardiac geometries and flow. This study introduces the first generation of anatomically accurate, 3D functional models of developing human heart. This platform enables precise tuning of microenvironmental factors, such as flow and geometry, thus allowing the study of normal developmental processes and underlying diseases.

Collaterals in portal hypertension: anatomy and clinical relevance

Portal hypertension is a key pathophysiology of chronic liver diseases typified with cirrhosis or noncirrhotic portal hypertension. The development of collateral vessels is a characteristic feature of impaired portal hemodynamics. The paraumbilical vein (PUV), left gastric vein (LGV), posterior gastric vein (PGV), short gastric vein (SGV), splenorenal shunt (SRS), and inferior mesenteric vein (IMV) are major collaterals, and there are some rare collaterals. The degree and hemodynamics of collateral may affect the portal venous circulation and may compensate for the balance between inflow and outflow volume of the liver. Additionally, the development of collateral shows a relation with the liver function reserve and clinical manifestations such as esophageal varices (EV), gastric varices, rectal varices and the other ectopic varices, hepatic encephalopathy, and prognosis. Furthermore, there may be an interrelationship in the development between different collaterals, showing additional influences on the clinical presentations. Thus, the assessment of collaterals may enhance the understanding of the underlying pathophysiology of the condition of patients with portal hypertension. This review article concluded that each collateral has a specific function depending on the anatomy and hemodynamics and is linked with the relative clinical presentation in patients with portal hypertension. Imaging modalities may be essential for the detection, grading and evaluation of the role of collaterals and may help to understand the pathophysiology of the patient condition. Further investigation in a large-scale study would elucidate the basic and clinical significance of collaterals in patients with portal hypertension and may provide information on how to manage them to improve the prognosis as well as quality of life.

Abnormal pulmonary flow is associated with impaired right ventricular coupling in patients with COPD

Cor Pulmonale or right ventricular (RV) dysfunction due to pulmonary disease is an expected complication of COPD resulting primarily from increased afterload mediated by chronic alveolar hypoxemia and resulting hypoxic pulmonary vasoconstriction. Early detection of elevated RV afterload has been previously demonstrated by visualization of abnormal flow patterns in the proximal pulmonary arteries. Prior analysis of helicity in the pulmonary arteries in pulmonary hypertension patients has demonstrated a strong association between helicity and increased RV afterload. However, these flow hemodynamics have yet to be fully explored in patients with COPD. We hypothesized that patients with COPD will have abnormal pulmonary flow as evaluated by 4D-Flow MRI and associated with RV function and pulmonary arterial stiffness. Patients with COPD (n = 15) (65 years ± 6) and controls (n = 10) (58 years ± 9) underwent 4D-Flow MRI to calculate helicity. The helicity was calculated in the main pulmonary artery (MPA) and along the RV outflow tract (RVOT)-MPA axis. Main pulmonary arterial stiffness was measured using the relative area change (RAC). We found COPD patients had decreased helicity relative to healthy controls in the MPA (19.4 ± 7.8vs 32.8 ± 15.9, P = 0.007) and reduced helicity along the RVOT-MPA axis (33.2 ± 9.0 vs 43.5 ± 8.3, P = 0.010). Our investigation indicates a strong association between helicity along the MPA-RV outflow tract axis and RV function and suggests that 4D-Flow MRI might be a sensitive tool in evaluating RV-pulmonary arterial coupling in COPD.

False lumen pressure estimation in type B aortic dissection using 4D flow cardiovascular magnetic resonance: comparisons with aortic growth

BACKGROUND

Chronic type B aortic dissection (TBAD) is associated with poor long-term outcome, and accurate risk stratification tools remain lacking. Pressurization of the false lumen (FL) has been recognized as central in promoting aortic growth. Several surrogate imaging-based metrics have been proposed to assess FL hemodynamics; however, their relationship to enlarging aortic dimensions remains unclear. We investigated the association between aortic growth and three cardiovascular magnetic resonance (CMR)-derived metrics of FL pressurization: false lumen ejection fraction (FLEF), maximum systolic deceleration rate (MSDR), and FL relative pressure (FL ΔPmax).

METHODS

CMR/CMR angiography was performed in 12 patients with chronic dissection of the descending thoracoabdominal aorta, including contrast-enhanced CMR angiography and time-resolved three-dimensional phase-contrast CMR (4D Flow CMR). Aortic growth rate was calculated as the change in maximal aortic diameter between baseline and follow-up imaging studies over the time interval, with patients categorized as having either 'stable' (< 3 mm/year) or 'enlarging' (≥ 3 mm/year) growth. Three metrics relating to FL pressurization were defined as: (1) FLEF: the ratio between retrograde and antegrade flow at the TBAD entry tear, (2) MSDR: the absolute difference between maximum and minimum systolic acceleration in the proximal FL, and (3) FL ΔPmax: the difference in absolute pressure between aortic root and distal FL.

RESULTS

FLEF was higher in enlarging TBAD (49.0 ± 17.9% vs. 10.0 ± 11.9%, p = 0.002), whereas FL ΔPmax was lower (32.2 ± 10.8 vs. 57.2 ± 12.5 mmHg/m, p = 0.017). MSDR and conventional anatomic variables did not differ significantly between groups. FLEF showed positive (r = 0.78, p = 0.003) correlation with aortic growth rate whereas FL ΔPmax showed negative correlation (r = - 0.64, p = 0.026). FLEF and FL ΔPmax remained as independent predictors of aortic growth rate after adjusting for baseline aortic diameter.

CONCLUSION

Comparative analysis of three 4D flow CMR metrics of TBAD FL pressurization demonstrated that those that focusing on retrograde flow (FLEF) and relative pressure (FL ΔPmax) independently correlated with growth and differentiated patients with enlarging and stable descending aortic dissections. These results emphasize the highly variable nature of aortic hemodynamics in TBAD patients, and suggest that 4D Flow CMR derived metrics of FL pressurization may be useful to separate patients at highest and lowest risk for progressive aortic growth and complications.

In vitro evaluation of cerebrospinal fluid velocity measurement in type I Chiari malformation: repeatability, reproducibility, and agreement using 2D phase contrast and 4D flow MRI

Background

Phase contrast magnetic resonance imaging, PC MRI, is a valuable tool allowing for non-invasive quantification of CSF dynamics, but has lacked adoption in clinical practice for Chiari malformation diagnostics. To improve these diagnostic practices, a better understanding of PC MRI based measurement agreement, repeatability, and reproducibility of CSF dynamics is needed.

Methods

An anatomically realistic in vitro subject specific model of a Chiari malformation patient was scanned three times at five different scanning centers using 2D PC MRI and 4D Flow techniques to quantify intra-scanner repeatability, inter-scanner reproducibility, and agreement between imaging modalities. Peak systolic CSF velocities were measured at nine axial planes using 2D PC MRI, which were then compared to 4D Flow peak systolic velocity measurements extracted at those exact axial positions along the model.

Results

Comparison of measurement results showed good overall agreement of CSF velocity detection between 2D PC MRI and 4D Flow (p = 0.86), fair intra-scanner repeatability (confidence intervals ± 1.5 cm/s), and poor inter-scanner reproducibility. On average, 4D Flow measurements had a larger variability than 2D PC MRI measurements (standard deviations 1.83 and 1.04 cm/s, respectively).

Conclusion

Agreement, repeatability, and reproducibility of 2D PC MRI and 4D Flow detection of peak CSF velocities was quantified using a patient-specific in vitro model of Chiari malformation. In combination, the greatest factor leading to measurement inconsistency was determined to be a lack of reproducibility between different MRI centers. Overall, these findings may help lead to better understanding for application of 2D PC MRI and 4D Flow techniques as diagnostic tools for CSF dynamics quantification in Chiari malformation and related diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12987-021-00246-3.

An LDV based method to quantify the error of PC-MRI derived Wall Shear Stress measurement

Wall Shear Stress (WSS) has been demonstrated to be a biomarker of the development of atherosclerosis. In vivo assessment of WSS is still challenging, but 4D Flow MRI represents a promising tool to provide 3D velocity data from which WSS can be calculated. In this study, a system based on Laser Doppler Velocimetry (LDV) was developed to validate new improvements of 4D Flow MRI acquisitions and derived WSS computing. A hydraulic circuit was manufactured to allow both 4D Flow MRI and LDV velocity measurements. WSS profiles were calculated with one 2D and one 3D method. Results indicated an excellent agreement between MRI and LDV velocity data, and thus the set-up enabled the evaluation of the improved performances of 3D with respect to the 2D-WSS computation method. To provide a concrete example of the efficacy of this method, the influence of the spatial resolution of MRI data on derived 3D-WSS profiles was investigated. This investigation showed that, with acquisition times compatible with standard clinical conditions, a refined MRI resolution does not improve WSS assessment, if the impact of noise is unreduced. This study represents a reliable basis to validate with LDV WSS calculation methods based on 4D Flow MRI.

Partial Anomalous Pulmonary Venous Connection Repair: Customized Approach and Outcomes

Alternative options for the correction of partial anomalous pulmonary venous connection (PAPVC) have been proposed. Each can be associated with variable risk for dysrhythmias, caval or pulmonary venous (PV) obstruction. A selective customized strategy to address PAPVC taking into account atrial shunt (AS) and growth potential was pursued. Between September 2014 and August 2018 21 PAPVC patients were identified. Two levels of reference determined the chosen repair strategy; azygous vein (AzV) and cavoatrial junction (CAJ). Six (Group-A) with PAPVC entering SVC cephalad to AV underwent a combined in-situ cavoatrial autologous reconstruction with atrial appendage advancement flap (CARAF). PAPVC entering caudally to AzV (Group-B) underwent alternative repair (caval division/Warden-type or intraatrial rerouting) (n = 15). Age was 8.3 (IQR:4.2-18.5) years for Group-A (vs 11.9; IQR:8.8-34.7 in Group-B) (p = 0.07). In Group-A 5(83%) had AS (vs 12[80%] Group-B; p = 0.9). None had left SVC in Group-A (vs 1 in Group-B; p = 0.9). Preoperative advanced imaging and echocardiographic hemodynamic evaluation was undertaken. Follow-up was complete (median 2.9; IQR:1.2-3.7 years). Freedom from atrial dysrhythmias, caval or PV obstruction was assessed. There were no early or late deaths. ICU and hospital length of stay were 1.8 ± 1.1 and 3.2 ± 0.5 days, respectively. No atrial dysrhythmias occurred postoperatively in Group-A (vs 1 in Group-B; p = 0.9). No permanent pacemaker was implanted. All patients remained in normal sinus rhythm. There were no early or late caval/PV obstruction. A customized approach reserves the advantages of each technique tailored to patient's needs. Expanding surgical capacity with favorable outlook for all PAPVC variations, irrespective of association with AS, can maximize efficiency and reproducibility paired with the lowest morbidity.

Multiscale modeling of blood flow to assess neurological complications in patients supported by venoarterial extracorporeal membrane oxygenation

Computational blood flow models in large arteries elucidate valuable relationships between cardiovascular diseases and hemodynamics, leading to improvements in treatment planning and clinical decision making. One such application with potential to benefit from simulation is venoarterial extracorporeal membrane oxygenation (VA-ECMO), a support system for patients with cardiopulmonary failure. VA-ECMO patients develop high rates of neurological complications, partially due to abnormal blood flow throughout the vasculature from the VA-ECMO system. To better understand these hemodynamic changes, it is important to resolve complex local flow parameters derived from three-dimensional (3D) fluid dynamics while also capturing the impact of VA-ECMO support throughout the systemic arterial system. As high-resolution 3D simulations of the arterial network remain computationally expensive and intractable for large studies, a validated, multiscale model is needed to compute both global effects and high-fidelity local hemodynamics. In this work, we developed and demonstrated a framework to model hemodynamics in VA-ECMO patients using coupled 3D and one-dimensional (1D) models (1D→3D). We demonstrated the ability of these multiscale models to simulate complex flow patterns in specific regions of interest while capturing bulk flow throughout the systemic arterial system. We compared 1D, 3D, and 1D→3D coupled models and found that multiscale models were able to sufficiently capture both global and local hemodynamics in the cerebral arteries and aorta in VA-ECMO patients. This study is the first to develop and compare 1D, 3D, and 1D→ 3D coupled models on the larger arterial system scale in VA-ECMO patients, with potential use for other large scale applications.

Fetal whole heart blood flow imaging using 4D cine MRI

Prenatal detection of congenital heart disease facilitates the opportunity for potentially life-saving care immediately after the baby is born. Echocardiography is routinely used for screening of morphological malformations, but functional measurements of blood flow are scarcely used in fetal echocardiography due to technical assumptions and issues of reliability. Magnetic resonance imaging (MRI) is readily used for quantification of abnormal blood flow in adult hearts, however, existing in utero approaches are compromised by spontaneous fetal motion. Here, we present and validate a novel method of MRI velocity-encoding combined with a motion-robust reconstruction framework for four-dimensional visualization and quantification of blood flow in the human fetal heart and major vessels. We demonstrate simultaneous 4D visualization of the anatomy and circulation, which we use to quantify flow rates through various major vessels. The framework introduced here could enable new clinical opportunities for assessment of the fetal cardiovascular system in both health and disease.

In Vitro Assessment of Flow Variability in an Intracranial Aneurysm Model Using 4D Flow MRI and Tomographic PIV

Aneurysm rupture has been suggested to be related to aneurysm geometry, morphology, and complex flow activity; therefore, understanding aneurysm-specific hemodynamics is crucial. 4D Flow MRI has been shown to be a feasible tool for assessing hemodynamics in intracranial aneurysms with high spatial resolution. However, it requires averaging over multiple heartbeats and cannot account for cycle-to-cycle hemodynamics variations. This study aimed to assess cycle-to-cycle flow dynamics variations in a patient-specific intracranial aneurysm model using tomographic particle image velocimetry (tomo-PIV) at a high image rate under pulsatile flow conditions. Time-resolved and time-averaged velocity flow fields within the aneurysm sac and estimations of wall shear stress (WSS) were compared with those from 4D Flow MRI. A one-way ANOVA showed a significant difference between cardiac cycles (p value < 0.0001); however, differences were not significant after PIV temporal and spatial resolution was matched to that of MRI (p value 0.9727). This comparison showed the spatial resolution to be the main contributor to assess cycle-to-cycle variability. Furthermore, the comparison with 4D Flow MRI between velocity components, streamlines, and estimated WSS showed good qualitative and quantitative agreement. This study showed the feasibility of patient-specific in-vitro experiments using tomo-PIV to assess 4D Flow MRI with high repeatability in the measurements.

4D flow imaging of the thoracic aorta: is there an added clinical value?

Four-dimensional (4D) flow MRI has emerged as a powerful non-invasive technique in cardiovascular imaging, enabling to analyse in vivo complex flow dynamics models by quantifying flow parameters and derived features. Deep knowledge of aortic flow dynamics is fundamental to better understand how abnormal flow patterns may promote or worsen vascular diseases. In the perspective of an increasingly personalized and preventive medicine, growing interest is focused on identifying those quantitative functional features which are early predictive markers of pathological evolution. The thoracic aorta and its spectrum of diseases, as the first area of application and development of 4D flow MRI and supported by an extensive experimental validation, represents the ideal model to introduce this technique into daily clinical practice. The purpose of this review is to describe the impact of 4D flow MRI in the assessment of the thoracic aorta and its most common affecting diseases, providing an overview of the actual clinical applications and describing the potential role of derived advanced hemodynamic measures in tailoring follow-up and treatment.

Detecting Aortic Valve-Induced Abnormal Flow with Seismocardiography and Cardiac MRI

Cardiac MRI (CMR) techniques offer non-invasive visualizations of cardiac morphology and function. However, imaging can be time-consuming and complex. Seismocardiography (SCG) measures physical vibrations transmitted through the chest from the beating heart and pulsatile blood flow. SCG signals can be acquired quickly and easily, with inexpensive electronics. This study investigates relationships between CMR metrics of function and SCG signal features. Same-day CMR and SCG data were collected from 28 healthy adults and 6 subjects with aortic valve disease history. Correlation testing and statistical median/decile calculations were performed with data from the healthy cohort. MR-quantified flow and function parameters in the healthy cohort correlated with particular SCG energy levels, such as peak aortic velocity with low-frequency SCG (coefficient 0.43, significance 0.02) and peak flow with high-frequency SCG (coefficient 0.40, significance 0.03). Valve disease-induced flow abnormalities in patients were visualized with MRI, and corresponding abnormalities in SCG signals were identified. This investigation found significant cross-modality correlations in cardiac function metrics and SCG signals features from healthy subjects. Additionally, through comparison to normative ranges from healthy subjects, it observed correspondences between pathological flow and abnormal SCG. This may support development of an easy clinical test used to identify potential aortic flow abnormalities.

Automated Quantitative Extraction and Analysis of 4D flow Patterns in the Ascending Aorta: An intraindividual comparison at 1.5 T and 3 T

4D flow MRI enables quantitative assessment of helical flow. Current methods are susceptible to noise. To evaluate helical flow patterns in healthy volunteers and patients with bicuspid aortic valves (BAV) at 1.5 T and 3 T using pressure-based helix-extraction in 4D flow MRI. Two intraindividual 4D flow MRI examinations were performed at 1.5 T and 3 T in ten healthy volunteers (5 females, 32 ± 3 years) and 2 patients with BAV using different acceleration techniques (kt-GRAPPA and centra). Several new quantitative parameters for the evaluation of volumes [ml], lengths [mm] as well as temporal parameters [ms] of helical flow were introduced and analyzed using the software tool Bloodline. We found good correlations between measurements in volunteers at 1.5 T and 3 T regarding helical flow volumes (R = 0.98) and temporal existence (R = 0.99) of helices in the ascending aorta. Furthermore, we found significantly larger (11.7 vs. 77.6 ml) and longer lasting (317 vs. 769 ms) helices in patients with BAV than in volunteers. The assessed parameters do not depend on the magnetic field strength used for the acquisition. The technique of pressure-based extraction of 4D flow MRI pattern is suitable for differentiation of normal and pathological flow.

Altered 4-D magnetic resonance imaging flow characteristics in complex congenital aortic arch repair

BACKGROUND

Interrupted aortic arch (IAA) is a rare but severe congenital abnormality often associated with bicuspid aortic valve (BAV). Complex re-interventions are often needed despite surgical advances, but the impact of aortic hemodynamics in repaired patients is unknown.

OBJECTIVE

Investigate effect of IAA repairs on aortic hemodynamics, wall shear stress and flow derangements via 4-D flow MRI.

MATERIALS AND METHODS

We retrospectively analyzed age- and gender-matched cohorts (IAA [n=6], BAV alone [n=6], controls [n=6]) undergoing cardiac MRI including 4-D flow. Aortic dimensions were measured from standard MR angiography. We quantified peak systolic velocities, regurgitant fractions and wall shear stress in the ascending aorta (AAo), transverse arch and descending aorta (DAo) from 4-D flow, and we graded helix/vortex flow patterns from 3-D blood flow visualization.

RESULTS

Children and young adults with IAA had a wide range of arch dimensions, peak systolic velocities, regurgitant fractions and flow grades. Peak transverse arch systolic velocities were higher in patients with IAA versus controls (P=0.02). Flow derangements in the AAo were found in patients with IAA (median grade=2, 5/6 patients, P=0.04) and BAV (median grade=3, 5/6 patients, P=0.03) versus controls. Flow derangements in the DAo were only seen in patients with IAA (median grade=1, 5/6 patients, P=0.04), and 5/6 people with IAA had helical flow in head and neck vessels. Wall shear stress was increased in people with IAA along the superior transverse arch and proximal DAo versus controls (P=0.02).

CONCLUSION

Complex congenital aortic arch repairs can change aortic hemodynamics. Associated cardiac defects can further alter findings. Studies are warranted to investigate clinical implications in larger cohorts.