https://orcid.org/0000-0003-1020-0818
Figures
Abstract
Background
Coronavirus disease 2019 (COVID-19) patients reportedly have high bicarbonate concentration. However, its relationship to the disease progression are obscure.
Methods
In this two-center retrospective study, we included COVID-19 patients with moderate severity between March 2020 and May 2021. We classified patients into three groups according to bicarbonate concentrations: high (>27 mEq/L), normal (21 to 27 mEq/L), and low (<21 mEq/L). The primary outcome was the time to clinical worsening defined by the requirement of intubation or death during 90 days. We evaluated high or low bicarbonate concentration during the clinical course related to the primary outcome using multivariable Cox proportional hazard models.
Results
Of the 60 participants (median age 72 years), 60% were men. Participants were classified into high (13 patients), normal (30 patients), and low (17 patients) groups. Clinical worsening occurred in 54% of patients in the high group, 23% in the normal group, and 65% in the low group. Both high and low groups were associated with a higher clinical worsening rate: HR, 3.02 (95% CI, 1.05 to 8.63) in the high group; 3.49 (95% CI: 1.33 to 9.12) in the low group.
Citation: Sada K-e, Yamamoto R, Yano A, Miyauchi A, Kawamura M, Ito H (2022) Bicarbonate concentration as a predictor of prognosis in moderately severe COVID-19 patients: A multicenter retrospective study. PLoS ONE 17(6): e0270141. https://doi.org/10.1371/journal.pone.0270141
Editor: Tatsuo Shimosawa, International University of Health and Welfare, School of Medicine, JAPAN
Received: February 5, 2022; Accepted: June 3, 2022; Published: June 24, 2022
Copyright: © 2022 Sada et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Data cannot be shared publicly because the data contain potentially identifying or sensitive patient information. However, data are available from the Ethics Committee of Kochi Medical School (contact via is21@kochi-u.ac.jp) for researchers who meet the criteria for access to confidential data.
Funding: The authors received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Since December 2019, coronavirus disease 2019 (COVID-19) caused by SARS- CoV-2 has rapidly spread throughout the world and is still life-threatening despite the hard work of researchers.
High bicarbonate concentration was reportedly noted in COVID-19 patients with metabolic alkalosis [1]. Other studies have reported that hypokalemia is common among patients with COVID-19 [2, 3]. High bicarbonate concentration, metabolic alkalosis, and hypokalemia in COVID-19 patients are assumed to be due to the activation of the renin-angiotensin-aldosterone (RAA) system via the downregulation of angiotensin-converting enzyme 2 (ACE2) by SARS-CoV-2 [4]. There is a difference in the reports of two observational studies in terms of the prevalence of metabolic alkalosis and hypokalemia [1, 2]. Because patients with severe COVID-19 often exhibit metabolic acidosis and hyperkalemia resulting from multiple organ damage, the prevalence of high bicarbonate concentration, metabolic alkalosis, and hypokalemia might change according to COVID-19 severity.
The association between high bicarbonate concentration and disease progression is not fully elucidated. One report showed no difference in mortality between patients with and without metabolic alkalosis [1]. However, metabolic acidosis and low bicarbonate concentration are caused by multiple organ failure and related to high mortality in severe COVID-19 patients, so patients with metabolic acidosis and low bicarbonate concentration should be evaluated separately [5–8]. In addition, one case report showed that metabolic alkalosis developed with the worsening of COVID-19 severity in a patient. Thus, it is difficult to cross-sectionally evaluate the association between metabolic alkalosis and mortality without considering patients’ clinical courses. One report evaluated the temporal change in bicarbonate concentration [9]. No difference in bicarbonate concentration at baseline was found between survivors and non-survivors, and a lower bicarbonate concentration was found at the last visit in non-survivors. Because the sequential time course of bicarbonate concentration was not evaluated and the disease severity at baseline was not taken into consideration in that study, the importance of a high bicarbonate concentration remains unclear.
The aim of our study was to investigate the sequential time course of bicarbonate levels with respect to changes in disease severity in COVID-19 patients with moderate disease severity.
Materials and methods
Study design and setting
We conducted a two-center retrospective cohort study at Kochi Health Sciences Center and Kochi Prefectural Hata-Kenmin Hospital. This study was performed following the Strengthening the Reporting of Observational studies in Epidemiology (STROBE) guidelines for reporting [10].
Study population
We included patients with COVID-19 infection who were admitted to Kochi Health Sciences Center and Kochi Prefectural Hata-Kenmin Hospital from March 1, 2020, to May 31, 2021. All patients were over the age of 18 years and had SARS-CoV-2 infection confirmed by either reverse transcriptase-polymerase chain reaction or antigen test on respiratory tract samples and were started on oxygen therapy (moderate severity). The day when oxygen therapy was started was considered as the first day of inclusion. Patients were excluded from the study if any of the following were applicable on the day of inclusion: the need for maintenance dialysis due to end-stage kidney disease; blood gases were not obtained; intubated at other hospitals.
Data collection
From a review of electronic health records, we collected data such as age, sex, body mass index (BMI), comorbidity, medication history, vital sign, laboratory test, treatment limitations (limitations in providing life-sustaining therapies such as mechanical ventilation, cardiopulmonary resuscitation, and extracorporeal membranous oxygenation), and sequential organ failure assessment (SOFA) score. We collected serum pH, bicarbonate, and potassium concentrations that were obtained from venous or arterial blood gas analysis. Recent studies have suggested that the relationship between pH and PaCO2 concentration obtained with venous blood gas and arterial blood gas sampling could allow venous blood to be used instead of arterial blood in analyses [11]. The PaCO2 concentrations were collected only from arterial blood gas data. All blood gas analyses were performed using either an “ABL800 FLEX” or “ABL90FLEX” (Radiometer Medical ApS, Copenhagen, Denmark). Venous or arterial pH and PaCO2 concentration were measured using selective electrodes. Bicarbonate concentration was calculated from pH and PaCO2 using the Henderson-Hasselbalch equation. Potassium concentration was determined with direct potentiometry using ion-selective electrodes.
Bicarbonate concentration categories
The patients were classified into three categories according to the bicarbonate concentration based on previous literature: high (>27 mEq/L), normal (21 to 27 mEq/L), and low (<21 mEq/L) [12–15]. Patients with bicarbonate concentration exceeding 27 mEq/L at least once within 7 days after inclusion were assigned to the high bicarbonate group while patients with bicarbonate concentration that decreased below 21 mEq/L at least once were assigned to the low bicarbonate group.
Outcome measurement
The primary outcome was time to clinical worsening defined by the requirement of intubation or death within 90 days. Secondary outcomes included time to intubation, time to death, and time to clinical worsening within 28 days. Observations were censored if the event of interest did not occur, i.e., the patient was discharged from the hospital or did not visit the outpatient clinic after discharge.
Statistical methods
Patient characteristics were described as the mean (standard deviation [SD]) or median (interquartile range [IQR]) as appropriate. Survival curves were plotted using the Kaplan-Meier method to compare the event probability at different points of time and to compare the three bicarbonate concentration groups. A log-rank test was applied to find the statistical significance among the three groups. As the primary analysis, Cox proportional hazard models were used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for the association between bicarbonate concentration categories and clinical worsening within 90 days; the incidence of clinical worsening in the normal bicarbonate group was compared with those in the high bicarbonate and low bicarbonate groups. Multivariable analysis was performed after adjustment for treatment limitation.
We applied the same analyses as for the primary outcome to assess the association among three bicarbonate concentration groups and the following secondary outcomes: intubation within 90 days, 90-day mortality, and clinical worsening within 28 days.
To perform sensitivity analyses, we tested several Cox proportional hazard models to assess the robustness of the primary analysis. First, we adjusted the following covariates in the four models: age, sex, treatment limitation, and respiratory sub-SOFA score in model 1; age, sex, treatment limitation, respiratory sub-SOFA score, and renal sub-SOFA score in model 2; age, sex, treatment limitation, and PaO2/FiO2 ratio in model 3; age, sex, treatment limitation, respiratory sub-SOFA score, chronic pulmonary disease, hypertension, diabetes, and BMI in model 4. Second, we used the Cox regression model using penalized spline to evaluate the association between bicarbonate concentration as continuous value and clinical worsening within 90 days. GAM is an extension of the generalized linear model, where the predictors are related to the outcome via a smooth, possibly non-linear, function [16]. We adjusted with treatment limitation with a complete case analysis.
A p-value <0.05 was considered to indicate statistical significance. The analyses were performed using R software, version 4.0.3 (The R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/).
Compliance with ethical standards
This study was conducted in accordance with the Declaration of Helsinki and the Ethical Guidelines for Medical and Health Research Involving Human Subjects in Japan. This study was approved by the Ethics Committee of Kochi Medical School (2021–007), and the Ethics Committee of Kochi Prefectural Hata-Kenmin Hospital. Patient data were anonymized and de-identified before the analysis. According to the Ethical Guidelines for Medical and Health Research Involving Human Subjects in Japan, the need to obtain written informed consent was waived due to the retrospective nature of the study.
Results
Patient characteristics
From March 1, 2020, to May 31, 2021, 112 patients with COVID-19 were admitted to the hospitals and 70 patients met the inclusion criteria. Of these, 10 patients were excluded. Finally, 60 patients were included in the analysis (Fig 1).
ESKD, end-stage kidney disease; HD, hemodialysis.
The median age of the included patients was 72 (IQR [64, 78]) years and 60% were males. The median BMI was 25.0 kg/m2 (IQR [22.3, 28.7]). All patients were Japanese and 57% of patients had hypertension. The patients were classified into high bicarbonate (13 patients), normal bicarbonate (30 patients), and low bicarbonate (17 patients) groups according to bicarbonate concentrations. The median day from disease onset to inclusion was 8 (IQR [6, 9]) days. The median total SOFA score was 3 (IQR [2, 4]) and the median respiratory sub-SOFA score was 2 (IQR [2, 2]). Other patient characteristics at inclusion were summarized in Table 1.
Blood gas data
The trend of pH, bicarbonate, PaCO2, and potassium concentrations among the three bicarbonate concentration groups were shown in Fig 2. There was a tendency for the pH to decrease in the high bicarbonate group after day 3, but there was no significant change in the pH trend among the three groups (Fig 2a). The mean (±SD) bicarbonate concentration on days 1, 3, and 7 were 25.1 ± 3.0 mEq/L, 29.0 ± 5.2, and 33.5 ± 4.5 in the high bicarbonate group while 19.4 ± 3.8 mEq/L, 21.1 ± 1.6 mEq/L, and 21.9 ± 3.4 mEq/L in the low bicarbonate group (Fig 2b). PaCO2 tended to increase in the high bicarbonate group (Fig 2c). After day 3, there was a slight but not significant tendency for potassium to increase in the low bicarbonate group (Fig 2d).
Association between bicarbonate concentration and the primary outcome
Clinical worsening occurred in 54% (7/13) of patients in the high bicarbonate group, 23% (7/30) in the normal bicarbonate group, and 65% (11/17) in the low bicarbonate group, respectively (Table 2). Fig 3 depicts the time to clinical worsening among the three bicarbonate concentration groups. The median times to clinical worsening among the three groups were significantly different; 11 (IQR [6, 41]) days in the high bicarbonate group, 16 (IQR [11, 25]) days in the normal bicarbonate, and 11 (IQR [7, 13]) days in the low bicarbonate, respectively; log-rank test, p = 0.012). The incidence of death and mechanical ventilation was higher in the high and low bicarbonate groups compared to the normal bicarbonate group (Table 2).
High bicarbonate (>27 mEq/L); normal bicarbonate (21 to 27 mEq/L); low bicarbonate (<21 mEq/L).
In the unadjusted analysis, high bicarbonate concentration and low bicarbonate concentration were associated with increased clinical worsening within 90 days (unadjusted HR 2.98 [95% CI 1.04 to 8.53], p = 0.04 in the high bicarbonate group, unadjusted HR: 3.80 [95% CI 1.46 to 9.89], p = 0.006 in the low bicarbonate group). After adjusting for treatment limitation, this association remained (adjusted HR: 3.02 [95% CI 1.05 to 8.63], p = 0.04 in the high bicarbonate group; adjusted HR: 3.49 [95% CI 1.33 to 9.12], p = 0.01 in the low bicarbonate group; Table 2).
Secondary outcomes
The high and low bicarbonate groups had higher intubation within 90 days compared to the normal group (Table 2). For death within 90 days, the low bicarbonate group was higher than the normal group. When comparing the high bicarbonate group with the normal bicarbonate group, there was a trend toward more deaths in the high bicarbonate group, but the difference was not significant (Table 2). Compared with the normal bicarbonate group, clinical worsening within 28 days was significantly higher in the high bicarbonate group and the low bicarbonate group (Table 2).
Sensitivity analysis
The association between the bicarbonate categories and clinical worsening within 90 days remained similar in various multivariable Cox regression analyses (Table 3). We also performed a sensitivity analysis using continuous serum bicarbonate concentration as the exposure. Multivariable-Cox regression model using penalized spline, with the same adjustments as in the primary model, showed that both high and low bicarbonate concentrations were associated with an increased risk of clinical worsening within 90 days (Fig 4).
The black solid line indicates log hazard ratio (HR) and the dotted lines indicate standard error (SE).
Discussion
In our study, we evaluated the time course of bicarbonate concentrations in COVID-19 patients with moderate severity. High bicarbonate concentrations were noted in 19% of the included patients during the 7 days after inclusion. Of the patients who experienced high bicarbonate concentrations, 54% developed severe conditions.
Bicarbonate concentrations in blood might increase before the development of severe disease status. A recent report showed metabolic alkalosis and high bicarbonate levels were common among COVID-19 patients [1]. It is speculated that the downregulation of ACE2 by the activation of the RAA system leads to high bicarbonate concentrations and metabolic alkalosis. However, one other report found no increase in bicarbonate concentrations in COVID-19 patients [9]. One possible reason for this discrepancy is the disease severity of the enrolled patients. Although the former report did not show the outcomes in detail, the latter report included 23% of non-survivors. The patients with severe COVID-19 often exhibited metabolic acidosis and low bicarbonate concentrations caused by multiple organ damage, which is probably the reason why the bicarbonate concentrations were not high in the latter study. There was also a big difference in the prevalence of hypokalemia between the two aforementioned studies. One study reported hypokalemia in only 9% of the patients, [1] whereas 55% of the patients in the other study had hypokalemia [2]. These differences might be related not only to disease severity but also to disease time courses. In our study, potassium concentrations gradually decreased in the high bicarbonate concentration group up to day 5. Therefore, the prevalence of hypokalemia may change according to the timing of the evaluation.
The patients with high bicarbonate concentrations showed worse prognosis compared to those with low bicarbonate concentrations. There are two possible reasons why elevated bicarbonate is associated with poor prognosis. First, elevated bicarbonate levels may be an early predictor of respiratory acidosis. Increased PaCO2 and respiratory acidosis contribute to the severity of COVID-19. The PaCO2 increased on Day 4 in the high bicarbonate group but only one patient met the definition of respiratory acidosis (pH<7.38 and PaCO2>42mmHg). Higher bicarbonate may predict a subsequent rise in PaCO2 and respiratory acidosis. Another reason is that elevated bicarbonate levels may reflect overactivity of the RAA system in patients with worsening COVID-19. Our results showed that PaCO2 levels compensatorily increased without lowering pH in the high bicarbonate concentration group, which may indicate that high bicarbonate concentrations were related to metabolic alkalosis. The previous study [1] found no difference in mortality between COVID-19 patients with and without metabolic alkalosis, but the study did not separate the patients with metabolic acidosis. Metabolic acidosis is a well-known risk factor for mortality in COVID-19 patients [5]. In our study, we could confirm the worst prognosis in the patients with low bicarbonate concentrations, which indicated that the patients exhibited metabolic acidosis. In the course of the activation of the RAA system, angiotensin II, via activation of angiotensin type 1a receptor, reportedly promotes inflammatory responses in COVID-19 patients [17]. Therefore, the high bicarbonate concentrations noted in our study population possibly reflects the hyperactivation of the RAA system in patients with worsening COVID-19.
Our patient had a slightly higher pH on Day 1 (Fig 2a), but whether this was due to disease specificity is not known because it was not compared to non-COVID-19 patients in our study. Although evidence is limited, higher pH has also been reported in patients with H1N1 influenza [18, 19]. Even before hypotension occurs, fever and sepsis can lead to an increase in pH via hyperventilation [20]. Therefore, the high pH on Day 1 may not be attributed to the disease specificity of COVID-19.
There are some limitations to our study. First, bicarbonate concentrations were evaluated using both venous and arterial blood analysis. The utility of venous blood gas should be evaluated further although it is a less invasive procedure than arterial blood gas analysis. However, previous meta-analysis studies indicated that venous and arterial bicarbonate levels are reasonably close [11, 21], so we believe that venous blood gas analysis could also be useful for bicarbonate monitoring in COVID-19 patients. Second, the bicarbonate concentration obtained at baseline may not be a useful predictor of prognosis. Since a few patients showed high bicarbonate concentrations at baseline, we cannot exclude the requirement for sequential monitoring of patients with normal bicarbonate concentrations at baseline. The third was the use of a Cox model adjusted for six covariates for fewer outcomes. However, the standard errors of the regression coefficients were small, and we believe that we have achieved a stable estimation.
Conclusions
High bicarbonate concentrations during the clinical course in COVID-19 patients with moderate disease status were related to a worse prognosis. Sequential monitoring of bicarbonate concentrations may be useful to predict the prognosis of COVID-19 patients.
Acknowledgments
The authors would like to thank Mrs. Kimiko Shirotake for her assistance in data management and all medical staff members who cooperated in data collection.
References
- 1. Rood J, Davids R, Le Roux A, Du Plessis M, Parker A, Allwood BW, et al. Metabolic alkalosis in hospitalised COVID-19 patients: A window to the pathogenesis? S Afr Med J. 2020;110(11):13109. pmid:33403973
- 2. Chen D, Li X, Song Q, Hu C, Su F, Dai J, et al. Assessment of Hypokalemia and Clinical Characteristics in Patients With Coronavirus Disease 2019 in Wenzhou, China. JAMA Netw Open. 2020;3(6):e2011122. pmid:32525548
- 3. Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med. 2020;382(18):1708–20. pmid:32109013
- 4. Wiese OJ, Allwood BW, Zemlin AE. COVID-19 and the renin-angiotensin system (RAS): A spark that sets the forest alight? Med Hypotheses. 2020;144:110231. pmid:33254538
- 5. Choron RL, Butts CA, Bargoud C, Krumrei NJ, Teichman AL, Schroeder ME, et al. Fever in the ICU: A Predictor of Mortality in Mechanically Ventilated COVID-19 Patients. J Intensive Care Med. 2021;36(4):484–93. pmid:33317374
- 6. Kafan S, Tadbir Vajargah K, Sheikhvatan M, Tabrizi G, Salimzadeh A, Montazeri M, et al. Predicting Risk Score for Mechanical Ventilation in Hospitalized Adult Patients Suffering from COVID-19. Anesth Pain Med. 2021;11(2):e112424. pmid:34336617
- 7. Türkay Kunt A, Kozaci N, Torun E. Mortality Predictors in Patients Diagnosed with COVID-19 in the Emergency Department: ECG, Laboratory and CT. Medicina (Kaunas). 2021;57(6). pmid:34204209
- 8. Asghar MS, Haider Kazmi SJ, Khan NA, Akram M, Hassan M, Rasheed U, et al. Poor Prognostic Biochemical Markers Predicting Fatalities Caused by COVID-19: A Retrospective Observational Study From a Developing Country. Cureus. 2020;12(8):e9575. pmid:32913691
- 9. Ouyang SM, Zhu HQ, Xie YN, Zou ZS, Zuo HM, Rao YW, et al. Temporal changes in laboratory markers of survivors and non-survivors of adult inpatients with COVID-19. BMC Infect Dis. 2020;20(1):952. pmid:33308159
- 10. Von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: Guidelines for Reporting Observational Studies. Annals of Internal Medicine. 2007;147(8):573. pmid:17938396
- 11. Byrne AL, Bennett M, Chatterji R, Symons R, Pace NL, Thomas PS. Peripheral venous and arterial blood gas analysis in adults: are they comparable? A systematic review and meta-analysis. Respirology. 2014;19(2):168–75. pmid:24383789
- 12.
Theodore AC. Arterial blood gases. In: TW P, editor. UpToDate. UpToDate, Waltham, MA. (Accessed on March 29, 2022.)2022.
- 13.
Michael Emmett BFP. Simple and mixed acid-base disorders. In: Post T, editor. UpToDate. UpToDate, Waltham, MA, (Accessed on March 29, 2022).
- 14.
Longenecker JC, Nelson TR. High-yield Acid-base: Lippincott Williams & Wilkins; 2007.
- 15.
Adrogué H, Madias N. Acid-Base Disorders and Their Treatment. 2005.
- 16.
Wood SN. Generalized additive models: an introduction with R: CRC press; 2017.
- 17. Magrone T, Magrone M, Jirillo E. Focus on Receptors for Coronaviruses with Special Reference to Angiotensin- Converting Enzyme 2 as a Potential Drug Target—A Perspective. Endocr Metab Immune Disord Drug Targets. 2020;20(6):807–11. pmid:32338224
- 18. Courouble P, Geukens P, Laarbaui F, Beauloye C, Van Caenegem O, Jacquet LM. Adult respiratory distress syndrome caused by 2009 H1N1 influenza during pregnancy: success of ECMO for both the mother and the child. J Extra Corpor Technol. 2011;43(2):75–8. pmid:21848176
- 19. Yassin Z, Saadat M, Abtahi H, Foroushani AR, Peiman S. Prognostic value of on admission arterial PCO2 in hospitalized patients with community-acquired pneumonia. Journal of Thoracic Disease. 2016;8(10):2765–71. pmid:27867552
- 20. Elisaf M, Theodorou J, Pappas H, Siamopoulos KC. Acid-base and electrolyte abnormalities in febrile patients with bacteraemia. Eur J Med. 1993;2(7):404–7. pmid:8258028
- 21. Bloom BM, Grundlingh J, Bestwick JP, Harris T. The role of venous blood gas in the emergency department: a systematic review and meta-analysis. Eur J Emerg Med. 2014;21(2):81–8. pmid:23903783