Martin Girard1,2,3*
Marie-Hélène Roy Cardinal4Michaël Chassé2,5Sébastien Garneau1
Yiorgos Alexandros Cavayas6
Guy Cloutier4,7
André Y. Denault8- 1Department of Anesthesiology, University of Montreal Hospital, Montréal, QC, Canada
- 2Division of Critical Care, Department of Medicine, University of Montreal Hospital, Montréal, QC, Canada
- 3University of Montreal Hospital Research Center, Montréal, QC, Canada
- 4Laboratory of Biorheology and Medical Ultrasonics, University of Montreal Hospital Research Center, Montréal, QC, Canada
- 5Department of Medicine, University of Montreal, Montréal, QC, Canada
- 6Department of Medicine, Sacré-Coeur Hospital of Montréal, Montréal, QC, Canada
- 7Department of Radiology, Radio-Oncology and Nuclear Medicine, Institute of Biomedical Engineering, University of Montreal, Montréal, QC, Canada
- 8Department of Anesthesiology, Montreal Heart Institute, Montréal, QC, Canada
Background: Mechanical ventilation is a common therapy in operating rooms and intensive care units. When ill-adapted, it can lead to ventilator-induced lung injury (VILI), which is associated with poor outcomes. Excessive regional pulmonary strain is thought to be a major mechanism responsible for VILI. Scarce bedside methods exist to measure regional pulmonary strain. We propose a novel way to measure regional pleural strain using ultrasound elastography. The objective of this study was to assess the feasibility and reliability of pleural strain measurement by ultrasound elastography and to determine if elastography parameters would correlate with varying tidal volumes.
Methods: A single-blind randomized crossover proof of concept study was conducted July to October 2017 at a tertiary care referral center. Ten patients requiring general anesthesia for elective surgery were recruited. After induction, patients received tidal volumes of 6, 8, 10, and 12 mL.kg–1 in random order, while pleural ultrasound cineloops were acquired at 4 standardized locations. Ultrasound radiofrequency speckle tracking allowed computing various pleural translation, strain and shear components. We screened 6 elastography parameters (lateral translation, lateral absolute translation, lateral strain, lateral absolute strain, lateral absolute shear and Von Mises Strain) to identify those with the best dose-response with tidal volumes using linear mixed effect models. Goodness-of-fit was assessed by the coefficient of determination. Intraobserver, interobserver and test-retest reliability were calculated using intraclass correlation coefficients.
Results: Analysis was possible in 90.7% of ultrasound cineloops. Lateral absolute shear, lateral absolute strain and Von Mises strain varied significantly with tidal volume and offered the best dose-responses and data modeling fits. Point estimates for intraobserver reliability measures were excellent for all 3 parameters (0.94, 0.94, and 0.93, respectively). Point estimates for interobserver (0.84, 0.83, and 0.77, respectively) and test-retest (0.85, 0.82, and 0.76, respectively) reliability measures were good.
Conclusion: Strain imaging is feasible and reproducible. Future studies will have to investigate the clinical relevance of this novel imaging modality.
Clinical trial registration: www.Clinicaltrials.gov, identifier NCT03092557.
Introduction
While often a live-saving therapy and a necessity during general anesthesia, mechanical ventilation may promote lung damage when ill-adapted, a phenomenon called ventilator-induced lung injury (VILI) (1). Some patient populations have been shown to be at higher risk of developing VILI: patients suffering from the acute respiratory distress syndrome (2) and patients requiring one-lung ventilation (3).
Assessing mechanical properties of the lung has been suggested as a promising approach to detect conditions leading to VILI (4) with excessive pulmonary strain thought to be a critical factor in its development (5, 6). Unfortunately, pulmonary strain is non-trivial to measure at the bedside (7) and most techniques measure global pulmonary strain, an oversimplification of the complex lung physiology. Areas of higher regional pulmonary strain that cannot be assessed using global measures have been implicated in the development of pulmonary inflammation (8). Detecting these areas of higher regional pulmonary strain may enable identification of patients at higher risk of VILI who have seemingly benign global strain measures and yet may benefit from further strain-lowering interventions (9).
While magnetic resonance imaging (10), positron emission tomography (11) and thoracic computed tomography (9) are currently available in a research setting to measure regional pulmonary strain, these techniques are particularly burdensome, subject patients to the dangers of intra-hospital transport (12), expose patients to ionizing radiations (13) and haven’t been validated against a gold standard nor have they been shown to be associated with clinical outcome. Electrical impedance tomography is an increasingly popular monitoring modality in ventilated patients (14) and has the advantage of being a bedside and radiation-free imaging modality. Unfortunately, electrical impedance tomography cannot currently measure regional pulmonary strain and requires costly dedicated equipment that is currently unavailable in most operating rooms and intensive care units. Consequently, all of the above modalities have major drawbacks that limit their use and validity.
Lung ultrasonography, a bedside, precise, repeatable, easy to learn, and low-cost exam (15), is already used to monitor pulmonary aeration during mechanical ventilation (16). Preliminary works hint that it may be an interesting avenue to measure regional pleural strain as a surrogate for regional pulmonary strain (17, 18). Ultrasound (US) systems, now standard equipment in most operating rooms and intensive care, could make bedside regional pleural strain measurement commonplace. However, prospective feasibility data is lacking.
The objective of this study was to assess the feasibility and reliability (intraobserver, interobserver, and test-retest) of pleural strain measurements by ultrasound elastography. Moreover, we sought to determine if elastography indices of pleural translation, strain and shear would correlate with varying tidal volumes.
Materials and methods
This report was redacted following the CONSORT 2010 extension to randomized crossover trials statement (19).
Study design
This study was a pilot single-center, single-blind, randomized, four period crossover trial. A crossover design was chosen for this study to improve its power while dropouts and a potential carryover effect were not expected. The protocol was approved by the ethics committee of the Centre hospitalier de l’Université de Montréal (16.386) and registered at ClinicalTrials.gov (NCT03092557, registered March 28th, 2017). Written informed consent was obtained from all study participants.
Study population
Between July and October 2017, adult patients with healthy lungs who were scheduled to undergo an elective surgery under general anesthesia requiring endotracheal intubation and muscle relaxation were screened for inclusion. Patients were recruited at the Centre Hospitalier de l’Université de Montréal (Montréal, Canada), a tertiary care referral center. Patients with healthy lungs were defined as: no active or past history of smoking, no previous intrathoracic procedure, no known pulmonary disease, no oxygen requirement and metabolic equivalent task (METS) greater than or equal to 4. To provide optimal imaging conditions for this pilot study, obese patients (body mass index > 30 kg.m–2) were excluded.
Interventions
All patients were pre-oxygenated in the supine position with 100% oxygen for 3 min without any continuous positive airway pressure. While specific doses were not protocolized, general anesthesia induction was performed using standard doses of propofol and fentanyl. Rocuronium was used in all cases to facilitate tracheal intubation. Anesthesia was maintained with desflurane or sevoflurane. All patients were ventilated with Datex-Ohmeda Aestiva 3000 machines (GE Healthcare, WI, United States) using volume-controlled ventilation, an inspired oxygen fraction (FiO2) of 40–50%, a respiratory rate of 12 min–1, an inspiratory to expiratory ratio of 1:2, no inspiratory pause and a positive end-expiratory pressure of 6 cm H2O.
After anesthesia induction, patients were administered tidal volumes of 6, 8, 10, and 12 mL.kg–1 predicted body weight (20) in random order (Supplementary Figure 1). For each tidal volume, mean expired tidal volume of 3 consecutive breaths were collected and the pleura was imaged at 4 predetermined anatomical locations: Left and right 3rd intercostal space at the mid-clavicular line, and left and right 8th intercostal space at the posterior axillary line. The correct intercostal spaces were identified by sliding the ultrasound transducer from the clavicle downwards and visually counting rib spaces. For each tidal volume and anatomical location tested, 3 ultrasound radiofrequency cineloops at a frame rate of 30 Hz were acquired over 3 separate respiratory cycles. All cineloops were saved to digital format for offline analysis. Interobserver reliability was assessed by repeating cineloop acquisition by a second blinded observer (SG) for the 10 mL.kg–1 predicted body weight tidal volume. Immediately after, test-retest reliability was also assessed by repeating cineloop acquisition by the first observer (MG) for the 10 mL.kg–1 predicted body weight tidal volume.
Lung ultrasonography
Lung ultrasonography was performed by experienced lung echographists (MG and SG for repeated measures with 8 and 1 years of experience) using a Terason T3000cv scanner (Teratech Corporation, Burlington, MA) and a 12 MHz transducer (probe #12L5). Initial depth of field was 4 cm and adjusted as needed to position the pleura between half to three-quarters of the screen. A single focal zone was placed nearest to the pleura. With the marker pointing toward the head, the probe was oriented perpendicular to the ribs with the pleura as horizontal as possible.
Elastography
B-mode images were reconstructed from radiofrequency data (Supplementary Figure 2A). For each cineloop, the pleura was segmented manually on a single frame (Supplementary Figure 2B). With the segmented pleura forming the upper boundary, a region of interest (ROI) of a fixed depth of 2 mm was defined (Supplementary Figure 2C). The geometry of the ROI was automatically adapted and tracked throughout the respiratory cycle (21) or, if inadequate, simply copied over from frame to frame. Tracking was considered adequate when the pleural line remained within the ROI at all-time as assessed visually by an experienced lung echographist (MG). Elastography parameters were computed on all cineloops with a visible pleura and an adequately tracked ROI. We used the Lagrangian speckle model estimator to compute tissue translation, strain and shear values (Supplementary Figure 3) (22). We defined translations as rigid displacements produced by lung sliding, strain as expansion or contraction of the pleura from tidal volume administration and shear as the angular deformation.
Elastography images and mechanical parameters were computed within the ROI over consecutive frames. We used an implementation of the estimator integrated into a commercial imaging platform (Visual, Object Research Systems, Montréal, Canada) (23). This implementation computes the various elastography parameters using the raw radiofrequency data and not on the reconstructed B-mode images. As such, computed elastography parameters are not affected by gain or image post-processing techniques. Axial and lateral elastography components were determined (“axial” indicating the direction along the US beam and “lateral” indicating the direction perpendicular to it). Considering the planar nature of the pleura and its perpendicular orientation with respect to the US beam, we restricted our analysis to lateral translation, strain and shear components along with the Von Mises strain, a combination of bidimensional strain and shear components. Six elastography parameters were computed per cineloop of a given patient (Figure 1 and Table 1).
Figure 1. Calculating lateral strain and lateral absolute strain values. (A) Instantaneous strain values are computed in all sub-ROIs between consecutive frames of a cineloop. (B) By averaging all instantaneous sub-ROI strain values in a single frame, instantaneous strain values for the whole ROI are plotted for all frames of the cineloop. (C) The summation of instantaneous strain values produces the cumulative strain of the pleura. Lateral strain is the range of the cumulative lateral strain experienced by the lung in the ROI. (D) On the other hand, by averaging all absolute sub-ROI instantaneous strain values in a single frame, instantaneous absolute strain values for the whole ROI are plotted for all frames of the cineloop. (E) The summation of the instantaneous absolute strain values produces the cumulative absolute strain of the pleura. Lateral absolute strain is the range of the cumulative lateral absolute strain experienced by the lung in the ROI.
Table 1. Description of elastography parameters.
Randomization and blinding
Simple randomization was performed by a statistician not otherwise implicated in this study using a computer random number generator [R v3.4.0, R Core Team, (24); blockrand package v1.3, Snow, (25)]. Concealment was ensured by the use of sequentially numbered sealed and opaque envelopes. Both observers were not blinded to tidal volume during lung cineloop acquisition and during image analysis.
Objectives
The primary objective was the feasibility of pleural strain measurement using ultrasound elastography. Secondary objectives were: (1) Verifying dose-response with varying tidal volumes and model fit for elastography parameters; (2) measuring reliability (intraobserver, interobserver and test-retest) of elastography parameters.
Outcomes
The primary outcome was feasibility defined as the percentage of all cineloops mandated by the protocol on which we computed elastography parameters. Secondary outcomes were: (1) Estimated elastography parameter slopes for tidal volume and marginal and conditional coefficients of determination (R2); (2) intraclass correlation coefficients (ICC) for intraobserver, interobserver and test-retest reliability values of elastography parameters.
Statistical analysis
We enrolled a convenience sample of 10 patients. For the primary outcome, with 72 potential cineloops per patient, our margin of error is 2.2% for a 95% confidence level and an expected 90% feasibility.
For the first secondary outcomes, computation results of all 6 elastography parameters were modeled. All continuous dependent and independent variables were centered and reduced. This allowed direct comparison of slope estimates from the various models. To simplify models, all repeated measurements performed at a tidal volume of 10 mL.kg–1 used to measure interobserver and test-retest reliability values were not included in the models. Linear mixed-effect models were used with elastography parameters as dependent variables and tidal volume, side of measurement (left/right) and gravity dependence of measurement (dependent/non-dependent) as independent variables. Because of the repeated (and thus correlated) nature of measurements per patient and per anatomical location, a random intercept per patient was included in the model as well as a one per anatomical location as a nested grouping factor, whereas triplicate measurements used to measure intraobserver reliability were averaged. Interaction between tidal volume and gravity dependence was included in the model (9). Model assumptions were verified. All six slope estimates were tested for significance. Using Bonferroni’s adjustment, a p-value of 0.008 (0.05/6) was considered significant. All other analyses are considered exploratory. To identify elastography parameters with the best dose-response with tidal volumes, absolute values of estimated slopes for significant parameters were ordered. Parameters with highest absolute values of estimated slopes with non-overlapping 95% confidence intervals were selected. Goodness of fit was assessed by the marginal and conditional R2 (26). To ensure robustness of results, a sensitivity analysis was performed by using the same linear mixed-effect models described above but using a smaller subset of cineloops with better imaging quality. Cineloops with better imaging quality were defined as having a thin and clearly defined pleural line and perfect ROI tracking.
For the second secondary outcome, intraobserver, interobserver and test-retest reliability measures were calculated using ICC type (2,1) (27). As outlined above, cineloops were acquired by different observers while the analysis process was performed by the same one. Bootstrap was performed as a second sensitivity analysis to calculate 95% confidence intervals using 10,000 iterations and the bias-corrected and accelerated method (28). As suggested, ICC values less than 0.5 indicate poor reliability, values between 0.5 and 0.75 indicate moderate reliability, values between 0.75 and 0.9 indicate good reliability and values greater than 0.9 indicate excellent reliability (27). Bland-Altman plots were also performed for intraobserver, interobserver and test-retest reliability. We computed mean bias and 95% limits of agreements for interobserver and test-retest reliability measures (29) but not for intraobserver reliability measure as we could not identify a generally accepted methodology that takes into account triplicate measurements for a single observer.
Results are expressed as mean ± standard deviation or median and interquartile range (25–75%) as appropriate. No imputation for missing values was performed. Statistical analyses were performed using R [v3.4.0, R Core Team, (24)].
Results
Twelve patients were assessed for eligibility. Two patients were excluded post-enrollment (Supplementary Figure 4): one change in anesthetic plan and one consent withdrawal. Patients’ baseline characteristics are summarized in Table 2.
Table 2. Patient characteristics.
With our protocol specifying that 72 cineloops would be acquired for each of the 10 patients randomized, 720 elastograms were planned to be computed upon completion of cineloop acquisition. Sixty-two cineloops (8.6%) weren’t subjected to segmentation and tracking: 4 cineloops were not acquired because of a manipulation error, 7 cineloops were incomplete because of a technical issue with the echograph, and 51 cineloops didn’t show sufficient pleura throughout the respiratory cycle. This last problem was only encountered for cineloops acquired over the left 3rd intercostal space at the mid-clavicular line where the heart was mostly seen. Lastly, elastograms were not computed for 6 cineloops because of inadequate tracking of the ROI throughout the respiratory cycle. Elastogram computation was thus performed on 652 (90.6%) cineloops.
Estimates for fixed effects of all 6 models can be found in Table 3. We observed a significant linear increase in 4 elastography parameters with increasing tidal volume: lateral absolute shear, lateral absolute strain, Von Mises strain and lateral absolute translation (Figure 2). With overlapping 95% confidence intervals, no significant elastography parameter was superior to the other in its ability to capture an increase in tidal volume (Figure 3). Amongst them, lateral absolute shear, lateral absolute strain and Von Mises strain were more precisely predicted by tidal volume compared to lateral absolute translation, as shown by higher goodness of fit (Table 3). The first sensitivity analysis excluding suboptimal images yielded similar results (Supplementary Table 1 and Supplementary Figure 5). When side and gravity dependence of measurements were considered, only gravity dependence seemed to be associated with our results (Table 3). Dependent lung zones had initially lower values for all significant strain elastography parameters, but pleural strain measurements in dependent lung zones increased more rapidly with increasing tidal volume although this was only significantly with the Von Mises strain parameter (Table 3 and Figures 2, 3).
Table 3. Modeled elastography parameters.
Figure 2. Regression lines, 95% confidence bands and individual data points for all elastography parameters across the various tidal volumes stratified by gravity dependence. (A) Lateral translation. (B) Lateral absolute translation. (C) Lateral strain. (D) Lateral absolute strain. (E) Lateral absolute shear. (F) Von Mises strain. PBW, predicted body weight.
Figure 3. Slope estimates for elastography parameters in increasing order stratified by gravity dependence. Significant parameters are identified by an asterix.
Measured intraobserver reliability was excellent for lateral absolute shear, lateral absolute strain and Von Mises strain. Interobserver and test-retest measured reliabilities were found to be moderate to good from Von Mises strain and moderate to excellent for lateral absolute shear and lateral absolute strain (Figure 4 and Supplementary Table 2). In our sensitivity analysis, bootstrapping yielded identical estimates (Supplementary Table 3). Bland-Altman plots, mean bias and 95% limits of agreements can be found in Supplementary Table 4 and in Supplementary Figure 6).
Figure 4. Intraclass correlation coefficients and 95% confidence intervals for intraobserver, interobserver and test-retest reliability measures for elastography parameters. Intraclass correlation coefficients in the red panel indicate poor reliability, values in the orange panel indicate moderate reliability, values in the yellow panel indicate good reliability and values in the green panel indicate excellent reliability.
Discussion
We showed preliminary evidence that computing regional pleural translation, strain and shear components is feasible in over 90% of cineloops. Lateral absolute shear, lateral absolute strain and Von Mises strain were more informative in regards to administered tidal volumes than lateral translation, lateral absolute translation and lateral strain. Measured intraobserver reliability was excellent for all significant parameters while interobserver and test-retest measured reliability were found to be moderate to excellent. This supports our hypothesis that lung ultrasound could be used to measure regional pleural strain in mechanically ventilated patients.
While others have previously described ultrasonographically-measured regional pleural strain in humans (17, 18), this study significantly extends previous work by quantifying the impact of tidal volume on regional pleural strain. Despite the absence of a gold standard precluding any firm conclusion about the utility of each elastography parameter, we described an expected dose-response relationship between various strain measurements and administered tidal volumes.
At first glance, it is somewhat surprising that a 2 mm thick ROI delineated from an image generated by a 12 MHz probe could capture physiological signal from the pleura, a structure only 30–40 microns thick (30). With the topmost part of the ROI aligned with the pleura, the image contained within the ROI is actually mainly composed of reverberation artifacts originating from the pleura and other soft tissue/ultrasound beam interfaces (31). These artifacts are most easily seen when contrasting normal lung sliding in M-mode and their disappearance when a pneumothorax occurs (Supplementary Figure 7). We hypothesize that reverberation artifacts in the ROI convey anatomical information on the pleura thus allowing our algorithm to compute elastography parameters. Another possibility would be that imperfect ROI tracking leads to intercoastal muscle being included in the strain calculation. While it is technically possible to measure intercostal muscle strain (32), we do not believe this is the case with our results as patients were paralyzed for induction of general anesthesia.
One important finding is that lateral strain had a poor dose-response with tidal volume despite a time plot of lateral strain that seemed to track delivered tidal volume in many patients (Figure 1C). Perhaps explaining the poor performance of lateral strain compared to lateral absolute strain, we observed small areas of negative strain during lung inflation (Figure 1A). Supporting this hypothesis, zones of local deflation during inspiration have been described in animals (Supplementary Figure 8) (33, 34). Analogous to how lateral strain and lateral absolute strain are calculated, global pulmonary compliance was correlated with the sum of positive values of local compliance, measured by computed tomography, but not with the sum of positive and negative values of local compliance (34). While the exact causes for this phenomenon are still debated, tidal recruitment with gas redistribution may contribute. Discrepancies between lateral strain and lateral absolute strain could signal tidal recruitment and ongoing atelectrauma.
Our study has some important strengths and limitations. First, our results are physiologically plausible with lateral absolute shear, lateral absolute strain and Von Mises strain increasing with increasing tidal volumes. Second, our two sensitivity analyses demonstrate the robustness of results. Third, by performing a randomized crossover study, we ruled out the possibility of an order effect. Fourth, we demonstrated the good to excellent reliability of these measurements. Given the proof-of-concept nature of this study, several limitations should be considered. First, the limited sample size restricts our ability to obtain precise estimates of pleural strain in human. Second, no gold standard for pulmonary strain was available for comparison. While helium dilution or nitrogen-washout techniques could have been used to measure end-expiratory lung volume at each step of our protocol, these methods cannot account for tidal recruitment or overdistension (35). Third, whether regional pleural strain is a good surrogate for regional pulmonary strain remains to be determined. Alveolar deflation kinetics in a small animal model suggest that differences between the lung’s periphery and core regions exist and that these may be altered by disease states (36, 37). On the other hand, lung lesions in another VILI animal model developed predominantly in the subpleural regions (38). Further work will need to be done to determine the clinical usefulness of lung US-measured regional pleural strain. Fourth, flow, pressure and volume curves were not recorded simultaneously with each image acquisition. As such, the precise beginning and end of each respiratory cycle wasn’t exactly known which may have affected our results. The exact tidal volume administered while acquiring cineloops may also have been slightly different from the tidal volume recorded beforehand. Fifth, physiological lung sliding and rib shadows preclude measuring the same pleural surface throughout the respiratory cycle. While we computed strain between two consecutive frames of a given cineloop, the pleura imaged in the ROI at end-inspiration was likely different from the pleura in the ROI at end-expiration. As such, strain measures likely represent average strain values from a pulmonary neighborhood. Sixth, the elastography parameters calculated from the Lagrangian speckle model estimator reflect dynamic pleural strain (pleural deformation from regional tidal volume/gas volume at set positive end-expiratory pressure) and not total pleural strain (pleural deformation from total gas volume/gas volume at functional residual capacity). While experimental evidence suggests dynamic strain is the most important strain component leading to the development of VILI (39–41), the clinical usefulness of lung US-measured regional pleural strain will need to be verified in future studies. Seventh, patients included in this proof-of-concept study were highly selected as we excluded patients with any lung pathology and obese patients. While our first sensitivity analysis that excluded suboptimal images did not show different results (Supplementary Table 1 and Supplementary Figure 5), our results’ external validity remains to be demonstrated and future studies will have to include obese patients and patients with various lung pathologies. Finally, like all ultrasonographic exams, lung ultrasonography is operator dependent. Fortunately, because of lung ultrasonography’s shallow learning curve, it is possible to train a clinician rapidly. As few as 10 supervised exams have been found to be necessary when assessing simple pathologies (42) although interpreting more complicated quantitative aspects of lung US requires a 2-month period and 25 supervised exams (43).
In conclusion, measuring regional pleural strain by ultrasound elastography is feasible. A significant dose-response relationship between tidal volume and lateral absolute shear, lateral absolute strain and Von Mises strain was observed, supporting the hypothesis that regional pleural strain can be measured at the bedside with non-invasive lung ultrasound. Further work will be required to compare elastography parameters to gold-standard strain measurements and how they can be used to tailor mechanical ventilation in the operating room and the intensive care unit to improve patient outcomes.
Data availability statement
Due to national regulations in the Province of Quebec (Canada), health medical data cannot be made available publicly. However, access to the research dataset is possible for research purpose after appropriate privacy agreements between research parties have been completed. Data access requests may be sent to the corresponding author MG, martin.girard@umontreal.ca, or directly to the CHUM REB (ethique.recherche.chum@ssss.gouv.qc.ca). The R code will be available upon request to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Comité d’éthique de la recherche du CHUM. The patients/participants provided their written informed consent to participate in this study.
Author contributions
MG, YC, GC, and AD designed the study. MG and SG collected the data. MG, M-HRC, MC, GC, and AD analyzed the data. MG wrote the manuscript. All authors have read and approved the final manuscript.
Funding
Funding for the present work was secured through internal departmental funds.
Acknowledgments
We thank Monique Ruel, RN, for her valuable assistance and dedication and Brian Grondin-Beaudoin, MD, for his help with preliminary testing procedures.
Conflict of Interest
MG is a paid consultant for the point-of-care ultrasonography group of GE Healthcare. GC has an active commercial license with Rhéolution Inc. (Montréal, Canada) and a license option with Siemens Healthcare. AD reported non-financial educational material support from CAE Healthcare, research equipment grants from Edwards and was on Masimo’s speaker bureau.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2022.935482/full#supplementary-material
References
1. Slutsky AS, Ranieri VM. Ventilator-induced lung injury. N Engl J Med. (2013) 369:2126–36. doi: 10.1056/NEJMra1208707
2. Del Sorbo L, Goligher EC, McAuley DF, Rubenfeld GD, Brochard LJ, Gattinoni L, et al. Mechanical ventilation in adults with acute respiratory distress syndrome. Summary of the experimental evidence for the clinical practice guideline. Ann Am Thorac Soc. (2017) 14:S261–70.s doi: 10.1513/AnnalsATS.201704-345OT
3. Lohser J, Slinger P. Lung injury after one-lung ventilation: A review of the pathophysiologic mechanisms affecting the ventilated and the collapsed lung. Anesth Analg. (2015) 121:302–18. doi: 10.1213/ANE.0000000000000808
4. Modesto I Alapont V, Aguar Carrascosa M. Medina Villanueva A. Stress, strain and mechanical power: Is material science the answer to prevent ventilator induced lung injury? Med Intensiva. (2019) 43:165–75. doi: 10.1016/j.medin.2018.06.008
5. Gattinoni L, Marini JJ, Pesenti A, Quintel M, Mancebo J, Brochard L. The “baby lung” became an adult. Intensive Care Med. (2016) 42:663–73. doi: 10.1007/s00134-015-4200-8
6. Marini JJ, Rocco PRM, Gattinoni L. Static and dynamic contributors to ventilator-induced lung injury in clinical practice. pressure, energy, and power. Am J Respir Crit Care Med. (2020) 201:767–74. doi: 10.1164/rccm.201908-1545CI
7. Blankman P, Hasan D, Bikker IG, Gommers D. Lung stress and strain calculations in mechanically ventilated patients in the intensive care unit. Acta Anaesth Scand. (2016) 60:69–78. doi: 10.1111/aas.12589
8. Wellman TJ, Winkler T, Costa ELV, Musch G, Harris RS, Zheng H, et al. Effect of local tidal lung strain on inflammation in normal and lipopolysaccharide-exposed sheep*. Crit Care Med. (2014) 42:e491–500. doi: 10.1097/CCM.0000000000000346
9. Paula LF, Wellman TJ, Winkler T, Spieth PM, Güldner A, Venegas JG, et al. Regional tidal lung strain in mechanically ventilated normal lungs. J Appl Physiol. (2016) 121:1335–47. doi: 10.1152/japplphysiol.00861.2015
10. Tustison NJ, Awate SP, Cai J, Altes TA, Miller GW, de Lange EE, et al. Pulmonary kinematics from tagged hyperpolarized helium-3 MRI. J Magn Reson Imaging. (2010) 31:1236–41. doi: 10.1002/jmri.22137
11. Retamal J, Hurtado D, Villarroel N, Bruhn A, Bugedo G, Amato MBP, et al. Does regional lung strain correlate with regional inflammation in acute respiratory distress syndrome during nonprotective ventilation? An experimental porcine study. Crit Care Med. (2018) 46:e591–9. doi: 10.1097/CCM.0000000000003072
12. Schwebel C, Clec’h C, Magne S, Minet C, Garrouste-Orgeas M, Bonadona A, et al. Safety of intrahospital transport in ventilated critically Ill patients: A multicenter cohort study*. Crit Care Med. (2013) 41:1919–28. doi: 10.1097/CCM.0b013e31828a3bbd
13. Krishnan S, Moghekar A, Duggal A, Yella J, Narechania S, Ramachandran V, et al. Radiation exposure in the medical ICU: Predictors and characteristics. Chest. (2018) 153:1160–8. doi: 10.1016/j.chest.2018.01.019
14. Pereira SM, Tucci MR, Morais CCA, Simões CM, Tonelotto BFF, Pompeo MS, et al. Individual positive end-expiratory pressure settings optimize intraoperative mechanical ventilation and reduce postoperative atelectasis. Anesthesiology. (2018) 129:1070–81. doi: 10.1097/ALN.0000000000002435
15. Mongodi S, De Luca D, Colombo A, Stella A, Santangelo E, Corradi F, et al. Quantitative lung ultrasound: Technical aspects and clinical applications. Anesthesiology. (2021) 134:949–65. doi: 10.1097/ALN.0000000000003757
16. Généreux V, Chassé M, Girard F, Massicotte N, Chartrand-Lefebvre C, Girard M. Effects of positive end-expiratory pressure/recruitment manoeuvres compared with zero end-expiratory pressure on atelectasis during open gynaecological surgery as assessed by ultrasonography: A randomised controlled trial. Br J Anaesth. (2020) 124:101–9. doi: 10.1016/j.bja.2019.09.040
17. Rubin JM, Horowitz JC, Sisson TH, Kim K, Ortiz LA, Hamilton JD. Ultrasound strain measurements for evaluating local pulmonary ventilation. Ultrasound Med Biol. (2016) 42:2525–31. doi: 10.1016/j.ultrasmedbio.2016.05.020
18. Duclos G, Bobbia X, Markarian T, Muller L, Cheyssac C, Castillon S, et al. Speckle tracking quantification of lung sliding for the diagnosis of pneumothorax: A multicentric observational study. Intensive Care Med. (2019) 45:1 212–8. doi: 10.1007/s00134-019-05710-1
19. Dwan K, Li T, Altman DG, Elbourne D. CONSORT 2010 statement: Extension to randomised crossover trials. BMJ Br Med J Publ Group. (2019) 366:l4378. doi: 10.1136/bmj.l4378
20. Acute Respiratory Distress Syndrome Network, Brower RG, Matthay MA, Morris A, Schoenfeld D, Thompson BT, et al. Ventilation with lower tidal volumes as compared with traditional tidal volumes for acute lung injury and the acute respiratory distress syndrome. N Engl J Med. (2000) 342:1301–8. doi: 10.1056/NEJM200005043421801
21. Destrempes F, Meunier J, Giroux M-F, Soulez G, Cloutier G. Segmentation of plaques in sequences of ultrasonic B-mode images of carotid arteries based on motion estimation and a Bayesian model. IEEE Trans Biomed Eng. (2011) 58:2202–11. doi: 10.1109/TBME.2011.2127476
22. Maurice RL, Daronat M, Ohayon J, Stoyanova E, Foster FS, Cloutier G. Non-invasive high-frequency vascular ultrasound elastography. Phys Med Biol. (2005) 50:1611–28. doi: 10.1088/0031-9155/50/7/020
23. Roy Cardinal M-H, Heusinkveld MHG, Qin Z, Lopata RGP, Naim C, Soulez G, et al. Carotid artery plaque vulnerability assessment using noninvasive ultrasound elastography: Validation with MRI. AJR Am J Roentgenol. (2017) 209:142–51. doi: 10.2214/AJR.16.17176
24. R Core Team. R: A Language And Environment For Statistical Computing. Vienna: R Foundation for Statistical Computing (2017).
25. Snow G. Blockrand: Randomization for Block Random Clinical Trials_. R Package Version 1.3. (2013). Available online at: https://CRAN.R-project.org/package=blockrand
26. Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. O’Hara RB, editor. Methods Ecol Evol. (2013) 4:133–42.
27. Koo TK, Li MYA. Guideline of selecting and reporting intraclass correlation coefficients for reliability research. J Chiropr Med. (2016) 15:155–63. doi: 10.1016/j.jcm.2016.02.012
29. Zou GY. Confidence interval estimation for the Bland-Altman limsits of agreement with multiple observations per individual. Stat Methods Med Res. (2013) 22:630–42. doi: 10.1177/0962280211402548
30. Lai-Fook SJ. Pleural mechanics and fluid exchange. Physiol Rev. (2004) 84:385–410. doi: 10.1152/physrev.00026.2003
31. Soldati G, Demi M, Smargiassi A, Inchingolo R, Demi L. The role of ultrasound lung artifacts in the diagnosis of respiratory diseases. Expert Rev Respir Med. (2019) 13:163–72. doi: 10.1080/17476348.2019.1565997
32. Rubin JM, Horowitz JC, Sisson TH, Kim K, Ortiz LA, Hamilton JD Ultrasound strain measurements for evaluating local pulmonary ventilation. Proceedings of the 2015 IEEE International Ultrasonics Symposium (IUS), Taipei. (2015). doi: 10.1109/ULTSYM.2015.0181
33. Schiller HJ, Steinberg J, Halter J, McCann U, DaSilva M, Gatto LA, et al. Alveolar inflation during generation of a quasi-static pressure/volume curve in the acutely injured lung. Crit Care Med. (2003) 31:1126–33.s doi: 10.1097/01.CCM.0000059997.90832.29
34. Perchiazzi G, Rylander C, Derosa S, Pellegrini M, Pitagora L, Polieri D, et al. Regional distribution of lung compliance by image analysis of computed tomograms. Respir Physiol Neurobiol. (2014) 201:60–70. doi: 10.1016/j.resp.2014.07.001
35. Chiumello D, Marino A, Brioni M, Cigada I, Menga F, Colombo A, et al. Lung recruitment assessed by respiratory mechanics and computed tomography in patients with acute respiratory distress syndrome. What is the relationship? Am J Respir Crit Care Med. (2016) 193:1254–63. doi: 10.1164/rccm.201507-1413OC
36. Scaramuzzo G, Broche L, Pellegrini M, Porra L, Derosa S, Tannoia AP, et al. The effect of positive end-expiratory pressure on lung micromechanics assessed by synchrotron radiation computed tomography in an animal model of ARDS. J Clin Med. (2019) 8:1117. doi: 10.3390/jcm8081117
37. Scaramuzzo G, Broche L, Pellegrini M, Porra L, Derosa S, Tannoia AP, et al. Regional behavior of airspaces during positive pressure reduction assessed by synchrotron radiation computed tomography. Front Physiol. (2019) 10:719. doi: 10.3389/fphys.2019.00719
38. Cressoni M, Chiurazzi C, Gotti M, Amini M, Brioni M, Algieri I, et al. Lung inhomogeneities and time course of ventilator-induced mechanical injuries. Anesthesiology. (2015) 123:618–27. doi: 10.1097/ALN.0000000000000727
39. Protti A, Andreis DT, Monti M, Santini A, Sparacino CC, Langer T, et al. Lung stress and strain during mechanical ventilation: Any difference between statics and dynamics?*. Crit Care Med. (2013) 41:1046–55. doi: 10.1097/CCM.0b013e31827417a6
40. Tschumperlin DJ, Oswari J, Margulies AS. Deformation-induced injury of alveolar epithelial cells. Effect of frequency, duration, and amplitude. Am J Respir Crit Care Med. (2000) 162:357–62. doi: 10.1164/ajrccm.162.2.9807003
41. Yen S, Preissner M, Bennett E, Dubsky S, Carnibella R, O’Toole R, et al. The link between regional tidal stretch and lung injury during mechanical ventilation. Am J Respir Cell Mol Biol. (2019) 60:569–77. doi: 10.1165/rcmb.2018-0143OC
42. See KC, Ong V, Wong SH, Leanda R, Santos J, Taculod J, et al. Lung ultrasound training: Curriculum implementation and learning trajectory among respiratory therapists. Intensive Care Med. (2016) 42:63–71. doi: 10.1007/s00134-015-4102-9
Keywords: mechanical ventilalion, ventilator-induced lung injury, general anesthesia, lung imaging, pulmonary strain, ultrasound elastography
Citation: Girard M, Roy Cardinal M-H, Chassé M, Garneau S, Cavayas YA, Cloutier G and Denault AY (2022) Regional pleural strain measurements during mechanical ventilation using ultrasound elastography: A randomized, crossover, proof of concept physiologic study. Front. Med. 9:935482. doi: 10.3389/fmed.2022.935482
Received: 04 May 2022; Accepted: 22 August 2022;
Published: 15 September 2022.
Edited by:
Paola Pierucci, Azienda Ospedaliero Universitaria Consorziale Policlinico di Bari, ItalyReviewed by:
Penny Andrews, R Adams Cowley Shock Trauma Center, United StatesEwan Goligher, University of Toronto, Canada
Gaetano Scaramuzzo, University of Ferrara, Italy
Copyright © 2022 Girard, Roy Cardinal, Chassé, Garneau, Cavayas, Cloutier and Denault. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Martin Girard, martin.girard@umontreal.ca