Background
Sickle cell disease (SCD) is a common inherited hemoglobinopathy. Adults with SCD manifest both systolic and diastolic cardiac dysfunction, though the age of onset of dysfunction has not been defined. Left ventricular (LV) rotational mechanics have not been studied in children with SCD. The aim of this study was to investigate whether cardiac rotational mechanics differed between children with SCD and age-matched controls.
Methods
Basal and apical LV short-axis images were acquired prospectively in 213 patients with SCD (mean age, 14.1 ± 2.6 years) and 49 controls (mean age, 13.3 ± 2.8 years) from the Muhimbili Sickle Cohort in Dar es Salaam, Tanzania. The magnitude of basal and apical rotation, net twist angle, torsion, and untwist rate were obtained by two-dimensional speckle-tracking. The timing of events was normalized to aortic valve closure.
Results
Mean basal rotation was significantly lower in patients with SCD compared with controls ( P = .012), although no difference was observed in apical rotation ( P = .37). No statistically significant differences in torsion or net twist angle were detected. Rotation rate at the apex ( P = .001) and base ( P = .0004) were significantly slower in subjects with SCD compared with controls. Mean peak untwisting rate was also significantly slower in patients with SCD ( P = .006). No associations were found between hemoglobin concentration and apical rotation, basal rotation, net twist, and torsion.
Conclusion
This study demonstrates alterations in LV rotational mechanics in children with SCD, including lower basal rotation, peak differential twist, and untwist rate. These abnormalities denote subclinical changes in LV systolic and diastolic performance in children with SCD. Future work may reveal an association between rotational metrics and long-term patient outcomes.
Sickle cell disease (SCD) is a common hemoglobinopathy, the sequelae of which can be life threatening. Patients with SCD are at risk for left ventricular (LV) systolic dysfunction. Previous studies have also demonstrated that this population is also susceptible to diastolic dysfunction, and some suggest that this diastolic dysfunction plays a role in the development of pulmonary hypertension. One proposed mechanism for the development of cardiac dysfunction is microvascular occlusion in the coronary circulation as a result of red blood cell sickling, resulting in chronic, indolent myocardial ischemia. In addition, chronic anemia leads to a high cardiac output state and LV volume overload. Over time, this may cause LV dilation, which may also adversely affect LV function.
Conventional echocardiographic measures of LV function, such as shortening fraction and ejection fraction, have not been demonstrably different when patients with SCD and unaffected individuals have been compared ; however, decreased LV peak systolic strain in subjects with SCD suggests that conventional methods of evaluating LV function by measuring chamber size, such as shortening fraction, are not precise enough to detect early dysfunction. LV rotation, net twist angle (twist), and torsion, a means of quantifying the rotational movement of the myocardium during systole, have been studied in other disease states and have demonstrated an early decline in LV systolic performance. LV torsion has not been previously studied in pediatric patients with SCD. We hypothesized that there would be significant differences in rotational metrics between individuals with SCD and healthy age-matched controls. The aim of this study was to evaluate cardiac rotational mechanics in patients with SCD, compared with age-matched, unaffected control subjects.
Methods
A nested case-control study was undertaken, including study participants with SCD and unaffected control patients from the Muhimbili Sickle Cohort (MSC) in Dar es Salaam, Tanzania.
Study Population
Cases were recruited prospectively from patients enrolled in the MSC and controls from children who had attended the clinic for sickle screening and were determined not to have SCD. Inclusion criteria included age between 9 and 19 years and, for patients with SCD, sickle status confirmed by hemoglobin (Hb) electrophoresis and high-performance liquid chromatography, while controls without SCD had the HbAS or HbAA phenotype by electrophoresis. Blood pressure was measured twice after 5 min of rest using an automated digital device (PRO 400 V2 Dinamap; GE Healthcare, Little Chalfont, United Kingdom) and the mean calculated. Axial temperature, height, and weight were recorded, and a clinical examination was performed to exclude acute complications. Exclusion criteria included history of blood transfusion within the previous 2 months, symptoms of sickle crisis within the previous 2 weeks, febrile illness within the 7 days before echocardiography, fever or signs of acute illness on the day of echocardiography, and known or echocardiographically identified hemodynamically significant congenital heart disease.
The study was approved in the institutional review board of Children’s Hospital Colorado (Reference No. COMIRB protocol #10-0030) and by the Muhimbili University of Health & Allied Sciences Research and Publications Committee (Reference No. MU/RP/AEC/Vol.XIIV/01). Written informed consent was obtained from all study participants.
Definitions
Rotation was defined as circumferential myocardial motion around the centroid, or long axis, of the ventricle. By convention, rotation is said to be positive in the counterclockwise direction when viewed from the apex of the left ventricle and vice versa. As depicted in Figure 1 , a basal short-axis plane used to determine twist was obtained just below the mitral valve annulus, where muscle appears in the plane of insonation throughout the cardiac cycle. The apical short-axis plane used was defined as a short-axis plane just above the apex, where the LV lumen can be seen throughout the cardiac cycle. Net twist angle was defined as the instantaneous difference in rotation between the apex and the base at the time of aortic valve closure, or end-systole. Peak differential twist is the maximum difference in rotation between the apex and base, independent of time. Torsion was defined as net twist divided by the length of the left ventricle at end-diastole to account for heart size. Untwisting is defined as the difference between apical recoil and basal recoil. The timing of these events is described both in milliseconds and as a percentage of systole, with the onset of systole defined as the initial deflection of the Q wave on the surface electrocardiogram and end-systole defined as aortic valve closure.
Echocardiography
For each patient, measurements of LV chamber size and fractional shortening were made from short-axis basal images. From the apical window, LV peak systolic mitral annular velocity (LV S′), LV peak early diastolic mitral annular velocity (LV E′), early diastolic mitral inflow (LV E), duration of late diastolic mitral inflow (mitral valve a-wave duration), and duration of the late diastolic a wave in pulmonary vein inflow (pulmonary vein a-wave duration) were made. The LV myocardial performance index was also calculated from mitral inflow and aortic outflow spectral Doppler patterns.
Short-axis images of the base and apex of the left ventricle were obtained using the VividQ platform (GE Healthcare) and an M4S, 5S, or 7S probe, depending on patient characteristics and image quality. Frame rates of >75 and <125 frames/sec were accepted. An acquisition of five consecutive beats was obtained. An apical view was also obtained for the purposes of measuring end-diastolic LV length and performing spectral Doppler measurement of mitral inflow and aortic outflow. End-diastolic LV length was measured from the apical four-chamber view as the distance from the mitral valve to the apex of the left ventricle.
Commercially available software (EchoPAC version 113, revision 0.4; GE Healthcare) was used to perform speckle-tracking analysis. A single cardiac cycle was isolated for the base, and a region of interest was identified. The magnitude and timing of basal rotation was calculated. The process was repeated for the short-axis apical image. Net twist was calculated as the difference between apical and basal rotation. Torsion was calculated by dividing the net twist by LV end-diastolic length. This process was repeated two additional times and an average over three cardiac cycles was calculated.
Laboratory Data
Hb phenotype was determined by alkaline Hb electrophoresis (Helena; Sunderland, United Kingdom) and by quantification of Hb fractions, including HbF, by high-performance liquid chromatography (β-Thalassemia Short Program, VARIANT analyzer; BioRad, Hercules, CA). In subjects with SCD, the most recent steady-state Hb, total bilirubin, conjugated bilirubin, and lactate dehydrogenase levels were obtained from MSC records (steady state defined as absence of pain or fever, malaria rapid test negative, and with no recorded hospital admission or blood transfusion within 30 days). Full blood counts were performed using an automated cell counter (Pentra 60; Horiba ABX, Kyoto, Japan). Biochemical tests were performed using an automated chemistry analyzer (Cobas Mira [Roche, New York, NY] or Abbott Architect [Abbott Diagnostics, New York, NY]).
Statistical Analysis
SAS version 9.2 (SAS Institute Inc, Cary, NC) was used to perform statistical analyses. Continuous variables are presented as mean ± SD or median (interquartile range), and categorical variables are presented as proportions. The Student independent two-sample t test was used to compare mean values between two patient groups. Simple linear regression was used to evaluate associations between independent variables and a continuous, dependent variable. Pearson correlation was used to evaluate associations between measures of rotational mechanics and Hb or laboratory measures of hemolysis in the children with SCD. A P value < .05 was determined to be the threshold for statistical significance. Analyses were compared when the five children with HbSB+ were included or excluded from the patients with SCD and no material differences were observed, and therefore all reported analyses include both the HbSS and HbSB+ cases.
Interobserver and Intraobserver Variability
Eleven consecutive patients were evaluated separately by two independent observers (H.H.H. and G.A.), who evaluated net twist, torsion, twist rate, and peak untwisting rate in a blinded fashion. To reflect the methodology used in this study, each measurement was made on what the investigator felt was the best cardiac cycle among the five that were acquired. A single observer (H.H.H.) performed serial analysis of the same 11 patients 2 months later using the same methodology to determine intraobserver variability. Analysis included calculation of intraclass correlation coefficients of variation.
Results
Three hundred subjects were recruited into the study. On the day of the echocardiographic assessment, 22 subjects were excluded before echocardiography because they were unwell, while seven subjects were excluded after echocardiography because of minor structural cardiac abnormalities. Four additional subjects were excluded after speckle-tracking analysis because of poor image quality. Thus a total of 262 participants were included in the present analysis, consisting of 213 patients and 49 unaffected controls (see Figure 2 ).
Baseline demographic and conventional echocardiographic data are summarized in Table 1 . There was no significant difference in age or gender composition between the two groups. Patients with SCD had higher heart rates ( P = .0042), lower mean Hb levels ( P < .0001), and lower systolic ( P < .0001) and diastolic ( P = .0002) blood pressures. The SCD group also had lower body mass index Z scores ( P < .0001) and lower height-for-age Z scores ( P < .0001). The SCD group had significantly greater LV dilation in both the long- and short-axis dimensions, as assessed by LV end-diastolic dimension Z score and LV length ( P < .0001 for both comparisons). The LV sphericity index was significantly lower in the SCD group ( P = .006), indicating that the left ventricle more closely approximated a sphere, which has a sphericity index of 1. No difference was observed in fractional shortening ( P = .16) or in tissue Doppler measurements of LV function. However, the LV myocardial performance index was significantly different in the SCD group (0.336 ± 0.12 vs 0.368 ± 0.12, P = .03).
Demographic variable | Controls ( n = 49) | Patients with SCD ( n = 218) | P | ||
---|---|---|---|---|---|
n | Mean ± SD or count (%) | n | Mean ± SD or count (%) | ||
Age (y) | 49 | 13.3 ± 2.8 | 218 | 14.1 ± 2.6 | .07 |
Male gender | 49 | 26 (52%) | 218 | 122 (55%) | .20 |
Heart rate (beats/min) | 49 | 80.5 ± 16.4 | 218 | 86.1 ± 11.1 | .0042 |
Hb (g/dL) | 49 | 11.6 ± 1.7 | 218 | 7.2 ± 1.1 | <.0001 |
Systolic blood pressure (mm Hg) | 49 | 116 ± 13.2 | 218 | 109 ± 10.2 | <.0001 |
Diastolic blood pressure (mm Hg) | 49 | 68 ± 8.4 | 218 | 64 ± 6.0 | .0002 |
Body mass index Z score | 49 | −0.57 ± 1.29 | 218 | −1.5 ± 1.2 | <.0001 |
Height-for-age Z score | 49 | −0.96 ± 1.3 | 218 | −2.26 ± 1.14 | <.0001 |
LV dimension | |||||
LV end-diastolic diameter Z score | 49 | -0.42 ± 1.1 | 217 | 2.25 ± 1.4 | <.0001 |
LV end-diastolic length (cm) | 49 | 6.7 ± 0.7 | 217 | 7.3 ± 0.6 | <.0001 |
LV sphericity index | 49 | 1.62 ± 0.2 | 216 | 1.53 ± 0.1 | .0006 |
LV function | |||||
Fractional shortening (%) | 49 | 37.91 ± 4.67 | 217 | 38.14 ± 4.26 | .16 |
LV myocardial performance index | 47 | 0.34 ± 0.11 | 205 | 0.37 ± 0.12 | .03 |
LV S′ velocity (cm/sec) | 46 | 10.18 ± 2.08 | 213 | 10.01 ± 1.84 | .39 |
LV E′ velocity (cm/sec) | 47 | 19.86 ± 3.48 | 210 | 19.93 ± 3.38 | .50 |
Mitral E-wave velocity (m/sec) | 49 | 1.00 ± 0.17 | 216 | 1.07 ± 0.17 | .50 |
Mitral valve inflow E wave to E′ ratio | 47 | 5.15 ± 1.24 | 209 | 5.54 ± 1.51 | .34 |
A-wave duration ratio (mitral inflow/pulmonary A reversal) | 20 | 1.22 ± 0.26 | 77 | 1.26 ± 0.40 | .35 |
The results for intraobserver and interobserver variability are summarized in Table 2 . There was good agreement between both serial measurements by the same investigator and blinded measurements by different, sequential observers.
Variable | Intraclass correlation coefficient |
---|---|
Interobserver | |
Torsion | 0.95 |
Net twist | 0.95 |
Torsion rate | 0.75 |
Untwist rate | 0.77 |
Intraobserver | |
Torsion | 0.98 |
Net twist | 0.98 |
Torsion rate | 0.81 |
Untwist rate | 0.90 |
Table 3 reports the rotational mechanics of patients with SCD compared with controls without SCD. Mean peak apical rotation was not different between the two groups, but peak basal rotation was significantly lower in the SCD group (−3.88 ± 2.3° vs −3.13 ± 2.5°). There was no significant difference in net twist, but controls had significantly higher peak differential twist (9.55 ± 3.93° vs 8.2 ± 3.25°) and peak normalized differential twist (1.42 ± 0.56°/cm vs 1.13 ± 0.46°/cm). There was no significant difference in mean torsion, but controls had a significantly higher torsion rate (12.06 ± 4.5°/sec/cm vs 9.53 ± 3.5°/sec/cm). Peak untwist rate was significantly slower in subjects with SCD compared with controls ( P = .006).
Variable | Controls ( n = 49) | Patients with SCD ( n = 218) | P | ||
---|---|---|---|---|---|
n | Mean ± SD | n | Mean ± SD | ||
Apical rotation (°) | 48 | 5.6 ± 2.7 | 200 | 5.0 ± 2.4 | .37 |
Basal rotation (°) | 47 | −3.9 ± 2.3 | 208 | −3.1 ± 2.5 | .012 |
Net twist (°) | 46 | 7.7 ± 4.1 | 198 | 6.9 ± 3.8 | .09 |
Peak differential twist (°) | 47 | 9.5 ± 3.9 | 198 | 8.2 ± 3.2 | .016 |
Peak normalized differential twist (°/cm) | 47 | 1.4 ± 0.6 | 198 | 1.1 ± 0.5 | .017 |
Torsion (°/cm) | 46 | 1.1 ± 0.6 | 198 | 0.9 ± 0.5 | .12 |
Torsion rate (°/sec/cm) | 46 | 12.1 ± 4.5 | 198 | 9.5 ± 3.5 | .0002 |
Untwist rate (°/sec) | 46 | −94.5 ± 37.1 | 198 | −83.0 ± 32.2 | .006 |
Normalized untwist rate (°/sec/cm) | 46 | −14.2 ± 5.8 | 198 | −11.4 ± 4.6 | .004 |
Time to peak apical rotation (% systole) | 48 | 86.4 ± 22.0 | 200 | 89.6 ± 25.9 | .23 |
Time to peak basal rotation (% systole) | 47 | 125.9 ± 26.2 | 208 | 116.6 ± 19.7 | .53 |
Time to peak torsion (% systole) | 46 | 102.6 ± 7.9 | 198 | 105.9 ± 9.9 | .16 |
Time to peak untwist rate (% systole) | 46 | 129.1 ± 10.2 | 198 | 131.3 ± 10.3 | .16 |
Apical rotation rate (°/sec) | 48 | 70.7 ± 28.3 | 200 | 57.7 ± 23.0 | .00127 |
Normalized apical rotation rate (°/sec/cm) | 48 | 10.6 ± 4.4 | 200 | 7.9 ± 3.3 | .00065 |
Basal rotation rate (°/sec) | 47 | −73.8 ± 24.4 | 208 | −64.90 ± 21.1 | .00048 |
Normalized basal rotation rate (°/sec/cm) | 47 | −11.0 ± 3.7 | 206 | −8.95 ± 3.0 | .00039 |
Net twist rate (°/sec) | 46 | 80.5 ± 30.7 | 199 | 69.12 ± 24.8 | .0004 |