key: cord-0764594-hp5sa483
authors: Chomistek, Andrea K.; Liang, Caihua; Doherty, Michael C.; Clifford, C. Robin; Ogilvie, Rachel P.; Gately, Robert V.; Song, Jennifer N.; Enger, Cheryl; Lin, Nancy D.; Wang, Florence T.; Seeger, John D.
title: Predictors of critical care, mechanical ventilation, and mortality among hospitalized patients with COVID-19 in an electronic health record database
date: 2022-04-29
journal: BMC Infect Dis
DOI: 10.1186/s12879-022-07383-6
sha: 16a47ecac91cc4a023397135485604aac198d13a
doc_id: 764594
cord_uid: hp5sa483

BACKGROUND: There are limited data on risk factors for serious outcomes and death from COVID-19 among patients representative of the U.S. population. The objective of this study was to determine risk factors for critical care, ventilation, and death among hospitalized patients with COVID-19. METHODS: This was a cohort study using data from Optum’s longitudinal COVID-19 electronic health record database derived from a network of healthcare provider organizations across the US. The study included patients with confirmed COVID-19 (presence of ICD-10-CM code U07.1 and/or positive SARS-CoV-2 test) between January 2020 and November 2020. Patient characteristics and clinical variables at start of hospitalization were evaluated for their association with subsequent serious outcomes (critical care, mechanical ventilation, and death) using odds ratios (OR) and 95% confidence intervals (CI) from logistic regression, adjusted for demographic variables. RESULTS: Among 56,996 hospitalized COVID-19 patients (49.5% male and 72.4% ≥ 50 years), 11,967 received critical care, 9136 received mechanical ventilation, and 8526 died. The median duration of hospitalization was 6 days (IQR: 4, 11), and this was longer among patients that experienced an outcome: 11 days (IQR: 6, 19) for critical care, 15 days (IQR: 8, 24) for mechanical ventilation, and 10 days (IQR: 5, 17) for death. Dyspnea and hypoxemia were the most prevalent symptoms and both were associated with serious outcomes in adjusted models. Additionally, temperature, C-reactive protein, ferritin, lactate dehydrogenase, D-dimer, and oxygen saturation measured during hospitalization were predictors of serious outcomes as were several in-hospital diagnoses. The strongest associations were observed for acute respiratory failure (critical care: OR, 6.30; 95% CI, 5.99–6.63; ventilation: OR, 8.55; 95% CI, 8.02–9.11; death: OR, 3.36; 95% CI, 3.17–3.55) and sepsis (critical care: OR, 4.59; 95% CI, 4.39–4.81; ventilation: OR, 5.26; 95% CI, 5.00–5.53; death: OR, 4.14; 95% CI, 3.92–4.38). Treatment with angiotensin-converting enzyme inhibitors/angiotensin receptor blockers during hospitalization were inversely associated with death (OR, 0.57; 95% CI, 0.54–0.61). CONCLUSIONS: We identified several clinical characteristics associated with receipt of critical care, mechanical ventilation, and death among COVID-19 patients. Future studies into the mechanisms that lead to severe COVID-19 disease are warranted. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12879-022-07383-6.

Since coronavirus disease 2019 (COVID-19) first emerged in China in December 2019, there have been over 270 million confirmed cases and 5.3 million deaths worldwide due to COVID-19 as of 15 December 2021, according to the World Health Organization [1] . The United States (US) leads the world in cases (49.8 million) and deaths (792,371) from COVID-19 and these numbers are expected to continue to rise through the start of 2022.

There is significant heterogeneity in the clinical presentation of COVID-19 infection, ranging from patients who are asymptomatic to those with severe disease [2] [3] [4] . It is important to determine predictors of serious outcomes as patients may decline rapidly after initially presenting with mild symptoms [5] . Identifying predictors of serious outcomes may enable clinicians to deliver appropriate care to patients early as well as inform interventions to reduce risk of death [6] .

The serious outcomes of COVID-19 (e.g., intensive care unit admission, receipt of mechanical ventilation, death) and their preceding risk factors have been identified previously [7] [8] [9] [10] . Common factors associated with progression to serious disease include age, male sex, obesity, and comorbid diseases, including diabetes and renal disease. Additionally, it has been recognized that biomarkers, such as C-reactive protein (CRP) and D-dimer, may be associated with serious outcomes. However, studies from early in the pandemic were small and sought to identify the strongest predictors of serious disease and death from a broad set of variables. Further, some of these studies were hospital-based case series and may not be representative of all patients hospitalized with COVID-19 in the United States.

The purpose of this study was to apply an exploratory, data-driven approach to the identification of potential risk factors for serious outcomes among patients with COVID-19 in order to inform clinicians and researchers of characteristics that may be integral to identifying high risk patients. Thus, the objective was to determine demographic and clinical predictors associated with critical care, mechanical ventilation, and death among hospitalized COVID-19 patients in a large electronic health record (EHR) database that is representative of a geographically diverse U.S. population.

This was a retrospective cohort study that included patients confirmed with COVID-19 infection between January 2020 and November 2020. Among these patients, a subset of patients hospitalized with COVID-19 was identified from those with an inpatient health care encounter. Confirmation of COVID-19 infection was based on presence of a specific International Classification of Diseases, 10 th Revision, Clinical Modification (ICD-10-CM) diagnosis code (U07.1) and/or a positive severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral test. The date of confirmed infection was the earlier of the date of diagnosis or the date of a positive test. The cohort entry date was the earliest date that the patient met both of the following criteria: confirmed COVID-19 infection and admission to the hospital. For patients who were hospitalized prior to contracting COVID-19, the date of cohort entry was the date of confirmed infection. For patients who were confirmed to have COVID-19 before they were admitted to the hospital, the date of cohort entry was date of hospital admission. This approach for assigning cohort entry date allows for the description of clinical characteristics at the time patients were first hospitalized with COVID-19. For clinical characteristics other than death, patients were followed from cohort entry to discharge or 30 days after cohort entry, whichever came first. Deaths were identified in all follow-up available, including during and after hospitalization.

Patients were identified from Optum's longitudinal COVID-19 EHR database derived from a network of healthcare provider organizations across the U.S. The COVID-19 EHR database consists of a subset of patients from Optum's EHR database, which represents a geographically diverse U.S. patient population with more than 85 million patients from 2007 through 2019. Optum's EHR database includes data collected from tens of thousands of providers and hundreds of hospitals representing more than 60 electronic medical record (EMR)-based provider/hospital networks across the U.S. This database incorporates clinical and medical administrative data from both inpatient and ambulatory Conclusions: We identified several clinical characteristics associated with receipt of critical care, mechanical ventilation, and death among COVID-19 patients. Future studies into the mechanisms that lead to severe COVID-19 disease are warranted. Keywords: Coronavirus, SARS-CoV-2, Risk factors, Death, Biomarkers, Treatments EMRs, practice management systems, and numerous other internal systems. Information is processed from across the continuum of care, including acute inpatient stays and outpatient visits. The COVID-19 data captures diagnostics specific to the COVID-19 patient during initial presentation at hospital admission, acute illness, and convalescence with over 500 mapped labs and bedside observations, including COVID-19 specific testing. The data are incorporated into the underlying database on a biweekly basis, allowing for near real-time analysis and assessment of the COVID-19 clinical landscape. The database is certified as de-identified by an independent statistical expert following Health Insurance Portability and Accountability Act statistical de-identification rules.

Demographic characteristics were assessed on the date of cohort entry. Comorbidities and medication use were assessed in the 21 days prior to cohort entry. Comorbidities were identified by ICD-10-CM diagnosis code with a diagnosis status of "history of " and medications were mapped according to the Anatomical Therapeutic Chemical (ATC) classification scheme.

Additionally, vital signs, laboratory results, symptoms, diagnoses, and treatments during hospitalization were assessed. For vital signs and laboratory results, the first measurement on or after the date of cohort entry was selected if the patient had multiple measurements. Symptoms and diagnoses were identified by ICD-10-CM diagnosis codes.

The primary outcomes of interest were receipt of critical care, mechanical ventilation, and death. Patients were classified according to the presence or absence of an outcome, and outcome groups were not mutually exclusive.

Receipt of critical care was identified using Current Procedural Terminology, 4th Edition (CPT-4) codes. Receipt of mechanical ventilation was identified by CPT-4 and International Classification of Diseases, 10th Revision, Procedural Coding System (ICD-10-PCS) codes and included intubation, mechanical ventilation, and extracorporeal membrane oxygenation (ECMO). Receipt of critical care and mechanical ventilation were ascertained during hospitalization.

Death was ascertained via linkage to the Social Security Administration's Death Master File and/or the presence of a death indicator in the structured EHR data within all data available during and after hospitalization.

All analyses were conducted using SAS version 9.4 (SAS Institute Inc, Cary, NC). Baseline characteristics were examined overall and according to outcome (critical care, mechanical ventilation, and death). Categorical variables were summarized using frequency and percent while continuous variables were summarized using median and interquartile range (IQR).

For the association analyses, vital signs and laboratory values were transformed into categorical variables. Dichotomous variables were created based on clinicallyrelevant cutpoints. Additionally, categories based on quintiles were generated to examine the shape of the dose-response relationship of vital signs and laboratory values with each outcome. Quintiles were determined based on the distribution of each biomarker among hospitalized patients overall. Tests for linear trend were computed by using the medians of the quintiles modeled as a continuous variable.

Logistic regression models were used to estimate unadjusted and adjusted odds ratios (OR) and corresponding 95% confidence intervals (CI) for associations between the covariates and each outcome. Adjusted models included age, gender, region, race, and week of cohort entry. Week of cohort entry was included as a covariate to adjust for any potential changes in patient characteristics or treatments over time.

A total of 56,996 hospitalized patients with COVID-19 between January 2020 and November 2020 were identified ( Fig. 1 ). Of these, 11,967 (21.0%) received critical care, 9136 (16.0%) received mechanical ventilation, Statins and angiotensin-converting enzyme inhibitors (ACEs)/angiotensin receptor blockers (ARBs) were the most prevalent patient-reported medications (26.2% and 20.2%, respectively, overall). Table 2 presents the distributions of vital signs, laboratory values, symptoms, diagnoses, and treatments received during hospitalization among COVID-19 patients. The median duration of hospitalization overall was 6 days (IQR: 4, 11). The duration was longer among patients that experienced an outcome: 11 days (IQR: 6, 19) for critical care, 15 days (IQR: 8, 24) for mechanical ventilation, and 10 days (IQR: 5, 17) for death. Among hospitalized patients overall, 6.6% had a temperature > 38 degrees Celsius and 10.1% had an oxygen saturation < 90%. Markers of inflammation and coagulation were elevated for many patients, including 85.7% of patients with C-reactive protein (CRP) > 10 mg/L and 76.8% of patients with D-dimer > 250 ng/mL (DDU).

The most common symptoms during hospitalization among COVID-19 patients were dyspnea (28.0%) and hypoxemia (25.8%) ( Table 2 ). Prevalence of these symptoms was higher among patients who received critical care (37.4% and 33.2%, respectively), those who received ventilation (41.4% and 35.4%, respectively), and those who died (34.8% and 31.9%, respectively). Relatedly, the most prevalent diagnoses during hospitalization among patients overall were pneumonia (59.6%) and acute respiratory failure (44.8%). These diagnoses were even more common among those with an outcome, with the highest prevalence observed among patients that received mechanical ventilation (acute respiratory failure, 85.7%; pneumonia, 84.7%). Anticoagulants were the most prevalent treatment, received by 79.8% of hospitalized patients with COVID-19 overall during their hospitalization. Among other treatments, 15.5% of patients received chloroquine or hydroxychloroquine, 29.3% received dexamethasone, and 22.8% received remdesivir. Figure 2 presents the associations between baseline characteristics and receipt of critical care, ventilation, and death adjusted for age, gender, region, race, and week of cohort entry. Unadjusted ORs for the associations between covariates and outcomes are provided in Additional file 1: Tables S1-S3. Age was associated with all 3 outcomes, particularly death; the OR for patients ≥ 80 years of age compared to those 50-59 years of age was 7.61 (95% CI, 6.96-8.32) ( Fig. 2A) . Female patients were less likely to experience any of the 3 outcomes compared to males. Regarding race, hospitalized patients that received critical care, ventilation, or died were less likely to be African American compared to Caucasian (critical care: OR, 0.82; 95% CI, 0.77-0.87; mechanical ventilation: OR, 0.94; 95% CI, 0.89-1.00; death: OR, 0.83; 95% CI, 0.77-0.88) (Fig. 2B) .

After adjusting for demographic variables, several comorbidities at baseline were associated with serious outcomes (Fig. 2C ). Comorbidities associated with higher odds of all 3 outcomes included diabetes, obesity, chronic obstructive pulmonary disease, coronary artery disease, congestive heart failure, kidney disease, and liver disease. Asthma, hypertension, and cancer showed varying associations with each outcome. In the unadjusted model, hypertension was positively associated with death (OR, 1.37; 95% CI, 1.30-1.43; Additional file 1: Table S3 ). However, once adjusted for demographic variables, there was an inverse association between hypertension and death (OR, 0.85; 95% CI, 0.81-0.89). Asthma was associated with higher odds of critical care (OR, 1.21; 95% CI, 1.12-1.30) and ventilation (OR, 1.33; 95% CI, 1.23-1.44), but lower odds of death (OR, 0.87; 95% CI, 0.79-0.97).

Among select patient-reported medications at baseline, use of statins, corticosteroids, and PPIs was positively associated with receipt of critical care and ventilation, but not death (Fig. 2D) . Statins, ACEs/ARBs, and NSAIDs showed inverse associations with death. Figure 3 presents the adjusted ORs for the associations between clinical characteristics during hospitalization and receipt of critical care, ventilation, and death among hospitalized patients with COVID-19. Measurements of temperature, CRP, ferritin, lactate dehydrogenase (LDH), and D-dimer that exceeded the clinically-relevant cutpoint were significantly associated with an increased risk of all 3 serious outcomes in adjusted models (Fig. 3A) . Likewise, measurements of oxygen saturation and platelets that were less than the clinically-relevant cutpoint were also significantly associated with an increased risk of all 3 serious outcomes. D-dimer and LDH measured during hospitalization had the highest adjusted ORs for the association with death: OR, 2.95 (95% CI, 2.59-3.35) for D-dimer > 250 ng/mL (DDU) and OR, 2.81 (95% CI, 2.61-3.02) for LDH > 280 U/L. Fibrinogen > 400 mg/dL was associated with a lower risk of all 3 outcomes.

The associations between each biomarker in quintiles and serious outcomes are shown in Fig. 4 . For CRP, D-dimer, ferritin, and LDH, the relationships appeared linear for all 3 outcomes (p values for linear trend < 0.0001, Additional file 1: Tables S1-S3). For temperature, oxygen saturation, and platelets, the relationships were less linear, although many of the p values for linear trend were < 0.0001. For fibrinogen, the for ventilation, and 6.07 (95% CI, 5.21-7.08) for death (Additional file 1: Tables S1-S3).

Dyspnea and hypoxemia were positively associated with all 3 serious outcomes in adjusted models, with both showing the strongest association with receipt of ventilation (dyspnea: OR, 1.68; 95% CI, 1.59-1.76; hypoxemia: OR, 1.44; 95% CI, 1.37-1.52) (Fig. 3B ). In contrast, cough was associated with lower odds of all 3 outcomes (critical care: OR, 0.87; 95% CI, 0.81-0.92; mechanical ventilation: OR, 0.80; 95% CI, 0.74-0.85; death: OR, 0.70; 95% CAD, coronary artery disease; CHF, congestive heart failure; ACE, angiotensin-converting enzyme; ARB, angiotensin II receptor blocker; NSAIDS, non-steroidal anti-inflammatory drugs; PPIs, proton-pump inhibitors. Logistic regression models were used to estimate odds ratios (OR) and corresponding 95% confidence intervals (CI) adjusted for age, gender, region, race, and week of cohort entry CI, 0.65-0.76). Patients with nausea and vomiting as well as malaise also had lower odds of death. All selected in-hospital diagnoses showed positive associations with each of the 3 serious outcomes (Fig. 3C) 

With the exception ACEs/ARBs, all of the selected treatments were associated with higher odds of serious outcomes (Fig. 3D) . The highest ORs were observed for the associations between treatments and receipt of ventilation, ranging from 1.87 (95% CI, 1.68-2.08) for antivirals to 4.02 (95% CI, 3.71-4.36) for immunosuppressants. Treatment with ACEs/ARBs was inversely associated with receipt of critical care (OR, 0.84; 95% CI, 0.80-0.89) and death (OR, 0.57; 95% CI, 0.54-0.61).

In this study, we identified and described patients who experienced a serious outcome (critical care, mechanical ventilation, and death) among 56,996 hospitalized patients with COVID-19 within a large, EHR database. We conducted an evaluation of the association of many , oxygen saturation; PLT, platelet count; CRP, C-reactive protein; LDH, lactate dehydrogenase; MI, myocardial infarction; ACE, angiotensin-converting enzyme inhibitors; ARB, angiotensin II receptor blocker Logistic regression models were used to estimate odds ratios (OR) and corresponding 95% confidence intervals (CI) adjusted for age, gender, region, race, and week of cohort entry demographic and clinical characteristics with these outcomes in order to identify potential signals for experiencing a serious outcome. As was observed in prior studies [7] [8] [9] [10] , older age and male gender were associated with higher risk of serious outcomes, along with comorbidities such as diabetes and kidney disease. We also observed associations with clinical characteristics measured during hospitalization, including several laboratory markers, symptoms, and diagnoses.

In this study of hospitalized patients with COVID-19, we found that African-Americans and women were at lower risk of experiencing a serious outcome compared to Caucasians and men, respectively. Since the start of the pandemic, African-Americans have been more likely than Caucasians to contract and be hospitalized with COVID-19 [11] . Nonetheless, evidence suggests that once hospitalized, they do not have a higher risk of adverse outcomes [10, [12] [13] [14] . In contrast, there does not appear to be sex difference in number of confirmed cases of COVID-19, but the death rate has been higher in men than women [15] . The reason for the sex difference in rate of mortality among patients with COVID-19 remains unknown, but it has been suggested that the mechanism involves a combination of biological and psychosocial factors [16] .

Hypertension was the most prevalent comorbidity among hospitalized COVID-19 patients in this study and its prevalence was even higher among patients with a serious outcome. Hypertension was positively associated with receipt of critical care and mechanical ventilation, but inversely associated with mortality after adjusting for age and other demographics. Findings in this study are consistent with previous studies that have found hypertension to be common among adults diagnosed with COVID-19, but not associated with mortality after adjusting for covariates [14, 17] . Thus, while there may be overrepresentation of hypertension among adults with COVID-19, it appears this association may be confounded by age and other covariates, and possibly affected by treatment.

In the current study, treatment with ACEs/ARBs during hospitalization was associated with lower odds of critical care and mortality. It is recommended that patients who are prescribed ACEs/ARBs for cardiovascular disease continue taking these medications if hospitalized with COVID-19 [18, 19] . With the exception of ACEs/ARBs, we observed that many treatments received during hospitalization were associated with higher odds of receipt of critical care, mechanical ventilation, or death. An explanation for this finding may be that most medications, particularly those that were investigational, were only recommended for use among patients with severe disease. For example, remdesivir and dexamethasone are recommended for hospitalized COVID-19 patients that require supplemental oxygen [19] . As of 15 December 2021, remdesivir is the only FDA-approved drug for the treatment of COVID-19, although emergency use authorizations have been issued for multiple anti-SARS-CoV-2 monoclonal antibody products (i.e., bamlanivimab plus etesevimab, casirivimab plus imdevimab, and sotrovimab) [19] . Other medications are under investigation.

Several vital signs and laboratory results were associated with serious outcomes in COVID-19 patients in this study, with the strongest associations observed for LDH, D-dimer, and CRP. LDH is an enzyme found within cells in almost all organ systems [20] . Its levels rise when the body's tissues are damaged. Recent studies have found that high levels of LDH may be predictive of COVID-19 severity and death [17, 20] . Similarly, increased D-dimer, an indicator of coagulopathy, has been linked to higher risk of mortality in COVID-19 patients [8, 21] . CRP is an inflammatory marker and has been shown to be elevated in patients with severe COVID-19 disease [9, 12, 22] . The findings in the current study are consistent with these smaller studies and, taken together, suggest these laboratory markers, if measured soon after admission, may help clinicians triage patients who may be at higher risk of progression to severe disease.

The current study was based on an analysis of EHR data, which are valuable for the examination of clinical health outcomes and treatment patterns. Nonetheless, all EHR databases have inherent limitations because the data are collected for the purpose of clinical patient management, not research. Unlike in clinical trials, where the collection of clinical and laboratory measures is standardized, the Optum EHR includes real-world clinical data obtained from multiple medical and laboratory settings used for patient care. Because data are not collected in a systematic way, clinical measurements (e.g., vital signs and laboratory results) were not available for all patients. Additionally, the presence of a diagnosis code in EHR data may not represent the actual presence of disease, as the diagnosis code may be incorrectly coded, or included as rule-out criteria rather than actual disease. We assumed the absence of a diagnosis code meant the patient did not have the disease. This assumption may be a reason why we observed an inverse association between cough and serious outcomes; in a patient who is very ill, severe symptoms like dyspnea and hypoxemia may be more likely to be recorded than minor symptoms such as cough. Furthermore, it is possible that some comorbidities and medications may not have been captured as health care encounters with medical providers who do not contract with Optum would not be observed.

The median duration of hospitalization among patients who died was 10 days. However, deaths were identified within all data available, not only during hospitalization. As such, it is possible that duration of hospitalization may have been shorter among patients that died during hospitalization. An additional limitation of the EHR database is the data lag at the time of data extraction, which likely resulted in an underestimation of the number of deaths. Finally, residual confounding is a concern as we only adjusted for demographic covariates.

In summary, we utilized an exploratory, data-driven approach to identify many clinical characteristics that were associated with receipt of critical care, mechanical ventilation, and death among patients hospitalized with COVID-19. As more continues to be learned about COVID-19 by clinicians and researchers, future studies should move toward causal inference and focus on identifying the etiologic factors and mechanisms responsible for some patients experiencing more severe COVID-19 disease.

COVID-19) Dashboard

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

Clinical characteristics of coronavirus disease 2019 in China

Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City Area

Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan, China

Precision medicine and public health: new challenges for effective and sustainable health

Development and validation of a clinical risk score to predict the occurrence of critical illness in hospitalized patients with COVID-19

convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year • At BMC, research is always in progress. Learn more biomedcentral.com/submissions Ready to submit your research Ready to submit your research ? Choose BMC

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

Characteristics and predictors of death among 4035 consecutively hospitalized patients with COVID-19 in Spain

Risk factors for intensive care unit admission and in-hospital mortality among hospitalized adults identified through the U.S. coronavirus disease 2019 (COVID-19)-associated hospitalization surveillance network (COVID-NET)

COVID-19 Hospitalization and Death by Race/Ethnicity

Patient trajectories among persons hospitalized for COVID-19: a cohort study

Hospitalization and mortality among black patients and white patients with Covid-19

Clinical characteristics and morbidity associated with coronavirus disease 2019 in a series of patients in metropolitan detroit

The COVID-19 Sex-Disaggregated Data Tracker

Men and COVID-19: a biopsychosocial approach to understanding sex differences in mortality and recommendations for practice and policy interventions

Prediction model and risk scores of ICU admission and mortality in COVID-19

Continuing versus suspending ACE inhibitors and ARBs: impact of adverse outcomes in hospitalized patients with COVID-19-the BRACE CORONA Trial

COVID-19) Treatment Guidelines

Lactate dehydrogenase levels predict coronavirus disease 2019 (COVID-19) severity and mortality: a pooled analysis

Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia

COVID-19-associated hyperinflammation and escalation of patient care: a retrospective longitudinal cohort study

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

The authors would like to acknowledge Katherine Reed for her technical assistance in preparing the manuscript.

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12879-022-07383-6.

Author contributions AKC, CL, NDL, FTW, JDS contributed to the study concept and design. NDL, FTW, and JDS acquired the data. AKC and CL drafted the manuscript. MCD, CRC, RPO, RVG, JNS, CE, NDL, FTW, and JDS critically revised the manuscript for important intellectual content. All authors contributed to the analysis and interpretation of data. All authors read and approved the final manuscript.

This work was funded by Optum Epidemiology.

The de-identified database used for the current study is not publicly available, but is available from Optum through a data license agreement. More information can be found at the following website: https:// www. optum. com/ busin ess/ solut ions/ life-scien ces/ real-world-data/ ehr-data. html.

Ethics approval and consent to participate Not applicable because this study utilized a commercially available, de-identified database. No administrative permissions were required for this study as this database has been certified as de-identified by an independent statistical expert following Health Insurance Portability and Accountability Act statistical de-identification rules and managed according to Optum ® customer data use agreements. As no study team members had access to patient identifiers linked to the data, review by an ethics committee or institutional review board was not required, nor was consent to participate.

Not applicable.