ORIGINAL RESEARCH article
Impact of Different Tidal Volume Levels at Low Mechanical Power on Ventilator-Induced Lung Injury in Rats
- 1Laboratory of Pulmonary Investigation, Carlos Chagas Filho Institute of Biophysics, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil
- 2Laboratory of Cellular and Molecular Physiology, Carlos Chagas Filho Institute of Biophysics, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil
- 3Department of Pathology, School of Medicine, University of São Paulo, São Paulo, Brazil
- 4Pulmonary Engineering Group, Department of Anesthesiology and Intensive Care Therapy, University Hospital Carl Gustav Carus, Dresden University of Technology, Dresden, Germany
- 5Department of Surgical Sciences and Integrated Diagnostics, San Martino Policlinico Hospital, IRCCS for Oncology, University of Genoa, Genoa, Italy
Tidal volume (VT) has been considered the main determinant of ventilator-induced lung injury (VILI). Recently, experimental studies have suggested that mechanical power transferred from the ventilator to the lungs is the promoter of VILI. We hypothesized that, as long as mechanical power is kept below a safe threshold, high VT should not be injurious. The present study aimed to investigate the impact of different VT levels and respiratory rates (RR) on lung function, diffuse alveolar damage (DAD), alveolar ultrastructure, and expression of genes related to inflammation [interleukin (IL)-6], alveolar stretch (amphiregulin), epithelial [club cell secretory protein (CC)16] and endothelial [intercellular adhesion molecule (ICAM)-1] cell injury, and extracellular matrix damage [syndecan-1, decorin, and metalloproteinase (MMP)-9] in experimental acute respiratory distress syndrome (ARDS) under low-power mechanical ventilation. Twenty-eight Wistar rats received Escherichia coli lipopolysaccharide intratracheally. After 24 h, 21 animals were randomly assigned to ventilation (2 h) with low mechanical power at three different VT levels (n = 7/group): (1) VT = 6 mL/kg and RR adjusted to normocapnia; (2) VT = 13 mL/kg; and 3) VT = 22 mL/kg. In the second and third groups, RR was adjusted to yield low mechanical power comparable to that of the first group. Mechanical power was calculated as [(ΔP,L2/Est,L)/2]× RR (ΔP,L = transpulmonary driving pressure, Est,L = static lung elastance). Seven rats were not mechanically ventilated (NV) and were used for molecular biology analysis. Mechanical power was comparable among groups, while VT gradually increased. ΔP,L and mechanical energy were higher in VT = 22 mL/kg than VT = 6 mL/kg and VT = 13 mL/kg (p < 0.001 for both). Accordingly, DAD score increased in VT = 22 mL/kg compared to VT = 6 mL/kg and VT = 13 mL/kg [23(18.5–24.75) vs. 16(12–17.75) and 16(13.25–18), p < 0.05, respectively]. VT = 22 mL/kg was associated with higher IL-6, amphiregulin, CC16, MMP-9, and syndecan-1 mRNA expression and lower decorin expression than VT = 6 mL/kg. Multiple linear regression analyses indicated that VT was able to predict changes in IL-6 and CC16, whereas ΔP,L predicted pHa, oxygenation, amphiregulin, and syndecan-1 expression. In the model of ARDS used herein, even at low mechanical power, high VT resulted in VILI. VT control seems to be more important than RR control to mitigate VILI.
Mechanical ventilation with high tidal volume (VT) can promote ventilator-induced lung injury (VILI) (Tremblay and Slutsky, 2006). In patients with the acute respiratory distress syndrome (ARDS), VT level has been recognized as a major risk factor for organ failure and death (Brower et al., 2000; Putensen et al., 2009), and the use of low VT (4–8 mL/kg) represents the cornerstone of protective mechanical ventilation (Bellani et al., 2016; Fan et al., 2017). Additionally, the respiratory system driving pressure (ΔP,RS) is associated with increased lung inflammation (Bellani et al., 2011) and predicts mortality rate in ARDS (Amato et al., 2015). Respiratory rate (RR) (Rich et al., 2003) as well as inspiratory (Fujita et al., 2007) and expiratory (Schumann et al., 2014) airflow have also been implicated in VILI. The knowledge that different respiratory variables may be injurious has led to the concept that mechanical energy or power transferred from the ventilator to the lungs may be determinants of VILI pathogenesis (Cressoni et al., 2016). According to this concept, we hypothesized that, as long as mechanical power is kept below a safe threshold, high VT levels should not be injurious to the lungs.
Thus, the present study aimed to investigate the impact of different VT levels and RR on lung function, diffuse alveolar damage (DAD), alveolar ultrastructure, and expression of genes related to inflammation [interleukin (IL)-6], alveolar stretch (amphiregulin), epithelial [club cell secretory protein (CC)16] and endothelial [intercellular adhesion molecule (ICAM)-1] cell injury, and extracellular matrix damage [syndecan-1, decorin, and metalloproteinase (MMP)-9] in experimental ARDS under low-power mechanical ventilation.
Materials and Methods
This study was approved by the Animal Care and Use Committee (CEUA: 064-17) of the Health Sciences Center, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil (chair: Prof. M. Frajblat). All animals received humane care in compliance with the “Principles of Laboratory Animal Care” formulated by the National Society for Medical Research and the U.S. National Academy of Sciences Guide for the Care and Use of Laboratory Animals. The present study followed the ARRIVE guidelines for reporting of animal research (Kilkenny et al., 2010).
Animal Preparation and Experimental Protocol
Twenty-eight Wistar rats (weight 383 ± 80 g) were anaesthetized by inhalation of sevoflurane 2% (Sevorane®; Cristália, Itapira, SP, Brazil), and intratracheal (i.t.) instillation of Escherichia coli lipopolysaccharide (O55:B5, LPS Ultrapure, Invivogen), 200 μg suspended in saline solution to a total volume of 200 μL (Samary et al., 2016a) to induce acute lung inflammation. After 24 h, animals were premedicated intraperitoneally (i.p.) with 10 mg/kg diazepam (Compaz, Cristália, Itapira, SP, Brazil), followed by 100 mg/kg ketamine (Ketamin-S+, Cristália, Itapira, SP, Brazil) and 2 mg/kg midazolam (Dormicum, União Química, São Paulo, SP, Brazil). After local anesthesia with 2% lidocaine (0.4 mL), a midline neck incision and tracheostomy were performed. Seven rats were used for molecular biology analysis and were not mechanically ventilated (nonventilated, NV).
An intravenous (i.v.) catheter (Jelco 24G, Becton, Dickinson and Company, New Jersey, NJ, USA) was inserted into the tail vein, and anesthesia induced and maintained with midazolam (2 mg/kg/h) and ketamine (50 mg/kg/h). Additionally, 10 mL/kg/h Ringer's lactate (B. Braun, Crissier, Switzerland) was administered i.v. Gelafundin® 4% (B. Braun, São Gonçalo, RJ, Brazil) was administered intravenously (in 0.5-mL increments) to maintain mean arterial pressure (MAP) >60 mmHg. A second catheter (18G, Arrow International, USA) was then placed in the right internal carotid artery for blood sampling and gas analysis (Radiometer ABL80 FLEX, Copenhagen NV, Denmark), as well as monitoring of MAP (Networked Multiparameter Veterinary Monitor LifeWindow 6,000 V; Digicare Animal Health, Boynton Beach, FL, USA). A 30-cm-long water-filled catheter (PE-205, Becton, Dickinson and Company) with side holes at the tip, connected to a differential pressure transducer (UT-PL-400, SCIREQ, Montreal, QC, Canada), was used to measure the esophageal pressure (Pes). The catheter was passed into the stomach and then slowly returned into the esophagus; its proper positioning was assessed with the “occlusion test” (Baydur et al., 1982). Heart rate (HR), MAP, and rectal temperature were continuously monitored (Networked Multiparameter Veterinary Monitor LifeWindow 6,000V, Digicare Animal Health, Florida, USA). Body temperature was maintained at 37.5 ± 1°C using a heating bed.
Animals were paralyzed with pancuronium bromide (2 mg/kg, i.v.), and their lungs mechanically ventilated (Servo-i, MAQUET, Solna, Sweden) in volume-controlled mode (VCV) with constant inspiratory airflow, VT = 6 mL/kg, RR to maintain normocapnia (PaCO2 = 35–45 mmHg) (around 70 bpm), positive end-expiratory pressure (PEEP) = 3 cmH2O, FiO2 = 1.0, and an inspiratory-expiratory ratio of 1:2 (Figure 1). Arterial blood gases (300 μl) were determined using a Radiometer ABL80 FLEX analyzer (Copenhagen NV, Denmark) and respiratory system mechanics were assessed (BASELINE). FiO2 was reduced to 0.4 to prevent possible iatrogenic effects. Rats were then assigned to three different tidal volumes with low mechanical power (n = 7/group): (1) VT = 6 mL/kg and RR adjusted to normocapnia; (2) VT = 13 mL/kg; (3) VT = 22 mL/kg. In the second and third groups, RR was adjusted to yield mechanical power comparable to that of the first group. The VT level of 6 mL/kg combined with RR to lead to normocapnia (Samary et al., 2015) corresponds to the lower limit of protective mechanical ventilation, and has been shown to lead to low mechanical power (Samary et al., 2016b) in rats. The VT sizes of 13 and 22 mL/kg were chosen because they have been shown to promote VILI when mechanical power is not kept within a safe limit (Samary et al., 2016b). Settings were maintained for 2 h. Respiratory system mechanics were measured at 5 min (INITIAL), 60 min, and at the end of the experiment (FINAL). At FINAL, arterial blood gases were analyzed, heparin (1000 IU) was injected into the tail vein. The trachea was then clamped at the same airway pressure at end-expiration, at PEEP = 3 cmH2O, and all animals were killed by overdose of sodium thiopental (60 mg/kg i.v.) followed by transection of the abdominal aorta and vena cava. Lungs were then extracted for histology and molecular biology analysis.
Figure 1. Timeline representation of the experimental protocol. FDA, functional data acquisition; BGA, blood gas analysis; FiO2, fraction of inspired oxygen; LPS, Escherichia coli lipopolysaccharide; i.t., intratracheally; PEEP, positive end-expiratory pressure; PV curve, pressure-volume curve; RR, respiratory rate; VCV, volume-controlled ventilation; VT, tidal volume. Lung mechanics were assessed every 15 min.
Data Acquisition and Processing
A pneumotachograph (internal diameter = 1.5 mm, length = 4.2 cm, distance between side ports = 2.1 cm) was connected to the tracheal cannula for airflow (V′) measurements (Mortola and Noworaj, 1983). The pressure gradient across the pneumotachograph was determined using a SCIREQ differential pressure transducer (UT-PDP-02, SCIREQ, Montreal, QC, Canada). The flow resistance of the equipment, tracheal cannula included, was constant up to flow rates of 26 mL.s−1 and amounted to 0.12 cmH2O.mL−1.s −1. The equipment dead space was 0.3 mL. Airflow, airway pressure (Paw, measured at the tip of the tracheal cannula), and Pes were continuously recorded throughout the experiments with a computer running custom-made software written in LabVIEW (National Instruments, Austin, TX) (Samary et al., 2015; Spieth et al., 2015). Briefly, VT was calculated by digital integration of the airflow signal obtained from a custom-made pneumotachograph (Mortola and Noworaj, 1983) that was connected to the Y-piece of the ventilator tubing, while RR was calculated from the Pes swings as the frequency per minute of each type of breathing cycle.
Respiratory system and lung mechanics were calculated by occluding the airways at end-inspiration for 5 s until a respiratory system plateau pressure (Pplat,RS) and transpulmonary plateau pressure (Pplat,L), respectively, were reached (every 15 min) (Samary et al., 2015). Respiratory system driving pressure (ΔP,RS) was calculated as the difference between Pplat,RS (post-end-inspiratory pause) and PEEP. Transpulmonary pressure (P,L) was calculated as the difference between the Paw and Pes, whereas transpulmonary driving pressure (ΔP,L) was the difference between the transpulmonary pressure during end-inspiration (post-inspiratory pause) and end-expiration. The mechanical energy (Energy,L) was calculated based on the equation described by Guerin et al. (2016) and Marini and Jaber (2016) (simplified formula), as Energy, L = Δ/Est,L, where Est,L is the static lung elastance. Energy, L = Δ/(ΔP,L/VT) = ΔP,LxVT, which is the area of rectangle (see Supplementary Figure 1). Therefore, one must compute the area of the rectangle and divide the result by 2. Values were converted to mJ and multiplied by RR to obtain the mechanical power in mJ/min. This equation estimates the elastic work, the major component of the total work, without taking into account resistive properties and PEEP, unlike the equation proposed by Gattinoni et al. (2016). During the experiment, mechanical power was computed every 15 min to keep it at a low level. After finishing the experiments, all mechanical data were analyzed offline. At every 15 min of acquired data, a 3-to-5-min period of mechanical ventilation was assessed, and the values obtained during each of these periods were averaged for analysis.
All signals were amplified in a four-channel signal conditioner (SC-24, SCIREQ, Montreal, QC, Canada), and sampled at 200 Hz with a 12-bit analog-to-digital converter (National Instruments; Austin, Texas, USA). Mechanical data were computed offline by a routine written in MATLAB (Version R2007a; The Mathworks Inc., Natick, Massachusetts, USA).
The lungs and heart were removed en bloc. The left lung was frozen in liquid nitrogen and immersed in formaldehyde solution (4%), embedded in paraffin, cut longitudinally in the central zone by means of a microtome into three slices, each 4 μm thick, and stained with hematoxylin–eosin for histological analysis (Samary et al., 2015; Padilha et al., 2016). Photomicrographs at magnifications of ×100, ×200, and ×400 were obtained from eight nonoverlapping fields of view per section using a light microscope (Olympus BX51, Olympus Latin America Inc., Brazil). DAD was quantified using a weighted scoring system by two investigators (V.M. and V.L.C.) blinded to group assignment and independently, as described elsewhere (Kiss et al., 2016). Briefly, scores of 0–4 were used to represent edema, atelectasis, and overdistension, with 0 standing for no effect and 4 for maximum severity. Additionally, the extent of each scored characteristic per field of view was determined on a scale of 0–4, with 0 standing for no visible evidence and 4 for complete involvement. Scores were calculated as the product of severity and extent of each feature, on a range of 0–16. The cumulative DAD score was calculated as the sum of each score characteristic and ranged from 0 to 48.
Transmission Electron Microscopy
Three slices (2 × 2 × 2 mm) were cut from three different segments of the right lung and fixed in 2.5% glutaraldehyde in 0.1M Sorensen's Sodium-Potassium Phosphate Buffer. They were then stained en bloc in aqueous uranyl acetate for 24 h. After this procedure, samples were dehydrated in acetone and embedded in Araldite® resin. Ultrathin sections were obtained with a diamond knife on an LKB Ultratome microtome (Leica, Deerfield, IL, USA), placed on a 200 mesh copper grid (Lab Research Industries, Burlington, VT, USA), and counterstained with uranyl acetate and lead citrate. Micrographs were taken at 80 kV with a JEOL 100cx 100KW transmission electron microscope (Philips, Munich, Germany). For each lung electron microscopy image (20 fields per animal), type II epithelial cell damage, extracellular matrix (ECM) injury, and endothelial cell damage were graded on a five-point, semiquantitative, severity-based scoring system as follows: 0 = normal lung parenchyma, 1–4 = changes in 1–25%, 26–50%, 51–75%, and 76–100% of examined tissue, respectively (Samary et al., 2015). Electron microscopy analyses were performed by two investigators (V.M. and V.L.C.) blinded to group assignment.
Quantitative real-time reverse transcription polymerase chain reaction (PCR) was performed by an investigator (C.L.S.) blinded to group assignment in order to measure biomarkers associated with inflammation [interleukin (IL)-6], mechanical pulmonary stretch (amphiregulin), epithelial cell damage [club cell secretory protein 16 (CC16)], endothelial cell damage [intercellular adhesion molecular (ICAM)-1], extracellular matrix damage (syndecan-1 and decorin), and metalloproteinase (MMP)-9. The primers used are described in the Supplementary Table 1. Central slices of the right lung were cut, collected in cryotubes, flash-frozen by immersion in liquid nitrogen, and stored at −80 °C. Total RNA was extracted from frozen tissues using the RNeasy Plus Mini Kit (Qiagen, Hilden, Germany), following the manufacturer's recommendations. RNA concentrations were measured by spectrophotometry in a Nanodrop ND-1000 system (Thermo Scientific, Wilmington, DE, USA). First-strand cDNA was synthesized from total RNA using a Quantitec reverse transcription kit (Qiagen, Hilden, Germany). Relative mRNA levels were measured with a SYBR green detection system in an ABI 7500 real-time PCR system (Applied Biosystems, Foster City, California, USA). Samples were run in triplicate. For each sample, the expression of each gene was normalized to the acidic ribosomal phosphoprotein P0 (36B4) housekeeping gene (Akamine et al., 2007) and expressed as fold change relative to respective NV animals, using the 2−ΔΔCt method, where ΔCt = Ct (target gene) – Ct (reference gene) (Schmittgen and Livak, 2008).
The sample size was calculated to allow detection of the lowest, but still significant, changes in IL-6 between ventilatory strategies, which was expected between VT = 6 and 13 mL/kg and derived from a previous publication of our group (Samary et al., 2015). A sample size of 7 animals per group would provide the appropriate power (1 – β = 0.8) to identify significant differences in IL-6 (adjusted α = 0.025 for two comparisons), taking into account an effect size d = 2.0, a two-sided t-test, and a sample size ratio = 1 (G*Power 22.214.171.124, University of Düsseldorf, Düsseldorf, Germany).
Data were tested for normality using the Kolmogorov-Smirnov test with Lilliefors' correction, while the Levene median test was used to evaluate the homogeneity of variances. If both conditions were satisfied, time-dependent differences among groups were determined with general linear models followed by Bonferroni tests. Molecular biology, DAD, and electron microscopy variables were assessed with the Kruskal-Wallis test followed by Dunn's test. Multiple linear regression analyses between pHa, PaCO2, PaO2/FiO2, DAD score, IL-6, amphiregulin, CC16, decorin, syndecan-1, and MMP-9 (dependent variables, known as outcome variables) and VT, ΔP,L, and RR (independent variables, known as predictor variables), which represent the mathematical components of mechanical power, were performed. Parametric data were expressed as mean ± SD, while nonparametric data were expressed as median (interquartile range). General linear models and multiple linear regression analyses were carried out in SPSS 11.5.1 (SPSS Inc., Chicago, IL, USA). All the other tests were carried out in GraphPad Prism 6.00 (GraphPad Software, La Jolla, CA, USA). Significance was established at p < 0.05.
Respiratory parameters and blood gas exchange at BASELINE were comparable among groups (see Supplementary Table 2). Mean arterial pressure was >70 mmHg throughout the experiments in all groups. At FINAL, no significant differences among groups were observed in the volume of fluids required to keep MAP >70 mmHg (VT = 6 mL/kg: 12 ± 6 mL, VT = 13 mL/kg: 12 ± 1 mL, VT = 22 mL/kg: 14 ± 5 mL). Oxygenation was higher in VT = 13 mL/kg compared to VT = 6 mL/kg (p = 0.03) and VT = 22 mL/kg (p = 0.02). PaCO2 was higher whereas pHa was lower in VT = 22 mL/kg compared to VT = 13 mL/kg and VT = 6 mL/kg (p = 0.001 and p = 0.03, respectively) (Figure 2).
Figure 2. (A) PaO2/FiO2, (B) PaCO2, and (C) arterial pH at FINAL. The following groups were analyzed: (1) VT = 6 mL/kg and RR adjusted to normocapnia; (2) VT = 13 mL/kg; (3) VT = 22 mL/kg. In the second and third groups, RR was adjusted to yield mechanical power comparable to that in the first group. PaO2/FiO2 and PaCO2 values are medians and interquartile ranges of 7 rats in each group. pH values are mean (SD) of 7 animals/group. *vs. 6 mL/kg (p < 0.05); #vs. 13 mL/kg (p < 0.05).
With the progressive increase in VT, RR reduced significantly in order to keep comparable mechanical power among groups (Table 1). Airflow reduced significantly in VT = 22 mL/kg compared to VT = 6 and VT = 13 mL/kg. Pplat,RS, ΔP,RS, Pplat,L, ΔP,L and Energy,L were higher in VT = 22 mL/kg in comparison to VT = 6 and 13 mL/kg. Pplat,RS and ΔP,RS were higher in VT = 13 compared to 6 mL/kg, while no differences were observed in Pplat,L, ΔP,L, or Energy,L between them.
Figure 3 depicts light microscopy images of representative animals. VT = 22 mL/kg resulted in a higher DAD score compared to VT = 6 mL/kg, mainly due to overdistension. Electron microscopy images of one representative animal per group are shown in Figure 4. VT = 22 mL/kg increased damage to type II epithelial cells and endothelial cells compared to VT = 6 mL/kg. VT = 22 mL/kg also resulted in greater injury to extracellular matrix (ECM) than VT = 6 mL/kg and VT = 13 mL/kg. Other transmission electron microscopy images detailing damage to type II epithelial cells, ECM, and endothelial cells are depicted in Supplementary Figure 2.
Figure 3. Photomicrographs (light microscopy) of lung parenchyma stained with hematoxylin and eosin. The following groups were analyzed: (1) VT = 6 mL/kg and RR adjusted to normocapnia (A); (2) VT = 13 mL/kg (B); (3) VT = 22 mL/kg (C). In the second and third groups, RR was adjusted to yield mechanical power comparable to that in the first group. Asterisks show alveolar/interstitial edema; arrows indicate alveolar collapse, and double arrows, areas of alveolar overdistension. Photomicrographs are representative of data obtained from lung sections of seven animals (original magnification, 200×). Bars = 100 μm. Cumulative diffuse alveolar damage (DAD) score representing injury from alveolar/interstitial edema, atelectasis, and alveolar overdistension. Values are median (interquartile range) of 7 animals/group. *vs. VT = 6 mL/kg (p < 0.05).
Figure 4. Transmission electron microscopy of alveolar architecture. Photomicrographs are representative of data obtained from lung sections of 7 animals in each group. The following groups were analyzed: (1) VT = 6 mL/kg and RR adjusted to normocapnia (A); (2) VT = 13 mL/kg (B); (3) VT = 22 mL/kg (C). In the second and third groups, RR was adjusted to yield mechanical power comparable to that in the first group. Asterisk: extracellular matrix (ECM), PII: Type II epithelial cell, Cap: capillary, LB, lamellar body. PII damage is visible in all ARDS groups, regardless of tidal volume. LB decreased with increasing tidal volume. Moreover, the ECM was characterized by damage to collagen/elastic fibers (asterisk). Endothelial cells were less damaged at VT = 6 mL/kg and VT = 13 mL/kg. Semi-quantitative electron microscopy analysis of lung tissue. A five-point, semi-quantitative, severity-based scoring system was used: 0 = normal; or damage to 1 = 1–25%; 2 = 26–50%; 3 = 51–75%; and 4 = 76–100% of examined tissue. Values are median (interquartile range) of 7 animals/group. *vs. VT = 6 mL/kg (p < 0.05); #vs. VT = 13 mL/kg (p < 0.05).
IL-6, amphiregulin, and CC16 expressions were higher in VT = 22 mL/kg compared to VT = 6 mL/kg (Figure 5). Syndecan expression was higher in VT = 22 mL/kg, while decorin expression was lower, in comparison to VT = 6 mL/kg. MMP-9 was higher in VT = 22 mL/kg than VT = 6 mL/kg and 13 mL/kg (Figure 6).
Figure 5. Expression of biomarkers associated with inflammation, alveolar stretch, and epithelial and endothelial cell damage. Real-time polymerase chain reaction analysis of markers of inflammation [interleukin (IL)-6], alveolar stretch (amphiregulin), and epithelial [club cell protein (CC)16] and endothelial cell damage [intercellular adhesion molecule (ICAM)-1]. Relative gene expression was calculated as the ratio of average gene expression levels to reference gene (36B4) expression and presented as fold change relative to non-ventilated animals (NV). Values are medians and interquartile ranges of 7 rats in each group. Comparisons among all groups were done by the Kruskal-Wallis test followed by Dunn's test (p < 0.05). *vs. VT = 6 mL/kg (p < 0.05).
Figure 6. Expression of biomarkers associated with extracellular matrix damage and metalloproteinase. Real-time polymerase chain reaction analysis of markers of extracellular matrix damage (syndecan-1 and decorin) and of metalloproteinase (MMP)-9. Relative gene expression was calculated as the ratio of average gene expression levels to reference gene (36B4) expression and presented as fold change relative to non-ventilated animals (NV). Values are medians and interquartile ranges of 7 rats in each group. Comparisons among all groups were done by the Kruskal-Wallis test followed by Dunn's test (p < 0.05). *vs. VT = 6 mL/kg (p < 0.05). #vs. VT = 13 mL/kg (p < 0.05).
On multiple linear regression analyses, VT best predicted changes in IL-6 and CC16, while ΔP,L best predicted variance in pHa, PaO2/FiO2, and amphiregulin (Table 2). Syndecan was predicted by VT and ΔP,L. RR only predicted changes in decorin. The r2-values from multiple linear regression, when components of Energy,L (VT, ΔP,L, and RR) were used, were compared to r2-values from the linear regression performed using Energy,L which encompasses all of these components. Interestingly, the isolated components of Energy,L, mainly VT, predicted changes in IL-6 better than the Energy,L construct did (r2 = 0.71 vs. 0.19, respectively), while DAD was comparable (r2 = 0.46 vs. 0.47) (Table 2; see Supplementary Table 3).
In the rat model of mild-to-moderate ARDS used herein, we found that, at low mechanical power, higher VT resulted in increased PaCO2, DAD scores, and gene expression of mediators associated with inflammation (IL-6), alveolar stretch (amphiregulin), damage to epithelial cells (CC16), and extracellular matrix (MMP-9 and syndecan-1). Furthermore, only VT was able to predict changes in IL-6 and CC16, whereas ΔP,L predicted changes in pHa, oxygenation, amphiregulin, and syndecan-1 expression. Thus, even when mechanical power was low, high VT resulted in VILI. Changes in PaCO2 might be considered as a parameter to monitor for VILI development even at low mechanical power. Our findings suggest that VT is a better determinant of lung injury than mechanical power or RR.
To the best of our knowledge, this was the first study to investigate the impact of different VT levels, under the same low mechanical power, on lung function, DAD score, ultrastructural analysis of pulmonary parenchyma and gene expression of different biomarkers of VILI in experimental ARDS. We calculated the mechanical power by using a simplified equation to allow possible extrapolation of results to the clinical scenario (Marini and Jaber, 2016). PEEP was kept constant throughout the experiment; thus, different VT levels led only to dynamic energy changes. A particular strength of our study was the tight control of low mechanical power and the prospective nature of the investigation. Low mechanical power was ensured by ventilating animals with low VT levels (protective mechanical ventilation) and setting the RR to keep gas exchange and pHa within physiologic limits. Thus, minute ventilation was kept through the modulation of VT and RR. Lung damage was evaluated by light and electron microscopy and partitioned into its components, such as epithelial, endothelial, and extracellular matrix injury. Additionally, the severity of lung injury (mild to moderate) induced in this model corresponds to that observed in approximately two-thirds of ARDS cases (Bellani et al., 2016).
Recently, it was suggested that the various potential causes of VILI might be unified in a single variable, such as mechanical power. This variable is measured by the pressure–volume loop and can be computed from its components, such as VT, ΔP,L, airflow, PEEP, and RR. Thus, the relative contributions of each of these component variables to mechanical power as a determinant of VILI can be evaluated. Previous theoretical (Gattinoni et al., 2016) and experimental (Cressoni et al., 2016) studies suggested that lung damage is the result of the energy and power delivered to the lungs. The quantification of mechanical power at bedside and its interpretation might be helpful to guide optimization of settings for “safe” mechanical ventilation. The main hypothesis of the present study was that VILI at low mechanical power would be independent of delivered VT. In fact, low mechanical power may be achieved by setting lower VT and higher RR or vice-versa.
A previous theoretical study (Gattinoni et al., 2016) hypothesized that the main determinants of energy and power would be VT, ΔP,L, and airflow, followed by RR. Our findings suggest that, at low mechanical power, higher VT plays a greater role than RR as a determinant VILI, yielding higher respiratory system and lung Pplat and ΔP. RR has been recognized as an early potential factor underlying lung injury (Hotchkiss et al., 2000), mainly when associated with increased power (Cressoni et al., 2016).
Increased DAD scores and damage to the epithelium, extracellular matrix, and endothelial cells were observed only at VT = 22 mL/kg. Therefore, at low mechanical power, there was probably a threshold of inspiratory stress and dynamic strain after which the lung epithelium, endothelium, and ECM started to sustain further injury. Lung morphology changes were followed by increased expression of IL-6, a pro-inflammatory mediator, and biomarkers of epithelial and extracellular matrix damage. This suggests that the proportional role of VT in promoting lung damage is greater than predicted by the mechanical power model, which, in turn, is likely explained by the fact that high VT and the transpulmonary driving pressure used herein may have exceeded the safe limits of lung elasticity (Protti et al., 2011; Cressoni et al., 2015).
A previous study investigated the role of mechanical power in VILI in healthy piglets (Cressoni et al., 2016). Two experiments were attempted: VT was kept constant and RR increased progressively; and VT was increased while keeping the RR constant. Lung damage was evaluated through the increase in lung weight, assessed by CT scanning. They showed that lung edema increased while oxygenation decreased proportional to RR. We are unable to compare our results with those reported (Cressoni et al., 2016), due to major methodological differences: in the present study, lung damage was induced by endotoxin, mechanical power was kept low and constant while different combinations of VT and RR were trialed, and VILI was evaluated not by lung weight but by histological analysis under light and electron microscopy, as well as by gene expression of biomarkers associated with inflammation, alveolar stretch, epithelial and ECM damage. To date, no other experimental studies have assessed the impact of power on VILI in ARDS.
The damage caused by higher VT is not limited to epithelial and endothelial cells; it also affects ECM components. In this line, higher stress and strain may cause displacement or rupture of the ECM (Moriondo et al., 2007; Pelosi and Negrini, 2008; Protti et al., 2015). In the presence of high VT (high dynamic strain), MMP-9 activity has been shown to increase significantly (Gonzalez-Lopez et al., 2011; Jain et al., 2017). Accordingly, in our study, MMP-9 expression was higher in VT = 22 mL/kg compared to VT = 13 mL/kg and VT = 6 mL/kg. MMP-9 cleaves collagen, but also increases neutrophil influx, contributing to the ongoing inflammatory process (Gaggar and Weathington, 2016). Gene expression of syndecan-1, a cell surface heparan sulfate proteoglycan primarily expressed in the lung epithelium and involved in regulation of lung repair, remodeling (Brauer et al., 2016), and inflammation (Hayashida et al., 2009), was higher in VT = 22 mL/kg compared to VT = 6 mL/kg. Since the lung epithelium is damaged by high VT and syndecan-1 is expressed to regulate lung inflammation, we may postulate that syndecan-1 may also contribute to lung inflammation in VILI. Both VT and ΔP,L predicted changes in syndecan-1 expression, highlighting the sensitivity of this biomarker to detect changes in ECM homeostasis. In line with this result, decorin, which has been associated with inhibition of cell proliferation (Iacob et al., 2008) and collagen synthesis (Jahanyar et al., 2007), was lower in VT = 22 mL/kg compared to VT = 6 mL/kg. These findings suggest that, at low mechanical power, high VT increased ECM fragmentation (Moriondo et al., 2007; Pelosi and Negrini, 2008), which may promote lung inflammation. In this context, gene expression of IL-6, and early surrogate of inflammation and VILI, was higher in the presence of higher VT. The higher gene expression of amphiregulin at VT = 22 mL/kg might be attributed to overdistension (Dolinay et al., 2006). CC16, a protein mainly produced and secreted by club cells in the distal respiratory or terminal bronchioles, has been proposed as a biomarker of lung epithelial injury (Broeckaert et al., 2000) and its gene expression was also increased in VT = 22 mL/kg. According to multiple linear regression analyses, changes in VT better predicted changes in IL-6 and CC16 mRNA expressions. This can be explained by the importance of VT in promoting lung injury under the same low mechanical power. Changes in ΔP,L better reflected changes in amphiregulin mRNA expression, which has been observed in cell culture under cyclic deformations (Tschumperlin et al., 2000).
Even at low mechanical power, PaCO2 and end-expiratory alveolar overdistension also increased at higher VT. Since overall minute ventilation was comparable among groups, we hypothesize that alveolar dead space probably increased with lung injury. Previous studies showed that increased PaCO2 and dead space correlated with mortality (Nuckton et al., 2002), remodeling of lung parenchyma in ARDS (Gattinoni et al., 1993), and alveolar septal retraction due to high surface tension (Wilson, 1982). This suggests that monitoring changes in PaCO2 and dead space might be helpful to optimize protective mechanical ventilation settings, even at low mechanical power.
Possible Clinical Implications
Our data reinforce the concept that mechanical power per se might not be a useful parameter to guide optimization protective mechanical ventilation. To minimize VILI, VT seems to be the main ventilator parameter which should be controlled, irrespective of RR. In fact, the hypothesis that low mechanical power might be a major target to minimize VILI could erroneously result in acceptance of higher VT. Thus, even during mechanical ventilation at low mechanical power, inappropriately high VT can still promote lung injury. We thus insist that, even when using low RR and ΔP,L, VT should always be kept in low protective range. Certainly, this claim remains to be confirmed in clinical observational studies and randomized controlled trials.
This study has some limitations that should be considered. First, experimental ARDS was induced by endotoxin, and our results cannot be extrapolated to other models or to human ARDS. Second, as the observation time was only 2 h, we did not assess protein levels of biological markers. However, this duration of mechanical ventilation was sufficient to identify gene expression of markers of interest (Samary et al., 2015). Third, a simplified equation was chosen in order to facilitate its routine use in the clinical setting. This equation estimates the major component of the total work that is elastic work, without taking into account resistive properties and PEEP.
In the experimental model of mild-to-moderate ARDS used herein, at low mechanical power, high VT was associated with VILI. To mitigate VILI, VT control seems more important than RR control.
LM and PS participated in the design of the study, carried out the experiments, performed data analyses, and drafted the manuscript; AT, CS RS contributed to the study design and carried out the experiments; MF, VM, and VC carried out the experiments and lung morphology; LM performed analyses of lung mechanics; CS and MM carried out the molecular biology analyses and contributed to the manuscript; PP and MdA contributed to the study design, supervised the entire project, and helped write the manuscript; PR and PS contributed to the study design, supervised the experimental work and statistical analysis, wrote the manuscript, and supervised the entire project. All authors read and approved the final manuscript.
This work was supported by grants from the Carlos Chagas Filho Rio de Janeiro State Research Foundation (FAPERJ) (grant number E-26/103.118/2), Rio de Janeiro, Brazil; and the Brazilian Council for Scientific and Technological Development (CNPq) (grant number 469716/2014-2), Brasilia, Brazil.
Conflict of Interest Statement
The 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.
Assistance with the study: We thank Mr. Andre Benedito da Silva for animal care, Mrs. Arlete Fernandes for her help with microscopy, and Mrs. Moira Elizabeth Schottler and Mr. Filippe Vasconcellos for their assistance in editing the manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00318/full#supplementary-material
Akamine, R., Yamamoto, T., Watanabe, M., Yamazaki, N., Kataoka, M., Ishikawa, M., et al. (2007). Usefulness of the 5′ region of the cDNA encoding acidic ribosomal phosphoprotein P0 conserved among rats, mice, and humans as a standard probe for gene expression analysis in different tissues and animal species. J. Biochem. Biophys. Methods 70, 481–486. doi: 10.1016/j.jbbm.2006.11.008
Amato, M. B., Meade, M. O., Slutsky, A. S., Brochard, L., Costa, E. L., Schoenfeld, D. A., et al. (2015). Driving pressure and survival in the acute respiratory distress syndrome. N. Engl. J. Med. 372, 747–755. doi: 10.1056/NEJMsa1410639
Bellani, G., Guerra, L., Musch, G., Zanella, A., Patroniti, N., Mauri, T., et al. (2011). Lung regional metabolic activity and gas volume changes induced by tidal ventilation in patients with acute lung injury. Am. J. Respir. Crit. Care Med. 183, 1193–1199. doi: 10.1164/rccm.201008-1318OC
Bellani, G., Laffey, J. G., Pham, T., Fan, E., Brochard, L., Esteban, A., et al. (2016). Epidemiology, patterns of care, and mortality for patients with acute respiratory distress syndrome in intensive care units in 50 countries. JAMA 315, 788–800. doi: 10.1001/jama.2016.0291
Brauer, R., Ge, L., Schlesinger, S. Y., Birkland, T. P., Huang, Y., Parimon, T., et al. (2016). Syndecan-1 attenuates lung injury during influenza infection by potentiating c-met signaling to suppress epithelial apoptosis. Am. J. Respir. Crit. Care Med. 194, 333–344. doi: 10.1164/rccm.201509-1878OC
Broeckaert, F., Clippe, A., Knoops, B., Hermans, C., and Bernard, A. (2000). Clara cell secretory protein (CC16): features as a peripheral lung biomarker. Ann. N.Y. Acad. Sci. 923, 68–77. doi: 10.1111/j.1749-6632.2000.tb05520.x
Brower, R. G., Matthay, M. A., Morris, A., Schoenfeld, D., Thompson, B. T., and Wheeler, A. (2000). 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. 342, 1301–1308. doi: 10.1056/NEJM200005043421801
Cressoni, M., Chiurazzi, C., Gotti, M., Amini, M., Brioni, M., Algieri, I., et al. (2015). Lung inhomogeneities and time course of ventilator-induced mechanical injuries. Anesthesiology 123, 618–627. doi: 10.1097/ALN.0000000000000727
Cressoni, M., Gotti, M., Chiurazzi, C., Massari, D., Algieri, I., Amini, M., et al. (2016). Mechanical power and development of ventilator-induced lung injury. Anesthesiology 124, 1100–1108. doi: 10.1097/ALN.0000000000001056
Dolinay, T., Kaminski, N., Felgendreher, M., Kim, H. P., Reynolds, P., Watkins, S. C., et al. (2006). Gene expression profiling of target genes in ventilator-induced lung injury. Physiol. Genomics 26, 68–75. doi: 10.1152/physiolgenomics.00110.2005
Fan, E., Del Sorbo, L., Goligher, E. C., Hodgson, C. L., Munshi, L., Walkey, A. J., et al. (2017). An Official American Thoracic Society/European Society of Intensive Care Medicine/Society of Critical Care Medicine Clinical Practice Guideline: mechanical ventilation in adult patients with Acute Respiratory Distress Syndrome. Am. J. Respir. Crit. Care Med. 195, 1253–1263. doi: 10.1164/rccm.201703-0548ST
Gattinoni, L., D'Andrea, L., Pelosi, P., Vitale, G., Pesenti, A., and Fumagalli, R. (1993). Regional effects and mechanism of positive end-expiratory pressure in early adult respiratory distress syndrome. JAMA 269, 2122–2127. doi: 10.1001/jama.1993.03500160092039
Gattinoni, L., Tonetti, T., Cressoni, M., Cadringher, P., Herrmann, P., Moerer, O., et al. (2016). Ventilator-related causes of lung injury: the mechanical power. Intensive Care Med. 42, 1567–1575. doi: 10.1007/s00134-016-4505-2
González-López, A., Astudillo, A., García-Prieto, E., Fernández-García, M. S., López-Vázquez, A., Batalla-Solís, E., et al. (2011). Inflammation and matrix remodeling during repair of ventilator-induced lung injury. Am. J. Physiol. Lung cell. Mol. Physiol. 301, L500–L509. doi: 10.1152/ajplung.00010.2011
Guerin, C., Papazian, L., Reignier, J., Ayzac, L., Loundou, A., and Forel, J. M. (2016). Effect of driving pressure on mortality in ARDS patients during lung protective mechanical ventilation in two randomized controlled trials. Crit. Care 20:384. doi: 10.1186/s13054-016-1556-2
Hayashida, A., Bartlett, A. H., Foster, T. J., and Park, P. W. (2009). Staphylococcus aureus beta-toxin induces lung injury through syndecan-1. Am. J. Pathol. 174, 509–518. doi: 10.2353/ajpath.2009.080394
Hotchkiss, J. R. Jr., Blanch, L., Murias, G., Adams, A. B., Olson, D. A., Wangensteen, O. D., et al. (2000). Effects of decreased respiratory frequency on ventilator-induced lung injury. Am. J. Respir. Crit. Care Med. 161, 463–468. doi: 10.1164/ajrccm.161.2.9811008
Iacob, D., Cai, J., Tsonis, M., Babwah, A., Chakraborty, C., Bhattacharjee, R. N., et al. (2008). Decorin-mediated inhibition of proliferation and migration of the human trophoblast via different tyrosine kinase receptors. Endocrinology 149, 6187–6197. doi: 10.1210/en.2008-0780
Jahanyar, J., Joyce, D. L., Southard, R. E., Loebe, M., Noon, G. P., Koerner, M. M., et al. (2007). Decorin-mediated transforming growth factor-beta inhibition ameliorates adverse cardiac remodeling. J. Heart Lung Transplant. 26, 34–40. doi: 10.1016/j.healun.2006.10.005
Jain, S. V., Kollisch-Singule, M., Satalin, J., Searles, Q., Dombert, L., Abdel-Razek, O., et al. (2017). The role of high airway pressure and dynamic strain on ventilator-induced lung injury in a heterogeneous acute lung injury model. Intensive Care Med. Exp. 5:25. doi: 10.1186/s40635-017-0138-1
Kilkenny, C., Browne, W. J., Cuthill, I. C., Emerson, M., and Altman, D. G. (2010). Improving bioscience research reporting: the ARRIVE guidelines for reporting animal research. PLoS Biol. 8:e1000412. doi: 10.1371/journal.pbio.1000412
Kiss, T., Silva, P. L., Huhle, R., Moraes, L., Santos, R. S., Felix, N. S., et al. (2016). Comparison of different degrees of variability in tidal volume to prevent deterioration of respiratory system elastance in experimental acute lung inflammation. Br. J. Anaesth. 116, 708–715. doi: 10.1093/bja/aew093
Moriondo, A., Pelosi, P., Passi, A., Viola, M., Marcozzi, C., Severgnini, P., et al. (2007). Proteoglycan fragmentation and respiratory mechanics in mechanically ventilated healthy rats. J. Appl. Physiol. 103, 747–756. doi: 10.1152/japplphysiol.00056.2007
Nuckton, T. J., Alonso, J. A., Kallet, R. H., Daniel, B. M., Pittet, J. F., Eisner, M. D., et al. (2002). Pulmonary dead-space fraction as a risk factor for death in the acute respiratory distress syndrome. N. Engl. J. Med. 346, 1281–1286. doi: 10.1056/NEJMoa012835
Padilha, G. A., Horta, L. F., Moraes, L., Braga, C. L., Oliveira, M. V., Santos, C. L., et al. (2016). Comparison between effects of pressure support and pressure-controlled ventilation on lung and diaphragmatic damage in experimental emphysema. Intensive Care Med. Exp. 4:35. doi: 10.1186/s40635-016-0107-0
Protti, A., Andreis, D. T., Milesi, M., Iapichino, G. E., Monti, M., Comini, B., et al. (2015). Lung anatomy, energy load, and ventilator-induced lung injury. Intensive Care Med. Exp. 3:34. doi: 10.1186/s40635-015-0070-1
Protti, A., Cressoni, M., Santini, A., Langer, T., Mietto, C., Febres, D., et al. (2011). Lung stress and strain during mechanical ventilation: any safe threshold? Am. J. Respir. Crit. Care Med. 183, 1354–1362. doi: 10.1164/rccm.201010-1757OC
Putensen, C., Theuerkauf, N., Zinserling, J., Wrigge, H., and Pelosi, P. (2009). Meta-analysis: ventilation strategies and outcomes of the acute respiratory distress syndrome and acute lung injury. Ann. Intern. Med. 151, 566–576. doi: 10.7326/0003-4819-151-8-200910200-00011
Samary, C. S., Moraes, L., Santos, C. L., Huhle, R., Santos, R. S., Ornellas, D. S., et al. (2016a). Lung functional and biologic responses to variable ventilation in experimental pulmonary and extrapulmonary acute respiratory distress syndrome. Crit. Care Med. 44, e553–e562. doi: 10.1097/CCM.0000000000001611
Samary, C. S., Santos, R. S., Santos, C. L., Felix, N. S., Bentes, M., Barboza, T., et al. (2015). Biological impact of transpulmonary driving pressure in experimental acute respiratory distress syndrome. Anesthesiology 123, 423–433. doi: 10.1097/ALN.0000000000000716
Samary, C. S., Silva, P. L., Gama de Abreu, M., Pelosi, P., and Rocco, P. R. (2016b). Ventilator-induced lung injury: power to the mechanical power. Anesthesiology 125, 1070–1071. doi: 10.1097/ALN.0000000000001297
Schumann, S., Goebel, U., Haberstroh, J., Vimlati, L., Schneider, M., Lichtwarck-Aschoff, M., et al. (2014). Determination of respiratory system mechanics during inspiration and expiration by FLow-controlled EXpiration (FLEX): a pilot study in anesthetized pigs. Minerva Anestesiol. 80, 19–28.
Spieth, P. M., Silva, P. L., Garcia, C. S., Ornellas, D. S., Samary, C. S., Moraes, L., et al. (2015). Modulation of stress versus time product during mechanical ventilation influences inflammation as well as alveolar epithelial and endothelial response in rats. Anesthesiology 122, 106–116. doi: 10.1097/ALN.0000000000000415
Tschumperlin, D. J., Oswari, J., and Margulies, A. S. (2000). Deformation-induced injury of alveolar epithelial cells. Effect of frequency, duration, and amplitude. Am. J. Respir. Crit. Care Med. 162, 357–362. doi: 10.1164/ajrccm.162.2.9807003
Keywords: mechanical power, energy, extracellular matrix, driving pressure, respiratory rate, tidal volume
Citation: Moraes L, Silva PL, Thompson A, Santos CL, Santos RS, Fernandes MVS, Morales MM, Martins V, Capelozzi VL, de Abreu MG, Pelosi P and Rocco PRM (2018) Impact of Different Tidal Volume Levels at Low Mechanical Power on Ventilator-Induced Lung Injury in Rats. Front. Physiol. 9:318. doi: 10.3389/fphys.2018.00318
Received: 26 January 2018; Accepted: 15 March 2018;
Published: 04 April 2018.
Edited by:Yu Ru Kou, National Yang-Ming University, Taiwan
Reviewed by:Bradford Julian Smith, University of Colorado Denver, United States
Matteo Maria Pecchiari, Università degli Studi di Milano, Italy
Copyright © 2018 Moraes, Silva, Thompson, Santos, Santos, Fernandes, Morales, Martins, Capelozzi, de Abreu, Pelosi and Rocco. 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 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: Patricia R. M. Rocco, firstname.lastname@example.org
†These authors have contributed equally to this work.