- Open Access
Blood flow dynamics with four-dimensional flow cardiovascular magnetic resonance in patients with aortic stenosis before and after transcatheter aortic valve replacement
Journal of Cardiovascular Magnetic Resonance volume 23, Article number: 81 (2021)
Pre- and post-procedural hemodynamic changes which could affect adverse outcomes in aortic stenosis (AS) patients who undergo transcatheter aortic valve replacement (TAVR) have not been well investigated. Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) enables accurate analysis of blood flow dynamics such as flow velocity, flow pattern, wall shear stress (WSS), and energy loss (EL). We sought to examine the changes in blood flow dynamics of patients with severe AS who underwent TAVR.
We examined 32 consecutive severe AS patients who underwent TAVR between May 2018 and June 2019 (17 men, 82 ± 5 years, median left ventricular ejection fraction 61%, 6 self-expanding valve), after excluding those without CMR because of a contraindication or inadequate imaging from the analyses. We analyzed blood flow patterns, WSS and EL in the ascending aorta (AAo), and those changes before and after TAVR using 4D flow CMR.
After TAVR, semi-quantified helical flow in the AAo was significantly decreased (1.4 ± 0.6 vs. 1.9 ± 0.8, P = 0.002), whereas vortical flow and eccentricity showed no significant changes. WSS along the ascending aortic circumference was significantly decreased in the left (P = 0.038) and left anterior (P = 0.033) wall at the basal level, right posterior (P = 0.011) and left (P = 0.010) wall at the middle level, and right (P = 0.012), left posterior (P = 0.019) and left anterior (P = 0.028) wall at the upper level. EL in the AAo was significantly decreased (15.6 [10.8–25.1 vs. 25.8 [18.6–36.2]] mW, P = 0.012). Furthermore, a significant negative correlation was observed between EL and effective orifice area index after TAVR (r = − 0.38, P = 0.034).
In severe AS patients undergoing TAVR, 4D flow CMR demonstrates that TAVR improves blood flow dynamics, especially when a larger effective orifice area index is obtained.
Transcatheter aortic valve replacement (TAVR) has emerged as a minimally invasive treatment option for patients with severe symptomatic aortic stenosis (AS) who cannot tolerate surgical aortic valve replacement (SAVR) and those at intermediate–high surgical risk, and is expanding to younger, lower-surgical-risk patients [1,2,3].
In patients with AS, chronic pressure overload evolves to left ventricular (LV) remodeling, characterized by cardiac muscle hypertrophy and interstitial collagen deposition, ultimately leading to systolic and/or diastolic dysfunction [4, 5]. While it is known that patients who undergo TAVR can improve diastolic dysfunction and undergo LV reverse remodeling, residual diastolic dysfunction and lack of LV reverse remodeling after TAVR may adversely affect the prognosis . The ascending aorta (AAo), as an immediate buffer of LV ejection, plays important roles in determining LV afterload . In fact, increased arterial stiffness and hemodynamic abnormalities in the AAo increase LV oxygen demand and cardiac work, leading to exacerbation of heart failure (HF) [8, 9]. Thus, non-invasive pre- and post-procedural hemodynamic assessment in the AAo may help evaluate the effects of TAVR.
Time-resolved three-dimensional (3D) phase-contrast cardiovascular magnetic resonance (CMR), also known as four-dimensional (4D) flow CMR, is a novel blood flow imaging technique that allows accurate visualization and quantification of aortic blood flow dynamics . Notably, 4D flow CMR can also calculate wall shear stress (WSS), which is the force acting in a tangential direction to the blood vessel lumen surface due to blood flow, and energy loss (EL), which is the kinetic energy lost in blood flow due to frictional force. Several studies have demonstrated various blood flow patterns assessed by this modality in healthy subjects [11, 12], patients with AS , and patients post-SAVR or -TAVR [14,15,16]. Recently, we have reported a successful case of assessment of blood flow and WSS in the AAo before and after TAVR using 4D flow CMR . Blood flow dynamic parameters obtained with 4D flow CMR, such as WSS and EL in the AAo, are expected to be an important indicator of LV afterload, which would be related to impaired LV expansion and progression of LV remodeling . However, it is uncertain how blood flow dynamics are impaired in patients with AS and how TAVR alters blood flow dynamic parameters obtained by 4D flow CMR. Furthermore, it is unclear which parameters related to blood flow dynamics in the AAo could be improved after TAVR. Accordingly, we sought to investigate the changes in blood flow dynamics before and after TAVR using 4D flow CMR to clarify their functional significance.
This was a single-center, observational, prospective study that included consecutive patients who underwent TAVR with a diagnosis of severe symptomatic AS according to the updated guidelines . The study protocol was approved by the Ethics Committee of Hokkaido University Hospital (018-0223 and 019-0090). The investigation conformed with the principles outlined in the Declaration of Helsinki. All patients gave informed consent to participate in the study.
We originally screened 48 consecutive patients who underwent TAVR between May 2018 and June 2019. All patients enrolled in this study were diagnosed with severe AS. They met any of the following criteria on transthoracic echocardiography: a peak aortic velocity (Vmax) ≥ 4.0 m/s, a mean transaortic pressure gradient (mean PG) ≥ 40 mmHg, an aortic valve area (AVA) ≤ 1.0 cm2, or an aortic valve area index (AVAI) ≤ 0.6 cm2/m2. Among patients with severe AS, those satisfying any of Vmax ≥ 5.0 m/s, mean PG ≥ 60 mmHg, AVA ≤ 0.6 cm2, or AVAI ≤ 0.4 cm2/m2 criteria were called very severe AS . Of these, patients without CMR because of a contraindication (n = 9) to CMR and those with inadequate imaging unsuitable for analysis because of artifacts (n = 5) were excluded. Two patients declined to participate in this study. Ultimately, 32 patients were included in this study (Fig. 1). In addition, 12 control subjects without significant AS or aortic regurgitation were recruited to compare the blood flow dynamic parameters obtained by 4D flow CMR to the pre- or post-TAVR groups.
All cases were discussed at a meeting of the multidisciplinary Heart Team and were indicated to undergo TAVR as recommended by current guidelines . TAVR was performed through a transfemoral, transapical or transsubclavian route depending on pre-procedural vascular assessment. In most patients, TAVR was performed with transesophageal echocardiographic guidance under general anesthesia. Pre- and post-dilatation was performed at the operator’s discretion.
Echocardiography was performed within 2 weeks before and after TAVR, respectively. LV systolic or diastolic function were assessed by echocardiography before and after TAVR: LV end-diastolic dimension (LVEDD), LV end-systolic dimension (LVESD), LV ejection fraction (LVEF), the ratio between early and late diastolic transmitral flow velocity (E/A), the ratio of maximal early diastolic filling wave velocity to maximal early diastolic myocardial velocity (E/e’), and left atrial (LA) volume index (LAVI). LVEF was measured from the apical 4- and 2-chamber images using the biplane method of disks. LA volume was measured from standard apical 4-chamber views at end-systole just before mitral valve opening. The biplane method of disks was used to calculate LA volume. LV mass (LVM) was calculated according to the following formula:
LAVI and LVM index (LVMI) were calculated by dividing LA volume or LV mass by body surface area of patients, respectively. Relative wall thickness (RWT) was defined as 2 times inferolateral wall thickness divided by the LVEDD. LV remodeling was assessed based on LVEDD, LVMI and RWT.
After TAVR, the effective orifice area (EOA) was estimated using the LV outflow tract diameter and the LV outflow tract time velocity integral measured immediately proximal to the stent, as recommended previously . The effective orifice area index (EOAI) was calculated with the EOA of the prosthesis divided by the body surface area of patients.
Cardiovascular magnetic resonance
CMR imaging using a 3 T CMR scanner (Achieva TX, Philips Healthcare, Best, The Netherlands) with a 32-channel phased-array receiver torso-cardiac coil was performed in patients before TAVR (median interval, 24 [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39] days) and after TAVR (median interval, 6 [4,5,6] days). Retrospective electrocardiogram gating and a respiratory navigator placed on the lung–liver interface (7 mm window size) were used. In all patients, LV and aortic blood flow dynamics were evaluated using 4D flow CMR measuring 3D blood flow velocities with full volumetric coverage of the heart and the thoracic aorta. 4D flow CMR data were acquired in a sagittal oblique 3D data including the entire heart and thoracic aorta without contrast agent. We acquired 4D flow sequence with 3D spatial encoding. We scanned 4 or 5 mm slice thickness without overlapping using 3D spatial encoding. After data acquisition, we reconstructed images with slice gap 2 or 2.5 mm. The scan parameters were: echo time = 1.73 ms, repetition time = 3.2 ms, flip angle α = 10°, field of view 400 × 400 mm, matrix 256 × 229, in-plane spatial resolution 1.56 × 1.75 mm2, slice thickness 4 mm and slice overlapping 2 mm, or slice thickness 5 mm and slice overlapping 2.5 mm (depending on the patient's body size), temporal resolution 12 phases/cardiac cycle, k-space segmentation factor 6, Sensitivity encoding factor R = 3. Velocity encoding timing was TR-interleaved. Partial k-space coverage methods and k–t undersampling were not used. Velocity encoding (VENC) was set to 2.5–8.0 m/s individually based on the peak blood flow velocity in the AAo with a secured margin. The acquisition time of 4D flow CMR was approximately 8–20 min.
Maxwell correction, eddy current correction, and gradient non-linearity correction methods were performed by the CMR scanner. Acquired data including phase-contrast 3-axis cine images, magnitude images and steady-state free procession cine images were processed using commercially available software (iTFlow, Cardio Flow Design Inc., Tokyo, Japan), which visualized the cardiovascular geometry and blood flow [22, 23]. Moreover, WSS and EL were calculated using the software .
The blood flow pattern in the AAo was evaluated as previously described . Three analysis planes were positioned perpendicular to the aortic wall at the level of the sinotubular junction (Slice 1), mid AAo (Slice 2), and proximal to the brachiocephalic trunk (Slice 3) (Fig. 2a).
Blood flow pattern in the AAo was evaluated using streamlines, which is the instantaneous velocity field at a specific temporal phase. Two readers (one cardiologist and one radiologist) simultaneously observed the streamlines of each phase from the patient's left front. The systolic blood flow pattern was semi-quantitatively evaluated as three grades for vortical flow and helical flow, respectively. Discordant cases were evaluated by a third reader (a radiologist experienced in cardiovascular imaging). Vortical flow was defined as revolving particles around a point within the vessel with a rotation direction deviating by more than 90° from the physiological flow direction . Helical flow was defined as regional fluid circulation along the longitudinal axis of the vessel, thereby creating a corkscrew-like motion . Vortex and helix strengths were graded in three categories: 1 = none (none or almost none), 2 = moderate (obvious. between 1 and 3), 3 = marked (mainstream). The detailed distinction between moderate and marked strengths was as follows. Vortical flow from the brachiocephalic artery bifurcation to the sinotubular junction was determined as marked, and smaller vortical flow was determined as moderate. Helical flow with laminar flow was classified as moderate, and helical flow without laminar flow was classified as marked. Examples of grading for vortical and helical flow are shown in Fig. 2c.
Peak velocity blood flow eccentricity in the mid AAo (Slice 2) during systole was semi-quantitatively evaluated as three grades: 1 = none (if high velocity systolic flow was centrally focused, occupying the majority of the vessel lumen), 2 = mild (if high velocity systolic flow occupied between one- and two-thirds of the vessel lumen), 3 = marked (if high velocity systolic flow occupied one-third or less of the vessel lumen) (Fig. 2d) .
Mean and maximum blood flow velocity and WSS in the AAo were calculated using the above-mentioned software tool (iTFlow, Cardio Flow Design Inc.) as described previously . In brief, the surface of the segmented aorta at the moment of peak systole is represented by connected triangles, and the WSS is defined as the spatial gradient of the 3D velocity vector in each edge point, perpendicular to the vessel wall. This algorithm has been established previously by Potters et al. . The WSS was calculated from the formula below (Formula 1).
μ: blood viscosity (μ = 0.004 Pa s), v: blood flow velocity, y: distance from the wall.
First, an anatomical segmentation of the AAo is obtained using iTFlow software. The AAo was defined as the aortic segment between the sinotubular junction and the origin of the brachiocephalic trunk before TAVR. After TAVR, the AAo was defined as the segment between the first circumferential area of the AAo not susceptible to metal-induced artifacts and the origin of the brachiocephalic trunk. Subsequently, we determined the 3D WSS over the complete AAo and displayed it on a color-coded map. Peak WSS (WSS peak) used the highest WSS value in this color-coded map of the entire AAo. As the average WSS (WSS average), the mean value of WSS of the entire AAo in the phase in which the peak was recorded is used (Fig. 2b). Furthermore, to investigate changes in WSS in each region, calculation of WSS were performed for 12 segments along the aortic circumference for each analysis plane (Slice 1 to Slice 3) (Fig. 2a and e) . By specifying the position of a mesh sphere whose size can be changed, it is possible to calculate WSS on any surface in contact with the sphere. The WSS value in each segment was the average of the calculated values at multiple contact points (31.9 ± 11.7) between the aorta reconstructed in 3D and the sphere. The assessment of the WSS is shown in Fig. 2a, b, and e.
EL was calculated from the spatial velocity gradient of blood flow and blood viscosity according to the formula below (Formula 2) [22, 25]. EL across the region of interest (left ventricle or AAo) was calculated for each of the 12 phases/cardiac cycles and averaged for systolic and diastolic phases, respectively. Regions of interest were selected for each slice of the coronal image. The heart chamber from the mitral valve to the aortic valve was defined as the region of the LV, and the aorta from the aortic valve to the brachiocephalic artery was defined as the region of the AAo.
μ: blood viscosity (μ = 0.004 Pa s), x: horizontal direction of the phase image, u: horizontal direction component of the blood velocity vector.
Continuous variables are presented as the mean ± standard deviation when normally distributed, and as the median and interquartile range when not normally distributed. Comparisons between the TAVR group and the control group were performed by the Mann–Whitney U-test for continuous variables and by the chi-squared test for dichotomous variables. The changes in the blood flow dynamic parameters after TAVR were evaluated by the Wilcoxon signed-rank test. Inter-reader agreement for blood flow patterns was assessed using quadratic weighted kappa statistics (along with their standard errors). Linear regression models were constructed to assess the relationship between WSS average or WSS peak in the AAo and AVAI before TAVR (or EOAI after TAVR). The same analysis was performed to confirm the relationship between EL and AVAI before TAVR (or EOAI after TAVR). Another linear regression model was constructed to confirm the relationship between WSS or EL in the AAo and the parameters of the LV function or remodeling, as well as the relationship between vortical flow, helical flow, blood flow eccentricity, WSS or EL and aortic size and pre- and post-TAVR aortic pressure gradient. All tests were two tailed, and a P value < 0.05 was considered statistically significant. All analyses were performed using Stata MP64 (version 15, Stata Corp, College Station, Texas, USA).
The baseline characteristics of all 32 subjects (82 ± 5 years, 17 (53%) men) are shown in Table 1. Mean transaortic PG was 53 ± 21 mmHg and AVAI was 0.5 ± 0.1 cm2/m2. Median LVEF was 61% (interquartile range 54–70%), and mean stroke volume index was 47 ± 12 ml/m2. Twenty-eight (88%) patients underwent transfemoral TAVR, and a balloon-expanding valve was used in 26 (81%) patients. Eighteen (56%) were classified as very severe AS. There was no significant difference in hemodynamic changes between severe AS and very severe AS groups (Additional file 5: Table S1).
The changes in the echocardiographic parameters before and after TAVR are shown in Table 2. EOAI was significantly increased, and the peak and mean transaortic PG were significantly decreased after TAVR. Conversely, there were no significant changes in the echocardiographic data for assessment of LV systolic or diastolic function and remodeling. None of our patients had severe prosthesis–patient mismatch (EOAI < 0.65 cm2/m2), or paravalvular leakage of 3 degrees or more, which have an adverse prognostic effect after TAVR.
The characteristics of the control group are shown in Additional file 5: Table S2; there were no significant differences in age and sex between the TAVR group and the control group. Cardiac function was comparable between the groups, with similar LV dimensions and LVEF. However, the TAVR group had a higher LVMI than the control group.
Blood flow patterns and flow eccentricity in ascending aorta
Two representative cases which showed significant improvement in vortical or helical blood flow in the AAo following TAVR procedure are demonstrated in Fig. 3 and supplementary videos (see Additional files: 1, 2, 3, and 4 [Videos S1, S2, S3, and S4]).
Of the 64 examinations, the number of discordant cases was 16 for vortex, 11 for helix, and 12 for eccentricity. The quadratic weighted kappa coefficients for inter-observer agreement for blood flow patterns were as follows: k = 0.654; standard error = 0.097 for vortex; k = 0.751, standard error = 0.098 for helix; and k = 0.627, standard error = 0.111 for eccentricity. The reproducibility of each parameter for the blood flow pattern was substantial.
Semi-quantified flow eccentricity values were significantly higher in the pre- and post-TAVR groups compared to the control group. Helical blood flow values were higher in the pre-TAVR group compared to the control group; however, there was no significant difference after TAVR. Furthermore, there was no significant difference in the degree of vortical blood flow between the pre- or post-TAVR groups and the control group (Additional file 5: Table S3).
After TAVR, semi-quantified helical blood flow in the AAo significantly decreased, whereas vortical blood flow and flow eccentricity showed no significant changes (Fig. 4).
Wall shear stress in ascending aorta
WSS average and WSS peak in the entire AAo were significantly higher in the pre- and post-TAVR groups compared to the control group (Additional file 5: Table S3).
The changes of WSS average and WSS peak in the entire AAo after TAVR are shown in Fig. 5. WSS average and WSS peak significantly decreased from 6.7 [6.1–8.4] Pa to 6.0 [5.4–7.0] Pa (P = 0.009) and from 52.0 [45.6–62.2] Pa to 47.5 [38.4–53.1] Pa (P = 0.023), respectively. The relationship between WSS at pre- or post-TAVR and AVAI or EOAI is shown in Additional file 5: Table S4. There was no significant correlation between WSS average or WSS peak and AVAI before TAVR. There was also no significant correlation between WSS average or WSS peak and EOAI after TAVR.
Figure 6 shows peak WSS in 12 segments along the aortic circumference for analysis in plane slices 1 to 3 before and after TAVR. WSS was significantly decreased in the left and left anterior wall at the basal level, right posterior and left wall at the middle level, and right, left posterior and left anterior wall at the upper level.
Energy loss in ascending aorta and left ventricle
The systolic EL in the AAo was significantly higher in the pre-TAVR group compared to the control group; however, there was no significant difference after TAVR (Additional file 5: Table S3).
The changes of EL in the AAo and LV before and after TAVR is shown in Fig. 7. EL in the AAo during the systolic phase was significantly decreased (15.6 [10.8–25.1 vs. 25.8 [18.6–36.2]] mW, P = 0.012), while that in the LV was unchanged after TAVR. Furthermore, a significant negative correlation was observed between EL in the AAo during the systolic phase and EOAI after TAVR (r = − 0.38, P = 0.034, Fig. 8). However, there was no significant correlation between EL in the AAo during the systolic phase and AVAI before TAVR (Additional file 5: Table S4). A significant positive correlation was observed between systolic EL in the AAo and mean transaortic pressure gradient after TAVR, and between diastolic EL in the AAo and the AAo diameter after TAVR.
Differences between the self- and the balloon-expanding bioprostheses
Six subjects received self-expanding TAVR bioprostheses. There were no significant differences on blood flow dynamic parameters in 4D flow CMR after TAVR between self-expanding and balloon-expanding bioprostheses (n = 26) (Additional file 5: Table S5). However, the self-expanding bioprosthesis group had significantly greater decrease in WSS average after TAVR than the balloon-expanding bioprosthesis group. In addition, systolic EL in the AAo tended to be decreased in the self-expanding bioprosthesis group compared with the balloon-expanding bioprosthesis group (Additional file 5: Table S6).
This study is an initial report evaluating the changes of blood flow dynamics using 4D flow CMR before and after TAVR. The major findings of this study were that (1) helical blood flow, WSS and EL in the AAo were significantly decreased, and (2) EL in the AAo was negatively correlated with EOAI after TAVR.
Using 3D phase-contrast cine CMR, Kilner et al. observed that helical blood flow develops in the AAo and extends towards the aortic arch in healthy subjects . Segadal and Matre also observed a bidirectional velocity profile in late systole and early diastole, featuring retrograde velocity (i.e., vortical blood flow) along the inner wall and sustained antegrade flow along the outer wall . Previous studies have demonstrated the role of helical blood flow in the AAo, showing that relatively coherent swirling of blood might avoid excessive dissipation of energy by limiting flow instability in arteries [27, 28]. Although helical blood flow plays a positive physiological role in enhancing blood flow transport, the flow pattern changes when aortic morphological changes occur, such as aortic enlargement or AS [13, 29]. Vortical flow also affects the elasticity of the aorta, protruding aortic atheroma, and diastolic coronary flow. Notably, helical blood flow, vortical blood flow, and blood eccentricity assessed by 4D flow CMR became stronger as the disease progressed in patients with AS . Furthermore, a higher grade of helical blood flow would lead to higher EL and WSS distribution in the AAo .
In this study, the degree of helical blood flow was higher in patients with AS than in the non-AS control group; however, there was no significant difference in vortical blood flow. Furthermore, we found a significant decrease in helical blood flow, but not in vortical blood flow, after TAVR. Generally, as the diameter of the AAo and age increase, vortical blood flow in the AAo increases . The loss of arterial wall elasticity and the large difference in blood flow velocity between the outside and inside due to dilation of the AAo are thought to cause an increase of vortical blood flow. The AAo diameter of the control subjects was larger than that of the patients with AS but might not produce a significant difference in vortical blood flow. This suggests that vortical blood flow may be greatly affected by the properties of the AAo. On the other hand, helical blood flow also developed in the non-dilated AAo  and its grade was significantly associated with the severity of AS resulting from aortic valve calcific fusion and reduced mobility . Therefore, helical rotation in blood flow rather than vortical flow might be reduced after TAVR.
In this study, the eccentricity observed before TAVR remained considerable after TAVR. Florian et al. reported that blood flow eccentricity on 4D flow CMR remains strong even after stented bioprosthesis SAVR . This finding may reflect the blood flow dynamic consequences caused by the rigid opening mechanism and the asymmetric orifices. This may be explained by differences in the opening mechanism and in the sensibility to co-factors such as the angle between the open cusp and the annular plane . Furthermore, Trauzeddel et al. also reported that patients who undergo TAVR show marked eccentricity compared to healthy controls . The characteristics of the implanted stented bioprosthesis may explain the presence of eccentricity after TAVR. Furthermore, it is known that the angle and height of placement for the annular plane differs depending on the type of TAVR. Although eccentricity may vary, we did not have sufficient sample size for a subgroup analysis to determine the differences in blood flow dynamics between balloon- and self-expanding bioprostheses.
Wall shear stress
Disorganized blood flow patterns such as helical and vortical flow due to abnormal valve opening resulted in markedly altered regional and global WSS in the AAo . Furthermore, elevated WSS in the AAo has been attributed to changes in endothelial function and vascular remodeling in patients with a bicuspid aortic valve and degenerative AS [34, 35]. Our patients showed regional and global decreases in WSS after TAVR, but it did not reach the level of non-AS control subjects. TAVR reduced WSS in the right posterior wall of the middle AAo and the anterior wall of the distal AAo, which are strongly aggravated by increased counterclockwise helical blood flow that passes through the stenotic aortic valve. This finding is consistent with the findings of von Knobelsdorff-Brenkenhoff et al. . They reported a decrease in helical blood flow and eccentricity after TAVR, which altered WSS localization and reduced global WSS in the AAo, suggesting that TAVR decreases WSS and reduces the risk of future aortic disease such as aortic dissection and aneurysm. In contrast, Farag et al. demonstrated that patients who underwent TAVR showed an increase in WSS in the AAo when compared to age- and gender- matched controls . Importantly, WSS decreased after TAVR but did not reach a normal level in our study population. The mechanisms for higher WSS in TAVR patients compared to control subjects could be explained as follows. First, TAVR prosthetic valves are implanted inside a calcified, native aortic valve annulus. This causes the TAVR valve to have a smaller EOAI than that of control subjects. Indeed, in our study, EOAI after TAVR was 1.1 ± 0.2 cm2/m2, which was smaller than the normal EOAI of 1.8 cm2/m2. This may provide an explanation as to why WSS post-TAVR did not recover to control level. Second, LV hypertrophy progresses to overcome the pressure load of the stenotic aortic valve. Therefore, stroke volume is higher in patients who undergo TAVR than in healthy subjects .
Our study showed that patients with AS had significantly higher systolic EL in the AAo than non-AS controls. Furthermore, systolic EL in the AAo was significantly reduced after TAVR. Using 4D flow CMR to conduct non-invasive estimations of EL might be useful for assessing hemodynamics in patients with AS, because it more closely reflects the amount of LV energy that is lost during systole due to the obstruction derived from the calcified aortic valve. EL in the AAo is a key indicator of LV afterload, which can worsen LV function . In fact, it is higher in patients with AS than in healthy subjects, suggesting that high EL reflects an increased LV afterload due to abnormal blood flow patterns in the AAo [18, 32]. These findings indicate that TAVR provides efficient blood flow dynamics and can reduce LV afterload.
Interestingly, we also found that EL had a significant negative correlation with EOAI after TAVR. However, there was no significant correlation between EL and AVAI before TAVR in patients with AS. One possible explanation is that most patients had an AVAI < 0.6 cm2/m2 before TAVR, and the range of AVAI values was small. The patient prosthesis mismatch is described as the mismatch between the hemodynamics derived from a prosthetic valve and the patient’s requirements for adequate cardiac output. Patient prosthesis mismatch is defined by the EOAI, and approximately 37% of patients after TAVR have moderate or severe patient prosthesis mismatch . Severe patient prosthesis mismatch after TAVR leads to subsequent adverse events such as worsening of HF or death . These adverse events may be associated with increased EL, which is a result of a small EOAI and consistent abnormal blood flow patterns after TAVR. In patients with severe AS who undergo TAVR, 4D flow CMR helps to visualize and assess the changes of blood flow dynamics.
Abnormal blood flow pattern in the AAo increases LV afterload as energy is dissipated by frictional losses associated with various blood flow phenomena and the aortic wall [18, 32]. As a consequence, accurate assessment of blood flow patterns in patients undergoing TAVR may allow estimation of the risk of TAVR-related adverse events represented by worsening HF due to increased LV afterload. It is noteworthy that TAVR prosthetic valves have been reported to obtain a larger EOAI than those in SAVR. Furthermore, self-expanding valves could achieve a larger EOAI than balloon-expanding valves, and it is necessary to carefully consider which valve to use for patients with a narrow annulus [3, 37]. Obtaining a larger EOAI may lead to an improved prognosis by avoiding the increase in EL. Moreover, Bahlmann et al. reported that EL was superior to EOA in predicting adverse clinical outcomes in patients with moderate to severe AS . Therefore, further studies investigating whether impaired blood flow dynamics including flow patterns, WSS and EL assessed by 4D flow CMR are associated with adverse events after TAVR, or these assessments are useful for selecting the type of prosthetic valve in AS patients undergoing TAVR, are warranted.
There are several potential limitations that should be acknowledged. First, the spatial and temporal resolution of 4D flow CMR in this study was relatively low. In the setting of spatial and temporal resolution, there is a trade-off between scan time and the accuracy of parameters such as the flow rate, WSS, and oscillatory shear index estimation . Decreasing the spatial and temporal resolution in order to shorten the scan time affects the accuracy of flow quantification and visualization adversely  and leads to underestimation of WSS [11, 39, 40]. As the voxel size increases, the accuracy of WSS estimates decreases . Since we favored the accuracy of WSS, we used an optimized, small-size in-plane spatial resolution (1.6 × 1.8 mm2). To minimize the influence of anisotropic voxels, we set the slice gap to 2–2.5 mm. In our preliminary study with young healthy volunteers (n = 5), the scan time was much longer for 20 phases (12.3 ± 1.9 min) than for 12 phases (6.6 ± 0.6 min). In elderly patients with severe heart disease, longer scan duration was impractical, and we decided to reduce the temporal resolution. The lower temporal resolution might imply the loss of possible information contained in other systolic phases . However, in this study, all acquisitions were made with the same imaging parameters and analyzed with the same methodology both pre- and post-TAVR. Since our study focused on changes in parameters between pre- and post-TAVR, underestimation of WSS might be compensated. Compressed sensing (CS), which exploits the inherent compressibility of CMR data, has been combined with parallel imaging to achieve even higher acceleration rates, enabling 4D flow CMR scan times to become clinically feasible [43, 44]. Ma et al. applied CS to 4D flow imaging of the thoracic aorta and achieved a scan time of 2 min . However, CS 4D flow CMR has been reported to underestimate peak flow and velocity. Further investigation using CS-based protocols is needed. Second, since our CMR unit did not allow multiple VENC settings, we set a single VENC value according to the velocity in the head/foot direction. However, this setting could measure the majority of high velocity compared to anterior/posterior or right/left directions. Third, the Reynolds number could provide more information on the impact of turbulent flow in the aorta from the stenotic aortic valve and how it affects the calculation of the parameters using 4D flow CMR. However, we were unable to actually calculate the Reynolds number due to the inability to measure aortic inlet diameter by 4D flow CMR and to accurately evaluate flow velocity after TAVR due to signal loss caused by metal artifacts of the implanted valve. Therefore, we calculated the substitute Reynolds number for each case using the diameter of the AAo on computed tomography, mean aortic velocity on echocardiography, and blood kinematic viscosity ν = 3.3 × 10−6. Substitute Reynolds number decreased significantly from 9741 [8992–12,574) before TAVR to 4358 [3688–4790] after TAVR (P < 0.001). From these numbers, because the WSS calculation was done in a turbulent state in both pre- and post-TAVR, we think that the comparison of WSS is valid. Furthermore, because the WSS is calculated by directly deriving the blood flow velocity gradient from the flow velocity near the wall in the software used for this study, the calculation is applicable in both laminar and turbulent flows. Generally, in 4D flow CMR, turbulent flow decreases the signal and causes an underestimation of the flow velocity . Due to the limitations of the spatial resolution of 4D flow CMR, the reliability of the absolute value of WSS is controversial, however, we believe that it is meaningful in comparing relative changes before and after the intervention. We recognize that the failure to assess the actual effect of the Reynolds number on the changes of blood flow dynamics in 4D flow CMR is a major limitation of the present study. Fourth, although the blood flow patterns were evaluated after thorough confirmation of the evaluation methods by two radiologists and one cardiologist, these evaluations may have differed among them. This is considered to be a limitation of the visual semi-quantitative method. In this study, those assessments were reproducible (but not yet well standardized) and susceptible to observer bias. Fifth, we did not have sufficient sample size for a subgroup analysis to determine the differences in blood flow dynamics between balloon- and self-expanding bioprostheses. Although our findings showed that the self-expanding bioprostheses might reduce the WSS and EL in the AAo more than the balloon-expanding bioprostheses, further studies are warranted to confirm these findings. Finally, the number of adverse events and the sample size were too small to evaluate changes in the LV function or remodeling. Therefore, we were unable to confirm the relationship between impaired blood flow dynamics (obtained from 4D flow CMR) and worse clinical outcomes. Accordingly, further larger-scale studies are necessary to confirm the prognostic significance of the changes in blood flow dynamics in patients who undergo TAVR.
In severe AS patients who undergo TAVR, 4D flow CMR is useful for visualizing and assessing the changes of blood flow dynamics. TAVR improves blood flow dynamics especially when a larger EOAI is obtained.
Availability of data and materials
The datasets used during the current study are available from the corresponding author on reasonable request.
Aortic valve area
Aortic valve area index
Cardiovascular magnetic resonance
Effective orifice area
Effective orifice area index
Left atrium/left atrial
Left atrial volume index
Left ventricle/left ventricular
Left ventricular end-diastolic dimension
Left ventricular ejection fraction
Left ventricular end-systolic dimension
Left ventricular mass
Left ventricular mass index
Relative wall thickness
Surgical aortic valve replacement
Transcatheter aortic valve replacement
- Vmax :
Wall shear stress
Nishimura RA, Otto CM, Bonow RO, Carabello BA, Erwin JP 3rd, Fleisher LA, et al. 2017 AHA/ACC Focused Update of the 2014 AHA/ACC Guideline for the Management of Patients With Valvular Heart Disease: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. Circulation. 2017;135:e1159–95.
Mack MJ, Leon MB, Thourani VH, Makkar R, Kodali SK, Russo M, et al. Transcatheter aortic-valve replacement with a balloon-expandable valve in low-risk patients. N Engl J Med. 2019;380:1695–705.
Popma JJ, Deeb GM, Yakubov SJ, Mumtaz M, Gada H, O’Hair D, et al. Transcatheter aortic-valve replacement with a self-expanding valve in low-risk patients. N Engl J Med. 2019;380:1706–15.
Hatani T, Kitai T, Murai R, Kim K, Ehara N, Kobori A, et al. Associations of residual left ventricular and left atrial remodeling with clinical outcomes in patients after aortic valve replacement for severe aortic stenosis. J Cardiol. 2016;68:241–7.
Barone-Rochette G, Pierard S, De Meester de Ravenstein C, Seldrum S, Melchior J, Maes F, et al. Prognostic significance of LGE by CMR in aortic stenosis patients undergoing valve replacement. J Am Coll Cardiol. 2014;64:144–54.
Genereux P, Head SJ, Hahn R, Daneault B, Kodali S, Williams MR, et al. Paravalvular leak after transcatheter aortic valve replacement: the new Achilles’ heel? A comprehensive review of the literature. J Am Coll Cardiol. 2013;61:1125–36.
Chirinos JA. Ventricular-arterial coupling: invasive and non-invasive assessment. Artery Res. 2013;7:2–14.
Chen CH, Nakayama M, Nevo E, Fetics BJ, Maughan WL, Kass DA. Coupled systolic-ventricular and vascular stiffening with age: implications for pressure regulation and cardiac reserve in the elderly. J Am Coll Cardiol. 1998;32:1221–7.
Redheuil A, Kachenoura N, Bollache E, Yu WC, Opdahl A, Decesare A, et al. Left ventricular and proximal aorta coupling in magnetic resonance imaging: aging together? Am J Physiol Heart Circ Physiol. 2019;317:H300–7.
Itatani K, Miyazaki S, Furusawa T, Numata S, Yamazaki S, Morimoto K, et al. New imaging tools in cardiovascular medicine: computational fluid dynamics and 4D flow MRI. Gen Thorac Cardiovasc Surg. 2017;65:611–21.
Stalder AF, Russe MF, Frydrychowicz A, Bock J, Hennig J, Markl M. Quantitative 2D and 3D phase contrast MRI: optimized analysis of blood flow and vessel wall parameters. Magn Reson Med. 2008;60:1218–31.
Kilner PJ, Yang GZ, Mohiaddin RH, Firmin DN, Longmore DB. Helical and retrograde secondary flow patterns in the aortic arch studied by three-directional magnetic resonance velocity mapping. Circulation. 1993;88:2235–47.
von Knobelsdorff-Brenkenhoff F, Karunaharamoorthy A, Trauzeddel RF, Barker AJ, Blaszczyk E, Markl M, et al. Evaluation of aortic blood flow and wall shear stress in aortic stenosis and its association with left ventricular remodeling. Circ Cardiovasc Imaging. 2016;9:e004038.
von Knobelsdorff-Brenkenhoff F, Trauzeddel RF, Barker AJ, Gruettner H, Markl M, Schulz-Menger J. Blood flow characteristics in the ascending aorta after aortic valve replacement–a pilot study using 4D-flow MRI. Int J Cardiol. 2014;170:426–33.
Trauzeddel RF, Lobe U, Barker AJ, Gelsinger C, Butter C, Markl M, et al. Blood flow characteristics in the ascending aorta after TAVI compared to surgical aortic valve replacement. Int J Cardiovasc Imaging. 2016;32:461–7.
Farag ES, Vendrik J, van Ooij P, Poortvliet QL, van Kesteren F, Wollersheim LW, et al. Transcatheter aortic valve replacement alters ascending aortic blood flow and wall shear stress patterns: a 4D flow MRI comparison with age-matched, elderly controls. Eur Radiol. 2019;29:1444–51.
Komoriyama H, Tsuneta S, Oyama-Manabe N, Kamiya K, Nagai T. Four-dimensional flow magnetic resonance imaging visualizes significant changes in flow pattern and wall shear stress in the ascending aorta after transcatheter aortic valve implantation in a patient with severe aortic stenosis. Eur Heart J Cardiovasc Imaging. 2020;21:21.
Barker AJ, van Ooij P, Bandi K, Garcia J, Albaghdadi M, McCarthy P, et al. Viscous energy loss in the presence of abnormal aortic flow. Magn Reson Med. 2014;72:620–8.
Nishimura RA, Otto CM, Bonow RO, Carabello BA, Erwin JP 3rd, Guyton RA, et al. 2014 AHA/ACC guideline for the management of patients with valvular heart disease: a report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines. J Thorac Cardiovasc Surg. 2014;148:e1–132.
Tribouilloy C, Rusinaru D, Bohbot Y, Maréchaux S, Vanoverschelde JL, Enriquez-Sarano M. How should very severe aortic stenosis be defined in asymptomatic individuals? J Am Heart Assoc. 2019;8:011724.
Zoghbi WA, Chambers JB, Dumesnil JG, Foster E, Gottdiener JS, Grayburn PA, et al. Recommendations for evaluation of prosthetic valves with echocardiography and doppler ultrasound: a report From the American Society of Echocardiography’s Guidelines and Standards Committee and the Task Force on Prosthetic Valves, developed in conjunction with the American College of Cardiology Cardiovascular Imaging Committee, Cardiac Imaging Committee of the American Heart Association, the European Association of Echocardiography, a registered branch of the European Society of Cardiology, the Japanese Society of Echocardiography and the Canadian Society of Echocardiography, endorsed by the American College of Cardiology Foundation, American Heart Association, European Association of Echocardiography, a registered branch of the European Society of Cardiology, the Japanese Society of Echocardiography, and Canadian Society of Echocardiography. J Am Soc Echocardiogr. 2009;22:975–1014.
Miyazaki S, Itatani K, Furusawa T, Nishino T, Sugiyama M, Takehara Y, et al. Validation of numerical simulation methods in aortic arch using 4D Flow MRI. Heart Vessels. 2017;32:1032–44.
Takei Y, Itatani K, Miyazaki S, Shibasaki I, Fukuda H. Four-dimensional flow magnetic resonance imaging analysis before and after thoracic endovascular aortic repair of chronic type B aortic dissection. Interact Cardiovasc Thorac Surg. 2019;28:413–20.
Potters WV, van Ooij P, Marquering H, vanBavel E, Nederveen AJ. Volumetric arterial wall shear stress calculation based on cine phase contrast MRI. J Magn Reson Imaging. 2015;41:505–16.
Stugaard M, Koriyama H, Katsuki K, Masuda K, Asanuma T, Takeda Y, et al. Energy loss in the left ventricle obtained by vector flow mapping as a new quantitative measure of severity of aortic regurgitation: a combined experimental and clinical study. Eur Heart J Cardiovasc Imaging. 2015;16:723–30.
Segadal L, Matre K. Blood velocity distribution in the human ascending aorta. Circulation. 1987;76:90–100.
Morbiducci U, Ponzini R, Rizzo G, Cadioli M, Esposito A, Montevecchi FM, et al. Mechanistic insight into the physiological relevance of helical blood flow in the human aorta: an in vivo study. Biomech Model Mechanobiol. 2011;10:339–55.
Stonebridge PA, Buckley C, Thompson A, Dick J, Hunter G, Chudek JA, et al. Non spiral and spiral (helical) flow patterns in stenosis. In vitro observations using spin and gradient echo magnetic resonance imaging (MRI) and computational fluid dynamic modeling. Int Angiol. 2004;23:276–83.
Markl M, Draney MT, Hope MD, Levin JM, Chan FP, Alley MT, et al. Time-resolved 3-dimensional velocity mapping in the thoracic aorta: visualization of 3-directional blood flow patterns in healthy volunteers and patients. J Comput Assist Tomogr. 2004;28:459–68.
Ha H, Kim GB, Kweon J, Lee SJ, Kim YH, Kim N, et al. The influence of the aortic valve angle on the hemodynamic features of the thoracic aorta. Sci Rep. 2016;6:32316.
Frydrychowicz A, Berger A, Munoz Del Rio A, Russe MF, Bock J, Harloff A, et al. Interdependencies of aortic arch secondary flow patterns, geometry, and age analysed by 4-dimensional phase contrast magnetic resonance imaging at 3 Tesla. Eur Radiol. 2012;22:1122–30.
Garcia J, Barker AJ, Markl M. The role of imaging of flow patterns by 4D flow MRI in aortic stenosis. JACC Cardiovasc Imaging. 2019;12:252–66.
Yoganathan AP, He Z, Casey JS. Fluid mechanics of heart valves. Annu Rev Biomed Eng. 2004;6:331–62.
Hope MD, Hope TA, Crook SE, Ordovas KG, Urbania TH, Alley MT, et al. 4D flow CMR in assessment of valve-related ascending aortic disease. JACC Cardiovasc Imaging. 2011;4:781–7.
Farag ES, van Ooij P, Planken RN, Dukker KCP, de Heer F, Bouma BJ, et al. Aortic valve stenosis and aortic diameters determine the extent of increased wall shear stress in bicuspid aortic valve disease. J Magn Reson Imaging. 2018;48:522–30.
Herrmann HC, Daneshvar SA, Fonarow GC, Stebbins A, Vemulapalli S, Desai ND, et al. Prosthesis-patient mismatch in patients undergoing transcatheter aortic valve replacement: from the STS/ACC TVT Registry. J Am Coll Cardiol. 2018;72:2701–11.
Hahn RT, Pibarot P, Stewart WJ, Weissman NJ, Gopalakrishnan D, Keane MG, et al. Comparison of transcatheter and surgical aortic valve replacement in severe aortic stenosis: a longitudinal study of echocardiography parameters in cohort A of the PARTNER trial (placement of aortic transcatheter valves). J Am Coll Cardiol. 2013;61:2514–21.
Bahlmann E, Gerdts E, Cramariuc D, Gohlke-Baerwolf C, Nienaber CA, Wachtell K, et al. Prognostic value of energy loss index in asymptomatic aortic stenosis. Circulation. 2013;127:1149–56.
Cibis M, Potters WV, Gijsen FJ, Marquering H, van Ooij P, vanBavel E, et al. The effect of spatial and temporal resolution of cine phase contrast MRI on wall shear stress and oscillatory shear index assessment. PLoS ONE. 2016;11:e0163316.
Dyverfeldt P, Bissell M, Barker AJ, Bolger AF, Carlhall CJ, Ebbers T, et al. 4D flow cardiovascular magnetic resonance consensus statement. J Cardiovasc Magn Reson. 2015;17:72.
Petersson S, Dyverfeldt P, Ebbers T. Assessment of the accuracy of MRI wall shear stress estimation using numerical simulations. J Magn Reson Imaging. 2012;36:128–38.
Rodríguez-Palomares JF, Dux-Santoy L, Guala A, Kale R, Maldonado G, Teixidó-Turà G, et al. Aortic flow patterns and wall shear stress maps by 4D-flow cardiovascular magnetic resonance in the assessment of aortic dilatation in bicuspid aortic valve disease. J Cardiovasc Magn Reson. 2018;20:28.
Tariq U, Hsiao A, Alley M, Zhang T, Lustig M, Vasanawala SS. Venous and arterial flow quantification are equally accurate and precise with parallel imaging compressed sensing 4D phase contrast MRI. J Magn Reson Imaging. 2013;37:1419–26.
Rich A, Potter LC, Jin N, Liu Y, Simonetti OP, Ahmad R. A Bayesian approach for 4D flow imaging of aortic valve in a single breath-hold. Magn Reson Med. 2019;81:811–24.
Ma LE, Markl M, Chow K, Huh H, Forman C, Vali A, et al. Aortic 4D flow MRI in 2 minutes using compressed sensing, respiratory controlled adaptive k-space reordering, and inline reconstruction. Magn Reson Med. 2019;81:3675–90.
Kadbi M, Negahdar M, Cha JW, Traughber M, Martin P, Stoddard MF, et al. 4D UTE flow: a phase-contrast MRI technique for assessment and visualization of stenotic flows. Magn Reson Med. 2015;73:939–50.
The authors are grateful to Kinya Ishizaka for his support with the analyses of 4D flow CMR.
This work was supported by grants from the Japan Society for the Promotion of Science, Tokyo, Japan (KAKENHI No. JP19K17510, Dr. Kamiya), the Takeda Science Foundation (Dr. Nagai), and a Bayer Scholarship for Cardiovascular Research, the Japan Cardiovascular Research Foundation (Dr. Nagai).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Hokkaido University Hospital (018-0223 and 019-0090). All subjects provided informed consent for research participation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Video S1. In a patient, marked vortical flow was observed before TAVR (MP4 675 KB)
Additional file 2: Video S2. Vortical flow decreased after TAVR. In another patient
Additional file 3: Video S3. Marked helical flow was observed before TAVR
Additional file 4: Video S4. Helical flow disappeared after TAVR. TAVR indicates transcatheter aortic valve replacement
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Komoriyama, H., Kamiya, K., Nagai, T. et al. Blood flow dynamics with four-dimensional flow cardiovascular magnetic resonance in patients with aortic stenosis before and after transcatheter aortic valve replacement. J Cardiovasc Magn Reson 23, 81 (2021). https://doi.org/10.1186/s12968-021-00771-y
- Four-dimensional flow cardiovascular magnetic resonance
- Aortic stenosis
- Transcatheter aortic valve replacement
- Blood flow dynamics