key: cord-349980-x1h5dhn9
authors: Ge, Huiqing; Pan, Qing; Zhou, Yong; Xu, Peifeng; Zhang, Lingwei; Zhang, Junli; Yi, Jun; Yang, Changming; Zhou, Yuhan; Liu, Limin; Zhang, Zhongheng
title: Lung Mechanics of Mechanically Ventilated Patients With COVID-19: Analytics With High-Granularity Ventilator Waveform Data
date: 2020-08-21
journal: Front Med (Lausanne)
DOI: 10.3389/fmed.2020.00541
sha: 
doc_id: 349980
cord_uid: x1h5dhn9

Background: Lung mechanics during invasive mechanical ventilation (IMV) for both prognostic and therapeutic implications; however, the full trajectory lung mechanics has never been described for novel coronavirus disease 2019 (COVID-19) patients requiring IMV. The study aimed to describe the full trajectory of lung mechanics of mechanically ventilated COVID-19 patients. The clinical and ventilator setting that can influence patient-ventilator asynchrony (PVA) and compliance were explored. Post-extubation spirometry test was performed to assess the pulmonary function after COVID-19 induced ARDS. Methods: This was a retrospective study conducted in a tertiary care hospital. All patients with IMV due to COVID-19 induced ARDS were included. High-granularity ventilator waveforms were analyzed with deep learning algorithm to obtain PVAs. Asynchrony index (AI) was calculated as the number of asynchronous events divided by the number of ventilator cycles and wasted efforts. Mortality was recorded as the vital status on hospital discharge. Results: A total of 3,923,450 respiratory cycles in 2,778 h were analyzed (average: 24 cycles/min) for seven patients. Higher plateau pressure (Coefficient: −0.90; 95% CI: −1.02 to −0.78) and neuromuscular blockades (Coefficient: −6.54; 95% CI: −9.92 to −3.16) were associated with lower AI. Survivors showed increasing compliance over time, whereas non-survivors showed persistently low compliance. Recruitment maneuver was not able to improve lung compliance. Patients were on supine position in 1,422 h (51%), followed by prone positioning (499 h, 18%), left positioning (453 h, 16%), and right positioning (404 h, 15%). As compared with supine positioning, prone positioning was associated with 2.31 ml/cmH(2)O (95% CI: 1.75 to 2.86; p < 0.001) increase in lung compliance. Spirometry tests showed that pulmonary functions were reduced to one third of the predicted values after extubation. Conclusions: The study for the first time described full trajectory of lung mechanics of patients with COVID-19. The result showed that prone positioning was associated with improved compliance; higher plateau pressure and use of neuromuscular blockades were associated with lower risk of AI.

The novel coronavirus disease 2019 (COVID- 19) imposes an important and urgent threat to global health (1, 2) . A substantial proportion of COVID-19 cases will develop severe acute respiratory distress syndrome (ARDS) that requires invasive mechanical ventilation (IMV). The mortality rate of such patients can be as high as 40% (3), depending on comorbidities and the available medical resources. Mechanical ventilation is an important strategy to treat such patients; and lung mechanics can have both prognostic and therapeutic implications. Lung compliance is an important mechanical parameter that should be monitored during IMV. For example, lung recruitment maneuver (RM) has been used to improve lung compliance in severe ARDS (4) . There is also evidence in general ARDS population that poor lung compliance without improvement during IMV is associated with poor clinical outcome (5) . Patient ventilator asynchrony (PVA) is another important parameter that should be stressed during IMV. Risk factors of PVA has been widely investigated, including hours of the day, use of sedatives, ventilation mode and tidal volume (6, 7) . While several studies showed that PVA was associated with clinical outcome, others did not (8, 9) . There is preliminary opinion suggesting that lung mechanics of COVID-19 induced ARDS can be quite different from general ARDS (10) . However, there is no empirical data on the lung mechanics in COVID-19 patients on IMV. Furthermore, previous studies are limited in several aspects. First, there is no continuous pulmonary mechanics evaluation, including the response of lung recruitment during IMV, all events during prone ventilation. Second, most techniques for the detection of PVA and other parameters requires physical presence of an expert physician at the bedside and is thus only feasible during short periods (11) (12) (13) . In addition, most studies explored risk factors for PVA in a fixed-time model (14) . In reality, both risk factors and PVA and compliance were time-varying (15) .

In order to make this gap end, the purpose of the study were 4-folds: (1) to describe the lung mechanics of COVID-19 patients by analyzing high-granularity ventilator waveform data; (2) to explore whether the lung compliance can be influenced by clinical factors, such as recruitment maneuver (RM) and body positioning; (3) to identify risk factors for PVA during IMV in COVID-19 patients; and (4) To describe post-extubation lung functions for survivors with spirometry test.

Abbreviations: AI, asynchrony index; WOB, work of breathing; PEEP, positive end expiratory pressure; DT, delayed triggering; IEE, ineffective effort during expiration; IQR, interquartile range; COVID-19, coronavirus disease 2019; PVA, patient-ventilator asynchrony; ARDS, acute respiratory distress syndrome; IMV, invasive mechanical ventilation.

The study was conducted in the First People's hospital of Jingmen. Clinical data and ventilator wave data were retrospectively collected. All ventilator parameters were collected as longitudinally in hourly basis using a ventilator information system (RespCare TM , ZhiRuiSi Tech. Co., Ltd., Hangzhou, China). The impact of RM and positioning on lung compliance was explored in mixed linear model. The study was approved by the ethics committee of the First People's hospital of Jingmen (Approval number: 202002007) and the ethics committee of Sir Run Run Shaw hospital (20200407-32). Individual patient data were de-identified before analysis. Informed consent was waived as determined by the IRB due to retrospective nature of the study design.

All COVID-19 patients treated with IMV were included for analysis. COVID-19 was confirmed by one of the following criteria: (1) novel coronavirus nucleic acid was positive as confirmed by real time (RT)-PCT in respiratory or blood specimen; and (2) genetic sequencing showed highly homogenous sequence with the known novel coronavirus (16) . For adults with COVID-19 and acute hypoxemic respiratory failure despite conventional oxygen therapy (<92%), we would start using high-flow nasal cannula (HFNC) or non-invasive ventilation (NIV). If the condition further deteriorated and the oxygenation saturation could not be maintained above 92% with HFNC or NIV, IMV would be started (17) . Patients were excluded if (1) they were younger than 18 years old; (2) patients with do-not-resuscitate order and (3) with terminally ill disease; (4) patients with incomplete record of waveform data.

Demographic data including age and sex were collected as time-fixed data. Hospital mortality was obtained on discharge. Pulmonary functions including forced vital capacity (FVC), forced expiratory volume (FEV1), FEV1/FVC ratio, peak expiratory flow (PEF), Peak inspiratory flow (PIF), maximal inspiratory pressure (MIP), maximal expiratory pressure (MEP) were measured for hospital survivors.

Ventilator parameters including lung compliance, measured PEEP, plateau pressure, tidal volume, work of breathing (WOB), and peak flow rate were measured based on pressure and flow waveforms. Details of the measurement approaches are described in the ESM.

Interventions including RM, positioning, sedatives and neuromuscular blockades were recorded in our analysis. Date and time of these interventions used to match to a period when ventilator parameters and lung mechanics were recorded. The body position was recorded as one of supine, right, left and prone positions at a specific time. Non-supine position was applied during daytime, and the specific positioning (prone, right or left) was determine at the discretion of the attending physician and respiratory therapist depending on the improvement in oxygenation. Prone positioning was applied for at least 10 h one day. RM could be accurately identified from ventilator waves as those with more than 30 cmH 2 O sustained inflation maintained for at least 30 s, the upper limit pressure was 45 cmH 2 O (18).

We developed an interpretable deep learning approach to detect double triggering (DT) and ineffective inspiratory effort during expiration (IEE). Individual deep learning models were developed under all ventilation modes. Under each ventilation mode, two models were established for detecting DT and IEE. Each model uses the raw ventilator waveforms (airway pressure and flow) as input for a binary classification (PVA or non-PVA). It is also capable of explaining the classification by highlighting the segments that contributes mostly to the results. Datasets were annotated by a group of clinical professionals for training and validating the models based on our previously proposed approach (19) . The accuracy reached above 95% for both types of PVA in all the ventilation modes. Asynchrony index (AI) was calculated as the number of asynchronous events divided by the number of ventilator cycles and wasted efforts (14) . Details of the algorithm development is described in the Electronic Supplemental Material.

Ventilator parameters were described for each individual patient by median and interquartile range (IQR) (20) . Temporal trends of ventilator parameters were visualized with scatter plots and described with Locally Weighted Scatterplot Smoothing (LOWESS) curves (21) . These curves were drawn for each individual patient and survivors and non-survivors were denoted with different colors. Risk factors for IEE and DT were explored with mixed negative binomial regression models, which was a generalization of the Poisson regression allowing for the conditional variance exceeds the conditional mean (22) . Random-effects was allowed for intercepts to account for between-subject variance. Predictors of IEE and DT included compliance, plateau pressure, PEEP, TV, respiratory rate, peak flow rate, WOB, sedatives, and neuromuscular blockades. We reported relative risk (RR) for the risk estimate associated with a unit change of these predictors. Risk factors for AI was explored with mixed linear effect model because the response variable AI was in linear scale. We reported coefficient and 95% confidence interval (CI) to represent how AI increased with a unit change in predictors. Factors that can influence lung compliance was explored with a mixed-effects linear model. Factors including age, sex, RM, PEEP, AI, and body position were included in the model. All statistical analyses were performed with RStudio (Version 1.1.463) . A two-tailed p < 0.05 were considered as statistical significance.

A total of 7 patients with full record of ventilator waveforms were included for analysis. There was no excluded patient due to predefined exclusion criteria. Four patients died and three survived to hospital discharge ( ; p < 0.001] than non-survivors. Survivors were more likely to adopt prone position than non-survivors (21 vs. 16%; p < 0.001). All RM was performed in non-survivors. More neuromuscular blockades were used in non-survivors ( Table 2) . Patients were on supine position in 1,422 h (51%), followed by prone positioning (499 h, 18%), left positioning (453 h, 16%), and right positioning (404 h, 15%). Survivors showed increasing compliance over time, whereas non-survivors showed persistently low compliance ( Figure 1A) . Plateau pressure, PEEP and tidal volume are shown in Figures 1B-D. WOB and respiratory rate did not show difference between survivors and non-survivors in temporal pattern (Figures 1E,F) . Temporal trends of PVA were not different between survivors and nonsurvivors (Figure 2 ).

Risk factors for PVA (IEE and DT) were investigated in the mixed negative binomial regression models. Higher plateau pressure (RR: 0.945; 95% CI: 0.934-0.956; p < 0.001) and respiratory rate (RR: 0.963; 95% CI: 0.951-0.976; p < 0.001) was associated with less IEE. However, greater tidal volume and WOB were associated with more IEE. In contrast to IEE, higher respiratory rate was associated with increased risk of DT (RR: 1.066; 95% CI: 1.054-1.078; p < 0.001). Higher plateau pressure (Coefficient: −0.90; 95% CI: −1.02 to −0.78) and neuromuscular blockades (Coefficient: −6.54; 95% CI: −9.92 to −3.16) were associated with lower AI. Sedatives had no significant impact on PVAs ( Table 3) . 

In multivariable mixed-effects linear model, we found two variables were significantly associated with lung compliance. Each 1 cmH 2 O increase in PEEP was associated 0.27 ml/cmH 2 O decrease in lung compliance (95% CI: −0.36 to −0.18; p < 0.001). As compared with supine positioning, prone positioning was associated with 2. 31 

Spirometry tests were performed in survivors at day 8, 11, and 13 after extubation. It showed that FVC was consistently decreased for the three measurements. FEV1/FVC was decreased in patient 3 (0.73 at day 8 and 0.707 at day 11); but was preserved in 

The study integrated high-granularity ventilator waveform data with clinical variables to describe the temporal change of lung mechanics of critically ill patients with COVID-19. At the time of intubation, the lung compliance was similar in survivors and non-survivors; but the survivors showed gradually improved compliance. Prone positioning is effective in improve lung compliance. Two types of PVA, IEE, and DT, were identified with deep learning algorithm. Higher plateau pressure and use of muscular relaxant were associated with lower risk of PVAs. Spirometry tests showed that pulmonary functions were significantly compromised after recovery from COVID-19 induced ARDS. Long-term follow up for the change of pulmonary functions would be relevant. Although the lung compliance was similar at the time of intubation, survivors showed gradual improvement in lung compliance, while non-survivors showed persistently low lung compliance. This is consistent with other studies that lung compliance was an independent predictor of mortality (5, 23, 24) . An important finding in our study was that RM was not effective in improving lung compliance, which is in contrast to findings from general ARDS patients. Although the effect of RM on mortality was conflicting in general ARDS, it has been consistently reported to be able to improve lung compliance (25) (26) (27) . For example, Kung and colleagues observed that the respiratory system compliance was significantly higher in the RM group from day 1 to 7 (25) . There is evidence that direct/pulmonary ARDS is more responsive to RM than indirect/extrapulmonary ARDS. While only 21% patients with lower percentage of recruitable lung were caused by pulmonary ARDS, 51% patients with higher percentage of recruitable lung caused by pulmonary ARDS (p = 0.01) (28) . Thus, COVID-19 induced ARDS is pulmonary ARDS but is less responsive to RM as shown in our study. The second reasons may be due to the fact that we only employed sustained inflation RM. Since there are many types of RM, it is largely unknown whether other types of RM can be effective in improve lung compliance in COVID-19 patients. Finally, the ARDS in COVID-19 may be due to viral, bacterial, or any kind of lung insults. Thus, the RM should not be able to demonstrate the benefits in this group of patients.

PVAs are commonly observed in patients with IMV, especially those with protective ventilation strategy. Our study observed that AI was 4.95% (IQR: 1.69-18.93) in overall observed hours. Non-survivors had more AI than survivors, indicating AI is a risk factor for mortality, which was consistent with other studies (29) .

Ventilator parameters can have differing effects on different PVA types. For example, while higher respiratory rate was associated with lower risk of IEE, it was associated with higher risk of DT. Use of neuromuscular blockades was associated with lower risk of both IEE and DT. However, we did not observe significant effect of sedatives on AI. Other studies have shown that Propofol or other sedatives can reduce AI (13, 30) . It is not surprising to observe that neuromuscular blockades are associated with significantly reduced risk of PVAs.

Post-extubation pulmonary function has never been reported for COVID-19 patients. Our results indicated that pulmonary functions can be significantly compromised in a short period. The FVC is reduced to one third of the predicted value. Other pulmonary function parameters, such as PEF and MIP were also reduced by one third of the predicted value. Boucher and coworkers observed that the pulmonary function can be significantly compromised in pediatric ARDS in short follow-up period (31) . In adult patients with general ARDS, the FVC can recover to 3.34 ± 0.77 and 3.78 ± 1.11 L at 1 and 6 months follow up (32), which is significantly higher than that in our study. However, since we did not obtain the long term follow up data, it is largely unknown whether COVID-19 can have long-term effect on pulmonary functions.

Several limitations should be acknowledged in the study. First, the sample size was limited, which prohibited patient-level analysis. The effect of prone-positioning on mortality outcome could not be analyzed with sufficient statistical power. Thus, further large-scale studies are needed to validate our findings. However, our data is rich with high-granularity waveform data, which allows for patient-hour analysis for epidemiological analysis. Second, we only developed deep learning algorithms for identifying two types of PVA. There are other types of PVAs, such as reverse triggering and short/long cycling. However, these analyses are not applicable in pressure-controlled ventilation and pleural pressure is required for reverse triggering (33) . Third, the impact of sedative on PVAs were estimated without the dosing of sedatives. We only recorded the use of sedatives as a binary variable. Such treatment would lose some information but is easy to interpret because different sedatives imposes challenge to standardize the dose. Finally, the pulmonary function was measured in a short period of time; long-term follow up data may provide important information for critically ill COVID-19 patients.

In conclusion, the study for the first time described full trajectory of lung mechanics of patients with COVID-19. The result showed that prone positioning was associated with improved compliance; higher plateau pressure and use of neuromuscular blockades were associated with lower risk of AI. RM was not associated with improvement on compliance.

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

The study was approved by the ethics committee of the First People's Hospital of Jingmen (Approval number: 202002007) and the ethics committee of Sir Run Run Shaw hospital (20200407-32). Informed consent was waived as determined by the IRB due to retrospective nature of the study design.

Clinical characteristics of coronavirus disease 2019 in China

Epidemiologic characteristics of early cases with 2019 novel coronavirus (2019-nCoV) disease in Republic of Korea

Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study

Effect of lung recruitment maneuver on oxygenation, physiological parameters and mortality in acute respiratory distress syndrome patients: a systematic review and meta-analysis

Estimated dead space fraction and the ventilatory ratio are associated with mortality in early ARDS

Identifying and managing patient-ventilator asynchrony: an international survey

Patient-ventilator asynchrony during daily interruption of sedation versus no sedation protocol

Asynchronies during mechanical ventilation are associated with mortality

Clusters of ineffective efforts during mechanical ventilation: impact on outcome

Covid-19 does not lead to a "typical" acute respiratory distress syndrome

Defining patient-ventilator asynchrony severity according to recurrence

Patient-ventilator asynchronies during mechanical ventilation: current knowledge and research priorities

Effects of dexmedetomidine and propofol on patient-ventilator interaction in difficult-to-wean, mechanically ventilated patients: a prospective, open-label, randomised, multicentre study

Predictors of asynchronies during assisted ventilation and its impact on clinical outcomes: the EPISYNC cohort study

Time-varying covariates and coefficients in Cox regression models

A rapid advice guideline for the diagnosis and treatment of 2019 novel coronavirus (2019-nCoV) infected pneumonia (standard version)

Surviving sepsis campaign: guidelines on the management of critically ill adults with coronavirus disease 2019 (COVID-19)

Optimal duration of a sustained inflation recruitment maneuver in ARDS patients

Detection of patientventilator asynchrony from mechanical ventilation waveforms using a twolayer long short-term memory neural network

Univariate description and bivariate statistical inference: the first step delving into data

LOWESS: a program for smoothing scatterplots by robust locally weighted regression

A flexible mixed-effect negative binomial regression model for detecting unusual increases in MRI lesion counts in individual multiple sclerosis patients

Impact of lung compliance on neurological outcome in patients with acute respiratory distress syndrome following out-of-hospital cardiac arrest

Driving pressure is associated with outcome during assisted ventilation in acute respiratory distress syndrome

Effects of stepwise lung recruitment maneuvers in patients with early acute respiratory distress syndrome: a prospective, randomized, controlled trial

Lung recruitment maneuvers for adult patients with acute respiratory distress syndrome. A systematic review and meta-analysis

Variation of poorly ventilated lung units (silent spaces) measured by electrical impedance tomography to dynamically assess recruitment

Lung recruitment in patients with the acute respiratory distress syndrome

Patient-ventilator asynchrony in adult critically ill patients

Dexmedetomidine vs midazolam or propofol for sedation during prolonged mechanical ventilation: two randomized controlled trials

Post-discharge respiratory outcomes of children with acute respiratory distress syndrome

Driving pressure and long-term outcomes in moderate/severe acute respiratory distress syndrome

Mechanical ventilation-induced reverse-triggered breaths: a frequently unrecognized form of neuromechanical coupling

HG, QP, and ZZ analyzed and interpreted the results and drafted the manuscript. JZ, CY, and LZ handled the FINNAKI data. All authors took part in designing the study, revised the manuscript critically for important intellectual content, read, and approved the final manuscript. 

The authors would like thank Biyun Lai and Jie Pan from College of Information Engineering, Zhejiang University of Technology for preparing the data, ZhiRuiSi Tech. Co., Ltd. for donating the RespCareTM system and Mr. Bing Li from ZhiRuiSi for deploying the system on the scene to help the current research work.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed. 2020.00541/full#supplementary-material