- 1School of Biomedical Engineering, Universidad de Valparaíso, Valparaíso, Chile
- 2Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile
- 3Millennium Institute for Intelligent Healthcare Engineering, iHEALTH, Santiago, Chile
- 4Millennium Nucleus in Cardiovascular Magnetic Resonance, Cardio MR, Santiago, Chile
- 5Department of Electrical Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
- 6Department of Cardiology, Hospital Universitari Vall d'Hebron, CIBER-CV, Vall d'Hebron Institut de Recerca (VHIR), Barcelona, Spain
- 7Department of Mechanical Engineering, Universidad Técnica Federico Santa María, Santiago, Chile
- 8Department of Structural and Geotechnical Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
- 9Institute for Biological and Medical Engineering, Schools of Engineering, Medicine and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile
- 10Department of Radiology, Schools of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile
Background and Purpose: Prognostic models based on cardiovascular hemodynamic parameters may bring new information for an early assessment of patients with bicuspid aortic valve (BAV), playing a key role in reducing the long-term risk of cardiovascular events. This work quantifies several three-dimensional hemodynamic parameters in different patients with BAV and ranks their relationships with aortic diameter.
Materials and Methods: Using 4D-flow CMR data of 74 patients with BAV (49 right-left and 25 right-non-coronary) and 48 healthy volunteers, aortic 3D maps of seventeen 17 different hemodynamic parameters were quantified along the thoracic aorta. Patients with BAV were divided into two morphotype categories, BAV-Non-AAoD (where we include 18 non-dilated patients and 7 root-dilated patients) and BAV-AAoD (where we include the 49 patients with dilatation of the ascending aorta). Differences between volunteers and patients were evaluated using MANOVA with Pillai's trace statistic, Mann–Whitney U test, ROC curves, and minimum redundancy maximum relevance algorithm. Spearman's correlation was used to correlate the dilation with each hemodynamic parameter.
Results: The flow eccentricity, backward velocity, velocity angle, regurgitation fraction, circumferential wall shear stress, axial vorticity, and axial circulation allowed to discriminate between volunteers and patients with BAV, even in the absence of dilation. In patients with BAV, the diameter presented a strong correlation (> |+/−0.7|) with the forward velocity and velocity angle, and a good correlation (> |+/−0.5|) with regurgitation fraction, wall shear stress, wall shear stress axial, and vorticity, also for morphotypes and phenotypes, some of them are correlated with the diameter. The velocity angle proved to be an excellent biomarker in the differentiation between volunteers and patients with BAV, BAV morphotypes, and BAV phenotypes, with an area under the curve bigger than 0.90, and higher predictor important scores.
Conclusions: Through the application of a novel 3D quantification method, hemodynamic parameters related to flow direction, such as flow eccentricity, velocity angle, and regurgitation fraction, presented the best relationships with a local diameter and effectively differentiated patients with BAV from healthy volunteers.
Introduction
Bicuspid aortic valve (BAV) is the most common congenital heart defect (1), with an estimated prevalence between 1 and 2%, increasing in white population and men (3:1) (2). Ascending aorta (AAo) dilation is presented in approximately 50% of patients with BAV, and it is associated with the increased risk of aortic dissection, rupture, and sudden death (3, 4). The most common BAV leaflet fusion phenotypes involve the right-left cusps (BAV-RL), and right-non-coronary cups (BAV-RN), with the prevalence of around 80 and 17% (4, 5), respectively. Genetic, histological, mechanical, and hemodynamic factors related to aortopathy in patients with BAV are still poorly understood (6–9). Currently, preventive aortic surgery is indicated when the diameter of the AAo is larger than 50 [mm] in patients with any of the following risk factors: aortic coarctation, systemic hypertension, a family history of aortic dissection, or rapid aortic growth (>3–5 mm/year) in experienced hands (10). However, these indications are debatable since some aortic events occur with aortic diameters below the suggested threshold (11). Therefore, there is a need for a better understanding of the mechanisms that influence the progression of these structural changes, which may allow the development of prognostic models for risk assessment, the indication of surgical correction, and pre- and post-operative monitoring (12, 13).
Cardiovascular hemodynamic parameters quantified using 4D-flow cardiac magnetic resonance (CMR) are emerging as the essential biomarkers in the early diagnosis of cardiovascular diseases, bringing new insights about complex flows as in patients with BAV (6, 12, 14). Recent studies have provided strong evidence that altered blood flow hemodynamics and wall shear stress (WSS) in the AAo of patients with BAV are associated with histological and proteolytic changes in the aortic wall, which may induce aortic remodeling (6, 9, 15). Moreover, altered WSS has been related to aortic wall disruption (12, 15), and the separated axial (WSSA) and circumferential (WSSC) components of WSS with abnormal flow eccentricity and leaflet fusion phenotype and extent (6, 16–19). Furthermore, several other hemodynamic parameters, such as flow eccentricity, circulation, vorticity, and helicity density, have also been reported in patients with BAV. However, they have mainly been assessed in a limited number of 2D planes (14, 20) with possible unappreciation of important aspects. Some studies have shown the application of three-dimensional (3D) WSS and its association with the valvular dysfunction (21) and the elastic fiber thinning (15), but the WSS quantification has been most often limited to its magnitude, neglecting its axial and circumferential components, at least in 3D applications. Other studies have shown in 3D the relationship between the absolute local normalized helicity (22) and energy loss (23) and the aortic dilation in patients with BAV. Considering many abnormal flow descriptors proposed, it is of utmost importance to identify those deserving special attention for the follow-up of patients with BAV.
A comprehensive methodology can help to identify those parameters related to aortic dilation in patients with BAV. Previously, we have developed a seamless computational framework to obtain several 3D quantitative parameters, which have been validated in phantoms and different cohorts of patients including aortic dissection (24, 25) and transposition of the great arteries (26, 27). This study aimed to compare quantitative 3D hemodynamic parameters between healthy volunteers (HVs) and patients with BAV and their relationships with aortic dilation in clinically relevant subgroups of patients with BAV. We hypothesize that there are differences in hemodynamic parameters between clinically relevant patients with BAV subgroups (morphotypes and phenotypes), such as those with a non-dilated ascending aorta (BAV-Non-AAoD), which includes non-dilated (BAV-NonD) and root dilated (BAV-RootD), with dilation of the AAo (BAV-AAoD), BAV-RL, BAV-RN, and the group of volunteers.
Methods
Study Population
A total of 49 patients with BAV-RL and 25 BAV-RN with AAo diameters ≤ 50 mm and no severe valvular disease (aortic regurgitation ≤ III, maximum aortic valve velocity <3 m/s by echocardiography) were consecutively and prospectively recruited between June 2014 and December 2015 at the Hospital Universitari Vall d'Hebron (Barcelona, Spain). Inclusion criteria were as follows: age >18 years, no connective tissue disorders, no aortic coarctation or other congenital heart diseases, no previous aortic surgery or aortic valve replacement, and no contraindication for CMR. A total of 48 HVs matched for sex, weight, height, body surface area (BSA), and stroke volume were also included. The local ethics committee approved the study, and informed consent was obtained from all participants.
Cardiovascular Magnetic Resonance Protocol
4D-flow CMR data were obtained in a clinical GE 1.5T Signa scanner (GE Healthcare, Waukesha, WI, USA) using the Vastly undersampled Isotropic Projection Reconstruction (VIPR) technique (28, 29). The volumetric acquisition included the entire thoracic aorta and was performed with retrospective ECG-gating during free-breathing and without administration of an endovenous contrast agent. Acquisitions were made using the following parameters: field of view (FOV) 400 mm × 400 mm × 400 mm, voxel size of 2.5 mm × 2.5 mm × 2.5 mm, flip angle 8°, repetition time 4.2–6.4 ms, echo time 1.9–3.7 ms, and velocity encoding 200 cm/s. This dataset was reconstructed with a temporal resolution that ranged between 21–32 ms.
Aortic Diameters and Aortic Morphotype
The two-dimensional balanced steady-state free precession (bSSFP) cine CMR images were used to assess BAV morphotype and aortic diameters as described (16). Briefly, the three aortic root cusp-to-commissure diameters were measured using double-oblique cine images at the level of the aortic root at end-diastole, and the maximum value was retained for the analysis. Similarly, AAo diameter was measured at end-diastole by double-oblique cine CMR at the level of the pulmonary artery bifurcation. To determine the presence of aortic root or ascending dilation, aortic diameters were adjusted with a logarithm transformation to set the z-score for both sinuses (zsinus) and AAo (zAAo) accounting for sex, age, and BSA as described by Campens et al. (30). Using a z-score cutoff value for the definition of aortic dilation of two standard errors of the estimate, patients were categorized according to the aorta segment predominantly or exclusively involved in dilation according to Della Corte's classification (31). Thus, patients were classified as non-dilated (zsinus ≤ 2 and zAAo ≤ 2), root-dilated morphotype (zsinus>2 and zsinus>zAAo), or AAo-dilated morphotype (zAAo>2 and zAAo>zsinus). Given the small number of patients (seven patients) with rood-dilated morphotype that we have, the groups of patients with BAV were divided only into two morphotypes categories, BAV-Non-AAoD (where we include 18 non-dilated patients and seven root-dilated patients) and BAV-AAoD (where we include the 49 patients with dilatation of the ascending aorta).
3D Quantification of Hemodynamic Parameters and Diameter
The hemodynamic parameters and the diameters in the entire thoracic aorta were obtained based on a finite-element method described previously (32–36) (Figure 1A). To provide a ground for location-specific comparisons, 16 regions (Figure 1D) were semi-automatically selected. Those regions were delimited by anatomical landmarks, regions 1 to 4 = AAo (between the Valsalva level and the brachiocephalic trunk), regions 5 to 8 = aortic arch (AArch, between the brachiocephalic trunk and the isthmus level), regions 9 to 12 = proximal descending aorta (pDAo, between the isthmus level and the level of Valsalva), and regions 13 to 16 = distal descending aorta (dDAo, between the Valsalva level and the diaphragmatic level).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/885338/fcvm-09-885338-HTML/image_m/fcvm-09-885338-g001.jpg)
Figure 1. Summary of the quantification process. (A) From the 4D-flow CMR acquisition, a semiautomatic segmentation of the aorta was generated and transformed into a tetrahedral mesh. Then, 4D-flow CMR velocity values were interpolated to each node of the mesh, and (B,C) 3D maps of 17 hemodynamics parameters and the diameter were calculated using finite elements. (D) A total of 16 different regions of the aorta were created to compare the data between volunteers and patients.
In each region, the mean value of 17 hemodynamic parameters and the aortic diameter were calculated (Figures 1B,C) for each cardiac phase. Peak systolic values were retained for comparison [averaged at peak systole using 1 time-frame before and two time-frames after to reduce noise in the data (16)], except for regurgitation fraction (35) and oscillatory shear index (OSI) calculated along the entire cardiac cycle. The detailed information about the quantification of each parameter is provided in Supplementary Figures S1–S8.
Statistical Analysis
The global statistical differences between healthy volunteers and patients with BAV were evaluated using a multivariate analysis of variance (MANOVA), with Pillai's trace statistic, whereas the local statistical differences were evaluated using the Mann–Whitney U test. To determine which parameters were more relevant to classify BAV morphotypes and phenotypes with respect to healthy volunteers, the receiver operating characteristic (ROC) curves and the minimum redundancy maximum relevance (MRMR) algorithm were used. Moreover, the Spearman's correlation coefficient was used to find the relationship between the hemodynamic parameters and the maximum AAo diameter for all patients with BAV together and for each dilation morphotype and fusion phenotype separately. The Spearman's correlation values (S) were discretized in three different groups: good (0.5 < |S| < =0.7), strong (0.7 < |S| < =0.9), and excellent (0.9 < |S| < =1). The statistical analysis was performed using SPSS Statistics (version 25.0 IBM SPSS, Chicago, IL).
Results
Demographics
A total of 74 patients with BAV (43 men, age 45[37–61] years) and 48 HV (29 men, age 41[30–49] years) were included in the study. Demographic and clinical variables of patients with BAV and HV are shown in Table 1. HVs were matched to patients with BAV in terms of sex, weight, height, BSA, and stroke volume. RL fusion phenotype was present in 66% of patients with BAV. According to the dilation morphotype, 34% patients with BAV (7 BAV root-dilated and 18 BAV non-dilated, or 19 BAV-RL and 6 BAV-RN) were classified as non-ascending aorta dilated (BAV-Non-AAoD) and 66% patients with BAV (19 BAV-RN and 30 BAV-RL) as AAo-dilated (BAV-AAoD).
Hemodynamic Parameters and Diameter vs. Aortic Morphotypes and Phenotypes
The three-dimensional maps of the different parameters are shown in Figure 2 for one representative volunteer and two representative patients (one BAV-Non-AAoD and one BAV-AAoD). The phenotype was not included in this figure due to the limitation of selecting a representative morphotype case in each phenotype group. Differences between volunteers and patients with BAV in the AAo and part of the AArch are visible for several parameters: diameter, eccentricity, backward velocity, velocity angle, regurgitation fraction, WSSC, and axial circulation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/885338/fcvm-09-885338-HTML/image_m/fcvm-09-885338-g002.jpg)
Figure 2. Three-dimensional maps of a representative HV and patients with BAV. Columns represent each analyzed parameter (Figure 1C). First row shows a HV, second row a non-dilated BAV (BAV-Non-AAoD) patient, and third row a patient with BAV with ascending aorta dilation (BAV-AAoD).
From the multivariate analysis test MANOVA, we obtained p-values from Pillai's trace of <0.001 for the HV vs. BAV-All, HV vs. BAV-Non-AAoD, HV vs. BAV-AAoD, HV vs. BAV-RL, and HV vs. BAV-RN. This test showed significant differences between the HV and all BAV groups in all parameters tested. In Figure 3, the comparison between the different BAV groups and HV is shown for each parameter (rows) and aortic region (columns). Most of the analyzed parameters showed significantly higher values in patients with BAV than HV in the AAo and the proximal AArch (regions 1 to 6). In contrast, in the DAo, most of the significant differences resulted from lower values in BAV groups. Eccentricity, backward velocity, velocity angle, regurgitation fraction, WSSC, vorticity, axial vorticity, and axial circulation were consistently greater in patients with BAV than in HV, even in the absence of clinically significant dilation. All parameters (except velocity) showed significant differences in more than one segment for almost all analyzed cases. Velocity, forward velocity, WSS, and WSSA showed lower values in patients with BAV compared to HV.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/885338/fcvm-09-885338-HTML/image_m/fcvm-09-885338-g003.jpg)
Figure 3. Statistical differences between groups: (A) HV vs. all BAV; (B) HV vs. BAV-RL; (C) HV vs. BAV-RN; (D) HV vs. BAV-Non-AAoD; (E) HV vs. BAV-AAoD. Yellow marks highlight statistically significant differences between compared groups (p-value <0.05). The minus sign in yellow boxes indicates the value was lower in BAV than in volunteers.
We also observed similar results for fusion phenotypes (Figures 3B,C) compared to HV. We found that the significant differences in patients with BAV-RN are more concentrated between the ascending aorta and aortic arch, in comparison with BAV-RL, probably, this is influenced by the dilation of the aorta at the proximal part of the aortic arch in patients with BAV-RL, and the descending aorta of these patients show only a few parameters with significant differences. The velocity angle and regurgitation fraction were significantly different for most regions in BAV-RN and BAV-RL. Nevertheless, BAV-RL showed more differences for parameters related to the flow turbulence (parameters from 12 to 17) in the descending aorta with lower values for patients with BAV.
In Supplementary Figure S9, we evaluated the differences between morphotypes and phenotypes groups. BAV-AAoD presented significant differences compared to the BAV-Non-AAoD for velocity, forward velocity, velocity angle, WSS, WSSA, and kinetic energy in the entire aorta. Other parameters such as backward velocity, regurgitation fraction, and vorticity are mostly significant in the ascending aorta. Interestingly, some parameters, which are related to the rotation of the flow, are only significant in the proximal part of the aortic arch, as axial circulation, helicity density, viscous dissipation, and energy loss. Comparing BAV-RL and RN phenotypes groups, we found significant differences in velocity angle, axial vorticity, and axial circulation in more than one region.
Table 2 shows the mean and standard deviation values of each parameter in the AAo (averaged between the regions 1 to 4) for patients with HV and BAV. Patients with BAV presented higher values in most of the parameters in comparison with HV group, whereas velocity, forward velocity, WSS, WSSA, and kinetic energy were lower in patients with BAV. In the AArch, pDAo, and dDAo, the velocity, forward velocity, WSS, WSSA, and the parameters related to the turbulence such as vorticity, axial vorticity, axial circulation, helicity density, viscous dissipation, energy loss, and kinetic energy were lower than in volunteers (Supplementary Tables S1–S3).
Analysis of ROC Curves and MRMR Algorithm
To identify the capacity of each parameter to differentiate between BAV and HV, receiver operating characteristic (ROC) curves were calculated. In the AAo, the area under the ROC curves (Supplementary data online Supplementary Table S4) was >0.8 for diameter, eccentricity, backward velocity, velocity angle, regurgitation fraction, WSSC, axial vorticity, and axial circulation (refer to Figure 4). Of note, the velocity angle was the best performing parameter, showing a specificity close to 99% to classify BAV-AAoD from volunteers, which was reduced to 90% in patients with BAV-Non-AAoD. Similar performances were obtained for eccentricity, backward velocity, WSSC, and axial circulation. In the other aortic sections (AArch, pDAo, and dDAo), only one parameter showed an area under the ROC curves >0.8 (Supplementary data online Supplementary Tables S5–S7), the most characteristic being the regurgitation fraction in the AArch sections.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/885338/fcvm-09-885338-HTML/image_m/fcvm-09-885338-g004.jpg)
Figure 4. ROC curves for the AAo; (A) HV vs. all BAV; (B) HV vs. BAV-RL; (C) HV vs. BAV-RN; (D) HV vs. BAV-Non-AAoD; (E) HV vs. BAV-AAoD. The parameters of the legend are the nine parameters with bigger area under the curve (AUC) values.
In Supplementary Figure S10, we show the ROC curves for the AAo between morphotypes and phenotypes groups. For patients with BAV-RN vs. BAV-RL, we obtained AUC values lower than 0.63, being the velocity angle the greater value. These results are in concordance with the few significant differences found in Supplementary Figure S9A. Comparing BAV-Non-AAoD vs. BAV-AAoD, excluding diameter and forward velocity, whose values are attributed to the greater area of the segment, the parameter with the highest AUC (0.77) was also the velocity angle.
The MRMR algorithm was applied in each section (AAo, AArch, pDAo, and dDAo), to assess similarities with the ROC curves. In the AAo, the velocity angle, axial circulation, regurgitation fraction, diameter, eccentricity, circumferential WSS, backward velocity, forward velocity, and axial vorticity were the parameters with a greater importance predictor score for most of the cases (HV vs. BAV-All, HV vs. BAV-RL, HV vs. BAV-RN, and HV vs. BAV-AAoD). For HV vs. BAV-Non-AAoD, the velocity angle was the parameter with a greater importance predictor score (Supplementary data online Supplementary Figure S11; Supplementary Table S8). Additionally, in Supplementary data online Supplementary Tables S9–11, we show the predictor important scores for the AArch, pDAo, and dDAo, respectively.
The MRMR was also evaluated between morphotypes and phenotypes only for the AAo (Supplementary Figure S12). The comparison of BAV-RN and BAV-RL shows only three parameters with a higher importance prediction score (vorticity, axial vorticity, and axial circulation), and two of them are in concordance with the parameter found in the Supplementary Figure S9A. For BAV-Non-AAoD and BAV-AAoD, we found that the diameter, backward velocity, WSS, forward velocity, regurgitation fraction, velocity angle, vorticity, and axial vorticity were the parameters with a higher importance prediction score.
Correlation Matrices of the Hemodynamic Parameters
In Figure 5, Spearman's correlation coefficient values between the diameter and each hemodynamic parameter in the AAo (the value was averaged between regions 1 to 4) are shown for patients with BAV. Considering all patients with AV together, the aortic diameter presented a strong correlation with forward velocity and velocity angle, the regurgitation fraction, WSS, WSSA, and vorticity presented a good correlation with the diameter. For patients with BAV-RN showed a good correlation for forward velocity and velocity angle. Patients with BAV-RL showed more parameters correlated with the diameter that the other cases, the forward velocity, velocity angle, and WSSA, which are correlated strongly with the diameter, the velocity, backward velocity, regurgitation fraction, WSS, vorticity, axial vorticity, and kinetic energy, have a good correlation with the diameter in this group. For patients with BAV-AAoD, we found a good correlation between the diameter and the forward velocity, WSSA, and vorticity. Finally, patients with BAV-Non-AAoD showed a good correlation for the regurgitation fraction and WSSA.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/885338/fcvm-09-885338-HTML/image_m/fcvm-09-885338-g005.jpg)
Figure 5. Spearman's correlation coefficient values for patients with BAV, BAV-RN, BAV-RL, BAV-AAoD, and BAV-Non-AAoD are reported. These values were grouped as good correlation (red box, Spearman's correlation between |0.5| and |0.7|), strong correlation (blue box, Spearman's correlation between |0.7| and |0.9|), and excellent correlation (green box, Spearman's correlation above |0.9|). The minus sign indicates a negative correlation.
Discussion
One of the strengths of our study is that it provides a comprehensive analysis of a novel quantification tool to obtain several 3D aortic maps of different hemodynamic parameters in a cohort of patients with BAV. Moreover, we presented an extensive comparison with HV, identifying a subset of high-performing parameters. Several hemodynamic parameters showed significant differences between HV and patients with BAV in the AAo. Among them, eccentricity, backward velocity, velocity angle, regurgitation fraction, WSSC, axial vorticity, and axial circulation showed the most marked differences, even in patients without aortic dilation.
In patients with BAV with AAo dilation, eccentricity, backward velocity, velocity angle, regurgitation fraction, WSSC, axial vorticity, and axial circulation showed statistical differences in the AAo compared to volunteers, with higher values in patients with BAV. Moreover, the velocity, forward velocity, WSS, and WSSA also showed statistical differences between patients with BAV-AAoD and volunteers, but lower values characterized the v with BAV. The WSSC was the only present in the AAo of patients with BAV with AAo dilation and may thus arise from dilation per se (17, 37). Some of these parameters have been previously studied separately from a qualitative point of view (38, 39). Those studies have concluded that the altered jet direction (i.e., flow eccentricity) is one of the major contributors to aortic dilation in patients with BAV and related to axial circulation (refer to Figure 6). Flow eccentricity in patients with BAV (40) may alter the integrity of the aortic wall and promote dilation of the AAo (7, 41).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/885338/fcvm-09-885338-HTML/image_m/fcvm-09-885338-g006.jpg)
Figure 6. (A) A patient with maximum values of diameter and regurgitation fraction; (B) a patient with maximum WSSC and axial circulation; (C) a patient with positive axial circulation; and (D) one case with negative axial circulation.
Wall shear stress has also been extensively studied in patients with BAV (6, 9, 15–17, 21, 39, 42, 43). In this study, we found statistically significant differences for the magnitude of WSS between HV and patients with BAV with AAo dilated and all patients together. Interestingly, the WSS was not significant in patients with BAV-Non-AAoD. The WSSC showed larger values in BAV and significant differences in the AAo (Figure 3) in comparison with HV, even in non-dilated patients. Similar results have been reported in the literature (9, 44, 45).
To the best of our knowledge, the velocity angle and the axial circulation in 3D domain have not been previously studied in a large cohort of patients with BAV. The velocity angle represents the deviation of the flow velocity at each point in the vessel's lumen concerning their axial direction calculated with the Laplace equation (32) (refer to Supplementary Figure S4). The velocity angle also showed statistical differences between patients with BAV and HV. Other studies reported similar results when analyzing the direction of flow at the level of the aortic valve (1, 46–48). The circulation represents the integral of vorticity with respect to the cross-sectional area of the vessel, in our case, this vorticity was the axial vorticity at each point of the level set generated with the Laplace equation (32), and for this reason, we call it axial circulation (refer to Supplementary Figure S7 and the reference 36). This parameter showed statistical differences between BAV groups and HV, with higher values in the ascending aorta, bigger than 200 cm2/s for patients with BAV and <30 cm2/s for HV.
For phenotypes, BAV-RN showed a more significant difference in the proximal part of the aortic arch, probably influenced by the dilation of the vessel in this location (17). In the descending aorta, the BAV-RL showed more significant differences for parameters related to the flow turbulence such as vorticity, axial vorticity, axial circulation, helicity density, viscous dissipation, energy loss, and kinetic energy in the descending aorta with lower values for patients with BAV. Finally, the velocity angle and regurgitation fraction were the most significant for all regions.
Receiver operating characteristic curve analysis showed that in the AAo (Figure 4; Supplementary Table 4), eccentricity, backward velocity, velocity angle, regurgitation fraction, WSSC, axial vorticity, and axial circulation were the best performing parameters in the differentiation between volunteers and patients with BAV, with values under the curve bigger than 0.8. Interestingly, this value was equal to 0.96 for the velocity angle with a sensitivity and specificity of 90%, the largest among all parameters. We validated these results with the application of the MRMR classification algorithm, where the velocity angle, axial circulation, regurgitation fraction, diameter, eccentricity, circumferential WSS, backward velocity, forward velocity, and axial vorticity are the parameters with a bigger importance predictor score for most of the cases (refer to Supplementary Figures S11, S12; Supplementary Tables S8–S11). One interesting parameter analyzed in this study is the axial circulation, which can provide important information on the helical rotation of the blood flow in the entire vessel. Figure 6 shows one representative patient with BAV with positive (men 27 years, Figure 6C) and one with negative (men 76 years, Figure 6D) rotations of the blood flow in the ascending aorta.
Although many parameters were found to be different between HV and patients with BAV, only a few of them showed a good correlation with the diameter of the aorta. The diameter correlated with the forward velocity, WSSA, and vorticity in patients with AAo dilation. There was a strong inverse correlation between the diameter and the forward velocity for all BAV and patients with BAV-RL, and also, there was a strong correlation between the diameter and the velocity angle for the same groups. BAV-RL showed more parameters correlated with the diameter than all other cases.
Limitations
The main limitation of this study is its cross-sectional nature, which implies that no causal relationship could be inferred. Longitudinal studies are needed to confirm the present results. Another limitation relates to the limited sample size for the group of BAV-RootD, which was included in the non-dilated group, because statistical tests could not be performed on only seven patients. Furthermore, averaging parameters along the circumference of each region can induce sub-estimation of local value. Moreover, the movement of the aorta along the cardiac cycle was not considered in this study since technical limitations in 4D-flow CMR acquisitions (poor contrast and low signal-to-noise ratio) make it difficult to obtain a time-resolved segmentation of the aorta. For that reason, we only analyze the peak systolic cardiac phase. Another limitation of the proposed methodology is that the Laplace approach can be evaluated only for one vessel of interest without branches. Nonetheless, this limitation may also be present in centerline-based methods. The analysis of a large number of parameters can be overwhelming, and for that reason in future contributions, we may generate some artificial intelligence (AI) algorithms that can help to classify groups and categories of these parameters and thus be able to deliver only the parameters that are relevant for a particular pathology, but for this purpose, balanced groups will be needed.
Conclusion
Through a novel 3D quantification method, a limited number of hemodynamic parameters that differentiate between healthy volunteers and patients with BAV were unveiled and ranked. Those related to the direction of blood flow (forward velocity, velocity angle, regurgitation fraction, and WSSA) presented the best relationships with local diameter; however, eccentricity, backward velocity, velocity angle, regurgitation fraction, WSSC, axial vorticity, and axial circulation are the parameters that accurately differentiate between patients with BAV and volunteers. Longitudinal studies are required to assess whether these parameters are the significant contributors to aortic dilation in patients with BAV.
Data Availability Statement
The data analyzed in this study is subject to the following licenses/restrictions. The data sets generated during the current study will be made available with prior permission from the partner institutions, the toolbox used to process de data is free available on (https://github.com/JulioSoteloParraguez/4D-Flow-Matlab-Toolbox). Requests to access these datasets should be directed to JS, julio.sotelo@uv.cl.
Ethics Statement
The studies involving human participants were reviewed and approved by the Local Ethics Committee (Hospital Universitari Vall d'Hebron). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JS: responsible for the study concepts, design, and analysis. AG, LD-S, AR-M, AE, and JR-P: data acquisition and clinical interpretation of the data. JS, PF, HM, JM, DH, and SU: technical interpretation of the data. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by ANID – Millennium Science Initiative Program – ICN2021_004 and ANID – Millennium Science Initiative Program – NCN17_129, CONICYT-FONDECYT Postdoctorado #3170737, ANID – FONDECYT Postdoctorado #3220266, ANID Ph. D. Scholarship 21170592, ANID FONDECYT de Iniciación en Investigación #11200481, ANID FONDECYT #1181057, ANID Ph. D. Scholarship 21180391, the Spanish Society of Cardiology (SEC/FEC-INV-CLI 20/015) and the Biomedical Research Networking Center on Cardiovascular Diseases (CIBERCV). AG has received funding from the Spanish Ministry of Science, Innovation and Universities (IJC2018-037349-I).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank the Millennium Science Initiative of the Ministry of Economy, Development and Tourism, grant Nucleus for Cardiovascular Magnetic Resonance. This work was funded by the ANID – Millennium Science Initiative Program – ICN2021_004. We are also grateful to Biomedical Imaging Center at Pontificia Universidad Católica de Chile and Hospital Universitari Vall d'Hebron for support this research. Franco P. would like to thank CONICYT – PCHA/Doctorado-Nacional/2018-21180391.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.885338/full#supplementary-material
References
1. Kang JW, Song HG, Yang DH, Baek S, Kim DH, Song JM, et al. Association bsetween bicuspid aortic valve phenotype and patterns of valvular dysfunction and bicuspid aortopathy: comprehensive evaluation using MDCT and echocardiography. JACC Cardiovasc Imaging. (2013) 6:150–61. doi: 10.1016/j.jcmg.2012.11.007
2. Itagaki S, Chikwe JP, Chiang YP, Egorova NN, Adams DH. Long-term risk for aortic complications after aortic valve replacement in patients with bicuspid aortic valve vs. Marfan syndrome. J Am Coll Cardiol. (2015) 65:2363–9. doi: 10.1016/j.jacc.2015.03.575
3. Michelena HI, Prakash SK, Della Corte A, Bissell MM, Anavekar N, Mathieu P, et al. Bicuspid aortic valve: identifying knowledge gaps and rising to the challenge from the International Bicuspid Aortic Valve Consortium (BAVCon). Circulation. (2014) 129:2691–704. doi: 10.1161/CIRCULATIONAHA.113.007851
4. Evangelista A, Gallego P, Calvo-Iglesias F, Bermejo J, Robledo-Carmona J, Sánchez V, et al. Anatomical and clinical predictors of valve dysfunction and aortic dilation in bicuspid aortic valve disease. Heart. (2018) 104:566–73. doi: 10.1136/heartjnl-2017-311560
5. Schaefer BM, Lewin MB, Stout KK, Gill E, Prueitt A, Byers PH, et al. The bicuspid aortic valve: an integrated phenotypic classification of leaflet morphology and aortic root shape. Heart. (2008) 94:1634–8. doi: 10.1136/hrt.2007.132092
6. Mahadevia R, Barker AJ, Schnell S, Entezari P, Kansal P, Fedak PW, et al. Bicuspid aortic cusp fusion morphology alters aortic three-dimensional outflow patterns, wall shear stress, and expression of aortopathy. Circulation. (2014) 129:673–82. doi: 10.1161/CIRCULATIONAHA.113.003026
7. Guala A, Rodriguez-Palomares J, Dux-Santoy L, Teixido-Tura G, Maldonado G, Galian L, et al. Influence of aortic dilation on the regional aortic stiffness of bicuspid aortic valve assessed by 4-Dimensional flow cardiac magnetic resonance: comparison with Marfan syndrome and degenerative aortic aneurysm. JACC Cardiovasc Imaging. (2019) 12:1020–9. doi: 10.1016/j.jcmg.2018.03.017
8. Galian-Gay L, Carro Hevia A, Teixido-Turà G, Rodríguez Palomares J, Gutiérrez-Moreno L, Maldonado G, et al. Familial clustering of bicuspid aortic valve and its relationship with aortic dilation in first-degree relatives. Heart. (2019) 105:603–8. doi: 10.1136/heartjnl-2018-313802
9. Guzzardi DG, Barker AJ, van Ooij P, Malaisrie SC, Puthumana JJ, Belke DD, et al. Valve-Related Hemodynamics Mediate Human Bicuspid Aortopathy: Insights From Wall Shear Stress Mapping. J Am Coll Cardiol. (2015) 66:892–900. doi: 10.1016/j.jacc.2015.06.1310
10. Borger MA, Fedak PWM, Stephens EH, Gleason TG, Girdauskas E, Ikonomidis JS, et al. The American Association for Thoracic Surgery consensus guidelines on bicuspid aortic valve-related aortopathy: Full online-only version. J Thorac Cardiovasc Surg. (2018) 156:e41–74. doi: 10.1016/j.jtcvs.2017.10.161
11. Pape LA, Tsai TT, Isselbacher EM, Oh JK, O'gara PT, Evangelista A, et al. Aortic diameter >or = 5.5 cm is not a good predictor of type A aortic dissection: observations from the International Registry of Acute Aortic Dissection (IRAD). Circulation. (2007) 116:1120–7. doi: 10.1161/CIRCULATIONAHA.107.702720
12. Cavalcante JL, Lalude OO, Schoenhagen P, Lerakis S. Cardiovascular magnetic resonance imaging for structural and valvular heart disease interventions. JACC Cardiovasc Interv. (2016) 9:399–425. doi: 10.1016/j.jcin.2015.11.031
13. Whitlock MC, Hundley WG. Noninvasive imaging of flow and vascular function in disease of the aorta. JACC Cardiovasc Imaging. (2015) 8:1094–106. doi: 10.1016/j.jcmg.2015.08.001
14. Dyverfeldt P, Bissell M, Barker AJ, Bolger AF, Carlhäll CJ, Ebbers T, et al. 4D flow cardiovascular magnetic resonance consensus statement. J Cardiovasc Magn Reson. (2015) 17:72. doi: 10.1186/s12968-015-0174-5
15. Bollache E, Guzzardi DG, Sattari S, Olsen KE, Di Martino ES, Malaisrie SC, et al. Aortic valve-mediated wall shear stress is heterogeneous and predicts regional aortic elastic fiber thinning in bicuspid aortic valve-associated aortopathy. J Thorac Cardiovasc Surg. (2018) 156:2112-−20.e2. doi: 10.1016/j.jtcvs.2018.05.095
16. 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 dilation in bicuspid aortic valve disease. J Cardiovasc Magn Reson. (2018) 20:28. doi: 10.1186/s12968-018-0451-1
17. Dux-Santoy L, Guala A, Teixidó-Turà G, Ruiz-Muñoz A, Maldonado G, Villalva N, et al. Increased rotational flow in the proximal aortic arch is associated with its dilation in bicuspid aortic valve disease. Eur Heart J Cardiovasc Imaging. (2019) 20:1407–17. doi: 10.1093/ehjci/jez046
18. Guala A, Rodriguez-Palomares J, Galian-Gay L, Teixido-Tura G, Johnson KM, Wieben O, et al. Partial aortic valve leaflet fusion is related to deleterious alteration of proximal aorta hemodynamics. Circulation. (2019) 139:2707–9. doi: 10.1161/CIRCULATIONAHA.119.039693
19. Guala A, Evangelista A, Teixido-Tura G, La Mura L, Dux-Santoy L, Ruiz-Muñoz et al. Leaflet fusion length is associated with aortic dilation and flow alterations in non-dysfunctional bicuspid aortic valve. Eur Radiol. (2021) 31:9262-9272. doi: 10.1007/s00330-021-08016-3
20. Burris NS, Hope MD. 4D flow MRI applications for aortic disease. Magn Reson Imaging Clin N Am. (2015) 23:15–23. doi: 10.1016/j.mric.2014.08.006
21. Rahman O, Scott M, Bollache E, Suwa K, Collins J, Carr J, et al. Interval changes in aortic peak velocity and wall shear stress in patients with bicuspid aortic valve disease. Int J Cardiovasc Imaging. (2019) 35:1925–34. doi: 10.1007/s10554-019-01632-7
22. Garcia J, Barker AJ, Collins JD, Carr JC, Markl M. Volumetric quantification of absolute local normalized helicity in patients with bicuspid aortic valve and aortic dilation. Magn Reson Med. (2017) 78:689–701. doi: 10.1002/mrm.26387
23. 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. doi: 10.1002/mrm.24962
24. David Marlevi, Julio Sotelo, Ross Grogan-Kaylor, et al. False lumen pressure estimation in type B aortic dissection using 4D flow cardiovascular magnetic resonance: comparisons with aortic growth. J Cardiovasc Magn Reason. (2021) 23:51. doi: 10.1186/s12968-021-00741-4
25. Burris NS, Nordsletten DA, Sotelo JA, et al. False lumen ejection fraction predicts growth in type B aortic dissection: preliminary results. Eur J Cardio-Thorac Surgs. (2020) 57:896–903. doi: 10.1093/ejcts/ezz343
26. Sotelo J, Valverde I, Martins D, et al. Impact of aortic arch curvature in flow hemodynamics in patients with transposition of the great arteries after arterial switch operation. Eur Heart J Cardiovasc Imaging. (2022) 23:402–11. doi: 10.1093/ehjci/jeaa416
27. Warmerdam E, van Assen HC, Sotelo J, Grotenhuis HB. Abnormal vortex formation in the right pulmonary artery after the arterial switch operation. Eur Heart J Case Rep. (2021) 5:ytaa549. doi: 10.1093/ehjcr/ytaa549
28. Gu T, Korosec FR, Block WF, Fain SB, Turk Q, Lum D, et al. PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. AJNR Am J Neuroradiol. (2005) 26:743–9.
29. Bieging ET, Frydrychowicz A, Wentland A, Landgraf BR, Johnson KM, Wieben O, et al. In vivo three-dimensional MR wall shear stress estimation in ascending aortic dilatation. J Magn Reson Imaging. (2011) 33:589–97. doi: 10.1002/jmri.22485
30. Campens L, Demulier L, De Groote K, Vandekerckhove K, De Wolf D, Roman MJ, et al. Reference values for echocardiographic assessment of the diameter of the aortic root and ascending aorta spanning all age categories. Am J Cardiol. (2014) 114:914–20. doi: 10.1016/j.amjcard.2014.06.024
31. Della Corte A, Bancone C, Dialetto G, Covino FE, Manduca S, Montibello MV, et al. The ascending aorta with bicuspid aortic valve: a phenotypic classification with potential prognostic significance. Eur J Cardiothorac Surg. (2014) 46:240–7. doi: 10.1093/ejcts/ezt621
32. Sotelo J, Dux-Santoy L, Guala A, Rodríguez-Palomares J, Evangelista A, Sing-Long C, et al. 3D axial and circumferential wall shear stress from 4D flow MRI data using a finite element method and a laplacian approach. Magn Reson Med. (2018) 79:2816–23. doi: 10.1002/mrm.26927
33. Sotelo J, Urbina J, Valverde I, Mura J, Tejos C, Irarrazaval P, et al. Three-dimensional quantification of vorticity and helicity from 3D cine PC-MRI using finite-element interpolations. Magn Reson Med. (2018) 79:541–53. doi: 10.1002/mrm.26687
34. Sotelo J, Urbina J, Valverde I, Tejos C, Irarrazaval P, Andia ME, et al. 3D quantification of wall shear stress and oscillatory shear index using a finite-element method in 3D CINE PC-MRI data of the thoracic aorta. IEEE Trans Med Imaging. (2016) 35:1475-87. doi: 10.1109/TMI.2016.2517406
35. Sotelo J, Bächler P, Urbina J, Crelier G, Toro L, Ferreiro M, et al. Quantification of pulmonary regurgitation in patients with repaired Tetralogy of Fallot by 2D phase-contrast MRI: Differences between the standard method of velocity averaging and a pixel-wise analysis. JRSM Cardiovasc Dis. (2017) 6:2048004017731986 doi: 10.1177/2048004017731986
36. Sotelo J, Bissell MM, Jiang Y, Mella H, Mura J, Uribe S. Three-dimensional quantification of circulation using finite-element methods in four-dimensional flow MR data of the thoracic aorta. Magn Reson Med. (2022) 87:1036–45. doi: 10.1002/mrm.29004
37. Bürk J, Blanke P, Stankovic Z, Barker A, Russe M, Geiger J, et al. Evaluation of 3D blood flow patterns and wall shear stress in the normal and dilated thoracic aorta using flow-sensitive 4D CMR. J Cardiovasc Magn Reson. (2012) 14:84. doi: 10.1186/1532-429X-14-84
38. Schnell S, Smith DA, Barker AJ, Entezari P, Honarmand AR, Carr ML, et al. Altered aortic shape in bicuspid aortic valve relatives influences blood flow patterns. Eur Heart J Cardiovasc Imaging. (2016) 17:1239–47. doi: 10.1093/ehjci/jew149
39. Bissell MM, Hess AT, Biasiolli L, Glaze SJ, Loudon M, Pitcher A, et al. Aortic dilation in bicuspid aortic valve disease: flow pattern is a major contributor and differs with valve fusion type. Circ Cardiovasc Imaging. (2013) 6:499–507. doi: 10.1161/CIRCIMAGING.113.000528
40. Sigovan M, Hope MD, Dyverfeldt P, Saloner D. Comparison of four-dimensional flow parameters for quantification of flow eccentricity in the ascending aorta. J Magn Reson Imaging. (2011) 34:1226–30. doi: 10.1002/jmri.22800
41. Burris NS, Dyverfeldt P, Hope MD. Ascending Aortic Stiffness with Bicuspid Aortic Valve is Variable and Not Predicted by Conventional Parameters in Young Patients. J Heart Valve Dis. (2016) 25:270–80. doi: 10.1186/1532-429X-17-S1-Q80
42. Dux-Santoy L, Guala A, Sotelo J, Uribe S, Teixidó-Turà G, Ruiz-Muñoz A, et al. Low and oscillatory wall shear stress is not related to aortic dilation in patients with bicuspid aortic valve: a time-resolved 3-dimensional phase-contrast magnetic resonance imaging study. Arterioscler Thromb Vasc Biol. (2020) 40:e10–20. doi: 10.1161/ATVBAHA.119.313636
43. van Ooij P, Garcia J, Potters WV, Malaisrie SC, Collins JD, Carr JC, et al. Age-related changes in aortic 3D blood flow velocities and wall shear stress: implications for the identification of altered hemodynamics in patients with aortic valve disease. J Magn Reson Imaging. (2016) 43:1239–49. doi: 10.1002/jmri.25081
44. Uretsky S, Gillam LD. Nature versus nurture in bicuspid aortic valve aortopathy: more evidence that altered hemodynamics may play a role. Circulation. (2014) 129:622–4. doi: 10.1161/CIRCULATIONAHA.113.007282
45. van Ooij P, Markl M, Collins JD, Carr JC, Rigsby C, Bonow RO, et al. Aortic valve stenosis alters expression of regional aortic wall shear stress: new insights from a 4-dimensional flow magnetic resonance imaging study of 571 subjects. J Am Heart Assoc. (2017) 6:e005959. doi: 10.1161/JAHA.117.005959
46. 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. doi: 10.1038/srep32316
47. den Reijer PM, Sallee D, van der Velden P, Zaaijer ER, Parks WJ, Ramamurthy S, et al. Hemodynamic predictors of aortic dilatation in bicuspid aortic valve by velocity-encoded cardiovascular magnetic resonance. J Cardiovasc Magn Reson. (2010) 12:4. doi: 10.1186/1532-429X-12-4
Keywords: 4D flow CMR, finite elements, hemodynamics parameters, Bicuspid aortic valve, congenital heart disease, aneurysm, magnetic resonance imaging (MRI), vascular disease
Citation: Sotelo J, Franco P, Guala A, Dux-Santoy L, Ruiz-Muñoz A, Evangelista A, Mella H, Mura J, Hurtado DE, Rodríguez-Palomares JF and Uribe S (2022) Fully Three-Dimensional Hemodynamic Characterization of Altered Blood Flow in Bicuspid Aortic Valve Patients With Respect to Aortic Dilatation: A Finite Element Approach. Front. Cardiovasc. Med. 9:885338. doi: 10.3389/fcvm.2022.885338
Received: 28 February 2022; Accepted: 22 April 2022;
Published: 18 May 2022.
Edited by:
Chengcheng Zhu, University of Washington, United StatesCopyright © 2022 Sotelo, Franco, Guala, Dux-Santoy, Ruiz-Muñoz, Evangelista, Mella, Mura, Hurtado, Rodríguez-Palomares and Uribe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sergio Uribe, suribe@uc.cl