Skip to main content

Late effects of pediatric hematopoietic stem cell transplantation on left ventricular function, aortic stiffness and myocardial tissue characteristics



Pediatric hematopoietic stem cell transplantation (HSCT) recipients are at increased risk of cardiovascular disease later in life. As HSCT survival has significantly improved, with a growing number of HSCT indications, tailored screening strategies for HSCT-related late effects are warranted. Little is known regarding the value of cardiovascular magnetic resonance (CMR) for early identification of high-risk patients after HSCT, before symptomatic cardiovascular disease manifests. This study aimed to assess CMR-derived left ventricular (LV) systolic and diastolic function, aortic stiffness and myocardial tissue characteristics in young adults who received HSCT during childhood.


Sixteen patients (22.1 ± 1.5 years) treated with HSCT during childhood and 16 healthy controls (22.1 ± 1.8 years) underwent 3 T CMR. LV systolic and diastolic function were measured as LV ejection fraction (LVEF), the ratio of transmitral early and late peak filling rate (E/A), the estimated LV filling pressure (E/Ea) and global longitudinal and circumferential systolic strain and diastolic strain rates, using balanced steady-state free precession cine CMR and 2D velocity-encoded CMR over the mitral valve. Aortic stiffness, myocardial fibrosis and steatosis were assessed with 2D velocity-encoded CMR, native T1 mapping and proton CMR spectroscopy (1H-CMRS), respectively.


In the patient compared to the control group, E/Ea (9.92 ± 3.42 vs. 7.24 ± 2.29, P = 0.004) was higher, LVEF (54 ± 6% vs. 58 ± 5%, P = 0.055) and global longitudinal strain (GLS) ( -20.7 ± 3.5% vs. -22.9 ± 3.0%, P = 0.063) tended to be lower, while aortic pulse wave velocity (4.40 ± 0.26 vs. 4.29 ± 0.29 m/s, P = 0.29), native T1 (1211 ± 36 vs. 1227 ± 28 ms, P = 0.16) and myocardial triglyceride content (0.47 ± 0.18 vs. 0.50 ± 0.13%, P = 0.202) were comparable. There were no differences between patients and controls in E/A (2.76 ± 0.92 vs. 2.97 ± 0.91, P = 0.60) and diastolic strain rates.


In young adults who received HSCT during childhood, LV diastolic function was decreased (higher estimated LV filling pressure) and LV systolic function (LVEF and GLS) tended to be reduced as compared to healthy controls, whereas no concomitant differences were found in aortic stiffness and myocardial tissue characteristics. When using CMR, assessment of LV diastolic function in particular is important for early detection of patients at risk of HSCT-related cardiovascular disease, which may warrant closer surveillance.


Hematopoietic stem cell transplantation (HSCT) recipients are exposed to several pre-transplant and/or HSCT-related therapies which may increase the risk of cardiovascular disease [1, 2]. As HSCT survival has significantly improved over the last decades [3, 4], with an increasing number of HSCT indications for both malignant and non-malignant disease [5], targeted follow-up strategies for the HSCT population are needed [6]. Recently, several international working groups have been established aimed at a greater understanding of the late effects including arterial disease and cardiac dysfunction [7, 8]. Most of the available studies on late effects involve HSCT in adults. Adequate screening in young HSCT recipients is even more challenging [9].

According to current guidelines [10], pediatric HSCT recipients who have an increased susceptibility to complications later in life, based on pre-existing comorbidities, pre-transplant exposures, the HSCT preparative regimen, post-transplant complications such as graft-versus-host-disease, or relapse of the primary disease, are selected for patient-specific follow-up programs. Identification of imaging markers which indicate subclinical disease would be supportive in the detection of high-risk patient groups for closer monitoring or targeted therapy [6]. In general, clinical follow-up after cardiotoxic exposures comprises echocardiography for left ventricular (LV) systolic and diastolic function. Less is known regarding the value of cardiovascular magnetic resonance (CMR) for screening of late effects [11]. CMR may be suited for comprehensive evaluation of subclinical deteriorations within the cardiovascular system after HSCT, that may be present before overt LV functional abnormalities arise.

Pre-transplant or HSCT-related cardiotoxic exposures may cause endothelial damage leading to increased aortic stiffness and myocyte cell death with reactive interstitial fibrosis [1, 2]. In addition, the immunosuppressive therapies in allogeneic HSCT increase the susceptibility to developing the metabolic syndrome at young age [12]. Aortic stiffening is known to occur in relation to normal, physiological ageing [13], but will be more progressive in response to hypertension, dyslipidemia or hyperglycemia [14] and possibly due to iron overload [15]. Increased aortic stiffness induces LV concentric remodeling and is recognized as an independent predictor for cardiovascular events [16]. Myocardial diffuse fibrosis in different types of cardiomyopathies is considered to reflect subclinical disease before cardiac dysfunction becomes manifest [17]. The metabolic derangements among allogeneic HSCT recipients may predispose to myocardial steatosis. In individuals with the metabolic syndrome, myocardial steatosis has been associated with LV remodeling [18].

We hypothesize that, when using CMR, subclinical deteriorations in LV function, aortic stiffness and/or myocardial tissue characteristics can be detected in young adults who received HSCT during childhood. Therefore, CMR may be suitable for early identification of patients at increased risk of developing HSCT-related cardiovascular disease. Accordingly, the aim of this study is to assess CMR-derived LV systolic and diastolic function, aortic stiffness and myocardial fibrosis and steatosis in young adults who have received pediatric HSCT and to compare these measures with those in healthy controls in the same age range.


Study population

The patient group consisted of young adults (18 to 25 years old), who received HSCT for malignant or non-malignant disease during childhood. Patients were recruited from the outpatient clinic for screening and treatment of late effects of childhood cancer and/or HSCT of the Internal Medicine Department (Leiden University Medical Centre, the Netherlands). Healthy controls were recruited by local advertising in Leiden University, the Netherlands. The control group was in the same age range as the patients and was sex-matched. Laboratory measures in the patient group were performed based on clinical indication and were typically measured within 1 year prior or after CMR examination. No blood samples were drawn in the healthy control group.

CMR acquisition

The study participants underwent 3 T CMR (Ingenia, Philips Healthcare, Best, the Netherlands), with a dStream Torso anterior coil and a FlexCoverage posterior coil in the table top, resulting in up to 32 coil elements for signal reception. The protocol consisted of standard electrocardiographic (ECG)-triggered two-, three- and four-chamber and short-axis cine balanced steady-state free precession (bSSFP) CMR and ECG-gated gradient-echo 2D velocity-encoded CMR over the mitral valve to quantify LV structure, systolic and diastolic function; 2D velocity-encoded CMR transecting the aortic arch and abdominal aorta to derive aortic stiffness; cardiac native T1 mapping to assess diffuse fibrosis and proton cardiovascular magnetic resonance spectroscopy (1H-CMRS) to measure the myocardial triglyceride content. For standardization of the measurement of myocardial triglyceride content, all participants were asked to fast for 6 h and the CMR examinations were scheduled at fixed times (evenings). No CMR contrast material was used.

For the bSSFP cines, typical field-of-view (FOV) was 350 × 350 mm2 (long-axis) and 400 × 352 mm2 (short-axis), acquired voxel size 2.0 × 1.6 mm2 (long-axis) and 1.5 × 1.5 mm2 (short-axis), slice thickness 8 mm, echo/repetition time (TE/TR) 1.5/3.0 ms, flip angle 45°, number of phases 30 (long-axis) and 35 (short-axis). For short-axis bSSFP cine, the complete LV was imaged, using 14–16 slices, depending on the size of the LV, with 2 slices per breath-hold. Free-breathing through-plane 2D velocity-encoded (venc 100 cm/s) CMR was acquired perpendicular to the mitral valve, with FOV 350 × 350 mm2, acquired voxel size 2.5 × 2.5mm2, slice thickness 8 mm, flip angle 10°, TE/TR 2.9/4.6 ms, number of phases 40.

Aortic pulse wave velocity (PWV) was quantified as a measure of aortic stiffness. For aortic PWV, first, a double-oblique sagittal scout view of the aorta was obtained. Then, two free-breathing through-plane 2D velocity-encoded CMR scans were acquired, one transecting the ascending aorta (venc 150 cm/s) and one transecting the abdominal aorta, above the aortic bifurcation (VENC 100 cm/s). Typical imaging parameters were: FOV 350 × 282 mm2, slice thickness 8 mm, acquired voxel 2.8 × 2.8 mm2, flip angle 20°, TE/TR 2.5/4.4 ms, temporal resolution 10 ms.

Native T1 mapping was acquired in the mid-ventricular short-axis slice, with breath-holding, using the 5s(3s)3s modified Look-Locker inversion recovery (MOLLI) scheme. Typical imaging parameters were: FOV 350 × 300 mm2, slice thickness 8 mm, acquired voxel size 2.1 × 2.1 mm2, flip angle 20°, TE/TR 1.1/2.3 ms, SENSitivity Encoding (SENSE) factor 2. 1H-MRS was performed as described previously [19]. In summary, a voxel of 40 x 15 x 25 mm3 was placed in the interventricular septum. For the acquisition with and without water suppression, 48 and 6 signal averages were obtained, respectively. ECG-triggering was used to acquire 1H-CMRS at 200 ms after the R-wave, and a respiratory navigator, tracking the lung-liver interface, for acquisition at end-expiration. A high permittivity pad was placed on the chest for improved signal-to-noise ratio.

Image analysis

CMR data were analyzed using MASS Research Software V2016-EXP (Leiden University Medical Center, the Netherlands) for LV structure and function and aortic PWV, custom-made software for further analysis of the aortic velocity-time curves [20], Medis Suite 3.0 (Medis Medical Imaging systems, Leiden, the Netherlands) for LV systolic and diastolic strain (QStrain 2.0) and native T1 (QMap 2.2.18), and the Java-based magnetic resonance user interface (jMRUI v5.0; MRUI Consortium) for 1H-CMRS. The image analysis was blinded to patient or healthy control status.

For LV mass and volumes, the endocardial and epicardial LV borders were manually outlined in the end-diastolic and end-systolic phase, with exclusion of the LV papillary muscles. For the feature tracking based strain calculations, the manually annotated endocardial LV borders were automatically tracked throughout the cardiac cycle. Global longitudinal strain (GLS), global longitudinal peak systolic strain rate (GLSR-S) and global longitudinal early peak diastolic strain rate (GLSR-E) were calculated based on the two-, three- and four-chamber cine images. Global circumferential strain (GCS), global circumferential peak systolic strain rate (GCSR-S) and global circumferential early peak diastolic strain rate (GCSR-E) were extracted from the mid-ventricular short-axis cine slice (Fig. 1). The ratio of the transmitral early and late peak filling rate (E/A ratio) and the transmitral early peak maximum velocity were derived from the 2D velocity-encoded scans (Fig. 1), as described previously [21]. The transmitral filling rate was measured after correction for the through-plane background velocity of the LV myocardial wall. Furthermore, the early diastolic mitral septal tissue velocity (Ea) was extracted from the four-chamber cine images. Subsequently, the estimated LV filling pressure (ratio of transmitral early peak maximum velocity without through-plane myocardial motion correction and early diastolic mitral septal tissue velocity) was calculated (E/Ea ratio).

Fig. 1
figure 1

A 22-year-old man, who was transplanted for a non-malignant bone marrow failure disorder at the age of 8 years. a the transmitral E/A ratio (early peak filling rate/late peak filling rate) was measured using 2D velocity-encoded CMR (left panel) and Ea (early peak diastolic mitral septal tissue velocity) was derived from the four-chamber long-axis relaxation (right panel). LV filling pressure was estimated by the ratio of the transmitral early peak maximum velocity and the early peak diastolic mitral septal tissue velocity. b From two, three- and four-chamber and mid-ventricular short-axis cine CMR (left panel), the longitudinal and circumferential strain and strain rate curves were extracted (right panel). The myocardial features at the endocardial borders (red dots), which were automatically tracked throughout the cardiac cycle (green lines), were manually annotated in the end-diastolic and end-systolic phase. GLS: global longitudinal strain; GCS: global circumferential strain; GLSR-S: global longitudinal peak systolic strain rate; GCSR-E: global circumferential early peak diastolic strain rate

Aortic PWV was calculated by dividing the distance between ascending and abdominal aorta by the transit time of the onset of the systolic velocity wave front (Fig. 2), as described previously [20]. In short, the aortic path length was measured manually along the aortic centerline on the double-oblique sagittal aorta scout scan. The onset of the systolic wave front was automatically determined from the resulting velocity graph by the intersection point of the constant horizontal diastolic velocity and upslope of the systolic wave front, modeled by linear regression (using the velocity values between 20 and 80% of the total range) along the upslope. T1 maps were constructed after manual in-plane motion correction of the T1 images, using a pixel-wise, mono-exponential three-parameter fit for the T1 relaxation curve. For the measurement of native T1, a region-of-interest in the mid-ventricular septum was drawn (Fig. 2). Myocardial triglyceride content was expressed as the percentage of triglyceride methyl (at 0.9 ppm) and triglyceride methylene (at 1.3 ppm) relative to the sum of the triglyceride signal and the unsuppressed water signal (at 4.7 ppm) (Fig. 2) [22].

Fig. 2
figure 2

The same patient as in Fig. 1 is presented. a Aortic pulse wave velocity was calculated from through-plane 2D velocity-encoded CMR transecting the ascending aorta (red) and the abdominal aorta, above the aortic bifurcation (green) (left panel), according to: aortic pulse wave velocity = ∆x/∆t, with ∆x: the distance between the ascending and abdominal aorta (yellow dotted line) and ∆t: transit time of the onset of the systolic velocity wave front (black arrow) (right panel). b Proton-cardioavascular magnetic resonance spectroscopy (1H-CMRS) was used to measure the myocardial triglyceride content. The voxel of interest was placed in the mid-ventricular septum (yellow box) (left panel). Myocardial triglyceride content was calculated as Tg-(CH2)n and Tg-CH3 relative to the sum of the triglyceride and the unsuppressed water signal (not shown). Triglycerides were measured using the water-suppressed spectrum (right panel). c Native T1 was measured in the mid-ventricular septal segments in short-axis view (black, dotted region of interest)


Statistical analyses were performed in SPSS 23 (International Business Machines, Armonk, New York, USA). Normal distribution was checked using the Shapiro-Wilk test. Differences between groups were tested for statistical significance using the Student’s t-test or the Mann-Whitney U test for normally and non-normally distributed variables, respectively. Patients with prior HSCT vs. controls were compared and additional analyses were performed for patients with HSCT for malignancies vs. controls. Levels for statistical significance were set at P < 0.05 and all tests were two-sided. Additionally, we performed a post-hoc power analysis (alfa 0.05, two-sided).


In total, sixteen patients (22.1 ± 1.5 years, 11/16 (69%) men) and sixteen healthy controls (22.1 ± 1.8 years, 11/16 (69%) men) were included and analyzed. In one patient, 1H-CMRS data were not analyzed due to insufficient quality. Also, in one patient 2D velocity-encoded CMR over the mitral valve was not acquired due to imaging time constraints. Height and body surface area (BSA) in the patient group were lower compared to the healthy control group. The anthropometric characteristics are presented in Table 1.

Table 1 Demographic and anthropometric characteristics

In ten patients, HSCT was indicated for a malignant disorder. Seven patients were treated with anthracyclines prior to HSCT, eight received total body irradiation (typically unfractionated 7.5 Gy or 2 × 6 Gy) and fourteen were given high-dose cyclophosphamide (> 1 g/m2) for HSCT conditioning. Five patients had a ferritin level above 250 μmol/L. In two of them, cardiac T2* was assessed upon clinical indication. In both, the T2* values were not suggestive of myocardial iron deposition. Clinical and biochemical patients characteristics are presented in Tables 2 and 3.

Table 2 Clinical patient characteristics
Table 3 Biochemical patient characteristics

Patients as compared to the healthy controls had a higher E/Ea ratio (9.92 ± 3.42 vs. 7.24 ± 2.29, P = 0.004), while the E/A ratio (2.76 ± 0.92 vs. 2.97 ± 0.91, P = 0.599) and diastolic strain rates were comparable. There was a trend towards a lower LV ejection fraction (LVEF) (54 ± 6 vs. 58 ± 5%, P = 0.055) and lower GLS (-20.7 ± 3.5 vs. -22.9 ± 3.0%, P = 0.063), whereas GCS was preserved (-23.2 ± 3.6 vs. -23.9 ± 3.5%, P = 0.587). Stroke volume (83 ± 15 vs. 101 ± 17 mL, P = 0.003) but also cardiac output (5.2 ± 1.0 vs. 6.5 ± 1.0 L, P = 0.001) and cardiac index (2.9 ± 0.6 vs. 3.4 ± 0.5 L/m2, P = 0.021) were lower. In contrast, aortic PWV (4.40 ± 0.26 vs. 4.29 ± 0.29 m/s, P = 0.288), LV concentricity (0.62 ± 0.10 vs. 0.61 ± 0.08 g/mL, P = 0.867), native T1 (1211 ± 36 vs. 1227 ± 28 ms, P = 0.158) and myocardial triglyceride content (0.47 ± 0.18 vs. 0.50 ± 0.13%, P = 0.202) were comparable. CMR results are presented in Table 4 and Fig. 3.

Table 4 CMR parameters
Fig. 3
figure 3

CMR measurements with means and 95% confidence intervals and median and interquartile ranges for normally and non-normally distributed data, respectively. Despite a non-significantly lower left ventricular (LV) systolic function (a-c) and lower LV diastolic function (d-e) as indicated by the increased estimated LV filling pressure, the CMR parameters for aortic stiffness (f) and LV structure (g) and myocardial tissue characteristics (h-i) were comparable for the patients who received HSCT and the healthy controls. Abbreviations as in Table 4

When comparing the subgroup with HSCT for malignancies with controls, the higher E/Ea persisted (10.33 ± 3.55 vs. 7.24 ± 2.29, P = 0.013). Additionally, LVEF (52 ± 4 vs. 58 ± 5%, P = 0.008), GLS (− 19.8 ± 3.3 vs. -22.9 ± 3.0%, P = 0.022) and GCSR-S (-1.13 ± 0.15 vs. -1.32 ± 0.26 1/s, P = 0.047) were lower. For the HSCT recipients treated for malignancies, aortic PWV was: 4.39 ± 0.18 m/s (P = 0.36 vs. controls), LV concentricity: 0.63 ± 0.10 g/mL (P = 0.73), native T1: 1208 ± 36 ms (P = 0.138) and myocardial triglyceride content: 0.53 ± 0.21 (P = 0.69).

Post-hoc power analysis

We had 0.90 power to detect a statistically significant difference if the HSCT population would have had an aortic PWV of 4.64 m/s or higher (+ 0.33 m/s higher as compared to the healthy controls), myocardial triglyceride content of 0.69% or higher (+ 0.19%) and native T1 of 1268 ms or higher (+ 39 ms). These differences in aortic PWV, myocardial triglyceride content and native T1 are comparable to approximately 3, 20 and 35 years of normal aging, respectively [23,24,25].


Our study showed that young adults with a history of pediatric HSCT, as compared to healthy controls in the same age range, have subclinical impairments in LV diastolic function and tend to have a lower LV systolic function, whereas aortic stiffness and myocardial tissue characteristics are comparable. Our results indicate that CMR-derived LV diastolic parameters, rather than aortic PWV, cardiac native T1 or myocardial triglyceride content, may be early markers of cardiovascular deterioration after HSCT.

LV function

The estimated LV filling pressure was found to be increased after childhood HSCT, which is indicative of reduced LV diastolic function. When LV myocardial relaxation progressively deteriorates, LV filling pressure becomes increased to compensate for the impaired LV diastolic filling [26]. Accordingly, in the patients who received HSCT during childhood, the E/A ratio was comparable as in the controls, whereas LV filling pressure was higher. In contrast to the impairments in LV diastolic function, the decrease in LV systolic function in the HSCT group compared to the control group was non-significant. This observation is in keeping with prior longitudinal studies using echocardiography, which showed that diastolic dysfunction precedes systolic dysfunction in response to cardiotoxic exposures [27, 28]. Anthracycline-related cardiomyopathy is considered to result from free radical formation, the consequent mitochondrial dysfunction and myofibrillar disarray, and eventual myocyte necrosis; in turn, with the ongoing remodeling of the injured heart, functional impairments may progressively develop [2]. Furthermore, irradiation may cause microcirculatory damage, leading to myocardial ischemia and cell death [2]. However, despite the reduced LV diastolic function, there was no interstitial fibrosis based on native T1 mapping. Furthermore, aortic stiffening or myocardial steatosis did not appear contributing factors to the abnormalities in LV diastolic function. Hence, when using CMR, LV diastolic parameters in particular may be important in the monitoring of HSCT-related cardiovascular deteriorations.

In our study, the trend for lower LVEF and significantly lower cardiac index in the pediatric HSCT recipients compared to the controls seemed to be related to impairments in GLS rather than GCS. GLS and GCS, respectively, can be interpreted as measures of the contractility of the cardiomyocytes in oblique orientation in the subendocardium and those in circumferential arrangement in the mesomyocardium [29, 30]. Accordingly, our results may suggest that the HSCT-related therapies predominately affect the subendocardial myocardium, and not the midwall. In the subgroup with a history of a malignant disorder, longitudinal and circumferential systolic strain-derived parameters were lower as compared to healthy controls; this is consistent with previous studies which assessed LV systolic strain parameters in patients receiving anthracyclines [31, 32].

We hypothesized that aortic stiffening or altered myocardial tissue characteristics would be early markers of cardiovascular disturbances after HSCT, preceding manifest LV dysfunction. Of interest, as chemotherapy is considered to disrupt mitochondrial function and irradiation has been demonstrated to induce microcirculatory damage [2], several other CMR-derived parameters may have potential for the early detection of cardiac effects of HSCT, for example phosphorus CMRS (31P-CMRS) for the assessment of myocardial energetics and perfusion imaging [33].

Aortic stiffness

No differences were observed in aortic stiffness between the HSCT study population and the healthy controls. Based on previous studies, we expected to find a significantly higher aortic stiffness, at least in the patients with HSCT for a malignant disorder, because of the pre-transplant cardiotoxic exposures. Anthracycline and/or chest irradiation in the treatment of the primary cancer may cause endothelial damage due to the generation of reactive oxygen species [34, 35]. Also, the conditioning in allogeneic HSCT, including total body irradiation and/or chemotherapy, is considered to initiate disruption of the endothelium [36]. Subsequently, the structural changes in the vascular matrix and the disturbed endothelial function may increase the vascular tone of the arterial wall [2]. In addition, metabolic disturbances caused by the immune dysregulation and iron overload are well recognized adverse effects of HSCT [12, 37]. There is extensive evidence for an elevated aortic stiffness due to hypertension, dyslipidemia and hyperglycemia in the general population [14]. The role of iron overload in endothelial dysfunction is less evident, although iron chelation in coronary artery disease has been shown to improve endothelium-dependent vasodilation [38].

Previous longitudinal studies in patients with breast cancer, lymphoma or leukemia showed that aortic PWV increases by approximately 1.5 to 2-fold upon anthracyclines in the first 4 to 6 months, also for low or moderate dosages, compared to the pre-treatment measurements [39, 40]. Another study measured a higher carotid artery stiffness in children who had received allogeneic HSCT, for malignant or non-malignant disease, even though this study population showed only minor, subclinical metabolic derangements, and no other cardiovascular impairments were detected [41]. However, the increased aortic stiffness due to cardiotoxic exposures may decrease after the cardiotoxic therapy has been discontinued, as was shown in breast cancer patients who were followed before and 1, 4 and 14 months after anthracycline or trastuzumab chemotherapy [42]. We measured aortic PWV after a mean time of 14.8 ± 5.0 years after HSCT. Possibly, there may have been changes in aortic PWV acutely after pediatric HSCT, which may have normalized after several years.

Myocardial tissue characteristics

There were no differences in myocardial tissue characteristics and LV concentricity between the HSCT study population and the healthy controls, which otherwise would have been indicative of LV remodeling. To our knowledge, this is the first study in which native T1, as a measure of interstitial fibrosis, and myocardial triglyceride content, as a measure of metabolic remodeling, have been assessed in young adults with prior HSCT. It should be noted that diffuse fibrosis can be estimated based on calculating native T1 or the extracellular volume fraction (ECV); the latter requires CMR contrast administration [17]. Recent CMR studies have demonstrated that both native T1 and ECV are increased in middle-aged patients with prior anthracycline therapy [43, 44]. Other studies in adolescents exposed to anthracyclines, however, reported native T1 and/or ECV values within the normal range [45, 46]; nonetheless, a correlation was found between the cumulative anthracycline dose and native T1 and ECV [46]. In contrast to these observations suggestive of diffuse fibrosis, focal fibrosis, visualized with late gadolinium enhancement CMR, is uncommon after anthracycline exposure, even when systolic function is subnormal [44,45,46,47]. In our HSCT study population as compared to the controls, we did not find a higher native T1, even though a lower LVEF was measured. Also, when analyzing the patients who received HSCT for a malignant indication separately, no difference with the control group was found. Therefore, in contrast to middle-aged patients with prior anthracycline therapy as described in other studies, we may conclude that there is no myocardial fibrosis, at least not substantial, in young adults with a history of childhood HSCT as compared to healthy controls.

From previous studies it is known that the metabolic syndrome, presumably due to changes in myocardial fatty acid uptake and/or oxidation, is associated with myocardial steatosis [18, 48]. Although allografting increases the risk of early development of the metabolic syndrome [12], our HSCT study population aged 18–25 years had on average no obesity, impaired fasting glucose or hypertriglyceridemia. Possibly, myocardial steatosis may arise after HSCT at middle-age, when metabolic disorders may have developed, while myocardial triglyceride content remains within the normal ranges at younger age. Based on our findings, we may conclude that myocardial triglyceride content does not represent an early marker of cardiometabolic disease after pediatric HSCT in the young adult population.


Strength of this study is the comprehensive CMR evaluation, including measures of both LV systolic and diastolic function, vascular function and myocardial tissue characteristics. CMR rather than ultrasound is the most accurate non-invasive modality to quantify aortic stiffness [49]. Diffuse fibrosis was measured using native T1 mapping, which is a non-contrast technique; contrast-enhanced CMR can be considered less suitable for screening programs. Cardiac 1H-CMRS is the gold standard for non-invasive assessment of myocardial steatosis [48]. In clinical practice, echocardiography is commonly used for follow-up after cardiotoxic therapies [11]. Compared to echocardiography, little is known regarding the value of CMR for cardiovascular screening after HSCT.


Our study did not comprise CMR before and shortly after HSCT, but had a cross-sectional design because of ethical considerations. Therefore, we cannot exclude temporary deteriorations in, for example, aortic stiffness. Our study population was too small to evaluate the patient group with HSCT for non-malignant disease separately. The observed abnormalities in LV diastolic function and the non-significantly lower LV systolic function seemed to be driven by the patients who had received HSCT for a malignant disease. However, based on our study we cannot rule out late LV functional impairment after HSCT for non-malignant disease. Also, we were not able to assess the correlations between the different HSCT conditioning regimens, graft-versus-host-disease or iron overload and cardiovascular measures. However, this study was aimed at identifying CMR parameters which may show deteriorations before overt cardiovascular disease develops, and not to identify HSCT-related factors which may be helpful in late cardiovascular disease risk stratification. Echocardiography was not applied in the present study, which would have allowed for comparison of CMR with echocardiography for the early detection of cardiac effects of HSCT.

As a measure of interstitial fibrosis, ECV may be preferable. The native T1 value may be affected by several processes other than edema or increased myocardial collagen in the extracellular compartment [17]. Hence, based on the native T1 measurements, we may not be able to exclude the presence of diffuse fibrosis in post-HSCT patients. However, as all patients were clinically monitored for cardiac iron deposition and as myocardial steatosis was ruled out based on the 1H-MRS measurements, it is unlikely that fibrosis related increases in native T1 might have been cancelled out by myocardial iron or fat related native T1 decreases.


In young adults who received HSCT during childhood, LV diastolic function was decreased (higher estimated LV filling pressure) and LV systolic function (LVEF and GLS) tended to be reduced as compared to healthy controls, whereas no concomitant differences were found in aortic stiffness (aortic PWV) and myocardial tissue characteristics (native T1 and myocardial triglyceride content). Therefore, when using CMR, the assessment of LV diastolic function in particular is important for early detection of patients at risk of HSCT-related cardiovascular disease. Further research, including longitudinal CMR measurements, is needed to show the predictive value of subclinical LV functional abnormalities for the development of symptomatic cardiovascular disease after HSCT. Also, comparative studies of imaging modalities should reveal whether CMR-derived LV function has additive value for screening programs when performed next to the current echocardiography-based follow-up of the cardiovascular late effects of HSCT. For example, CMR may be used for the selection of patients who require frequent follow-up by standard echocardiography.



Proton cardiovascular magnetic resonance spectroscopy


Body mass index


Body surface area


Balanced steady-state free precession


Cardiovascular magnetic resonance

E/A ratio:

Ratio of the transmitral early and late peak filling rate


Estimated left ventricular filling pressure


Early diastolic mitral septal tissue velocity




Extracellular volume


Ejection fraction




Global circumferential strain


Global circumferential early peak diastolic strain rate


Global circumferential peak systolic strain rate


Global longitudinal strain


Global longitudinal early peak diastolic strain rate


Global longitudinal peak systolic strain rate


Hematopoietic stem cell transplantation


Left ventricle/left ventricular


Left ventricular ejection fraction


Modified Look-Locker inversion recovery


Pulse wave velocity


SENSitivity Encoding


Echo time


Repetition time




  1. Uderzo C, Pillon M, Corti P, Tridello G, Tana F, Zintl F, et al. Impact of cumulative anthracycline dose, preparative regimen and chronic graft-versus-host disease on pulmonary and cardiac function in children 5 years after allogeneic hematopoietic stem cell transplantation: a prospective evaluation on behalf of the EBMT pediatric diseases and late effects working parties. Bone Marrow Transplant. 2007;39(11):667–75.

    Article  CAS  Google Scholar 

  2. Lipshultz SE, Adams MJ, Colan SD, Constine LS, Herman EH, Hsu DT, et al. Long-term cardiovascular toxicity in children, adolescents, and young adults who receive cancer therapy: pathophysiology, course, monitoring, management, prevention, and research directions: a scientific statement from the American Heart Association. Circulation. 2013;128(17):1927–95.

    Article  Google Scholar 

  3. Pulte D, Gondos A, Brenner H. Trends in 5- and 10-year survival after diagnosis with childhood hematologic malignancies in the United States, 1990-2004. J Natl Cancer Inst. 2008;100(18):1301–9.

    Article  Google Scholar 

  4. Chima RS, Daniels RC, Kim MO, Li D, Wheeler DS, Davies SM, et al. Improved outcomes for stem cell transplant recipients requiring pediatric intensive care. Pediatr Crit Care Med. 2012;13(6):e336–42.

    Article  Google Scholar 

  5. Niederwieser D, Baldomero H, Szer J, Gratwohl M, Aljurf M, Atsuta Y, et al. Hematopoietic stem cell transplantation activity worldwide in 2012 and a SWOT analysis of the worldwide network for blood and marrow transplantation group including the global survey. Bone Marrow Transplant. 2016;51(6):778–85.

    Article  CAS  Google Scholar 

  6. Battiwalla M, Tichelli A, Majhail NS. Long-term survivorship after hematopoietic cell transplantation: roadmap for research and care. Biol Blood Marrow Transplant. 2017;23(2):184–92.

    Article  Google Scholar 

  7. Battiwalla M, Hashmi S, Majhail N, Pavletic S, Savani BN, Shelburne N. National Institutes of Health hematopoietic cell transplantation late effects initiative: developing recommendations to improve survivorship and long-term outcomes. Biol Blood Marrow Transplant. 2017;23(1):6–9.

    Article  Google Scholar 

  8. Armenian SH, Chemaitilly W, Chen M, Chow EJ, Duncan CN, Jones LW, et al. National Institutes of Health hematopoietic cell transplantation late effects initiative: the cardiovascular disease and associated risk factors working group report. Biol Blood Marrow Transplant. 2017;23(2):201–10.

    Article  Google Scholar 

  9. Pulsipher MA, Skinner R, McDonald GB, Hingorani S, Armenian SH, Cooke KR, et al. National Cancer Institute, National Heart, lung and blood institute/pediatric blood and marrow transplantation consortium first international consensus conference on late effects after pediatric hematopoietic cell transplantation: the need for pediatric-specific long-term follow-up guidelines. Biol Blood Marrow Transplant. 2012;18(3):334–47.

    Article  Google Scholar 

  10. Chow EJ, Anderson L, Baker KS, Bhatia S, Guilcher GM, Huang JT, et al. Late effects surveillance recommendations among survivors of childhood hematopoietic cell transplantation: a Children’s oncology group report. Biol Blood Marrow Transplant. 2016;22(5):782–95.

    Article  Google Scholar 

  11. Shankar SM, Marina N, Hudson MM, Hodgson DC, Adams MJ, Landier W, et al. Monitoring for cardiovascular disease in survivors of childhood cancer: report from the cardiovascular disease task force of the Children’s oncology group. Pediatrics. 2008;121(2):e387–96.

    Article  Google Scholar 

  12. Taskinen M, Saarinen-Pihkala UM, Hovi L, Lipsanen-Nyman M. Impaired glucose tolerance and dyslipidaemia as late effects after bone-marrow transplantation in childhood. Lancet. 2000;356(9234):993–7.

    Article  CAS  Google Scholar 

  13. McEniery CM, Yasmin HIR, Qasem A, Wilkinson IB, Cockcroft JR, et al. Normal vascular aging: differential effects on wave reflection and aortic pulse wave velocity: the Anglo-Cardiff Collaborative Trial (ACCT). J Am Coll Cardiol. 2005;46(9):1753–60.

    Article  Google Scholar 

  14. Terentes-Printzios D, Vlachopoulos C, Xaplanteris P, Ioakeimidis N, Aznaouridis K, Baou K, et al. Cardiovascular risk factors accelerate progression of vascular aging in the general population: results from the CRAVE study (cardiovascular risk factors affecting vascular age). Hypertension. 2017;70(5):1057–64.

    Article  CAS  Google Scholar 

  15. Cheung YF, Chan GC, Ha SY. Arterial stiffness and endothelial function in patients with beta-thalassemia major. Circulation. 2002;106(20):2561–6.

    Article  CAS  Google Scholar 

  16. Ben-Shlomo Y, Spears M, Boustred C, May M, Anderson SG, Benjamin EJ, et al. Aortic pulse wave velocity improves cardiovascular event prediction: an individual participant meta-analysis of prospective observational data from 17,635 subjects. J Am Coll Cardiol. 2014;63(7):636–46.

    Article  Google Scholar 

  17. Puntmann VO, Peker E, Chandrashekhar Y, Nagel E. T1 mapping in characterizing myocardial disease: a comprehensive review. Circ Res. 2016;119(2):277–99.

    Article  CAS  Google Scholar 

  18. Rayner JJ, Banerjee R, Holloway CJ, Lewis AJM, Peterzan MA, Francis JM, et al. The relative contribution of metabolic and structural abnormalities to diastolic dysfunction in obesity. Int J Obes. 2018;42(3):441-7.

  19. de Heer P, Bizino MB, Versluis MJ, Webb AG, Lamb HJ. Improved cardiac proton magnetic resonance spectroscopy at 3 T using high permittivity pads. Investig Radiol. 2016;51(2):134–8.

    Article  Google Scholar 

  20. Grotenhuis HB, Westenberg JJ, Steendijk P, van der Geest RJ, Ottenkamp J, Bax JJ, et al. Validation and reproducibility of aortic pulse wave velocity as assessed with velocity-encoded MRI. J Magn Reson Imaging. 2009;30(3):521–6.

    Article  Google Scholar 

  21. Brandts A, Bertini M, van Dijk EJ, Delgado V, Marsan NA, van der Geest RJ, et al. Left ventricular diastolic function assessment from three-dimensional three-directional velocity-encoded MRI with retrospective valve tracking. J Magn Reson Imaging. 2011;33(2):312–9.

    Article  Google Scholar 

  22. de Heer P, Bizino MB, Lamb HJ, Webb AG. Parameter optimization for reproducible cardiac (1) H-MR spectroscopy at 3 tesla. J Magn Reson Imaging. 2016;44(5):1151–8.

    Article  Google Scholar 

  23. Ohyama Y, Teixido-Tura G, Ambale-Venkatesh B, Noda C, Chugh AR, Liu CY, et al. Ten-year longitudinal change in aortic stiffness assessed by cardiac MRI in the second half of the human lifespan: the multi-ethnic study of atherosclerosis. Eur Heart J Cardiovasc Imaging. 2016;17(9):1044–53.

    Article  Google Scholar 

  24. van der Meer RW, Rijzewijk LJ, Diamant M, Hammer S, Schar M, Bax JJ, et al. The ageing male heart: myocardial triglyceride content as independent predictor of diastolic function. Eur Heart J. 2008;29(12):1516–22.

    Article  Google Scholar 

  25. Roy C, Slimani A, de Meester C, Amzulescu M, Pasquet A, Vancraeynest D, et al. Age and sex corrected normal reference values of T1, T2 T2* and ECV in healthy subjects at 3T CMR. J Cardiovasc Magn Reson. 2017;19(1):72.

    Article  Google Scholar 

  26. Zile MR, Baicu CF, Gaasch WH. Diastolic heart failure--abnormalities in active relaxation and passive stiffness of the left ventricle. N Engl J Med. 2004;350(19):1953–9.

    Article  CAS  Google Scholar 

  27. Tassan-Mangina S, Codorean D, Metivier M, Costa B, Himberlin C, Jouannaud C, et al. Tissue Doppler imaging and conventional echocardiography after anthracycline treatment in adults: early and late alterations of left ventricular function during a prospective study. Eur J Echocardiogr. 2006;7(2):141–6.

    Article  Google Scholar 

  28. Nagy AC, Cserep Z, Tolnay E, Nagykalnai T, Forster T. Early diagnosis of chemotherapy-induced cardiomyopathy: a prospective tissue Doppler imaging study. Pathol Oncol Res. 2008;14(1):69–77.

    Article  Google Scholar 

  29. Ingels NB Jr. Myocardial fiber architecture and left ventricular function. Technol Health Care. 1997;5(1–2):45–52.

    PubMed  Google Scholar 

  30. Dorri F, Niederer PF, Lunkenheimer PP, Anderson RH. The architecture of the left ventricular myocytes relative to left ventricular systolic function. Eur J Cardiothorac Surg. 2010;37(2):384–92.

    PubMed  Google Scholar 

  31. Ong G, Brezden-Masley C, Dhir V, Deva DP, Chan KKW, Chow CM, et al. Myocardial strain imaging by cardiac magnetic resonance for detection of subclinical myocardial dysfunction in breast cancer patients receiving trastuzumab and chemotherapy. Int J Cardiol. 2018.

  32. Lunning MA, Kutty S, Rome ET, Li L, Padiyath A, Loberiza F, et al. Cardiac magnetic resonance imaging for the assessment of the myocardium after doxorubicin-based chemotherapy. Am J Clin Oncol. 2015;38(4):377–81.

    Article  CAS  Google Scholar 

  33. Hudsmith LE, Neubauer S. Magnetic resonance spectroscopy in myocardial disease. JACC Cardiovasc Imaging. 2009;2(1):87–96.

    Article  Google Scholar 

  34. Keltai K, Cervenak L, Mako V, Doleschall Z, Zsary A, Karadi I. Doxorubicin selectively suppresses mRNA expression and production of endothelin-1 in endothelial cells. Vasc Pharmacol. 2010;53(5–6):209–14.

    Article  CAS  Google Scholar 

  35. Soucy KG, Attarzadeh DO, Ramachandran R, Soucy PA, Romer LH, Shoukas AA, et al. Single exposure to radiation produces early anti-angiogenic effects in mouse aorta. Radiat Environ Biophys. 2010;49(3):397–404.

    Article  CAS  Google Scholar 

  36. Woywodt A, Scheer J, Hambach L, Buchholz S, Ganser A, Haller H, et al. Circulating endothelial cells as a marker of endothelial damage in allogeneic hematopoietic stem cell transplantation. Blood. 2004;103(9):3603–5.

    Article  CAS  Google Scholar 

  37. de Witte T. The role of iron in patients after bone marrow transplantation. Blood Rev. 2008;22(Suppl 2):S22–8.

    Article  Google Scholar 

  38. Duffy SJ, Biegelsen ES, Holbrook M, Russell JD, Gokce N, Keaney JF Jr, et al. Iron chelation improves endothelial function in patients with coronary artery disease. Circulation. 2001;103(23):2799–804.

    Article  CAS  Google Scholar 

  39. Chaosuwannakit N, D'Agostino R Jr, Hamilton CA, Lane KS, Ntim WO, Lawrence J, et al. Aortic stiffness increases upon receipt of anthracycline chemotherapy. J Clin Oncol. 2010;28(1):166–72.

    Article  CAS  Google Scholar 

  40. Drafts BC, Twomley KM, D'Agostino R Jr, Lawrence J, Avis N, Ellis LR, et al. Low to moderate dose anthracycline-based chemotherapy is associated with early noninvasive imaging evidence of subclinical cardiovascular disease. JACC Cardiovasc Imaging. 2013;6(8):877–85.

    Article  Google Scholar 

  41. Turanlahti MI, Taskinen M, Saarinen-Pihkala U, Jokinen EV. Time-related arterial changes after allogeneic hematopoietic stem cell transplantation in children. Pediatr Res. 2013;73(6):777–82.

    Article  Google Scholar 

  42. Grover S, Lou PW, Bradbrook C, Cheong K, Kotasek D, Leong DP, et al. Early and late changes in markers of aortic stiffness with breast cancer therapy. Intern Med J. 2015;45(2):140–7.

    Article  CAS  Google Scholar 

  43. Jordan JH, Vasu S, Morgan TM, D'Agostino RB Jr, Melendez GC, Hamilton CA, et al. Anthracycline-associated T1 mapping characteristics are elevated independent of the presence of cardiovascular comorbidities in Cancer survivors. Circ Cardiovasc Imaging. 2016;9(8).

  44. Neilan TG, Coelho-Filho OR, Shah RV, Feng JH, Pena-Herrera D, Mandry D, et al. Myocardial extracellular volume by cardiac magnetic resonance imaging in patients treated with anthracycline-based chemotherapy. Am J Cardiol. 2013;111(5):717–22.

    Article  CAS  Google Scholar 

  45. Toro-Salazar OH, Gillan E, O'Loughlin MT, Burke GS, Ferranti J, Stainsby J, et al. Occult cardiotoxicity in childhood cancer survivors exposed to anthracycline therapy. Circ Cardiovasc Imaging. 2013;6(6):873–80.

    Article  Google Scholar 

  46. Tham EB, Haykowsky MJ, Chow K, Spavor M, Kaneko S, Khoo NS, et al. Diffuse myocardial fibrosis by T1-mapping in children with subclinical anthracycline cardiotoxicity: relationship to exercise capacity, cumulative dose and remodeling. J Cardiovasc Magn Reson. 2013;15:48.

    Article  Google Scholar 

  47. Ylanen K, Poutanen T, Savikurki-Heikkila P, Rinta-Kiikka I, Eerola A, Vettenranta K. Cardiac magnetic resonance imaging in the evaluation of the late effects of anthracyclines among long-term survivors of childhood cancer. J Am Coll Cardiol. 2013;61(14):1539–47.

    Article  Google Scholar 

  48. van der Meer RW, Lamb HJ, Smit JW, de Roos A. MR imaging evaluation of cardiovascular risk in metabolic syndrome. Radiology. 2012;264(1):21–37.

    Article  Google Scholar 

  49. Cavalcante JL, Lima JA, Redheuil A, Al-Mallah MH. Aortic stiffness: current understanding and future directions. J Am Coll Cardiol. 2011;57(14):1511–22.

    Article  Google Scholar 

Download references


Not applicable.


Not applicable.

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



EHMP drafted the manuscript. AAWR, HJL and ML were major contributors for writing the manuscript. AAWR and ACL designed the study. JJMW and HJL developed the scanning protocol. ML, DB and EHMP were responsible for patient and volunteer recruitment. EHMP acquired the data. EHMP, RJvdG, QT, JJMW and HJL were involved in the data analysis. All authors critically reviewed and approved the manuscript.

Corresponding author

Correspondence to Elisabeth H. M. Paiman.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the local institutional review board (Leiden University Medical Center, the Netherlands) and complies with the Declaration of Helsinki. All participants provided written informed consent prior to enrolment into the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Paiman, E.H.M., Louwerens, M., Bresters, D. et al. Late effects of pediatric hematopoietic stem cell transplantation on left ventricular function, aortic stiffness and myocardial tissue characteristics. J Cardiovasc Magn Reson 21, 6 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: