The use of strain, strain rate, and displacement by 2D speckle tracking for assessment of systolic left ventricular function in goats: applicability and influence of general anesthesia
Cardiovascular Ultrasound volume 13, Article number: 11 (2015)
Assessment of left ventricular (LV) systolic function can be achieved by conventional echocardiographic methods, but quantification of contractility, regional myocardial function, and ventricular synchrony is challenging. The goal of this study was to investigate the applicability of two-dimensional speckle tracking (2DST) to characterize segmental and global wall motion for assessment of LV function and LV synchrony in healthy goats. We aimed to describe the techniques, report normal values of a variety of 2DST indices, and determine the influence of general anesthesia.
Prospective study on 22 healthy female Saanen goats (3.7 ± 1.1 y, 60.2 ± 10.5 kg [mean ± SD]). All goats underwent two transthoracic echocardiographic examinations, the first standing and unsedated and the second 7.4 ± 3.5 days later during isoflurane anesthesia and positioned in sternal recumbency. Data analyses were performed offline, blinded, and in random order. Left ventricular longitudinal, radial and circumferential strain and strain rate as well as longitudinal and radial displacement were measured using 2DST methods. Summary statistics were generated and differences of 2DST variables between myocardial segments and treatments (i.e., awake vs. anesthetized) were assessed statistically (alpha level=0.05).
Echocardiographic analyses by 2DST were feasible in all goats and at both time points. Longitudinal systolic strain, strain rate and displacement followed a gradient from apex to base. Absolute systolic strain was generally lower and strain rate was higher in awake goats compared to anesthetized goats. Circumferential and radial indices did not consistently follow a segmental pattern. Generally, peak strain occurred later in anesthetized goats compared to awake goats. General anesthesia did not significantly influence LV synchrony.
2SDT is a valid method for non-invasive characterization of LV wall motion in awake and anesthetized goats. The results of this study add to the understanding of LV mechanical function, aid in the diagnosis of global and segmental LV systolic dysfunction, and will be useful for future cardiovascular studies in this species. However, effects of anesthesia and species-specific characteristics should be considered when goats are used as animal models for human disease.
The echocardiographic assessment of myocardial function is a pivotal part of every clinical or experimental cardiologic examination, but still presents a considerable challenge for every cardiologist. Traditionally, myocardial function has been assessed by use of two-dimensional (2D), M-Mode, and Doppler echocardiography [1,2]. Although these methods are well established, they serve to assess overall systolic function, but do not provide specific data for quantification of myocardial contractility, regional myocardial function, or ventricular synchrony. At best, conventional echocardiographic methods can be used for subjective visual assessment of regional wall motion or manual tracking of myocardial movement, which requires considerable experience and is highly operator dependent [3-6].
Echocardiographic strain and strain rate imaging using two-dimensional speckle tracking (2DST) has been advocated for assessment of regional and global ventricular function beyond the conventional echocardiographic approach [7-9]. The high sensitivity of this method for early detection of myocardial dysfunction or dyssynchrony related to coronary artery disease, myocardial ischemia, myocardial infarction, or heart failure render this non-invasive diagnostic method suitable for clinical use in humans . It has been used in a variety of species, including humans, dogs, cats, horses, and pigs [4,10-16]. The goat is commonly used as an animal model for the study of human cardiovascular disease [17-20]. Therefore, novel echocardiographic techniques are also pertinent to this species, despite the fact that natural cardiac disease is of minor relevance in goats . However, to our knowledge there are no comprehensive studies investigating the use of strain and strain rate by 2DST in goats.
Therefore, the goal of this study was to investigate the applicability of 2DST to characterize segmental and global wall motion for assessment of LV function and LV synchrony in healthy Saanen goats. We hypothesized that 2DST can be applied to characterize myocardial function in standing, awake goats as well as in anesthetized goats during general anesthesia. We aimed to describe the techniques, report normal values of a variety of 2DST indices, and determine the influence of general anesthesia. We further intended to explore the agreement of the findings of this study to other studies conducted in people and in other animal species to identify similarities and disparities with regard to the physiology of LV mechanics. The results of this study provide fundamental information on the use of 2DST for assessment of LV mechanics in goats and will be useful for future cardiovascular studies in this species.
22 female Saanen goats aged 3.7 ± 1.1 [mean ± SD] years and with a body weight of 60.2 ± 10.5 kg were studied prospectively. All goats were considered healthy based upon physical examination, cardiac auscultation, and routine echocardiographic examination .
None of the goats received medications during the 2 weeks preceding entry into the study. The goats were acclimatized to the hospital for one week before the study. They were housed indoors, kept on straw and had free access to water and hay. Animal experiments were carried out in accordance with the Swiss law on animal protection. All animal studies were approved by the district veterinary office of the Canton of Zurich.
All goats underwent a complete echocardiographic examination while standing in a quiet room, unsedated and restrained by an experienced handler. A second echocardiographic examination within 7.4 ± 3.5 days of the first examination was performed immediately following a computed tomography scan conducted within the scope of another investigation, with the goats in general anesthesia and positioned in sternal recumbency. The mean time from induction of anesthesia to the start of the echocardiographic examination was 68 (52–103) min [mean (range)].
The animals were fasted for 24 h and deprived of water for 2 h prior to anaesthesia. A 14 G/1.88 in catheter (BD Angiocath, Becton Dickinson, Allschwil, Switzerland) was placed aseptically into the right jugular vein. All animals were premedicated intravenously with 0.1 mg/kg xylazine (Rompun 2%, Provet, Lyssach, Switzerland), diluted in 20 mL of saline (NaCl 0.9%, Braun Medical, Sempach, Switzerland) and given over 5 min with a syringe pump (Syramed μ6000, Arcomed, Regensdorf, Switzerland). Anaesthesia was then induced with either 3 mg/kg racemic ketamine (Narketan 10%, Vetoquinol, Ittigen, Switzerland) or 1.5 mg/kg S-ketamine (Keta-S 6%, Dr. Graeub AG, Bern, Switzerland) i.v. given by manual injection.
Once the animals were recumbent, they were positioned sternally, intubated, and connected to an anesthesia machine. Generally, the goats were allowed to breathe spontaneously throughout the entire duration of anesthesia. Mechanical ventilation was only applied when arterial blood gas analysis revealed a PaCO2 > 50 mmHg. In this case, intermittent positive pressure ventilation was used with a tidal volume of 10-15 mL/kg at a rate of 10-15 breaths per minute. Anesthesia was maintained with isoflurane (IsoFlo®, Abbott, Baar, Switzerland) delivered in oxygen and air via a semi-closed circle absorption system with a flow of 60 mL/kg/min during the first 10 min and a flow of 35 mL/kg/min thereafter. The vaporizer was initially set to 2.5% and after 15 min it was reduced and held between 1.5 and 2% to maintain an end-tidal isoflurane concentration of 1.1%. Lactated Ringer’s solution (Ringer-Laktat, Fresenius Kabi, Stans, Switzerland) was administered i.v. at a dose rate of 10 mL/kg/h. During anesthesia, the goats were kept in sternal recumbency.
Each goat was instrumented with an electrocardiogram, a pulse oximeter and an end-tidal gas analyzer. A portable multiparameter monitor (Datex-Ohmeda AS/3 Compact Monitor; Anandic, Schaffhausen, Switzerland) continuously displayed the physiologic monitoring data during the anesthetic period. The ear artery was cannulated for intraarterial blood pressure monitoring and blood sampling for blood gas analyses. Heart rate (HR), respiratory rate, saturation of arterial oxygen, direct arterial blood pressures, end-tidal CO2, end-tidal isoflurane, and rectal temperature were monitored. Blood gas analysis was performed at 10 minutes after anesthesia induction and then every 30 minutes for the duration of anesthesia. Following anesthesia, the goats were placed sternally in a padded box. The endotracheal tube was removed once swallowing reflex returned.
Transthoracic echocardiography was performed using a high-end digital echocardiograph (GE Vivid 7 Dimension, BTO6, GE Medical Systems, Glattbrugg, Switzerland) with a phased array transducer (M4S, GE Medical Systems, Glattbrugg, Switzerland) at a frequency of 1.9/4.0 MHz (octave harmonics) . A single lead electrocardiogram (lead I) was recorded simultaneously. The lead position was consistent among animals and between repeat echocardiograms. The aortic valve was imaged in 2D grayscale mode in a right parasternal long-axis view of the left ventricular outflow tract (LVOT) and in M-mode in a right parasternal short-axis view, respectively, for subsequent determination of the time of aortic valve closure. The left ventricle (LV) was imaged in 2D mode using a right parasternal four-chamber view optimized to obtain an image of the entire LV at its largest dimensions, including the LV apex. Subsequently, the LV was imaged in three short-axis views, namely at the level of the apex, at the level of the papillary muscles, and at the level of the chordae tendineae. Imaging depth and sector width, respectively, were adjusted to achieve a frame rate > 50 frames/second in 2D imaging mode. Three representative, non-consecutive cardiac cycles were recorded in each view and stored as cine-loops in digital raw data format. Data analyses were performed offline, blinded, and in random order, using a dedicated software package (EchoPAC Software version 6.1.2, GE Medical Systems, Glattbrugg, Switzerland). Three cardiac cycles were analyzed for each imaging plane. For each measured or calculated variable, the average of the three measurements was reported.
The time to aortic valve closure was measured manually (tAVCm). It was defined as the time interval between the peak of the electrocardiographic R wave and the closure point of the aortic valve identified on an M-mode recording of aortic valve motion. In cases, in which the closure point of the aortic valve was not clearly identifiable on the M-mode recording, anatomical M-mode was applied to the recorded two-dimensional echo loop of the LVOT, and cursor placement was adjusted to identify the point of aortic valve closure. The corresponding heart rate (HRAVCm) for each cycle was calculated as 60,000/R-R interval.
The 2D speckle tracking (2DST) analyses were performed using the 2D Strain module of the analysis software (EchoPAC Software version 6.1.2, GE Medical Systems, Glattbrugg, Switzerland). The 2DST variables were measured as follows: (1) The appropriate long-axis or short-axis image was selected and the Q-Analysis module was started. (2) A single heart cycle was selected by moving the left and the right cursor, respectively, to the peak of the R waves of the electrocardiogram. Each cycle was identified by the image number, the time when the image was recorded, and the first and the last frame of the cycle. The R-R interval was measured and the frame rate was recorded. (3) Subsequently, the 2D Strain module was started. The long-axis (subsequently called LAX) grayscale loops of the LV were analyzed using the “4CH” option, and the short-axis grayscale loops of the LV were analyzed using the “SAX-AP” (apical level, subsequently called SAX-AP), “SAX-PM” (papillary muscle level, subsequently called SAX-PM), and “SAX-MV” option (chordal level, subsequently called SAX-CH), respectively. Once the option had been selected, a region of interest (ROI) was determined by tracing the endocardial border of the LV at end-systole. For the long-axis images, tracing started at the septal mitral valve (MV) annulus and ended at the lateral MV annulus; for the short-axis images, tracing was started at mid-septum and proceeded in a clockwise direction. The papillary muscles were not included in the tracings. The ROI width was adjusted, so that the entire myocardial thickness was covered throughout the cardiac cycle. Subsequently the speckle tracking analysis was started. The software algorithm automatically divided the myocardium into 6 segments, performed the speckle tracking analysis, and provided confirmation of adequate tracking for each segment (Additional file 1: Figure S5). The segments were preselected by the software based on regional wall motion analysis standards applied to human patients and were not adjusted for use in the goat. Hence, apart from slight individual variations, the short-axis segments termed “Sept”, “AntSept” and “Ant” generally depict the interventricular septum, whereas “Lat”, “Post”, and “Inf” depict the cranial, lateral, and caudal LV free wall.
The quality of the tracking was visually assessed by the operator during motion playback. If necessary, the line tracing of the endocardium was adjusted and the speckle tracking analysis was repeated until adequate tracking was achieved. If adequate tracking was not possible despite repeated adjustments of the ROI, another cardiac cycle was chosen for analysis. If the software failed to adequately track a segment even after repeated manual tracing of multiple cycles, the loop was excluded from analysis.
Six curve profiles were obtained corresponding to the average of each myocardial segment (Additional file 1: Figure S6). For the long-axis images (4CH), the following variables were reported by the software: Longitudinal strain (εL), longitudinal strain rate (SRL), longitudinal displacement (DL), and transverse displacement (DT). For the short-axis images (SAX-AP, SAX-PM, SAX-CH), the following variables were reported by the software: Circumferential strain (εC), circumferential strain rate (SRC), radial strain (εR), radial strain rate (SRR), radial displacement (DR). The time of aortic valve closure was automatically calculated and displayed by the software (AVCa, Additional file 1: Figure S6).
The measurements for strain, strain rate, and displacement, respectively, were performed on the “Results” screen of the 2D Strain software module. Automated detection of the peak values was verified on the graphical display and corrected as necessary. The mean of the 6 segmental measurements of each variable was calculated to obtain indices of averaged strain, strain rate, and displacement. Strain measurements for each segment included longitudinal, radial, and circumferential peak strain (peak independent of aortic valve closure; termed εL, εR, and εC) as well as peak systolic strain (peak prior to or at the time of aortic valve closure; termed εL-sys, εR-sys, and εC-sys). Where the highest strain value occurred before or at the time of aortic valve closure, peak strain (ε) was identical with peak systolic strain (εsys). Where 2 peaks for strain were present and εsys was followed by a 2nd, higher peak after aortic valve closure, the 2nd peak was considered ε. Where only 1 peak occurring after aortic valve closure was present, εsys was defined as strain at the time of aortic valve closure and the 2nd peak was termed ε. Post-systolic motion was diagnosed where the highest strain value occurred after automatically determined aortic valve closure (AVCa) [7,9,13,14].
Strain rate measurements for each segment included peak systolic strain rate (SRL-sys, SRR-sys, and SRC-sys), peak early-diastolic strain rate (SRL-E, SRR-E, and SRC-E) and peak late-diastolic strain rate (SRL-A, SRR-A, and SRC-A). Displacement measurements for each segment included peak displacement (DL, DT, DR). The R-R interval was measured for each cycle. The time to automatically determined aortic valve closure (tAVCa) was calculated as the time interval between the peak of the electrocardiographic R wave and AVCa displayed on the “Trace” screen. The corresponding heart rate (HRAVCa) for each cycle was calculated as 60,000/R-R interval.
The time interval from the electrocardiographic R wave to longitudinal, circumferential, and radial peak strain (tεL, tεC, tεR) of each segment was measured in each cycle. The synchrony time index (STIε), a measure of myocardial dyssynchrony, was calculated as the difference in tε from the earliest to the latest segment [12,13].
Data analysis and statistics
All statistical and graphical analyses were performed using standard computer software (Microsoft Office Excel 2003, Microsoft Corporation, Redmond, WA; SigmaStat v3.5, SPSS Inc, Chicago, IL; GraphPad Prism v5.00 for Windows, GraphPad Software, San Diego California USA). The agreement between tAVCa and tAVCm was compared by paired t test and Bland-Altman statistics. The corresponding heart rates (HRAVCm and HRAVCa) were compared by paired t test. For segmental 2DST indices, 2-way repeated-measures analysis of variance was used to detect differences between segments and treatment (i.e., awake vs. anesthetized). When the F test indicated significant differences, all pairwise multiple comparisons were performed using the Holm-Sidak post hoc test. The effect of general anesthesia on averaged 2DST indices and on STIε was assessed by paired t tests, reporting the 95% confidence intervals for the difference of means. Validity of the normality assumption was confirmed by assessment of normal probability plots of the residuals. The level of significance was set at P=0.05.
Feasibility and quality of recordings
Echocardiographic analyses by 2DST were feasible in all 22 goats and at both time points (i.e., awake, anesthetized). Accidentally, in one goat the LAX and in another one the SAX-CH images had not been recorded.
The frame rate of the recordings ranged between 70.2 and 103.8 frames per seconds (85.1 ± 2.6 [mean ± SD]). Tracking was considered inaccurate in 1,869 of total 31,512 measured segments (5.93%) based on automated verification and visual assessment of tracking by the operator (LAX, 178 segments; SAX-AP, 488 segments; SAX-PM, 724 segments; SAX-CH, 479 segments).
AVC and HR
There was no significant difference between tAVCa and tAVCm except for measurements obtained at the apical short-axis level during anesthesia and at the chordal short-axis level in awake and anesthetized goats (Table 1). Generally, tAVCa occurred later at the chordal level compared to the papillary muscle and the apical level, respectively. Mean bias was lowest for long-axis analyses and short-axis analyses at the papillary muscle level.
No significant differences were identified between HRAVCa and HRAVCm in any of the views and treatments.
The occurrence of post-systolic motion is summarized in Figure 1. Peak longitudinal strain occurred before or at AVC in the apical segments of most animals, whereas post-systolic motion, with tεL occurring after AVC, was more common in mid and basal segments. In circumferential direction, the occurrence of post-systolic motion differed depending on imaging plane and myocardial segments. Finally, in radial direction, post-systolic motion was seen in the majority of animals in all three imaging planes.
Averaged strain and strain rate
Heart rate, arterial blood pressure, averaged strain and strain rate values in awake and anesthetized goats, respectively, are presented in Table 2. The difference in heart rate between awake and anesthetized goats did not reach statistical significance. Absolute strain values of awake goats were generally lower compared to those in anesthetized goats, but not all differences reached statistical significance. Conversely, absolute systolic strain rate values in awake goats were all significantly higher than those recorded under anesthesia, except SRR-sys at the SAX-AP level. No clear trend was detected for early-diastolic strain rate, whereas absolute values for late-diastolic strain rate were generally lower during anesthesia compared to the awake state.
Absolute values for peak systolic strain (εsys) were significantly lower than those for peak strain (ε), except for circumferential strain in awake goats recorded at the SAX-AP level (Table 3).
Segmental strain, strain rate, and displacement
Segmental analyses are presented in Figures 2 and 3 and in Additional file 1: Table S5. In LAX (Figure 2), there was a gradient in εL from apex to the base, with the largest strain at the apex and the lowest at the base. Furthermore, anesthesia significantly influenced εL independent of segment, with absolute values for εL being higher during anesthesia than in awake goats. Similarly, SRL-sys significantly differed between segments (highest strain rate in apical segments) and between treatments (higher absolute strain rate in awake compared to anesthetized goats, independent of segment). SRL-E significantly differed between segments (highest strain rate in apical segments) but not between treatments. SRL-A significantly differed between segments (without a clear gradient) and between treatments (higher strain rate in awake compared to anesthetized goats, independent of segment). Longitudinal displacement differed significantly between segments (with the largest displacement in basal segments and the smallest in apical segments) and between treatments (with varying effect of anesthesia depending on segment). Transverse displacement differed significantly between segments (with an apparent gradient from lateral basal over apical to septal basal segments) and between treatments (higher displacement in anesthetized compared to awake goats).
In SAX views (Figure 3), εC significantly differed between segments and treatments (with varying effect of anesthesia depending on segment) and with an apparent gradient in anesthetized goats at PM and CH levels. At PM and CH levels, εR significantly differed between segments (without a clear gradient) but not between treatments. DR significantly differed between segments (with highest values in LVFW segments) and between treatments. SRC-sys differed significantly between segments and treatments (with effect of anesthesia depending on segment, except at AP level). SRC-E differed between segments at all levels (without an obvious gradient), with a segmental effect of anesthesia at PM and CH levels. At AP level, SRC-A differed significantly between segments (without an obvious gradient), with varying segmental effects of anesthesia. SRR-sys significantly differed between segments (without an obvious gradient), with a segmental effect of anesthesia at the PM and CH level. SRR-E differed significantly between segments at the PM and CH, but not the AP level. Finally, absolute SRR-A was significantly higher in the awake state at the AP and PM level, and it differed between segments at the PM and the CH level.
Segmental timing of peak strain
Segmental timing of peak strain values (tεL, tεC, tεR) are presented in Figure 4A-G and in Additional file 1: Table S6. Generally, peak strain occurred significantly later in anesthetized goats compared to awake goats. Longitudinal peak strain occurred slightly earlier in apical segments compared to mid and basal segments. Circumferential peak strain generally occurred later in segments “Ant” and “Lat” in the anesthetized goats (although not all SAX-levels reached the level of significance), whereas segmental differences in radial peak strain did not follow an obvious pattern.
Synchrony time index
There were no significant differences between synchrony time indices between anesthetized and awake goats (Table 4).
The present study shows that 2DST is applicable to characterize LV wall motion in goats in an awake, standing position and in sternal recumbency under general anesthesia.
Aortic valve closure and post-systolic motion
The time interval between the electrocardiographic R wave and AVC reflects the duration of the electromechanical systole. Since positioning of the ECG leads might influence QRS conformation and thereby the ability to depict the true onset of electrical events in the ventricle, consistent ECG lead placement is important. Automatically determined tAVCa may be less susceptible to variability compared to manually measured tAVCm using M-mode recordings, because the aortic valve is moving relative to the M-mode cursor line during the cardiac cycle . Furthermore, automated timing of aortic valve closure on each cardiac cycle may in fact be superior to an averaged estimation of tAVCm, as the latter cannot account for variations of heart rate between different cycles . However, relative bias between tAVCa and tAVCm can differ depending on imaging plane. The findings of this study are consistent with previous findings in horses , indicating that tAVCa occurs later at the chordal level compared to the papillary muscle and the apical level, respectively, and that mean bias is small in the long-axis plane (Table 1). However, even the largest relative mean bias (for awake recordings at the chordal level in a short-axis view) was less than 10% and might not be clinically relevant. Therefore, we suggest that automated tAVCa can be readily used for routine 2DST analyses in goats.
The timing of AVC may become important for identification of post-systolic myocardial motion, which is reported to occur in > 30% of myocardial segments in normal human subjects, but may be pathologic if there is a concomitant reduction in systolic strain, especially if the post-systolic thickening is marked [7,9,24]. In previous studies in healthy horses, post-systolic myocardial motion was present in 15 of 16 horses, affecting 31–96% of myocardial segments [13,15,16]. Similarly, the high prevalence of post-systolic motion in this study suggests that this phenomenon might be physiological in goats (Figure 1), with absolute peak strain values (independent of AVC) generally being higher than peak systolic strain values occurring before or at AVC (Table 3). However, no attempts were made to grade the degree of post-systolic motion or to investigate its causes or its clinical relevance. Therefore, it remains uncertain whether the definition of post-systolic motion used in this study would be clinically useful and relevant.
Averaged strain and strain rate
Myocardial function during the heart cycle is determined by preload, contractility, afterload, and heart rate, which are all influenced by autonomic traffic [25,26]. Traditional echocardiographic ejection phase indices of systolic LV function generally do not reflect contractility per se, since they are also influenced by loading, rate, and rhythm. Peak strain correlates well with stroke volume and ejection fraction and is therefore also influenced by changes in loading conditions, contractility, and heart rate [6,8]. In this study, heart rate during anesthesia was slightly lower compared to heart rate in awake goats, although the difference was not statistically significant. However, the relationship between strain and conventional indices of LV function was not assessed and the study design did not allow detailed investigation of the influence of loading conditions and heart rate. Longitudinal and circumferential strain measurements were generally higher in anesthetized compared to awake goats, consistent with an overall increase in systolic LV performance. These findings were in agreement with previously reported findings obtained from the same population of goats, indicating that M-mode-based and 2DE-based LV ejection phase indices and peak systolic LV wall thickness were significantly higher under general anesthesia .
Peak systolic strain rate seems to be more resistant against changes in heart rate and loading conditions and more closely reflects myocardial contractility [5,7,8]. In this study, the lower peak systolic strain rate values in anesthetized compared to awake goats therefore suggest depressed myocardial contractility, which could be explained by negative inotropic drug effects and decreased sympathoadrenergic activity in anesthetized goats. However, overall LV systolic performance appeared to be increased (see above). This can be attributed to the sum of interacting effects of altered loading, contractility, and heart rate caused by anesthetic drugs, fluid therapy, and mechanical ventilation .
Overall, longitudinal and circumferential peak strain and strain rate appear most sensitive to detect differences between awake and anesthetized goats, suggesting that these variables best be used for assessment of LV systolic function in goats.
Strain rate measurements can also be applied for assessment of LV diastolic function. Early-diastolic strain rate (SRE) is determined by active LV relaxation, LV compliance, and filling pressures . Late-diastolic strain rate (SRA) is related to active atrial contraction (determined by atrial contractility, preload, and afterload) and LV compliance at end-diastole . In this study, SRE was not significantly different in anesthetized compared to awake goats, suggesting that LV relaxation and filling was not markedly altered by anesthesia or that 2DST is not sensitive enough to detect slight alterations. Conversely, SRA was higher in awake goats in longitudinal as well as radial (apical and papillary muscle level) and circumferential (apical level) direction, consistent with a slight depression of active atrial contraction in the anesthetized state. This was in agreement with previous findings obtained on the same population of goats, indicating that active left-atrial contraction was slightly lower in anesthetized goats .
Segmental strain, strain rate and displacement
Segmental 2DST indices represent regional contractile function and are used for detection of abnormal regional wall motion patterns and ventricular dyssynchrony and in cardiac resynchronization therapy in humans .
In this study, longitudinal strain, systolic strain rate, and early-diastolic strain rate decreased gradually from apex to base (Figure 2). This suggests that the apical myocardial segments of the LV deform more dynamically compared to the mid-wall and basal segments, which is in agreement with findings in previous experimental studies and in human and horse studies [14,16,28-30]. The gradual increase in longitudinal displacement from apex to base (Figure 2E) is consistent with the fact that the mitral annulus is pulled down toward the apex during the ejection phase [16,28].
In SAX imaging planes, segmental differences and effects of anesthesia were generally more pronounced in the circumferential compared to the radial direction. The higher radial strain, strain rate, and displacement in the “Lat”, “Post”, “Inf”, and “Sept” segments in the SAX-PM and SAX-CH planes indicate more pronounced systolic thickening and motion of the LV free wall compared to the septum. These findings also correspond to a previous study in horses .
Overall, effects of anesthesia corresponded to those described above for averaged strain and strain rate values, suggesting improved overall systolic LV function in face of depressed myocardial contractility. Effects of anesthesia on radial displacement were in agreement with previously reported 2DE and M-mode indices of LV systolic function in the same population of goats .
The timing of segmental peaks strain describes the propagation of mechanical activation. In agreement with previous studies [28,31], maximum longitudinal LV deformation in this study occurred significantly earlier at the apex compared to mid ventricular and basal segments (Figure 4A), consistent with the fact that myocardial depolarization begins near the apical septum and spreads toward the base .
In SAX imaging planes, timing of circumferential peak strain values showed more segmental variation compared to radial strain values (Figure 4B-G), possibly related to the helical arrangement of the myofibers. Cranial and lateral segments (“Ant”, “Lat”) generally showed a delayed peak activation compared to other segments. This has also been observed in horses .
During anesthesia, peak strain generally was delayed, suggesting that ejection time was prolonged. This can be attributed to the anesthesia-induced decrease in heart rate (Table 1), combined with altered contractility and loading conditions . No attempts were made to correct segmental peak timing for differences in heart rate.
There is no single synchrony index that is preferred for assessment of ventricular dyssynchrony in people. Studies suggest that combining dyssynchrony data from different methods may be of additive value . In this study, we chose to investigate a synchrony time index that can easily be calculated based on the maximum differences in timing of peak strain [12,27]. Based on the STIs, ventricular synchrony appeared higher in radial compared to longitudinal and circumferential direction and was not influenced by general anesthesia. However, other studies in horses [13,14] and humans [9,34] indicated that reliability of the STI was insufficient for clinical use. Therefore, the effect of anesthesia on ventricular synchrony in goats may not be conclusively assessed based on the present data and further studies are needed to investigate the clinical relevance and the best diagnostic approach to LV dyssynchrony in goats.
Acquisition of 2DE images for assessment of longitudinal motion of the LV in large animals such as adult goats is limited by the fact that apical long-axis views cannot be obtained because of anatomical constraints. While the 2DST software algorithm used in this study is able to correctly track longitudinal myocardial motion independent of image orientation (personal communication, manufacturer’s application specialist), spatial resolution is slightly reduced when the longitudinal axis of the heart is perpendicular to the ultrasound beam. Therefore it is important to standardize imaging planes when comparing studies between individuals or over time and when comparing individual measurements to normal reference intervals.
Since only female goats were included in this study, the influence of sex on 2DST variables could not be assessed. Regression analyses were not able to detect any significant relationship of 2DST variables to age and body weight, respectively (data not shown), but the study population was relatively homogenous and the range of available ages and body weights was narrow.
A true limitation of this study is the lack consideration of alterations in blood pressures between the awake and the anesthetized state, since loading conditions might significantly influence strain variables. However, the current study setting did not allow invasive measurements of arterial blood pressures in awake goats.
Also, this study was not designed to comprehensively assess the repeatability and reproducibility of recordings and measurements. Therefore, while significant differences between segments and treatments, respectively, can be established on a population level, the diagnostic validity for detection of subtle changes in strain and strain rate in individual goats in a clinical setting is unknown to date.
In conclusion, 2SDT is a valid method for non-invasive characterization of LV wall motion in awake and anesthetized goats. In conjunction with conventional 2D, M-mode, and Doppler echocardiography, it may add to the understanding of LV mechanical function and may aid in the diagnosis of global and segmental LV systolic dysfunction in goats. However, effects of anesthesia and species-specific characteristics should be considered when goats are used as animal models for human disease. Furthermore, future studies are required to assess the reliability of 2DST measurements when used repeatedly in individual animals over time.
Two-dimensional speckle tracking
Four chamber view
Interventricular septal segment in short axis view
Interventricular septal segment in short axis view
Aortic valve closure (automatically determined)
Aortic valve closure (manually determined)
Circumferential peak systolic strain
Peak systolic strain
Longitudinal peak systolic strain
Radial peak systolic strain
Instantaneous heart rate (for AVCa)
Instantaneous heart rate (for AVCm)
Caudal LV free wall segment
Cranial LV free wall segment
Left ventricle or left ventricular
Left ventricular free wall
Left ventricular outflow tract
Mean arterial blood pressure
Lateral LV free wall segment
Region of interest
Short axis at apical level
Short axis at chordal level
Short axis at papillary muscle level
Interventricular septal segment in short axis view
Circumferential strain rate
Circumferential peak late-diastolic strain rate
Circumferential peak early-diastolic strain rate
Circumferential peak systolic strain rate
Longitudinal strain rate
Longitudinal peak late-diastolic strain rate
Longitudinal peak early-diastolic strain rate
Longitudinal peak systolic strain rate
Radial strain rate
Radial peak late-diastolic strain rate
Radial peak early-diastolic strain rate
Radial peak systolic strain rate
Synchrony time index
Time interval from the electrocardiographic R wave to circumferential peak strain
Time interval from the electrocardiographic R wave to longitudinal peak strain
Time interval from the electrocardiographic R wave to radial peak strain
Time of aortic valve closure (manually determined)
Time of aortic valve closure (automatically determined)
Tissue Doppler imaging
Boon JA. Evaluation of size, function, and hemodynamics. In: Veterinary Echocardiography. 2nd edition. Oxford: Wiley-Blackwell; 2010. p. 153–266.
Otto CM. Left and right ventricular systolic function. In: Textbook of Clinical Echocardiography. 3rd ed. Philadelphia: Elsevier Saunders; 2004. p. 131–65.
Teske AJ, Boeck BWD, Melman PG, Sieswerda GT, Doevendans PA, Cramer MJ. Echocardiographic quantification of myocardial function using tissue deformation imaging, a guide to image acquisition and analysis using tissue Doppler and speckle tracking. Cardiovasc Ultrasound. 2007;5:27.
Dandel M, Hetzer R. Echocardiographic strain and strain rate imaging—clinical applications. Int J Cardiol. 2009;132:11–24.
Edvardsen T, Gerber BL, Garot J, Bluemke DA, Lima JAC, Smiseth OA. Quantitative assessment of intrinsic regional myocardial deformation by Doppler strain rate echocardiography in humans: validation against three-dimensional tagged magnetic resonance imaging. Circulation. 2002;206:50–6.
Sutherland GR, Salvo GD, 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.
Marwick TH. Measurement of strain and strain rate by echocardiography: ready for prime time? J Am Coll Cardiol. 2006;47:1313–27.
Strain rate imaging: Cardiac deformation imaging by ultrasound/echocardiography–Tissue Doppler and speckle tracking. [http://folk.ntnu.no/stoylen/strainrate/] Last accessed: March 2015.
Abraham TP, Dimaano VL, Liang H-Y. Role of tissue Doppler and strain echocardiography in current clinical practice. Circulation. 2007;116:2597–609.
Lindqvist P, Borgström E, Gustafsson U, Mörner S, Henein MY. Asynchronous normal regional left ventricular function assessed by speckle tracking echocardiography Appearances can be deceptive. Int J Cardiol. 2009;134:195–200.
Cho GY, Marwick TH, Kim HS, Kim MK, Kyung-Soon H, Oh DJ. Global 2-dimensional strain as a new prognosticator in patients with heart failure. J Am Coll Cardiol. 2009;54:618–24.
Chetboul V, Serres F, Gouni V, Tissier R, Pouchelon JL. Radial strain and strain rate by two-dimensional speckle tracking echocardiography and the tissue velocity based technique in the dog. J Vet Cardiol. 2007;9:69–81.
Schwarzwald CC, Schober KE, Berli AS, Bonagura JD. Left ventricular radial and circumferential wall motion analysis in horses using strain, strain rate, and displacement by 2D speckle tracking. J Vet Intern Med. 2009;23:890–900.
Schefer KD, Bitschnau C, Weishaupt MA, Schwarzwald CC. Quantitative analysis of stress echocardiograms in healthy horses with 2-dimensional (2D) echocardiography, anatomical M-mode, tissue Doppler imaging, and 2D speckle tracking. J Vet Intern Med. 2010;24:918–31.
Decloedt A, Verheyen T, Sys S, De Clercq D, Van Loon G. Two-dimensional speckle tracking for quantification of left ventricular circumferential and radial wall motion in horses. Equine Vet J. 2013;45:47–55.
Decloedt A, Verheyen T, Sys S, De Clercq D, Van Loon G. Quantification of left ventricular longitudinal strain, strain rate, velocity, and displacement in healthy horses by 2-dimensional speckle tracking. J Vet Intern Med. 2011;25:330–8.
Smith MC, Sherman DM. Cardiovascular System. In: Goat Medicine. Ames. Iowa: Wiley-Blackwell; 2009. p. 319–38.
Greiser M, Neuberger HR, Harks E, El-Armouche A, Boknik P, De Haan S, et al. Distinct contractile and molecular differences between two goat models of atrial dysfunction: AV block-induced atrial dilatation and atrial fibrillation. J Mol Cell Cardiol. 2009;46:385–94.
Remes J, Van Brakel TJ, Bolotin G, Garber C, De Jong MM, Van Der Veen FH, et al. Persistent atrial fibrillation in a goat model of chronic left atrial overload. J Thorac Cardiovasc Surg. 2008;136:1005–11.
Van Brakel TJ, Hermans JJR, Accord RE, Schotten U, Smits JFM, Allessie MA, et al. Effects of intrapericardial sotalol and flecainide on transmural atrial electrophysiology and atrial fibrillation. J Cardiovasc Electrophysiol. 2009;20:207–15.
Hallowell GD, Potter TJ, Bowen IM. Reliability of quantitative echocardiography in adult sheep and goats. BMC Vet Res. 2012;8:181.
Steininger K, Berli AS, Jud R, Schwarzwald CC. Echocardiography in Saanen-goats: normal findings, reference intervals in awake goats, and the effect of general anesthesia. Schweiz Arch Tierheilkd. 2011;153:553–64.
Aase SA, Støylen A, Ingul CB, Frigstad S, Torp H. Automatic timing of aortic valve closure in apical tissue Doppler images. Ultrasound Med Biol. 2006;32:19–27.
Voigt JU, Lindenmeier G, Exner B, Regenfus M, Werner D, Reulbach U, et al. Incidence and characteristics of segmental postsystolic longitudinal shortening in normal, acutely ischemic, and scarred myocardium. J Am Soc Echocardiogr. 2003;16:415–23.
Opie LH. Heart Physiology: From Cell to Circulation. Philadelphia: Lippincott Williams & Wilkins; 2004.
Schwarzwald CC, Bonagura JD, Muir WW. The cardiovascular system. In: Equine anesthesia-Monitoring and Emergency Therapy. 2nd ed. St. Louis: Elsevier S; 2009. p. 37–100.
Delgado V, Ypenburg C, Van Bommel RJ, Tops LF, Mollema SA, Marsan NA, et al. Assessment of left ventricular dyssynchrony by speckle tracking strain imaging comparison between longitudinal, circumferential, and radial strain in cardiac resynchronization therapy. J Am Coll Cardiol. 2008;51:1944–52.
Sengupta PP, Krishnamoorthy VK, Korinek J, Narula J, Vannan MA, Lester SJ, et al. Left ventricular form and function revisited: applied translational science to cardiovascular ultrasound imaging. J Am Soc Echocardiogr. 2007;20:539–51.
Sengupta PP, Korinek J, Belohlavek M, Narula J, Vannan MA, Jahangir A, et al. Left ventricular structure and function: basic science in cardiac imaging. J Am Coll Cardiol. 2006;48:1988–2001.
Saito K, Okura H, Watanabe N, Hayashida A, Obase K, Imai K, et al. Comprehensive evaluation of left ventricular strain using speckle tracking echocardiography in normal adults: comparison of three-dimensional and two-dimensional approaches. J Am Soc Echocardiogr. 2009;22:1025–30.
Sengupta PP, Khandheria BK, Korinek J, Wang J, Jahangir A, Seward JB, et al. Apex-to-Base dispersion in regional timing of left ventricular shortening and lengthening. J Am Coll Cardiol. 2006;47:163–72.
Hamlin RL, Smith CR. Categorization of common domestic mammals based upon their ventricular activation process. Ann N Y Acad Sci. 1965;127:195–203.
Gorcsan 3rd J, Abraham T, Agler DA, Bax JJ, Derumeaux G, Grimm RA, et al. 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. 2008;21:191–213.
Voigt JU, Exner B, Schmiedehausen K, Huchzermeyer C, Reulbach U, Nixdorff U, et al. Strain-rate imaging during dobutamine stress echocardiography provides objective evidence of inducible ischemia. Circulation. 2003;107:2120–6.
We would like to thank Prof. U. Braun for putting the goats at our disposal and the veterinarians and grooms of the Farm Animal Department for animal care and animal handling.
Parts of this study have been presented as a poster presentation at the 27th annual ACVIM Forum, June 3–6, 2009, Montreal, Quebec, Canada.
The authors declare that they have no competing interests.
AJB was responsible for echocardiographic measurements, data analysis and statistics, and manuscript writing. RJ was responsible for conducting anesthesia and manuscript writing. KS was assisting with data collection and echocardiographic examinations. CCS was responsible for planning and supervising the study, for conducting all echocardiographic examinations, and assisting with data analysis, statistics, and manuscript writing. All authors read and approved the final manuscript.
Two-dimensional Speckle Tracking Analysis Tool. Figure S6. Two-dimensional Speckle Tracking Trace Screens. Table S5. Segmental 2DST analyses of left ventricular long-axis and short-axis recordings under awake and anesthetized conditions. Table S6. Segmental timing of peak strain in left ventricular long-axis and short-axis recordings under awake and anesthetized conditions.
About this article
Cite this article
Berli, AS.J., Jud Schefer, R., Steininger, K. et al. The use of strain, strain rate, and displacement by 2D speckle tracking for assessment of systolic left ventricular function in goats: applicability and influence of general anesthesia. Cardiovasc Ultrasound 13, 11 (2015). https://doi.org/10.1186/s12947-015-0005-8