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.

Echocardiographic haemodynamic monitoring in the context of HeartMate 3™ therapy: a systematic review

AIMS

While echocardiography remains essential within haemodynamic monitoring of durable mechanical circulatory support, previous echocardiographic guidelines are missing scientific evidence for the novel HeartMate 3™ (HM3) system. Accordingly, this review aims to summarize available echocardiographic evidence including HM3.

METHODS AND RESULTS

This systematic review adhered to the PRISMA 2020 guidelines. Searches were conducted during August 2023 across PubMed, Embase, and Google Scholar using specific echocardiographic terms combined with system identifiers. Study quality was assessed using the Newcastle-Ottawa Scale (NOS) for cohort studies and Critical Appraisal Instrument (PCAI) for cross-sectional studies. Nine studies met the inclusion criteria, of which eight cohort studies and one cross-sectional study. Aortic regurgitation (AR) prevalence at approximately 12 months of support exhibited heterogenicity (33.5% (Δ 33%)) in a limited number of studies (n = 3). Several studies (n = 5) demonstrated an increasing prevalence and severity of AR during HM3 support, generating moderate to high level of evidence. One AR study showed a higher cumulative incidence of death and heart failure (HF) readmission compared with those without significant AR, hazard ratio 3.42 (95% CI 1.48-8.76). A second study showed that a worsening AR group had significantly lower survival-free from HF readmission (59% vs. 89%, P = 0.023) with a hazard ratio of 5.18 (95% CI 1.07-25.0), while a third study did not reveal any differences in cardiac-related hospitalizations in the 12 months follow-up or non-cardiac-related hospitalization. Mitral regurgitation (MR) prevalence at approximately 12 months of support exhibited good consistency 15.0% (Δ 0.8%) in both included studies, which did not reveal any significant pattern of changing prevalence over time. Tricuspid regurgitation (TR) prevalence at approximately 12 months of support exhibited fair consistency 28.5% (Δ 8.3%) in a limited number of studies (n = 2); both studies showed a statistically un-confirmed trend of increased TR prevalence over time. The evidence of general prevalence of right ventricular dysfunction (RVD) was insufficient due to lack of studies.

CONCLUSIONS

There are few methodologically consistent studies with focus on long-term haemodynamic effects. Aortic regurgitation still seems to be a prevalent and potentially significant finding. The available evidence concerning right heart function is limited despite clinical relevance and potential prognostic value. Potential interventricular and haemodynamic interplay are identified as a white field for future research.

Regional and Global Aortic Pulse Wave Velocity in Patients with Abdominal Aortic Aneurysm

OBJECTIVE

Abdominal aortic aneurysm (AAA) is commonly defined as localised aortic dilatation with a diameter > 30 mm. The pathophysiology of AAA includes chronic inflammation and enzymatic degradation of elastin, possibly increasing aortic wall stiffness and pulse wave velocity (PWV). Whether aortic stiffness is more prominent in the abdominal aorta at the aneurysm site is not elucidated. The aim of this study was to evaluate global and regional aortic PWV in patients with AAA.

METHODS

Experimental study of local PWV in the thoracic descending and abdominal aorta in patients with AAA and matched controls. The study cohort comprised 25 patients with an AAA > 30 mm (range 36 - 70 mm, all male, age range 65 - 76 years) and 27 age and sex matched controls free of AAA. PWV was measured with applanation tonometry (carotid-femoral PWV, cfPWV) as well as a 4D flow MRI technique, assessing regional aortic PWV. Blood pressure and anthropometrics were measured.

RESULTS

Global aortic PWV was greater in men with an AAA than controls, both by MRI (AAA 8.9 ± 2.4 m/s vs. controls 7.1 ± 1.5 m/s; p = .007) and cfPWV (AAA 11.0 ± 2.1 m/s vs. controls 9.3 ± 2.3 m/s; p = .007). Regionally, PWV was greater in the abdominal aorta in the AAA group (AAA 7.0 ± 1.8 m/s vs. controls 5.8 ± 1.0 m/s; p = .022), but similar in the thoracic descending aorta (AAA 8.7 ± 3.2 m/s vs. controls 8.2 ± 2.4 m/s; p = .59). Furthermore, PWV was positively associated with indices of central adiposity both in men with AAA and controls.

CONCLUSION

PWV is higher in men with AAA compared with matched controls in the abdominal but not the thoracic descending aorta. Furthermore, aortic stiffness was linked with central fat deposition. It remains to be seen whether there is a causal link between AAA and increased regional aortic stiffness.

Comprehensive left atrial flow component analysis reveals abnormal flow patterns in paroxysmal atrial fibrillation

Left atrial (LA) blood flow plays an important role in diseases such as atrial fibrillation (AF) and atrial cardiomyopathy since alterations in the blood flow might lead to thrombus formation and stroke. Using traditional techniques, such as echocardiography, atrial flow velocities can be measured at the pulmonary veins and the mitral valve, but a comprehensive understanding of the three-dimensional atrial flow field is missing. Previously, ventricular flow has been analyzed using flow component analysis, revealing new insights into ventricular flow and function. Thus, the aim of this project was to develop a comprehensive flow component analysis method for the LA and explore its utility in 21 patients with paroxysmal atrial fibrillation compared with a control group of 8 participants. The flow field was derived from time-resolved CT acquired during sinus rhythm using computational fluid dynamics. Flow components were computed from particle tracking. We identified six atrial flow components: conduit, reservoir, delayed ejection, retained inflow, residual volume, and pulmonary vein backflow. It was shown that conduit flow, defined as blood entering and leaving the LA within the same diastolic phase, exists in most subjects. Although the volume of conduit and reservoir is similar in patients with paroxysmal AF in sinus rhythm and controls, the volume of the other components is increased in paroxysmal AF. Comprehensive quantification of LA flow using flow component analysis makes atrial blood flow quantifiable, thus facilitating investigation of mechanisms underlying atrial dysfunction and can increase understanding of atrial blood flow in disease progression and stroke risk.NEW & NOTEWORTHY We developed a new comprehensive approach to atrial blood component analysis that includes both conduit flow and residual volume and compared the flow components of atrial fibrillation (AF) patients in sinus rhythm with controls. Conduit and reservoir flow were similar between the groups, whereas components with longer residence time in the left atrium were increased in the AF group. This could add to the pathophysiological understanding of atrial diseases and possibly clinical management.

Corrigendum: Elevated atrial blood stasis in paroxysmal atrial fibrillation during sinus rhythm: a patient-specific computational fluid dynamics study

[This corrects the article DOI: 10.3389/fcvm.2023.1219021.].

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.

Novel insights into diminished cardiac reserve in non-obstructive hypertrophic cardiomyopathy from four-dimensional flow cardiac magnetic resonance component analysis

AIMS

Hypertrophic cardiomyopathy (HCM) is characterized by hypercontractility and diastolic dysfunction, which alter blood flow haemodynamics and are linked with increased risk of adverse clinical events. Four-dimensional flow cardiac magnetic resonance (4D-flow CMR) enables comprehensive characterization of ventricular blood flow patterns. We characterized flow component changes in non-obstructive HCM and assessed their relationship with phenotypic severity and sudden cardiac death (SCD) risk.

METHODS AND RESULTS

Fifty-one participants (37 non-obstructive HCM and 14 matched controls) underwent 4D-flow CMR. Left-ventricular (LV) end-diastolic volume was separated into four components: direct flow (blood transiting the ventricle within one cycle), retained inflow (blood entering the ventricle and retained for one cycle), delayed ejection flow (retained ventricular blood ejected during systole), and residual volume (ventricular blood retained for >two cycles). Flow component distribution and component end-diastolic kinetic energy/mL were estimated. HCM patients demonstrated greater direct flow proportions compared with controls (47.9 ± 9% vs. 39.4 ± 6%, P = 0.002), with reduction in other components. Direct flow proportions correlated with LV mass index (r = 0.40, P = 0.004), end-diastolic volume index (r = -0.40, P = 0.017), and SCD risk (r = 0.34, P = 0.039). In contrast to controls, in HCM, stroke volume decreased with increasing direct flow proportions, indicating diminished volumetric reserve. There was no difference in component end-diastolic kinetic energy/mL.

CONCLUSION

Non-obstructive HCM possesses a distinctive flow component distribution pattern characterised by greater direct flow proportions, and direct flow-stroke volume uncoupling indicative of diminished cardiac reserve. The correlation of direct flow proportion with phenotypic severity and SCD risk highlight its potential as a novel and sensitive haemodynamic measure of cardiovascular risk in HCM.

Elevated atrial blood stasis in paroxysmal atrial fibrillation during sinus rhythm: a patient-specific computational fluid dynamics study

Introduction

Atrial fibrillation (AF) is associated with an increased risk of stroke, often caused by thrombi that form in the left atrium (LA), and especially in the left atrial appendage (LAA). The underlying mechanism is not fully understood but is thought to be related to stagnant blood flow, which might be present despite sinus rhythm. However, measuring blood flow and stasis in the LAA is challenging due to its small size and low velocities. We aimed to compare the blood flow and stasis in the left atrium of paroxysmal AF patients with controls using computational fluid dynamics (CFD) simulations.

Methods

The CFD simulations were based on time-resolved computed tomography including the patient-specific cardiac motion. The pipeline allowed for analysis of 21 patients with paroxysmal AF and 8 controls. Stasis was estimated by computing the blood residence time.

Results and Discussion

Residence time was elevated in the AF group (p < 0.001). Linear regression analysis revealed that stasis was strongest associated with LA ejection ratio (p < 0.001, R2 = 0.68) and the ratio of LA volume and left ventricular stroke volume (p < 0.001, R2 = 0.81). Stroke risk due to LA thrombi could already be elevated in AF patients during sinus rhythm. In the future, patient specific CFD simulations may add to the assessment of this risk and support diagnosis and treatment.

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.

Enhancing students' understanding of cardiac physiology by using 4D visualization

Difficulties in achieving knowledge about physiology and anatomy of the beating heart highlight the challenges with more traditional pedagogical methods. Recent research regarding anatomy education has mainly focused on digital three-dimensional models. However, these pedagogical improvements may not be entirely applicable to cardiac anatomy and physiology due to the multidimensional complexity with moving anatomy and complex blood flow. The aim of this study was therefore to evaluate whether high quality time-resolved anatomical images combined with realistic blood flow simulations improve the understanding of cardiac structures and function. Three time-resolved datasets were acquired using time-resolved computed tomography and blood flow was computed using Computational Fluid Dynamics. The anatomical and blood flow information was combined and interactively visualized using volume rendering on an advanced stereo projection system. The setup was tested in interactive lectures for medical students. Ninety-seven students participated. Summative assessment of examinations showed significantly improved mean score (18.1 ± 4.5 vs 20.3 ± 4.9, p = 0.002). This improvement was driven by knowledge regarding myocardial hypertrophy and pressure-velocity differences over a stenotic valve. Additionally, a supplementary formative assessment showed significantly more agreeing answers than disagreeing answers (p < 0.001) when the participants subjectively evaluated the contribution of the visualizations to their education and knowledge. In conclusion, the use of simultaneous visualization of time-resolved anatomy data and simulated blood flow improved medical students' results, with a particular effect on understanding of cardiac physiology and these simulations may be useful educational tools for teaching complex anatomical and physiological concepts.

Assessment of transmitral and left atrial appendage flow rate from cardiac 4D-CT

BACKGROUND

Cardiac time-resolved CT (4D-CT) acquisitions provide high quality anatomical images of the heart. However, some cardiac diseases require assessment of blood flow in the heart. Diastolic dysfunction, for instance, is diagnosed by measuring the flow through the mitral valve (MV), while in atrial fibrillation, the flow through the left atrial appendage (LAA) indicates the risk for thrombus formation. Accurate validated techniques to extract this information from 4D-CT have been lacking, however.

METHODS

To measure the flow rate though the MV and the LAA from 4D-CT, we developed a motion tracking algorithm that performs a nonrigid deformation of the surface separating the blood pool from the myocardium. To improve the tracking of the LAA, this region was deformed separately from the left atrium and left ventricle. We compared the CT based flow with 4D flow and short axis MRI data from the same individual in 9 patients.

RESULTS

For the mitral valve flow, good agreement was found for the time span between the early and late diastolic peak flow (bias: <0.1 s). The ventricular stroke volume is similar compared to short-axis MRI (bias 3 ml). There are larger differences in the diastolic peak flow rates, with a larger bias for the early flow rate than the late flow rate. The peak LAA outflow rate measured with both modalities matches well (bias: -6 ml/s).

CONCLUSIONS

Overall, the developed algorithm provides accurate tracking of dynamic cardiac geometries resulting in similar flow rates at the MV and LAA compared to 4D flow MRI.

Corrigendum: A head-to-head comparison of myocardial strain by fast-strain encoding and feature tracking imaging in acute myocardial infarction

[This corrects the article DOI: 10.3389/fcvm.2022.949440.].

Non-invasive cardiovascular magnetic resonance assessment of pressure recovery distance after aortic valve stenosis

BACKGROUND

Decisions in the management of aortic stenosis are based on the peak pressure drop, captured by Doppler echocardiography, whereas gold standard catheterization measurements assess the net pressure drop but are limited by associated risks. The relationship between these two measurements, peak and net pressure drop, is dictated by the pressure recovery along the ascending aorta which is mainly caused by turbulence energy dissipation. Currently, pressure recovery is considered to occur within the first 40-50 mm distally from the aortic valve, albeit there is inconsistency across interventionist centers on where/how to position the catheter to capture the net pressure drop.

METHODS

We developed a non-invasive method to assess the pressure recovery distance based on blood flow momentum via 4D Flow cardiovascular magnetic resonance (CMR). Multi-center acquisitions included physical flow phantoms with different stenotic valve configurations to validate this method, first against reference measurements and then against turbulent energy dissipation (respectively n = 8 and n = 28 acquisitions) and to investigate the relationship between peak and net pressure drops. Finally, we explored the potential errors of cardiac catheterisation pressure recordings as a result of neglecting the pressure recovery distance in a clinical bicuspid aortic valve (BAV) cohort of n = 32 patients.

RESULTS

In-vitro assessment of pressure recovery distance based on flow momentum achieved an average error of 1.8 ± 8.4 mm when compared to reference pressure sensors in the first phantom workbench. The momentum pressure recovery distance and the turbulent energy dissipation distance showed no statistical difference (mean difference of 2.8 ± 5.4 mm, R2 = 0.93) in the second phantom workbench. A linear correlation was observed between peak and net pressure drops, however, with strong dependences on the valvular morphology. Finally, in the BAV cohort the pressure recovery distance was 78.8 ± 34.3 mm from vena contracta, which is significantly longer than currently accepted in clinical practise (40-50 mm), and 37.5% of patients displayed a pressure recovery distance beyond the end of the ascending aorta.

CONCLUSION

The non-invasive assessment of the distance to pressure recovery is possible by tracking momentum via 4D Flow CMR. Recovery is not always complete at the ascending aorta, and catheterised recordings will overestimate the net pressure drop in those situations. There is a need to re-evaluate the methods that characterise the haemodynamic burden caused by aortic stenosis as currently clinically accepted pressure recovery distance is an underestimation.

Impact of dobutamine stress on diastolic energetic efficiency of healthy left ventricle: an in vivo kinetic energy analysis

The total kinetic energy (KE) of blood can be decomposed into mean KE (MKE) and turbulent KE (TKE), which are associated with the phase-averaged fluid velocity field and the instantaneous velocity fluctuations, respectively. The aim of this study was to explore the effects of pharmacologically induced stress on MKE and TKE in the left ventricle (LV) in a cohort of healthy volunteers. 4D Flow MRI data were acquired in eleven subjects at rest and after dobutamine infusion, at a heart rate that was ∼60% higher than the one in rest conditions. MKE and TKE were computed as volume integrals over the whole LV and as data mapped to functional LV flow components, i.e., direct flow, retained inflow, delayed ejection flow and residual volume. Diastolic MKE and TKE increased under stress, in particular at peak early filling and peak atrial contraction. Augmented LV inotropy and cardiac frequency also caused an increase in direct flow and retained inflow MKE and TKE. However, the TKE/KE ratio remained comparable between rest and stress conditions, suggesting that LV intracavitary fluid dynamics can adapt to stress conditions without altering the TKE to KE balance of the normal left ventricle at rest.

Automatic Time-Resolved Cardiovascular Segmentation of 4D Flow MRI Using Deep Learning

BACKGROUND

Segmenting the whole heart over the cardiac cycle in 4D flow MRI is a challenging and time-consuming process, as there is considerable motion and limited contrast between blood and tissue.

PURPOSE

To develop and evaluate a deep learning-based segmentation method to automatically segment the cardiac chambers and great thoracic vessels from 4D flow MRI.

STUDY TYPE

Retrospective.

SUBJECTS

A total of 205 subjects, including 40 healthy volunteers and 165 patients with a variety of cardiac disorders were included. Data were randomly divided into training (n = 144), validation (n = 20), and testing (n = 41) sets.

FIELD STRENGTH/SEQUENCE

A 3 T/time-resolved velocity encoded 3D gradient echo sequence (4D flow MRI).

ASSESSMENT

A 3D neural network based on the U-net architecture was trained to segment the four cardiac chambers, aorta, and pulmonary artery. The segmentations generated were compared to manually corrected atlas-based segmentations. End-diastolic (ED) and end-systolic (ES) volumes of the four cardiac chambers were calculated for both segmentations.

STATISTICAL TESTS

Dice score, Hausdorff distance, average surface distance, sensitivity, precision, and miss rate were used to measure segmentation accuracy. Bland-Altman analysis was used to evaluate agreement between volumetric parameters.

RESULTS

The following evaluation metrics were computed: mean Dice score (0.908 ± 0.023) (mean ± SD), Hausdorff distance (1.253 ± 0.293 mm), average surface distance (0.466 ± 0.136 mm), sensitivity (0.907 ± 0.032), precision (0.913 ± 0.028), and miss rate (0.093 ± 0.032). Bland-Altman analyses showed good agreement between volumetric parameters for all chambers. Limits of agreement as percentage of mean chamber volume (LoA%), left ventricular: 9.3%, 13.5%, left atrial: 12.4%, 16.9%, right ventricular: 9.9%, 15.6%, and right atrial: 18.7%, 14.4%; for ED and ES, respectively.

DATA CONCLUSION

The addition of this technique to the 4D flow MRI assessment pipeline could expedite and improve the utility of this type of acquisition in the clinical setting.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY: Stage 1.

Non-contrast myocardial perfusion in rest and exercise stress using systolic flow-sensitive alternating inversion recovery

OBJECTIVE

To evaluate systolic flow-sensitive alternating inversion recovery (FAIR) during rest and exercise stress using 2RR (two cardiac cycles) or 1RR intervals between inversion pulse and imaging.

MATERIALS AND METHODS

1RR and 2RR FAIR was implemented on a 3T scanner. Ten healthy subjects were scanned during rest and stress. Stress was performed using an in-bore ergometer. Heart rate, mean myocardial blood flow (MBF) and temporal signal-to-noise ratio (TSNR) were compared using paired t tests.

RESULTS

Mean heart rate during stress was higher than rest for 1RR FAIR (85.8 ± 13.7 bpm vs 63.3 ± 11.1 bpm; p < 0.01) and 2RR FAIR (83.8 ± 14.2 bpm vs 63.1 ± 10.6 bpm; p < 0.01). Mean stress MBF was higher than rest for 1RR FAIR (2.97 ± 0.76 ml/g/min vs 1.43 ± 0.6 ml/g/min; p < 0.01) and 2RR FAIR (2.8 ± 0.96 ml/g/min vs 1.22 ± 0.59 ml/g/min; p < 0.01). Resting mean MBF was higher for 1RR FAIR than 2RR FAIR (p < 0.05), but not during stress. TSNR was lower for stress compared to rest for 1RR FAIR (4.52 ± 2.54 vs 10.12 ± 3.69; p < 0.01) and 2RR FAIR (7.36 ± 3.78 vs 12.41 ± 5.12; p < 0.01). 2RR FAIR TSNR was higher than 1RR FAIR for rest (p < 0.05) and stress (p < 0.001).

DISCUSSION

We have demonstrated feasibility of systolic FAIR in rest and exercise stress. 2RR delay systolic FAIR enables non-contrast perfusion assessment during stress with relatively high TSNR.

Assessment of Blood Flow Transport in the Left Ventricle Using Ultrasound. Validation Against 4-D Flow Cardiac Magnetic Resonance

Four-dimensional flow cardiac magnetic resonance (CMR) is the reference technique for analyzing blood transport in the left ventricle (LV), but similar information can be obtained from ultrasound. We aimed to validate ultrasound-derived transport in a head-to-head comparison against 4D flow CMR. In five patients and two healthy volunteers, we obtained 2D + t and 3D + t (4D) flow fields in the LV using transthoracic echocardiography and CMR, respectively. We compartmentalized intraventricular blood flow into four fractions of end-diastolic volume: direct flow (DF), retained inflow (RI), delayed ejection flow (DEF) and residual volume (RV). Using ultrasound we also computed the properties of LV filling waves (percentage of LV penetration and percentage of LV volume carried by E/A waves) to determine their relationships with CMR transport. Agreement between both techniques for quantifying transport fractions was good for DF and RV (Ric [95% confidence interval]: 0.82 [0.33, 0.97] and 0.85 [0.41, 0.97], respectively) and moderate for RI and DEF (Ric= 0.47 [-0.29, 0.88] and 0.55 [-0.20, 0.90], respectively). Agreement between techniques to measure kinetic energy was variable. The amount of blood carried by the E-wave correlated with DF and RV (R = 0.75 and R = 0.63, respectively). Therefore, ultrasound is a suitable method for expanding the analysis of intraventricular flow transport in the clinical setting.

Ectopic fat is associated with cardiac remodeling—A comprehensive assessment of regional fat depots in type 2 diabetes using multi-parametric MRI

Background

Different regional depots of fat have distinct metabolic properties and may relate differently to adverse cardiac remodeling. We sought to quantify regional depots of body fat and to investigate their relationship to cardiac structure and function in Type 2 Diabetes (T2D) and controls.

Methods

From the SCAPIS cohort in Linköping, Sweden, we recruited 92 subjects (35% female, mean age 59.5 ± 4.6 years): 46 with T2D and 46 matched controls. In addition to the core SCAPIS data collection, participants underwent a comprehensive magnetic resonance imaging examination at 1.5 T for assessment of left ventricular (LV) structure and function (end-diastolic volume, mass, concentricity, ejection fraction), as well as regional body composition (liver proton density fat fraction, visceral adipose tissue, abdominal subcutaneous adipose tissue, thigh muscle fat infiltration, fat tissue-free thigh muscle volume and epicardial adipose tissue).

Results

Compared to the control group, the T2D group had increased: visceral adipose tissue volume index (P &lt; 0.001), liver fat percentage (P &lt; 0.001), thigh muscle fat infiltration percentage (P = 0.02), LV concentricity (P &lt; 0.001) and LV E/e'-ratio (P &lt; 0.001). In a multiple linear regression analysis, a negative association between liver fat percentage and LV mass (St Beta −0.23, P &lt; 0.05) as well as LV end-diastolic volume (St Beta −0.27, P &lt; 0.05) was found. Epicardial adipose tissue volume and abdominal subcutaneous adipose tissue volume index were the only parameters of fat associated with LV diastolic dysfunction (E/e'-ratio) (St Beta 0.24, P &lt; 0.05; St Beta 0.34, P &lt; 0.01, respectively). In a multivariate logistic regression analysis, only visceral adipose tissue volume index was significantly associated with T2D, with an odds ratio for T2D of 3.01 (95% CI 1.28–7.05, P &lt; 0.05) per L/m2 increase in visceral adipose tissue volume.

Conclusions

Ectopic fat is predominantly associated with cardiac remodeling, independently of type 2 diabetes. Intriguingly, liver fat appears to be related to LV structure independently of VAT, while epicardial fat is linked to impaired LV diastolic function. Visceral fat is associated with T2D independently of liver fat and abdominal subcutaneous adipose tissue.

A head-to-head comparison of myocardial strain by fast-strain encoding and feature tracking imaging in acute myocardial infarction

Background

Myocardial infarction (MI) is a major cause of heart failure. Left ventricular adverse remodeling is common post-MI. Several studies have demonstrated a correlation between reduced myocardial strain and the development of adverse remodeling. Cardiac magnetic resonance (CMR) with fast-strain encoding (fast-SENC) or feature tracking (FT) enables rapid assessment of myocardial deformation. The aim of this study was to establish a head-to-head comparison of fast-SENC and FT in post-ST-elevated myocardial infarction (STEMI) patients, with clinical 2D speckle tracking echocardiography (2DEcho) as a reference.

Methods

Thirty patients treated with primary percutaneous coronary intervention for STEMI were investigated. All participants underwent CMR examination with late gadolinium enhancement, cine-loop steady-state free precession, and fast-SENC imaging using a 1.5T scanner as well as a 2DEcho. Global longitudinal strain (GLS), segmental longitudinal strain (SLS), global circumferential strain (GCS), and segmental circumferential strain (SCS) were assessed along with the MI scar extent.

Results

The GCS measurements from fast-SENC and FT were nearly identical: the mean difference was 0.01 (2.5)% (95% CI - 0.92 to 0.95). For GLS, fast-SENC values were higher than FT, with a mean difference of 1.8 (1.4)% (95% CI 1.31-2.35). Tests of significance for GLS did not show any differences between the MR methods and 2DEcho. Average strain in the infarct-related artery (IRA) segments compared to the remote myocardium was significantly lower for the left anterior descending artery and right coronary artery culprits but not for the left circumflex artery culprits. Fast-SENC displayed a higher area under the curve for detecting infarcted segments than FT for both SCS and SLS.

Conclusion

GLS and GCS did not significantly differ between fast-SENC and FT. Both showed acceptable agreement with 2DEcho for longitudinal strain. Segments perfused by the IRA showed significantly reduced strain values compared to the remote myocardium. Fast-SENC presented a higher sensitivity and specificity for detecting infarcted segments than FT.

Four-dimensional flow cardiac magnetic resonance assessment of left ventricular diastolic function

Left ventricular diastolic dysfunction is a major cause of heart failure and carries a poor prognosis. Assessment of left ventricular diastolic function however remains challenging for both echocardiography and conventional phase contrast cardiac magnetic resonance. Amongst other limitations, both are restricted to measuring velocity in a single direction or plane, thereby compromising their ability to capture complex diastolic hemodynamics in health and disease. Time-resolved three-dimensional phase contrast cardiac magnetic resonance imaging with three-directional velocity encoding known as '4D flow CMR' is an emerging technology which allows retrospective measurement of velocity and by extension flow at any point in the acquired 3D data volume. With 4D flow CMR, complex aspects of blood flow and ventricular function can be studied throughout the cardiac cycle. 4D flow CMR can facilitate the visualization of functional blood flow components and flow vortices as well as the quantification of novel hemodynamic and functional parameters such as kinetic energy, relative pressure, energy loss and vorticity. In this review, we examine key concepts and novel markers of diastolic function obtained by flow pattern analysis using 4D flow CMR. We consolidate the existing evidence base to highlight the strengths and limitations of 4D flow CMR techniques in the surveillance and diagnosis of left ventricular diastolic dysfunction.

Turbulent Intensity of Blood Flow in the Healthy Aorta Increases With Dobutamine Stress and is Related to Cardiac Output

Introduction: The blood flow in the normal cardiovascular system is predominately laminar but operates close to the threshold to turbulence. Morphological distortions such as vascular and valvular stenosis can cause transition into turbulent blood flow, which in turn may cause damage to tissues in the cardiovascular system. A growing number of studies have used magnetic resonance imaging (MRI) to estimate the extent and degree of turbulent flow in different cardiovascular diseases. However, the way in which heart rate and inotropy affect turbulent flow has not been investigated. In this study we hypothesized that dobutamine stress would result in higher turbulence intensity in the healthy thoracic aorta.

Method: 4D flow MRI data were acquired in twelve healthy subjects at rest and with dobutamine, which was infused until the heart rate increased by 60% when compared to rest. A semi-automatic segmentation method was used to segment the thoracic aorta in the 4D flow MR images. Subsequently, flow velocity and several turbulent kinetic energy (TKE) parameters were calculated in the ascending aorta, aortic arch, descending aorta and whole thoracic aorta.

Results: With dobutamine infusion there was an increase in heart rate (66 ± 9 vs. 108 ± 13 bpm, p < 0.001) and stroke volume (88 ± 13 vs. 102 ± 25 ml, p < 0.01). Additionally, there was an increase in Peak Average velocity (0.7 ± 0.1 vs. 1.2 ± 0.2 m/s, p < 0.001, Peak Max velocity (1.3 ± 0.1 vs. 2.0 ± 0.2 m/s, p < 0.001), Peak Total TKE (2.9 ± 0.7 vs. 8.0 ± 2.2 mJ, p < 0.001), Peak Median TKE (36 ± 7 vs. 93 ± 24 J/m3, p = 0.002) and Peak Max TKE (176 ± 33 vs. 334 ± 69 J/m3, p < 0.001). The relation between cardiac output and Peak Total TKE in the whole thoracic aorta was very strong (R² = 0.90, p < 0.001).

Conclusion: TKE of blood flow in the healthy thoracic aorta increases with dobutamine stress and is strongly related to cardiac output. Quantification of such turbulence intensity parameters with cardiac stress may serve as a risk assessment of aortic disease development.

Wall shear stress and relative residence time as potential risk factors for abdominal aortic aneurysms in males: a 4D flow cardiovascular magnetic resonance case-control study

BACKGROUND

Abdominal aortic aneurysms (AAA) can lead to catastrophic events such as dissection or rupture, and are an expression of general aortic disease. Low wall shear stress (WSS), high oscillatory shear index (OSI), and high relative residence time (RRT) have been correlated against increased uptake of inflammatory markers in the vessel wall and may improve risk stratification of AAA. We sought to obtain a comprehensive view of WSS, OSI, and RRT in the whole aorta for patients with AAA and age-matched elderly controls and young normal controls.

METHODS

4D Flow cardiovascular magnetic resonance images of the whole aorta were acquired in 18 AAA patients (70.8 ± 3.4 years), 22 age-matched controls (71.4 ± 3.4 years), and 23 young subjects (23.3 ± 3.1 years), all males. Three-dimensional segmentations of the whole aorta were created for all timeframes using a semi-automatic approach. The aorta was divided into five segments: ascending aorta, arch, descending aorta, suprarenal and infrarenal abdominal aorta. For each segment, average values of peak WSS, OSI, and RRT were computed. Student's t-tests were used to compare values between the three cohorts (AAA patients vs elderly controls, and elderly controls vs young controls) where the data were normally distributed, and the non-parametric Wilcoxon rank sum tests were used otherwise.

RESULTS

AAA patients had lower peak WSS in the descending aorta as well as in the abdominal aorta compared to elderly controls (p ≤ 0.001), similar OSI, but higher RRT in the descending and abdominal aorta (p ≤ 0.001). Elderly controls had lower peak WSS compared to young controls throughout the aorta (p < 0.001), higher OSI in all segments except for the infrarenal aorta (p < 0.001), and higher RRT throughout the aorta, except the infrarenal aorta (p < 0.001).

CONCLUSIONS

This study provides novel insights into WSS, OSI, and RRT in patients with AAA in relation to normal ageing, highlighting how AAA patients have markedly abnormal hemodynamic stresses not only in the infrarenal, but in the entire aorta. Moreover, we identified RRT as a marker for abnormal AAA hemodynamics. Further investigations are needed to explore if RRT or other measures of hemodynamics stresses best predict AAA growth and/or rupture.

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.

Cardiac Kinetic Energy and Viscous Dissipation Rate From Radial Flow Data

Recent studies have correlated kinetic energy (KE) and viscous dissipation rate (VDR) in the left ventricle (LV) with heart health. These studies have relied on 4D-flow imaging or computational fluid dynamics modeling, which are able to measure, or compute, all 3 components (3C) of the blood flow velocity in 3 dimensional (3D) space. This richness of data is difficult to acquire clinically. Alternatively, color Doppler echocardiography (CDE) is more widespread clinically, but only measures a single radial component of velocity and typically only over a planar section. Because of this limitation, prior CDE-based studies have first reconstructed a second component of velocity in the measurement plane prior to evaluating VDR or KE. Herein, we propose 1C-based surrogates of KE and VDR that can be derived directly from the radial component of the flow velocity in the LV. Our results demonstrate that the proposed 1C-based surrogates of KE and VDR are generally as well-correlated with the true KE and VDR values as surrogates that use reconstructed 2C flow data. Moreover, the correlation of these 1C-based surrogates with the true values indicate that CDE (3D in particular) may be useful in evaluating these metrics in practice.

The transluminal attenuation gradient does not add diagnostic accuracy to coronary computed tomography

BACKGROUND

A method for improving the accuracy of coronary computed tomography angiography (CCTA) is highly sought after as it would help to avoid unnecessary invasive coronary angiographies. Measurement of the transluminal attenuation gradient (TAG) has been proposed as an alternative to other existing methods, i.e. CT perfusion and CT fractional flow reserve (FFR).

PURPOSE

To evaluate the incremental value of three types of TAG in high-pitch spiral CCTA with invasive FFR measurements as reference.

MATERIAL AND METHODS

TAG was measured using two semi-automatic methods and one manual method. A receiver operating characteristic (ROC) analysis was made to determine the usefulness of TAG alone as well as TAG combined with CCTA for detection of significant coronary artery stenoses defined by an invasive FFR value ≤0.80.

RESULTS

A total of 51 coronary vessels in 37 patients were included in this retrospective study. Hemodynamically significant stenoses were found in 13 vessels according to FFR. The ROC analysis TAG alone resulted in areas under the curve (AUCs) of 0.530 and 0.520 for the semi-automatic TAG and 0.557 for the manual TAG. TAG and CCTA combined resulted in AUCs of 0.567, 0.562 for semi-automatic TAG, and 0.569 for the manual TAG.

CONCLUSION

The results from our study showed no incremental value of TAG measured in single heartbeat CCTA in determining the severity of coronary artery stenosis degrees.

Evaluating the prevalence and severity of NAFLD in primary care: the EPSONIP study protocol

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) affects 20-30% of the general adult population. NAFLD patients with type 2 diabetes mellitus (T2DM) are at an increased risk of advanced fibrosis, which puts them at risk of cardiovascular complications, hepatocellular carcinoma, or liver failure. Liver biopsy is the gold standard for assessing hepatic fibrosis. However, its utility is inherently limited. Consequently, the prevalence and characteristics of T2DM patients with advanced fibrosis are unknown. Therefore, the purpose of the current study is to evaluate the prevalence and severity of NAFLD in patients with T2DM by recruiting participants from primary care, using the latest imaging modalities, to collect a cohort of well phenotyped patients.

METHODS

We will prospectively recruit 400 patients with T2DM using biomarkers to assess their status. Specifically, we will evaluate liver fat content using magnetic resonance imaging (MRI); hepatic fibrosis using MR elastography and vibration-controlled transient elastography; muscle composition and body fat distribution using water-fat separated whole body MRI; and cardiac function, structure, and tissue characteristics, using cardiovascular MRI.

DISCUSSION

We expect that the study will uncover potential mechanisms of advanced hepatic fibrosis in NAFLD and T2DM and equip the clinician with better diagnostic tools for the care of T2DM patients with NAFLD.

TRIAL REGISTRATION

Clinicaltrials.gov, identifier NCT03864510. Registered 6 March 2019, https://clinicaltrials.gov/ct2/show/NCT03864510 .

Using Deep Learning to Emulate the Use of an External Contrast Agent in Cardiovascular 4D Flow MRI

BACKGROUND

Although contrast agents would be beneficial, they are seldom used in four-dimensional (4D) flow magnetic resonance imaging (MRI) due to potential side effects and contraindications.

PURPOSE

To develop and evaluate a deep learning architecture to generate high blood-tissue contrast in noncontrast 4D flow MRI by emulating the use of an external contrast agent.

STUDY TYPE

Retrospective.

SUBJECTS

Of 222 data sets, 141 were used for neural network (NN) training (69 with and 72 without contrast agent). Evaluation was performed on the remaining 81 noncontrast data sets.

FIELD STRENGTH/SEQUENCES

Gradient echo or echo-planar 4D flow MRI at 1.5 T and 3 T.

ASSESSMENT

A cyclic generative adversarial NN was trained to perform image translation between noncontrast and contrast data. Evaluation was performed quantitatively using contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), structural similarity index (SSIM), mean squared error (MSE) of edges, and Dice coefficient of segmentations. Three observers performed a qualitative assessment of blood-tissue contrast, noise, presence of artifacts, and image structure visualization.

STATISTICAL TESTS

The Wilcoxon rank-sum test evaluated statistical significance. Kendall's concordance coefficient assessed interobserver agreement.

RESULTS

Contrast in the regions of interest (ROIs) in the NN enhanced images increased by 88%, CNR increased by 63%, and SNR improved by 48% (all P < 0.001). The SSIM was 0.82 ± 0.01, and the MSE of edges was 0.09 ± 0.01 (range [0,1]). Segmentations based on the generated images resulted in a Dice similarity increase of 15.25%. The observers managed to differentiate between contrast MR images and our results; however, they preferred the NN enhanced images in 76.7% of cases. This percentage increased to 93.3% for phase-contrast MR angiograms created from the NN enhanced data. Visual grading scores were blood-tissue contrast = 4.30 ± 0.74, noise = 3.12 ± 0.98, and presence of artifacts = 3.63 ± 0.76. Image structures within and without the ROIs resulted in scores of 3.42 ± 0.59 and 3.07 ± 0.71, respectively (P < 0.001).

DATA CONCLUSION

The proposed approach improves blood-tissue contrast in MR images and could be used to improve data quality, visualization, and postprocessing of cardiovascular 4D flow data. EVIDENCE LEVEL: 3 TECHNICAL EFFICACY: Stage 1.

Impact of prosthetic mitral valve orientation on the ventricular flow field: Comparison using patient-specific computational fluid dynamics

Significant mitral valve regurgitation creates progressive adverse remodeling of the left ventricle (LV). Replacement of the failing valve with a prosthesis generally improves patient outcomes but leaves the patient with non-physiological intracardiac flow patterns that might contribute to their future risk of thrombus formation and embolism. It has been suggested that the angular orientation of the implanted valve might modify the postoperative distortion of the intraventricular flow field. In this study, we investigated the effect of prosthetic valve orientation on LV flow patterns by using heart geometry from a patient with LV dysfunction and a competent native mitral valve to calculate intracardiac flow fields with computational fluid dynamics (CFD). Results were validated using in vivo 4D Flow MRI. The computed flow fields were compared to calculations following virtual implantation of a mechanical heart valve oriented in four different angles to assess the effect of leaflet position. Flow patterns were visualized in long- and short-axes and quantified with flow component analysis. In comparison to a native valve, valve implantation increased the proportion of the mitral inflow remaining in the basal region and further increased the residual volume in the apical area. Only slight changes due to valve orientation were observed. Using our numerical framework, we demonstrated quantitative changes in left ventricular blood flow due to prosthetic mitral replacement. This framework may be used to improve design of prosthetic heart valves and implantation procedures to minimize the potential for apical flow stasis, and potentially assist personalized treatment planning.

Clinical validation of three cardiovascular magnetic resonance techniques to measure strain and torsion in patients with suspected coronary artery disease

BACKGROUND

Several cardiovascular magnetic resonance (CMR) techniques can measure myocardial strain and torsion with high accuracy. The purpose of this study was to compare displacement encoding with stimulated echoes (DENSE), tagging and feature tracking (FT) for measuring circumferential and radial myocardial strain and myocardial torsion in order to assess myocardial function and infarct scar burden both at a global and at a segmental level.

METHOD

116 patients with a high likelihood of coronary artery disease (European SCORE > 15%) underwent CMR examination including cine images, tagging, DENSE and late gadolinium enhancement (LGE) in the short axis direction. In total, 97 patients had signs of myocardial disease and 19 had no abnormalities in terms of left ventricular (LV) wall mass index, LV ejection fraction, wall motion, LGE or a history of myocardial infarction. Thirty-four patients had myocardial infarct scar with a transmural LGE extent (transmurality) that exceeded 50% of the wall thickness in at least one segment. Global circumferential strain (GCS) and global radial strain (GRS) was analyzed using FT of cine loops, deformation of tag lines or DENSE displacement.

RESULTS

DENSE and tagging both showed high sensitivity (82% and 71%) at a specificity of 80% for the detection of segments with > 50% LGE transmurality, and receiver operating characteristics (ROC) analysis showed significantly higher area under the curve-values (AUC) for DENSE (0.87) than for tagging (0.83, p < 0.001) and FT (0.66, p = 0.003). GCS correlated with global LGE when determined with DENSE (r = 0.41), tagging (r = 0.37) and FT (r = 0.15). GRS had a low but significant negative correlation with LGE; DENSE r = - 0.10, FT r = - 0.07 and tagging r = - 0.16. Torsion from DENSE and tagging had a weak correlation (- 0.20 and - 0.22 respectively) with global LGE.

CONCLUSION

Circumferential strain from DENSE detected segments with > 50% scar with a higher AUC than strain determined from tagging and FT at a segmental level. GCS and torsion computed from DENSE and tagging showed similar correlation with global scar size, while when computed from FT, the correlation was lower.

Improved Efficiency of Intraventricular Blood Flow Transit Under Cardiac Stress: A 4D Flow Dobutamine CMR Study

Introduction: The effects of heart rate, inotropy, and lusitropy on multidimensional flow patterns and energetics within the human heart remain undefined. Recently, reduced volume and end-diastolic kinetic energy (KE) of the portion of left ventricular (LV) inflow passing directly to outflow, Direct flow (DF), have been shown to reflect inefficient LV pumping and to be a marker of LV dysfunction in heart failure patients. In this study, we hypothesized that increasing heart rate, inotropy, and lusitropy would result in an increased efficiency of intraventricular blood flow transit. Therefore, we sought to investigate LV 4D blood flow patterns and energetics with dobutamine infusion.

Methods: 4D flow and morphological cardiovascular magnetic resonance (CMR) data were acquired in twelve healthy subjects: at rest and with dobutamine infusion to achieve a target heart rate ~60% higher than the resting heart rate. A previously validated method was used for flow analysis: pathlines were emitted from the end-diastolic (ED) LV blood volume and traced forward and backward in time to separate four functional LV flow components. For each flow component, KE/mL blood volume at ED was calculated.

Results: With dobutamine infusion there was an increase in heart rate (64%, p < 0.001), systolic blood pressure (p = 0.02) and stroke volume (p = 0.01). Of the 4D flow parameters, the most efficient flow component (DF), increased its proportion of EDV (p < 0.001). The EDV proportion of Residual volume, the blood residing in the ventricle over at least two cardiac cycles, decreased (p < 0.001). The KE/mL at ED for all flow components increased (p < 0.001). DF had the largest absolute and relative increase while Residual volume had the smallest absolute and relative increase.

Conclusions: This study demonstrates that it is feasible to compare 4D flow patterns within the normal human heart at rest and with stress. At higher heart rate, inotropy and lusitropy, elicited by dobutamine infusion, the efficiency of intraventricular blood flow transit improves, as quantified by an increased relative volume and pre-systolic KE of the most efficient DF component of the LV volume. The change in these markers may allow a novel assessment of LV function and LV dysfunction over a range of stress.

In vitro experiments on ICOSA6 4D flow MRI measurement for the quantification of velocity and turbulence parameters

PURPOSE

To perform comprehensive in vitro experiments using six-directional icosahedral flow encoding (ICOSA6) 4D flow magnetic resonance imaging (MRI) under various scan conditions to analyze the robustness of velocity and turbulence quantification.

MATERIALS AND METHODS

In vitro flow phantoms with steady flow rates of 10 and 20 L/min were scanned using both conventional 4D flow MRI and ICOSA6. Experiments focused on comparisons between ICOSA6 and conventional four point (4P) methods, and the effects of contrast agents, velocity encoding range (Venc), and scan direction on velocity and turbulence quantification.

RESULTS

The results demonstrated that 1) ICOSA6 improves the velocity-to-noise ratio (VNR) of velocity estimation by 33% (on average) and results in similar turbulent kinetic energy (TKE) estimation as the 4P method. 2) Measurements with a contrast agent resulted in more than a 2.5 fold increase in average VNR. However, the improvement of total TKE quantification was not obvious. 3) TKE estimation was less affected by Venc and the scan direction, whereas turbulence production (TP) estimation was largely affected by these measurement conditions. The effects of Venc and scan direction accounted for less than 11.63% of TKE estimation, but up to 33.89% of TP estimation.

CONCLUSION

The ICOSA6 scheme is compatible with conventional 4D flow MRI for velocity and TKE measurement. Contrast agents are effective at increasing VNR, but not signal-to-noise ratio for TKE quantification. The effects of Venc and scan direction influence total TP more than total TKE.

Left Ventricular Flow Analysis

BACKGROUND

Cardiac remodeling, after a myocardial insult, often causes progression to heart failure. The relationship between alterations in left ventricular blood flow, including kinetic energy (KE), and remodeling is uncertain. We hypothesized that increasing derangements in left ventricular blood flow would relate to (1) conventional cardiac remodeling markers, (2) increased levels of biochemical remodeling markers, (3) altered cardiac energetics, and (4) worsening patient symptoms and functional capacity. Methods Thirty-four dilated cardiomyopathy patients, 30 ischemic cardiomyopathy patients, and 36 controls underwent magnetic resonance including 4-dimensional flow, BNP (brain-type natriuretic peptide) measurement, functional capacity assessment (6-minute walk test), and symptom quantification. A subgroup of dilated cardiomyopathy and control subjects underwent cardiac energetic assessment. Left ventricular flow was separated into 4 components: direct flow, retained inflow, delayed ejection flow, and residual volume. Average KE throughout the cardiac cycle was calculated.

RESULTS

Patients had reduced direct flow proportion and direct-flow average KE compared with controls ( P<0.0001). The residual volume proportion and residual volume average KE were increased in patients ( P<0.0001). Importantly, in a multiple linear regression model to predict the patient's 6-minute walk test, the independent predictors were age (β=-0.3015; P=0.019) and direct-flow average KE (β=0.280, P=0.035; R² model, 0.466, P=0.002). In contrast, neither ejection fraction nor left ventricular volumes were independently predictive.

CONCLUSIONS

This study demonstrates an independent predictive relationship between the direct-flow average KE and a prognostic measure of functional capacity. Intracardiac 4-dimensional flow parameters are novel biomarkers in heart failure and may provide additive value in monitoring new therapies and predicting prognosis.

Non-invasive estimation of relative pressure in turbulent flow using virtual work-energy

Vascular pressure differences are established risk markers for a number of cardiovascular diseases. Relative pressures are, however, often driven by turbulence-induced flow fluctuations, where conventional non-invasive methods may yield inaccurate results. Recently, we proposed a novel method for non-turbulent flows, νWERP, utilizing the concept of virtual work-energy to accurately probe relative pressure through complex branching vasculature. Here, we present an extension of this approach for turbulent flows: νWERP-t. We present a theoretical method derivation based on flow covariance, quantifying the impact of flow fluctuations on relative pressure. νWERP-t is tested on a set of in-vitro stenotic flow phantoms with data acquired by 4D flow MRI with six-directional flow encoding, as well as on a patient-specific in-silico model of an acute aortic dissection. Over all tests νWERP-t shows improved accuracy over alternative energy-based approaches, with excellent recovery of estimated relative pressures. In particular, the use of a guaranteed divergence-free virtual field improves accuracy in cases where turbulent flows skew the apparent divergence of the acquired field. With the original νWERP allowing for assessment of relative pressure into previously inaccessible vasculatures, the extended νWERP-t further enlarges the method's clinical scope, underlining its potential as a novel tool for assessing relative pressure in-vivo.

Data Quality and Optimal Background Correction Order of Respiratory-Gated k-Space Segmented Spoiled Gradient Echo (SGRE) and Echo Planar Imaging (EPI)-Based 4D Flow MRI

Background

A reduction in scan time of 4D Flow MRI would facilitate clinical application. A recent study indicates that echo‐planar imaging (EPI) 4D Flow MRI allows for a reduction in scan time and better data quality than the recommended k‐space segmented spoiled gradient echo (SGRE) sequence. It was argued that the poor data quality of SGRE was related to the nonrecommended absence of respiratory motion compensation. However, data quality can also be affected by the background offset compensation.

Purpose

To compare the data quality of respiratory motion‐compensated SGRE and EPI 4D Flow MRI and their dependence on background correction (BC) order.

Study Type

Retrospective.

Subjects

Eighteen healthy subjects (eight female, mean age 32 ± 5 years).

Field Strength and Sequence

1.5 T. [Correction added on July 26, 2019, after first online publication: The preceding field strength was corrected.] SGRE and EPI‐based 4D Flow MRI.

Assessment

Data quality was investigated visually and by comparing flows through the cardiac valves and aorta. Measurements were obtained from transvalvular flow and pathline analysis.

Statistical Tests

Linear regression and Bland–Altman analysis were used. Wilcoxon test was used for comparison of visual scoring. Student's t‐test was used for comparison of flow volumes.

Results

No significant difference was found by visual inspection (P = 0.08). Left ventricular (LV) flows were strongly and very strongly associated with SGRE and EPI, respectively (R² = 0.86–0.94 SGRE; 0.71–0.79 EPI, BC0–4). LV and right ventricular (RV) outflows and LV pathline flows were very strongly associated (R² = 0.93–0.95 SGRE; 0.88–0.91 EPI, R² = 0.91–0.95 SGRE; 0.91–0.93 EPI, BC1–4). EPI LV outflow was lower than the short‐axis‐based stroke volume. EPI RV outflow and proximal descending aortic flow were lower than SGREs.

Data Conclusion

Both sequences yielded good internal data consistency when an adequate background correction was applied. Second and first BC order were considered sufficient for transvalvular flow analysis in SGRE and EPI, respectively. Higher BC orders were preferred for particle tracing.

Level of Evidence 4

Technical Efficacy Stage 1

J. Magn. Reson. Imaging 2020;51:885–896.

Intracardiac 4D Flow MRI in Congenital Heart Disease: Recommendations on Behalf of the ISMRM Flow & Motion Study Group

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2019;50:677-681.

Post-cardioversion Improvement in LV Function Defined by 4D Flow Patterns and Energetics in Patients With Atrial Fibrillation

Background

Atrial fibrillation (AF) is a prevalent cause of cardiovascular morbidity, including thromboembolism and heart failure. Left ventricular dysfunction (LVD) detected in AF patients may be either precursor or consequence of the arrythmia. Successful cardioversion of chronic AF is often followed by a transient period of left atrial (LA) stunning, where depressed mechanical atrial contraction persists despite reinstitution of sinus rhythm. To determine if AF-associated LVD would improve with resolution of LA dysfunction, AF patients were examined immediately and 4 weeks after cardioversion to sinus rhythm. 4D flow cardiovascular magnetic resonance (CMR) assesses ventricular function according to the volumes and energetics of functional components of the LV volume. Previously, described 4D CMR markers of LVD include decreased volume and end-diastolic kinetic energy (KE) of the Direct flow, which is the portion of LV volume that passes directly from inflow to outflow in a single cycle. We hypothesize that impaired LV flow patterns and energetics will be found immediately after cardioversion during atrial stunning, and that those parameters will improve as atrial function returns.

Methods

Ten patients with a history of AF underwent CMR 2-3 h (Time-1) and 4 weeks (Time-2), following electrical cardioversion to sinus rhythm. 4D phase-contrast velocity data and morphological images were acquired at a 3T CMR system. Using a previously evaluated method, pathlines were emitted from the LV end diastolic volume (LVEDV) and traced forward and backward in time until end-systole. The LVEDV was automatically separated into four functional flow components whose volume and KE were calculated.

Results

Left atrial fractional area change increased over the follow-up period (P = 0.001), indicating recovery of LA mechanical function. LVEF increased between Time-1 and Time-2 (P = 0.003); LVEDVI did not change (P = 0.319). Over that interval, the ratios of Direct flow/LVEDV volume and KE increased (P = 0.001 and P = 0.003, respectively), while the ratios of Residual volume/LVEDV volume and KE decreased (P = 0.001 and P = 0.005, respectively).

Conclusion

Post-cardioversion recovery of LA function was associated with improvements in conventional and 4D CMR markers of LV function. Flow-specific measures demonstrate the negative but potentially reversible impact of LA dysfunction on volume and energetic aspects of LV function.

Pathophysiology behind prolonged whiplash associated disorders: study protocol for an experimental study

BACKGROUND

There is insufficient knowledge of pathophysiological parameters to understand the mechanism behind prolonged whiplash associated disorders (WAD), and it is not known whether or not changes can be restored by rehabilitation. The aims of the projects are to investigate imaging and molecular biomarkers, cervical kinaesthesia, postural sway and the association with pain, disability and other outcomes in individuals with longstanding WAD, before and after a neck-specific exercise intervention. Another aim is to compare individuals with WAD with healthy controls.

METHODS

Participants are a sub-group (n = 30) of individuals recruited from an ongoing randomized controlled study (RCT). Measurements in this experimental prospective study will be carried out at baseline (before intervention) and at a three month follow-up (end of physiotherapy intervention), and will include muscle structure and inflammation using magnetic resonance imaging (MRI), brain structure and function related to pain using functional MRI (fMRI), muscle function using ultrasonography, biomarkers using samples of blood and saliva, cervical kinaesthesia using the "butterfly test" and static balance test using an iPhone app. Association with other measures (self-reported and clinical measures) obtained in the RCT (e.g. background data, pain, disability, satisfaction with care, work ability, quality of life) may be investigated. Healthy volunteers matched for age and gender will be recruited as controls (n = 30).

DISCUSSION

The study results may contribute to the development of improved diagnostics and improved rehabilitation methods for WAD.

TRIAL REGISTRATION

Clinicaltrial.gov Protocol ID: NCT03664934, initial release 09/11/2018.

Non-invasive Assessment of Systolic and Diastolic Cardiac Function During Rest and Stress Conditions Using an Integrated Image-Modeling Approach

Background: The possibility of non-invasively assessing load-independent parameters characterizing cardiac function is of high clinical value. Typically, these parameters are assessed during resting conditions. However, for diagnostic purposes, the parameter behavior across a physiologically relevant range of heart rate and loads is more relevant than the isolated measurements performed at rest. This study sought to evaluate changes in non-invasive estimations of load-independent parameters of left-ventricular contraction and relaxation patterns at rest and during dobutamine stress. Methods: We applied a previously developed approach that combines non-invasive measurements with a physiologically-based, reduced-order model of the cardiovascular system to provide subject-specific estimates of parameters characterizing left ventricular function. In this model, the contractile state of the heart at each time point along the cardiac cycle is modeled using a time-varying elastance curve. Non-invasive data, including four-dimensional magnetic resonance imaging (4D Flow MRI) measurements, were acquired in nine subjects without a known heart disease at rest and during dobutamine stress. For each of the study subjects, we constructed two personalized models corresponding to the resting and the stress state. Results: Applying the modeling framework, we identified significant increases in the left ventricular contraction rate constant [from 1.5 ± 0.3 to 2 ± 0.5 (p = 0.038)] and relaxation constant [from 37.2 ± 6.9 to 46.1 ± 12 (p = 0.028)]. In addition, we found a significant decrease in the elastance diastolic time constant from 0.4 ± 0.04 s to 0.3 ± 0.03 s (p = 0.008). Conclusions: The integrated image-modeling approach allows the assessment of cardiovascular function given as model-based parameters. The agreement between the estimated parameter values and previously reported effects of dobutamine demonstrates the potential of the approach to assess advanced metrics of pathophysiology that are otherwise difficult to obtain non-invasively in clinical practice.

Impact of Pulmonary Venous Inflow on Cardiac Flow Simulations: Comparison with In Vivo 4D Flow MRI

Blood flow simulations are making their way into the clinic, and much attention is given to estimation of fractional flow reserve in coronary arteries. Intracardiac blood flow simulations also show promising results, and here the flow field is expected to depend on the pulmonary venous (PV) flow rates. In the absence of in vivo measurements, the distribution of the flow from the individual PVs is often unknown and typically assumed. Here, we performed intracardiac blood flow simulations based on time-resolved computed tomography on three patients, and investigated the effect of the distribution of PV flow rate on the flow field in the left atrium and ventricle. A design-of-experiment approach was used, where PV flow rates were varied in a systematic manner. In total 20 different simulations were performed per patient, and compared to in vivo 4D flow MRI measurements. Results were quantified by kinetic energy, mitral valve velocity profiles and root-mean-square errors of velocity. While large differences in atrial flow were found for varying PV inflow distributions, the effect on ventricular flow was negligible, due to a regularizing effect by mitral valve. Equal flow rate through all PVs most closely resembled in vivo measurements and is recommended in the absence of a priori knowledge.

Validation of pressure drop assessment using 4D flow MRI-based turbulence production in various shapes of aortic stenoses

PURPOSE

To validate pressure drop measurements using 4D flow MRI-based turbulence production in various shapes of stenotic stenoses.

METHODS

In vitro flow phantoms with seven different 3D-printed aortic valve geometries were constructed and scanned with 4D flow MRI with six-directional flow encoding (ICOSA6). The pressure drop through the valve was non-invasively predicted based on the simplified Bernoulli, the extended Bernoulli, the turbulence production, and the shear-scaling methods. Linear regression and agreement of the predictions with invasively measured pressure drop were analyzed.

RESULTS

All pressure drop predictions using 4D Flow MRI were linearly correlated to the true pressure drop but resulted in different regression slopes. The regression slope and 95% limits of agreement for the simplified Bernoulli method were 1.35 and 11.99 ± 21.72 mm Hg. The regression slope and 95% limits of agreement for the extended Bernoulli method were 1.02 and 0.74 ± 8.48 mm Hg. The regression slope and 95% limits of agreement for the turbulence production method were 0.89 and 0.96 ± 8.01 mm Hg. The shear-scaling method presented good correlation with an invasively measured pressure drop, but the regression slope varied between 0.36 and 1.00 depending on the shear-scaling coefficient.

CONCLUSION

The pressure drop assessment based on the turbulence production method agrees well with the extended Bernoulli method and invasively measured pressure drop in various shapes of the aortic valve. Turbulence-based pressure drop estimation can, as a complement to the conventional Bernoulli method, play a role in the assessment of valve diseases.

Test-retest variability of left ventricular 4D flow cardiovascular magnetic resonance measurements in healthy subjects

BACKGROUND

Quantification and visualisation of left ventricular (LV) blood flow is afforded by three-dimensional, time resolved phase contrast cardiovascular magnetic resonance (CMR 4D flow). However, few data exist upon the repeatability and variability of these parameters in a healthy population. We aimed to assess the repeatability and variability over time of LV 4D CMR flow measurements.

METHODS

Forty five controls underwent CMR 4D flow data acquisition. Of these, 10 underwent a second scan within the same visit (scan-rescan), 25 returned for a second visit (interval scan; median interval 52 days, IQR 28-57 days). The LV-end diastolic volume (EDV) was divided into four flow components: 1) Direct flow: inflow that passes directly to ejection; 2) Retained inflow: inflow that enters and resides within the LV; 3) Delayed ejection flow: starts within the LV and is ejected and 4) Residual volume: blood that resides within the LV for > 2 cardiac cycles. Each flow components' volume was related to the EDV (volume-ratio). The kinetic energy at end-diastole (ED) was measured and divided by the components' volume.

RESULTS

The dominant flow component in all 45 controls was the direct flow (volume ratio 38 ± 4%) followed by the residual volume (30 ± 4%), then delayed ejection flow (16 ± 3%) and retained inflow (16 ± 4%). The kinetic energy at ED for each component was direct flow (7.8 ± 3.0 microJ/ml), retained inflow (4.1 ± 2.0 microJ/ml), delayed ejection flow (6.3 ± 2.3 microJ/ml) and the residual volume (1.2 ± 0.5 microJ/ml). The coefficients of variation for the scan-rescan ranged from 2.5%-9.2% for the flow components' volume ratio and between 13.5%-17.7% for the kinetic energy. The interval scan results showed higher coefficients of variation with values from 6.2-16.1% for the flow components' volume ratio and 16.9-29.0% for the kinetic energy of the flow components.

CONCLUSION

LV flow components' volume and their associated kinetic energy values are repeatable and stable within a population over time. However, the variability of these measurements in individuals over time is greater than can be attributed to sources of error in the data acquisition and analysis, suggesting that additional physiological factors may influence LV flow measurements.

Longitudinal changes in myocardial T1 and T2 relaxation times related to diffuse myocardial fibrosis in aortic stenosis; before and after aortic valve replacement

BACKGROUND

Diffuse myocardial fibrosis is associated with adverse outcomes, although detection and quantification is challenging. Cardiac MR relaxation times mapping represents a promising imaging biomarker for diffuse myocardial fibrosis.

PURPOSE

To investigate whether relaxation times can detect longitudinal changes in myocardial tissue composition associated with diffuse fibrosis in patients with severe aortic stenosis (AS) before and after aortic valve replacement (AVR).

STUDY TYPE

Prospective longitudinal study.

POPULATION/SUBJECTS/PHANTOM/SPECIMEN/ANIMAL MODEL

Fifteen patients with severe AS.

FIELD STRENGTH/SEQUENCE

3T / 3(3)3(3)5-MOLLI, T2 -GraSE, and 3D-QALAS.

ASSESSMENT

Patients underwent MR examinations at three timepoints: before AVR, as well as 3 and 12 months after AVR. Data from each patient was analyzed in 16 myocardial segments.

STATISTICAL TESTS

The segment-wise T1 and T2 data were analyzed over time after surgery using linear mixed models for repeated measures analysis.

RESULTS

The results showed that T1 relaxation times were significantly (P < 0.05) shorter 3 and 12 months postoperative than preoperative and that the T2 relaxation times were significantly (P < 0.05) longer 3 and 12 months postoperative than preoperative for both 3D and 2D mapping methods. No significant changes were seen between 3 and 12 months postoperative for any of the methods (P = 0.06/0.19 for T1 with 3D-QALAS/MOLLI and P = 0.09/0.25 for T2 with 3D-QALAS/GraSE).

DATA CONCLUSION

We demonstrated that changes in myocardial relaxation times and thus tissue characteristics can be observed within 3 months after AVR surgery. The significant changes in relaxation times from preoperative examinations to the follow-up may be interpreted as a reduction of interstitial fibrosis in the left ventricular wall.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2018.

Age-Related Vascular Changes Affect Turbulence in Aortic Blood Flow

Turbulent blood flow is implicated in the pathogenesis of several aortic diseases but the extent and degree of turbulent blood flow in the normal aorta is unknown. We aimed to quantify the extent and degree of turbulece in the normal aorta and to assess whether age impacts the degree of turbulence. 22 young normal males (23.7 ± 3.0 y.o.) and 20 old normal males (70.9 ± 3.5 y.o.) were examined using four dimensional flow magnetic resonance imaging (4D Flow MRI) to quantify the turbulent kinetic energy (TKE), a measure of the intensity of turbulence, in the aorta. All healthy subjects developed turbulent flow in the aorta, with total TKE of 3–19 mJ. The overall degree of turbulence in the entire aorta was similar between the groups, although the old subjects had about 73% more total TKE in the ascending aorta compared to the young subjects (young = 3.7 ± 1.8 mJ, old = 6.4 ± 2.4 mJ, p < 0.001). This increase in ascending aorta TKE in old subjects was associated with age-related dilation of the ascending aorta which increases the volume available for turbulence development. Conversely, age-related dilation of the descending and abdominal aorta decreased the average flow velocity and suppressed the development of turbulence. In conclusion, turbulent blood flow develops in the aorta of normal subjects and is impacted by age-related geometric changes. Non-invasive assessment enables the determination of normal levels of turbulent flow in the aorta which is a prerequisite for understanding the role of turbulence in the pathophysiology of cardiovascular disease.