European Heart Journal Advance Access originally published online on October 24, 2007
European Heart Journal 2007 28(22):2738-2748; doi:10.1093/eurheartj/ehm443
Different effects of cardiac resynchronization therapy on left atrial function in patients with either idiopathic or ischaemic dilated cardiomyopathy: a two-dimensional speckle strain study
1 Monaldi Hospital, Second University of Naples, Via G. Martucci 35, 80121 Naples, Italy
2 Department of Internal Medicine and Public Health, University of L'Aquila, L'Aquila, Italy
3 Dipartimento di Internistica Clinica e Sperimentale, Sezione Scientifica di Diagnostica per Immagini, Second University of Naples, Naples, Italy
4 Department of Cardiothoracic and Respiratory Sciences, Second University of Naples, Naples, Italy
Received 28 February 2007; revised 23 July 2007; accepted 11 September 2007; online publish-ahead-of-print 24 October 2007.
* Corresponding author. Tel: +39 081 7618525; fax: +39 081 7145205. E-mail address: antonellodandrea{at}libero.it
| Abstract |
|---|
|
|
|---|
Aims: In dilated cardiomyopathy (DCM), attenuation of left atrial (LA) booster pump function has been observed, and attributed both to altered LA loading conditions owing to left ventricular (LV) diastolic dysfunction and to LA involvement in the myopathic process. The aim of the present study was to detect LA systolic dysfunction in DCM using speckle-tracking two-dimensional strain echocardiography (2DSE), and to assess the effects of cardiac resynchronization therapy (CRT) on LA myocardial strain during 6 month follow-up.
Methods and results: A total of 90 patients (aged, 52.4 ± 10.2 years) with either idiopathic (n = 47) or ischaemic (n = 43) DCM underwent standard Doppler echo and 2DSE analysis of atrial longitudinal strain in the basal segments of LA septum and LA lateral wall, and in LA roof. The two groups were comparable for clinical variables (NYHA class: III in 72.2%; IV in 27.8%). LV volumes, ejection fraction, stroke volume, and mitral valve effective regurgitant orifice were similar between the two groups. No significant differences were evidenced in Doppler transmitral inflow measurements. LA diameter and maximal volume were also similar between the two groups. Conversely, LA active emptying volume and fraction were both lower in patients with idiopathic DCM. Peak systolic myocardial atrial strain was significantly compromised in patients with idiopathic DCM compared with ischaemic DCM in all the analysed atrial segments (P < 0.001). At follow-up, 64 patients (71.1%) (37 idiopathic and 27 ischaemic) were responders, and 26 (28.9%) (10 idiopathic; 16 ischaemic) were non-responders to CRT (responder: decrease of LV end-systolic volume >15%). A significant improvement in LA systolic function was obtained only in patients with ischaemic DCM responders to CRT (P < 0.001). By multivariable analysis, in the overall population, it was found that ischaemic aetiology of DCM (ß-cofficient = 0.62; P < 0.0001) and positive response to CRT (ß-cofficient = 0.42; P < 0.01) were the only independent determinants of LA lateral wall systolic strain.
Conclusions: Two-dimensional strain represents a promising non-invasive technique to assess LA atrial myocardial function in patients with DCM. LA pump and reservoir function at baseline and after CRT are more depressed in idiopathic compared with ischaemic DCM patients. Future longitudinal studies are warranted to understand further the natural history of LA myocardial function, the extent of reversibility of LA dysfunction with CRT, and the possible prognostic impact of such indexes in patients with congestive heart failure.
Key Words: Heart failure Resynchronization therapy Idiopathic dilated cardiomyopathy Two-dimensionalstrain imaging Left atrium Left atrial function
Cardiac resynchronization therapy (CRT) has become an attractive therapeutic option for patients with end-stage chronic heart failure (HF), left ventricular (LV) dilatation, and left bundle branch block (LBBB) with wide QRS duration.1–2 Currently, patients are selected for CRT on the basis of clinical and ECG criteria and of standard LV echocardiographic indexes.3–4
Among various echocardiographic techniques, Doppler myocardial imaging (DMI) and strain rate imaging have gained their acceptance by virtue of their ability to define myocardial timing, contractility, and deformation in patients with LBBB and HF, and are highly feasible and easy-repeatable.5–12 Although DMI and Doppler strain measures have been used most often in clinical practice, they are limited by Doppler angle of incidence. A novel approach to quantify regional LV function from routine gray-scale 2D echocardiographic images, known as speckle-tracking two-dimensional strain echocardiography (2DSE), calculates myocardial strain independent of angle of incidence, and has been recently validated against sonomicrometry and tagged magnetic resonance imaging.13–16
Augmentation of left atrial (LA) systolic function and active LA contribution to LV filling (booster pump function) is a common finding in conditions associated with LV diastolic dysfunction.17 Conversely, in dilated cardiomyopathy (DCM), attenuation of LA booster pump function has been observed, and attributed both to altered LA loading conditions owing to LV diastolic dysfunction and to LA involvement in the myopathic process.18–19 By the use of standard echocardiography, previous authors have showed a LA systolic function more depressed in idiopathic compared with ischaemic DCM.20–22 To date, no previous report analysed the potential effects of CRT on LA myocardial function.
On these grounds, aims of the present study were: (i) to detect by speckle-tracking 2D strain underlying LA myocardial dysfunction in patients with either idiopathic or ischaemic DCM with comparable atrial mechanical overload; (ii) to assess possible different effects of CRT on LA myocardial strain during the 6 month follow-up.
| Methods |
|---|
|
|
|---|
Study population
From January 2005 to April 2006, 455 ambulatory HF patients with known LV systolic dysfunction were referred to our echocardiographic laboratory. Patients with QRS value of <120 ms (n = 92), atrial fibrillation (n = 69), inducible myocardial ischaemia (n = 25), poor echocardiographic window (n = 11) and mildly impaired LV ejection fraction (EF) (>40%) (n = 36) were ineligible for our study. Among the remaining 222 patients (120 with idiopathic DCM and 102 with ischaemic DCM), 90 patients were prospectively studied and selected for CRT for severe HF despite optimal pharmacological therapy, and were followed up for at least 6 months after the implantation.
The final study population, therefore, included 90 patients (mean age, 52.4 ± 10.2 years) with either idiopathic (47 patients) or ischaemic (43 patients) DCM. Inclusion criteria for CRT were: NYHA class III–IV refractory HF, LV end-diastolic diameter value of >55 mm, LV EF of <35%, QRS interval of >120 ms, and sinus rhythm.1–2 Exclusion criteria were: acute HF, coronary artery by-pass graft surgery, or myocardial infarction within the previous 3 months, valvular stenosis, and previous valve replacement or reconstruction. The study was approved by the local Ethics Committee. All gave their informed consent.
Since 2DSE represents a novel technique with few clinical applications, we also studied LA myocardial function in 25 age- and sex-matched subjects without detectable cardiovascular risk factor. Volunteer controls were all recruited in Naples (Italy), were selected from our departments of Cardiology amongst subjects investigated for work eligibility, and were examined in a single centre (Monaldi Hospital, Naples, Italy). None of the control subjects had cardiovascular structural or functional abnormalities and received any medication.
Implantation technique
Cardiac resynchronization therapy was initiated with implantation of a biventricular pacing system (CONTAK CD H115, CD II H119, RENEWAL H135, or III H170, H175, H177, Guidant Corp; InSync ICD 7272, Marquis 7277, or InSync ICD II Marquis 7289, Medtronic). Three transvenous pacing leads were inserted: one in right atrium and another in the high interventricular septum or in the right ventricular apex. A coronary sinus lead was positioned on the LV free wall through a coronary sinus tributary veins. LV leads were positioned as follows: posterior–lateral or lateral in 72 patients, and anterior or anterior–lateral in 18 patients.
To optimize the benefit from CRT, in the overall population of DCM patients, we performed an echo-guided positioning of the LV pacing lead, determining by use of DMI the site of latest mechanical activation.3 In addition, among the 43 patients with ischaemic cardiomyopathy, pre-implantation evaluation of myocardial thickness of the region targeted for LV pacing was assessed. As reported by others, we considered the presence of myocardial scar tissue in case of measurement of end-diastolic wall thickness (
0.6 cm).23 In particular, end-diastolic wall thickness was assessed at the centre of each myocardial segment from the leading endocardial edge to leading epicardial edge using the echocardiographic window that best identified the endocardial and epicardial borders. In our population, myocardial scar was observed in the anterior/antero-lateral position in 29 ischaemic patients, in the posterior/postero-lateral position in eight patients. In addition, in six ischaemic DCM patients, a normal myocardial thickness (>0.6 cm) was observed, and contrast-enhanced cardiac magnetic resonance was used to localize and measure the extent of the myocardial scar. As a consequence, in 27 of 29 patients with anterior/antero-lateral scar, LV lead was implanted in the posterio-lateral position, whereas in six of eight patients with posterior/postero-lateral scar LV lead was implanted in anterio-lateral position. Only in four patients, patient-tailored LV lead positioning was not possible as a result of the limited number of suitable branches of the coronary sinus.
The pacing leads were connected to a dual-chamber biventricular device with programmable inter-ventricular delay. The device was programmed in DDD mode and the atrio-ventricular (AV) delay was adjusted during simultaneous pacing. Optimal AV delay was determined first using Ritter's method24 followed by optimization of interventricular (VV) delay. According to Ritter et al.,24 the longest filling time without truncation of the A wave was chosen by means of pulsed Doppler analysis of transmitral flow. Readjustment of AV delay was done again after VV optimization. LV pre-excitation of 4–50 ms and right ventricular pre-excitation of 4–40 ms were tested in each patient. Post-optimization mean atrioventricular delay was 120 ± 55 ms, with a wide variability in optimized AV delay intervals (60–180 ms) among the patient cohort, whereas post-optimization mean VV delay was 35 ± 15 ms.
Study protocol
One week before and 6 months after implantation, all the patients underwent a clinical examination, 12-lead ECG, standard Doppler echo, and two-dimensional strain by Vivid 7 ultrasound system (GE Vingmed Ultrasound). The patients were considered as responders to CRT if LV end-systolic volume decreased by 15% and as non-responders in all other cases.2
Standard echocardiography
Standard Doppler echocardiography and 2DSE were performed with the subjects in partial left decubitus. A variable frequency phased-array transducer (2.5–3.5–4.0 MHz) was used for 2DSE and Doppler imaging. Doppler echocardiographic tracings were recorded on magneto-optical disk. All the measurements were analysed by two experienced readers blinded to the aetiology and responder status of DCM patients, on the average of
3 cardiac cycles. Stroke volume was obtained by LV outflow Doppler method as the product between outflow tract area and LV output time–velocity integral.25 LV EF was measured using a commercially available software program that applied modified Simpson's rule on the two- and four-chamber views. Tricuspid annular plane systolic excursion (TAPSE) was calculated as index of right ventricular global systolic function by the difference between end-diastolic and end-systolic measurement (in millimetres). The early (Ea) and late (Aa) diastolic annular velocities were measured at the lateral corner of the mitral annulus by pulsed DMI, in accordance with the method proposed by Nagueh et al.26 Mitral Ea, corrected for the influence of relaxation (i.e. the E/Ea ratio), was calculated to estimate LV filling pressures.
The proximal flow convergence (PFC) technique has been validated as a quantitative Doppler method to calculate regurgitant volume (RV) of flow and orifice area [effective regurgitant orifice (ERO)]. The regurgitant flow is measured as 2
r2 x Vr, where r is the radius of the hemispheric PFC region and Vr is the aliasing velocity. The following parameters are calculated: ERO = regurgitant flow/maximal regurgitant velocity, and RV = ERO x RTVI, where RTVI being the regurgitant time–velocity integral. At least three consecutive beats of sinus rhythm were measured and the average value was taken.27
According to American Society of Echocardiography guidelines, M-mode measurement of LA diameter was performed in parasternal long-axis view at end-systole, measured at the maximum dimension from the leading edge of the posterior wall of the aorta to the dominant line of posterior wall of LA.17 LA volumes were obtained as described previously.17,20–22 In brief, LA volumes were calculated using the biplane area–length method as follows: LA volume = 8 x (LA area in four-chamber view) x (LA area in two-chamber view)/3
(LA length). Left atrial length was defined as the longest line that could be drawn between the posterior LA wall and the mid-portion of the mitral annulus and was similar in the six-chamber views, which are perpendicular to each other. Left atrial maximal volume was measured at the point of mitral valve opening, whereas LA minimal volume was measured at the point of mitral valve closure. Left atrial volume at onset of atrial systole was considered the volume corresponding to the onset of the P wave in the simultaneously recorded ECG. All LA volume values were obtained as the mean of three consecutive beats and were corrected for body surface area. Left atrial systolic (active emptying) function was assessed with the use of the following parameters17,20–22:
|
| (1) |
|
| (2) |
Two-dimensional echocardiographic left atrial strain
Two-dimensional strain uses gray scale (B-Mode) sector image and is based on frame by frame tracking of small rectangular image blocks with stable speckle pattern.16 A minimum frame rate of 30 Hz was required for reliable operation of this program and frame rates of 30–90 Hz were used for routine gray-scale imaging.13–15 Apical four-chamber view value (the same of the LA volume measurement) was obtained using the same ultrasound system and probe used for standard echocardiography; end-systole was chosen as the single frame for the endocardial to epicardial region of interest to include maximal wall thickness for strain calculation. The Zoom/RES feature on the echocardiographic machine was used to improve the accuracy of atrial measurements. A circular region of interest was traced on the endocardial cavity interface of the apical four-chamber view at the LA systole (minimum cavity area) using point-and-click approach. Then a second larger concentric circle was automatically generated that was near the epicardium with a default width of 15 mm. The region of interest then included the entire LA myocardial wall, and a click feature increased or decreased the width of the two circles for thicker or thinner walls, respectively. The tracking algorithm followed the endocardium from this single frame throughout the cardiac cycle. The image-processing algorithm automatically subdivided the region of interest into blocks of
20–40 pixels containing stable patterns of speckles. Subsequent frames were then analysed automatically by searching for the new location of each of the blocks with correlation criteria and the sum of the absolute differences. The location shift of these acoustic markers from frame to frame, which represents tissue movement, provides the spatial and temporal data used to calculate velocity vectors. Temporal alterations in this stable speckle patterns are identified as moving further apart or closer together, and a series of regional strain vectors are calculated as change in length/initial length. Accordingly, for atrial longitudinal strain, myocardial thickening was represented with a positive value, colour coded as red; myocardial thinning was represented with a negative value, colour coded as blue; and then these were superimposed to conventional two-dimensional images.28–30 The software then automatically divided the image into six standard segments and provided an automated tracking score, similar to statistical standard deviation, as feedback of the stability of the regional speckle tracking, ranging from 1.0 to 3.0 in arbitrary units. A tracking score value of <2.5 was determined as acceptable as previously described, and slight adjustments were made to the placement of the region of interest in regions with greater SDs to attempt to improve the tracking stability.16
The tracking process and conversion to Lagrangian strains were performed offline using a dedicated software (EchoPAQ PC 2D strain, GE Healthcare, Milwakee). The analysis was performed for atrial longitudinal strain from the four-chamber apical view for the basal segments of LA septum and LA lateral wall, and for LA roof. Atrial radial strain was not calculated from the parasternal short- and long-axis views because the atrial wall is too thin to be properly analysed in these views, as demonstrated by our previous pilot studies.31–32 Continuous care was taken to keep the sample volume out of the pulmonary veins and of oval fossa.
Statistical methods
All the analyses were performed using a commercially available SPSS software (SPSS, Rel 11.0 2002, Chicago, SPSS Inc.). Variables are presented as mean ± SD. Two-tailed t-test for paired and unpaired data, and analysis of variance by Newman–Keuls post hoc test for multiple comparisons were used to assess changes between groups. Linear regression analyses and partial correlation test by Pearsons method were performed to assess univariate relations.
To identify significant, independent determinants of LA myocardial strain in patients with DCM, their individual association with clinical relevant and echocardiographic variables was assessed by multivariable Cox regression analysis. The following variables were included into the analysis: clinical data (age, sex, BSA, mean blood pressure, and etiology of HF), standard echocardiographic indexes (LV volumes, LA volumes, Doppler transmitral inflow measurements, mitral valve ERO, and response to CRT). These variables were selected according to their clinical relevance and potential impact on LA function, as shown by earlier studies.17 Variable selection was performed in the multivariable Cox regression as an interactive stepwise backward elimination method, each time excluding the one variable with the highest P-value according to Walds statistics. The assumption of linearity was checked graphically by studying the smoothed martingale residuals from the null model plotted against the covariate variables.33 The linearity assumptions were satisfied. The Hosmer–Lemeshow goodness-of-fit test was used to check that the model adequately fit the data.34 The model also underwent bootstrap validation (200 runs). To decrease the inflation of the Type 1 error rate because of multiple testing, the statistical significance was defined as two-sided P-value of <0.01.
Reproducibility of 2DSE measurements was determined in all the subjects. Inter- and intra-observer variability was determined using both Pearson's bivariate two-tailed correlations and Bland–Altman analysis. Relation coefficients, 95% confidence limits, and percentage errors were reported.
| Results |
|---|
|
|
|---|
Study population
Clinical features of the 90 patients enrolled in the study are shown in Table 1. At the time of CRT system implantation, 72.2% patients were in NYHA functional class III and 27.8% in NYHA class IV. QRS duration ranged from 120 to 220 ms (mean, 134.8 ± 28 ms). The two groups were comparable for most of the clinical variables. However, patients with ischaemic DCM were more frequently smokers and diabetics.
|
Standard echocardiography at baseline
Left ventricular volumes, EF, stroke volume, as well as mitral valve ERO and right ventricular TAPSE were comparable between the two groups. No significant differences were evidenced in Doppler transmitral inflow measurements and in E/Ea ratio (Table 2).
|
Left atrial dimension and function at baseline
Left atrial diameter and maximal volume were similar between the two groups. Conversely, LA active emptying volume and fraction were both lower in patients with idiopathic DCM (Table 3).
|
Overall, atrial speckle tracking was possible in 96.5% of 270 attempted segments from the 90 subjects with technically adequate images, with only 3.5% of segments eliminated with tracking variation scores of >2.5. Overall tracking variation scores were <2.0 in 77% segments.
Peak systolic myocardial atrial strain was significantly compromised in patients with idiopathic DCM compared with ischaemic DCM in all the analysed atrial segments (Table 3; Figures 1 and 2). Comparing standard echo indexes with 2DSE values, we found a significant linear correlation between LA emptying fraction peak and systolic myocardial atrial strain in LA lateral wall (r = 0.71, P = 0.001), atrial septum (r = 0.68, P = 0.002), and LA roof (r = 0.63, P = 0.005).
|
|
Twenty-five age- and sex-matched healthy subjects (15 males; mean age: 52.1 ± 4.2 years), without history of cardiovascular disease, with normal ECG (mean QRS duration, 81 ± 12 ms), and normal EF (mean EF, 67.2 ± 8.1%) were selected as control group for LA myocardial function analysis. The mean values of LA peak systolic strain (in percentage) in the control group were significantly higher than those in HF patients in all the analysed atrial segments (Table 3; Figure 3).
|
Effects of 6 month cardiac resynchronization therapy
All the patients had echocardiographic studies 6 months after CRT to determine long-term response. When LV end-systolic volume was used as the principal outcome marker for response to CRT, 64 patients (71.1%) (37 with idiopathic DCM and 27 with ischaemic DCM) were long-term responders, and 26 (28.9%) (10 with idiopathic DCM; 16 with ischaemic DCM) were non-responders.
Overall, significant changes in the entire CRT group were observed. A significant increase in LV EF and LV stroke volume was detected; however, mitral regurgitation decreased. In addition, standard and 2DSE indexes of LA systolic function increased. Individual responses, however, were variable.
In particular, among the responders to CRT, patients with ischaemic DCM showed a significant improvement in both LA active emptying volume and fraction, as well as in LA peak systolic strain in all the analysed myocardial segments (Table 4). Significant inverse correlation between changes in myocardial LA lateral strain and both changes in mitral valve ERO (r = –0.58, P < 0.001) and changes in transmitral E/A ratio (r = –0.52, P < 0.005) were observed. Conversely, in patients with idiopathic DCM, no significant differences were found before and after CRT for both standard and 2DSE indexes of LA systolic function (Table 4).
|
On the other hand, in patients non-responders to CRT, no improvement in LA function was evidenced both in patients with idiopathic and with ischaemic aetiology of DCM (Table 5). In addition, no significant correlation was found between changes in LA myocardial strain and changes in either mitral valve ERO (r = –0.2, P = 0.23) and in transmitral E/A ratio (r = –0.22, P = 0.25).
|
By multivariable analysis, after adjusting for potential determinants such as age, heart rate, BSA, mean blood pressure, LV and LA volumes, in the overall population of patients, ischaemic aetiology of DCM (ß-coefficient = 0.62; P < 0.0001) and positive response to CRT (ß-coefficient = 0.48; P < 0.001) are emerged as the only independent determinants of LA lateral wall systolic strain after CRT.
Reproducibility of speckle-tracking two-dimensional strain echocardiography left atrial strain
Intra-observer variability
Pearson's correlations—lateral wall: r = 0.87, P < 0.00001; atrial septum: r = 0.88, P < 0.00001; atrial roof: r = 0.85, P < 0.00001.
Bland–Altman analysis—lateral wall (95% CI ±1.8; percent error 3.3%); atrial septum (95% CI ±1.2; percent error 3.1%); atrial roof (95% CI ±1.6; percent error 3.2%).
Inter-observer variability
Pearson's correlations—lateral wall: r = 0.85; P < 0.00001; atrial septum: r = 0.86; P < 0.00001; atrial roof: r = 0.84; P < 0.00001.
Bland–Altman analysis—lateral wall (95% CI ± 1.9; percent error 3.4%), atrial septum (95% CI ±1.5; percent error 3.5%), atrial roof (95% CI ±1.7; percent error 3.6%).
| Discussion |
|---|
|
|
|---|
The results of the present study demonstrate the usefulness of 2DSE in analysing LA myocardial function in patients with either idiopathic or ischaemic DCM undergoing CRT.
The main findings of our study are: (i) 2DSE detected at baseline more impaired LA systolic function in patients with idiopathic compared with ischaemic DCM; (ii) a significant improvement in LA systolic function was obtained only in patients with ischaemic DCM responders to CRT; (iii) ischaemic aetiology of DCM and positive response to CRT were powerful independent determinants of LA myocardial function after CRT.
The advantages and perspectives of speckle-tracking two-dimensional strain echocardiography in the study of left atrial function
Previous reports have documented that echocardiography with its new applications such as DMI and strain imaging may represent useful techniques to better select candidate to CRT.3,8–12
In particular, myocardial strain is a dimensionless index of change in myocardial length in response to an applied force, and is expressed as fractional or percent change. This technique has a theoretic advantage over DMI in that it is relatively immune to cardiac translational motion and tethering allowing differentiation between active systolic contraction and passive motion.10 In a previous report, we have documented for the first time the feasibility of Doppler strain in analysing myocardial deformation properties of both atria during recent-onset lone atrial fibrillation, and atrial myocardial strain measurements as predictive of sinus rhythm maintenance.31
However, in patients with DCM and large and spherical atria and ventricles, Doppler strain estimation with its angle dependency represents an important limitation as a result of poor alignment between Doppler beam and myocardial wall.10 Conversely, in our study, we used 2DSE, a novel approach to quantify regional myocardial deformation within a scan plane that is inherently two-dimensional and independent of interrogation angle as it tracks speckle patterns (acoustic markers) within serial B-Mode sector scans.28–30 A recent study in patients with DCM has demonstrated that asynchrony indexes derived from the novel speckle-tracking radial 2DSE predict response to CRT.16 To the best of our knowledge, the present paper is the first report analysing the LA function in patients with DCM at baseline and after CRT using speckle-tracking 2DSE.
Left atrial function in idiopathic vs. ischaemic dilated cardiomyopathy
Normal LA function consists of three components. Left atrium is: (i) a contractile chamber (booster pump) that actively empties immediately before the end of LV diastole and establishes final LV end-diastolic volume; (ii) a reservoir that stores pulmonary venous return during LV systole and isovolumic relaxation after the closure and before the opening of the mitral valve; (iii) a conduit that empties its contents into the LV down a pressure gradient after the mitral valve opens and during LV diastole.17
The LA is therefore exposed during ventricular diastole to the pressures of the LV. With increased stiffness or non-compliance of the LV, LA pressure rises to maintain adequate LV filling, and the increased atrial wall tension leads to chamber dilatation and stretch of the atrial myocardium.17 As a result, early in HF, LA pump function is augmented but LA stiffness increases and work mismatch occurs. With further progression of LV dysfunction, LA pump function decreases as a result of increased afterload imposed on the LA myocardium. However, in idiopathic DCM, although loading conditions are same, a more depressed LA booster pump function at rest has been observed compared with ischaemic patients, and attributed both to altered LA overload and to LA larger involvement in the myopathic process.20–22,35
Several authors have recently demonstrated that LA peak systolic myocardial atrial strain, which is LA passive stretching during LV systole, could be used as an index of LA reservoir function.31–32,36–39
In our study protocol, both LA booster pump function (i.e. LA active emptying fraction) and LA passive deformation indexes (i.e. peak systolic myocardial atrial strain) in patients with DCM were significantly more impaired in patients with idiopathic DCM compared with patients with ischaemic DCM in all the analysed atrial segments.
Since heart rate, LA volumes, transmitral flow velocity profile, LV isovolumic relaxation time, LV stroke volume, and mitral valve ERO in our population were similar in the two cardiomyopathy groups, it is unlikely that the differences between the two groups were only due to differences in LA loading conditions. Therefore, it is reasonable to assume that in addition to LV diastolic dysfunction, depressed LA contractility and deformation resulting from LA involvement in the myopathic process might have contributed to the LA myocardial dysfunction in our patients with idiopathic DCM. This hypothesis is further supported by previous studies demonstrating: (i) higher degree of LA fibrotic change in autopsied heart of patients with idiopathic DCM compared with patients with old myocardial infarction35; (ii) more impaired LA function in patients with idiopathic DCM than in those with LV dysfunction caused by valvular abnormalities, despite comparable mechanical loads40; (iii) reduced LA pump function in idiopathic vs. ischaemic DCM despite a similar LA size and tension at onset of LA systole20; (iv) after dobutamine infusion, absence of significant changes in the LA active emptying fraction in idiopathic DCM, in contrast to significant increase in these indexes in ischaemic DCM.21–22
Left atrial myocardial function after 6 month cardiac resynchronization therapy
Few reports have previously investigated the long-term effect of CRT on LA volumes and function by standard echocardiography. In a first study by Vural et al.41 in 20 patients with DCM (13 idiopathic/seven ischaemic), CRT resulted in atrial reverse remodelling (reduced both maximal and minimal volume), increase in LA total emptying fraction, and reduction in both frequency and intensity of LA spontaneous echo contrast. In a subsequent study of the same group by transoesophageal echocardiography, in 18 patients (10 idiopathic/eight ischaemic), treatment of HF by CRT resulted with marked improvements in LA appendage function and increases pulmonary venous systolic velocity.42
In our study, a significant improvement in LA myocardial function after CRT was obtained only in patients with ischaemic DCM responders to CRT. Such reversible LA dysfunction in ischaemic patients may further confirm that the initial LA dysfunction was mainly because of LA afterload mismatch rather than intrinsic LA disease. In fact, in our patients, changes in lateral LA strain were strongly associated with changes in mitral valve ERO and in Doppler transmitral A wave after CRT. Therefore, in these patients, enhanced LA electromechanical function after CRT may result in improved LV filling, contribute to the haemodynamic improvement observed with CRT, and also has implications for reduction in the propensity for atrial fibrillation.17
Conversely, in idiopathic DCM, our findings detected a reduced LA systolic reserve as well as an irreversible impairment of LA deformation even after improvement of loading conditions during CRT. The lack of association between changes in LA strain and improvement in diastolic measurements may confirm the presence in idiopathic DCM patients of a larger involvement of LA structure and function in the myopathic process.35 Since in congestive HF the loss of LA contribution to LV filling adversely affects prognosis, and a recent report has used LA strain analysis for demonstrating beneficial effects of renin–angiotensin system inhibitors on LA reservoir function,39 our results may have important clinical implications on the evaluation long-term effects of pharmacological therapy and of CRT in patients with DCM.
| Study limitations |
|---|
|
|
|---|
Our study has several limitations. First, even though the biplane area–length method is a better way to assess LA dimensions compared with the M-mode evaluation,17 it is necessary to underline that the volumes obtained with this method are estimations derived from two-dimensional datasets, and not true measurements of a three-dimensional structure.
Atrial strain is influenced by loading conditions; thus, the comparisons of atrial deformation properties between idiopathic and ischaemic patients may have been influenced by preload and by different atrial wall stress.10 However, even though our study protocol did not include any invasive haemodynamic measurement, our patients with DCM presented similar mean heart rates, LA and LV volumes, and Doppler indexes of global diastolic function. As a result, this confounding factor should not detract from the value of atrial strain analysis.
A technical limitation is that speckle-tracking 2DSE is dependent on frame rate as well as image resolution. Low frame rate results in the speckle pattern changing too much from frame to frame, which prevents the precise characterization of regional myocardial motion and impacts the overall temporal resolution of the regional strain map. In contrast, increase in the frame rate reduces the scan line density, thereby reducing the image resolution.28–30 Frame rate in our setting ranged from 35 to 90 frame/s; this value is lower than frame rate available with Doppler strain, however indexes of LA function used in this study did not rely on difference in the timing of contraction.
As previously specified, this is the first report analysing LA function in patients with DCM using speckle-tracking 2DSE. In particular, we used EchoPAC PC program (GE Healthcare) for two-dimensional strain analysis. Since a definite 2DSE software for LA two-dimensional strain has not yet been provided, we applied a program of 2DSE for LV to analyse also LA strain. The feasibility and the reproducibility of LA strain patterns and measurements were good. However, in future studies, changes in the software may be needed in order to improve the tracking ability of the speckle-tracking system for LA functional study.
Our analysis was limited by data collection from the apical four-chamber view as this was the only image plane that provided reliable information on longitudinal atrial deformation in a sufficient numbers of patients. It would be desirable to obtain information on left atrial deformation from the complete atrial cavity. However, in the present study, atrial radial strain was not calculated from the parasternal short- and long-axis views because the atrial wall is too thin to be properly analysed in these views, as demonstrated by our previous studies.31–32
In our patients with ischaemic DCM, end-diastolic myocardial thickness of the region targeted for LV pacing was measured by standard echocardiography. This is a rude but simple and largely available method to detect myocardial scar tissue in ischaemic patients, as previously showed by Cwajg JM et al.,23 who reported that end-diastolic wall thickness can predict in patients with suspected myocardial hibernation recovery of function similar to T1-201 scintigraphy. Recent reports have documented that delayed contrast-enhanced magnetic resonance imaging, positron emission tomography, and gated SPECT are highly effective and more quantitative methods to define myocardial viability.43–45 However, in our study, we limited the use of contrast-enhanced cardiac magnetic resonance imaging to patients with ischaemic DCM and normal myocardial thickness. Although these techniques may serve as gold standard, high costs and limited availability reduce their routine use in our country in selecting candidates for CRT.
| Conclusions |
|---|
|
|
|---|
This study demonstrates for the first time that a novel speckle-tracking algorithm applied to routine gray-scale two-dimensional images represents a promising and feasible non-invasive and easy-repeatable technique to assess LA atrial myocardial function in patients with DCM. Left atrial pump function and passive deformation both at baseline and after CRT are more impaired in idiopathic compared with ischaemic DCM. Since LA enlargement has been proposed as a barometer of diastolic burden and a predictor of common cardiovascular outcomes such as atrial fibrillation, stroke, and cardiovascular death,17 future longitudinal studies are warranted to understand further the natural history of LA myocardial function, the extent of reversibility of LA dysfunction with medical therapy and/or CRT, and the possible long-term impact of such changes on outcomes in patients with congestive HF.
| Supplementary material |
|---|
|
|
|---|
Supplementary material is available at European Heart Journal online.
| Acknowledgements |
|---|
|
|
|---|
The authors are grateful to Mr Marco Maglione (GE Medical Systems, Italy) for excellent technical support to the study protocol. The authors would also like to express their gratitude for Ms Gemma D'Andrea and Mrs Francesca Romana D'Andrea for being a valuable source of ideas and motivation as well.
Conflict of interest: none declared.
| References |
|---|
|
|
|---|
- Abraham WT, Fisher WG, Smith AL, Delurgio DB, Leon AR, Loh E, Kocovic DZ, Packer M, Clavell AL, Hayes DL, Ellestad M, Trupp RJ, Underwood J, Pickering F, Truex C, McAtee P, Messenger J, MIRACLE Study Group. Multicenter InSync randomized clinical evaluation. Cardiac resynchronization in chronic heart failure. N Engl J Med (2002) 346:1845–1852.
[Abstract/Free Full Text] - Cleland JG, Daubert JC, Erdmann E, Freemantle N, Gras D, Kappenberger L, Tavazzi L, Cardiac Resynchronization-Heart Failure (CARE-HF) Study Investigators. The effect of cardiac resynchronization on morbidity and mortality in heart failure. N Engl J Med (2005) 352:1539–1549.
[Abstract/Free Full Text] - Bax JJ, Abraham T, Barold SS, Breithardt OA, Fung JW, Garrigue S, Gorcsan J III, 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: part 1—issues before device implantation. J Am Coll Cardiol (2005) 20:2153–2167.
- Hawkins NM, Petrie MC, MacDonald MR, Hogg KJ, McMurray JJ. Selecting patients for cardiac resynchronization therapy: electrical or mechanical dyssynchrony? Eur Heart J (2006) 27:1270–1281.
[Abstract/Free Full Text] - D'Andrea A, Ducceschi V, Caso P, Galderisi M, Mercurio B, Liccardo B, Sarubbi B, Scherillo M, Cotrufo M, Calabro R. Usefulness of Doppler tissue imaging for the assessment of right and left ventricular myocardial function in patients with dual-chamber pacing. Int J Cardiol (2001) 81:75–83.[CrossRef][Web of Science][Medline]
- D'Andrea A, Caso P, Galderisi M, Ducceschi V, Scherillo M, Sarubbi B, Mininni N, Calabro R. Ventricular interdependence in patients with dual-chamber pacing: a Doppler tissue imaging study. Echocardiography (2002) 19:289–297.[CrossRef][Web of Science][Medline]
- D'Andrea A, Caso P, Cuomo S, Scarafile R, Salerno G, Limongelli G, Di Salvo G, Severino S, Ascione L, Calabro P, Romano M, Romano G, Santangelo L, Maiello C, Cotrufo M, Calabro R. Effect of dynamic myocardial dyssynchrony on mitral regurgitation during supine bicycle exercise stress echocardiography in patients with idiopathic dilated cardiomyopathy and narrow QRS. Eur Heart J (2007) 28:1004–1011.
[Abstract/Free Full Text] - Yu C-M, Fung JW-H, Zhang Q, Chan C-K, Chan Y-S, Lin H, Kum LCC, Kong S-L, Zhang Y, Sanderson JE. Tissue Doppler imaging is superior to strain rate imaging and postsystolic shortening on the prediction of reverse remodeling in both ischemic and nonischemic heart failure after cardiac resynchronization therapy. Circulation (2004) 110:66–73.
[Abstract/Free Full Text] - D'Andrea A, Severino S, Caso P, Fusco A, Lo Piccolo R, Liccardo B, Forni A, Di Salvo G, Scherillo M, Mininni N, Calabro R. Prognostic value of intra-left ventricular electromechanical asynchrony in patients with hypertrophic cardiomyopathy. Eur Heart J (2006) 27:1311–1318.
[Abstract/Free Full Text] - Sutherland GR, Di Salvo G, Claus P, D'hooge J, Bijnens B. Strain and strain rate imaging: a new clinical approach to quantifying regional myocardial function. J Am Soc Echocardiogr (2004) 17:788–802.[CrossRef][Web of Science][Medline]
- Yu CM, Zhang Q, Wing-Hong Fung J, Chi-Kin Chan H, Chan YS, Wai-Kwok Yip G, Kong SL, Lin H, Zhang Y, Sanderson JE. A novel tool to assess systolic asynchrony and identify responders of cardiac resynchronization therapy by tissue synchronization imaging. J Am Coll Cardiol (2005) 45:677–684.
[Abstract/Free Full Text] - Mele D, Pasanisi G, Capasso F, De Simone A, Morales MA, Poggio D, Capucci A, Tabacchi G, Sallusti L, Ferrari R. Left intraventricular myocardial deformation dyssynchrony identifies responders to cardiac resynchronization therapy in patients with heart failure. Eur Heart J (2006) 27:1070–1078.
[Abstract/Free Full Text] - Langeland S, D'hooge J, Wouters PF, Leather HA, Claus P, Bijnens B, Sutherland GR. Experimental validation of a new ultrasound method for the simultaneous assessment of radial and longitudinal myocardial deformation independent of insonation angle. Circulation (2005) 112:2157–2162.
[Abstract/Free Full Text] - Notomi Y, Lysyansky P, Setser RM, Shiota T, Popovic ZB, Martin-Miklovic MG, Weaver JA, Oryszak SJ, Greenberg NL, White RD, Thomas JD. Measurement of ventricular torsion by two-dimensional ultrasound speckle tracking imaging. J Am Coll Cardiol (2005) 45:2034–2041.
[Abstract/Free Full Text] - Amundsen BH, Helle-Valle T, Edvardsen T, Torp H, Crosby J, Lyseggen E, Stoylen A, Ihlen H, Lima JA, Smiseth OA, Slordahl SA. Noninvasive myocardial strain measurement by speckle tracking echocardiography: validation against sonomicrometry and tagged magnetic resonance imaging. J Am Coll Cardiol (2006) 47:789–793.
[Abstract/Free Full Text] - Suffoletto MS, Dohi K, Cannesson M, Saba S, Gorcsan J 3rd. Novel speckle-tracking radial strain from routine black-and-white echocardiographic images to quantify dyssynchrony and predict response to cardiac resynchronization therapy. Circulation (2006) 113:960–968.
[Abstract/Free Full Text] - Abhayaratna WP, Seward JB, Appleton CP, Douglas PS, Oh JK, Tajik AJ, Tsang TS. Left atrial size: physiologic determinants and clinical applications. J Am Coll Cardiol (2006) 47:2357–2363.
[Abstract/Free Full Text] - Triposkiadis F, Pitsavos C, Boudoulas H, Trikas A, Toutouzas P. Left atrial myopathy in idiopathic dilated cardiomyopathy. Am Heart J (1994) 128:308–315.[CrossRef][Web of Science][Medline]
- Trikas A, Triposkiadis F, Pitsavos C, Tentolouris K, Kyriakidis M, Gialafos J, Toutouzas P. Relation of left atrial volume and systolic function to the hormonal response in idiopathic dilated cardiomyopathy. Int J Cardiol (1994) 47:139–143.[CrossRef][Web of Science][Medline]
- Triposkiadis F, Moyssakis I, Hadjinikolaou L, Makris T, Zioris H, Hatzizaharias A, Kyriakidis M. Left atrial systolic function is depressed in idiopathic and preserved in ischemic dilated cardiomyopathy. Eur J Clin Invest (1999) 29:905–912.[CrossRef][Web of Science][Medline]
- Moyssakis I, Papadopoulos DP, Kelepeshis G, Gialafos E, Votteas V, Triposkiadis F. Left atrial systolic reserve in idiopathic vs. ischaemic-dilated cardiomyopathy. Eur J Clin Invest (2005) 35:355–361.[CrossRef][Web of Science][Medline]
- Paraskevaidis I, Doduras T, Adamopoulos S, Kremastinos D. Left atrial functional reserve in patients with nonischaemic dilated cardiomyopathy. Chest (2002) 122:1340–1347.[CrossRef][Web of Science][Medline]
- Cwajg JM, Cwajg E, Nagueh SF, He ZX, Qureshi U, Olmos LI, Quinones MA, Verani MS, Winters WL, Zoghbi WA. End-diastolic wall thickness as a predictor of recovery of function in myocardial hibernation: relation to rest-redistribution T1-201 tomography and dobutamine stress echocardiography. J Am Coll Cardiol (2000) 35:1152–1161.
[Abstract/Free Full Text] - Ritter P, Padeletti L, Gillio-Meina L, Gaggini G. Determination of the optimal atrioventricular delay in DDD pacing. Comparison between echo and peak endocardial acceleration measurements. Europace (1999) 1:126–130.
[Abstract/Free Full Text] - Dubin J, Wallerson DC, Cody RJ, Devereux RB. Comparative accuracy of Doppler echocardiographic methods for clinical stroke volume determination. Am Heart J (1990) 120:116–123.[CrossRef][Web of Science][Medline]
- Nagueh SF, Sun H, Kopelen HA, Middleton KJ, Khoury DS. Hemodynamic determinants of the mitral annulus diastolic velocities by tissue Doppler. J Am Coll Cardiol (2001) 37:278–285.
[Abstract/Free Full Text] - Enriquez-Sarano M, Miller F, Hayes S, Bailey K, Tajik A, Seward J. Effective mitral regurgitant orifice area: clinical use and pitfalls of the proximal isovelocity surface area method. J Am Coll Cardiol (1995) 25:703–709.[Abstract]
- Liang HY, Cauduro S, Pellikka P, Wang J, Urheim S, Yang EH, Rihal C, Belohlavek M, Khandheria B, Miller FA, Abraham TP. Usefulness of two-dimensional speckle strain for evaluation of left ventricular diastolic deformation in patients with coronary artery disease. Am J Cardiol (2006) 98:1581–1586.[CrossRef][Web of Science][Medline]
- Cho GY, Chan J, Leano R, Strudwick M, Marwick TH. Comparison of two-dimensional speckle and tissue velocity based strain and validation with harmonic phase magnetic resonance imaging. Am J Cardiol (2006) 97:1661–1666.[CrossRef][Web of Science][Medline]
- Ashraf M, Li XK, Young MT, Jensen AJ, Pemberton J, Hui L, Lysyansky P, Friedman Z, Park B, Sahn DJ. Delineation of cardiac twist by a sonographically based 2-dimensional strain analysis method: an in vitro validation study. J Ultrasound Med (2006) 25:1193–1198.
[Abstract/Free Full Text] - Di Salvo G, Caso P, Lo Piccolo R, Fusco A, Martiniello AR, Russo MG, D'Onofrio A, Severino S, Calabro P, Pacileo G, Mininni N, Calabro R. Atrial myocardial deformation properties predict maintenance of sinus rhythm after external cardioversion of recent-onset lone atrial fibrillation: a color Doppler myocardial imaging and transthoracic and transesophageal echocardiographic study. Circulation (2005) 112:387–395.
[Abstract/Free Full Text] - Di Salvo G, Drago M, Pacileo G, Rea A, Carrozza M, Santoro G, Bigazzi MC, Caso P, Russo MG, Carminati M, Calabro R. Atrial function after surgical and percutaneous closure of atrial septal defect: a strain rate imaging study. J Am Soc Echocardiogr (2005) 18:930–933.[CrossRef][Web of Science][Medline]
- Grambsch PM, Therneau TM, Fleming TR. Diagnostic plots to reveal functional form for covariates in multiplicative intensity models. Biometrics (1995) 51:1469–1482.[CrossRef][Web of Science][Medline]
- Weitzen S, Lapane KL, Toledano AY, Hume AL, Mor V. Principles for modeling propensity scores in medical research: a systematic literature review. Pharmacoepidemiol Drug Saf (2004) 13:855–857.[CrossRef][Web of Science][Medline]
- Ohtani K, Yutani C, Nagata S, Koretsune Y, Hori M, Kamada T. High prevalence of atrial fibrosis in patients with dilated cardiomyopathy. J Am Coll Cardiol (1995) 25:1162–1169.[Abstract]
- Inaba Y, Yuda S, Kobayashi N, Hashimoto A, Uno K, Nakata T, Tsuchihashi K, Miura T, Ura N, Shimamoto K. Strain rate imaging for noninvasive functional quantification of the left atrium: comparative studies in controls and patients with atrial fibrillation. J Am Soc Echocardiogr (2005) 18:729–736.[CrossRef][Web of Science][Medline]
- Sirbu C, Herbots L, D'hooge J, Claus P, Marciniak A, Langeland T, Bijnens B, Rademakers FE, Sutherland GR. Feasibility of strain and strain rate imaging for the assessment of regional left atrial deformation: a study in normal subjects. Eur J Echocardiogr (2006) 7:199–208.
[Abstract/Free Full Text] - Wang T, Wang M, Fung JWH, Yip GWK, Zhang Y, Ho PPY, Tse DMK, Yu CM, Sanderson JE. Atrial strain rate echocardiography can predict success or failure of cardioversion for atrial fibrillation: a combined transthoracic tissue Doppler and transesophageal imaging study. Int J Cardiol (2007) 114:202–209.[CrossRef][Web of Science][Medline]
- Kokubu N, Yuda S, Tsuchihashi K, Hashimoto A, Nakata T, Miura T, Ura N, Nagao K, Tsuzuki M, Wakanayashi C, Shimamoto K. Noninvasive assessment of left atrial function by strain rate imaging in patients with hypertension: a possible beneficial effect of Renin–Angiotensin system inhibition on left atrial function. Hypertens Res (2007) 30:13–21.[CrossRef][Web of Science][Medline]
- Triposkiadis F, Pitsavos C, Boudoulas H, Trikas A, Kallikazaros I, Stefanadis C. Left atrial function in valvular aortic stenosis. J Heart Valve Dis (1993) 2:104–113.[Medline]
- Vural A, Agacdiken A, Ural D, Sahin T, Kozdag G, Kahraman G, Ural E, Akbas H, Suzer K, Komsuoglu B. Effect of cardiac resynchronization therapy on left atrial reverse remodeling and spontaneous echo contrast. Tohoku J Exp Med (2004) 202:143–153.[CrossRef][Web of Science][Medline]
- Vural A, Agacdiken A, Ural D, Sahin T, Kozdag G, Kahraman G, Ural E, Akbas H, Suzer K, Komsuoglu B. Effect of cardiac resynchronization therapy on left atrial appendage function and pulmonary venous flow pattern. Int J Cardiol (2005) 102:103–109.[CrossRef][Web of Science][Medline]
- Bleeker GB, Kaandorp TA, Lamb HJ, Boersma E, Steendijk P, de Roos A, van der Wall EE, Schalij MJ, Bax JJ. Effect of posterolateral scar tissue on clinical and echocardiographic improvement after cardiac resynchronization therapy. Circulation (2006) 113:969–976.
[Abstract/Free Full Text] - De Winter O, Van de Veire N, Van Heuverswijn F, Van Pottelberge G, Gillebert TC, De Sutter J. Relationship between QRS duration, left ventricular volumes and prevalence of nonviability in patients with coronary artery disease and severe left ventricular dysfunction. Eur J Heart Fail (2006) 8:275–277.[CrossRef][Web of Science][Medline]
- Ypenburg C, Schalij MJ, Bleeker GB, Steendijk P, Boersma E, Dibbets-Schneider P, Stokkel MP, van der Wall EE, Bax JJ. Extent of viability to predict response to cardiac resynchronization therapy in ischemic heart failure patients. J Nucl Med (2006) 47:1565–1570.
[Abstract/Free Full Text]
This article has been cited by other articles:
![]() |
S Eshoo, A C Boyd, D L Ross, T H Marwick, and L Thomas Strain rate evaluation of phasic atrial function in hypertension Heart, July 15, 2009; 95(14): 1184 - 1191. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Caso, R. Ancona, G. Di Salvo, S. Comenale Pinto, M. Macrino, V. Di Palma, A. D'Andrea, A.R. Martiniello, S. Severino, and R. Calabro Atrial reservoir function by strain rate imaging in asymptomatic mitral stenosis: prognostic value at 3 year follow-up Eur J Echocardiogr, May 13, 2009; (2009) jep058v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Kurt, J. Wang, G. Torre-Amione, and S. F. Nagueh Left Atrial Function in Diastolic Heart Failure Circ Cardiovasc Imaging, January 1, 2009; 2(1): 10 - 15. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. D'Andrea, P. Caso, R. Scarafile, L. Riegler, G. Salerno, F. Castaldo, R. Gravino, R. Cocchia, L. Del Viscovo, G. Limongelli, et al. Effects of global longitudinal strain and total scar burden on response to cardiac resynchronization therapy in patients with ischaemic dilated cardiomyopathy Eur J Heart Fail, January 1, 2009; 11(1): 58 - 67. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Jarnert, A. Melcher, K. Caidahl, H. Persson, L. Ryden, and M. J. Eriksson Left atrial velocity vector imaging for the detection and quantification of left ventricular diastolic function in type 2 diabetes Eur J Heart Fail, November 1, 2008; 10(11): 1080 - 1087. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Valzania, M. J. Eriksson, G. Boriani, and F. Gadler Cardiac resynchronization therapy during rest and exercise: comparison of two optimization methods Europace, October 1, 2008; 10(10): 1161 - 1169. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||







