Objectives. DOPPLER-CIP aims to determine the optimal noninvasive parameters (myocardial function, perfusion, ventricular blood flow, cell integrity) and methodology (ergometry, echocardiography, scintigraphy, MRI) in a given ischemic substrate that best predicts the impact of an intervention (or the lack thereof) on adverse morphological ventricular remodeling and functional recovery. Moreover, the relative predictive value of each of these parameters, in respect to the cost of extracting this information in order to enable optimization of cost-effectiveness for improved health care, will be determined by this project. Design. DOPPLER-CIP is a multi-center registry study. All patients with ischemic heart disease included in this study undergo at least two noninvasive stress imaging examinations at baseline. The presence/or absence of left ventricular (LV) remodeling will be assessed after a follow-up of 2 years, during which all cardiac events will be registered. Results. 676 patients were included. Currently, baseline data analysis is almost finished and the follow-up is ongoing. Conclusions. After completion, DOPPLER-CIP will provide evidence-based guidelines toward the most effective use of cardiac imaging in the chronically ischemic heart disease patient. The study will generate information, knowledge, and insight into the new imaging methodologies and into the pathophysiology of chronic ischemic heart disease.
Wall shear stress (WSS) distribution in the human aorta is a highly interesting hemodynamic factor for atherosclerosis development. We present a magnetic resonance imaging (MRI) and computational fluid dynamics (CFD) based subject specific WSS estimation method and demonstrate it on a group of nine healthy volunteers (males age 23.6 ± 1.3 years). In all nine subjects, the aortic blood flow was simulated in a subject specific way, where the 3D segmented geometries and inflow profiles were obtained using MRI. No parameter settings were tailored using data from the nine subjects. Validation was performed by comparing CFD gained velocity with magnetic resonance imaging (MRI) velocity measurements. CFD and MRI velocity profiles were comparable, but the temporal variations of the differences during the cardiac cycle were significant. Spatio-temporal analyzes on the WSS distribution showed a strong subject specific influence. Subject specific models are decisive to estimate WSS distribution.
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 similar to 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.
Purpose: Eddy current induced velocity offsets are of concern for accuracy in cardiovascular magnetic resonance (CMR) volume flow quantification. However, currently known theoretical aspects of eddy current behavior have not led to effective guidelines for the optimization of flow quantification sequences. This study is aimed at identifying correlations between protocol parameters and the resulting velocity error in clinical CMR flow measurements in a multi-vendor study. Methods: Nine 1.5T scanners of three different types/vendors were studied. Measurements were performed on a large stationary phantom. Starting from a clinical breath-hold flow protocol, several protocol parameters were varied. Acquisitions were made in three clinically relevant orientations. Additionally, a time delay between the bipolar gradient and read-out, asymmetric versus symmetric velocity encoding, and gradient amplitude and slew rate were studied in adapted sequences as exploratory measurements beyond the protocol. Image analysis determined the worst-case offset for a typical great-vessel flow measurement. Results: The results showed a great variation in offset behavior among scanners (standard deviation among samples of 0.3, 0.4, and 0.9 cm/s for the three different scanner types), even for small changes in the protocol. Considering the absolute values, none of the tested protocol settings consistently reduced the velocity offsets below the critical level of 0.6 cm/s neither for all three orientations nor for all three scanner types. Using multilevel linear model analysis, oblique aortic and pulmonary slices showed systematic higher offsets than the transverse aortic slices (oblique aortic 0.6 cm/s, and pulmonary 1.8 cm/s higher than transverse aortic). The exploratory measurements beyond the protocol yielded some new leads for further sequence development towards reduction of velocity offsets; however those protocols were not always compatible with the time-constraints of breath-hold imaging and flow-related artefacts. Conclusions: This study showed that with current systems there was no generic protocol which resulted into acceptable flow offset values. Protocol optimization would have to be performed on a per scanner and per protocol basis. Proper optimization might make accurate (transverse) aortic flow quantification possible for most scanners. Pulmonary flow quantification would still need further (offline) correction.
A four-dimensional mapping (three spatial dimensions + time) of myocardial strain-rate would help to describe the mechanical properties of the myocardium, which affect important physiological factors such as the pumping performance of the ventricles. Strain-rate represents the local instantaneous deformation of the myocardium and can be calculated from the spatial gradients of the velocity field. Strain-rate has previously been calculated using one-dimensional (ultrasound) or two-dimensional (2-D) magnetic resonance imaging techniques. However, this assumes that myocardial motion only occurs in one direction or in one plane, respectively. This paper presents a method for calculation of the time-resolved three-dimensional (3-D) strain-rate tensor using velocity vector information in a 3-D spatial grid during the whole cardiac cycle. The strain-rate tensor provides full information of both magnitude and direction of the instantaneous deformation of the myocardium. A method for visualization of the full 3-D tensor is also suggested. The tensors are visualized using ellipsoids, which display the principal directions of strain-rate and the ratio between strain-rate magnitude in each direction. The presented method reveals the principal strain-rate directions without a priori knowledge of myocardial motion directions.
The complex three-dimensional structure of the heart muscle (myocardium) has anisotropic, non-linear and time-dependent mechanical properties. During the cardiac cycle, the myocardium undergoes large elastic deformations as a consequence of the active muscle contraction along the muscle fibers and their relaxation, respectively. A four-dimensional (4D) description (three spatial dimensions + time) of the mechanical properties of the myocardium would be of interest in the assessment of myocardial function. Time-resolved 3D phase contrast MRI makes it possible to quantify all three velocity components, which is necessary to as accurately as possible describe the velocities in the heart. The velocity data may be used for investigation of the deformation of the heart and calculation of strain in the myocardial wall. We present a method for estimation of myocardial kinematics using finite elements and 3D time-resolved phase contrast MRI.
Blood flow patterns are closely linked to the morphology and function of the cardiovascular system. These patterns reflect the exceptional adaptability of the cardiovascular system to maintain normal blood circulation under a wide range of workloads. Accurate retrieval and display of flow-related information remains a challenge because of the processes involved in mapping the flow velocity fields within specific chambers of the heart. We review the potentials and pitfalls of current approaches for blood flow visualization, with an emphasis on acquisition, display, and analysis of multidirectional flow. This document is divided into 3 sections. First, we provide a descriptive outline of the relevant concepts in cardiac fluid mechanics, including the emergence of rotation in flow and the variables that delineate vortical structures. Second, we elaborate on the main methods developed to image and visualize multidirectional cardiovascular flow, which are mainly based on cardiac magnetic resonance, ultrasound Doppler, and contrast particle imaging velocimetry, with recommendations for developing dedicated imaging protocols. Finally, we discuss the potential clinical applications and technical challenges with suggestions for further investigations.
While many methods exist for visualising scalar and vector data, visualisation of tensor data is still troublesome. We present a method for visualising second order tensors in three dimensions using a hybrid between direct volume rendering and glyph rendering.
An overview scalar field is created by using three-dimensional adaptive filtering of a scalar field containing noise. The filtering process is controlled by the tensor field to be visualised, creating patterns that characterise the tensor field. By combining direct volume rendering of the scalar field with standard glyph rendering methods for detailed tensor visualisation, a hybrid solution is created.
A combined volume and glyph renderer was implemented and tested with both synthetic tensors and strain-rate tensors from the human heart muscle, calculated from phase contrast magnetic resonance image data. A comprehensible result could be obtained, giving both an overview of the tensor field as well as detailed information on individual tensors.
Aim: The influences on the SNR of DENSE MRI of field strength, receiver coil sensitivity and choice of flip angle strategy have been previously investigated individually. In this study, all of these parameters have been investigated in the same setting, and a mutual comparison of their impact on SNR is presented.
Materials and methods: Ten healthy volunteers were imaged in a 1.5T and a 3T MRI system, using standard 5 or 6 channel cardiac coils as well as 32 channel coils, with four different excitation patterns. Variation of spatial coil sensitivity was assessed by regional SNR analysis.
Results: SNR ranging from 2.8 to 30.5 was found depending on the combination of excitation patterns, coil sensitivity and field strength. The SNR at 3T was 53 ± 26% higher than at 1.5T (p<0.001), whereas spatial differences of 59 ± 26% were found in the ventricle (p<0.001). 32 channel coils provided 52 ± 29% higher SNR compared to standard 5 or 6 channel coils (p<0.001). A fixed flip angle strategy provided an excess of 50% higher SNR in half of the imaged cardiac cycle compared to a sweeping flip angle strategy, and a single phase acquisition provided a six-fold increase of SNR compared to a cine acquisition.
Conclusion: The effect of field strength and receiver coil sensitivity influences the SNR with the same order of magnitude, whereas flip angle strategy can have a larger effect on SNR. Thus, careful choice of imaging hardware in combination with adaptation of the acquisition protocol is crucial in order to realize sufficient SNR in DENSE MRI.
A method to acquire multiple displacement encoded slices within asingle breath hold is presented. Efficiency is improved overconventional Displacement ENcoding with Stimulated Echoes (DENSE) without compromising image quality by read-out of multiple slices inthe same cardiac cycle, thus utilizing the position encoded stimulatedecho available in the whole heart. The method was evaluated bycomparing strain values obtained using the proposed method to strainvalues obtained by conventional separate breath hold single-sliceDENSE acquisitions. Good agreement (Lagrangian E2 strainbias=0.000, 95% limits of agreement ±0.04,root-mean-square-difference 0.02 (9.4% of the mean end-systolic E2)) was found between the methods, indicating that the proposedmethod can replace a multiple breath hold acquisition. Eliminating theneed for multiple breath holds reduces the risk of changes in breathhold positions or heart rate, results in higher patient comfort andfacilitates inclusion of DENSE in a clinical routine protocol.
In multiphase DENSE MRI, available signal decreases during the cardiac cycle due to excitation and T1-relaxation. Because of this, a widely used method is to vary the flip angle to yield constant SNR during the cardiac cycle. Unfortunately, this may lead to significant sacrifices of SNR in the early cardiac phases.
Time-resolved three-dimensional phase-contrast MRI is an important tool for physiological as well as clinical studies of blood flow in the heart and vessels. The application of the technique is, however, limited by the long scan times required. In this work, we investigate the feasibility of using spiral readouts to reduce the scan time of four-dimensional flow MRI without sacrificing quality. Three spiral approaches are presented and evaluated in vivo and in vitro against a conventional Cartesian acquisition. In vivo, the performance of each method was assessed in the thoracic aorta in 10 volunteers using pathline-based analysis and cardiac output analysis. Signal-to-noise ratio and background phase errors were investigated in vitro. Using spiral readouts, the scan times of a four-dimensional flow acquisition of the thoracic aorta could be reduced 23-fold, with no statistically significant difference in pathline validity or cardiac output. The shortened scan time improves the applicability of four-dimensional flow MRI, which may allow the technique to become a part of a clinical workflow for cardiovascular functional imaging.
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 (Pamp;lt;0.0001). The residual volume proportion and residual volume average KE were increased in patients (Pamp;lt;0.0001). Importantly, in a multiple linear regression model to predict the patients 6-minute walk test, the independent predictors were age (beta=-0.3015; P=0.019) and direct-flow average KE (beta=0.280, P=0.035; R-2 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.
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 amp;gt;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.
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-2 = 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.
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 similar to 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.
Aims: 4D flow magnetic resonance imaging (MRI) allows quantitative assessment of left ventricular (LV) function according to characteristics of the dynamic flow in the chamber. Marked abnormalities in flow components volume and kinetic energy (KE) have previously been demonstrated in moderately dilated and depressed LVs compared to healthy subjects. We hypothesized that these 4D flow-based measures would detect even subtle LV dysfunction and remodeling.
Methods and Results: We acquired 4D flow and morphological MRI data from 26 patients with chronic ischemic heart disease with New York Heart Association (NYHA) class I and II and with no to mild LV systolic dysfunction and remodeling, and from 10 healthy controls. A previously validated method was used to separate the LV end-diastolic volume (LVEDV) into functional components: direct flow, which passes directly to ejection, and non-ejecting flow, which remains in the LV for at least 1 cycle. The direct flow and non-ejecting flow proportions of end-diastolic volume and KE were assessed. The proportions of direct flow volume and KE fell with increasing LVEDV-index (LVEDVI) and LVESV-index (LVESVI) (direct flow volume r = -0.64 and r = -0.74, both Pamp;lt;0.001; direct flow KE r = -0.48, P = 0.013, and r = -0.56, P = 0.003). The proportions of non-ejecting flow volume and KE rose with increasing LVEDVI and LVESVI (non-ejecting flow volume: r = 0.67 and r = 0.76, both Pamp;lt;0.001; non-ejecting flow KE: r = 0.53, P = 0.005 and r = 0.52, P = 0.006). The proportion of direct flow volume correlated moderately to LVEF (r = 0.68, P amp;lt; 0.001) and was higher in a sub-group of patients with LVEDVI amp;gt; 74 ml/m(2) compared to patients with LVEDVI amp;lt; 74 ml/m(2) and controls (both Pamp;lt;0.05).
Conclusion: Direct flow volume and KE proportions diminish with increased LV volumes, while non-ejecting flow proportions increase. A decrease in direct flow volume and KE at end-diastole proposes that alterations in these novel 4D flow-specific markers may detect LV dysfunction even in subtle or subclinical LV remodeling.
Patient specific modelling of the blood flow through the human aorta is performed using computational fluid dynamics (CFD) and magnetic resonance imaging (MRI). Velocity patterns are compared between computer simulations and measurements. The workflow includes several steps: MRI measurement to obtain both geometry and velocity, an automatic levelset segmentation followed by meshing of the geometrical model and CFD setup to perform the simulations follwed by the actual simulations. The computational results agree well with the measured data.
Purpose
To correct for displacement artifacts in 3D phase contrast imaging.
Materials and Methods
A 3D phase contrast pulse sequence was modified so that displacements of velocity measurements were restricted to one direction. By applying a postprocessing method, displaced measurements could be traced back to their accurate positions. Flow studies were performed using a phantom that generated flow through a stenosis, directed oblique relative to the phase and frequency encoding directions. Velocity profiles and streamline visualization were used to compare displaced and corrected velocity data to a reference.
Results
Velocity profiles obtained from the original measurement showed skewed profiles due to the displacement artifact, both at close proximity to the orifice as well as further downstream. After correction, concordance with the reference improved considerably.
Conclusion
The displacement artifact, which restricts the accuracy of phase contrast measurements, can be corrected for using the proposed method. Correction of the phase contrast velocity data may improve the accuracy of subsequent flow analysis and visualization.
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. Students 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.
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.
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 5T. 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. Students 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-2 = 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-2 = 0.93-0.95 SGRE; 0.88-0.91 EPI, R-2 = 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. Technical Efficacy Stage 1 J. Magn. Reson. Imaging 2019.
The flow patterns in the human heart are complex and difficult to visualize using conventional two-dimensional (2D) modalities, whether they depict a single velocity component (Doppler echocardiography) or all three components in a few slices (2D phase contrast MRI). To avoid these shortcomings, a temporally resolved 3D phase contrast technique was used to derive data describing the intracardiac velocity fields in normal volunteers. The MRI data were corrected for phase shifts caused by eddy currents and concomitant gradient fields, with improvement in the accuracy of subsequent flow visualizations. Pathlines describing the blood pathways through the heart were generated from the temporally resolved velocity data, starting from user-specified locations and time frames. Flow trajectories were displayed as 3D particle traces, with simultaneous demonstration of morphologic 2D slices. This type of visualization is intuitive and interactive and may extend our understanding of dynamic and previously unrecognized patterns of intracardiac flow.
Purpose: To assess turbulent kinetic energy (TKE) within the left ventricle (LV) of healthy subjects using novel 4D flow MRI methods and to compare TKE values to those from a spectrum of patients with dilated cardiomyopathy (DCM).
Methods: 4D flow and morphological MRI-data were acquired in 11 healthy subjects and 9 patients with different degrees of diastolic dysfunction. TKELV was calculated within the LV at each diastolic time frame. At peak early (E) and late (A) diastolic filling, the TKELV was compared to transmitral peak velocity, LV diameter and mitral annular diameter.
Results: In the majority of all subjects, peaks in TKELV could be observed at E and A. Peak TKELV at E was not different between the groups, and correlated with mitral annular dimensions. Peak TKELV at A was higher in DCM patients compared to healthy subjects, and was related to LV diameter and transmitral velocity.
Conclusions: In normal LVs, TKE values are low. Values are highest during early diastole, and diminish with increasing LV size. In a heterogeneous group of DCM patients, late diastolic TKE values are higher than in healthy subjects. Kinetic energy loss due to elevated late diastolic TKE may reflect inefficient flow in dilated LVs.
The impact of left bundle branch block (LBBB) related mechanical dyssynchrony on left ventricular (LV) diastolic function remains unclear. 4D flow cardiovascular magnetic resonance (CMR) has provided reliable markers of LV dysfunction: reduced volume and kinetic energy (KE) of the portion of LV inflow which passes directly to outflow (Direct Flow) has been demonstrated in failing hearts compared to normal hearts. We sought to investigate the impact of mechanical dyssynchrony on diastolic function by comparing 4D flow in myopathic LVs with and without LBBB. CMR data were acquired at 3 T in 22 heart failure patients; 11 with LBBB and 11 without LBBB matched according to several demographic and clinical parameters. An established 4D flow analysis method was used to separate the LV end-diastolic (ED) volume into functional flow components based on the blood's timing and route through the heart cavities. While the Direct Flow volume was not different between the groups, the KE possessed at ED was lower in LBBB patients (P = 0.018). Direct Flow entering the LV during early diastolic filling possessed less KE at ED in LBBB patients compared to non-LBBB patients, whereas no intergroup difference was observed during late filling. Pre-systolic KE of LV Direct Flow was reduced in patients with LBBB compared to matched patients with normal conduction. These intriguing findings propose that 4D flow specific measures can serve as markers of LV mechanical dyssynchrony in heart failure patients, and could possibly be investigated as predictors of response to cardiac resynchronization therapy.
Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2019.
Purpose: To explore the use of MR-estimated turbulence quantities for the assessment of turbulent flow effects on the vessel wall. Methods: Numerical velocity data for two patient-derived models was obtained using computational fluid dynamics (CFD) for two physiological flow rates. The four-dimensional (4D) Flow MRI measurements were simulated at three different spatial resolutions and used to investigate the estimation of turbulent wall shear stress (tWSS) using the intravoxel standard deviation (IVSD) of velocity and turbulent kinetic energy (TKE) estimated near the vessel wall. Results: Accurate estimation of tWSS using the IVSD is limited by the spatial resolution achievable with 4D Flow MRI. TKE, estimated near the wall, has a strong linear relationship to the tWSS (mean R(2=)0.84). Near-wall TKE estimates from MR simulations have good agreement to CFD-derived ground truth (mean R-2=0.90). Maps of near-wall TKE have strong visual correspondence to tWSS. Conclusion: Near-wall estimation of TKE permits assessment of relative maps of tWSS, but direct estimation of tWSS is challenging due to limitations in spatial resolution. Assessment of tWSS and near-wall TKE may open new avenues for analysis of different pathologies. (C) 2016 International Society for Magnetic Resonance in Medicine
Purpose: To examine methods for visualizing and quantifying flow stasis in abdominal aortic aneurysms (AAA) using 4D Flow MRI. Methods: Three methods were investigated: conventional volumetric residence time (VRT), mean velocity analysis (MVA), and particle travel distance analysis (TDA). First, ideal 4D Flow MRI data was generated using numerical simulations and used as a platform to explore the effects of noise and background phase-offset errors, both of which are common 4D Flow MRI artifacts. Error-free results were compared to noise or offset affected results using linear regression. Subsequently, 4D Flow MRI data for thirteen (13) subjects with AAA was acquired and used to compare the stasis quantification methods against conventional flow visualization. Results: VRT (R-2 = 0.69) was more sensitive to noise than MVA (R-2 = 0.98) and TDA (R-2 = 0.99) at typical noncontrast signal-to-noise ratio levels (SNR = 20). VRT (R-2 = 0.14) was more sensitive to background phase-offsets than MVA (R-2 = 0.99) and TDA (R-2 = 0.96) when considering a 95% effective background phase-offset correction. Qualitatively, TDA outperformed MVA (Wilcoxon p amp;lt; 0.005, mean score improvement 1.6/5), and had good agreement (median score 4/5) with flow visualizations. Conclusion: Flow stasis can be quantitatively assessed using 4D Flow MRI. While conventional residence time calculations fail due to error accumulation as a result of imperfect measured velocity fields, methods that do not require lengthy particle tracking perform better. MVA and TDA are less sensitive to measurement errors, and TDA generates results most similar to those obtained using conventional flow visualization.