ORIGINAL RESEARCH article
High Wall Shear Stress can Predict Wall Degradation in Ascending Aortic Aneurysms: An Integrated Biomechanics Study
- 1Department of Surgery and Cancer, Imperial College London, London, United Kingdom
- 2Department of Chemical Engineering, Imperial College London, London, United Kingdom
- 3Department of Bioengineering, Imperial College London, London, United Kingdom
- 4Politecnico di Milano, Milan, Italy
- 5MRC London Institute of Medical Sciences, Imperial College London, London, United Kingdom
- 6Barts Heart Centre, London, United Kingdom
Background: Blood flow patterns can alter material properties of ascending thoracic aortic aneurysms (ATAA) via vascular wall remodeling. This study examines the relationship between wall shear stress (WSS) obtained from image-based computational modelling with tissue-derived mechanical and microstructural properties of the ATAA wall using segmental analysis.
Methods: Ten patients undergoing surgery for ATAA were recruited. Exclusions: bicuspid aortopathy, connective tissue disease. All patients had pre-operative 4-dimensional flow magnetic resonance imaging (4D-MRI), allowing for patient-specific computational fluid dynamics (CFD) analysis and anatomically precise WSS mapping of ATAA regions (6–12 segments per patient). ATAA samples were obtained from surgery and subjected to region-specific tensile and peel testing (matched to WSS segments). Computational pathology was used to characterize elastin/collagen abundance and smooth muscle cell (SMC) count.
Results: Elevated values of WSS were predictive of: reduced wall thickness [coef −0.0489, 95% CI (−0.0905, −0.00727), p = 0.022] and dissection energy function (longitudinal) [−15,0, 95% CI (−33.00, −2.98), p = 0.048]. High WSS values also predicted higher ultimate tensile strength [coef 0.136, 95% CI (0 0.001, 0.270), p = 0.048]. Additionally, elevated WSS also predicted a reduction in elastin levels [coef −0.276, 95% (CI −0.531, −0.020), p = 0.035] and lower SMC count ([oef −6.19, 95% CI (−11.41, −0.98), p = 0.021]. WSS was found to have no effect on collagen abundance or circumferential mechanical properties.
Conclusions: Our study suggests an association between elevated WSS values and aortic wall degradation in ATAA disease. Further studies might help identify threshold values to predict acute aortic events.
Ascending thoracic aortic aneurysm (ATAA) is a permanent and irreversible dilatation of the thoracic aorta. Many patients remain asymptomatic until acute presentation with rupture or dissection (AD), which is associated with a 50% early mortality rate (Howard et al., 2013). Death from aortic aneurysm-related emergencies occurs at a rate of 2.4/100,000, and is one of the most common causes of death amongst conditions requiring emergency surgery in high-income countries (Stewart et al., 2014). Survivors of acute events often need repeat intervention and have high rates of stroke and renal failure. This has a huge impact on their quality-of-life, and carries significant societal burden (Tsai et al., 2006).
The evolution in our understanding of the disease from large observational studies has led to the reliance on ATAA diameter as the primary predictor of future AD. Intervention is recommended in patients who have a maximal aortic diameter ≥55 mm, with lower thresholds (≥50 mm) for those with bicuspid aortic valves (BAV) or connective tissue disease (Hiratzka et al., 2010; Erbel et al., 2014). Interventions in ATAA can be complex and not without risk, involving open surgery (particularly for ascending and arch aneurysms), endovascular stenting, or a combination of both. This drives the need for an accurate prognostic prediction prior to subjecting patients to treatment that involves moderate risk. However, isolated ATAA diameter measurements are inadequate predictors of acute aortic events. International registry data (>4,400 dissection patients) has highlighted that 40% of dissections occur below 50 mm (and up to 60% below 55 mm) (Pape et al., 2007). In this light, clinicians at present are unable to accurately predict the prognosis of enlarged thoracic aortas from routine imaging. Indeed, diameter alone fails to take into account local flow patterns, aortic wall stresses, mechanical properties and wall thickness. A patient specific approach is required.
Despite having the potential to predict disease progression, biomechanical assessment of the aorta remains experimental with limited translation into clinical practice. Biomechanically, AD occurs when haemodynamic forces exceed the aortic wall strength leading to intimal tear and false lumen propagation. It is possible that long-term exposure to certain blood flow patterns lead to changes in aortic wall structure (i.e., remodeling) that predispose wall degradation and a higher risk of AD. The identification of which flow patterns correlate with changes in aortic wall structure might lead to a more rational and useful patient-specific predictor (Geiger et al., 2012; Condemi et al., 2017; Piatti et al., 2017; Bollache et al., 2018).
The degradation of extracellular matrix (ECM) structures in ATAA is well reported. Methods to quantify degeneration center around computational histology techniques to characterize microstructural elements, such as elastin and collagen. Studies employing these methods have linked abnormal WSS patterns to disruption of the aortic microarchitecture (Guzzardi et al., 2015; Bollache et al., 2018). However, studies that use flow-to-tissue spatial registration (that incorporate the whole aneurysm) have been limited.
This study aims to explore the relationship between haemodynamic parameters, material properties and composition of the aortic wall in ATAA using whole aneurysm samples, with a view to assess the predictive ability of flow on aortic wall degeneration. We hypothesize elevated WSS in ATAAs is associated with accelerated degenerative disease.
Materials and Methods
The study was ethically approved (17/NI/0160) by the Health Research Authority (HRA) in the United Kingdom and was sponsored by the Imperial College London Joint Research and Compliance Office, as defined under the sponsorship requirements of the Research Governance Framework (2005). The study was designed as a cohort study during the years 2018-2020, incorporating computational flow analysis, in-vitro vascular wall mechanical testing and microstructural quantification. Patients and public were involved in the design, conduct, reporting, and dissemination plans of our research via the London Aortic Mechanobiology Working Group.
A total of 10 patients undergoing proximal aortic surgery (either aortic root, ascending aorta, proximal arch replacements, or a combination of these) were recruited into this cohort study. The main exclusion criteria were: connective tissue disease (i.e. Marfans, Ehler-Danlos, Loeys-Dietz), bicuspid aortic valve and redo-aortic surgery. Patient characteristics were ascertained from clinic letters and pre-operative echocardiography. Recruited patients provided informed consent for participation in the study, which involved the additional MR imaging and provision of tissue sample from surgery.
Patients were scheduled to undergo cardiac gated magnetic resonance imaging (MRI) scanning at a time point prior to surgery (range 1–90 days) using a 3T MRI scanner (Siemens Healthcare, Erlangen, Germany). The addition of a time-resolved 3-dimensional sequence (4D-flow) visualized and measured temporal changes and flow-patterns throughout the whole volume of the thoracic aorta.
Aneurysm Tissue Characterization
Patient-matched aneurysm specimens were obtained en-bloc and acquired immediately after surgical excision in the operating theatre. Rectangular or dog-bone shaped subsections aligned in either the longitudinal or circumferential direction were punched out at specific locations (Figure 1A) and assessed for tensile mechanical properties, peel strength and strain inhomogeneities (Olchanyi et al., 2020). In particular, for each patient, the aneurysm region was unfolded at the outer wall and subdivided into up to six (depending on the size of the collected sample) circumferential anatomical regions to the right (R1, R2, R3) and left (L1, L2, L3) of the midline (Figure 1B). These were further divided into upper and lower regions, giving 6 to 12 segments in total per patient (Figure 1B). These segments were matched to the 6–12 segmental regions undergoing computational WSS mapping (Section 2.4; Figure 1C), using the anterior midline (widest region of the aneurysm region in the vertical plane) as the common frame of reference. Whole-aneurysm tissue specimens underwent thickness mapping (∼0.01 µm resolution using a bench top device - Litematic VL-50-B Mitutoyo Ltd.).
FIGURE 1. Matching flow to material properties. (A) Reference map used for collection of samples throughout the whole aortic specimen. Segments are classified using alphanumeric codes in which the first letter and number denote the position (R = right, L = left), the second letter denotes the sample orientation (C = circumferential, L = longitudinal), the third letter denotes the sample type (S = dogbone, P = rectangular) (B–C) WSS mapping on to aneurysmal wall surface, which is unraveled by division along the anterior. This is subdivided into six circumferential regions (R1-R3, L1-L3), with a horizontal dividing line (12 segments total). (D) Streamlines demonstrating haemodynamic flow derived from CFD within an anatomically (patient-specific) segmented aortic wall.
Sharp stencils were used to create 20 mm × 5 mm dogbone samples (for uniaxial testing). The dogbone shaped tissue provided larger shoulders for easier gripping, whereas the narrower gauge section in the center provided a smaller cross-sectional area for a focused region of deformation and failure. The samples for peel testing required a different shape: the width of the peeling sample needed to be kept constant throughout the peeling arc. As such rectangular stencils were used (5 mm × 20 mm) to obtain eight further samples from both halves of the aneurysm specimen (both circumferential and longitudinal) –these would be used for peel testing.
Uniaxial tensile tests were performed on all dog-bone shaped subsections using a Test Resources R-Series Controller with a 44 N load cell. All tests were conducted in an environmental chamber containing phosphate buffered saline (PBS), maintained at 37°C. Tensile force-elongation data were used to provide estimations of region specific ultimate tensile strength (UTS) and maximum tangential stiffness (MTS) as described by previous groups (Vorp et al., 2003). Peel testing allowed estimation of peel force (Fpeel) and dissection energy function (DEF) for rectangular subsections.
Computational Fluid Dynamics (CFD)
Image-based CFD modelling was used to evaluate pre-operative aneurysmal WSS distribution in a patient-specific and region-specific fashion. To this aim, the 3D thoracic aortic geometry, including the arch branches, of each patient was reconstructed from MRI data (bright blood). Reconstructed geometries (Figures 1C, D) started upstream of the aneurysm, in the aortic root between mid-sinus and STJ levels, and ended distally, in the descending thoracic aorta at the level of the pulmonary artery.
Patient-specific geometries were discretized using unstructured meshes (Supplementary Figure S1) with a tetrahedral core and 10 prismatic layers at the walls. Local mesh refinements were prescribed at the aneurysm wall and arch branches. Mesh refinement in the region of interest (i.e. aneurysm wall) was guided by flow features observed from 4D flow MRI data (i.e. a finer mesh was designed where higher velocity gradients were observed). Sensitivity analyses were conducted to ensure mesh-independent results. Final meshes consisted of 5.5–16.6 million elements (Pirola et al., 2019).
For each patient, 3D subject-specific velocity profiles of the aortic valve were incorporated at the computational model inlet from 4D flow MR images (Armour et al., 2021). This therefore took into account the patient-to-patient variation in the ascending aorta hemodynamics which arise due to differences in left ventricular outflow tract and aortic valve anatomy (Figure 1D). Model outlets, located distal to the aneurysm, consisted of the descending aorta outlet and the arch branches. 3-element Windkessel model was applied at all outlets; model parameters were tuned using patient-specific central mean aortic pressure and blood flow rates evaluated using 4D flow MRI data (Pirola et al., 2017; 2019). Briefly, the total resistance (RT) of the 3-element Windkessel model was calculated as RT = Pm/Qm (Les et al., 2010). Pm is the mean central pressure of the patient, which was measured ∼30 min before the MR scan using a BP Plus device (BP Plus, Uscom, Australia). Qm is the mean flow through the outlet. This was evaluated using 4D flow MR images. The proximal resistance was evaluated as R1 = ρc/A (Xiao et al., 2014, where ρ is the blood density, c is the pulse wave speed and A is the outlet cross-sectional area. The pulse wave speed was evaluated as c = 13.3/(2r)0.3 Reymond et al. (2009)), where r is the outlet radius. The distal resistance was evaluated as R2 = RT-R1 (LaDisa et al., 2011). Total vasculature compliance was calculated as C = τ/RT, where τ is the time-constant of the exponential diastolic pressure-fall (Xiao et al. (2014)). This was assumed equal to 1.79 and 1.92 s for normotensive and hypertensive subjects, respectively (Simon et al., 1979).
Simulations were run in ANSYS CFX (v 15.0). The aortic wall was assumed to be rigid with a no-slip condition. Blood was modelled as an incompressible Newtonian fluid, with a density of 1,060 kg m−3 and viscosity of 4·10−3 Pa s. High-resolution advection scheme and second-order Backward Euler scheme were employed for spatial and temporal discretization of the Navier-Stokes equations, respectively. A fixed timestep of 10−3 s was used, and the maximum RMS residual was set to 10−5 as a convergence criterion. The shear stress transport transitional (SST-Tran) model was used to account for the transitional nature of the aortic flow (Tan et al., 2009), with a 1% turbulence level prescribed as the inlet boundary condition. The patient-specific central diastolic pressure was used as the initial condition. Simulations were run for the number of cycles necessary to achieve periodicity. The period of the cardiac cycle of each patient was recorded during 4D flow MR acquisition and stored in the image header. Periodicity was considered as achieved when differences in pulse pressure and pressure maxima between two consecutive heartbeats were less than 3 and 1%, respectively, at each model outlet. Results from the last cycle only were analyzed. Before proceeding with result analysis, for each patient, predicted blood flow features were checked against 4D flow MR acquired features.
Computational result postprocessing workflow was designed to coincide with tissue processing for mechanical and histological testing: for each patient, the aneurysm region was unfolded at the outer wall and subdivided into segments corresponding to the aortic aneurysm tissue segments used to characterize the material properties (Figures 1B,C). Results were analyzed in ANSYS EnSight.
The primary derived measurement from CFD analysis was the magnitude of the wall shear stress (WSS) vector. For each anatomical region, several WSS-derived parameters were evaluated. These are summarized in Table 1. Briefly, WSS values are reported as maxima in time (WSSmax, WSSmean) and time-averaged (TAWSSmax, TAWSSmean, WSSTmeanMax). Mean and max subscripts refer to spatial mean and maxima, respectively. Tmean superscript refers to temporal mean. The spatial mean value was calculated over each segment (R1, R2, R3, L1, L2, L3, upper and lower).
TABLE 1. Calculated wall shear stress (WSS) parameters from computational fluid dynamics (CFD) analysis.
Full circumferential rings of tissue were obtained from the inferior border of the ATAA specimen (Figure 2A). Whole slide imaging was performed using a high-resolution digital optical system (Hamamatsu TM) and uploaded onto the digital processing software QuPath (https://qupath.github.io). Cross-sectional images of each ATAA ring were divided into six circumferential regions (Figure 2B), matching the right-to-left segments used in the CFD and mechanical analyses (Figures 1B, C). Using a pre-defined workflow (Figure 2C), microstructural density calculation of medial structural proteins was conducted, namely of elastin (from slides stained with Elastin Van Gieson) and collagen (slides stained with picrosirius red). Using the H&E stained slides (Figure 2D), the thresholding function was used to highlight the stained cells of the medial layer. Setting the range of particle size from 20 to 100 µm the total cell count was extracted and divided by the total area to obtain the density.
FIGURE 2. Computational pathology: regional histological analysis. (A) the origin of the histological specimen from the ATAA region. (B) Whole slide imaging demonstrating the division into six segmental regions. (C) Computational steps used to quantify elastin density (using EVG stained slide image), including RGB-stack image conversion, thresholding and measurement. Collagen quantification uses the same steps but for a PSR stained image. (D) displays the computational steps used to measure smooth muscle cell count (from an H&E-stained slide image) including particle analysis, using a size window of 20–100 μm.
Amongst the 10 recruited patients, a total of 371 ATAA anatomical samples had matched CFD, material properties and microstructural features, and were thus included in the analysis, being treated as separate data points. Of these samples, thickness was measured in 102 samples, 63 samples (29 longitudinal, 34 circumferential) were used for tensile testing, 63 samples (30 longitudinal, 33 circumferential) for peel testing, 143 for histological analyses (60, 35 and 48 for elastin, collagen and SMCs, respectively). All haemodynamic parameters were acquired as continuous variables, as were the histological and mechanical data. Where relevant, data were reported as means and standard deviations, being averaged across the whole patient cohort. All statistical analysis were conducted using STATA 13.0 (Stata Corp. College Station, TX, United States).
Univariable linear regression analysis was used to explore the effect of WSS on acquired material parameters (including mechanical and histological). These models tested the hypothesis of a flow-mediated degenerative process in the aortic wall giving rise to altered material properties. Results of regression analysis were reported as standardized beta-coefficients with 95% confidence intervals. The significance level was set at α = 0.05. Values of longitudinal Fpeel, longitudinal DEF, UTS and MTS were positively skewed; values of elastin abundance were negatively skewed. Therefore, logarithmic transformation with skewness correction was applied in STATA to these datasets.
Scatter plots and trendlines have been generated using Microsoft Excel, using log10 transformation for longitudinal DEF and UTS, and log10 (90-x) transformation for elastin abundance.
Multilevel Mixed Effects (Hierarchical) Linear Models
As the aortic segments arose from among 10 different subjects, multilevel mixed-effect linear regression models were further constructed to account for the hierarchical structure of the data points. The data arising from each aortic segment was nested into a clustered hierarchical structure using two levels in the model: 1) patient; and 2) orientation of the subsection (circumferential versus longitudinal) (Figure 3). Univariate regression was again carried out, assessing the effect of WSS separately on wall thickness, UTS, MTS and histological parameters. The significance level for all models was set at α = 0.05.
Of the cohort of 10 patients, the majority were male (8/10), and of Caucasian ethnicity (8/10). Patient age was (mean ± standard deviation) 63.9years ±6.6 (Table 2). Whilst none had syndromic disease, 3 reported a family history of aortic disease. The mean aneurysm diameter was 54.7 mm ± 7.5. Three patients had root aneurysms, two patients had arch aneurysms, with the remaining patients having isolated ascending aorta aneurysms. All patients had good ventricular function (mean left ventricular ejection fraction 57.2% ± 9.0) with 3 patients suffering from severe aortic regurgitation. Maxima of jet velocity and area-averaged velocity at the model inlet were 1.86 ± 0.93 m/s and 0.51 ± 0.13 m/s (mean ± standard deviation), respectively.
The patient-specific TAWSS distribution maps for all 10 patients are shown in Figure 4. The R1 segment (outer right wall) was the region yielding the highest values of WSS (Table 3), including temporal maximum WSS values (WSSmax and WSSmean) as well as time-averaged values (TAWSSmax and TAWSSmean). The highest value for WSSmax (24.98 ± 7.79 Pa) in the R1 region contrasted to the lowest value yielded from the inner curve (10.18 ± 4.14 Pa). When measuring the (peak in time) WSS averaged over space (WSSmean), the R1 region still yielded the highest value (11.68 ± 6.42 Pa). For TAWSSmax, R1, also the region of highest WSS values, was 4.85Pa ± 2.07, compared to the inner curve (2.67 ± 0.81 Pa). This indicated that as well as the asymmetric peak flow patterns reached in systole acting on the potentially fragile intima in ATAA, the sustained stress over the cardiac cycle, as displayed by the TAWSS, was also asymmetrical and affecting the outer curve more intensely.
FIGURE 4. Patient-specific TAWSS distribution maps in the ten patients with ATAA. Regions in red indicate TAWSS >3 Pa, which, on gross inspection, are localized to the aneurysmal aorta, and particularly to the outer curve.
TABLE 3. Mean wall shear stress (WSS) parameters (±standard deviation) in anatomical regions around the aneurysm (averaged over 10 patients). All values are in Pascals (Pa).
Comparing WSS With Aortic Material Properties
In order to test the hypothesis of flow mediated wall degeneration, WSS was compared with material properties of the aortic wall. Statistically significant linear regression analysis results are summarized in Table 4. Scatter plots and trendlines showing the relationship between WSS and aortic wall material properties are shown in Figure 5.
TABLE 4. Linear regression analysis: comparison of measurements of WSS per ATAA segment and tissue-derived parameters of corresponding segment from patient-specific excised tissue. circ = circumferential, long = longitudinal, SMC = smooth muscle cell. All results are statistically significant (p-value < 0.05).
FIGURE 5. Scatter plots and trendlines demonstrating the relationship between WSS parameters (x-axes) and aortic wall material properties. (A) aortic wall thickness (mm) vs TAWSSmax (Pa) and (B) WSSTmeanMax (Pa); (C) longitudinal dissection energy function (log10 transformed, Log DEFlong) vs TAWSSmean (Pa); (D) ultimate tensile strength (log10 transformed, Log UTS) vs TAWSSmax (Pa); (E) Elastin abundance (log10 (90-x) transformed) vs TAWSSmax (Pa); (F) Smooth muscle cell (SMC) count (\mm2) vs TAWSSmax (Pa) and (G) WSSTmeanMax (Pa).
The data on measured wall thickness was normally distributed. Linear regression analysis was conducted which showed a statistically significant influence (p < 0.05) of TAWSSmax and WSSTmeanMax parameters on wall thickness (Table 4). These yielded negative coefficients of variance, indicating that higher WSS values are associated with aortic wall thinning (Table 4; Figures 5A,B).
Peeling force (Fpeel) and dissection energy function (DEF) (as measured from tissue mechanical testing) provide a surrogate for the likelihood of dissection (i.e. separation of the medial layers) per region of interest. These were obtained from ATAA specimens in both the circumferential and longitudinal directions. Regression analysis (Table 4) comparing log transformed DEF values in the longitudinal orientation (DEFlong) with TAWSSMean yielded a statistically significant relationship [coef −0.211, 95% CI (−0.427, −0.062), p = 0.048]. This relationship displayed a negative coefficient result: higher WSS values were associated with lower DEFlong values (as can also be observed in Figure 5C). In contrast, in the circumferential direction Fpeel and DEF values showed an opposite trend, with higher values corresponding to higher WSS. However, results for Fpeel (both longitudinal and circumferential) and circumferential DEF were not statistically significant.
Tensile Properties: Ultimate Tensile Strength (UTS) and Maximum Tangential Stiffness (MTS)
The comparison of log-transformed UTS data to TAWSSmax revealed a statistically significant positive correlation (coef 0.136, 95% CI [0 0.001, 0.270], p = 0.048) (Table 4; Figure 5D). A comparison of log-transformed MTS data to WSS revealed a nonsignificant trend toward higher values of WSS leading to stiffer aortic tissue (coef 0.120, 95% CI [-0.009, 0.248], p = 0.068). Comparisons to other WSS measures failed to show statistical significance. These results together suggest that high WSS leads to stiffer aortic tissue with higher tensile strength, but the statistical significance is marginal.
Elastin abundance was found to be lower in areas with higher TAWSSMax [coef −0.276, 95% CI (−0.531, −0.020), p = 0.035] (Table 4; Figure 5E). An opposing, but statistically insignificant, trend was found for collagen abundance, i.e., higher TAWSS showed some association with areas of high collagen content. Higher WSS levels were also associated with lower counts of SMCs [TAWSSMax: coef −6.19, 95% CI (−11.41, −0.98), p = 0.021; WSSTmeanMax: coef −5.87, 95% CI (−10.12, −1.62), p = 0.008] (Table 4 and Figures 5F,G).
Multi-Level Mixed Effects Regression
TABLE 5. Results of multilevel mixed-effects linear model for the main outcome of tissue thickness. The fixed effects part of the model tested the influence of TAWSSmax on thickness, whilst the random effect part of the model tested the influence of the patient. From these results, the influence of WSS alone on thickness does not occur in isolation and the variance occurring at the patient level is an important influencing factor.
TABLE 6. Multilevel mixed effects linear regression: association between wall shear stress (WSSTmeanMax) and smooth muscle cell (SMC) count. The fixed effects part of the model tests the effect of WSSTmeanMax on SMC count. Random effects part of the model tests the effect of the data being nested with the patient domain.
The fixed effects part of the model (i.e. the effect of TAWSSMax on tissue thickness) was found to be non-significant (coef -0.016, 95% CI [-0.063, 0.031], p = 0.496) (Table 5). Thus, despite the linear relationship observed earlier between WSS and thickness, the data variance evident at the patient level remains significant.
Mechanical Properties: MTS and UTS.
In the case of UTS, the results trended towards a higher UTS in response to higher WSS, although this relationship was non-significant (p = 0.08). There was no evidence of influence of the patient or aortic tissue orientation on the WSS-UTS relationship. TAWSSmax was found to have no significant influence on MTS (p = 0.115).
Microstructural Features: SMC Count
TAWSSmax was found to have a persistently negative influence on SMC count [coef −4.87, 95% CI (−9.26, −0.48), p = 0.030]. The multi-level model found WSSTmeanMax to be a stronger fit to the SMC count data than patient variance (Table 6). Multi-level mixed effects regression however found no significant association between the evaluated WSS parameters and elastin/collagen.
Computational modelling of thoracic aortic disease has expanded the repertoire of aneurysm diagnostics in recent years. Whilst several in vivo parameters related to aortic wall mechanics can be obtained from such models, their association with aortic wall mechanobiology, and in the pathogenesis of ATAA remains poorly understood. Further developments in this field may allow disease severity, progression, and acute events to be predicted at an individual patient level. This is required as isolated size measurements of the aorta are inadequate (Pape et al., 2007). Existing studies have associated abnormal WSS in the ATAA with aneurysm growth and wall degeneration (Geiger et al., 2012; Guzzardi et al., 2015; Condemi et al., 2017; Bollache et al., 2018). This study builds on such work, utilizing state of the art imaging, modelling, and biomechanical methods. However, because of the availability of pathologic human aortic specimens and the ability to assess their local properties, we made a more concerted effort to divide the aortic anatomy into finer areas for analysis, specifically including up to six regions circumferentially and a further upper and lower division (up to 12 in total).
Elevated WSS on the Outer Curve
The WSS distribution in the ascending aorta is dominated by the curvature of the aortic arch, which forces blood emanating from the heart to change its direction. Fluid at the centre of the vessel is more difficult to displace as it is travelling at a higher velocity compared to that closer to the wall, so it is displaced to a greater degree (Caballero et al., 2013). Thus, blood is skewed towards the outer curvature of the bend.
WSS and Wall Thickness
Owing to the low temporo-spatial resolution of cross-sectional imaging and difficulty in tracking the motion of the aortic wall (<2.5 mm thickness), estimations of stress distribution within the aneurysm wall are reliant on important assumptions of aortic wall thickness (commonly assumed to be constant) and material properties (commonly not patient- or even disease-specific) (Martin et al., 2015; Gomez et al., 2020). Our results have identified a potential link between areas of high WSS and aortic wall thinning. Reduced wall thickness in areas of high WSS has been previously shown in femoral artery bifurcation by Kornet et al. (1999), who explained this finding with the link between low WSS, increased influx of substances into the aortic wall (through an increase in blood residence time) and increased release of vasoactive molecules potentially causing thickness increase of the intima-media layer.
These results represents a step towards the incorporation of wall thickness inferences from baseline imaging that can be incorporated into fluid-structure interaction models and estimations of wall stress distribution. These results also help to strengthen the flow-mediated degeneration hypothesis, whereby persistently elevated WSS throughout the cardiac cycle exposes mechanocytes within the aortic wall to prolonged stimuli and downstream potentially maladaptive remodeling (Humphrey et al., 2015). However, the conducted hierarchical analysis also suggests that data variance at the patient level remains significant.
WSS and Aortic Wall Mechanical Properties
Our analyses on the effect of WSS on wall mechanical properties produced two main findings. Firstly, higher WSS was associated with a reduced DEF in the longitudinal direction. This contrasted with circumferential DEF, which showed an insignificant increasing trend. This further exemplifies the anisotropic nature of the aortic wall and suggests that it extends to influence both tear direction and location in aneurysm dissection (Manopoulos et al., 2018). Secondly, the aortic wall had significantly higher values for UTS with higher WSS, indicating a reduced likelihood for wall rupture in response to elevated WSS. In addition, the aorta might be stiffer in response to high WSS, as suggested by elevated values for MTS. This result, however, was not statistically significant.
WSS and Microstructural Features: Implications for Mechanobiology
Higher WSS levels were found to be associated with a reduction in both elastin abundance and smooth muscle cell (SMC) count. Loss of elastin integrity and relative increase in collagen describes the reduced compliance and increased stiffness of the aorta seen in ageing (Wagenseil et al., 2009). Studies have shown that global increases in vascular structural stiffness reflect increased central pulse pressures and pulse wave velocities (Lacolley et al., 2009; Safar et al., 2010) that pathologically increase proximal aortic loading. Collagen deposition in many cases of ATAA has been observed to be higher, reflecting the compensatory fibrotic changes as a result of the disease process (Wågsäter et al., 2013; Meng et al., 2014). The microstructural response to elevated WSS therefore explains these changes (Humphrey et al., 2015), and describes the resulting aortic wall stiffness and delamination potential, thus increasing the likelihood of dissection (Phillippi et al., 2011).
In addition to endothelial cells and fibroblasts, SMCs are part of the repertoire of crucial mechano-sensing and -regulating cells in the aortic wall. These cells display adaptive remodeling in response to shear stress encountered by endothelial cells, which they detect with integrins, glyclocalyx, membrane microdomains, cytoskeleton, receptor tyrosine kinases and others. The direct link between WSS level and SMC found in the present study is likely to result from mechanodysregulation in response to elevated WSS, involving a disruption of cell matrix connections which are vital to aortic wall integrity (Michel et al., 2018). Dysfunctional mechanosensing can lead to cellular apoptosis and/or an atrophic remodeling response, thus disturbing the structural integrity of the aortic wall (Leung et al., 1976). Whilst there are a number of possible intracellular and matrix signaling pathways associated with the process that we have not tested for, they are likely to culminate in a final common pathway, leading to cellular loss and matrix degeneration.
Strengths and Limitations
This study benefits from a robust method of ATAA flow-to-tissue patient-specific association, which arises from segmental aneurysm analysis. Conducting the tissue characterization and CFD portions of work separately has helped reduce the risk of bias, that could result from basing aortic tissue acquisition on findings from flow analysis retrospectively (Guzzardi et al., 2015). Our statistical methods have aimed to appreciate the spread of data and made use of multilevel regression models, which have not been utilized in similar studies. Altogether this aims to improve the validity of the findings.
The study is limited by its small sample size. This may perhaps explain the lack of significance in some relationships. Larger studies would also help deal with potential confounders such patient covariates and valve function. In addition, CFD simulations were conducted under rigid wall assumption (i.e. aortic wall compliance was not taken into account). Including aortic wall compliance could further improve WSS estimation. However, this would significantly increase the computational time required to conduct patient-specific simulation further decreasing the likelihood of adoption of this technique into the clinic. Rigid wall is therefore a common assumption is several computational studies (Morbiducci et al., 2013) with a translational goal. In addition, we do not expect WSS results to be significantly affected by this assumption. Firstly, because our mechanical test results suggest increased stiffness in regions of enhanced WSS. Secondly, previous imaging studies have shown a reduction in aneurysm wall compliance when compared with healthy tissues. Larger sample sizes and the incorporation of computational methods to couple flow and wall material properties will form the basis of future work.
Our findings make a strong case for the co-localization of elevated WSS and patterns of medial degeneration, the hallmark of TAA disease. This is likely a result of the negative remodeling process of the aortic wall in response to chronic exposure to locally higher shear forces. Detailed wall shear stress analysis could predict areas of altered vascular wall mechanics and microstructural features in ascending aortic aneurysms. Presented findings further validate 4D-flow MRI and computational fluid dynamics as powerful tools for risk stratification in aneurysmal disease. This can improve precision in the timing and planning of intervention in at-risk patients.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The studies involving human participants were reviewed and approved by (17/NI/0160) the Health Research Authority (HRA) in the United Kingdom and was sponsored by the Imperial College London Joint Research and Compliance Office, as defined under the sponsorship requirements of the Research Governance Framework (2005). The patients/participants provided their written informed consent to participate in this study.
MS: conception, recruitment, analysis, statistics, writing, review; SP: conception, analysis, writing, review; SS: analysis, writing, review; SF: analysis, review; AR: review; OJ: writing, review; DO: writing, review; AO: review; JM: writing, review; XX: conception, review; TA: conception, review.
This work was supported by the NIHR Imperial College BRC (P69559) and the British Heart Foundation Centre for Research Excellence (Imperial College) (RE/18/4/34215). DPO is funded by the Medical Research Council.
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.
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.
The authors would like to thank additional surgeons who agreed to participate in this study (Royal Brompton and Harefield) Fabio De Robertis, Shahzad Raja, Toufan Bahrami, (Hammersmith Hospital) Jon Anderson, Andrew Chukwuemeka, (Barts Heart Centre) John Yap, Rakesh Uppal, Kulvinder Lall.The authors also thank staff at the Robert Steiner MRI Unit who conducted the 4D-flow MRI sequences, including Ben Statton, Marina Quinlan, Alaine Berry and Faiza Ahmed. We also thank the Cardiovascular Histopathology Unit who helped in the preparation and staining of pathology slides, including Pratibha Shah, Alex Bowman and Toyin Adefila-Ideozu.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2021.750656/full#supplementary-material
Supplementary Figure 1S | Example of unstructured mesh used for geometry discretization. The image shows details of the tetrahedral core and 10 prismatic layers at the wall (A) and of local mesh refinements. These were prescribed at the arch branches and the aneurysm wall, where a finer mesh was designed where higher velocity gradients were observed from 4D flow MRI data. Smooth transition between finer and coarser regions was ensured by prescribing small growth ratios and verified through mesh expansion factor (Ansys ICEM) (B).
Armour, C. H., Guo, B., Pirola, S., Saitta, S., Liu, Y., Dong, Z., et al. (2021). The Influence of Inlet Velocity Profile on Predicted Flow in Type B Aortic Dissection. Bio. Model. Mech. 20, 481–490. doi:10.1007/s10237-020-01395-4
Bieging, E. T., Frydrychowicz, A., Wentland, A., Landgraf, B. R., Johnson, K. M., Wieben, O., et al. (2011). In Vivo Three-Dimensional MR wall Shear Stress Estimation in Ascending Aortic Dilatation. J. Magn. Reson. Imaging 33 (3), 589–597. doi:10.1002/jmri.22485
Bollache, E., Guzzardi, D. G., Sattari, S., Olsen, K. E., Di Martino, E. S., Malaisrie, S. C., et al. (2018). 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. 156 (6), 2112–2120. doi:10.1016/j.jtcvs.2018.05.095
Bürk, J., Blanke, P., Stankovic, Z., Barker, A., Frydrychowicz, A., Russe, J., et al. (2012). 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. 14 (1), 84. doi:10.1186/1532-429x-14-84
Condemi, F., Campisi, S., Viallon, M., Troalen, T., Xuexin, G., Barker, A. J., et al. (2017). Fluid- and Biomechanical Analysis of Ascending Thoracic Aorta Aneurysm with Concomitant Aortic Insufficiency. Ann. Biomed. Eng. 45 (12), 2921–2932. doi:10.1007/s10439-017-1913-6
Erbel, R., Aboyans, V., Boileau, C., Bossone, E., Bartolomeo, R. Di., Eggebrecht, H., et al. (2014). 2014 ESC Guidelines on the Diagnosis and Treatment of Aortic Diseases. Eur. Heart J. 35 (41), 2873–2926. doi:10.1093/eurheartj/ehu281
Geiger, J., Markl, M., Herzer, L., Hirtler, D., Loeffelbein, F., and Stiller, B. (2012). Aortic Flow Patterns in Patients with Marfan Syndrome Assessed by Flow-Sensitive Four-Dimensional MRI. J. Magn. Reson. Imaging 35 (3), 594–600. doi:10.1002/jmri.23500
Gomez, A., Wang, Z., Xuan, Y., Wisneski, A. D., Hope, M. D., Saloner, D. A., et al. (2020). Wall Stress Distribution in Bicuspid Aortic Valve–Associated Ascending Thoracic Aortic Aneurysms. Ann. Thorac. Surg. 110 (3), 807–814. doi:10.1016/j.athoracsur.2019.12.035
Guzzardi, D. G., Barker, A. J., Van Ooij, P., Malaisrie, S. C., Puthumana, J. J., Belke, D. D., et al. (2015). Valve-Related Hemodynamics Mediate Human Bicuspid Aortopathy: Insights from Wall Shear Stress Mapping. J. Am. Coll. Cardiol. 66 (8), 892–900. doi:10.1016/j.jacc.2015.06.1310
Hiratzka, L. F., Bakris, G. L., Beckman, J. A., Bersin, R. M., Carr, V. F., Casey, D. E., et al. (2010). ACCF/AHA/AATS/ACR/ASA/SCA/SCAI/SIR/STS/SVM Guidelines for the Diagnosis and Management of Patients with Thoracic Aortic Disease. Circulation 121 (13). doi:10.1161/CIR.0b013e3181d4739e
Howard, D. P. J., Banerjee, A., Fairhead, J. F., Perkins, J., Silver, L. E., and Rothwell, P. M. (2013). Population-Based Study of Incidence and Outcome of Acute Aortic Dissection and Pre-morbid Risk-Factor Control: 10-Year Results from the Oxford Vascular Study. Circulation 127 (20), 2031–2037. doi:10.1161/circulationaha.112.000483
Humphrey, J. D., Schwartz, M. A., Tellides, G., and Milewicz, D. M. (2015). Role of Mechanotransduction in Vascular Biology: Focus on Thoracic Aortic Aneurysms and Dissections. Circ. Res. 116 (8), 1448–1461. doi:10.1161/circresaha.114.304936
Kornet, L., Hoeks, A. P. G., Lambregts, J., and Reneman, R. S. (1999). In the Femoral Artery Bifurcation, Differences in Mean wall Shear Stress within Subjects Are Associated with Different Intima-media Thicknesses. Arterioscler Thromb. Vasc. Biol. 19 (12), 2933–2939. doi:10.1161/01.atv.19.12.2933
LaDisa, J. F., Figueroa, C. A., Vignon-Clementel, I. E., Kim, H. J., Xiao, N., Ellwein, L. M., et al. (2011). Computational Simulations for Aortic Coarctation: Representative Results from a Sampling of Patients. J. Biomech. Eng. 133, 091008. doi:10.1115/1.4004996
Les, A. S., Shadden, S. C., Figueroa, C. A., Park, J. M., Tedesco, M. M., Herfkens, R. J., et al. (2010). Quantification of Hemodynamics in Abdominal Aortic Aneurysms during Rest and Exercise Using Magnetic Resonance Imaging and Computational Fluid Dynamics. Ann. Biomed. Eng. 38, 1288–1313. doi:10.1007/s10439-010-9949-x
Leung, D. Y. M., Glagov, S., and Mathews, M. B. (1976). Cyclic Stretching Stimulates Synthesis of Matrix Components by Arterial Smooth Muscle Cells In Vitro. Science 191, 475–477. doi:10.1126/science.128820
Manopoulos, C., Karathanasis, I., Kouerinis, I., Angouras, D. C., Lazaris, A., Tsangaris, S., et al. Identification of Regional/layer Differences in Failure Properties and Thickness as Important Biomechanical Factors Responsible for the Initiation of Aortic Dissections. J. Biomech. 2018. 80:102–110. doi:10.1016/j.jbiomech.2018.08.024
Martin, C., Sun, W., and Elefteriades, J. (2015). Patient-specific Finite Element Analysis of Ascending Aorta Aneurysms. Am. J. Physiol. - Hear. Circ. Physiol. 308 (10), H1306–H1316. doi:10.1152/ajpheart.00908.2014
Meng, Y., Tian, C., Liu, L., Wang, L., and Chang, Q. (2014). Elevated Expression of Connective Tissue Growth Factor, Osteopontin and Increased Collagen Content in Human Ascending Thoracic Aortic Aneurysms. Vascular 22 (1), 20–27. doi:10.1177/1708538112472282
Michel, J. B., Jondeau, G., and Milewicz, Di. M. (2018). From Genetics to Response to Injury: Vascular Smooth Muscle Cells in Aneurysms and Dissections of the Ascending Aorta. Cardiovasc. Res. 114 (4), 578–589. doi:10.1093/cvr/cvy006
Morbiducci, U., Ponzini, R., Gallo, D., Bignardi, C., and Rizzo, G. (2013). Inflow Boundary Conditions for Image-Based Computational Hemodynamics: Impact of Idealized versus Measured Velocity Profiles in the Human Aorta. J. Biomech. 46 (1), 102–109. doi:10.1016/j.jbiomech.2012.10.012
Olchanyi, M. D., Sadikov, A., Frattolin, J., Sasidharan, S., Yousuf Salmasi, M., Edgar, L. T., et al. (2020). Validation of Markerless Strain-Field Optical Tracking Approach for Soft Tissue Mechanical Assessment. J. Biomech. 116 (2021), 110196. doi:10.1016/j.jbiomech.2020.110196
Pape, L. A., Tsai, T. T., Isselbacher, E. M., Oh, J. K., O’Gara, P. T., Evangelista, A., et al. (2007). Aortic Diameter ≥5.5 Cm Is Not a Good Predictor of Type A Aortic Dissection: Observations from the International Registry of Acute Aortic Dissection (IRAD). Circulation 116 (10), 1120–1127. doi:10.1161/circulationaha.107.702720
Phillippi, J. A., Pasta, S., and Vorp, D. A. (2011). “Biomechanics and Pathobiology of Aortic Aneurysms,” in Biomechanics and Mechanobiology of Aneurysms (New York, NY, USA: Springer). doi:10.1007/8415_2011_84
Piatti, F., Sturla, F., Bissell, M. M., Pirola, S., Lombardi, M., Nesteruk, I., et al. (2017). 4D Flow Analysis of BAV-Related Fluid-Dynamic Alterations: Evidences of wall Shear Stress Alterations in Absence of Clinically-Relevant Aortic Anatomical Remodeling. Front. Physiol. 8, 441. doi:10.3389/fphys.2017.00441
Pirola, S., Cheng, Z., Jarral, O. A., O'Regan, D. P., Pepper, J. R., Athanasiou, T., et al. (2017). On the Choice of Outlet Boundary Conditions for Patient-specific Analysis of Aortic Flow Using Computational Fluid Dynamics. J. Biomech. 60, 15–21. doi:10.1016/j.jbiomech.2017.06.005
Pirola, S., Guo, B., Menichini, C., Saitta, S., Fu, W., Dong, Z., et al. (2019). 4-D Flow MRI-Based Computational Analysis of Blood Flow in Patient-specific Aortic Dissection. IEEE Trans. Biomed. Eng. 66 (12), 3411–3419. doi:10.1109/tbme.2019.2904885
Reymond, P., Merenda, F., Perren, F., Rüfenacht, D., and Stergiopulos, N. (2009). Validation of a One-Dimensional Model of the Systemic Arterial Tree. Am. J. Physiol.-Heart Circulatory Physiol. 297, H208–H222. doi:10.1152/ajpheart.00037.2009
Stewart, B., Khanduri, P., McCord, C., Ohene-Yeboah, M., Uranues, S., Vega Rivera, F., et al. (2014). Global Disease burden of Conditions Requiring Emergency Surgery. Br. J. Surg. 101, e9–22. doi:10.1002/bjs.9329
Simon, A. C., Safar, M. E., Levenson, J. A., London, G. M., Levy, B. I., Chau, N. P., et al. (1979). An Evaluation of Large Arteries Compliance in Man. Am. J. Physiol. 237 (5), H550–H554. doi:10.1152/ajpheart.1979.237.5.H550
Tan, F., Borghi, A., Mohiaddin, R., Wood, N., Thom, S., and Xu, Y. X. (2009). Analysis of Flow Patterns in a Patient-specific Thoracic Aortic Aneurysm Model. Comput. Structures 87, 680–690. doi:10.1016/j.compstruc.2008.09.007
Tsai, T. T., Evangelista, A., Nienaber, C. A., Trimarchi, S., Sechtem, U., Fattori, R., et al. (2006). Long-term Survival in Patients Presenting with Type A Acute Aortic Dissection: Insights from the International Registry of Acute Aortic Dissection (IRAD). Circulation 114 (1), 350–357. doi:10.1161/CIRCULATIONAHA.105.000497
Van, O. P., Potters, W. V., Nederveen, A. J., Allen, B. D., Collins, J., Carr, J., et al. (2015). A Methodology to Detect Abnormal Relative wall Shear Stress on the Full Surface of the Thoracic Aorta Using Four-Dimensional Flow MRI. Magn. Reson. Med. 73 (3), 1216–1227. doi:10.1002/mrm.25224
Vorp, D. A., Schiro, B. J., Ehrlich, M. P., Juvonen, T. S., Ergin, M. A., and Griffith, B. P. (2003). Effect of Aneurysm on the Tensile Strength and Biomechanical Behavior of the Ascending Thoracic Aorta. Ann. Thorac. Surg. 75 (4), 1210–1214. doi:10.1016/s0003-4975(02)04711-2
Wågsäter, D., Paloschi, V., Hanemaaijer, R., Hultenby, K., Bank, R. A., Franco-Cereceda, A., et al. (2013). Impaired Collagen Biosynthesis and Cross-Linking in Aorta of Patients with Bicuspid Aortic Valve. J. Am. Heart Assoc. 2 (1), e000034. doi:10.1161/JAHA.112.000034
Keywords: aortic surgery, aneurysm, computational fluid dynamics, CFD, magnetic resonance imaging, wall shear stress, computational pathology, vascular biomechanics
Citation: Salmasi MY, Pirola S, Sasidharan S, Fisichella SM, Redaelli A, Jarral OA, O’Regan DP, Oo AY, Moore JE, Xu XY and Athanasiou T (2021) High Wall Shear Stress can Predict Wall Degradation in Ascending Aortic Aneurysms: An Integrated Biomechanics Study. Front. Bioeng. Biotechnol. 9:750656. doi: 10.3389/fbioe.2021.750656
Received: 30 July 2021; Accepted: 24 September 2021;
Published: 18 October 2021.
Edited by:Jolanda Wentzel, Erasmus Medical Center, Netherlands
Reviewed by:Natalya Kizilova, Warsaw University of Technology, Poland
Harvey Ho, The University of Auckland, New Zealand
Copyright © 2021 Salmasi, Pirola, Sasidharan, Fisichella, Redaelli, Jarral, O’Regan, Oo, Moore, Xu and Athanasiou. 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: M. Yousuf Salmasi, firstname.lastname@example.org
†These authors have contributed equally to this work and share first authorship