Abstract

OBJECTIVES

Aortic stiffness is an independent cardiovascular risk factor. Cardiac magnetic resonance (CMR) allows evaluation of aortic elastic properties by different indexes such as distensibility, the maximum rate of systolic distension (MRSD) and pulse wave velocity (PWV). We sought to define age-dependent changes of indexes of elastic properties of the thoracic aorta in healthy subjects.

METHODS

We enrolled 85 healthy subjects (53 males) free of overt cardiovascular disease subdivided into 6 classes of age (from 15 to >60 years). Distensibility, MRSD and PWV were measured by the analysis of CMR images acquired using a 1.5 T clinical scanner.

RESULTS

MRSD and distensibility decreased progressively through the classes of age (P < 0.001) after an initial plateau between 20 and 30 years in males and 15 and 20 years in females. Pulse wave velocity increased progressively with the age (P < 0.001). Distensibility was related to body mass index (P = 0.002), surface area (P < 0.005), weight (P = 0.005) and to left ventricular parameters such as mass index (P < 0.001) and end-diastolic volume index (P = 0.002). MRSD was related to end-diastolic volume index (P < 0.001) but not to body parameters. PWV was not related to body and ventricular parameters.

CONCLUSIONS

This study confirmed that physiological ageing is associated with a progressive impairment of the elastic properties of the aortic wall. Results of this study may be useful for the early identification of subjects with impaired aortic wall properties providing referral values of elasticity indexes assessed by CMR in different classes of age.

INTRODUCTION

Arterial stiffness is an independent cardiovascular risk factor [13] and there is a strong relation between ageing and stiffening of the proximal aorta, one of the earliest manifestations of vascular ageing in healthy humans [4].

Functional assessment of thoracic aortic properties such as aortic stiffness by cardiac magnetic resonance (CMR) has been proposed and validated [46]. In previous studies, the evaluation of aortic elastic properties with CMR was performed in patients with Marfan disease [710] and in subjects with bicuspid aortic valve [11, 12]. The conventional parameters for the assessment of the aortic stiffness were regional distensibility and flow wave velocity [13]. Currently, aortic stiffness and flow wave velocity are under investigation as parameters for the indication for ascending aorta replacement. Recently, the maximum rate of systolic distension (MRSD), a new index for the evaluation of aortic wall properties, was demonstrated to be more accurate than distensibility and flow wave velocity to detect initial changes of aortic elastic properties in young subjects with bicuspid aortic valve [12].

The evaluation of the physiological changes of aortic wall properties in healthy subjects of different ages may be relevant in clinical practice. The aim of the present study was to investigate age-related changes of aortic elastic properties by CMR in a population of patients free of overt cardiovascular disease.

METHODS

Healthy subjects were enrolled using the following inclusion criteria: (i) age ≥15 years; (ii) absence of contraindications to CMR; (iii) no history of cardiovascular disease; (iv) no family history of aortic disease, including Marfan syndrome and connective tissue disease; (v) no family history of sudden death; (vi) normal physical examination; (vii) normal echocardiographic examination, with evidence of tricuspid aortic valve and (viii) absence of acute or chronic disease, including diabetes, hypertension (defined as diastolic blood pressure >90 mmHg, systolic blood pressure >140 mmHg or receiving treatment for hypertension), hyperlipidaemia (abnormal lipid levels or receiving lipid-lowering therapy) and obesity (defined as body mass index [BMI] ≥30). Finally, subjects with referred smoke habit or drug assumption were excluded. All the subjects were informed about the study protocol and provided written consent. The procedures followed were in accordance with institutional guidelines and the Declaration of Helsinki. Subjects were subdivided into six class of age: 15–20 years; >20–30 years; >30–40 years; >40–50 years; >50–60 years; >60 years. The enrolment was completed when the CMR protocol was concluded by at least 12 subjects per class of age. The final population included 97 subjects (58 males).

Cardiac magnetic resonance: acquisition protocol

CMR was performed with a dedicated 1.5 Tesla scanner (Signa Hdx, General Electrics Healthcare, Milwaukee, WI, USA) with an eight-channel cardiac phased array coil. Left ventricular (LV) volumes and mass were obtained by acquiring short-axis cine images using conventional approach.

The thoracic aorta was visualized by acquiring sagittal-oblique cine images parallel to the major aortic axis using a breath-hold, electrocardiography-triggered, steady-state free-precession (SSFP) pulse sequence with the following parameters: 400-mm field of view, 8-mm slice thickness, no gap, 1 NEX, 12 views per segment, TE/TR 1.6/3.2, flip angle 45°, matrix 224 × 224 and reconstruction matrix 256 × 256. The number of cardiac phases was set according to the heart rate to obtain an aortic wall excursion temporal resolution of ∼10–3 s. Cross-sectional cine SSFP images with the same parameters were acquired at different aortic levels: (i) at the aortic valve plane to confirm the absence of aortic valve disease and the presence of tricuspid aortic valve; (ii) at the proximal ascending aorta (5 mm above the sinotubular junction); (iii) in section of the descending aorta included in the same images of the proximal ascending aorta and (iv) at the level of the maximum diameter of the ascending aorta.

The central arterial pressure was non-invasively estimated as previously described [14].

Magnetic resonance imaging: acquisition protocol

The SSFP images for the assessment of ventricular function and wall excursion were transferred to a workstation (Advantage Windows 4.1 GE) and elaborated with commercially available software (Mass Plus and CV Flow, MR Analytical Software Systems, Version 4.0, Leiden, Netherlands). Ventricular volumes were indexed for body surface area (BSA). The diameters of the aortic root, sinotubular junction and ascending aorta (at the maximal dimension level) were measured from the cine SSFP images.

The maximal rate of systolic distension (MRSD) was recently described as a parameter useful in the evaluation of the elastic properties of the aortic wall using cine-magnetic resonance imaging [12]. MRSD was measured as follows: The cross-sectional area of the proximal ascending aorta (5 mm above the sinotubular junction) was measured in each cardiac phase, using a previously validated software (Tracking Tool) [15] and if necessary manually corrected. The cross-sectional area of the proximal ascending aorta was indexed for the maximal end systolic cross-sectional area and plotted against time (relative cross-sectional area/time curve, Fig. 1). In this curve, the MRSD of the ascending aorta cross-sectional area was measured as the maximal systolic upslope and expressed as percentile of the maximal area/10–3 s.
Males: Maximum rate of systolic distension (MRSD) of ascending and descending and distensibility of ascending and descending aortas by age decile (mean, 95% prediction). MRSD is expressed as percentage of maximal area/10–3 s. Distensibility is expressed as 10–3 mmHg–1.
Figure 1:

Males: Maximum rate of systolic distension (MRSD) of ascending and descending and distensibility of ascending and descending aortas by age decile (mean, 95% prediction). MRSD is expressed as percentage of maximal area/10–3 s. Distensibility is expressed as 10–3 mmHg–1.

Aortic wall distensibility in the ascending aorta was measured in SSFP images, using a method described in a previous study [16]. Briefly the relative change in luminal aortic area (in square millimetres) was measured as (Area max − Area min)/Area min. Maximal and minimal areas were measured automatically as above. Then, distensibility of ascending aorta in each subject was defined as follows: distensibility = relative change in luminal aortic area/cPP, where cPP is the estimated central pulse pressure. Distensibility was expressed as 10–3 mmHg–1.

The cross-sectional area/time curves generated for the ascending and descending aortas were used to measure the pulse wave velocity (PWV). PWV was defined as the time delay to reach the systolic half-peak of the cross-sectional area/time curve in ascending to descending aorta, divided by the length of the aortic arch between the two sections.

Statistical analysis

The Kolmogorov–Smirnov test was employed to assess normality of data distribution. A P-value of <0.05 was considered statistically significant. The values of clinical variables evaluated in the current study are expressed as the mean ± standard deviation (for normal distribution) or median and 25th–75th (interquartile range [IQR]) percentile. The group means of numerical data were analysed using a two-sample Student’s t-test. Differences among the groups were analysed using the Anova F-test when appropriate. The coefficient of determination for linear regression R2 was used.

RESULTS

Study population

Among the 97 subjects enrolled, 12 were excluded because of claustrophobia (n = 7) and poor image quality (n = 5). Then, the final population included 85 subjects (53 males, mean age 40, IQR percentile 25–52).

The subjects were distributed over classes of age as follows: 12 in ≥15–20 years; 12 in >20–30 years; 14 in >30–40 years; 13 in >40–50 years; 14 in >50–60 years; 12 in >60 years.

Patient general characteristics are shown in Tables 1 and 2

Table 1:

Population characteristics

MaleFemaleP-value
n5332
Age (years)36 ± 18.9540 ± 17.90.34
Height (cm)174.06 ± 10.49163.96 ± 7.9<0.001
Weight (kg)71.93 ± 16.0864.21 ± 140.021
BMI (kg/m2)23.81 ± 3.4723.89 ± 5.220.93
BSA1.87 ± 0.241.7 ± 0.190.001
LV mass index (g/m2)77.28 ± 15.6357.2 ± 11.78<0.001
Aortic arch length (mm)128.55 ± 37.48124.43 ± 21.770.57
SAP (mmHg)120 ± 17120 ± 151
DAP (mmHg)74 ± 1176 ± 90.38
MAP (mmHg)90 ± 1191 ± 100.67
Ascending aorta
 Maximum diameter (mm)29.2 ± 7.0225.51 ± 6.280.016
 Minimum diameter (mm)27.11 ± 7.122.16 ± 5.220.001
 Distensibility0.52 (0.27–1.06)0.59 (0.38–0.94)
 MRSD8.64 ± 3.828.02 ± 3.720.46
 PWV3.96 ± 2.43.84 ± 1.480.83
Descending aorta
 Maximum diameter (mm)24.13 ± 6.319.3 ± 6.520.001
 Minimum diameter (mm)22.08 ± 6.1817.53 ± 5.960.001
 Distensibility0.39 (0.28–0.57)0.47 (0.37–0.62)
 MRSD7.5 ± 3.777.5 ± 3.211
MaleFemaleP-value
n5332
Age (years)36 ± 18.9540 ± 17.90.34
Height (cm)174.06 ± 10.49163.96 ± 7.9<0.001
Weight (kg)71.93 ± 16.0864.21 ± 140.021
BMI (kg/m2)23.81 ± 3.4723.89 ± 5.220.93
BSA1.87 ± 0.241.7 ± 0.190.001
LV mass index (g/m2)77.28 ± 15.6357.2 ± 11.78<0.001
Aortic arch length (mm)128.55 ± 37.48124.43 ± 21.770.57
SAP (mmHg)120 ± 17120 ± 151
DAP (mmHg)74 ± 1176 ± 90.38
MAP (mmHg)90 ± 1191 ± 100.67
Ascending aorta
 Maximum diameter (mm)29.2 ± 7.0225.51 ± 6.280.016
 Minimum diameter (mm)27.11 ± 7.122.16 ± 5.220.001
 Distensibility0.52 (0.27–1.06)0.59 (0.38–0.94)
 MRSD8.64 ± 3.828.02 ± 3.720.46
 PWV3.96 ± 2.43.84 ± 1.480.83
Descending aorta
 Maximum diameter (mm)24.13 ± 6.319.3 ± 6.520.001
 Minimum diameter (mm)22.08 ± 6.1817.53 ± 5.960.001
 Distensibility0.39 (0.28–0.57)0.47 (0.37–0.62)
 MRSD7.5 ± 3.777.5 ± 3.211

Data are expressed as: mean ± standard deviation or median (IQR, interquartile range).

BMI: body mass index; BSA: body surface area; LV: left ventricular; SAP: systolic aortic pressure; DAP: diastolic aortic pressure; MAP: mean aortic pressure; MRSD: maximal rate of systolic distension; PWV: pulse wave velocity.

Table 1:

Population characteristics

MaleFemaleP-value
n5332
Age (years)36 ± 18.9540 ± 17.90.34
Height (cm)174.06 ± 10.49163.96 ± 7.9<0.001
Weight (kg)71.93 ± 16.0864.21 ± 140.021
BMI (kg/m2)23.81 ± 3.4723.89 ± 5.220.93
BSA1.87 ± 0.241.7 ± 0.190.001
LV mass index (g/m2)77.28 ± 15.6357.2 ± 11.78<0.001
Aortic arch length (mm)128.55 ± 37.48124.43 ± 21.770.57
SAP (mmHg)120 ± 17120 ± 151
DAP (mmHg)74 ± 1176 ± 90.38
MAP (mmHg)90 ± 1191 ± 100.67
Ascending aorta
 Maximum diameter (mm)29.2 ± 7.0225.51 ± 6.280.016
 Minimum diameter (mm)27.11 ± 7.122.16 ± 5.220.001
 Distensibility0.52 (0.27–1.06)0.59 (0.38–0.94)
 MRSD8.64 ± 3.828.02 ± 3.720.46
 PWV3.96 ± 2.43.84 ± 1.480.83
Descending aorta
 Maximum diameter (mm)24.13 ± 6.319.3 ± 6.520.001
 Minimum diameter (mm)22.08 ± 6.1817.53 ± 5.960.001
 Distensibility0.39 (0.28–0.57)0.47 (0.37–0.62)
 MRSD7.5 ± 3.777.5 ± 3.211
MaleFemaleP-value
n5332
Age (years)36 ± 18.9540 ± 17.90.34
Height (cm)174.06 ± 10.49163.96 ± 7.9<0.001
Weight (kg)71.93 ± 16.0864.21 ± 140.021
BMI (kg/m2)23.81 ± 3.4723.89 ± 5.220.93
BSA1.87 ± 0.241.7 ± 0.190.001
LV mass index (g/m2)77.28 ± 15.6357.2 ± 11.78<0.001
Aortic arch length (mm)128.55 ± 37.48124.43 ± 21.770.57
SAP (mmHg)120 ± 17120 ± 151
DAP (mmHg)74 ± 1176 ± 90.38
MAP (mmHg)90 ± 1191 ± 100.67
Ascending aorta
 Maximum diameter (mm)29.2 ± 7.0225.51 ± 6.280.016
 Minimum diameter (mm)27.11 ± 7.122.16 ± 5.220.001
 Distensibility0.52 (0.27–1.06)0.59 (0.38–0.94)
 MRSD8.64 ± 3.828.02 ± 3.720.46
 PWV3.96 ± 2.43.84 ± 1.480.83
Descending aorta
 Maximum diameter (mm)24.13 ± 6.319.3 ± 6.520.001
 Minimum diameter (mm)22.08 ± 6.1817.53 ± 5.960.001
 Distensibility0.39 (0.28–0.57)0.47 (0.37–0.62)
 MRSD7.5 ± 3.777.5 ± 3.211

Data are expressed as: mean ± standard deviation or median (IQR, interquartile range).

BMI: body mass index; BSA: body surface area; LV: left ventricular; SAP: systolic aortic pressure; DAP: diastolic aortic pressure; MAP: mean aortic pressure; MRSD: maximal rate of systolic distension; PWV: pulse wave velocity.

Table 2:

Baseline characteristics of normal subjects studied (mean ± standard deviation)

Age class15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Males
 Height (cm)174 ± 7.88181.42 ± 7.04175 ± 14.17175.12 ± 13.87174.62 ± 6.04172.16 ± 5.60
 Weight (kg)68.62 ± 9.5080.28 ± 8.2673.14 ± 15.8674.15 ± 21.1575.12 ± 9.0479.33 ± 12.07
 BSA (m2)1.82 ± 0.142 ± 0.131.88 ± 0.271.88 ± 0.361.9 ± 0.111.94 ± 0.17
 BMI (kg/m2)22.63 ± 3.2024.39 ± 2.1923.6 ± 3.0023.51 ± 4.3524.72 ± 3.4026.67 ± 2.86
 EDVi (ml/m2)86 ± 1089 ± 1186 ± 984 ± 881 ± 1177 ± 9
 LV mass index (ml/m2)76 ± 978 ± 880 ± 1078 ± 976 ± 873 ± 11
 Ejection fraction (%)65 ± 666 ± 561 ± 460 ± 661 ± 759 ± 10
Female
 Height (cm)169.50 ± 3.53169.75 ± 10.96159.60 ± 6.38167.80 ± 6.14163.66 ± 3.93162.20 ± 5.11
 Weight (kg)54.00 ± 2.8270.00 ± 16.1460.20 ± 8.2264.20 ± 13.0261.33 ± 5.4277.20 ± 17.16
 BSA (m2)1.59 ± 0.031.81 ± 0.211.63 ± 0.131.72 ± 0.201.67 ± 0.081.86 ± 0.19
 BMI (kg/m2)18.83 ± 1.7724.59 ± 6.9123.64 ± 2.9222.63 ± 3.1222.90 ± 1.9229.53 ± 7.4
 EDVi (ml/m2)80 ± 1182 ± 1079 ± 1180 ± 876 ± 1173 ± 12
 LV mass index (ml/m2)54 ± 1054 ± 1255 ± 956 ± 859 ± 1260 ± 11
 Ejection fraction (%)61 ± 663 ± 564 ± 663 ± 763 ± 962 ± 11
Age class15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Males
 Height (cm)174 ± 7.88181.42 ± 7.04175 ± 14.17175.12 ± 13.87174.62 ± 6.04172.16 ± 5.60
 Weight (kg)68.62 ± 9.5080.28 ± 8.2673.14 ± 15.8674.15 ± 21.1575.12 ± 9.0479.33 ± 12.07
 BSA (m2)1.82 ± 0.142 ± 0.131.88 ± 0.271.88 ± 0.361.9 ± 0.111.94 ± 0.17
 BMI (kg/m2)22.63 ± 3.2024.39 ± 2.1923.6 ± 3.0023.51 ± 4.3524.72 ± 3.4026.67 ± 2.86
 EDVi (ml/m2)86 ± 1089 ± 1186 ± 984 ± 881 ± 1177 ± 9
 LV mass index (ml/m2)76 ± 978 ± 880 ± 1078 ± 976 ± 873 ± 11
 Ejection fraction (%)65 ± 666 ± 561 ± 460 ± 661 ± 759 ± 10
Female
 Height (cm)169.50 ± 3.53169.75 ± 10.96159.60 ± 6.38167.80 ± 6.14163.66 ± 3.93162.20 ± 5.11
 Weight (kg)54.00 ± 2.8270.00 ± 16.1460.20 ± 8.2264.20 ± 13.0261.33 ± 5.4277.20 ± 17.16
 BSA (m2)1.59 ± 0.031.81 ± 0.211.63 ± 0.131.72 ± 0.201.67 ± 0.081.86 ± 0.19
 BMI (kg/m2)18.83 ± 1.7724.59 ± 6.9123.64 ± 2.9222.63 ± 3.1222.90 ± 1.9229.53 ± 7.4
 EDVi (ml/m2)80 ± 1182 ± 1079 ± 1180 ± 876 ± 1173 ± 12
 LV mass index (ml/m2)54 ± 1054 ± 1255 ± 956 ± 859 ± 1260 ± 11
 Ejection fraction (%)61 ± 663 ± 564 ± 663 ± 763 ± 962 ± 11

Data are expressed as: mean ± standard deviation.

BSA: body surface area; BMI: body mass index; EDVi: end-diastolic volume index; LV: left ventricular.

Table 2:

Baseline characteristics of normal subjects studied (mean ± standard deviation)

Age class15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Males
 Height (cm)174 ± 7.88181.42 ± 7.04175 ± 14.17175.12 ± 13.87174.62 ± 6.04172.16 ± 5.60
 Weight (kg)68.62 ± 9.5080.28 ± 8.2673.14 ± 15.8674.15 ± 21.1575.12 ± 9.0479.33 ± 12.07
 BSA (m2)1.82 ± 0.142 ± 0.131.88 ± 0.271.88 ± 0.361.9 ± 0.111.94 ± 0.17
 BMI (kg/m2)22.63 ± 3.2024.39 ± 2.1923.6 ± 3.0023.51 ± 4.3524.72 ± 3.4026.67 ± 2.86
 EDVi (ml/m2)86 ± 1089 ± 1186 ± 984 ± 881 ± 1177 ± 9
 LV mass index (ml/m2)76 ± 978 ± 880 ± 1078 ± 976 ± 873 ± 11
 Ejection fraction (%)65 ± 666 ± 561 ± 460 ± 661 ± 759 ± 10
Female
 Height (cm)169.50 ± 3.53169.75 ± 10.96159.60 ± 6.38167.80 ± 6.14163.66 ± 3.93162.20 ± 5.11
 Weight (kg)54.00 ± 2.8270.00 ± 16.1460.20 ± 8.2264.20 ± 13.0261.33 ± 5.4277.20 ± 17.16
 BSA (m2)1.59 ± 0.031.81 ± 0.211.63 ± 0.131.72 ± 0.201.67 ± 0.081.86 ± 0.19
 BMI (kg/m2)18.83 ± 1.7724.59 ± 6.9123.64 ± 2.9222.63 ± 3.1222.90 ± 1.9229.53 ± 7.4
 EDVi (ml/m2)80 ± 1182 ± 1079 ± 1180 ± 876 ± 1173 ± 12
 LV mass index (ml/m2)54 ± 1054 ± 1255 ± 956 ± 859 ± 1260 ± 11
 Ejection fraction (%)61 ± 663 ± 564 ± 663 ± 763 ± 962 ± 11
Age class15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Males
 Height (cm)174 ± 7.88181.42 ± 7.04175 ± 14.17175.12 ± 13.87174.62 ± 6.04172.16 ± 5.60
 Weight (kg)68.62 ± 9.5080.28 ± 8.2673.14 ± 15.8674.15 ± 21.1575.12 ± 9.0479.33 ± 12.07
 BSA (m2)1.82 ± 0.142 ± 0.131.88 ± 0.271.88 ± 0.361.9 ± 0.111.94 ± 0.17
 BMI (kg/m2)22.63 ± 3.2024.39 ± 2.1923.6 ± 3.0023.51 ± 4.3524.72 ± 3.4026.67 ± 2.86
 EDVi (ml/m2)86 ± 1089 ± 1186 ± 984 ± 881 ± 1177 ± 9
 LV mass index (ml/m2)76 ± 978 ± 880 ± 1078 ± 976 ± 873 ± 11
 Ejection fraction (%)65 ± 666 ± 561 ± 460 ± 661 ± 759 ± 10
Female
 Height (cm)169.50 ± 3.53169.75 ± 10.96159.60 ± 6.38167.80 ± 6.14163.66 ± 3.93162.20 ± 5.11
 Weight (kg)54.00 ± 2.8270.00 ± 16.1460.20 ± 8.2264.20 ± 13.0261.33 ± 5.4277.20 ± 17.16
 BSA (m2)1.59 ± 0.031.81 ± 0.211.63 ± 0.131.72 ± 0.201.67 ± 0.081.86 ± 0.19
 BMI (kg/m2)18.83 ± 1.7724.59 ± 6.9123.64 ± 2.9222.63 ± 3.1222.90 ± 1.9229.53 ± 7.4
 EDVi (ml/m2)80 ± 1182 ± 1079 ± 1180 ± 876 ± 1173 ± 12
 LV mass index (ml/m2)54 ± 1054 ± 1255 ± 956 ± 859 ± 1260 ± 11
 Ejection fraction (%)61 ± 663 ± 564 ± 663 ± 763 ± 962 ± 11

Data are expressed as: mean ± standard deviation.

BSA: body surface area; BMI: body mass index; EDVi: end-diastolic volume index; LV: left ventricular.

Relationship between ageing and aortic elastic properties

An inverse relationship was found between MRSD of the ascending aorta and age (R2 = 0.66, 95% confidence interval (CI) −0.18 to −0.10, P < 0.001) (Figs 1 and 2). On average, the MRSD of the ascending aorta decreased by 0.14 percentile of maximal area/10–3 s every year. MRSD was not related to height (P = 0.8), body weight (P = 0.7), BMI (P = 0.4), BSA (0.6) and to systolic and diastolic arterial pressures (respectively, P = 0.9 and P = 0.8).
Females: Maximum rate of systolic distension (MRSD) of ascending and descending and distensibility of ascending and descending aortas by age decile (mean, 95% prediction). MRSD is expressed as percentage of maximal area/10–3 s. Distensibility is expressed as 10–3 mmHg–1.
Figure 2:

Females: Maximum rate of systolic distension (MRSD) of ascending and descending and distensibility of ascending and descending aortas by age decile (mean, 95% prediction). MRSD is expressed as percentage of maximal area/10–3 s. Distensibility is expressed as 10–3 mmHg–1.

MRSD in ascending aorta was related to LV end-diastolic volume index (R2 = 0.44, 95% CI 0.04–0.12, P < 0.001), while it was not associated with ejection fraction (P = 0.6) and mass index (P = 0.7).

In Tables 3 and 4 are presented the values of MRSD in the ascending and the descending aortas in the different classes of age, respectively, in males and females. No significant gender differences in MRSD in the ascending and the descending aortas was found in all classes of age.

Table 3:

Elasticity indexes in classes of age in males (mean ± standard deviation, 95% CI)

Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD11.6 ± 2.3 (9.5–13.7)12.6 ± 3.5 (9.7–15.7)10.6 ± 2.3 (8.6–13.1)6.1 ± 2.4 (4.5–9.7)5.5 ± 1.3 (4.4–6.6)4.0 ± 1.4 (2.7–5.3)
 Distensibility8.1 ± 2.9 (6.5–12.7)7.0 ± 1.4 (5.5–8.4)5.4 ± 2.1 (3.6–7.2)4.1 ± 2.8 (1.9–7.2)3.6 ± 2.2 (1.8–6.4)2.1 ± 1.2 (0.9–3.2)
Descending aorta
 MRSD15.0 ± 6.4 (10.5–19.6)11.6 ± 4.1 (6.5–16.8)7.7 ± 2.3 (5.5–9.9)6.7 ± 2.3 (4.6–8.9)6.3 ± 2.3 (4.3–8.1)4.4 ± 0.9 (3.5–5.4)
 Distensibility6.1 ± 0.6 (6.0–11.5)5.2 ± 1.5 (3.2–8.2)4.5 ± 1.1 (3.2–5.5)3.8 ± 1.8 (2.0–7.4)3.8 ± 2.1 (2.5–8.0)2.3 ± 1.2 (1.1–3.5)
 PWV1.7 ± 0.8 (1.1–2.8)3.4 ± 1.1 (2.3–4.4)3.9 ± 1.5 (2.6–5.1)4.7 ± 1.1 (3.3–5.6)5.0 ± 2.0 (3.1–6.9)6.0 ± 2.1 (4.0–8.0)
Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD11.6 ± 2.3 (9.5–13.7)12.6 ± 3.5 (9.7–15.7)10.6 ± 2.3 (8.6–13.1)6.1 ± 2.4 (4.5–9.7)5.5 ± 1.3 (4.4–6.6)4.0 ± 1.4 (2.7–5.3)
 Distensibility8.1 ± 2.9 (6.5–12.7)7.0 ± 1.4 (5.5–8.4)5.4 ± 2.1 (3.6–7.2)4.1 ± 2.8 (1.9–7.2)3.6 ± 2.2 (1.8–6.4)2.1 ± 1.2 (0.9–3.2)
Descending aorta
 MRSD15.0 ± 6.4 (10.5–19.6)11.6 ± 4.1 (6.5–16.8)7.7 ± 2.3 (5.5–9.9)6.7 ± 2.3 (4.6–8.9)6.3 ± 2.3 (4.3–8.1)4.4 ± 0.9 (3.5–5.4)
 Distensibility6.1 ± 0.6 (6.0–11.5)5.2 ± 1.5 (3.2–8.2)4.5 ± 1.1 (3.2–5.5)3.8 ± 1.8 (2.0–7.4)3.8 ± 2.1 (2.5–8.0)2.3 ± 1.2 (1.1–3.5)
 PWV1.7 ± 0.8 (1.1–2.8)3.4 ± 1.1 (2.3–4.4)3.9 ± 1.5 (2.6–5.1)4.7 ± 1.1 (3.3–5.6)5.0 ± 2.0 (3.1–6.9)6.0 ± 2.1 (4.0–8.0)

Data are expressed as: mean ± standard deviation (95% CI).

MRSD: maximum rate of systolic distension; PWV: pulse wave velocity.

Table 3:

Elasticity indexes in classes of age in males (mean ± standard deviation, 95% CI)

Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD11.6 ± 2.3 (9.5–13.7)12.6 ± 3.5 (9.7–15.7)10.6 ± 2.3 (8.6–13.1)6.1 ± 2.4 (4.5–9.7)5.5 ± 1.3 (4.4–6.6)4.0 ± 1.4 (2.7–5.3)
 Distensibility8.1 ± 2.9 (6.5–12.7)7.0 ± 1.4 (5.5–8.4)5.4 ± 2.1 (3.6–7.2)4.1 ± 2.8 (1.9–7.2)3.6 ± 2.2 (1.8–6.4)2.1 ± 1.2 (0.9–3.2)
Descending aorta
 MRSD15.0 ± 6.4 (10.5–19.6)11.6 ± 4.1 (6.5–16.8)7.7 ± 2.3 (5.5–9.9)6.7 ± 2.3 (4.6–8.9)6.3 ± 2.3 (4.3–8.1)4.4 ± 0.9 (3.5–5.4)
 Distensibility6.1 ± 0.6 (6.0–11.5)5.2 ± 1.5 (3.2–8.2)4.5 ± 1.1 (3.2–5.5)3.8 ± 1.8 (2.0–7.4)3.8 ± 2.1 (2.5–8.0)2.3 ± 1.2 (1.1–3.5)
 PWV1.7 ± 0.8 (1.1–2.8)3.4 ± 1.1 (2.3–4.4)3.9 ± 1.5 (2.6–5.1)4.7 ± 1.1 (3.3–5.6)5.0 ± 2.0 (3.1–6.9)6.0 ± 2.1 (4.0–8.0)
Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD11.6 ± 2.3 (9.5–13.7)12.6 ± 3.5 (9.7–15.7)10.6 ± 2.3 (8.6–13.1)6.1 ± 2.4 (4.5–9.7)5.5 ± 1.3 (4.4–6.6)4.0 ± 1.4 (2.7–5.3)
 Distensibility8.1 ± 2.9 (6.5–12.7)7.0 ± 1.4 (5.5–8.4)5.4 ± 2.1 (3.6–7.2)4.1 ± 2.8 (1.9–7.2)3.6 ± 2.2 (1.8–6.4)2.1 ± 1.2 (0.9–3.2)
Descending aorta
 MRSD15.0 ± 6.4 (10.5–19.6)11.6 ± 4.1 (6.5–16.8)7.7 ± 2.3 (5.5–9.9)6.7 ± 2.3 (4.6–8.9)6.3 ± 2.3 (4.3–8.1)4.4 ± 0.9 (3.5–5.4)
 Distensibility6.1 ± 0.6 (6.0–11.5)5.2 ± 1.5 (3.2–8.2)4.5 ± 1.1 (3.2–5.5)3.8 ± 1.8 (2.0–7.4)3.8 ± 2.1 (2.5–8.0)2.3 ± 1.2 (1.1–3.5)
 PWV1.7 ± 0.8 (1.1–2.8)3.4 ± 1.1 (2.3–4.4)3.9 ± 1.5 (2.6–5.1)4.7 ± 1.1 (3.3–5.6)5.0 ± 2.0 (3.1–6.9)6.0 ± 2.1 (4.0–8.0)

Data are expressed as: mean ± standard deviation (95% CI).

MRSD: maximum rate of systolic distension; PWV: pulse wave velocity.

Table 4:

Elasticity indexes in classes of age in females (mean ± standard deviation, 95% CI)

Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD13.4 ± 1.9 (11.0–15.7)11.2 ± 3.8 (9.2–14.4)7.41 ± 1.9 (5.4–9.4)6.7 ± 1.4 (5.1–8.4)4.6 ± 1.7 (2.5–6.8)3.3 ± 1.99 (1.2–5.8)
 Distensibility13.2 ± 3.4 (8.3–16.3)8.4 ± 2.6 (6.0–11.3)6.2 ± 3.7 (4.2–10.8)5.4 ± 3.7 (3.8–7.9)5.1 ± 1.8 (1.6–6.0)1.91 ± 0.4 (1.16–2.7)
Descending aorta
 MRSD11.7 ± 1.3 (8.1–14.4)9.4 ± 1.7 (6.2–10.6)7.8 ± 2.3 (5.4–10.1)7.6 ± 1.9 (5.2–10.1)5.8 ± 2.8 (2.3–9.2)4.5 ± 2.2 (1.1–8.9)
 Distensibility11.0 ± 1.0 (6.0–16.0)6.0 ± 1.2 (4.0–8.0)5.0 ± 1.7 (4.0–7.0)5.0 ± 1.5 (3.0–6.0)4.0 ± 1.1 (2.0–6.0)3.0 ± 1.0 (1.6–4.2)
 PWV2.4 ± 1.4 (1.2–4.3)2.6 ± 1.4 (1.4–4.4)3.6 ± 1.4 (2.4–4.7)4.2 ± 1.1 (3.2–5.6)4.2 ± 1.1 (2.9–5.6)4.9 ± 0.6 (3.7–6.1)
Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD13.4 ± 1.9 (11.0–15.7)11.2 ± 3.8 (9.2–14.4)7.41 ± 1.9 (5.4–9.4)6.7 ± 1.4 (5.1–8.4)4.6 ± 1.7 (2.5–6.8)3.3 ± 1.99 (1.2–5.8)
 Distensibility13.2 ± 3.4 (8.3–16.3)8.4 ± 2.6 (6.0–11.3)6.2 ± 3.7 (4.2–10.8)5.4 ± 3.7 (3.8–7.9)5.1 ± 1.8 (1.6–6.0)1.91 ± 0.4 (1.16–2.7)
Descending aorta
 MRSD11.7 ± 1.3 (8.1–14.4)9.4 ± 1.7 (6.2–10.6)7.8 ± 2.3 (5.4–10.1)7.6 ± 1.9 (5.2–10.1)5.8 ± 2.8 (2.3–9.2)4.5 ± 2.2 (1.1–8.9)
 Distensibility11.0 ± 1.0 (6.0–16.0)6.0 ± 1.2 (4.0–8.0)5.0 ± 1.7 (4.0–7.0)5.0 ± 1.5 (3.0–6.0)4.0 ± 1.1 (2.0–6.0)3.0 ± 1.0 (1.6–4.2)
 PWV2.4 ± 1.4 (1.2–4.3)2.6 ± 1.4 (1.4–4.4)3.6 ± 1.4 (2.4–4.7)4.2 ± 1.1 (3.2–5.6)4.2 ± 1.1 (2.9–5.6)4.9 ± 0.6 (3.7–6.1)

Data are expressed as: mean ± standard deviation or median (IQR, inter-quartile range).

MRSD: maximum rate of systolic distension; PWV: pulse wave velocity.

Table 4:

Elasticity indexes in classes of age in females (mean ± standard deviation, 95% CI)

Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD13.4 ± 1.9 (11.0–15.7)11.2 ± 3.8 (9.2–14.4)7.41 ± 1.9 (5.4–9.4)6.7 ± 1.4 (5.1–8.4)4.6 ± 1.7 (2.5–6.8)3.3 ± 1.99 (1.2–5.8)
 Distensibility13.2 ± 3.4 (8.3–16.3)8.4 ± 2.6 (6.0–11.3)6.2 ± 3.7 (4.2–10.8)5.4 ± 3.7 (3.8–7.9)5.1 ± 1.8 (1.6–6.0)1.91 ± 0.4 (1.16–2.7)
Descending aorta
 MRSD11.7 ± 1.3 (8.1–14.4)9.4 ± 1.7 (6.2–10.6)7.8 ± 2.3 (5.4–10.1)7.6 ± 1.9 (5.2–10.1)5.8 ± 2.8 (2.3–9.2)4.5 ± 2.2 (1.1–8.9)
 Distensibility11.0 ± 1.0 (6.0–16.0)6.0 ± 1.2 (4.0–8.0)5.0 ± 1.7 (4.0–7.0)5.0 ± 1.5 (3.0–6.0)4.0 ± 1.1 (2.0–6.0)3.0 ± 1.0 (1.6–4.2)
 PWV2.4 ± 1.4 (1.2–4.3)2.6 ± 1.4 (1.4–4.4)3.6 ± 1.4 (2.4–4.7)4.2 ± 1.1 (3.2–5.6)4.2 ± 1.1 (2.9–5.6)4.9 ± 0.6 (3.7–6.1)
Class of age15–20 years21–30 years31–40 years41–50 years51–60 years>60 years
Ascending aorta
 MRSD13.4 ± 1.9 (11.0–15.7)11.2 ± 3.8 (9.2–14.4)7.41 ± 1.9 (5.4–9.4)6.7 ± 1.4 (5.1–8.4)4.6 ± 1.7 (2.5–6.8)3.3 ± 1.99 (1.2–5.8)
 Distensibility13.2 ± 3.4 (8.3–16.3)8.4 ± 2.6 (6.0–11.3)6.2 ± 3.7 (4.2–10.8)5.4 ± 3.7 (3.8–7.9)5.1 ± 1.8 (1.6–6.0)1.91 ± 0.4 (1.16–2.7)
Descending aorta
 MRSD11.7 ± 1.3 (8.1–14.4)9.4 ± 1.7 (6.2–10.6)7.8 ± 2.3 (5.4–10.1)7.6 ± 1.9 (5.2–10.1)5.8 ± 2.8 (2.3–9.2)4.5 ± 2.2 (1.1–8.9)
 Distensibility11.0 ± 1.0 (6.0–16.0)6.0 ± 1.2 (4.0–8.0)5.0 ± 1.7 (4.0–7.0)5.0 ± 1.5 (3.0–6.0)4.0 ± 1.1 (2.0–6.0)3.0 ± 1.0 (1.6–4.2)
 PWV2.4 ± 1.4 (1.2–4.3)2.6 ± 1.4 (1.4–4.4)3.6 ± 1.4 (2.4–4.7)4.2 ± 1.1 (3.2–5.6)4.2 ± 1.1 (2.9–5.6)4.9 ± 0.6 (3.7–6.1)

Data are expressed as: mean ± standard deviation or median (IQR, inter-quartile range).

MRSD: maximum rate of systolic distension; PWV: pulse wave velocity.

MRSD measured in both the ascending and the descending aortas demonstrated a significant decrease passing between the first three classes of age (<31 years) and other classes (P < 0.001). MRSD values increased from 10 to 30 years reaching a plateau at 20 years, and then decreased progressively.

Distensibility in the ascending aorta decreased with ageing (Figs 1 and 2). Subjects aged <31 showed higher distensibility than those aged >40, as evidenced in Table 3. Distensibility was related to weight (P = 0.005) and BSA (P < 0.005). It was inversely related to aortic arch length (P = 0.002) and directly related to LV mass index (R2 = 0.55, 95% CI 0.4–0.8, P < 0.001) and LV end-diastolic volume index (P = 0.002).

PWV was directly related to age (R2 = 0.49, 95% CI 0.03–0.08, P < 0.001) (Fig. 3), while it was not related to height (P = 0.92), weight (P = 0.24), BMI (P = 0.19) and BSA (P = 0.33). Similarly, PWV was not related to LV ventricular parameters as ejection fraction (P = 0.13), end-diastolic volume index (P = 0.32) and mass index (P = 0.43). A relation between PWV and the maximal diameter of the ascending aorta was found (R2 = 0.4, 95% CI 0.19–1.14, P < 0.007). As evidenced in tables, PWV increased progressively from the first class of age to the highest class of age.
Pulse wave velocity (PWV) by age decile (mean, 95% prediction) in males (left panel) and females (right panel). PWV is expressed as m/s.
Figure 3:

Pulse wave velocity (PWV) by age decile (mean, 95% prediction) in males (left panel) and females (right panel). PWV is expressed as m/s.

No relation was found between MRSD of both the ascending (P = 0.43) and the descending aortas (P = 0.31) and PWV. Instead, a relation was found between ascending aorta distensibility and PWV (P = 0.005). A relation was also found between MRSD of the ascending aorta and distensibility (R2 = 0.5, 95% CI 0.3–0.7, P < 0.001).

DISCUSSION

The main findings of this study may be reassumed as: (i) MRSD and distensibility showed a peak value between 21 and 30 years in males and 15 and 20 in females, then decreased progressively reaching the lowest values in subjects aged >60 years; (ii) PWV increased progressively with age; (iii) cross-sectional area of both the ascending and the descending aortas increased significantly from the lowest to the highest class of age.

MRSD was recently proposed as a new index of elastic properties assessing the compliance of the aortic wall during the systolic distension provoked by LV ejection. A lower MRSD value indicates that the kinetic energy of systolic ejection is transformed more slowly into potential energy. The imbalance between the uptake and release of energy by the aortic wall (Windkessel function), due to impaired elasticity, could slow the centripetal acceleration of the aortic wall, increasing aortic wall stress and progressively leading to dilatation of the aortic wall. In the current study, MRSD was independent of aortic size, arterial pressure, LV mass index and body parameters such as weight, height, BMI and BSA. Distensibility, instead, is directly related to blood pressure and body conformation parameters and aortic size. However, we found that also ascending aorta distensibility reached a peak at 20 years and then decreased progressively with the increase of age.

In this study, PWV was measured by the analysis of SSFP images in which cross sections of the ascending and descending aortas were depicted. The time of pulse wave to propagate from the ascending to the descending aortas was measured as the time gap to reach half of the maximal cross-sectional area between the ascending and the descending aortas. Then PWV was measured dividing the pulse wave propagation time by the length of aortic arch. Thus, the propagation of aortic wall deformation over time was directly evaluated. Unlike MRSD and distensibility, PWV increased progressively with age. The impairment of the elastic properties of the thoracic aorta in elder patients is concordant with the previous demonstration that ageing is accompanied by a destruction of the elastic vascular structure and that elastogenesis occurs throughout life but with constantly decreasing levels [17]. In this study, for every class of age, physiological range of values of aortic elasticity indexes has been determined. Then, these physiological ranges of values may be useful in clinical practice permitting the early identification of subjects with impaired aortic elastic properties. These subjects may be potentially at higher risk of development of ascending aorta aneurysm and acute aortic complications.

The association between impaired aortic elastic properties and some clinical conditions at higher risk of aortic complication was previously demonstrated, as in Marfan disease and in subjects with bicuspid aortic valve [15, 18].

Some study limitations should be mentioned. First, B-natriuretic peptide was not considered an exclusion criterion for the selection of the study population, and then we cannot completely rule out the presence of initial cardiac disease. However, we enrolled subjects with no family history of cardiac and aortic disease, no risk factors for coronary artery disease and with normal echocardiogram including diastolic function. Moreover, the evaluation of cardiac function by CMR demonstrated a normal systolic function, ventricular volumes and mass index in all the subjects using the referral normality values [19, 20]. Thus, the probability of including subjects affected by initial cardiac disease in the study population was very low.

Another limitation was that measurement of pressure was performed not contemporary to CMR acquisitions because we could not perform it inside the magnet. This may potentially alter the quantification of distensibility. However, the time gap between pressure measurement and CMR acquisition was very short (immediately after the end of the examinations) and in the same conditions. Moreover, we also used two indexes MRSD and PWV, which are independent of the measurement of blood pressure.

CONCLUSIONS

This study confirmed that physiological ageing is associated with a progressive impairment of the elastic properties of the aortic wall. Results of this study may be useful for the early identification of subjects with impaired aortic wall properties providing referral values of elasticity indexes assessed by CMR in different classes of age. CMR is a valuable tool for the assessment of aortic wall function.

Funding

This study was supported by research funding from the Fondazione G.Monasterio CNR-Regione Toscana.

Conflict of interest: none declared.

REFERENCES

1

Laurent
S
Cockcroft
J
Van Bortel
L
Expert consensus document on arterial stiffness: methodological issues and clinical applications
Eur Heart J
2006
27
2588
605

2

Mitchell
GF
Hwang
SJ
Vasan
RS
Larson
MG
Pencina
MJ
Hamburg
NM
et al.
Arterial stiffness and cardiovascular events: the Framingham Heart Study
Circulation
2010
121
505
11

3

Vlachopoulos
C
Aznaouridis
K
Stefanadis
C
Prediction of cardiovascular events and all-cause mortality with arterial stiffness: a systematic review and meta-analysis
J Am Coll Cardiol
2010
55
1318
27

4

Redheuil
A
Yu
WC
Wu
CO
Mousseaux
E
de Cesare
A
Yan
R
et al.
Reduced ascending aortic strain and distensibility: earliest manifestations of vascular aging in humans
Hypertension
2010
55
319
26

5

Rogers
WJ
Hu
YL
Coast
D
Vido
DA
Kramer
CM
Pyeritz
RE
et al.
Age-associated changes in regional aortic pulse wave velocity
J Am Coll Cardiol
2001
38
1123
9

6

Mohiaddin
RH
Firmin
DN
Longmore
DB
Age-related changes of human aortic flow wave velocity measured noninvasively by magnetic resonance imaging
J Appl Physiol
1993
74
492
7

7

Groenink
M
de Roos
A
Mulder
BJ
Verbeeten
B
Jr
Timmermans
J
Zwinderman
AH
et al.
Biophysical properties of the normal-sized aorta in patients with Marfan syndrome: evaluation with MR flow mapping
Radiology
2001
219
535
40

8

Hirata
K
Triposkiadis
F
Sparks
E
Bowen
J
Wooley
CF
Boudoulas
H
The Marfan syndrome: abnormal aortic elastic properties
J Am Coll Cardiol
1992
18
57
63

9

Jeremy
RW
Huang
H
Hwa
J
McCarron
H
Hughes
CF
Richards
JG
Relation between age, arterial distensibility, and aortic dilatation in the Marfan syndrome
Am J Cardiol
1994
74
369
73

10

Savolainen
A
Keto
P
Hekali
P
Nisula
L
Kaitila
I
Viitasalo
M
et al.
Aortic distensibility in children with the Marfan syndrome
Am J Cardiol
1992
70
691
3

11

Debl
K
Djavidani
B
Buchner
S
Poschenrieder
F
Schmid
FX
Kobuch
R
et al.
Dilatation of the ascending aorta in bicuspid aortic valve disease: a magnetic resonance imaging study
Clin Res Cardiol
2009
98
114
20

12

Aquaro
GD
Ait-Ali
L
Basso
ML
Lombardi
M
Pingitore
A
Festa
P
Elastic properties of aortic wall in patients with bicuspid aortic valve by magnetic resonance imaging
Am J Cardiol
2011
108
81
7

13

Nollen
GJ
Groenink
M
Tijssen
JG
Van Der Wall
EE
Mulder
BJ
Aortic stiffness and diameter predict progressive aortic dilatation in patients with Marfan syndrome
Eur Heart J
2004
25
1146
52

14

Chen
CH
Nevo
E
Fetics
B
Pak
PH
Yin
FC
Maughan
WL
et al.
Estimation of central aortic pressure waveform by mathematical transformation of radial tonometry pressure. Validation of generalized transfer function
Circulation
1997
95
1827
36

15

Lionetti
V
Guiducci
L
Simioniuc
A
Aquaro
GD
Simi
C
De Marchi
D
et al.
Mismatch between uniform increase in cardiac glucose uptake and regional contractile dysfunction in pacing-induced heart failure
Am J Physiol Heart Circ Physiol
2007
293
H2747
56

16

Grotenhuis
HB
Ottenkamp
J
Westenberg
JJ
Bax
JJ
Kroft
LJ
de Roos
A
Reduced aortic elasticity and dilatation are associated with aortic regurgitation and left ventricular hypertrophy in nonstenotic bicuspid aortic valve patients
J Am Coll Cardiol
2007
49
1660
5

17

Fritze
O
Romero
B
Schleicher
M
Jacob
MP
Oh
D
Starcher
B
et al.
Age-related changes in the elastic tissue of the human aorta
J Vasc Res
2012
49
77
86

18

Westenberg
JJ
Scholte
AJ
Vaskova
Z
van der Geest
RJ
Groenink
M
Labadie
G
et al.
Age-related and regional changes of aortic stiffness in the Marfan syndrome: assessment with velocity-encoded MRI
J Magn Reson Imaging
2011
34
526
31

19

Maceira
AM
Prasad
SK
Khan
M
Pennell
DJ
Normalized left ventricular systolic and diastolic function by steady state free precession cardiovascular magnetic resonance
J Cardiovasc Magn Reson
2006
8
417
26

20

Maceira
AM
Prasad
SK
Khan
M
Pennell
DJ
Reference right ventricular systolic and diastolic function normalized to age, gender and body surface area from steady-state free precession cardiovascular magnetic resonance
Eur Heart J
2006
27
2879
88

Author notes

G.D.A. and A.C. equally contributed to this study.