Abstract
BACKGROUND: The impact of mechanical ventilation parameters and management on outcomes of patients with coronavirus disease 2019 (COVID-19) ARDS is unclear.
METHODS: This multi-center observational study enrolled consecutive mechanically ventilated patients with COVID-19 ARDS admitted to one of 7 Korean ICUs between February 1, 2020–February 28, 2021. Patients who were age < 17 y or had missing ventilation parameters for the first 4 d of mechanical ventilation were excluded. Multivariate logistic regression was used to identify which strategies or ventilation parameters that were independently associated with ICU mortality.
RESULTS: Overall, 129 subjects (males, 60%) with a median (interquartile range) age of 69 (62–78) y were included. Neuromuscular blocker (NMB) use and prone positioning were applied to 76% and 16% of subjects, respectively. The ICU mortality rate was 37%. In the multivariate analysis, higher dynamic driving pressure (ΔP) values during the first 4 d of mechanical ventilation were associated with increased mortality (adjusted odds ratio 1.16 [95% CI 1.00–1.33], P = .046). NMB use was associated with decreased mortality (adjusted odds ratio 0.27 [95% CI 0.09–0.81], P = .02). The median tidal volume values during the first 4 d of mechanical ventilation and the ICU mortality rate were significantly lower in the NMB group than in the no NMB group. However, subjects who received NMB for ≥ 6 d (vs < 6 d) had higher ICU mortality rate.
CONCLUSIONS: In subjects with COVID-19 ARDS receiving mechanical ventilation, ΔP during the first 4 d of mechanical ventilation was independently associated with mortality. The short-term use of NMB facilitated lung-protective ventilation and was independently associated with decreased mortality.
- COVID-19
- driving pressure
- mechanical ventilation
- neuromuscular blockade
- respiratory distress syndrome
- SARS-CoV-2
Introduction
As of August 21, 2021, over 209 million confirmed cases of coronavirus disease 2019 (COVID-19) have been reported worldwide, and over 4 million related deaths have been confirmed (https://www.who.int/emergencies/diseases/novel-coronavirus-2019. Accessed August 21, 2021). Many of those who were hospitalized for COVID-19 have received mechanical ventilation due to severe pneumonia and ARDS.1,2 ARDS is characterized by the acute onset of noncardiogenic pulmonary edema and hypoxemia that could be aggravated by ventilator-induced lung injury.3 The mortality of patients with severe COVID-19 admitted to the ICU is high, ranging from 15–74%.4,5
Early reports demonstrated that COVID-19 ARDS may have a different physiology and can be categorized into 2 subphenotypes according to respiratory mechanics: “early” ARDS with high compliance and low recruitability and “late” ARDS with low compliance and high recruitability.6,7 Thus, the applicability of current evidence-based ARDS management, such as low tidal volume (VT), PEEP, and prone positioning, for patients with COVID-19 ARDS is controversial. Although numerous studies have described ventilation parameters and strategies in mechanically ventilated patients with COVID-19,8 their potential impact on patient outcomes is unclear. For instance, although neuromuscular blockers (NMBs) and prone positioning have been extensively studied in subjects with non-COVID-19 ARDS,9,10 studies on their impact on mortality of patients with COVID-19 are limited.
Therefore, this study aimed to investigate the mechanical ventilation settings and management strategies independently associated with ICU mortality in subjects with COVID-19 ARDS.
QUICK LOOK
Current Knowledge
Given the recent controversy concerning differences in physiology between coronavirus disease 2019 (COVID-19) and non-COVID-19 ARDS, numerous studies have described ventilation parameters and management in patients with COVID-19. However, there are limited data on the potential impact on mortality.
What This Paper Contributes to Our Knowledge
Higher dynamic driving pressure, peak inspiratory pressure, and mechanical power values during the first 4 d of mechanical ventilation were independently associated with ICU mortality in subjects with COVID-19 ARDS. PEEP, tidal volume, compliance, and ventilatory ratio were not. In subjects who received a neuromuscular blocker, tidal volume values during the first 4 d of mechanical ventilation and ICU mortality rate were significantly lower. Prone positioning was not significantly associated with decreased mortality.
Methods
Setting and Subjects
The Korean study group on the oxygenated or ventilated subjects with severe COVID-19 (KOVIDS) is a multi-center, observational cohort study that enrolled adult (≥ 17 y old) subjects with COVID-19 who were admitted to one of the ICUs at the 7 tertiary or referral hospitals in South Korea. The institutional review boards (IRBs) of the participating hospitals approved the study protocol, and the need for informed consent was waived due to the observational nature of the study (IRB of Chung-Ang University Hospital, No. 2103–009-19360). Subjects’ electronic medical records were reviewed by intensivists or research nurses trained in critical care, collecting de-identified subject data on demographics, comorbidities, dates of first symptom and COVID-19 diagnosis, oxygen support devices (nasal cannula, oxygen mask, high-flow nasal cannula, or mechanical ventilation), in-hospital medications and procedures, vital signs, laboratory and microbiological data, and outcomes according to a standardized case record form. Each investigator added the data collected from this form to an online database (Google Sheets, Google, Mountain View, California). Two independent investigators (WYK and MSB) screened the database for errors and verified data by visual inspection of distribution plots. Local investigators were contacted for any queries; then validated data were entered into the final database. The recommendations from the Surviving Sepsis Campaign guidelines for critically ill adults with COVID-19, including standard management of ARDS, were followed in all hospitals.11 PEEP/FIO2 titration was not standardized but was left at the discretion of the attending physician. This study complied with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement.
For the current study, all consecutive subjects with COVID-19 registered in the data set between February 1, 2020–February 28, 2021 in whom mechanical ventilation was applied for ARDS were included. Patients were excluded if they were age < 17 y, were not intubated, did not meet the diagnostic criteria for ARDS, had a do-not-intubate order, or had missing ventilation parameter values during the first 4 d of mechanical ventilation. The subjects were followed until hospital discharge, death, or June 17, 2021, which was the date of last chart review. Survival status in the 60 d following hospital admission was determined for all subjects.
Data Collection and Definitions
Baseline data collected at hospital admission included age, sex, body mass index, smoking status, comorbidities, date of first symptom, date of hospital admission, and laboratory findings. In-hospital admission data included the date of mechanical ventilation initiation, corticosteroid, and rescue therapies for ARDS (NMB, prone positioning, and extracorporeal membrane oxygenation [ECMO]). The severity of ARDS, Sequential Organ Failure Assessment (SOFA) score,12 and arterial blood gas data at the time of mechanical ventilation initiation were recorded. We also collected the daily mean values of the following ventilation parameters for the first 4 d of mechanical ventilation, with day 1 indicating the first day of intubation: breathing frequency, VT, PEEP, and peak inspiratory pressure (PIP). The parameters measured before mechanical ventilation optimization were not included.
Laboratory confirmation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was based on a positive result according to a real-time reverse transcription polymerase chain reaction assay from upper respiratory tract specimens.13 Immunosuppression was identified when the subject had malignancy, organ transplant, or immunosuppressive medications. ARDS was defined based on the Berlin criteria and classified into mild (PaO2/FIO2 200–300 mm Hg), moderate (PaO2/FIO2 100–200 mm Hg), and severe (PaO2/FIO2 ≤ 100 mm Hg).14 The use of NMB was defined as continuous administration within 24 h after intubation to facilitate mechanical ventilation or prone positioning. Prone positioning was defined as any episode of proning initiated within the first 2 d of mechanical ventilation. VT was expressed in mL/kg ideal body weight (IBW). Dynamic driving pressure (ΔP) was calculated as PIP, PEEP, and dynamic compliance as VT/ΔP.15 We used dynamic ΔP because plateau pressure was not measured for subjects who received pressure-control continuous mandatory ventilation. Mechanical power was defined as 0.098 × VT × breathing frequency × (PIP 1/2 × ΔP).16 Ventilatory ratio was calculated as (minute ventilation × PaCO2)/(IBW × 100 × 37.5).17
Outcome Measures
The primary outcome measure was ICU mortality. The secondary outcome measures were (1) ventilator weaning, defined as freedom from mechanical ventilation for ≥ 48 h; (2) ventilator-free days at day 28, defined as the number of d that the subject was alive and free from mechanical ventilation within 28 d after intubation; (3) hospital length of stay; (4) tracheostomy; (5) renal replacement therapy during ICU stay; (6) superinfection; and (7) hospital, 28-d, and 60-d mortality.
Statistical Analysis
Data were reported as the median (interquartile range [IQR]) or the mean ± SD for continuous variables and as the number (percentage) for categorical variables. Continuous variables were compared using Mann-Whitney U test, whereas categorical variables were compared using the chi-square test or Fisher exact test. Continuous variables were compared among more than 2 groups with the Kruskal-Wallis test. Ventilation parameters during the first 4 d of mechanical ventilation were described by line graphs. Kaplan-Meier survival curves until day 60 were generated and compared between groups using log-rank tests.
Multivariate logistic regression was performed to quantify the association of mechanical ventilation settings and ventilation management with ICU mortality, with adjustment for other covariates. Independent variables were selected among variables with P values of < .10 in the univariate analysis and potential confounders that were judged based on clinical expertise. The national health system was extremely pressured during the third wave of the epidemic in Korea (December 2020–January 2021). Therefore, the potential effect of changes in the treatment period was assessed by including hospital admission date (November 30, 2020, and earlier vs December 1, 2020, and later) in the model. For ventilation parameters, the mean values from the first 4 d of mechanical ventilation were considered in the model. As serum creatinine showed a high collinearity with the SOFA score at the start of mechanical ventilation, serum creatinine was excluded in the final model because it was already included in the SOFA renal subscore. Similarly, PEEP, PIP, ΔP, and mechanical power showed high collinearity (Table S1, see supplementary material related to this paper at http://rc.rcjournal.com), but only ΔP was included in the final model because previous studies demonstrated that ΔP is the optimal ventilation parameter for stratifying the risk of death among subjects with ARDS receiving mechanical ventilation.18,19 After undergoing final assessment, a multivariate model was constructed using stepwise backward selection. The model’s discrimination was assessed with C statistic (c = 0.82), and its calibration was assessed with Hosmer-Lemeshow test (chi-square = 4.93, P = .77). To assess the possible association of PEEP, PIP, and mechanical power with ICU mortality, we performed a secondary multivariate analysis replacing ΔP by those variables.
No imputation strategy was conducted because only few variables had missing data (13/5,516, 0.2%). There were no missing values for ventilation parameters. All statistical analyses were performed using IBM SPSS software (version 26.0; IBM, Armonk, New York). A 2-sided P value of < .05 was considered statistically significant.
Results
Subject Characteristics
Among the 1,575 patients evaluated for eligibility (Fig. 1), 1,446 patients were excluded due to no oxygen support (n = 1,077), conventional oxygen therapy only (n = 184), high-flow nasal cannula only (n = 138), and a do-not-intubate order (n = 47). Finally, 129 mechanically ventilated subjects were included in the study, and all met the criteria for ARDS. The median duration of follow-up was 27 d (IQR 20–44 d). The baseline subject characteristics are shown in Table 1. The median age, body mass index, and SOFA score at the start of mechanical ventilation were 69 y (IQR 62–78 y), 25.3 kg/m2 (IQR 22.5–27.8 kg/m2), and 9 (IQR 8–11), respectively. In total, 60% and 18% of the subjects were male and ever smokers, respectively. The most frequent comorbidities were hypertension (60%) and diabetes (41%). The median time from symptom onset to hospital admission was 5 d (IQR 2–9 d). Corticosteroid was used in 92%, NMB in 76%, prone positioning in 16%, and ECMO in 16%. ICU nonsurvivors were older, had higher SOFA score at the start of mechanical ventilation, had lower NMB use, and had lower pH and serum bicarbonate at the start of mechanical ventilation. There was no significant difference in the median time from symptom onset to hospital admission between survivors and nonsurvivors.
Overall, 19%, 43%, and 38% of subjects had mild, moderate, and severe ARDS on day 1 of mechanical ventilation, respectively (Table S2, see supplementary material related to this paper at http://rc.rcjournal.com). The median SOFA score at the start of mechanical ventilation was higher in subjects with severe ARDS. However, the degree of ARDS severity was not associated with significant differences in the use of corticosteroids, NMB, prone positioning, or ECMO.
Ventilation Parameters
The median time from hospital admission to mechanical ventilation initiation was 1 d (IQR 0–3 d), with no significant difference between survivors and nonsurvivors (Table 1). All subjects received pressure- or volume-control continuous mandatory ventilation. On day 1 of mechanical ventilation, the median VT, PEEP, PIP, ΔP, and mechanical power were 6.9 mL/kg IBW (IQR 6.2–8.1 mL/kg IBW), 10 cm H2O (IQR 9–12 cm H2O), 24 cm H2O (IQR 21–28 cm H2O), 14 cm H2O (IQR 11–16 cm H2O), and 26.0 J/min (IQR 21.3–32.5 J/min), respectively (Table 1). Dis-tributions of the median VT and PEEP were similar between survivors and nonsurvivors, but the median PIP, ΔP, and mechanical power were significantly higher in nonsurvivors than in survivors. In addition, the median compliance was lower and median ventilatory ratio was higher in nonsurvivors than in survivors. However, these findings did not reach statistical significance. The ventilation strategies (VT and PEEP) did not vary with the degree of ARDS severity (Table S2, see supplementary material related to this paper at http://rc.rcjournal.com). However, the median PIP, ΔP, mechanical power, and ventilatory ratio increased with the severity of ARDS.
During the first 4 d of mechanical ventilation, the median VT values of 6.8–7.3 mL/kg IBW were maintained in both survivors and nonsurvivors (Fig. 2). The differences in the median PEEP between survivors and nonsurvivors were not clinically relevant, although statistically significant differences were observed on day 3 and day 4 of mechanical ventilation. The median PIP and ΔP values were significantly higher in nonsurvivors than in survivors, despite both groups receiving median PIP values of < 25 cm H2O and median ΔP values of < 15 cm H2O. The median compliance values were significantly lower in nonsurvivors on day 3 and day 4. The median mechanical power values were significantly higher in nonsurvivors during the first 4 d of mechanical ventilation.
Clinical Outcomes
The rates of ICU mortality, ventilator weaning, and superinfection were 37%, 60%, and 53%, respectively. The median ventilator-free d at day 28 was 8.2 ± 9.7 d, and the median hospital length of stay was 27 d (IQR 20–44 d). Tracheostomy was performed in 33%, and 25% received renal replacement therapy during their ICU stay. The ICU mortality rates were 20% in mild, 40% in moderate, and 43% in severe ARDS (P = .13) (Table S2, see supplementary material related to this paper at http://rc.rcjournal.com). Other study outcomes were not generally affected by the degree of ARDS severity, except for the highest superinfection rate in subjects with severe ARDS.
Association of Ventilation Parameters and Management With ICU Mortality
After adjustment for confounding factors, the independent predictors of ICU mortality were a higher age (adjusted odds ratio 1.10 [95% CI 1.04–1.15], P = .001) and a higher ΔP (adjusted odds ratio 1.16 [95% CI 1.00–1.33], P = .046) (Table 2). NMB use was associated with decreased mortality (adjusted odds ratio 0.27 [95% CI 0.09–0.81], P = .02). Prone positioning, VT, compliance, and ventilatory ratio were not significantly associated with increased mortality. When ΔP was replaced by PEEP, PIP, or mechanical power, a higher age, no NMB use, a higher PIP, and a higher mechanical power were significantly associated with increased mortality (Table S3, see supplementary material related to this paper at http://rc.rcjournal.com).
Characteristics, Ventilation Parameters, and Outcomes According to NMB Use
The NMB group had shorter median time from hospital admission to initiation of mechanical ventilation and lower median serum creatinine than the no NMB group (Table S4, see supplementary material related to this paper at http://rc.rcjournal.com). All subjects in whom prone positioning was performed required NMB. During the first 4 d of mechanical ventilation, both groups received a median VT < 8.0 cm H2O, but the values were significantly lower in the NMB group (Fig. 3). The median PIP, ΔP, and mechanical power values were also lower in the NMB group, but the difference did not reach statistical significance. Among the clinical outcomes, the NMB group had lower rates of ICU, hospital, and 60-d mortalities; a higher rate of ventilator weaning; and a lower rate of tracheostomy (Table 3). The survival curves for the NMB and no NMB groups are shown in Figure 4. The median duration of NMB use was 5 d (IQR 4–9 d) and was significantly longer in nonsurvivors (Table 1). Subjects who received NMB for ≥ 6 d (vs < 6 d) had higher rate of ICU mortality, lower rate of ventilator weaning, lower ventilator-free d at day 28, and higher rate of tracheostomy (Table S5, see supplementary material related to this paper at http://rc.rcjournal.com).
Clinical Outcomes according to Prone Positioning and VT
Clinical outcomes did not differ when subjects were stratified according to prone positioning (yes vs no) (Table S6, see supplementary material related to this paper at http://rc.rcjournal.com) and VT (low [≤ 6 mL/kg IBW] vs intermediate [> 6 to ≤ 8 mL/kg IBW] vs high [> 8 mL/kg IBW]) (Table S7, see supplementary material related to this paper at http://rc.rcjournal.com).
Discussion
Although previous studies have evaluated the ventilation parameters and management strategies for subjects with COVID-19 receiving mechanical ventilation,8 there are limited data on their potential associations with mortality. This study revealed that higher ΔP, PIP, and mechanical power values during the first 4 d of mechanical ventilation were independently associated with ICU mortality. NMB use had a protective effect on mortality.
The 37% ICU mortality rate of subjects with COVID-19 requiring mechanical ventilation is in line with the 30–41% mortality rates reported in recent studies with similar subject characteristics.20-22 Meanwhile, the median hospital length of stay of 27 d was substantially longer than the median of 14 d in the LUNG SAFE study,23 which included subjects with non-COVID-19 ARDS. The ΔP is associated with mortality in patients with non-COVID-19 ARDS receiving mechanical ventilation.18 In our study, ΔP was similar to non-COVID-19 ARDS and was also associated with mortality in the multivariate analysis. The changes in ΔP were determined largely by PIP changes because the PEEP changes were minimal. Interestingly, ΔP and PIP were associated with mortality despite those values being within the limits of lung-protective ventilation in our subjects,11 even for nonsurvivors. Thus, our findings suggest that lower ΔP and PIP targets than recommended in the current guidelines might be beneficial in patients with COVID-19 ARDS. Mechanical power, which has been suggested as a risk factor for mortality in patients with ARDS,24 was very high (median 26.0 J/min) and was associated with mortality.
Meanwhile, higher PEEP values were not associated with reduced mortality in the multivariate analysis. In ARDS, high PEEP may improve oxygenation via alveolar recruitment, although this finding can be inconsistent in COVID-19 ARDS.25,26 The increases in PEEP levels could only improve survival through its effects on ΔP.18 Moreover, inappropriate PEEP may result in lung overdistention, increase in dead space, and intrapulmonary shunting.27 Some investigators have suggested that patients with COVID-19 ARDS with high compliance may be ventilated with higher VT (7–9 mL/kg IBW) and lower PEEP (< 10 cm H2O) than recommended.28 The median PEEP values in survivors on day 3 and day 4 of mechanical ventilation were < 10 cm H2O. However, the benefit of a lower PEEP in COVID-19 ARDS should be interpreted with caution because our subjects generally had low compliance.
The median compliance on day 1 of mechanical ventilation was 30.4 mL/cm H2O, similar to findings by recent studies on subjects with COVID-19 and non-COVID-19 ARDS.20,23,29 These observations contrast with the concept of a relatively high compliance specific to COVID-19 ARDS.6,7 The finding that compliance values were not associated with mortality is also consistent with previous studies in COVID-19 and non-COVID-19 ARDS.14,20,29,30 Ventilatory ratio is a surrogate marker for dead-space ventilation and is associated with an increased risk of mortality in patients with ARDS.17 A single-center study demonstrated that ventilatory ratio was associated with hospital mortality, after adjustment for different ventilation parameters, in patients with COVID-19 ARDS.31 However, a secondary analysis from the multi-center cohort revealed that ventilatory ratio could be a marker of disease severity but not an independent predictor of mortality.32 Our findings also showed that ventilatory ratio increased with ARDS severity but was not associated with ICU mortality.
In mechanically ventilated patients, NMB is used to optimize endotracheal intubation, facilitate mechanical ventilation, and improve safety of prone positioning.33 However, data regarding the efficacy and indications of NMB in patients with COVID-19 are limited. The 76% rate of NMB use in our study falls within the 72–88% rate reported in recent studies.20,22,29,34 These rates are strikingly higher than in non-COVID-19 ARDS.23 Unexpectedly, we found better clinical outcomes of mortality, ventilator weaning, and tracheostomy rates in the NMB group. There are several explanations for this phenomenon. First, patients with COVID-19 ARDS appear to have a high respiratory drive, which may lead to self-inflicted lung injury.35 NMB may have been used frequently to maintain lower VT and transpulmonary pressure. Indeed, the median VT, PIP, ΔP, and mechanical power values in the current study were lower in the NMB group. Second, the use of NMB may be associated with a lower risk of patient-ventilator asynchrony.36 Third, the energy requirements of patients with COVID-19 receiving mechanical ventilation are considerably high during the early period of ICU admission, and NMB may significantly decrease total energy expenditure.37 Fourth, due to the unusual pandemic situation, there may have been inadequate staffing and training to provide lung-protective ventilation without paralytics and daily monitoring of the benefit of NMB. This may be supported by the findings that the duration of NMB use (median 5 d) was longer than recommended, and protective ventilation was only marginally achieved in the no NMB group. Nonetheless, the protective effects of NMB in ARDS have been inconsistent,9,38 and other studies showed worse outcomes in patients with COVID-19 who received NMB.36,39 In addition, subjects who received prolonged NMB in the current study had worse clinical outcomes, perhaps due to disease progression or the occurrence of ICU-acquired weakness. Thus, our results should be considered as hypothesis generating and need to be confirmed in larger prospective studies.
In total, 16% of the subjects was placed in prone positioning, and this rate is relatively lower than the 61–76% rate in previous studies20,22,29 but comparable to that in the LUNG SAFE study.23 Resource limitations, expertise, and practice differences could explain these discrepancies. Studies evaluating the efficacy of prone positioning in patients with COVID-19 are lacking, although a recent study showed that early prone positioning reduced the risk of death in patients with COVID-19 on mechanical ventilation.40 Our study also showed lower ICU mortality and higher ventilator weaning rates in the prone group. However, these differences did not reach statistical significance. The timing of mechanical ventilation was not associated with mortality. One possible explanation is that the median time from hospital admission to initiation of mechanical ventilation was only 1 d. Further studies are required regarding the optimal timing of mechanical ventilation in COVID-19 ARDS.
The main strength of our study was the study design that assessed potential associations of various ventilation parameters and strategies with the most clinically relevant outcome of mortality. Unlike previous studies that only considered baseline ventilation parameters, serial changes in parameters during the first 4 d of mechanical ventilation may provide further information regarding the progress of COVID-19 ARDS. In addition, the daily mean values of mechanical ventilation settings were collected instead of point values to account for dynamic characteristic of mechanical ventilation.
However, the present study also had some limitations. First, residual and unmeasured confounding may have been included because of the retrospective design, and these may have biased the results. Moreover, a nonrandomized design precludes any inference of causality on associations between ΔP values and NMB use and mortality. Second, the relatively small sample size might have limited the power to detect the significant effects of prone positioning, VT, and PaO2/FIO2. Third, data from only the first 4 d of mechanical ventilation were included in the multivariate analysis; and thus, we could not determine whether ventilation parameters after day 4 of mechanical ventilation would affect mortality. However, mechanical ventilation parameters during this period are crucial influencing factors of the prognosis of patients with ARDS.19,41 Fourth, the mechanical ventilation settings and management were not standardized across centers; and the indications and timings of rescue therapies were unclear, which might have led to selection bias. However, the collected ventilation parameters demonstrated high similarity and adherence to current ARDS guidelines. Patient enrollment and clinical care might have differed between the centers, although KOVIDS provided a standardized case record form and encouraged the investigators to follow the guidelines for critically ill patients with COVID-19. Fifth, some of the subjects were admitted in the early pandemic period, and overwhelmed critical care services may have biased the results. Finally, other factors such as superinfection may have influenced mortality, although the difference was not significant between survivors and nonsurvivors. Larger randomized controlled trials may address these limitations.
Conclusions
In conclusion, higher ΔP, PIP, and mechanical power values during mechanical ventilation were independently associated with ICU mortality in subjects with COVID-19 ARDS. The short-term use of NMB was associated with decreased mortality in these subjects. Further studies are needed to establish the feasibility and efficacy of lung-protective ventilation using NMB during periods of limited resources.
Footnotes
- Correspondence: Won-Young Kim MD PhD, Department of Internal Medicine, Chung-Ang University Hospital, Chung-Ang University College of Medicine, 102 Heukseok-ro, Dongjak-gu, Seoul 06973, Republic of Korea. E-mail: wykim81{at}cau.ac.kr
The authors have disclosed no conflicts of interest.
This research was supported by a National Research Foundation of Korea grant funded by the Korea government (MSIP) (2019R1G1A1099373).
Supplementary material related to this paper is available at http://rc.rcjournal.com.
- Copyright © 2022 by Daedalus Enterprises