- Research article
- Open Access
Diaphragm assessment by two dimensional speckle tracking imaging in normal subjects
BMC Anesthesiology volume 16, Article number: 43 (2015)
Conventionally, ultrasonographic assessment of diaphragm contractility has involved measuring respiratory changes in diaphragm thickness (thickening fraction) using B-mode or caudal displacement with M-mode. Two-dimensional speckle-tracking has been increasingly used to assess muscle deformation (‘strain’) in echocardiography. We sought to determine in a pilot study if this technology could be utilized to analyze diaphragmatic contraction.
Fifty healthy adult volunteers with normal exercise capacity underwent ultrasound imaging. A linear array transducer was used for the assessment of diaphragm thickness, thickening fraction (TF), and strain in the right anterior axillary line at approximately the ninth intercostal space. A phased array transducer was applied subcostally for the assessment of diaphragm displacement on the right mid-clavicular line. Diaphragmatic images were recorded from the end of expiration through the end of inspiration at 60 % maximal inspiratory capacity. Diaphragm strain was analyzed off-line by speckle tracking imaging. Blinded inter- and intra-rater variability was tested in 10 cases.
Mean right diaphragm thickness at end-expiration (±SD: standard deviation) was 0.24 cm (±0.1), with TF of 45.1 % (±12) at 60 % peak inspiratory effort. Mean right diaphragm caudal displacement was 4.9 cm (±1). Mean right diaphragm strain was -40.3 % (±9). A moderate correlation was seen between longitudinal strain and TF (R2 0.44, p < 0.0001). A weak correlation was seen between strain and caudal displacement (R2 0.14, p < 0.01), and an even weaker correlation was seen between caudal displacement and TF (R2 0.1, p = 0.04). Age, gender, and body mass index were not significantly associated with right diaphragm strain or TF. Although inter- and intra-rater variability was overall good for TF, caudal displacement, and strain (inter-rater R2; 0.8, 0.9, and 0.7, respectively [p < 0.01], intra-rater R2; 0.9, 0.7, and 0.9, respectively [p < 0.01]), strain values did have a slightly lower inter-rater repeatability.
Diaphragmatic strain estimated by speckle tracking imaging was associated with conventional ultrasound measures of diaphragmatic function (TF and caudal displacement). Further clinical studies are warranted to investigate its clinical utility.
Measurement of transdiaphragmatic pressures has been considered the gold standard for the assessment of diaphragmatic function; however, it is rarely conducted due to its invasiveness . While several alternative modalities have been used, such as chest radiography, video fluoroscopy, pulmonary function testing, electromyography, and nerve conduction studies, their accuracy and reproducibility may be modest and some are invasive or cause radiation exposure . Recently, there has been growing interest in ultrasonographic assessment of the diaphragm due to its portability, lack of radiation exposure, non-invasiveness, and ease of repeatability. Previous studies have demonstrated the utility of diaphragm thickening fraction (TF) and caudal displacement with inspiration as markers of diaphragmatic function, measured via M-mode in the zone of diaphragm apposition and via two-dimensional (2D) B-mode imaging in the subcostal area, respectively [3–6] (Figs. 1 and 2). Measurements have been shown to be feasible, reproducible, and accurate [7–9]. Other studies have shown that a large variety of disorders can result in diaphragm thinning and dysfunction, which was associated with respiratory complications and prolonged weaning from mechanical ventilation [10, 11].
Measuring diaphragm thickening and displacement are indirect estimates of diaphragm muscle contraction as they do not assess ‘longitudinal’ muscle shortening, ie: in the ‘plane’ of muscle fiber motion. Speckle tracking analysis has the potential to describe this function. Speckle tracking is a novel analysis method used in echocardiography to determine myocardial fiber deformation . Ultrasound images are made up of different grey-scale pixels, and the speckle tracking software follows unique groups of these pixels (known as ‘kernels’) and measures their displacement and how different ‘kernels’ move in relation to one another (deformation). The degree of deformation is known as ‘strain’ and negative values indicate ‘kernels’ coming closer together. For example, a strain value of -40 % indicates local muscle fiber shortening of 40 %. The more negative a number, the greater the degree of deformation and the greater the contraction (Fig. 3).
In this pilot study, we sought to determine (1) the feasibility of speckle tracking in assessing diaphragm function, (2) how this compares to conventional methods of diaphragm function assessment by ultrasound: TF and caudal displacement, and (3) reproducibility of speckle tracking analysis of the diaphragm.
Fifty adult subjects were included in the study. Inclusion criteria comprised of healthy adults at age over 18 years with a normal exercise tolerance (defined as being able to climb three flights of stairs with relative ease). Exclusion criteria were previous history of neuromuscular disease, chronic respiratory issues (defined as moderate to severe restrictive or obstructive lung disease), and cardiopulmonary or abdominal surgery. Ultrasonography was conducted in subjects during their inspiratory efforts. Each inspiration effort was standardized by using an incentive spirometry device (Hudson RCI®, Voldyne 5000). Using a predictive normogram based on sex, age and height , peak inspiratory capacity was calculated for each individual. Study subjects were instructed to inspire to 60 % of their peak inspiratory capacity in one second. A 60 % inspiratory breath was chosen to limit descent of the lung into the image but to ensure adequate diaphragm movement. The study was approved by the Mayo Institutional Review Board (11-002744), and all subjects provided informed consent.
Ultrasound imaging and analysis
Imaging was performed using a commercially available ultrasound machine (Vivid E9, General Electric Healthcare, Milwaukee, WI) with a linear array transducer (2.5-8 MHz) and a phased array transducer (1.6-6 MHz) by S.O.: Australian Intensive Care specialist, board certified in standard and advanced echocardiography in America (Fellow of American Society of Echocardiography, FASE) and Australia (Diploma Diagnostic Ultrasound, DDU). Subjects were positioned semi-recumbent (45° head-up). The linear array transducer was first placed perpendicular to the angle of the lower ribs, with the transducer marker directed cephalad. Images of two adjacent ribs were seen at approximately the same level with the diaphragm running parallel to the transducer surface (Fig. 1). Image acquisition began in the right mid-axillary line above the costal margin, where the zone of apposition of the diaphragm could be visualized at approximately the 9th intercostal space (Additional file 1: Figure S1a). If the lung expansion obscured the optimal images during inspiration (‘lung curtain’ sign) , then the transducer was moved anteriorly along the rib space toward the anterior axillary line until no ‘lung curtain’ sign was seen (Additional file 1: Figure S1b). Only one focus level was chosen, and the depth was optimized in order to visualize the diaphragm as clearly as possible throughout inspiration with maximal frame rate. Attempts were made to enhance sonographic speckles in the diaphragm by slight transducer angulation.
M-mode imaging was conducted with the cursor line equidistant between the two ribs perpendicular to the diaphragm. Diaphragm thickness was measured, excluding the outer layers of pleura and peritoneum, from the inner-edge to the inner-edge at end-expiration and at end-inspiration (Fig. 2a). TF was calculated by: diaphragm diameter at end-inspiration minus diaphragm diameter at end-expiration divided by diaphragm diameter at end-expiration, expressed as a percentage: [(DDinsp – DDexp)/DDexp] x 100.
Subsequently, three-second clips were recorded for off-line speckle-tracking analysis using EchoPacs (GE Healthcare, Milwaukee, MI). Various assumptions that EchoPacs requires for standard speckle-tracking were contradicted, for example: absent ECG tracings, use of linear array transducer, shallow depth, abnormal tracing, lack of drift compensation, and optimized spatial and temporal smoothing. Apical four chamber longitudinal analysis was used, and the inner surface of the two hyperechogenic lines bordering the diaphragm muscle (the pleura and peritoneal surface) were traced and followed during an inspiratory breath to contain the diaphragm muscle within the barriers. The central portion of the region of interest was analyzed only (Figs. 3 and 4, and Additional file 2: Video 1).
Diaphragm caudal displacement was calculated by imaging with a phased array transducer in the right anterior to mid-clavicular line with the transducer aimed dorso-cranially (Additional file 3: Figure S2). M-mode imaging was used with the cursor line placed at the dome of the diaphragm. Diaphragm displacement was calculated as the maximal height change of the echogenic line on the M-mode image (Fig. 2b).
In order to assess inter-rater variability, 10 subjects’ images were analyzed by two different observers (S.O. and D.F.) blinded to each other’s assessment. Intra-rater variability was evaluated in 10 subjects by performing the same imaging and analysis on a different examination date (by S.O.).
Data analysis was performed using JMP software version 10.0.0 (SAS Institute Inc. North Carolina). Continuous variables were expressed as mean ± standard deviation (SD) or median ± interquartile range (IQR) if not normally distributed. Normality was tested by visual assessment and the Shapiro-Wilk test . Non-normally distributed data (eg: right diaphragm longitudinal strain) underwent logarithmic transformation. Pearson’s correlation coefficient was used to assess correlation between continuous variables and t-test for direct comparison. P values <0.05 were considered statistically significant, and all probability values were 2 sided. Reproducibility was expressed by interclass correlation coefficient , Pearson’s correlation, and Bland-Altman’s coefficient of repeatability (twice the SD of differences in repeated measurements) [17, 18].
A total of 50 subjects were prospectively enrolled for the study. Twenty-eight (56 %) were female with a median age of 37 years old (IQR: 30.2 to 39.8). The median body mass index (BMI) was 22.8 kg/m2 (IQR 20.4–24.9). Estimated inspiratory capacity based on age, sex, and height had a median of 2900 ml (IQR 2375 to 3100). Baseline characteristics of the subjects by sex are described in Table 1.
The mean diaphragm thickness at end-expiration, TF, caudal displacement, and longitudinal strain values are listed in Table 2. The mean right diaphragmatic thickness was 0.24 cm (±0.1) and 0.35 cm (±0.1) at the end expiration and inspiration, respectively. The mean TF was 45.1 % (±12). The mean caudal displacement was 4.9 cm (±1). The right diaphragmatic longitudinal strain was -40.3 % (±9). No sex-specific differences were seen in these values. Neither BMI nor estimated inspiratory capacity demonstrated a significant correlation with diaphragmatic thickness, TF, caudal displacement, or strain on the univariate linear regression analysis. There was a moderate association between right diaphragm longitudinal strain (logarithmic scale) and TF (R2 0.44, p < 0.0001). Although a weaker association with right diaphragm strain and caudal displacement was observed (R2 0.14, p < 0.01), its visual correlation was not as strong as one between TF and caudal displacement (R2 0.1, p < 0.0001) (see Fig. 5).
Inter-and intra-operator reliability was good (Table 3). The correlation between the measurements obtained by two operators was highly significant in the right diaphragm TF (ICC 0.9, R2 0.8, p < 0.001, coefficient of repeatability 2.6 % and 8 % respectively) as well as caudal displacement (ICC 0.9, R2 0.9, p < 0.001, coefficient of repeatability 14.5 %). The correlation in the right diaphragm longitudinal strain measurements was also significant (ICC 0.9, R2 0.7, p = 0.004), with a slightly higher coefficient of repeatability (24.3 %). A similar trend was seen in intra-operator variability. The correlation in measurements of TF, caudal displacement, and strain were highly significant (R2 0.94, 0.72, and 0.91, respectively).
Our study demonstrated that diaphragmatic assessment using speckle tracking with a linear array transducer was feasible and reproducible, albeit more variable than conventional assessment by TF with use of M-mode ultrasonography. This supports a recent study that reported 2D speckle tracking assessment using a low frequency phased array transducer placed on the right subcostal, mid-clavicular line . In this study, longitudinal muscle contraction of the diaphragm was estimated by assessing the entire dome of the diaphragm, whereas we examined the diaphragm in the zone of apposition with use of a high frequency linear array transducer. We found a mean right diaphragm longitudinal strain value of -40.3 % ±9 at a 60 % maximal inspiratory capacity. This corresponded with a mean TF of 45.1 % (±12), which was in keeping with data from a larger population of normal subjects . Mean right diaphragm thickness at end-expiration in our study was 0.24 mm ±0.1 and caudal displacement 4.9 cm (±1), which are also comparable to previous literature [6, 20].
Although there was a moderate correlation seen between right diaphragm longitudinal strain and TF, the correlation between strain and caudal displacement was weak. The significance of correlation, albeit not a strong one, suggests that at least speckle tracking may be reasonable method for analysing diaphragm function. We have to be mindful of the fact that strain and TF are measuring function in different planes and direction: for example diaphragm longitudinal strain could be considered to measure deformation in the ‘x plane’, while TF is a measure of deformation in the ‘y plane’. However, a modest correlation between strain and TF can be explained through the principle of conservation of volume, as previously proposed in the assessment of incompressible tissue . Recent studies suggested that diaphragm displacement is not as representative of respiratory effort or diaphragmatic function as diaphragm thickening [22, 23]. There was no strong correlation found between diaphragm displacement and indices of respiratory workload, nor was between displacement and TF . Interestingly, in our study, a poor correlation was seen between strain and caudal displacement.
A significant difference was seen between male and female subjects for BMI and 60 % inspiratory capacity; however, this difference was not seen with respect to diaphragm function. Our result may represent the hypothesis that not the absolute value of inspiratory capacity, but the proportion (or percentage, here 60 %) of the maximal inspiratory capacity is an important factor to dictate the needed diaphragmatic strength and contraction when adjusted by sex and body size.
The diaphragm acts as a primary muscle of ventilation, and its dysfunction can result in dyspnoea and complications such as difficultly in liberation from mechanical ventilation , yet it is not routinely evaluated or monitored . Ultrasound has been recognized as a new tool to examine the diaphragm non-invasively in various clinical settings [25, 26]. However, conventional methods, such as TF, or caudal displacement assessed by M-mode have limitations, such as angle dependence and translational error (where other areas of diaphragm move into the ‘line of sight’). 2D speckle tracking follows specific grey-scale pixels (known as ‘kernels’) in specified regions of interest and is less prone to the errors of traditional ultrasound methods as it is relatively angle independent. Also tracking can help differentiate passive movement from active contraction as the ‘kernels’ do not move closer together in passive motion, whereas they do with any active contraction . Our study may help to provide a novel quantifiable method of diaphragm function analysis: longitudinal strain of the diaphragm assessed by speckle tracking.
There are limitations to our study. Primarily, the reference ultrasound methods used in this study for measuring diaphragm function have not been validated in a large study or in a study using trans-diaphragmatic pressure as a gold standard [3, 28]. However, diaphragm TF and caudal displacement have been shown to be relevant in the clinical environment and are commonly used non-invasive tools . Secondly, the technology used in our study is not designed for diaphragm analysis. Several assumptions that the software makes which help it to be sensitive in determining myocardial deformation needed to be broken. Third, the speckle tracking software automatically generated three regions of interest when the diaphragm was traced. We only analysed the central region of interest as the peripheral sections did not have appropriate tracking. Fourth, the time taken to perform analysis was considerable relative to a clinical situation. This limits the current clinical utility of this method of assessment. However, as technology improves, speckle tracking may become easier and has the potential to add to non-invasive diaphragm function assessment. Finally, our pilot study enrolled a small number of patients and did not include subjects with abnormal diaphragmatic function. Sample size was difficult to statistically estimate due to lack of prior data. Further studies enrolling subjects with abnormal diaphragmatic function and/or critical illness, with comparison to electromyography or transdiaphragmatic pressure measurements may be warranted.
2D speckle tracking analysis of the right diaphragm is feasible and reproducible. It is associated with conventional ultrasound measures of diaphragmatic function (TF and caudal displacement). It provides unique information on longitudinal muscle deformation (strain) which conventional methods of ultrasound analysis, such as M-mode, are unable to provide. Although this form of analysis can take considerable time to perform, technological progress may help make this a useful clinical tool in the future.
2D, Two dimensional; DDexp, Diaphragm diameter end-expiration; DDinsp, Diaphragm diameter end-inspiration; ICU, Intensive Care Unit; IQR, Interquartile range; SD, Standard deviation; TF, Thickening fraction.
Zambon M, Cabrini L, Zangrillo A. Diaphragmatic Ultrasound in Critically Ill Patients. Annual update in Intensive Care and Emergency Medicine. 2013 2013:427–438.
McCool FD, Tzelepis GE. Dysfunction of the Diaphragm. N Engl J Med. 2012;366:932–42.
Wait JL, Nahormek PA, Yost WT, Rochester DP. Diaphragmatic thickness-lung volume relationship in vivo. J Appl Physiol. 1989;67:1560–8.
Summerhill EM, Angov N, Garber C, McCool FD. Respiratory muscle strength in the physically active elderly. Lung. 2007;185:315–20.
Summerhill EM, El-Sameed YA, Glidden TJ, McCool FD. Monitoring recovery from diaphragm paralysis with ultrasound. Chest. 2008;133:737–43.
Boon AJ, Harper CJ, Ghahfarokhi LS, Strommen JA, Watson JC, Sorenson EJ. Two-dimensional ultrasound imaging of the diaphragm: quantitative values in normal subjects. Muscle Nerve. 2013;47:884–9.
Harper CJ, Shahgholi L, Cieslak K, Hellyer NJ, Strommen JA, Boon AJ: Variability in Diaphragm Motion During Normal Breathing, Assessed With B Mode Ultrasound. J Orthop Sports Phys Ther. 2013
Gottesman E, McCool FD. Ultrasound evaluation of the paralyzed diaphragm. Am J Respir Crit Care Med. 1997;155:1570–4.
Matamis D, Soilemezi E, Tsagourias M, Akoumianaki E, Dimassi S, Boroli F, Richard J-CM, Brochard L. Sonographic evaluation of the diaphragm in critically ill patients. Technique and clinical applications. Intensive Care Med. 2013;39:801–10.
Jiang J-R, Tsai T-H, Jerng J-S, Yu C-J, Wu H-D, Yang P-C. Ultrasonographic evaluation of liver/spleen movements and extubation outcome. Chest. 2004;126:179–85.
Kim WY, Suh HJ, Hong S-B, Koh Y, Lim C-M. Diaphragm dysfunction assessed by ultrasonography: influence on weaning from mechanical ventilation. Crit Care Med. 2011;39:2627–30.
Huang SJ, Orde S. From speckle tracking echocardiography to torsion: research tool today, clinical practice tomorrow. Curr Opin Crit Care. 2013;19:250–7.
Polgar G, Weng TR. The functional development of the respiratory system from the period of gestation to adulthood. Am Rev Respir Dis. 1979;120:625–95.
Doerschug KC, Schmidt GA. Intensive care ultrasound: III. Lung and pleural ultrasound for the intensivist. Ann Am Thorac Soc. 2013;10:708–12.
Ghasemi A, Zahediasl S. Normality Tests for Statistical Analysis: A Guide for Non-Statisticians. Int J Endocrinol Metab. 2012;10:486–9.
Bland JM, Altman DG. Measurement error. BMJ. 1996;312:1654.
Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;1:307–10.
Rothwell PM. Analysis of agreement between measurements of continuous variables: general principles and lessons from studies of imaging of carotid stenosis. J Neurol. 2000;247:825–34.
Ye X, Xiao H, Bai W, Liang Y, Chen M, Zhang S. Two-dimensional strain ultrasound speckle tracking as a novel approach for the evaluation of right hemidiaphragmatic longitudinal deformation. Exp Ther Med. 2013;6:368–72.
Gerscovich EO, Cronan M, McGahan JP, Jain K, Jones CD, McDonald C. Ultrasonographic evaluation of diaphragmatic motion. J Ultrasound Med. 2001;20:597–604.
Heimdal A, Stoylen A, Torp H, Skjaerpe T. Real-time strain rate imaging of the left ventricle by ultrasound. J Am Soc Echocardiogr. 1998;11:1013–9.
Umbrello M, Formenti P, Longhi D, Galimberti A, Piva I, Pezzi A, Mistraletti G, Marini JJ, Iapichino G. Diaphragm ultrasound as indicator of respiratory effort in critically ill patients undergoing assisted mechanical ventilation: a pilot clinical study. Crit Care. 2015;19:161.
Vivier E, Mekontso Dessap A, Dimassi S, Vargas F, Lyazidi A, Thille AW, Brochard L. Diaphragm ultrasonography to estimate the work of breathing during non-invasive ventilation. Intensive Care Med. 2012;38:796–803.
Nason LK, Walker CM, McNeeley MF, Burivong W, Fligner CL, Godwin JD. Imaging of the diaphragm: anatomy and function. Radiographics. 2012;32:E51–70.
Noh DK, Lee JJ, You JH. Diaphragm breathing movement measurement using ultrasound and radiographic imaging: a concurrent validity. Biomed Mater Eng. 2014;24:947–52.
Boussuges A. Diaphragmatic Motion Studied by M-Mode Ultrasonography. Chest. 2009;135:391.
Seo Y, Ishizu T, Sakamaki F, et al. Analysis of the origin of cardiac wall motion that constitutes myocardial velocity-time curves in patients with left bundle branch block. J Am Soc Echocardiogr. 2009;22(4):331–6.
Cohn D, Benditt JO, Eveloff S, McCool FD. Diaphragm thickening during inspiration. J Appl Physiol. 1997;83:291–6.
The authors would like to thank Mr Andrew J. Zafft (R.EP.T., R.NCS.T.) and Tammy Green (MBA RDCS) for their help with subject recruitment. This publication was made possible by the Mayo Clinic CTSA through grant number UL1 TR000135 from the National Center for Advancing Translational Sciences (NCATS), a component of the National Institutes of Health (NIH).
This publication was made possible by the Mayo Clinic CTSA through grant number UL1 TR000135 from the National Center for Advancing Translational Sciences (NCATS), a component of the National Institutes of Health (NIH).
SRO, AJB, DGF, HV, and HS had access to the data, contributed to data acquisition and take responsibility for the integrity and accuracy of the data and analysis. SRO, AJB, HV and HS contributed to study conception and design, statistical analysis and preparation of the manuscript. DGF contributed to data analysis and interpretation, and drafting of the manuscript. All authors read and approved the final manuscript
The authors declare that they have no competing interests.
Availability of data and materials
The dataset supporting the conclusions on this article are available in the ‘figshare’ repository: https://figshare.com/articles/Diaphragm_strain_imaging_data_SO_9_5_2016_xls/3362593.
Consent for publication
Consent to publish was given by the models used in the supplemental figure images.
Ethics approval and consent to participate
The study was approved by the Mayo Institutional Review Board (11-002744), and all subjects provided informed consent to participate in the study.
Linear array transducer application for the assessment of diaphragm thickness and strain. (a) Image acquisition in the right mid-axillary line above the costal margin, where the zone of apposition of the diaphragm can be visualized at approximately the 9th intercostal space. (b) If the lung expansion obscures the optimal images during inspiration (‘lung curtain’ sign), then the transducer can be moved anteriorly along the intercostal space toward the anterior axillary line until no ‘lung curtain’ sign is seen. (PPTX 1272 kb)
Phased array transducer application for the assessment of diaphragm displacement. Diaphragm caudal displacement can be measured by imaging with a phased array transducer in the right anterior to mid-clavicular line with the transducer aimed dorso-cranially. (PPTX 648 kb)
About this article
Cite this article
Orde, S.R., Boon, A.J., Firth, D.G. et al. Diaphragm assessment by two dimensional speckle tracking imaging in normal subjects. BMC Anesthesiol 16, 43 (2015). https://doi.org/10.1186/s12871-016-0201-6
- Speckle tracking