Echocardiographic quantification of myocardial function using tissue deformation imaging, a guide to image acquisition and analysis using tissue Doppler and speckle tracking
Cardiovascular Ultrasound volume 5, Article number: 27 (2007)
Recent developments in the field of echocardiography have allowed the cardiologist to objectively quantify regional and global myocardial function. Regional deformation (strain) and deformation rate (strain-rate) can be calculated non-invasively in both the left and right ventricle, providing information on regional (dys-)function in a variety of clinical settings. Although this promising novel technique is increasingly applied in clinical and preclinical research, knowledge about the principles, limitations and technical issues of this technique is mandatory for reliable results and for implementation both in the clinical as well as the scientific field.
In this article, we aim to explain the fundamental concepts and potential clinical applicability of strain and strain-rate for both tissue Doppler imaging (TDI) derived and speckle tracking (2D-strain) derived deformation imaging. In addition, a step-by-step approach to image acquisition and post processing is proposed. Finally, clinical examples of deformation imaging in hypertrophic cardiomyopathy (HCM), cardiac resynchronization therapy (CRT) and arrhythmogenic right ventricular dysplasia/cardiomyopathy (ARVD/C) are presented.
There is an increasing need for diagnostic modalities able to objectively quantify myocardial function. Quantification of regional myocardial function with ultrasound is challenging on its own. Visual assessment of wall motion and thickening requires extensive training and remains highly subjective. In addition to this visual assessment on 2-dimensional B-mode images and new 3-dimensional recordings, "classic" M-mode echocardiography can be used to evaluate myocardial thickening. Unfortunately, M-mode can only be used in a limited number of segments since it has to be perpendicular to the investigated myocardial segment. New techniques such as anatomical M-mode have partly overcome the problems with insonation angle, at the expense of a reduced temporal resolution. Moreover, this method only provides uni-dimensional semi-quantitative assessment of myocardial thickening/thinning without information on longitudinal or circumferential function. Tissue velocity imaging along the long axis of the ventricle is a commonly used and well validated technique to quantify longitudinal ventricular function. Tissue velocities at the ventricular base represent the integral of myocardial shortening velocity from base to apex and therefore provide information on global, rather than on regional ventricular function.
Tissue deformation imaging is a recently introduced technique which enables the objective assessment of regional myocardial deformation assessed by ultrasound based strain and strain rate using myocardial Doppler data or B-mode images. This is a promising new technique to quantify regional left and right ventricular function and appears of added value in unmasking or unraveling cardiac pathology (table 1). Knowledge of its principles and limitations is mandatory for proper application and reliable interpretation of the results both in the clinical as well as the scientific setting. This review therefore aims to explain the basic concepts and potential clinical applicability of strain (ε) and strain-rate (SR) for both tissue Doppler derived and speckle tracking derived deformation imaging. In addition, a step-by-step approach of image acquisition and post processing will be discussed.
Following electro-mechanical activation, the myocardium deforms during systole due to sarcomere shortening. This active deformation causes a reduction in intracavitary size, resulting in the ejection of blood from the ventricle. In diastole the original ventricular geometry is restored due to active relaxation and passive filling following atrial contraction. Since myocardial tissue is virtually incompressible, the volume of the ventricular wall remains the same during the cardiac cycle and, thus, deforms in three dimensions. During systole there is three-dimensional deformation which can be expressed in three ventricular coordinates: a longitudinal shortening and a circumferential shortening and a radial thickening.
Myocardial deformation is expressed as a one dimensional parameter: strain (ε). It defines the total deformation during the cardiac cycle relative to the initial length at the onset of the cardiac cycle, and is expressed in percentages (figure 1). This implies that longitudinal and circumferential (systolic) shortening result in a negative strain and radial (systolic) thickening in a positive strain. When the initial length of the investigated myocardial segment is known during the cardiac cycle, the relative length change (strain) can be calculated throughout the entire cardiac cycle. The local end-systolic strain value reflects the regional ejection fraction and the global LV end-systolic strain reflects the LV ejection fraction  (Figure 2). When visualized in a graph (figure 3) the different phases of the cardiac cycle can usually be identified: during systole the strain values become more negative (S-wave) with the negative peak at the aortic valve closure, representing the maximal longitudinal myocardial shortening during contraction (or peak systolic strain). In diastole the strain values return towards zero (towards the original length of the analyzed myocardial segment at the onset of the cardiac cycle) in three phases: (1) the early, or rapid filling phase (E-wave) followed by (2) a plateau phase, or diastasis and finally (3) atrial filling (A-wave).
The speed at which the myocardial deformation occurs is the strain-rate (SR) and is expressed in s-1. SR depicts the change in strain over a period of time. Thus, when the myocardium shortens there is a negative SR and, the steeper the slope of the strain-curve, the higher the SR-values. The peak systolic strain rate correlates well to loading independent indices of contractility and hence provides valuable information on regional contractile function [3, 4]. Comparable to strain, when plotted in a graph (figure 4) the different cardiac phases can be recognized: in systole a negative deflection indicating myocardial shortening with a peak systolic SR at the steepest part of the strain-curve, and in diastole two positive deflections: the E- and A-wave and diastasis when no change in strain occurs.
To optimally assess regional myocardial function, both strain and SR need to be calculated since they provide complementary information: end systolic strain estimates ejection fraction and peak systolic SR is a measure of contractility. By analogy, when driving a car both the total distance of the journey as well as the speed of the car during the journey provide valuable information.
Tissue deformation imaging
Currently there are two different methods to calculate myocardial deformation using ultrasound: tissue Doppler derived strain and two-dimensional strain derived from B-mode images (speckle tracking). These two techniques will be briefly explained in the following section. Notably, differences in software packages may exist between manufacturers . The examples in this review are based on the commercially available software packages from GE Vingmed (GE Vingmed ultrasound, Horten, Norway).
Tissue Doppler strain
By color coded tissue Doppler images (TDI), the velocity relative to the transducer is recorded and can be calculated in each pixel (figure 5). In the longitudinal recordings, the gradual increase in velocity from the nearly stationary apex to a maximal velocity at the annular plane is called the velocity gradient and is generated by longitudinal shortening of the myocardium. Within a defined region of interest (ROI) the difference in velocities is used to calculate the SR using the velocity of the distal (v1) and the proximal (v2) point within the ROI and the distance between these two points (see formula, figure 4). Strain is calculated from the SR values by temporal integration.
From this formula it is also clear that if a myocardial segment is not deforming during systole, for example due to an infarction, there is no difference in velocities and thus a SR of zero, although the segment could show a displacement and velocity due to tethering of a functional neighboring segment. This translated velocity on the non-deforming segment could give the impression of normal function when assessing it by velocity imaging only, but this is unmasked when calculating the regional SR and ε.
TDI is one-dimensional, thus only deformation along the ultrasound beam can be derived from velocity. Therefore, only the longitudinal SR/ε in the apical views and radial SR/ε in the parasternal views can be calculated.
Several studies have established the robustness and potential useful role of this novel echocardiographic technique.
The accuracy of the measurement of velocities has recently been investigated using tissue mimicking gelatin blocks in two commercially available ultrasound systems . This in-vitro study reported a very high degree of accuracy for velocity measurements, although there was a small overestimation of ε and SR. In-vivo invasive studies, using sonomicrometry as a reference method, have confirmed the validity to calculate strain and strain-rate from tissue velocity over a wide range of strains in animal models [6, 7]. Nowadays, the reference method for non-invasive validation in humans is three-dimensional tagged magnetic resonance imaging, which relies on deformation of tag lines to assess strain in 3 dimensions. Doppler derived strain measurements in humans, both healthy controls and patients with acute myocardial infarction, under different conditions (at rest and during Dobutamine infusion) have proven to be accurate when comparing ultrasound with MRI .
Thus, tissue Doppler imaging is a well validated technique to calculate strain and SR over a wide range of conditions, indicating the feasibility for the usefulness of this technique in the clinical setting.
When evaluating regional myocardial function with tissue Doppler derived SR/ε, knowledge of the limitations of this technique is essential to ensure appropriate acquisition as well as correct post-processing since artifacts often mimic pathology. The most important and common artifacts will be briefly discussed in this section.
Since all calculations are derived from the measured velocities, it is essential to work with a data set which provides the true velocities of the myocardium. Reverberations, side lobes and drop out all interfere with obtaining the true velocities and disturb the data used to calculate SR and strain. Reverberations and side lobes are stationary artifacts of extracardiac structures projecting through the myocardium. The measured velocity within such an artifact is close to zero, resulting in 1) random noise within the artifact 2) an overestimation of negative SR/ε below the artifact, 3) a reversal directly above the artifact and 4) an underestimation above the reversal region (Figure 6). If, for example, the true v1 is 4 cm/s and v2 is 3 cm/s, but due to a stationary artifact the measured v2 is 0 cm/s. This means that the calculated Δv is 4 cm/s instead of the true 1 cm/s, overestimating the SR (see formula) and strain. Drop-out interferes in a similar manner in the calculation of SR, since no velocities are measured, within the drop out zone SR is close to zero, giving the impression of hypo- or akinesia.
Deviation of the insonation angle is another important limitation of the technique. If the ultrasound beam deviates from the direction the investigated myocardial segment moves, the measured velocity (vm) is lower than the actual myocardial velocity. This underestimation is dependent on the angle between the myocardial wall and the ultrasound beam, this implies that vm reduces if the angle increases. Angle correction however, is not an option. As explained earlier, the myocardium deforms in three dimensions, meaning that when measuring longitudinal deformation, one has to keep in mind the simultaneous radial deformation. When the ultrasound beam is parallel to the myocardial wall, the vm is the actual velocity and the vm of the transverse deformation is zero, since the movement is perpendicular to the ultrasound beam. When the angle changes, the contribution of the transverse velocity increases in the vm. For the resulting SR/ε this results in a combination of both a negative (longitudinal) and a positive (transverse) strain, severely disturbing the derived SR/ε patterns. For acceptable calculations, an angle deviation below the 20 to 15 degrees is mandatory.
Image acquisition and post-processing
As pointed out, TDI is a technique which has some limitations inherent to the technique, which could potentially disturb the derived parameters. Optimal image acquisition and knowledge of post-processing techniques are essential to minimize artifacts and, respectively, deal with them in post processing.
The ultrasound investigation starts with an optimal ECG signal with a clear definition of the QRS-complex and P-wave ensuring a consistent ECG-triggering. Pulsed or continuous wave measurements of the valves are needed to incorporate event timing in the off-line measurements. An optimal quality in the 2D gray scale image is mandatory for a high quality tissue Doppler recording. In the following section we present a step-by-step approach for image acquisition.
1) The myocardial wall needs to be optimally visualized with a clear delineation of myocardial tissue and extracardiac structures (e.g. suppression of the blood pool), a minimum of extracardial artifacts should be achieved and the angle of interrogation limited (deviation below 15–20 degrees). The latter is achieved by orienting the transducer parallel or perpendicular to the investigated wall.
2) TDI settings need to be adjusted for optimal post-processing. The velocity scale needs to be adjusted to avoid aliasing. Although a reduction in the velocity scale will result in lower frame rates, the lowest possible velocity scale will result in the highest spatial and temporal resolution combination. Further optimization of the temporal and spatial resolution is achieved by narrowing of the image sector until it is slightly wider than the investigated wall or wall segment. Wall by wall acquisition is essential to achieve adequate frame rates, for reliable SR calculations the optimal sampling rate is >180 frames/sec.
3) Three complete cardiac cycles are recorded digitally during breath hold (to minimize image translation) in the presence of a regular rhythm. This will allow filtering of average noise during post processing without losing resolution.
Longitudinal ε/SR can be derived from the apical three (posterior and anteroseptal wall), two (inferior and anterior wall) and four chamber views (lateral, septal and RV free wall). Radial strain can be derived from the parasternal recordings of the posterior wall. The measurement of septal wall radial strain often produces artifacts due to the morphological and functionally bilayered architecture (LV and RV fibers) of the interventricular septal wall . Circumferential strain can be calculated from the short axis recording (lateral and inferoseptal wall).
Knowledge and experience in the post-processing techniques are important to interpret the ε and SR graphs and to distinguish pathological deformation from artifacts. All commercial software programs are equipped with numerous settings to offer the possibility to optimize the calculations to derive SR from velocity. These settings are available to reduce the 'random-noise' in the SR and ε graphs. With the reduction of noise, the SR and ε graphs are smoothed with the consequence that potentially important temporal or regional information is lost. Therefore, a maximal signal to noise ratio should be achieved with a maximal spatial resolution. Averaging more than one cardiac cycle (reduction of random noise and the beat-to-beat variation), temporal (reducing the effective sampling rate) and spatial averaging and increasing the offset length (changing the Δx in figure 1) are alternative options to reduce noise, although, with careful tracking this is often unnecessary using default settings. Another important setting is the correction for drift in the strain graphs. Assuming that the length of the myocardium returns to its original length at the end of the cycle, the strain value should return to zero. Time integrating the strain-rate can result in drifting of the strain curve. Drift compensation could introduce an error in the strain estimates and should be taken into account when interpreting the strain graphs.
All derived parameters are calculated within a specific region of interest (ROI), therefore the ROI needs to follow myocardial motion during the cardiac cycle. A brief explanation on tracking is described in the following segment.
1) The ROI should be defined; this implies adjustment of both the offset length and ROI size. These parameters can vary and depend on the investigated parameter (radial vs. longitudinal strain) or the clinical question (since small defects appear normal in a large ROI or offset length).
2) The ROI should be placed in the myocardium, typically, three ROI locations are defined in each wall (basal-, mid- and apical segment).
3) The ROI should be tracked throughout the cardiac cycle to follow the myocardial movement. Special attention should be made to avoid stationary artifacts (e.g. reverberations) and the blood pool (see Additional File 1).
Additional file 1: ROI tracking in tissue Doppler imaging. The ROI is tracked (bottom left) as described in the text in the basal segment of the RV. Note the changes in the strain-graph during the tracking. After careful tracking, the ROI follows the myocardial motion. (AVI 10 MB)
4) Finally the quality of the tracking should be checked. Sub-optimal tracking of the ROI will result in drift in the strain-graphs and in a noisy SR-graph. The ROI position can be corrected after the identification of the noisy frames (to eliminate frame-to-frame artifact or blood-pool noise) which are best identified in the SR graph.
Numerous measurements and calculations can be performed after the tracking of the ROI. These can be divided into parameters concerning the timing and/or the magnitude of strain and SR. The most relevant measurements are presented in figure 7. For shortening strains, the peak systolic strain value is defined as the lowest value between aortic valve opening and closure. Post-systolic strain is defined as myocardial deformation (shortening for longitudinal and lengthening for radial function) after the aortic valve closure. From these values, the post-systolic index (PSI) can be calculated (figure 7). For the SR trace, comparable parameters can be obtained as in the velocity trace; peak systolic, peak early- and late diastolic SR and time to peak systolic SR, although the last measurement can be unreliable, due to the noisy aspect of the SR curves.
Speckle tracking/two-dimensional strain
2-dimensional strain imaging is a new technique which uses standard B-mode images for speckle tracking analysis, in which the speckled pattern (acoustic backscatter generated by the reflected ultrasound beam) is followed frame by frame. This speckle pattern is unique for each myocardial region and it is relatively stable throughout the cardiac cycle. The displacement of this speckled pattern is considered to follow myocardial movement and a change between speckles represents myocardial deformation (Figure 8). When tracking a defined region of speckles, a software algorithm follows the change in geometric position of this region, frame by frame, and extracts the displacement, velocity, strain and strain rate of a defined myocardial segment. In contrast to TDI derived parameters, speckle tracking is an angle independent technique as the movement of speckles can be followed in any direction. For the apical views this implies that not only longitudinal, but also transverse parameters can be calculated, which is not possible in TDI recordings. In short axis images both circumferential and radial parameters can be calculated for all myocardial segments. In addition to the circumferential deformation parameters, ventricular rotation and twist (additional parameters for LV function) can also be calculated using this technique [10, 11]. Since the deformation parameters can be calculated in two dimensions (while TDI derived parameters are one-dimensional), this technique is often referred to as two-dimensional strain echocardiography (2DSE).
Validation of early versions of 2DSE software in a tissue mimicking gelatin block revealed a good correlation compared to sonomicrometry, although values (both strain and strain-rate) were overestimated in the lower range of values . These findings were also found in an animal model (sonomicrometry in pigs with myocardial infarction), where 2DSE was found to slightly overestimate values compared to the reference test. A weaker correlation was found in the lower range of deformation values . Tracking is affected by dyskinetic segments (e.g. ischemia), causing the fiber orientation to change, thus affecting acoustic properties and speckle integrity. Using a new (two-stage) tracking algorithm, a better correlation was found between the 2DSE values and those obtained using sonomicrometry over the entire broad range of tested clinical relevant values . For radial and circumferential parameters, a good correlation was found using sonomicrometry in dogs .
The first study to determine the feasibility and accuracy of 2DSE in patients with myocardial infarction reported that of all segments, 98% in the control population and 80% in the patient group could be analyzed, findings which were reproduced in later studies . Importantly, 2DSE data are highly reproducible and analysis is affected by only small intraobserver (mean 4.4 (SD 1.6)%) and interobserver variability (7.3 (SD 2.5)%) . Values were significantly reduced in the infarcted segments, implying a reduction in systolic deformation. There was no significant difference between TDI and 2DSE parameters . In a direct comparison of TDI and 2DSE with MRI-tagging as a reference, comparable values were found for 2DSE and TDI in both normal and dysfunctional segments. For radial measurements, 2DSE was more reliable than TDI .
2DSE is generally more smoothed and, therefore, less influenced by artifacts than TDI derived strain (Figure 9/Additional File 2), but also less sensitive to detect small regions of pathology. Some other limitations of this technique exist. Despite noise reduction by increased smoothing and despite an anti-clutter filter, drop-out is still problematic and stationary reverberations are sometimes still tracked or interfere with the frame-by-frame tracking, resulting in drift or incorrect calculation of strain and strain-rate. In addition, a clear delineation of the endocardial border is also important for a reliable radial and transverse tracking. Since spatial resolution is lower in the transverse direction, tracking perpendicular to the ultrasound beam becomes less robust.
At present, the optimal frame rate for speckle tracking seems to be 50–70 frames per second (FPS), which is lower compared to TDI (>180 FPS): this could result in undersampling, especially in patients with tachycardia. Rapid events during the cardiac cycle (e.g. isovolumic phases) may disappear all together, and peak SR and velocity values may be reduced due to under sampling, especially in isovolumic phases and in early diastole. Higher frame rates could reduce the undersampling problem, although this will result in a reduction of spatial resolution and, consequently, less optimal ROI tracking . Low frame rates increase the spatial resolution, but introduce a new problem. Speckle-tracking software uses a frame-by-frame approach to follow the myocardial movement and searches each consecutive frame for a speckle pattern closely resembling and in close proximity to the reference frame. With a too low frame rate the speckle pattern could be outside the search area, again resulting in poor tracking.
Another limitation of 2SDE is the fact that calculated parameters are averaged over the myocardial segment when using the "results page" of the software program (see figure 10). For most pathologic conditions this is not an important limitation, although in conditions with small regions of myocardial dysfunction, such as early stages of hypertrophic cardiomyopathy or arrhythmogenic right ventricular dysplasia, the averaging could result in normal deformation parameters due to normal deformation properties of the adjacent myocardium.
Optimal tracking of the ROI is not only dependent on optimal image quality, but also on the implemented tracking algorithm. Software programs designed for speckle tracking are relatively new and are subjected to periodical improvements. Implementation of a new two-stage tracking algorithm (correcting initial 'rough', block matching estimates using optical flow techniques ) has indeed shown to improve tracking quality . However, different tracking algorithms potentially produce different results, therefore it should be kept in mind that a periodical update of the software package conceivably influences reference values.
Image acquisition and post-processing
The same basic principles for image acquisition as mentioned for TDI recordings apply for 2DSE: digitalized recordings of a minimum of one cardiac cycle are made during breath hold (to minimize through plane motion) with a stable ECG recording. Since the angle of interrogation does not influence strain and strain-rate calculations, the transducer can be placed off-axis to obtain the optimal gray scale image. However, misalignment in the third dimension (foreshortening in apical views and oblique transsections in parasternal views) cannot be corrected in post-processing and should be correct. We would like to state that optimal frame rates for 2DSE are still under discussion, but in our own experience: a frame rate of 50–80 FPS seems to result in the optimal tracking quality in the left ventricle long and short axis. For single wall recordings, a frame rate between 70–110 FPS will result in proper tracking .
Visualization of the myocardial wall needs to be optimal, with a clear delineation of myocardial tissue and extracardiac structures and avoidance of drop-out since this will result in unacceptable ROI tracking and drift. The optimal balance between temporal and spatial resolution is achieved when adjusting the image sector slightly wider than the interrogated wall. For small-angle acquisitions this has the advantage that frame-rate can be reduced into the 50–110 FPS range while greatly enhancing the lateral resolution and overall quality of the image. Moreover, near field clutter can be reduced with the implementation of dual-focus imaging when artifacts are prominent. For parasternal recordings (for radial and circumferential function on the short axis) small angle recordings are not possible, all myocardial segments need to be in the image field. For any recording, it should be kept in mind that all myocardial segments need to be visualized optimal, since poor tracking of any segment often affects tracking of the adjacent segments.
Comparable to TDI calculations, careful post-processing is mandatory to obtain reliable regional data, although ROI tracking is less influenced by artifacts and are easier to identify as inappropriate ROI tracking. ROI tracking in 2DSE is discussed in the following section.
1) The ROI is defined by tracing the endocardium in a still frame (default at end systole).
2) The ROI width is set to match the myocardial thickness. An automated software program calculates the frame-to-frame displacements of the speckle-pattern within the ROI throughout the cardiac cycle.
3) The resulting tracking quality is then scored as either acceptable or non-acceptable. Since the program follows any speckle pattern, tracking of a stationary artifact during the cardiac cycle is not uncommon, and often scored as appropriate tracked. Therefore, the tracking has to be visually checked and adjusted if necessary. Even the slightest failure of ROI tracking can result in drifting of the calculated strain curve.
For each myocardial segment, velocity, displacement, strain and SR are calculated in two dimensions: longitudinal and transverse parameters in the apical recordings and circumferential and radial parameters in the parasternal recordings (Figure 10). For the right ventricle, only longitudinal parameters can be reliably calculated due to the thin wall. Additionally, ventricular rotation and rotation-rate are calculated for the short axis recordings. With these parameters, left ventricular torsion can be calculated. Likewise for the TDI parameters, temporal and spatial smoothing can be adjusted.
Comparability between TDI and 2DSE
Both TDI and 2DSE derived parameters are accurate quantitative measures of local longitudinal myocardial deformation, and thus (dys)function and will general yield comparable values for local deformation and deformation rates in both the LV  and the RV . Variability in the measurements and technical factors however, make that minor differences between the calculated parameters should be anticipated when comparing values of TDI and 2DSE in a single patient.
Most of the observed differences between the two techniques are those inherent to the limitations described earlier. For example, TDI is not possible when the investigated wall is not optimally aligned, or in the presence of a stationary artifact, whereas 2DSE is less influence by these interferences. On the other hand, 2DSE has a relative low temporal resolution hindering tracking in the presence of high heart rates where also undersampling becomes an issue , this is not a problem with TDI, where frame-rates are higher (up to 250 FPS).
Practical differences between the two techniques are an additional factor when choosing between one or the other. Training requirements for both image acquisition and post-processing are equal, although TDI is generally perceived as more difficult to master. This is particularly the case for post-processing. In our experience, acquisition time is comparable between the two techniques. Post-processing time, however is significantly shorter with speckle tracking  with approximately 2 minutes for 2DSE and 11 minutes for TDI analysis of 16 segments from 3 apical views. This mainly results from a more rapid tracking and parameter extraction in 2DSE, which is automated. The number of segments which can be analyzed with either technique is comparable (Table 2, feasibility), although some studies reported a lower feasibility for 2DSE compared to TDI . The reproducibility of 2DSE is generally better in the LV than with TDI (Table 2). In the RV, this difference is not seen. Despite the advantages 2DSE offers compared to TDI, we feel that both techniques are of equal or complementary value in the quantification of myocardial function.
Clinical and experimental research on the potential role of deformation imaging in various clinical settings is currently an on going topic. The additional value of deformation imaging in echocardiography has been proposed in various pathological conditions although more research is needed  In addition, tissue deformation imaging has given us recent new insights into ventricular function, adaptation and mal-adaptation in response to pathology. In the following section we present some clinical examples (case reports) of tissue deformation imaging (either TDI or 2DSE) in various cardiac pathology, underlining the additional advantage in the field of echocardiography.
Hypertrophic (obstructive) cardiomyopathy (HCM) is a genetic disorder characterized by hypertrophy of the myocardium. The differentiation of HCM from hypertensive left ventricular hypertrophy (H-LVH) on the basis of morphological information obtained by conventional echocardiography is cumbersome. Myocardial disarray is a typical histopathological finding in HCM and provides a substrate for potential lethal ventricular arrhythmias. This disarray also induces functional abnormalities, which could be detected using tissue deformation imaging. A recent article found that this technique reliably distinguishes HCM from H-LVH with a cut-off value for peak systolic strain of -10.6% 
In the following example we present an asymptomatic first degree family member (male, 44 years) of a patient with HCM. This patient had the same mutation as the index patient. On conventional echocardiography a septal thickness of 12 mm (upper limit of normal) was measured with no wall motion abnormalities. Deformation imaging (TDI) unmasked a pathologic deformation in the basal segment of the septum with prestrech, a low peak systolic strain (-3.3%) and post-systolic shortening (Figure 11).
A potential clinical role of tissue deformation imaging is for patient selection for cardiac resynchronization therapy. The maximal delay between wall segments can be visualized using different parameters. The application of regional deformation imaging in this patient group has the advantage to be able to carefully calculate timing issues of regional deformation properties of the myocardium, independent of translated displacement induced by adjacent or opposite wall segments. In addition, LV ejection fraction can be estimated using global LV end-systolic strain (Figure 2).
In the following example a patient (male, 54 years) is presented fulfilling CRT criteria (LBBB with a QRS-width > 120 ms, NYHA-class III-IV and an ejection fraction of 22%) who was implanted with a CRT-device. Prior to implantation, we found a septal-to-lateral delay to peak strain value of 425 msec using 2DSE (Figure 12). A bulls-eye plot of peak systolic strain values is presented in figure 13. After 8 months follow up after device implanation, there was a significant improvement of ejection fraction (to 42%) a reduction in LV volume and a dramatic improvement of longitudinal deformation during systole (figure 14). 2DSE revealed a nearly synchronous longitudinal shortening with a septal-to-lateral delay of -60 msec (figure 15).
Right ventricular function
Although clinical and preclinical research of left ventricular function is broadly being explored, little is known about normal and pathologic right ventricular (RV) function and there is an urgent need for more insight in RV function by quantitative techniques  In our own experience, the quantification of regional RV function using either 2DSE or TDI is feasible and reproducible  This technique has a potential place in the assessment and follow up of patients with primary (see below) or secondary RV disease such as pulmonary artery disease .
The following case is a 36 year old woman evaluated for arrhythmogenic right ventricular dysplasia (ARVD), an inherited cardiomyopathy characterized by fibro-fatty myocardial replacement predominantly in the right ventricle. She has a first degree family member diagnosed with ARVD. No genetic defect was found in either. The RV was not dilated, no wall motion abnormalities could be detected nor any other typical ARVD findings on conventional echocardiographic examination . Tissue Doppler strain analysis as well as 2DSE revealed an abnormal deformation pattern in the basal segment (Figure 16/Additional File 3). Based on this finding, this patient scored a minor criterion (regional RV hypokinesia) fulfilling the diagnostic criteria for ARVD (one major and two minor criteria).
Several new developments are to be expected, not only to improve tracking quality, reduction noise and improve post processing time , but also to make this promising technique more clinically appealing. New automated programs for 2DSE ("automated function imaging"), for instance, allow rapid evaluation of left ventricular deformation in a 17 segment model and visualization of peak systolic strain and post systolic strain index in a bulls-eye figure.
Angle dependency and the impact of artifacts in TDI derived parameters and the low temporal resolution in 2DSE are important limiting factors in the technology. A combination of the two techniques (high quality, angle independent gray scale data and Doppler data with high temporal resolution) could enhance the robustness of the technique and reduce the post-processing time.
MRI-tagging is a technique to calculate deformation in three dimensions. As explained earlier, new echocardiographic techniques have "upgraded" the one-dimensional tissue Doppler to 2DSE. A next step in development is the calculation of strain in three dimensions. With the recent advances in three-dimension echocardiography, this could be possible in the near future. Currently, the low temporal resolution is a severe limitation which has to be overcome before accurate calculations of deformation parameters are feasible.
Tissue deformation imaging is a rapidly evolving technique able to objectively evaluate regional myocardial function. In this review we aimed to give insight into the basic concepts, image acquisition, data analysis and potential clinical applications. Although this quantification provides both the clinician as well as the researcher with detailed information on the myocardial function during the cardiac cycle, it should be kept in mind that, like all ultrasound techniques, image acquisition and quality remain a limiting factor. Although tissue deformation imaging is a promising technique, knowledge of image acquisition, post processing techniques and interpretation of strain and strain-rate graphs are essential to reliably differentiate myocardial pathology from artifacts both in the clinical setting, and for reliable and reproducible research. Its limitations notwithstanding, tissue deformation imaging is an accurate and robust techniques that may successfully fulfill the need for quantitative assessment of myocardial function. The application of this technique in the clinical setting is currently not routine, but the potential of tissue deformation imaging and the need for quantitative measurements of ventricular function should promote its introduction in various clinical settings in the near future.
Picano E, Lattanzi F, Orlandini A, Marini C, L'Abbate AE: Stress echocardiography and the human factor: the importance of being expert. Journal of the American College of Cardiology. 1991, 17: 666-669.
Hoffmann R, Lethen H, Marwick T, Arnese M, Fioretti P, Pingitore A, Picano E, Buck T, Erbel R, Flachskampf FA, Hanrath P: Analysis of interinstitutional observer agreement in interpretation of dobutamine stress echocardiograms. Journal of the American College of Cardiology. 1996, 27: 330-336. 10.1016/0735-1097(95)00483-1.
Weidemann F, Jamal F, Sutherland GR, Claus P, Kowalski M, Hatle L, De Scheerder I, Bijnens B, Rademakers FE: Myocardial function defined by strain rate and strain during alterations in inotropic states and heart rate. Am J Physiol Heart Circ Physiol. 2002, 283: H792-H799.
Greenberg NL, Firstenberg MS, Castro PL, Main M, Travaglini A, Odabashian JA, Drinko JK, Rodriguez LL, Thomas JD, Garcia MJ: Doppler-Derived Myocardial Systolic Strain Rate Is a Strong Index of Left Ventricular Contractility. Circulation. 2002, 105: 99-105. 10.1161/hc0102.101396.
Kjaergaard J, Korinek J, Belohlavek M, Oh JK, Sogaard P, Hassager C: Accuracy, Reproducibility, and Comparability of Doppler Tissue Imaging by Two High-end Ultrasound Systems. Journal of the American Society of Echocardiography. 2006, 19: 322-328. 10.1016/j.echo.2005.10.007.
Urheim S, Edvardsen T, Torp H, Angelsen B, Smiseth OA: Myocardial Strain by Doppler Echocardiography : Validation of a New Method to Quantify Regional Myocardial Function. Circulation. 2000, 102: 1158-1164.
Skulstad H, Urheim S, Edvardsen T, Andersen K, Lyseggen E, Vartdal T, Ihlen H, Smiseth OA: Grading of Myocardial Dysfunction by Tissue Doppler Echocardiography: A Comparison Between Velocity, Displacement, and Strain Imaging in Acute Ischemia. Journal of the American College of Cardiology. 2006, 47: 1672-1682. 10.1016/j.jacc.2006.01.051.
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, 106: 50-56. 10.1161/01.CIR.0000019907.77526.75.
Boettler P, Claus P, Herbots L, McLaughlin M, D'hooge J, Bijnens B, Ho SY, Kececioglu D, Sutherland GR: New aspects of the ventricular septum and its function: an echocardiographic study. Heart. 2005, 91: 1343-1348. 10.1136/hrt.2004.043422.
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. Journal of the American College of Cardiology. 2005, 45: 2034-2041. 10.1016/j.jacc.2005.02.082.
Helle-Valle T, Crosby J, Edvardsen T, Lyseggen E, Amundsen BH, Smith HJ, Rosen BD, Lima JAC, Torp H, Ihlen H, Smiseth OA: New Noninvasive Method for Assessment of Left Ventricular Rotation: Speckle Tracking Echocardiography. Circulation. 2005, 112: 3149-3156. 10.1161/CIRCULATIONAHA.104.531558.
Korinek J, Wang J, Sengupta PP, Miyazaki C, Kjaergaard J, McMahon E, Abraham TP, Belohlavek M: Two-Dimensional Strain-A Doppler-Independent Ultrasound Method for Quantitation of Regional Deformation: Validation In Vitro and In Vivo. Journal of the American Society of Echocardiography. 2005, 18: 1247-1253. 10.1016/j.echo.2005.03.024.
Korinek J, Kjaergaard J, Sengupta PP, Yoshifuku S, McMahon EM, Cha SS, Khandheria BK, Belohlavek M: High Spatial Resolution Speckle Tracking Improves Accuracy of 2-Dimensional Strain Measurements: An Update on a New Method in Functional Echocardiography. Journal of the American Society of Echocardiography. 2007, 20: 165-170. 10.1016/j.echo.2006.08.031.
Amundsen BH, Helle-Valle T, Edvardsen T, Torp H, Crosby J, Lyseggen E, Stoylen A, Ihlen H, Lima JAC, Smiseth OA, Slordahl SA: Noninvasive Myocardial Strain Measurement by Speckle Tracking Echocardiography: Validation Against Sonomicrometry and Tagged Magnetic Resonance Imaging. Journal of the American College of Cardiology. 2006, 47: 789-793. 10.1016/j.jacc.2005.10.040.
Becker M, Bilke E, Kuhl H, Katoh M, Kramann R, Franke A, Bucker A, Hanrath P, Hoffmann R: Analysis of myocardial deformation based on pixel tracking in two dimensional echocardiographic images enables quantitative assessment of regional left ventricular function. Heart. 2006, 92: 1102-1108. 10.1136/hrt.2005.077107.
Leitman M, Lysyansky P, Sidenko S, Shir V, Peleg E, Binenbaum M, Kaluski E, Krakover R, Vered Z: Two-dimensional strain-a novel software for real-time quantitative echocardiographic assessment of myocardial function. Journal of the American Society of Echocardiography. 2004, 17: 1021-1029. 10.1016/j.echo.2004.06.019.
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. The American Journal of Cardiology. 2006, 97: 1661-1666. 10.1016/j.amjcard.2005.12.063.
Teske AJ, De Boeck BWL, Olimulder M, Prakken NH, Doevendans PA, Cramer MJM: Echocardiographic assessment of regional right ventricular function. A head to head comparison between 2D-strain and tissue Doppler derived strain analysis. Journal of the American Society of Echocardiography. 2007, in press:
Behar V, Adam D, Lysyansky P, Friedman Z: Improving motion estimation by accounting for local image distortion. Ultrasonics. 2004, 43: 57-65. 10.1016/j.ultras.2004.02.022.
Hanekom L, Cho GY, Leano R, Jeffriess L, Marwick TH: Comparison of two-dimensional speckle and tissue Doppler strain measurement during dobutamine stress echocardiography: an angiographic correlation. Eur Heart J. 2007, 28 (14): 1765-1772. 10.1093/eurheartj/ehm188.
Ingul CB, Torp H, Aase SA, Berg S, Stoylen A, Slordahl SA: Automated Analysis of Strain Rate and Strain: Feasibility and Clinical Implications. Journal of the American Society of Echocardiography. 2005, 18: 411-418. 10.1016/j.echo.2005.01.032.
Marwick TH: Measurement of Strain and Strain Rate by Echocardiography: Ready for Prime Time?. Journal of the American College of Cardiology. 2006, 47: 1313-1327. 10.1016/j.jacc.2005.11.063.
Kato TS, Noda A, Izawa H, Yamada A, Obata K, Nagata K, Iwase M, Murohara T, Yokota M: Discrimination of Nonobstructive Hypertrophic Cardiomyopathy From Hypertensive Left Ventricular Hypertrophy on the Basis of Strain Rate Imaging by Tissue Doppler Ultrasonography. Circulation. 2004, 110: 3808-3814. 10.1161/01.CIR.0000150334.69355.00.
Voelkel NF, Quaife RA, Leinwand LA, Barst RJ, McGoon MD, Meldrum DR, Dupuis J, Long CS, Rubin LJ, Smart FW, Suzuki YJ, Gladwin M, Denholm EM, Gail DB: Right Ventricular Function and Failure: Report of a National Heart, Lung, and Blood Institute Working Group on Cellular and Molecular Mechanisms of Right Heart Failure. Circulation. 2006, 114: 1883-1891. 10.1161/CIRCULATIONAHA.106.632208.
Pirat B, McCulloch ML, Zoghbi WA: Evaluation of Global and Regional Right Ventricular Systolic Function in Patients With Pulmonary Hypertension Using a Novel Speckle Tracking Method. The American Journal of Cardiology. 2006, 98: 699-704. 10.1016/j.amjcard.2006.03.056.
Yoerger DM, Marcus F, Sherrill D, Calkins H, Towbin JA, Zareba W, Picard MH: Echocardiographic findings in patients meeting task force criteria for arrhythmogenic right ventricular dysplasia: New insights from the multidisciplinary study of right ventricular dysplasia. Journal of the American College of Cardiology. 2005, 45: 860-865. 10.1016/j.jacc.2004.10.070.
Bjork Ingul C, Stoylen A, Slordahl SA, Wiseth R, Burgess M, Marwick TH: Automated Analysis of Myocardial Deformation at Dobutamine Stress Echocardiography: An Angiographic Validation. Journal of the American College of Cardiology. 2007, 49: 1651-1659. 10.1016/j.jacc.2007.01.059.
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. The American Journal of Cardiology. 2006, 98: 1581-1586. 10.1016/j.amjcard.2006.07.038.
Kowalski M, Kukulski T, Jamal F, D'hooge J, Weidemann F, Rademakers F, Bijnens B, Hatle L, Sutherland GR: Can natural strain and strain rate quantify regional myocardial deformation? A study in healthy subjects. Ultrasound in Medicine & Biology. 2001, 27: 1087-1097. 10.1016/S0301-5629(01)00388-X.
De Boeck BWL, Cramer MJ, Loh P, Doevendans PAFM: Two-Dimensional Strain Imaging to Assess the Origin and Extent of Ventricular Preexcitation Associated With an Accessory Bypass. Circulation. 2006, 113: e835-e839. 10.1161/CIRCULATIONAHA.105.589762.
Borges AC, Knebel F, Eddicks S, Panda A, Schattke S, Witt C, Baumann G: Right Ventricular Function Assessed by Two-Dimensional Strain and Tissue Doppler Echocardiography in Patients With Pulmonary Arterial Hypertension and Effect of Vasodilator Therapy. The American Journal of Cardiology. 2006, 98: 530-534. 10.1016/j.amjcard.2006.02.060.
Ingul CB, Rozis E, Slordahl SA, Marwick TH: Incremental Value of Strain Rate Imaging to Wall Motion Analysis for Prediction of Outcome in Patients Undergoing Dobutamine Stress Echocardiography. Circulation. 2007, 115: 1252-1259.
The author(s) declare that they have no competing interests.
AT designed the review, participated to analyzed paper collection and drafted the paper.
BB and PM helped to draft the manuscript (2DSE section and Doppler section respectively).
GS and PD revised the article.
MC helped to draft the manuscript.
All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Impact of stationary artifact on 2DSE. Appropriate ROI in a 4-chamber view using 2DSE despite a prominent stationary artifact in the lateral wall. The ROI is color coded: yellow-red = negative SR, green = no SR and blue = positive SR. In the region of the artifact there is no reversal of SR values, all myocardial segments show identical color coding during each phase of the cardiac cycle. See figure 9 for additional information. (AVI 4 MB)
Additional file 3: Deformation imaging of the right ventricle in a patient suspected for ARVD using 2DSE. 2DSE parameter extraction in the right ventricle. A ROI is tracked along the RV free wall and follows its movement during the cardiac cycle. An overlying color map is seen on the ROI: red indicates normal deformation; pink indicates reduced deformation and blue indicates no deformation, or stretch. See figure 16 for additional information. (AVI 6 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Teske, A.J., De Boeck, B.W., Melman, P.G. et al. Echocardiographic quantification of myocardial function using tissue deformation imaging, a guide to image acquisition and analysis using tissue Doppler and speckle tracking. Cardiovasc Ultrasound 5, 27 (2007). https://doi.org/10.1186/1476-7120-5-27