Impact Factor 6.050 | CiteScore 7.8
More on impact ›


Front. Cardiovasc. Med., 01 June 2021 |

Hemodynamic Abnormalities in the Aorta of Turner Syndrome Girls

  • 1Department of Biomedical Engineering, University of Strathclyde, Glasgow, United Kingdom
  • 2Department of Radiology, Royal Hospital for Children, Glasgow, United Kingdom
  • 3Department of Paediatric Endocrinology, Royal Hospital for Children, Queen Elizabeth University Hospital, Glasgow, United Kingdom

Congenital abnormalities in girls and women with Turner syndrome (TS), alongside an underlying predisposition to obesity and hypertension, contribute to an increased risk of cardiovascular disease and ultimately reduced life expectancy. We observe that children with TS present a greater variance in aortic arch morphology than their healthy counterparts, and hypothesize that their hemodynamics is also different. In this study, computational fluid dynamic (CFD) simulations were performed for four TS girls, and three age-matched healthy girls, using patient-specific inlet boundary conditions, obtained from phase-contrast MRI data. The visualization of multidirectional blood flow revealed an increase in vortical flow in the arch, supra-aortic vessels, and descending aorta, and a correlation between the presence of aortic abnormalities and disturbed flow. Compared to the relatively homogeneous pattern of time-averaged wall shear stress (TAWSS) on the healthy aortae, a highly heterogeneous distribution with elevated TAWSS values was observed in the TS geometries. Visualization of further shear stress parameters, such as oscillatory shear index (OSI), normalized relative residence time (RRTn), and transverse WSS (transWSS), revealed dissimilar heterogeneity in the oscillatory and multidirectional nature of the aortic flow. Taking into account the young age of our TS cohort (average age 13 ± 2 years) and their obesity level (75% were obese or overweight), which is believed to accelerate the initiation and progression of endothelial dysfunction, these findings may be an indication of atherosclerotic disease manifesting earlier in life in TS patients. Age, obesity and aortic morphology may, therefore, play a key role in assessing cardiovascular risk in TS children.


Turner syndrome (TS) is a rare genetic disorder where the second sex chromosome in females is partially or completely absent and can affect all or only a percentage of cells (mosaicism). Although rare, TS is the most common chromosomal abnormality among females, affecting 1 in 2,500 live births based on epidemiological data from Europe, East Asia, and North America (13). The clinical characteristics of TS are highly variable, however, with congenital heart abnormalities estimated to occur in as many as half of individuals (4). These defects predominantly affect the left side of the heart, with the most commonly reported being an elongated transverse aortic arch (ETA) (5, 6). Defined as an increased distance between the second (LCCA) and third (LSA) supra-aortic branch origins, ETA is reported in 49% of TS adults (5, 6). Bicuspid aortic valve (BAV) is the second most common abnormality in TS, affecting ~30%, closely followed by aortic dilatation with a prevalence of 27% (5, 7, 8). Coarctation of the aorta (CoA) is reported in around 16% of TS females, most commonly at the site between the third supra-aortic branch and the descending aorta (5, 7, 8). Aberrant right subclavian artery (RSA) is an anatomical variation of the RSA which atypically originates from the arch as a separate fourth branch and has a prevalence of 8% in TS (6). These congenital abnormalities, alongside an underlying predisposition to obesity and hypertension, contribute to a greater risk of cardiovascular disease and ultimately reduced life expectancy in TS (9).

Hemodynamic factors have been linked to the initiation and development of cardiovascular disease for over a century (10). However, the exact nature of pro-atherogenic flow is uncertain with researchers proposing contradicting theories (11, 12). Fry suggested that high wall shear stress (WSS) preceded endothelial dysfunction, one of the early biological markers for atherosclerotic lesions that underlie most cardiovascular diseases (11). Soon after, Caro et al. suggested that high WSS regions are in fact spared from disease, with low WSS areas prone to develop atherosclerosis (12, 13). The oscillatory shear index (OSI) was then put forward by Ku et al. (14) to characterize regions of reversing flow, and today the combined low and/or oscillatory WSS theory is generally accepted as the biological mechanism for atherosclerosis. However, atherosclerotic lesions appear to depend on and vary with age (15): studies in human fetuses, newborns, and children have demonstrated the development of sudanophilic lesions downstream of branch ostia in the thoracic aorta (16); in young adults, lesions are observed laterally of such branch origins, while in middle age, upstream (17, 18); last, in older people, atherosclerotic lesions develop around the origins of intercostal arteries (19). Some aspects of such distinct lesion distributions have been described before hemodynamically and were correlated with strikingly varied WSS patterns calculated numerically around aortic branches (2022).

In recent decades, computational fluid dynamic (CFD) modeling has gradually replaced traditional experimental methods in the study of cardiovascular diseases, largely due to its ability to quantify variables not measurable in-vivo, particularly at a temporal and spatial resolution exceeding all other methodologies (23). Increasingly complex models of arterial vessels have improved our understanding of the relationship between anatomical and hemodynamic factors (2431). However, the accuracy of these models is dependent on several assumptions. The assumption of blood as an incompressible, Newtonian fluid is common in the case of larger arteries. While this assumption is not physiologically correct, qualitative wall shear stress patterns have proven to be comparable between Newtonian and non-Newtonian simulations (32). Multiple studies ignore the pulsatile nature of blood flow, assuming instead a steady-state condition (25, 26). This assumption is computationally-efficient when only considering the mean WSS result, but pulsatility is crucial for capturing time-dependent parameters such as instantaneous velocities and WSS, time-averaged wall shear stress (TAWSS), transverse WSS (transWSS), oscillatory shear index (OSI), and relative residence time (RRT) (33). In addition to the above assumptions, the computational result is highly sensitive to the assigned boundary conditions (BCs) (3236). Madhavan and Kemmerling (32) compared five different inlet velocity profiles on human aortae and found only small differences in the flow solution approximately two diameters downstream from the aortic inlet. In a similar study on mice, Van Doormaal et al. (37) used realistic MRI-derived aortic root velocity profiles and suggested that idealized inflow profiles should be avoided. Pirola et al. (28) recommended the use of a 3D inlet velocity profile for hemodynamic analysis of the ascending aorta and arch, but a 1D inlet velocity profile was acceptable for evaluating flow in the descending aorta. Outlet BCs impact a greater percentage of the solution domain, but often in-vivo data is not available, which is a hurdle in patient-specific simulations (32). In the absence of in-vivo data, most studies apply either an outflow boundary condition, in which a percentage of the total flow is specified at each outlet, or a simple downstream resistance or the Windkessel model, in which the resistance and capacitance of the downstream vasculature are modeled (38). However, implementation of the Windkessel model is particularly challenging for children, due to the lack of values for the Windkessel parameters in the literature, especially for children with TS. In the absence of patient information, Murray's Law or the splitting method, both established from the power law relationship between branch diameter and flow rate, could be used (39, 40). Several authors have favored Murray's law over the zero-pressure boundary condition, for capturing more physiologically relevant flow features (33, 39).

In Turner syndrome, deviations from an anatomically healthy aorta are common, and therefore changes in blood flow may exacerbate the risk of cardiovascular disease. Understanding the pathogenesis of the increased cardiovascular morbidity and mortality in Turner syndrome, and determining the contributions of atherosclerotic disease, hypertension, and obesity was made a key cardiovascular research priority by the American Heart Association (4). To this end, our research primarily aims at providing new evidence for the aortic arch hemodynamics in the TS children population, with categorically abnormal aortic morphologies. Comparison with three age- and gender-matched, anatomically healthy cases is also made to highlight the differences with normal anatomies. The hemodynamic environment of each patient was studied, and the relationship between aortic morphology and flow was analyzed. Moreover, this study provides further evidence on understanding the aortic flow development generally in children and attempts to find correlations between hemodynamics and clinical significance in TS patients.

Materials and Methods

Magnetic Resonance Imaging (MRI) and Patient Cohort

In this study, retrospective MRI scans were obtained from four (n = 4) girls (average age 13 ± 2 years) with karyotypically proven Turner syndrome (Table 1), attending the pediatric TS clinic at the Royal Hospital for Children, Queen Elizabeth University Hospital (RHC, QEUH). TS3 underwent anomalous pulmonary venous drainage repair ~6 years prior to MRI imaging, and TS4 underwent left congenital diaphragmatic hernia repair at birth (13 years before the MRI scan). Both surgical corrections were unrelated to the region of interest, being the proximal aorta.


Table 1. Biometric and anatomical data for the Turner syndrome (TS) girls (n = 4).

All cardiac imaging was performed between 2014 and 2018 on a 1.5 Tesla diagnostic MRI scanner (MAGNATOM Aera/Avanto, Siemens Healthcare, GmbH). The scans were acquired with both ECG and respiratory gating. Additional 2D time-resolved phase-contrast MRI (PC-MRI) data was acquired at the aortic root of each individual over a cardiac cycle, oriented axially and centered on the anatomical position of the pulmonary artery. Flow in the through-plane direction was measured with PC-MRI imaging using a velocity encoding (VENC) of 150–400 cm/s and TE = 2.66–3.33 ms, TR = 9.94 ms (slice thickness of 3.5–5 mm; FOV = 300 × 300 mm), for all data, resulting in a 256 x 256 pixel matrix and resolution of 1.17 × 1.17 × 3.5–5 mm.

Biometric data (height, weight and blood pressure, BP) was obtained within 1 month of the cardiac imaging data and body mass index (BMI) and BMI percentile were computed according to the 1990 British childhood standards (Table 1) (41). The cardiac output (CO) and heart rate (HR) were obtained from the PC-MRI data. Patient-specific volumetric image and geometrical data from three gender-matched (average age 13 ± 6.8 years) patients (n = 3), with no known aortic abnormalities, were purchased from the Vascular Model Repository ( The average CO and HR of the healthy group were 3.77 ± 1 l/min (mean ± SD), and 81.33 ± 22.59 bpm (mean ± SD), respectively. Even though detailed BMI information was not included as part of the purchased data, it is expected that all healthy cases were within the normal BMI range, in accordance with the standard of the repository.

Anatomical Reconstruction and Mesh Generation

The three-dimensional geometries of the TS group were segmented and reconstructed from the MRI data, from above the location of the aortic valve to the end of the thoracic aorta, including the brachiocephalic, left common carotid, and subclavian arteries, using the medical-imaging software ITK-SNAP ( In Figure 1 the anterior view of the reconstructed aortic models of the healthy (H1–H3) and Turner syndrome girls (TS1–TS4) is shown, with the geometries being in scale for direct comparison. The healthy controls had aortic diameters (D) of 17.53-22.48 mm at the sinotubular junction (Table 2), similar to the mean diameter of 17.5 mm (range = 11.1–26.4 mm) found in healthy children of a related age (n = 53, range = 2–20 years) (42). Ascending aorta diameters in women with TS are generally about 10% greater than in healthy women (43). In our cohort, the TS girls had aortic diameters of 18.99–33.77 mm (Table 2). The excessively large aortic diameter in TS1 was expected due to the dilated ascending aorta. Among the Turner syndrome models were aortic abnormalities commonly reported in the literature, such as ascending aorta dilatation (TS1), elongated transverse arch (TS2 and TS3), and an aberrant right subclavian artery (TS3), as stated in the patient clinical records. Arch morphology in TS4 was suggestive of a triangular-shaped “gothic” arch, although this was not recorded in the clinical records. Both TS1 and TS3 had a bicuspid aortic valve, while TS2 and TS4 had a normal functioning tricuspid aortic valve. TS girls 1, 2, and 4 exhibited the conventional anatomy of three supra-aortic vessels arising from the arch. In TS3, the right subclavian artery arose from the posterior arch, distal to the left common carotid artery (Figure 1, inset). There were no aortic abnormalities in the healthy cases (H1–H3).


Figure 1. Anterior view of the reconstructed aortic models from the (H1-H3) healthy and (TS1-TS4) Turner syndrome (TS) girls. RSA, Right subclavian artery; RCCA, right common carotid artery; LCCA, left common carotid artery; LSA, left subclavian artery. Inset: Superior view of TS3 to highlight the origin of the aberrant RSA. Geometries are in scale.


Table 2. Hemodynamic information for healthy (H1–H3) and Turner syndrome (TS1–TS4) patients.

After segmentation, the surface models were smoothed in Autodesk Meshmixer ( to reduce post-segmentation artifacts, and flow extensions were added normal to all boundary faces in VMTK ( An extension of half aortic diameter in length was added at each patient-specific inlet, which was fitted to a circular inlet of the same area. The domain was then discretized in STAR-CCM+ software (Siemens PLM, USA, using polyhedral elements for the internal mesh, and prismatic elements for the boundary layer (Figure 2). While tetrahedral meshing is a more common approach, often a larger number of elements are required to limit element skewness and achieve acceptable mesh quality. Polyhedral meshing, introduced in STAR-CCM+, has two major benefits resulting in better numerical stability: one, each individual element has multiple neighbors, and two, the elements are less sensitive to stretching (44, 45). As a result, an accurate solution can be achieved with a much lower cell count of polyhedral elements than tetrahedral. Quantitative grid convergence was investigated on three meshes of increasing density (each grid was further refined by a factor of 2) and assessed using the Grid Convergence Index (GCI=Fs(e/(rp-1)) where Fs is a safety factor, taken as 1.25, e the relative error, r the mesh refinement ratio, and p the order of convergence) (46). The results confirmed that wall shear stress accuracy was suitably achieved with 1.5 million polyhedral elements (GCI ≅ 11%). Mesh generation with polyhedral elements is a promising discretization approach to reduce the computational time of simulations while retaining accuracy.


Figure 2. STAR-CCM+ polyhedral mesh shown on a healthy (H1) model with zoomed views of the inlet mesh with a prism boundary layer (left) and the arch surface mesh (right).

To accurately resolve the flow adjacent to the wall, the boundary layer mesh was generated using a total of 5 layers, with the height of the first layer (Δy1) satisfying a y+ value of 1 (47), according to Equation (1). The y+ value describes a non-dimensional distance from the wall to the first element node and therefore characterizes near-wall mesh quality. Refining the near-wall mesh with an appropriate y+ value is crucial for resolving flow in the viscous sublayer of the boundary layer. The values of the first boundary layer height (Δy1) and total boundary layer height (Δyn = 5) are shown in Table 2.

Δy1=μy+ρ UT    (1)

where Δy1 is the height of the first layer, y+ is equal to 1, μ and ρ are the fluid viscosity and density, taken as 3.5e-3 Pa s and 1050 kg/m3, respectively (25), and UT is the friction velocity computed from Equation (2):

UT=τwρ    (2)

where τw is the wall shear stress, calculated according to Equation (3), Umean is the mean velocity value corresponding to Qmean (Table 2), and Cf is the skin friction coefficient calculated per Equation (4) (48).

τw=12CfρUmean2    (3)
Cf=2[(8Re)12+((2.457ln((Re7)0.9))16+(37530Re)16)1.5]1/12    (4)

The adequacy of the y+ value was further confirmed from Equations (1, 2), based on the directly calculated maximum WSS values, which led to a y+ value much smaller than 1.

Boundary Conditions

The flow exiting the aortic valve was segmented from the two-dimensional time-resolved PC-MRI data with the use of Medviso Segment software (, resulting in a series of time-dependent flow waveforms (Figure 3). The highest peak flow rate was found for TS2 and TS3, and the lowest peak flow rate for H1 (Figure 3A). To account for diameter and cardiac cycle variability, the normalized flow rates were calculated (Figure 3B), where two of the healthy girls demonstrated higher peak flow rates than the TS girls. The subject-specific volumetric flow waveforms (Figure 3A) were applied at the inlet boundary, which was defined at the location of the sinotubular junction.


Figure 3. PC-MRI derived (A) volumetric flow rate and (B) normalized waveforms at the aortic root during one cardiac cycle for healthy (H) and Turner syndrome (TS) girls. Insets: average data calculated from H1-H3 (black line) and TS1-TS4 (red line). Flow rate and time normalized by the mean flow rate and cardiac cycle period (see Table 2 for values).

The cardiac cycle period (T) and mean and peak flow rates (Qmean, Qpeak, respectively) were extracted from the location of the aortic valve from the PC-MRI data, while the Reynolds and Womersley numbers were calculated post-segmentation using the patient-specific inlet diameter (D) and the cardiac cycle period (Table 2). The mean and peak Reynolds number (Re = 4Qρ/(πDμ) in TS patients varied, although the mean value for the TS group (Remean=1296) was greater than the healthy group (Remean = 1,182). Available Repeak data in the literature for TS children reported a range of 3,980–6,560 in the ascending aorta, meaning that the average value for our cohort (Remean = 4,568) was at the lower end of this range (49). The Womersley number is a non-dimensional expression of the pulsatile nature of blood flow, and was calculated based on Equation (5):

Wo=D2 2 π f ρμ    (5)

where f is the frequency of the cardiac cycle (s−1). The expected Womersley number in the abdominal aorta of a young healthy adult is 13 (50). At large Womersley numbers (above 10), the shape of the velocity profile is relatively flat or plug-like, with the maximum velocity no longer at the center (51). The computed Womersley number (Table 2) for the healthy aortae (Wo = 15–16) is within the range of TS2-TS4 (Wo = 15–18). For TS1, with the largest aortic diameter, the Womersley number is the highest with a value of 28.

The percentage of the total flow rate distributed to each outlet was calculated using Murray's law due to the shortage of patient-specific or literature values for children. In arterial bifurcations, Murray's law states that the flow is proportional to the diameter of that vessel raised to a power, n (39). For the right subclavian branch, it can be expressed as:


where QRSA, QRCCA, QLCCA, QLSA, QDescAo are the flow rates, and DRSA, DRCCA, DLCCA, DLSA, DDescAo the diameters at the corresponding vessels. This relationship has been used in several human CFD studies, and is well-known with an exponential power of 3, hence the name “Murray's cube law.” However, several authors have shown that a power of 2 is more valid in the case of larger vessels such as the aorta (5254). Using a power of 2, the calculated combined flow percentage to the aortic branches was 34–59% for the healthy aortae and 32–55% for the TS aortae, with the remaining flow distributed to the descending aorta (Table 3).


Table 3. Outlet flow percentages for each individual case calculated using Murray's Law [exponential power of 2, Equation (6)].

In line with other CFD studies in the aortic arch, the arterial wall was presumed non-deformable and the no-slip BC was assigned (24, 25, 27, 28, 30, 31).

Numerical Method

Throughout this study, the fluid was considered incompressible and Newtonian, with a constant density and viscosity. The fluid was governed by the time-dependent Navier-Stokes equations:

·u=0    (7)
ρut+ρ(u·)u=-p +μ2u    (8)

where u is the velocity vector and p the pressure.

Flow simulations were performed in the open-source software, OpenFOAM© (, version 6), using the combined pressure-implicit split-operator and semi-implicit method for pressure-linked equations (PIMPLE) solver for incompressible, transient flow. The flow was computed using the wall-adapted local eddy-viscosity (WALE) large eddy simulation (LES) model, with temporal and spatial discretization performed using second order accurate schemes (backward Euler and central differencing, respectively). As the peak Reynolds numbers (Table 2) are in the transitional to turbulent range, an investigation was made with a pulsatile laminar model, a k-omega SST model, and the LES (WALE) model for the same mesh (TS2). The wall shear stress results were qualitatively similar for all three models, and marginally different quantitatively (0.14% difference in the integral of the WSS between the LES and laminar models, and 6% difference between the LES and k-omega SST models). However, the LES model more accurately captured the flow features, especially during systole and early diastole. The motivation for the use of the LES turbulence model was further reinforced by other published studies on blood flow (55, 56). Blood flow was simulated until time-periodicity was reached at approximately five cardiac cycles with a time step (t) of 1 × 10−3 s (satisfying mean Courant number < 1). Residual control for the convergence criteria was set to 1 × 10−5 (57) for both pressure and velocity.

Shear Stress Parameters

The characterization of shear stress during pulsatile flow is commonly described with hemodynamic parameters such as the instantaneous wall shear stress (WSS) at various time points in the cardiac cycle, the time-averaged WSS (TAWSS), the oscillatory shear index (OSI), the relative residence time (RRT) and the transverse WSS (transWSS) (Equations 9–12, respectively). Such shear stress parameters can be used to identify areas where flow departs from a laminar, unidirectional pattern.

TAWSS=1T0T|τw| dt    (9)
OSI=12(1-|1T0T τw dt|1T0T|τw| dt)    (10)
RRT=[(1-2 OSI) TAWSS]-1    (11)
transWSS=1T0T|τw.  (nx0Tτwdt|0Tτwdt|)|dt    (12)

where |τw| is the magnitude of the wall shear stress vector, and n is the surface normal In this study, the TAWSS, OSI, RRT and transWSS were taken over the fifth cardiac cycle, and both the WSS and TAWSS were normalized with respect to the mean WSS at the inlet for each individual case. The OSI describes the degree of oscillatory flow ranging from zero, representing unidirectional flow, to 0.5, representing reversing flow with no mean shear direction (58). The RRT provides information on the residence time of flow particles in close proximity to the wall and is elevated in regions of low magnitude and high oscillatory WSS (59). RRT was normalized with respect to the surface integral average for the patient-specific aortic geometries. The transWSS quantifies multidirectional flow, with low values indicating that flow remains approximately parallel to a single axis throughout the cardiac cycle and high values indicating large changes in flow direction (60).


Blood Flow in the Aorta

Blood flow was visualized in the geometries of Figure 1, using a combination of 3D streamlines (Figure 4) at three time points in the cardiac cycle, peak velocity (t1), maximum deceleration (t2), and mid-diastole (t3), and 2D and 3D vector-fields (Figure 5) at peak velocity (t1). Maximum deceleration was defined as the minimum rate of change of the velocity with respect to time.


Figure 4. Velocity streamlines in the aortic arch of the healthy (H1-H3) and Turner syndrome (TS1-TS4) girls at (A) peak velocity (t1), (B) maximum deceleration (t2), and (C) mid-diastole (t3), colored by non-dimensional velocity magnitude that is normalized according to the average inlet velocity (Umean), derived from patient PC-MRI data (anterior view). Note that the color legends in (B) and (C) were shifted compared to (A) to enhance visualization. For interpretation of the colored legends, please refer to the online version of the paper.


Figure 5. (Left) Through-plane velocity profiles and (Right) contours of through-plane velocity overlayed by vectors of in-plane velocity on seven cross-sections α-α′ to η-η′ along the aorta [locations shown on the 3D healthy (H1) model]. Contours colored by non-dimensional axial velocity at peak velocity for the healthy (H1–H3) and Turner syndrome (TS1–TS4) girls. Cross-sections are oriented looking downstream, with the top and bottom edges corresponding to the anterior and posterior sides of the aorta, respectively, and the left and right points as shown on the left. Cross-sections are not to scale.

Velocity Streamlines

In Figure 4A, the non-dimensional velocity streamlines in the aortic arch were visualized at peak velocity (t1). In all three healthy cases (H1-H3), the flow at peak systole was laminar throughout the aortic arch, the descending aorta, and the three major branches arising from the arch. A similar flow pattern was observed for TS patients 1 and 4. In TS2 and TS3, the streamlines at peak systole were, for the most part, laminar throughout the aorta, except at the entrance to the LSA in both geometries, and the aberrant RSA in TS3, where small zones of recirculation were formed. In the ascending aorta of all cases, the velocity magnitude was greater at the inner wall. The majority of flow within the arch of H1–H2, and to some extent in H3 and TS2, was of lower magnitude than in the ascending aorta. In all TS aortae, the velocity magnitude for the majority of the arch was similar to, or greater than, the flow in the ascending aorta. In the majority of TS aortic models, flow velocity in the descending aorta was high, with the exception of TS2.

In Figure 4B, the non-dimensional velocity streamlines in the aortic arch were visualized at maximum deceleration (t2). In the healthy aortae, flow was relatively laminar with the exception of the entrance to the LCCA and LSA branches in H3. In contrast, a complex flow pattern with significant secondary flows was observed for all TS aortic models. In the ascending aorta of TS girls 2–4, the flow at maximum deceleration was laminar. However, in the TS girl with ascending aortic dilatation (TS1), the flow was highly disturbed in this region with significantly low velocity values. In the aortic arch of TS3–TS4, vortical flow developed along the lesser curvature and extended to the entrance of the descending aorta, while in TS1 and TS2 the region of vortical flow filled the majority of the aortic arch. In the descending aorta of the TS girls, with the exception of TS1, vortical flow was present at the proximal wall with undisturbed streamlines at the distal wall, with the inverse true for TS1. Unlike in the healthy geometries, the pattern of flow entering the TS aortic branches was unpredictable and highly disturbed. Strong vortical flow patterns were observed in the brachiocephalic branch of TS1, and recirculation regions at the proximal wall of TS2–TS4. The left common carotid artery of TS1 and TS3, and the left subclavian branch of TS1–TS3 further exhibit disturbed flow.

In Figure 4C, the non-dimensional velocity streamlines in the aortic arch were visualized at mid-diastole (t3). For the healthy aortae, the predominantly laminar flow pattern seen at maximum deceleration (Figure 4B) was replaced with strong, slow-moving vortical motion throughout the entire vessel at mid-diastole (Figure 4C). For TS1-TS2, there was little visible difference in the flow patterns between t2 (Figure 4B) and t3 (Figure 4C), albeit with lower velocity values and a slight increase in vortical flow in the ascending aorta of TS2 and descending aorta of TS1. For TS3, flow throughout all regions of the aorta and aortic branches was clearly more disturbed in mid-diastole than in systole. Similar observations could be made for TS4, with the greatest flow disturbance seen in the ascending aorta and arch.

Through-Plane and In-Plane Velocities

Furthermore, the through-plane and in-plane velocities were calculated for seven cross-sectional slices (α-α′ to η-η′, Figure 5) of all models, at peak velocity (t1). The 3D through-plane velocity profiles were shown only for a healthy (H1) case, in the corresponding geometry (left side of Figure 5), while the cross-sections, colored by through-plane velocity contours and overlaid by in-plane velocity streamlines, were shown to the right side of Figure 5 for all models. The cross-sectional slices were considered at the same relative locations for all cases, perpendicular to the aortic centerline and relative to the individual model inlet diameter, D: slice α-α′ was taken 1D upstream from the inlet; slice β-β′ was assumed midway between the LCCA and LSA branches; γ-γ′, 1D downstream from the LSA (TS3: from the LCCA); δ-δ′, 2D downstream from the LSA (TS3: from the LCCA); ε-ε′, 0.5D upstream from the brachiocephalic junction; ζ-ζ′, 0.5D upstream from the LCCA origin; and η-η′, 0.5D upstream from the LSA ostium. The slices in Figure 5 are oriented so that the top and bottom edges correspond to the anterior and posterior sides of the aorta, respectively, while the left and right points correspond to the greater and lesser curvature of the arch, respectively, for slices α-α′ to δ-δ′, and to the outer and inner walls of the branches, for slices ε-ε′ to η-η′. That is, all cross-sections are oriented looking downstream. In all figures, the velocity values were normalized according to the corresponding mean inlet velocity magnitude.

As the velocity increases to a maximum during peak systole (Figure 5), blood flow accelerates along the curvature of the arch, with a tendency to skew toward the inner wall of the ascending aorta (α′) as seen in TS1 and TS3, and to a lesser extent in H1 and TS4. In H2 and TS2 the flow was skewed anterolaterally, and in H3 the flow was uniform. In slice β-β′, the flow was skewed laterally or posterolaterally for all healthy aortae, posteriorly for TS1 and TS3, laterally for TS4, and posteromedially for TS2. At the entrance to the descending aorta, flow in the healthy aortae was only slightly skewed: toward the outer curvature wall (γ) in H1 and H3, and toward the inner curvature wall (γ′) in H2. In the TS aortae, flow was more visible skewed, specifically toward the outer curvature wall in TS1 and TS3, and the anterior wall in TS2. In the remainder of the descending aorta (δ-δ′), the flow field was well-structured with little asymmetry, except in TS2 where a region of higher velocity flow was located near the anteromedial wall. For slices ε-ε′ to η-η′, the individual branch anatomy of each aorta influenced the axial and transverse flow fields. When asymmetry of the axial flow was present, it tended to be skewed toward a region of the lateral wall, as seen in the brachiocephalic branch (ε-ε′) of TS2 and TS4, the left common carotid branch (ζ-ζ′) of H3 and TS2, and the left subclavian branch (η-η′) of H1. In the left common carotid branch (ζ-ζ′) of TS3, the flow velocity was greatest posterior to the vessel center.

Figure 5 also revealed further information on the presence of secondary flows at peak velocity that were not visible in Figure 4A. In the ascending aorta (α-α′), even though the flow was laminar for all cases presented (Figure 4A), the direction of secondary flow differed among the patients, although it was never in the direction of the posterior wall. By the time the flow reached the distal end of the aortic arch (β-β′), the secondary flow direction changed for all geometries, except for H3 and TS1, while recirculatory flow was present along the anterior, lateral, or anterolateral walls in TS2–TS4, respectively. Secondary flow in the first descending aortic slice (γ-γ′) of the healthy aortae was present as a single recirculation zone at the posterolateral wall of H3 exclusively. The same region of recirculating flow was present in the first descending aortic slice of TS1. In TS2–TS4, secondary flow in slice γ-γ′ presented as a pair of counter-rotating vortices at the posterior wall of TS2, the posterolateral wall of TS3, and the anterolateral wall of TS4. An additional region of recirculating flow was present at the anterior wall of TS3. In slice δ-δ′, taken further downstream in the descending aorta, secondary flow direction changed from slice γ-γ′ in all geometries except for H2, but remained laminar in all. In the majority of healthy aortic branches, excluding the BCA (ε-ε′) and LSA (η-η′) in H2 and H3 respectively, in-plane velocities were laminar. The same was true for all aortic branches of TS1 and TS4. In TS2, flow recirculation was present in the brachiocephalic branch (ε-ε′) near the anterior wall, and in the left subclavian branch (η-η′) near the anteromedial wall. A single recirculation region was present between the vessel center and medial wall in the LCCA branch (ζ-ζ′), and between the vessel center and anterior wall in the LSA branch (η-η′) of TS3.

Other Hemodynamic Metrics

The distribution of instantaneous normalized wall shear stress (WSSn) at peak systole (t1) and maximum deceleration (t2) are shown, respectively, in Figures 6A,B. Figure 6C displays the normalized time-averaged wall shear stress (TAWSSn) for all models. Additionally, the oscillatory shear index (OSI), normalized relative residence time (RRTn), and transverse WSS (transWSS) are presented in Figures 7A–C. A further analysis on the regional TAWSSn is provided in Figure 8.


Figure 6. (A,B) Instantaneous normalized wall shear stress (WSSn), and (C) normalized time-averaged wall shear stress (TAWSSn) distributions shown (anterior view) for the (H1-H3) healthy, and (TS1-T4) Turner syndrome cases. (A) Peak systole and (B) maximum deceleration. WSS and TAWSS were normalized with respect to the mean WSS at the inlet for each individual case.


Figure 7. (A) Oscillatory shear index (OSI), (B) normalized relative residence time (RRTn), based on the surface integral average for each individual case, and (C) transverse wall shear stress (transWSS) distributions shown (anterior view) for the (H1-H3) healthy and (TS1-T4) Turner syndrome cases.


Figure 8. Normalized time-averaged wall shear stress (TAWSSn) values at seven cross-sections along the aorta of the healthy (black, dotted lines, as average of H1–H3) and Turner syndrome girls (red, solid lines, as average of TS1–TS4). Cross-sections are located as shown on the 3D model of a healthy (H1) case. Standard deviation shown as error bars at each point.

The WSSn distribution differed throughout the cardiac cycle for all models but was most visible at peak systole (Figure 6A) due to the lower values at maximum deceleration (Figure 6B). At peak systole (Figure 6A), each TS aortic geometry presented a highly heterogeneous wall shear stress pattern compared to the less diversified distribution on the aortic wall of the healthy cases. The WSSn pattern at maximum deceleration (Figure 6B) was less clear, with values <20 throughout the majority of the H2, H3, and TS2-TS4 aorta. The peak systolic WSSn had evidently the greatest influence throughout the cardiac cycle, as the TAWSSn (Figure 6C) results were to a large degree similar in character to the WSSn patterns at t1 (Figure 6A). In the healthy (H1-H3) geometries the majority of the aortic body was exposed to relatively low TAWSSn (0–5) values (Figure 6C), with additional higher TAWSSn regions at the inner ascending aortic wall, and the proximal entrance to the descending aorta in H1. In general for H1–H3, along the greater arch curvature TAWSSn was highest upstream of the branch junctions and lowest downstream. Along the lesser arch curvature, TAWSSn was lowest at the arch entrance and highest at the proximal entrance to the descending aorta. TAWSSn on the three healthy branching arteries was lowest on the proximal walls. In the descending aorta, low TAWSSn was observed on the downstream proximal wall and higher TAWSSn on the distal wall. Of the regions identified on the healthy aortae with the lowest TAWSSn (the inner arch curvature wall, the proximal walls of the three branches, downstream of branch junctions, and the downstream proximal wall of the descending aorta), all were accompanied by elevated OSI and RRTn values (Figures 7A,B). The transWSS (Figure 7C), which provides information on the multi-directionality of WSS, was greater in H1 than H2-H3, with the exception of the brachiocephalic branch in H2. This indicates that the direction of flow in H1 is not parallel to a single axis and experiences fluctuations.

The relatively homogenous distribution of TAWSSn on the healthy aortae was replaced with higher TAWSSn values and highly heterogeneous patterns on the walls of the Turner syndrome geometries (Figure 6C). In the dilated aorta of TS1, low TAWSSn values (0–5) were concentrated at the proximal wall of the ascending aorta, downstream of the three branch junctions, and at the proximal wall of the left subclavian artery. Irregular low shear stresses were also found along the proximal wall of the descending aorta among high shear stress regions. Very high TAWSSn values dominated the three aortic branches and almost circumferentially at two locations along the aorta, the first being between the left common carotid and left subclavian branch, and the second at the transition from the arch to the descending aorta. OSI values were highest at the proximal walls of the ascending and descending aorta, and downstream of the three branch junctions, more specifically at the anterior side. As low magnitude and high oscillatory wall shear stress leads to an increase in the residence time of the blood adjacent to the arterial wall, particularly elevated RRTn values were observed at the ascending and descending aorta, and the distal wall of the brachiocephalic branch for TS1 (Figure 7B). The magnitude of transWSS in TS1 was overall very low, with higher values concentrated to the lesser curvature of the proximal arch and the aortic branches (Figure 7C). Case TS2 was defined by moderate-high TAWSSn and transWSS values (Figures 6C, 7C), extending from the distal wall of the ascending aorta and throughout the majority of the elongated arch. The lowest TAWSSn (Figure 6C), accompanied by high OSI (Figure 7A) and RRTn (Figure 7B), and low transWSS (Figure 7C), was found along the proximal wall of the descending aorta. TS3 had a very heterogeneous time-averaged wall shear stress pattern throughout the entire geometry, with the highest values concentrated to the outer curvature of the arch, near branch junctions, and at the transition from the arch to the descending aorta, as seen in TS1 and TS2. As expected, the OSI and RRTn distributions on the surface of TS3 were also highly heterogeneous, with elevated values dominating the descending aorta (Figures 7A,B). TransWSS (Figure 7C) was highest along the anterior wall of the ascending aorta and arch, as was the case for TS2. TAWSSn heterogeneity in TS4 was concentrated to the outer walls of the ascending aorta, arch, and descending aorta, with low TAWSSn and elevated OSI and RRTn transcending from the lateral wall of the ascending aorta, through the lesser curvature of the arch, to the proximal wall of the descending aorta. The transWSS in TS4 was of lower magnitude than that of TS2–TS3, although the highest values were again found along the anterior wall of the ascending aorta and arch.

A detailed, comparative analysis of the TAWSS and TAWSSn was performed at eight circumferential points of the selected aortic slices for the Turner syndrome group, taken as an average, and compared with those of the healthy group. This investigation revealed distinctly higher TAWSSn values at all eight points of all cross-sections for the TS group (Figure 8). The only regions of the aorta that did not fit this trend were at points ε and ε′ of the brachiocephalic trunk. Cross-sections ζ-ζ′, and η-η′ in the left common carotid and left subclavian branches showed the greatest TAWSSn variation between the healthy and TS groups. When comparing the dimensional TAWSS and non-dimensional TAWSSn values averaged over each analysis plane for both groups (Table 4), the TS group had higher average values at every aortic location for both TAWSS and TAWSSn, with the exception of the brachiocephalic branch (ε-ε′) for TAWSS.


Table 4. Time-averaged wall shear stress (TAWSS), expressed in Pascals (Pa), and normalized time-averaged wall shear stress (TAWSSn) averaged over each analysis plane (α-α′ to η-η') for the healthy (n = 3) and the Turner syndrome groups (n = 4).


In our study, patient-specific blood flow simulations were performed through the aortae of four Turner syndrome girls and three healthy girls. In healthy aortae, blood flow patterns range greatly throughout the systolic period (61, 62). The TS girls included in this study had on average larger aortic root diameters (23.66 ± 6.93 mm) and higher peak flow rates (269 ± 49 ml/s, mean ± SD) than the healthy group (19.90 ± 2.48 mm and 241 ± 55 ml/s respectively, mean ± SD). At peak systole, the velocity streamlines revealed the highest flow velocities in the descending aorta of the TS geometries with the largest diameters (TS1 and TS3), but no vortical flow, except at some small recirculation areas in the distal branches of TS3 (Figure 4A). The evolution of helical flow throughout peak to late systole is influenced by the curvature and non-planarity of the arch (63). The visualization of multidirectional blood flow at maximum deceleration (Figure 4B), revealed increased vortical flow in the arch, supra-aortic branches, and descending aorta of the TS patients compared to the healthy controls. This is likely due to the anatomical variants of the TS aortae, rather than the transition from peak systolic flow to diastolic flow, as the peak deceleration value was higher for the healthy group (−13.73 m/s2 ± 1.52) than the TS group (−7.98 ± 1.40 m/s2, mean ± SD).

Atherosclerosis has a non-uniform distribution within the arterial system and shows a predilection for arterial curvatures and branch junctions. However, it is yet unclear why the endothelium is at times prone or immune to disease, particularly with increasing age. One of our key findings in this study was that compared to the relatively homogeneous distribution of time-averaged wall shear stress on the healthy aortae, there was a highly heterogeneous pattern on the walls of the Turner syndrome geometries. Taking into account the young age of our TS cohort (average age 13 ± 2 years) and their obesity level (75% of the TS girls examined here were obese or overweight, Table 1), which is believed to accelerate the initiation and progression of endothelial dysfunction (64), the higher heterogeneity of WSS distribution found in TS girls may be an indication of atherosclerotic disease manifesting earlier in life in these patients.

Wall shear stress is a biomechanical force, predetermined by fluid flow and arterial geometry, and a key player in the pathogenesis of atherosclerosis (65). Aortic sites where blood flow departs from a laminar, unidirectional pattern, and wall shear stresses are low and/or oscillatory, are often associated with locations of atherosclerosis development (27, 63). However, this link has primarily been established in the literature for adults (27, 63). Studies on other species of younger age (e.g., mice and rabbits) do not necessarily follow the low and/or oscillatory WSS theory (22, 37) and the underlying flow mechanisms related to children are not well-known in general. In our study, there appeared to be a correlation between the presence of aortic abnormalities, as seen in the TS group, and overall elevated TAWSSn. Higher TAWSS values and increased vortical flows may indicate a less atheroprotective environment in the young TS patient cohort studied here, thus not showing complete consistency with the above theory. However, the time-averaged WSS that we used here as part of our conclusions may not be the best indicator of disease predisposition, since atherogenesis may also be associated with the time that endothelial cells are exposed to WSS, rather than TAWSS, during the cardiac cycle (37). Additionally, atherogenesis may not depend only on low and/or oscillatory WSS. Mass transfer of molecules such as Low Density Lipoproteins (LDL) or Nitric Oxide (NO) between the blood and the endothelium may also play a role in the mechanism of atherogenesis, even though they might be related to WSS (37).

Comparison With Other Studies

As far as we are aware, there are only three published studies on the aortic flow of TS patients. In the 2014 study by Chen et al. (25), a steady flow rate of 4.5l/min was simulated through three patient-specific Turner syndrome aortae (without patient age information), with aortic abnormalities common to two aortic geometries (TS1 and TS2) in our study. Although our study includes more complex simulations, accounting for patient-specific transient flow, similarities can be made with Chen et al. (25) regarding the flow patterns: weak secondary flow in the ascending aorta (Figure 5 slice α-α') became stronger throughout the transverse aorta (Figure 5 slice β-β'); vortical flow was present in the arch of the aortic model with ETA (TS2) (Figure 5 slice β-β'); and the flow in the descending aorta was remarkably different between TS models (Figure 5 slices γ-γ', δ-δ'). While we report a large variation in the flow patterns of the descending aorta among TS1–TS4, we further found much greater asymmetry in the through-plane velocities and stronger secondary flow, especially at maximum deceleration, than Chen et al. (25). In the same study, the authors determined three locations with low WSS (<0.5 Pa): (1) the brachiocephalic artery, (2) the inner side of the aortic arch exit, and (3) the entire inner wall of the descending aorta, which is in agreement with our results of low TAWSSn values, across all healthy and TS aortae.

In a more recent study by Prahl Wittberg et al. (49), a non-patient-specific flow rate of 5L/min was simulated through four different TS geometries using a non-Newtonian flow model (no patient age information provided). The authors found recirculation zones in all TS geometries, specifically at a cross-section in the descending aorta, which was most prominent (in length and width) in the aorta with ETA. This is in good agreement with Figure 5, where in-plane vectors revealed secondary flow in the descending aorta (slice γ-γ') of all TS geometries, as well as H3. While a prominent region of recirculating flow was seen in TS2 with ETA, this region was of a similar length and width to that seen in TS3 and TS4 (Figure 5). In the arch of a TS geometry with ETA, Prahl Wittberg et al. (49) reported a single recirculation zone, as did our results for TS2 at peak systole (t1) (Figure 5 slice β-β'). At the entrance to the left subclavian branch in the ETA model, Prahl Wittberg et al. (49) visualized recirculating flow, which we also observed anteromedially to the vessel center in TS2 (Figure 5 slice η-η'). The low-velocity recirculatory regions identified by Prahl Wittberg et al. (49) were accompanied by low TAWSS (<0.5 Pa) and high OSI values. Additionally, the authors identified high TAWSS (3 Pa) on the walls of the supra-aortic branches in their normal, dilated, and ETA geometries, and very high TAWSS (5–10 Pa) in the arch of the dilated and ETA models (49). We found similar TAWSS values on the brachiocephalic walls of both the healthy (3.50 ± 0.86 Pa) and TS groups (2.52 ± 0.86 Pa), on the left common carotid walls of the TS group (3.83 ± 0.98 Pa) (Table 4), and on the arch cross-section of TS2 (5.28 ± 1.73 Pa).

The 4D-flow MRI study in the aorta of both healthy and Turner syndrome girls by Arnold et al. (66) found significant differences between the two cohorts (mean age patient group 16 ± 5 years for TS, 17 ± 4 years for healthy controls). Specifically, during early and late systole, and early diastole, helical flow was increased in the ascending and descending regions of TS individuals with significantly larger diameters. Peak systolic velocities did not vary significantly between the TS and control group. In our study, the TS geometries with the largest aortic diameters were TS1 and TS3, for which we demonstrated a similar trend, particularly at mid-diastole (Figure 4C). Additionally, Arnold et al. (66) reported decreased peak systolic WSS throughout the body of the aorta in TS girls, especially at larger diameter regions, when compared to healthy controls. Compared with Figure 6A, we observed a similar pattern on the ascending aortic wall of TS1, but not for TS2-TS4 which had only marginally larger diameters than the healthy controls. Also, our regional TAWSS results (Table 4, Figure 8), averaged over each cross-section, and extended to the supra-aortic branches, concluded that TS girls generally have higher TAWSS and TAWSSn than their healthy counterpart.

Study Limitations

There are some limitations in this study considering the assumptions made in sections Boundary Conditions and Numerical Method. The assumptions of the arterial wall as rigid and blood as Newtonian are reasonable as the deformation of the arterial wall is relatively small in the proximal aorta and the effects of non-Newtonian blood can be ignored in larger vessels (67). In this study, the applied inlet boundary condition is the volumetric waveform and so information on the spatial profile of the inlet velocity is lost. This was because the 3D velocity profile was not available for the healthy cases. The Murray's law, utilized here to predict the outflow percentages, is based on fully developed laminar Poiseuille flow which is not the case for larger arteries where turbulence occurs (68, 69). However, in the case of absent patient data at the model outlets, and a lack of Windkessel parameter values from the literature, especially for children with Turner syndrome, Murray's law was considered the most appropriate method for the outlet boundary conditions. Furthermore, it was difficult to validate the results of this study since the patient data was obtained retrospectively. Prospective MRI data obtained alongside Doppler ultrasound measurements of blood velocity would allow the validation of the numerical results. Nonetheless, our numerical methodology has been used extensively in the past by our group and other researchers for biomedical engineering problems increasing our confidence in the obtained results (20, 21, 37, 70). Finally, due to the limited availability of retrospective MRI data, especially for children, this study examined only a small cohort of TS patients at a specific age range, which undoubtedly limits the statistical significance of our results.

Clinical Significance

Computational fluid dynamics can be a very powerful tool in cardiovascular medicine. The evaluation of hemodynamic parameters of patient-specific aortic models that are challenging to measure in vivo can improve our understanding of cardiovascular disease processes, thus enhancing diagnostic capabilities, and progressing toward patient-specific precision medicine. TS girls and women face a lifelong battle with a broad spectrum of cardiovascular concerns, from congenital heart abnormalities to an increased risk of hypertension, ultimately reducing life expectancy. The current management of cardiovascular conditions in TS is the same as in the general population, due to a lack of understanding of the developmental origins of the cardiovascular manifestations seen specifically in TS (4). The key clinical question is whether any of the cardiovascular risk in TS patients is modifiable. Currently, known modifiable risk factors are hypertension and obesity and therefore the treatment of hypertension and avoidance of being overweight and obese is necessary for this group of patients. This begs the question - would this alter the flow dynamics more favorably too? In this study, the anatomical abnormalities observed in the aortae of young TS girls, the majority (75%) of which were obese or overweight, were accompanied by abnormal flow patterns and highly non-uniform distribution of wall shear stresses, which may promote the development of cardiovascular diseases earlier in life. The type of analysis presented in this study could be used clinically to predict patients at higher risk and therefore be more pro-active in lifestyle measures. Aortic dilatation, although rare, is more common in TS patients and is catastrophic, associated with high mortality. Other risk factors are bicuspid aortic valves, hypertension and previous aortic surgeries, as well as pregnancy. It is still very difficult to predict even if aortic size index, a criterion to estimate the probability for aortic dissection in TS patients, has been evaluated correctly (71). Pregnancy in some recommendations is contraindicated in all of TS due to this risk. CFD modeling could help in this direction (72). Considering the excess of morbidity and mortality, the early diagnosis of cardiovascular changes associated with Turner syndrome is essential, and given the advantages of CFD in monitoring these changes, this method should be used alongside the standard Doppler echocardiography and magnetic resonance imaging in the clinical assessment of these patients.


Personalized flow in seven patient-specific aortic geometries (four TS, three healthy) was investigated in a comprehensive investigation of Turner syndrome children, using computational fluid dynamic methods. The visualization of multidirectional blood flow revealed increased vortical flow in the arch, supra-aortic branches, and descending aorta of the TS girls compared to the healthy controls, but no significant difference was found in the ascending aorta. This study found that young Turner syndrome patients, at about 13 years of age, exhibit markedly elevated TAWSS values when compared to age-and gender-matched healthy controls, suggesting a correlation with the presence of aortic abnormalities. The Turner syndrome girls, 75% of whom were obese or overweight, also display a highly heterogeneous pattern of TAWSS on the aorta, compared to the relatively homogeneous distribution of the healthy aortae, which may be an indication of atherosclerotic disease manifesting earlier in life in these patients. These results are not fully consistent with the low and/or oscillatory WSS theory of atherogenesis, which nonetheless may depend also on other factors, e.g., age, endothelial cells' exposure time to WSS and not TAWSS, and mass transfer of LDL and NO molecules to the endothelium, as explained here. For TS children, age, obesity and aortic morphology may, therefore, play a key role in assessing cardiovascular risk. Further investigations are required on the role of mechanobiological factors on atherogenesis in healthy and diseased children, and the hemodynamics in different age groups of TS patients in particular.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by UK Health Research Authority, London–Westminster Research Ethics Committee (IRAS Project ID: 252866, REC Reference: 18/LO/2052). Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author Contributions

All roles according to CRediT (contributor roles taxonomy). LJ: methodology, software, formal analysis, writing—original draft, and visualization. AK and LJ: conceptualization, investigation, data curation, writing—review and editing, and project administration. AK, AM, RA, and PH: resources and writing—review and editing. AK: supervision and funding acquisition. All authors contributed to the article and approved the submitted version.


This work was supported in part by the UK Research and Innovation (UKRI) Engineering and Physical Sciences Research Council (EPSRC) Award Ref. 2104390 through the University of Strathclyde Research Studentship Scheme (SRSS) Student Excellence Awards (SEA) Project No 1619, the UKRI Natural Environment Research Council (NERC) Award Ref. NE/T014113/1, and the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 749185. Funding for the open access publication fee was provided by UKRI through the above EPSRC and NERC projects.

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.


The authors greatly acknowledge the support from the University of Strathclyde, the Royal Hospital for Children (Glasgow, UK), and the Queen Elizabeth University Hospital (Glasgow, UK). Results were obtained using the ARCHIE-WeSt High Performance Computer ( based at the University of Strathclyde.


TS, Turner syndrome; ETA, Elongated transverse arch; BAV, Bicuspid aortic valve; CoA, Coarctation of the aorta; BCA, Brachiocephalic artery; RSA, Right subclavian artery; RCCA, Right common carotid artery; LCCA, Left common carotid artery; LSA, Left subclavian artery; WSS, Wall shear stress; TAWSS, Time-averaged wall shear stress; transWSS, transverse wall shear stress; OSI, Oscillatory shear index; RRT, Relative residence time.


1. Bondy CA. Clinical practice guideline: care of girls and women with turner syndrome: a guideline of the turner syndrome study group. J Clin Endocrinol Metab. (2007) 92:10–25. doi: 10.1210/jc.2006-1374

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Cockwell A, MacKenzie M, Youings S, Jacobs P. A cytogenetic and molecular study of a series of 45, X fetuses and their parents. J Med Genet. (1991) 28:151–5. doi: 10.1136/jmg.28.3.151

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Nielsen J, Wohlert M. Chromosome abnormalities found among 34910 newborn children: results from a 13-year incidence study in Århus, Denmark. Hum Genet. (1991) 87:81–3. doi: 10.1007/BF01213097

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Silberbach M, Roos-Hesselink JW, Andersen NH, Braverman AC, Brown N, Collins RT, et al. Cardiovascular health in turner syndrome: a scientific statement from the american heart association. Circ Genomic Precis Med. (2018) 11:e000048. doi: 10.1161/HCG.0000000000000048

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Gutmark-Little I, Backeljauw PF. Cardiac magnetic resonance imaging in turner syndrome. Clin Endocrinol. (2013) 78:646–658. doi: 10.1111/cen.12157

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ho VB, Bakalov VK, Cooley M, Van PL, Hood MN, Burklow TR, et al. Major vascular anomalies in turner syndrome: prevalence and magnetic resonance angiographic features. Circulation. (2004) 110:1694–700. doi: 10.1161/01.CIR.0000142290.35842.B0

PubMed Abstract | CrossRef Full Text | Google Scholar

7. de Castro AVB, Okoshi K, Ribeiro SM, Barbosa MF, Mattos PF, Pagliare L, et al. Cardiovascular assessment of patients with ullrich-turner's syndrome on doppler echocardiography and magnetic resonance imaging. Arq Bras Cardiol. (2002) 78:51–8. doi: 10.1590/S0066-782X2002000100005

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kim HK, Gottliebson W, Hor K, Backeljauw P, Gutmark-Little I, Salisbury SR, et al. Cardiovascular anomalies in turner syndrome: spectrum, prevalence, and cardiac mRI findings in a pediatric and young adult population. Am J Roentgenol. (2011) 196:454–60. doi: 10.2214/AJR.10.4973

PubMed Abstract | CrossRef Full Text | Google Scholar

9. De Groote K, Demulier L, De Backer J, De Wolf D, De Schepper J, T'Sjoen G, et al. Arterial hypertension in turner syndrome: a reviewof the literature and a practical approach for diagnosis and treatment. J Hypertens. (2015) 33:1342–51. doi: 10.1097/HJH.0000000000000599

CrossRef Full Text | Google Scholar

10. von Rindfleisch GE, Kloman WC, Miles FT. A Text-Book of Pathological Histology. An Introduction to the Study of Pathological Anatomy. London: Lindsay and Blakiston (1872). doi: 10.1097/00000441-187207000-00017

CrossRef Full Text | Google Scholar

11. Fry DL. Certain chemorheologic considerations regarding the blood vascular interface with particular reference to coronary artery disease. Circulation. (1969) 40:IV-38–IV-57. doi: 10.1161/01.CIR.40.5S4.IV-38

CrossRef Full Text | Google Scholar

12. Caro CG, Fitz-Gerald JM, Schroter RC. Atheroma and arterial wall shear. Observation, correlation and proposal of a shear dependent mass transfer mechanism for atherogenesis. Proc R Soc London Ser B Biol Sci. (1971) 177:109–59. doi: 10.1098/rspb.1971.0019

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Caro CG, Pedley TJ, Schroter RC, Seed WA, Parker KH. The Mechanics of the Circulation, 2nd Edn. Cambridge: Cambridge University Press (2011).

Google Scholar

14. Ku DN, Giddens DP, Zarins CK, Glagov S. Pulsatile flow and atherosclerosis in the human carotid bifurcation. Positive correlation between plaque location and low and oscillating shear stress. Arteriosclerosis. (1985) 5:293–302. doi: 10.1161/01.ATV.5.3.293

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Weinberg PD. Disease patterns at arterial branches and their relation to flow. Biorheology. (2002) 39:533–537. Available online at: (accessed February 20, 2021).

PubMed Abstract | Google Scholar

16. Sinzinger H, Silberbauer K, Auerswald W. Paradoxical relaxation of arterial strips induced by vasoconstrictive agents. J Vasc Res. (1980) 17:44–52. doi: 10.1159/000158233

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Sloop GD, Perret RS, Brahney JS, Oalmann M. A description of two morphologic patterns of aortic fatty streaks, and a hypothesis of their pathogenesis. Atherosclerosis. (1998) 141:153–60. doi: 10.1016/S0021-9150(98)00167-1

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Svindland A, Walløe L. Distribution pattern of sudanophilic plaques in the descending thoracic and proximal abdominal human aorta. Atherosclerosis. (1985) 57:219–24. doi: 10.1016/0021-9150(85)90035-8

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Simms MH, Welfe JHN, Campbell WB. Arterial Disease. Oxford: Blackwell Scientific Publications (1988). doi: 10.1002/bjs.1800750545

CrossRef Full Text | Google Scholar

20. Kazakidi A, Sherwin SJ, Weinberg PD. Effect of reynolds number and flow division on patterns of haemodynamic wall shear stress near branch points in the descending thoracic aorta. J R Soc Interface. (2009) 6:539–48. doi: 10.1098/rsif.2008.0323

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Kazakidi A, Plata AM, Sherwin SJ, Weinberg PD. Effect of reverse flow on the pattern of wall shear stress near arterial branches. J R Soc Interface. (2011) 8:1594–603. doi: 10.1098/rsif.2011.0108

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Vincent PE, Plata AM, Hunt AAE, Weinberg PD, Sherwin SJ. Blood flow in the rabbit aortic arch and descending thoracic aorta. J R Soc Interface. (2011) 8:1708–19. doi: 10.1098/rsif.2011.0116

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Morris PD, Narracott A, Von Tengg-Kobligk H, Soto DAS, Hsiao S, Lungu A, et al. Computational fluid dynamics modelling in cardiovascular medicine. Heart. (2016) 102:18–28. doi: 10.1136/heartjnl-2015-308044

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Antiga L. Patient-Specific Modeling of Geometry and Blood Flow in Large Arteries. Milan: Politec di Milano. (2002).

Google Scholar

25. Chen J, Gutmark E, Mylavarapu G, Backeljauw PF, Gutmark-Little I. Numerical investigation of mass transport through patient-specific deformed aortae. J Biomech. (2014) 47:544–52. doi: 10.1016/j.jbiomech.2013.10.031

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Liu X, Pu F, Fan Y, Deng X, Li D, Li S. A numerical study on the flow of blood and the transport of lDL in the human aorta: the physiological significance of the helical flow in the aortic arch. Am J Physiol Hear Circ Physiol. (2009) 297:H163–70. doi: 10.1152/ajpheart.00266.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Numata S, Itatani K, Kanda K, Doi K, Yamazaki S, Morimoto K, et al. Blood flow analysis of the aortic arch using computational fluid dynamics. Eur J Cardio-thoracic Surg. (2016) 49:1578–85. doi: 10.1093/ejcts/ezv459

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Pirola S, Jarral OA, O'Regan DP, Asimakopoulos G, Anderson JR, Pepper JR, et al. Computational study of aortic hemodynamics for patients with an abnormal aortic valve: the importance of secondary flow at the ascending aorta inlet. APL Bioeng. (2018) 2:026101. doi: 10.1063/1.5011960

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Sundström E, Jonnagiri R, Gutmark-Little I, Gutmark E, Critser P, Taylor MD, et al. Effects of normal variation in the rotational position of the aortic root on hemodynamics and tissue biomechanics of the thoracic aorta. Cardiovasc Eng Technol. (2020) 11:47–58. doi: 10.1007/s13239-019-00441-2

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Youssefi P, Gomez A, Arthurs C, Sharma R, Jahangiri M, Figueroa CA. Impact of patient-specific inflow velocity profile on hemodynamics of the thoracic aorta. J Biomech Eng. (2018) 140:1–14. doi: 10.1115/1.4037857

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zhu Y, Chen R, Juan YH, Li H, Wang J, Yu Z, et al. Clinical validation and assessment of aortic hemodynamics using computational fluid dynamics simulations from computed tomography angiography. Biomed Eng Online. (2018) 17:53. doi: 10.1186/s12938-018-0485-5

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Madhavan S, Kemmerling EMC. The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow. Biomed Eng Online. (2018) 17:66. doi: 10.1186/s12938-018-0497-1

PubMed Abstract | CrossRef Full Text | Google Scholar

33. McElroy M, Keshmiri A. Impact of using conventional inlet/Outlet boundary conditions on haemodynamic metrics in a subject-Specific rabbit aorta. Proc Inst Mech Eng Part H J Eng Med. (2018) 232:103–13. doi: 10.1177/0954411917699237

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Gallo D, De Santis G, Negri F, Tresoldi D, Ponzini R, Massai D, et al. On the use of in vivo measured flow rates as boundary conditions for image-based hemodynamic models of the human aorta: implications for indicators of abnormal flow. Ann Biomed Eng. (2012) 40:729–41. doi: 10.1007/s10439-011-0431-1

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Morbiducci U, Gallo D, Massai D, Consolo F, Ponzini R, Antiga L, et al. Outflow conditions for image-based hemodynamic models of the carotid bifurcation: implications for indicators of abnormal flow. J Biomech Eng. (2010) 132:091005. doi: 10.1115/1.4001886

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Morbiducci U, Ponzini R, Gallo D, Bignardi C, Rizzo G. Inflow boundary conditions for image-based computational hemodynamics: impact of idealized versus measured velocity profiles in the human aorta. J Biomech. (2013) 46:102–9. doi: 10.1016/j.jbiomech.2012.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Van Doormaal MA, Kazakidi A, Wylezinska M, Hunt A, Tremoleda JL, Protti A, et al. Haemodynamics in the mouse aortic arch computed from mRI-derived velocities at the aortic root. J R Soc Interface. (2012) 9:2834–44. doi: 10.1098/rsif.2012.0295

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Pirola S, Cheng Z, Jarral OA, O'Regan DP, Pepper JR, Athanasiou T, et al. On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics. J Biomech. (2017) 60:15–21. doi: 10.1016/j.jbiomech.2017.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Chnafa C, Brina O, Pereira VM, Steinman DA. Better than nothing: a rational approach for minimizing the impact of outflow strategy on cerebrovascular simulations. Am J Neuroradiol. (2018) 39:337–43. doi: 10.3174/ajnr.A5484

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Murray CD. The physiological principle of minimum work: I. The vascular system and the cost of blood volume. Proc Natl Acad Sci USA. (1926) 12:207–14. doi: 10.1073/pnas.12.3.207

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Freeman JV, Cole TJ, Chinn S, Jones PRM, White EM, Preece MA. Cross sectional stature and weight reference curves for the UK, 1990. Arch Dis Child. (1995) 73:17–24. doi: 10.1136/adc.73.1.17

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Kaiser T, Kellenberger CJ, Albisetti M, Bergsträsser E, Buechel ERV. Normal values for aortic diameters in children and adolescents - assessment in vivo by contrast-enhanced cMR-angiography. J Cardiovasc Magn Reson. (2008) 10:56. doi: 10.1186/1532-429X-10-56

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Ostberg JE, Donald AE, Halcox JPJ, Storry C, McCarthy C, Conway GS. Vasculopathy in turner syndrome: arterial dilatation and intimal thickening without endothelial dysfunction. J Clin Endocrinol Metab. (2005) 90:5161–6. doi: 10.1210/jc.2005-0677

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Febina J, Sikkandar MY, Sudharsan NM. Wall shear stress estimation of thoracic aortic aneurysm using computational fluid dynamics. Comput Math Methods Med. (2018) 2018:7126532. doi: 10.1155/2018/7126532

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Sosnowski M, Krzywanski J, Gnatowska R. Polyhedral meshing as an innovative approach to computational domain discretization of a cyclone in a fluidized bed CLC unit. E3S Web Conf. (2017) 14:01027. doi: 10.1051/e3sconf/20171401027

CrossRef Full Text | Google Scholar

46. Celik I, Ghia U, Roache PJ, Freitas CJ. Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. J Fluids Eng. (2008) 130:078001. doi: 10.1115/1.2960953

CrossRef Full Text | Google Scholar

47. Lewandowska N, Mosiȩżny J. Meshing strategy for bifurcation arteries in the context of blood flow simulation accuracy. E3S Web Conf. (2019) 128:02003. doi: 10.1051/e3sconf/201912802003

CrossRef Full Text | Google Scholar

48. Churchill SW. Friction factor equation spans all fluid flow regimes. Chem Eng J. (1977) 84:91–92. Available online at: (accessed February 20, 2021).

Google Scholar

49. Prahl Wittberg L, van Wyk S, Fuchs L, Gutmark E, Backeljauw P, Gutmark-Little I. Effects of aortic irregularities on blood flow. Biomech Model Mechanobiol. (2016) 15:345–60. doi: 10.1007/s10237-015-0692-y

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Ponzini R, Vergara C, Rizzo G, Veneziani A, Roghi A, Vanzulli A, et al. Womersley number-based estimates of blood flow rate in doppler analysis: in vivo validation by means of phase-contrast MRI. IEEE Trans Biomed Eng. (2010) 57:1807–15. doi: 10.1109/TBME.2010.2046484

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Loudon C, Tordesillas A. The use of the dimensionless womersley number to characterize the unsteady nature of internal flow. J Theor Biol. (1998) 191:63–78. doi: 10.1006/jtbi.1997.0564

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Cheng C, Helderman F, Tempel D, Segers D, Hierck B, Poelmann R, et al. Large variations in absolute wall shear stress levels within one species and between species. Atherosclerosis. (2007) 195:225–35. doi: 10.1016/j.atherosclerosis.2006.11.019

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Reneman RS, Hoeks APG. Wall shear stress as measured in vivo: consequences for the design of the arterial system. Med Biol Eng Comput. (2008) 46:499–507. doi: 10.1007/s11517-008-0330-2

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Zamir M, Sinclair P, Wonnacott TH. Relation between diameter and flow in major branches of the arch of the aorta. J Biomech. (1992) 25:1303–10. doi: 10.1016/0021-9290(92)90285-9

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Lantz J, Karlsson M. Large eddy simulation of lDL surface concentration in a subject specific human aorta. J Biomech. (2012) 45:537–42. doi: 10.1016/j.jbiomech.2011.11.039

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Chnafa C, Mendez S, Nicoud F. Image-based simulations show important flow fluctuations in a normal left ventricle: what could be the implications? Ann Biomed Eng. (2016) 44:3346–58. doi: 10.1007/s10439-016-1614-6

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Khan MO, Valen-Sendstad K, Steinman DA. Narrowing the expertise gap for predicting intracranial aneurysm hemodynamics: impact of solver numerics versus mesh and time-step resolution. Am J Neuroradiol. (2015) 36:1310–6. doi: 10.3174/ajnr.A4263

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Chen Z, Yu H, Shi Y, Zhu M, Wang Y, Hu X, et al. Vascular remodelling relates to an elevated oscillatory shear index and relative residence time in spontaneously hypertensive rats. Sci Rep. (2017) 7:2007. doi: 10.1038/s41598-017-01906-x

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Riccardello GJ, Shastri DN, Changa AR, Thomas KG, Roman M, Prestigiacomo CJ, et al. Influence of relative residence time on side-wall aneurysm inception. Clin Neurosurg. (2018) 83:574–81. doi: 10.1093/neuros/nyx433

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Peiffer V, Sherwin SJ, Weinberg PD. Computation in the rabbit aorta of a new metric – the transverse wall shear stress – to quantify the multidirectional character of disturbed blood flow. J Biomech. (2013) 46:2651–8. doi: 10.1016/j.jbiomech.2013.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Kilner PJ, Yang GZ, Mohiaddin RH, Firmin DN, Longmore DB. Helical and retrograde secondary flow patterns in the aortic arch studied by three-directional magnetic resonance velocity mapping. Circulation. (1993) 88:2235–47. doi: 10.1161/01.CIR.88.5.2235

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Laín S, Caballero AD. Simulación transitoria de la dinámica del flujo sanguíneo en la aorta torácica. Ing Investig. (2017) 37:92–101.

Google Scholar

63. Morbiducci U, Ponzini R, Rizzo G, Cadioli M, Esposito A, Montevecchi FM, et al. Mechanistic insight into the physiological relevance of helical blood flow in the human aorta: an in vivo study. Biomech Model Mechanobiol. (2011) 10:339–55. doi: 10.1007/s10237-010-0238-2

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Berenson GS, Srinivasan SR, Nicklas TA. Atherosclerosis: a nutritional disease of childhood. Am J Cardiol. (1988) 82:22T−29T. doi: 10.1016/s0002-9149(98)00719-x

CrossRef Full Text | Google Scholar

65. Cunningham KS, Gotlieb AI. The role of shear stress in the pathogenesis of atherosclerosis. Lab Investig. (2005) 85:9–23. doi: 10.1038/labinvest.3700215

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Arnold R, Neu M, Hirtler D, Gimpel C, Markl M, Geiger J. Magnetic resonance imaging 4-D flow-based analysis of aortic hemodynamics in turner syndrome. Pediatr Radiol. (2017) 47:382–90. doi: 10.1007/s00247-016-3767-8

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Wittek A, Karatolios K, Fritzen CP, Bereiter-Hahn J, Schieffer B, Moosdorf R, et al. Cyclic three-dimensional wall motion of the human ascending and abdominal aorta characterized by time-resolved three-dimensional ultrasound speckle tracking. Biomech Model Mechanobiol. (2016) 15:1375–88. doi: 10.1007/s10237-016-0769-2

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Moreau B, Mauroy B. Murray's law revisited: quémada's fluid model and fractal trees. J Rheol. (2015) 59:1419–30. doi: 10.1122/1.4934240

CrossRef Full Text | Google Scholar

69. Uylings HBM. Optimization of diameters and bifurcation angles in lung and vascular tree structures. Bull Math Biol. (1977) 39:509–20. doi: 10.1016/S0092-8240(77)80054-2

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Boumpouli M, Danton MHD, Gourlay T, Kazakidi A. Blood flow simulations in the pulmonary bifurcation in relation to adult patients with repaired tetralogy of fallot. Med Eng Phys. (2020) 85:123–38. doi: 10.1016/j.medengphy.2020.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Nijs J, Gelsomino S, Lucà F, Parise O, Maessen JG, La Meir M. Unreliability of aortic size index to predict risk of aortic dissection in a patient with turner syndrome. World J Cardiol. (2014) 6:349. doi: 10.4330/wjc.v6.i5.349

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Johnston L, Allen R, Hall-Barrientos P, Mason A, Kazakidi A. Computational haemodynamics in turner syndrome patient-specific aortae with PC-MRI obtained boundary conditions. Heart J. (2020) 106:A10. doi: 10.1136/heartjnl-2020-SCF.23

CrossRef Full Text | Google Scholar

Keywords: Turner syndrome, computational fluid dynamics, hemodynamics, cardiovascular disease, pediatric medicine, patient-specific, disturbed flow, atherosclerosis

Citation: Johnston L, Allen R, Hall Barrientos P, Mason A and Kazakidi A (2021) Hemodynamic Abnormalities in the Aorta of Turner Syndrome Girls. Front. Cardiovasc. Med. 8:670841. doi: 10.3389/fcvm.2021.670841

Received: 22 February 2021; Accepted: 06 May 2021;
Published: 01 June 2021.

Edited by:

Umberto Morbiducci, Politecnico di Torino, Italy

Reviewed by:

Michele Conti, University of Pavia, Italy
Bogdan Ene-Iordache, Istituto di Ricerche Farmacologiche Mario Negri (IRCCS), Italy

Copyright © 2021 Johnston, Allen, Hall Barrientos, Mason and Kazakidi. 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: Asimina Kazakidi,