Mechanical discoordination rather than dyssynchrony ... - CiteSeerX

8 downloads 0 Views 1MB Size Report
May 30, 2008 - Helm RH, Leclercq C, Faris OP, Ozturk C, McVeigh E, Lardo AC,. Kass DA. Cardiac dyssynchrony analysis using circumferential versus.
Am J Physiol Heart Circ Physiol 295: H640 –H646, 2008. First published May 30, 2008; doi:10.1152/ajpheart.00106.2008.

Mechanical discoordination rather than dyssynchrony predicts reverse remodeling upon cardiac resynchronization Borut Kirn,1,3 Annemieke Jansen,2 Frank Bracke,2 Berry van Gelder,2 Theo Arts,3 and Frits W. Prinzen4 1

Institute of Physiology, School of Medicine, University of Ljubljana, Slovenia; 2Department of Cardiology, Catharina Hospital, Eindhoven; and 3Department of Biomedical Engineering, and 4Department of Physiology, Cardiovascular Research Institute Maastricht, Maastricht University, Maastricht, The Netherlands

Submitted 31 January 2008; accepted in final form 27 May 2008

Kirn B, Jansen A, Bracke F, van Gelder B, Arts T, Prinzen FW. Mechanical discoordination rather than dyssynchrony predicts reverse remodeling upon cardiac resynchronization. Am J Physiol Heart Circ Physiol 295: H640 –H646, 2008. First published May 30, 2008; doi:10.1152/ajpheart.00106.2008.—By current guidelines a considerable part of the patients selected for cardiac resynchronization therapy (CRT) do not respond to the therapy. We hypothesized that mechanical discoordination [opposite strain within the left ventricular (LV) wall] predicts reversal of LV remodeling upon CRT better than mechanical dyssynchrony. MRI tagging images were acquired in CRT candidates (n ⫽ 19) and in healthy control subjects (n ⫽ 9). Circumferential strain (εcc) was determined in 160 regions. From εcc signals we derived 1) an index of mechanical discoordination [internal stretch fraction (ISF), defined as the ratio of stretch to shortening during ejection] and 2) indexes of mechanical dyssynchrony: the 10 –90% width of time to onset of shortening, time to peak shortening, and end-systolic strain. LV end-diastolic volume (LVEDV), end-systolic volume (LVESV), and ejection fraction (LVEF) were determined before and after 3 mo of CRT. Responders were defined as those patients in whom LVESV decreased by ⬎15%. In responders (n ⫽ 10), CRT increased LVEF and decreased LVEDV and LVESV (11 ⫾ 6%, 21 ⫾ 16%, and 30 ⫾ 16%, respectively) significantly more (P ⬍ 0.05) than in nonresponders (1 ⫾ 6%, 3 ⫾ 4%, and 5 ⫾ 10%, respectively). Among mechanical indexes, only ISF was different between responders and nonresponders (0.53 ⫾ 0.25 vs. 0.31 ⫾ 0.16; P ⬍ 0.05). In patients with ISF ⬎0.4 (n ⫽ 10), LVESV decreased by 31 ⫾ 18% vs. 5 ⫾ 11% in patients with ISF ⬍0.4 (P ⬍ 0.05). We conclude that mechanical discoordination, as estimated from ISF, is a better predictor of reverse remodeling after CRT than differences in time to onset and time to peak shortening. Therefore, discoordination rather than dyssynchrony appears to reflect the reserve contractile capacity that can be recruited by CRT. resynchronization therapy; cardiac mechanics

cardiac resynchronization therapy (CRT) has been shown to improve cardiac function, clinical status, and survival significantly compared with placebo control subjects. However, of all patients receiving CRT ⬃30% do not respond clinically to this therapy (15–17) and an even larger percentage of patients do not show reverse left ventricular (LV) remodeling (reversal of LV dilatation) (19, 27, 28). Some studies suggested that mechanical dyssynchrony can improve prediction of response to CRT on top of wide QRS complex (15). However, an increasing number of studies, including the recent PROSPECT multicenter trial, seem to disprove this idea (5). We hypothesized that the imperfect

IN LARGE CLINICAL TRIALS

Address for reprint requests and other correspondence: B. Kirn, Univ. of Maastricht, Dept. of Biomedical Engineering, Cardiovascular Research Inst. Maastricht, Maastricht, The Netherlands (e-mail: [email protected]). H640

prediction of CRT from mechanical dyssynchrony is due to the fact that pump function is dependent on coordination of contraction rather than on regional differences in timing of onset or peak of contraction. To investigate this hypothesis we developed a novel measure of discoordination, internal stretch fraction (ISF), defined as the amount of stretch relative to the amount of shortening during ventricular ejection. ISF and conventional dyssynchrony measures were determined with MRI-myocardial tagging in a group of CRT patients, and the baseline values of the discoordination and dyssynchrony indexes were related to the long-term outcome (reverse remodeling) of CRT in these patients. METHODS

Protocol: Patients and Control Subjects A group of 19 patients (age 67 ⫾ 8, male/female ⫽ 12/7), with [ischemic cardiomyopathy (ICM), n ⫽ 9] and without [dilated cardiomyopathy (DCM), n ⫽ 10] a history of ischemic disease and 9 healthy volunteers (control subjects) underwent measurements with tagged MRI. In patients, tagged MRI and echocardiography measurements (Philips Medical Systems) were performed 1 wk before the start of CRT and echocardiography was repeated 3 mo after the start to determine a change in LV end-diastolic volume (⌬LVEDV), endsystolic volume (⌬LVESV), and ejection fraction (⌬LVEF) in the apical four- and two-chamber views. The degree of mitral regurgitation (MR) was assessed by color jet area and as midsystolic percent jet area relative to left atrial size in the apical four-chamber view. The biventricular pacing leads (InSync 8042, Medtronic) at coronary sinus, right ventricle, and right atrium were positioned in a posterior or posterolateral branch, in the apical or midseptal region, and in the atrial appendage, respectively. The AV and VV delays were optimized as described previously (10, 11). Maximum rate of rise of LV pressure (dplv/dtmax) was measured with a micromanometer introduced into the LV (RADI Medical Systems) (11). The ethics committee of the Catharina Hospital approved the study protocol, and all subjects gave informed consent. MRI Scans MRI images were acquired (Philips Medical Systems, 1.5 T) with a 500-ms acquisition period started 20 ms after the ECG trigger. In five parallel short-axis cross sections MR tagging image series (⌬t ⫽ 20 ms) were acquired. The MR imaging parameters for each slice were as follows: echo time 10 ms, slice thickness 8 mm, field of view 213 mm, image size 256 ⫻ 256 pixels, and spatial modulation of the magnetization (SPAMM). Line tags were applied with intertag distance 5 mm, tag width 2.5 mm, and number of phases 16. The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

0363-6135/08 $8.00 Copyright © 2008 the American Physiological Society

http://www.ajpheart.org

H641

INTERNAL STRETCH FRACTION

Strain Measurements Tagged MR images were processed off-line to obtain midwall circumferential strain (εcc) in 32 sectors in each slice with custom software written in Matlab 7.0 (7, 23). In the LV midwall the fibers are directed predominantly in the circumferential direction. Each strain signal εcc(t) was filtered with a low-pass filter with a cutoff frequency of 12.5 s⫺1. Data Analysis Index of mechanical discoordination: internal stretch fraction. ISF was determined as the amount of stretch relative to the amount of shortening during the ejection phase. The “amount of stretch” was calculated as the product of the number of regions with stretch and their average amount of stretch, integrated over the ejection phase. A similar definition was used for the “amount of shortening.” The values of these variables are depicted as the size of the black and white areas under the curves in Figs. 1C, 2C, and 3C. To calculate ISF, first from each strain signal εcc(t) (Figs. 1A, 2A, 3A) a strain rate signal e(t) was derived by time differentiation (Figs. 1B, 2B, 3B). Then the strain rate signals were split into a lengthening signal ep(t) and a shortening signal en(t) with e共t兲 ⫽

dεcc共t兲 ep共t兲 ⫽ max共0,e共t兲兲 en共t兲 ⫽ min共0,e共t兲兲 dt

(1)

Fig. 2. Measured data and explanation of ISF concept for representative responder patient. For figure explanation see Fig. 1.

Averaging ep(t) and en(t) over all 160 signals resulted in average lengthening e៮ p(t) and average shortening e៮ n(t), respectively (Figs. 1B, 2B, 3B). Integration of these signals in the time interval [T1, T2] (Fig. 1C) results in a strain εp, representing the amount of stretch, and a strain εn, representing the amount of shortening: T2

T2



ε p ⫽ e៮ p共t兲dt



εn ⫽ e៮ n共t兲dt

T1

(2)

T1

We calculated ISF during the ejection phase; therefore, T1 was the start (begin ejection, BE) and T2 the end (end ejection, EE) of the ejection phase. BE and EE were determined from the pattern of the mean εcc(t) signal, as shown in Fig. 4. ISF was then calculated as ISF ⫽ ⫺

Fig. 1. Measured data and explanation of internal stretch fraction (ISF) concept for representative healthy control subject. A: tracings of circumferential strain (εcc) in all 160 regions. B: tracings of strain rate (thin gray lines) and the time course of the average lengthening [៮ep(t), solid black line] and average shortening [៮en(t), dotted black line]. C: the amount of stretch (εp) and shortening (εn) during ejection phase are represented with the area below the curve e៮ p(t) (black area) and above the curve e៮ n(t) (white area), respectively. BE, begin ejection; EE, end ejection. AJP-Heart Circ Physiol • VOL

εp εn

(3)

Equations 1–3 were implemented numerically by discretization to the frame rate. Indexes of mechanical dyssynchrony. From each strain signal three parameters were determined. 1) Onset of shortening time (Tonset) was determined by the first positive-to-negative zero crossing of strain rate e(t) as defined in Eq. 1. 2) Time to the first peak of shortening (Tpeak) was determined by the first negative-to-positive zero crossing of e(t) after Tonset. 3) End-systolic strain (εES) was determined as the value of εcc(t) at t ⫽ tES. To increase the precision of determination of Tonset and Tpeak, the values of εcc(t) time derivatives between two points were linearly interpolated. This approach yields Tonset variability within ⫾5 ms (26).

295 • AUGUST 2008 •

www.ajpheart.org

H642

INTERNAL STRETCH FRACTION

Figure 5 presents the percent decrease in LVEDV and LVESV and the change in LVEF after 3 mo of CRT and the percent increase of dplv/dtmax acutely upon CRT. Responders had significantly larger decrease in LVEDV and LVESV (%⌬LVEDV, %⌬LVESV, respectively) and a larger increase in LVEF, but the acute increase in dplv/dtmax (%⌬dplv/dtmax) was similar in both groups. Strain Analysis

Fig. 3. Measured data and explanation of ISF concept for representative nonresponder patient. For figure explanation see Fig. 1.

Subsequently, WTonset, WTpeak, and WεES were determined as the width of the interval between the 10th and 90th percentiles of Tonset, Tpeak, and εES, respectively. The width of these intervals reflects the degree of the mechanical dyssynchrony (25). WTonset and WTpeak are measures of intraventricular mechanical dyssynchrony, which are the equivalents of those determined with tissue Doppler velocities (4, 15, 28, 29). WεES was determined as well, because variation in εES is also increased within dyssynchronous ventricles (24). Statistics. Patients were divided post hoc into two groups, responders and nonresponders, defined as subjects with ⌬LVESV ⬎15% and ⱕ15%, respectively. Subsequently, all data were pooled within control, responder, and nonresponder groups. All data are presented as means ⫾ SD. Differences between groups were evaluated with oneway ANOVA; P ⬍ 0.05 was considered statistically significant. Changes over time within a group were analyzed by paired t-test.

Figures 1–3 show examples of strain analysis in three individuals. Within the LV wall of the healthy control subject εcc(t) signals were similar (Fig. 1A). In contrast, within the LV wall of the two patients (Figs. 2A and 3A), the strain signals were unalike. The amount of stretch [៮ep(t), solid line] and shortening [៮en(t)] is shown in Figs. 1B, 2B, and 3B. In the control subject, during ejection ៮ep(t) was negligible and ៮en(t) was large and the opposite was true during the filling phase. In contrast, in the patients ៮ep(t) was substantial, especially during early ejection in the responder patient. Figure 6 shows examples of the regional distribution of ISF, Tonset, and Tpeak. In the healthy control subject maps are almost uniform. In both the responder and nonresponder patients the yellow-red areas in the Tonset and Tpeak maps indicate regions with delayed onset and peak shortening, respectively. In both patients, maps of Tonset and Tpeak appear similar. In contrast, only in the responder patient is a region with high values of ISF observed, located predominantly in the anterior septum. Reverse Remodeling After 3 Months of CRT Is Associated with High ISF Both patient groups had significantly elevated values of WTonset, WTpeak, WεES, and ISF compared with the control group (Fig. 7). ISF values were significantly different between the responders and nonresponders. However, the two patient groups did not significantly differ with regard to any of the

RESULTS

Clinical Variables All patients survived during follow-up. By the criterion of a decrease in LVESV ⬎15% (see Statistics), 10 of 19 patients were responders and 9 were nonresponders. Of the responders 8 of 10 had DCM, and of the nonresponders 7 of 9 had ICM. Absolute values of clinical and echocardiographic variables are summarized in Table 1. Baseline characteristics were similar between the responder and nonresponder groups. Both groups improved clinically, as indicated by the decrease in New York Heart Association (NYHA) class. However, 3 mo of CRT only significantly increased 6-min walk test, cardiac index, and LVEF and significantly decreased LVESV in the responder group (Table 1). AJP-Heart Circ Physiol • VOL

Fig. 4. Definition of time points of BE and EE from mean εcc(t) signal. Tangent a meets the signal at maximum negative strain rate. BE and EE are defined in the intersection of tangent a and the tangents through the maximum (b) and the minimum (c), respectively.

295 • AUGUST 2008 •

www.ajpheart.org

H643

INTERNAL STRETCH FRACTION

Table 1. Clinical and echocardiographic variables before resynchronization (baseline) and after 3 mo of CRT Responders (n ⫽ 10)

New York Heart Association class 6-min walk test, ms LV ejection fraction, % LV end-diastolic volume, ml LV end-systolic volume, ml Mitral regurgitation/LA Cardiac index, l 䡠 min⫺1 䡠 m⫺2 dplv/dtmax, mmHg/s

Nonresponders (n ⫽ 9)

Baseline

3 mo

Baseline

3 mo

2.9⫾0.4 318⫾196 26⫾10 244⫾113 189⫾111 0.30⫾0.26 1.8⫾0.5 806⫾208

1.9⫾0.6* 389⫾119* 37⫾14* 192⫾96* 132⫾90* 0.18⫾0.20* 2.1⫾0.3* 1,021⫾244*

3.1⫾0.3 426⫾62 24⫾7 242⫾103 187⫾92 0.34⫾0.27 1.7⫾0.5 725⫾143

2.5⫾0.7* 458⫾62 25⫾8 234⫾96* 178⫾81 0.26⫾0.23* 1.7⫾0.3 916⫾121*

Values are means ⫾ SD for n patients. At baseline there was no difference between responder and nonresponder groups. CRT, cardiac resynchronization therapy; LV, left ventricle; LA, left atrium size (mm); dplv/dtmax, maximum rate of rise of LV pressure. *P ⬍ 0.05 vs. baseline.

other mechanical indexes or QRS duration or LVESV at baseline. Figure 8 presents the individual values of the indexes of reverse remodeling (%⌬LVEDV, %⌬LVESV, and ⌬LVEF) as a function of ISF and WTonset. In Fig. 8, top left, it is shown that all but two responders had an ISF ⬎0.4 and all but two nonresponders had an ISF ⬍0.4. When the patients were divided into groups with ISF ⬎0.4 and ⱕ0.4, the resulting groups with a high ISF (0.62 ⫾ 0.14, n ⫽ 10) had significantly more ventricular reverse remodeling than those with a low ISF (0.22 ⫾ 0.10; P ⬍ 0.05): %⌬LVESV ⫽ 31 ⫾ 18% vs. 5 ⫾ 11%, %⌬LVEDV ⫽ 20 ⫾ 16% vs. 3 ⫾ 5%, and ⌬LVEF ⫽ 10 ⫾ 6% vs. 1 ⫾ 6%, respectively (all P ⬍ 0.05). Figure 8 also shows that most ICM patients have lower ISF than most DCM patients, but there is some overlap between the data points of these two disease categories. In contrast, the value of WTonset was not associated with any of the reverse remodeling parameters, as shown in Fig. 8, right. When the patients were divided in groups with WTonset greater than and less than or equal to the median of WTonset (0.08 s) (WTonset: 0.09 ⫾ 0.01 s vs. 0.06 ⫾ 0.01 s, respectively; P ⬍ 0.05), no difference between groups with high (n ⫽ 10) and low WTonset was found: %⌬LVESV ⫽ 20 ⫾ 21% vs. 17 ⫾

Fig. 5. Left: % decrease in left ventricular (LV) end-diastolic volume (%⌬LVEDV) and end-systolic volume (%⌬LVESV) and increase of LV ejection fraction (⌬LVEF) after 3 mo of cardiac resynchronization therapy (CRT). Right: % increase of maximum rate of rise of LV pressure (%⌬dplv/ dtmax) acutely after pacemaker implant. For the sake of presentation a decrease in LVESV and LVEDV is depicted as positive (#P ⬍ 0.05 between responder and nonresponder groups). AJP-Heart Circ Physiol • VOL

19%, %⌬LVEDV ⫽ 12 ⫾ 16% vs. 12 ⫾ 14%, and ⌬LVEF ⫽ 8 ⫾ 8% vs. 4 ⫾ 8%, respectively. DISCUSSION

The present study demonstrates that the baseline value of the discoordination index ISF provides a prediction of reversal of LV remodeling that is considerably better than that from conventional indexes of mechanical dyssynchrony or from the acute increase in contractility (dplv/dtmax). Therefore, these data support the idea that an index of discoordination can contribute importantly to a good prediction of the volume response to CRT. Moreover, the lack of such prediction by the indexes of mechanical dyssynchrony may at least partly explain the limited added value of indexes of echocardiographically determined mechanical dyssynchrony to predict response to CRT (5). Reduction of dilatation indicates structural recovery at the organ level. Absence of an increase in LVEF or a decrease in LVESV does not exclude a beneficial effect of CRT on the clinical status of the patient. This can be derived from the significant decrease in NYHA class in the responders and nonresponders in our study. However, the NYHA class is

Fig. 6. Maps of ISF, time to onset of shortening (Tonset), and time to first peak shortening (Tpeak) for a healthy subject (control) and 2 patients (nonresponder and responder). The maps are shown as bull’s-eye plots, where the LV wall is represented as if it were looked at from the apex (Ap), which is represented with the inner circle, with the base (B) being the outer contour. The septum (S) and posterior (P), lateral (L) and anterior (A) walls are indicated. For the sake of presentation, smoothing was applied. Time scales in maps of Tonset and Tpeak are given as a difference from the average.

295 • AUGUST 2008 •

www.ajpheart.org

H644

INTERNAL STRETCH FRACTION

Fig. 7. Baseline values of various indexes for the control, nonresponder, and responder groups. WTonset, WTpeak, WεES, width of interval between 10th and 90th percentiles for Tonset, Tpeak, and strain at end systole (εES); QRS, QRS duration; ESV, end-systolic volume. *P ⬍ 0.05 between control subjects and patients; #P ⬍ 0.05 between responder and nonresponder patients.

subject to a placebo effect (1), and the most pronounced improvements, including improved survival, prevail in patients with reduction of dilatation (27). The mechanical discoordination may have two origins: asynchrony of activation as a result of conduction disorder and/or regional deterioration of contractility as a result of ischemia or infarction (12). Of the two only the former can be compensated by CRT. The infarcted myocardium forms scar, which is less compliant than viable nonactivated tissue and hence at raised pressure stretches less. The scar tissue cannot be recruited by CRT and does not contribute considerably to ISF value.

disease, heart failure) and may therefore explain the low ISF in some of the DCM patients. Clearly, future studies are required to investigate the cause of low ISF values in asynchronous hearts in more detail. A theoretical advantage of ISF over conventional measures of mechanical dyssynchrony is that it takes into account regional dispersion of strain throughout the entire ejection phase, rather than at a single time point in a strain signal. A practical advantage of this difference is also that when using ISF one does not rely on proper detection of the onset or peak of shortening, which can be problematic in cases of gradual changes, noisy signals, and low sampling rates. In our study strain was determined in 160 regions within 5 short-axis slices of the LV wall, covering at least 70% of that wall. Nowadays, myocardial strains can also be determined with speckle tracking, and such strains can, theoretically, be used equally well for ISF analysis. The ISF maps in Fig. 6 show that the major dispersion of ISF occurs in the septumto-lateral wall direction, so ISF could be reliably determined with a few well-chosen echocardiographic sections. Besides discoordination in ventricular mechanics during ejection phase there is also discoordination during filling phase. In Figs. 2A and 3A, we can see both as opposite strain slopes during each phase. We determined the amount of discoordination during filling phase with ISFd, which is a conceptual equivalent to ISF. ISFd was calculated as a ratio between the amount of shortening and stretch (inversely to ISF). ISFd is not as predictive for the CRT outcome as ISF. Although there are differences between responders and nonre-

The ISF Concept The observation that ISF is related to the increase in LVEF supports the idea that reduction of opposing systolic strain patterns leads to more effective ejection of blood into the arterial circulation. It is plausible that such a better ejection helps to let the LV operate at a smaller volume, thereby allowing a decrease of LV dilatation in the long run. Figure 8 shows that CRT response can be predicted from the ISF value but, however, also fairly well from the disease categories DCM and ICM alone. Preliminary results from animal experiments in our laboratory show that ISF increases because of creation of left bundle branch block (LBBB) but not because of chronic infarction. Therefore, a low ISF may indicate the absence of LBBB (even despite a wide QRS complex). An alternative explanation would be that in asynchronous ventricles large opposing strains are reduced because of increased myocardial stiffness. The source of this increased stiffness may be scar tissue resulting from infarctions. The amount of scar varies largely among ICM patients, some having no scar at all. Other sources of decreased myocardial compliance may be generalized fibrosis, as a result of myocardial remodeling, and altered properties of titin (14). Such generalized changes have been reported for various abnormal cardiac overload conditions (hypertension, valvular AJP-Heart Circ Physiol • VOL

Fig. 8. Measured % change of LVESV (%⌬LVESV) and LVEDV (%⌬LVEDV) and change of LVEF (⌬LVEF) after 3 mo of cardiac resynchronization therapy shown as a function of ISF and WTonset. If %⌬LVESV ⬎15% the patients were defined as responders (F); otherwise they were nonresponders (cross). Ischemic cardiomyopathy (ICM) patients are indicated by an open circle in the center. For the sake of presentation a decrease in LVESV and LVEDV is depicted as positive.

295 • AUGUST 2008 •

www.ajpheart.org

INTERNAL STRETCH FRACTION

sponders in ISFd (0.54 ⫾ 0.22 vs. 0.40 ⫾ 0.16, respectively) they are not significant. ISF is conceptually related to internal flow fraction (IFF), estimated from measurements of longitudinal segmental volume change by the conductance catheter technique. IFF was shown to be increased in heart failure patients with conduction disorders and reduced during CRT (20, 21), but a relation with long-term response to CRT has not been published until now. It is the advantage of ISF that it can be measured noninvasively and that septal-to-lateral wall differences are taken into account. The TUS (or CURE) index also includes information about coordination (18). This index is calculated with Fourier analysis of strain in LV regions in space. The TUS index reflects the relative first-order Fourier power within the LV wall. In doing so, TUS assumes a sine wave spatial variation in strains and it does not express the amount of stretch. ISF is calculated without any assumption on distribution of strain in space and time. Other studies have employed a vector analysis of mechanical dyssynchrony (9, 13, 25). This reflects the main orientation of the major component of dyssynchrony, but again does not reflect the amount of stretch. Mechanical Dyssynchrony A large number of single-center studies have reported that septal-to-lateral delay and/or the standard deviation of the time to peak velocity values of 12 segments have a high sensitivity/ specificity to predict CRT responders (3). However, the recent PROSPECT multicenter study (5) and several smaller singlecenter studies show that the prediction of responders is moderate and hardly better than that by QRS width alone (2, 3, 6, 18). On the basis of these findings, a recent group of experts declared in a consensus statement on behalf of the American Society of Echocardiography that the decision for treating a patient with CRT should not be based on or withheld because of results of an echocardiographic Doppler dyssynchrony study (8). An extensive discussion about the possible explanations of the current confusion in this area is beyond the scope of this article. Relevant to the results from the present study is that it has been argued that the poor performance of the indexes would, in part, be related to limitations of tissue Doppler imaging (TDI) and other echocardiographic techniques and potential subjective aspects in the analysis (6). In the present study we determined mechanical dyssynchrony with MRI tagging, the gold standard of myocardial deformation, and using completely objective analysis. The only manual, subjective, part in the entire analysis was the manual contouring of endocardial and epicardial borders of the LV wall in order to create the region of interest. Moreover, we used MRI tagging measurement covering a large part of the LV wall and with a high spatial resolution, making it unlikely that important information was missed. Despite the use of this state-of-the-art technique for measuring myocardial strains, the indexes related to a single point in time (time of onset or peak shortening, value of strain at end systole) did not differentiate between responders and nonresponders. This observation suggests that the disappointing results of various recent echocardiographic studies on prediction of CRT response may, in part, be exAJP-Heart Circ Physiol • VOL

H645

plained by the use of timing parameters and that potentially the use of indexes of dyscoordination may improve prediction of CRT response. Indexes of Mechanical Dyssynchrony We used WT measures based on 10th–90th percentile intervals, which is a rather sophisticated approach compared with more traditional parameters like septal-to-lateral delay. It can be seen in Fig. 6 that there is quite some variation in Tpeak and Tonset along a hypothetical ultrasound beam in the basal septal and lateral wall. Accordingly, the value of the septum-to-lateral wall delay is heavily dependent on the regions chosen and is therefore less likely to predict CRT response. We considered it fair to the dyssynchrony measures to use a measure that takes into account the variation in all regions, similar to the case with ISF. Mitral Regurgitation and ISF In this study we determined ISF only during the ejection phase, hence avoiding isovolumic contraction and relaxation phases where most of the effect of MR happens. In case of MR during the ejection phase, this would reduce the value of ISF, since shortening does not discriminate between ejection into the aorta and atrium. In exploring the best possible criteria, we also investigated whether calculation of ISF during the entire systole, including the isovolumic phases, would improve prediction of CRT response, and this was not the case. Predictive Value of dplv/dtmax In various scientific studies as well as in some clinical practice dplv/dtmax is being used to test the acute hemodynamic effect of CRT (11). While this information is most likely relevant, there is, to our knowledge, no publication showing the correlation between the acute hemodynamic effect and long-term reverse remodeling. In contrast, already in 2001 Stellbrink et al. (22) reported that acute hemodynamic improvement did not predict long-term reverse remodeling. Therefore, the finding in our small group of patients that the acute hemodynamic response does not predict long-term reverse remodeling is not necessarily in conflict with other data in the literature. It should, however, be kept in mind that this study was by no means designed to investigate this question. Obviously, larger studies are required to investigate this question in more detail. Critiques on the Study Design We chose to study the feasibility of the ISF analysis by using the best technique to measure myocardial strains, the relatively expensive and complex MRI tagging technique. As a consequence, the number of patients included in the study was relatively small. Therefore, only differences between the two patient groups can be compared, as opposed to determining specificity and sensitivity. Despite the small groups, ISF proved to be consistently different between responders and nonresponders, whereas conventional indexes of mechanical dyssynchrony did not differentiate at all between responders and nonresponders. The data in Fig. 8 show that the presence of ICM was also predictive of nonresponse to CRT. To what extent ISF is

295 • AUGUST 2008 •

www.ajpheart.org

H646

INTERNAL STRETCH FRACTION

coupled to ICM and DCM and whether ISF is a better predictor than ICM/DCM should be proven in a larger study, which is currently being performed. The patients included in this study can be regarded as typical CRT patients because of their degree of heart failure (NYHA class III) and QRS width (⬃160 ms), and all showed a clear increase in dplv/dtmax immediately on the start of CRT. It remains to be explored whether ISF can also be useful in the selection of patients with more moderate degrees of heart failure.

12. 13. 14.

15.

Conclusions All indexes of dyssynchrony and discoordination (WTonset, WTpeak, WεES, and ISF) discriminate normal from asynchronous hearts. However, the discoordination index ISF is the only one of these parameters that discriminates between patients in whom CRT does or does not reverses LV dilatation. Discoordination rather than dyssynchrony appears to reflect the reserve contractile capacity that can be recruited by CRT.

16.

17.

18.

REFERENCES 1. Abraham WT. Cardiac resynchronization therapy. Prog Cardiovasc Dis 48: 232–238, 2006. 2. Achilli A, Peraldo C, Sassara M, Orazi S, Bianchi S, Laurenzi F, Donati R, Perego GB, Spampinato A, Valsecchi S, Denaro A, Puglisi A. Prediction of response to cardiac resynchronization therapy: the selection of candidates for CRT (SCART) study. Pacing Clin Electrophysiol 29, Suppl 2: S11–S19, 2006. 3. Bax JJ, Abraham T, Barold SS, Breithardt OA, Fung JW, Garrigue S, Gorcsan J 3rd, Hayes DL, Kass DA, Knuuti J, Leclercq C, Linde C, Mark DB, Monaghan MJ, Nihoyannopoulos P, Schalij MJ, Stellbrink C, Yu CM. Cardiac resynchronization therapy. 1. Issues before device implantation. J Am Coll Cardiol 46: 2153–2167, 2005. 4. Breithardt OA, Stellbrink C, Herbots L, Claus P, Sinha AM, Bijnens B, Hanrath P, Sutherland GR. Cardiac resynchronization therapy can reverse abnormal myocardial strain distribution in patients with heart failure and left bundle branch block. J Am Coll Cardiol 42: 486 – 494, 2003. 5. Cleland JG, Abdellah AT, Khaleva O, Coletta AP, Clark AL. Clinical trials update from the European Society of Cardiology Congress 2007: 3CPO, ALOFT, PROSPECT and statins for heart failure. Eur J Heart Fail 9: 1070 –1073, 2007. 6. De Boeck BWL, Meine M, Leenders GE, Teske AJ, van Wessel H, Kirkels JH, Prinzen FW, Doevendans PA, Cramer MJ. Practical and conceptual limitations of tissue Doppler imaging to predict reverse remodelling in cardiac resynchronisation therapy. Eur J Heart Fail 10: 281–290, 2008. 7. Delhaas T, Kotte J, van der Toorn A, Snoep G, Prinzen FW, Arts T. Increase in left ventricular torsion-to-shortening ratio in children with valvular aortic stenosis. Magn Reson Med 51: 135–139, 2004. 8. Gorcsan J 3rd, Abraham T, Agler DA, Bax JJ, Derumeaux G, Grimm RA, Martin R, Steinberg JS, Sutton MS, Yu CM. Echocardiography for cardiac resynchronization therapy: recommendations for performance and reporting—a report from the American Society of Echocardiography Dyssynchrony Writing Group endorsed by the Heart Rhythm Society. J Am Soc Echocardiogr 21: 191–213, 2008. 9. Helm RH, Leclercq C, Faris OP, Ozturk C, McVeigh E, Lardo AC, Kass DA. Cardiac dyssynchrony analysis using circumferential versus longitudinal strain: implications for assessing cardiac resynchronization. Circulation 111: 2760 –2767, 2005. 10. Jansen AH, Bracke F, van Dantzig JM, Meijer A, Korsten EH, Peels KH, van Hemel NM. Optimization of pulsed wave tissue Doppler to predict left ventricular reverse remodeling after cardiac resynchronization therapy. J Am Soc Echocardiogr 19: 185–191, 2006. 11. Jansen AH, Bracke FA, van Dantzig JM, Meijer A, van der Voort PH, Aarnoudse W, van Gelder BM, Peels KH. Correlation of echo-Doppler optimization of atrioventricular delay in cardiac resynchronization therapy with invasive hemodynamics in patients with heart failure secondary to AJP-Heart Circ Physiol • VOL

19.

20.

21.

22.

23.

24. 25.

26. 27.

28.

29.

ischemic or idiopathic dilated cardiomyopathy. Am J Cardiol 97: 552–557, 2006. Kass DA. An epidemic of dyssynchrony: but what does it mean? J Am Coll Cardiol 51: 12–17, 2008. Lardo AC, Abraham TP, Kass DA. Magnetic resonance imaging assessment of ventricular dyssynchrony: current and emerging concepts. J Am Coll Cardiol 46: 2223–2228, 2005. Nagueh SF, Shah G, Wu Y, Torre-Amione G, King NM, Lahmers S, Witt CC, Becker K, Labeit S, Granzier HL. Altered titin expression, myocardial stiffness, and left ventricular function in patients with dilated cardiomyopathy. Circulation 110: 155–162, 2004. Penicka M, Bartunek J, De Bruyne B, Vanderheyden M, Goethals M, De Zutter M, Brugada P, Geelen P. Improvement of left ventricular function after cardiac resynchronization therapy is predicted by tissue Doppler imaging echocardiography. Circulation 109: 978 –983, 2004. Pitzalis MV, Iacoviello M, Romito R, Guida P, De Tommasi E, Luzzi G, Anaclerio M, Forleo C, Rizzon P. Ventricular asynchrony predicts a better outcome in patients with chronic heart failure receiving cardiac resynchronization therapy. J Am Coll Cardiol 45: 65– 69, 2005. Pitzalis MV, Iacoviello M, Romito R, Massari F, Rizzon B, Luzzi G, Guida P, Andriani A, Mastropasqua F, Rizzon P. Cardiac resynchronization therapy tailored by echocardiographic evaluation of ventricular asynchrony. J Am Coll Cardiol 40: 1615–1622, 2002. Soliman OI, Theuns DA, Geleijnse ML, Anwar AM, Nemes A, Caliskan K, Vletter WB, Jordaens LJ, Cate FJ. Spectral pulsed-wave tissue Doppler imaging lateral-to-septal delay fails to predict clinical or echocardiographic outcome after cardiac resynchronization therapy. Europace 9: 113–118, 2007. St John Sutton MG, Plappert T, Abraham WT, Smith AL, DeLurgio DB, Leon AR, Loh E, Kocovic DZ, Fisher WG, Ellestad M, Messenger J, Kruger K, Hilpisch KE, Hill MR. Effect of cardiac resynchronization therapy on left ventricular size and function in chronic heart failure. Circulation 107: 1985–1990, 2003. Steendijk P, Tulner SA, Bax JJ, Oemrawsingh PV, Bleeker GB, van Erven L, Putter H, Verwey HF, van der Wall EE, Schalij MJ. Hemodynamic effects of long-term cardiac resynchronization therapy: analysis by pressure-volume loops. Circulation 113: 1295–1304, 2006. Steendijk P, Tulner SA, Schreuder JJ, Bax JJ, van Erven L, van der Wall EE, Dion RA, Schalij MJ, Baan J. Quantification of left ventricular mechanical dyssynchrony by conductance catheter in heart failure patients. Am J Physiol Heart Circ Physiol 286: H723–H730, 2004. Stellbrink C, Breithardt OA, Franke A, Sack S, Bakker P, Auricchio A, Pochet T, Salo R, Kramer A, Spinelli J. Impact of cardiac resynchronization therapy using hemodynamically optimized pacing on left ventricular remodeling in patients with congestive heart failure and ventricular conduction disturbances. J Am Coll Cardiol 38: 1957–1965, 2001. Van Der Toorn A, Barenbrug P, Snoep G, Van Der Veen FH, Delhaas T, Prinzen FW, Maessen J, Arts T. Transmural gradients of cardiac myofiber shortening in aortic valve stenosis patients using MRI tagging. Am J Physiol Heart Circ Physiol 283: H1609 –H1615, 2002. Vernooy K, Verbeek XA, Peschar M, Crijns HJ, Arts T, Cornelussen RN, Prinzen FW. Left bundle branch block induces ventricular remodelling and functional septal hypoperfusion. Eur Heart J 26: 91–98, 2005. Wyman BT, Hunter WC, Prinzen FW, Faris OP, McVeigh ER. Effects of single- and biventricular pacing on temporal and spatial dynamics of ventricular contraction. Am J Physiol Heart Circ Physiol 282: H372– H379, 2002. Wyman BT, Hunter WC, Prinzen FW, McVeigh ER. Mapping propagation of mechanical activation in the paced heart with MRI tagging. Am J Physiol Heart Circ Physiol 276: H881–H891, 1999. Yu CM, Bleeker GB, Fung JW, Schalij MJ, Zhang Q, van der Wall EE, Chan YS, Kong SL, Bax JJ. Left ventricular reverse remodeling but not clinical improvement predicts long-term survival after cardiac resynchronization therapy. Circulation 112: 1580 –1586, 2005. Yu CM, Chau E, Sanderson JE, Fan K, Tang MO, Fung WH, Lin H, Kong SL, Lam YM, Hill MR, Lau CP. Tissue Doppler echocardiographic evidence of reverse remodeling and improved synchronicity by simultaneously delaying regional contraction after biventricular pacing therapy in heart failure. Circulation 105: 438 – 445, 2002. Zwanenburg JJ, Gotte MJ, Marcus JT, Kuijer JP, Knaapen P, Heethaar RM, van Rossum AC. Propagation of onset and peak time of myocardial shortening in time of myocardial shortening in ischemic versus nonischemic cardiomyopathy: assessment by magnetic resonance imaging myocardial tagging. J Am Coll Cardiol 46: 2215–2222, 2005.

295 • AUGUST 2008 •

www.ajpheart.org