Abstract
Objectives We aimed to describe trends in the incidence of adverse outcomes among patients who tested positive for SARS-CoV-2 between February and September 2020 within a national healthcare system.
Setting US Veterans Affairs national healthcare system.
Participants Enrollees in the VA healthcare system who tested positive for SARS-CoV-2 between 2/28/2020 and 9/30/2020 (n=55,952).
Outcomes Death, hospitalization, intensive care unit (ICU) admission and mechanical ventilation within 30 days of testing positive.The incidence of these outcomes was examined among patients infected each month and trends were evaluated using an interrupted time-series analysis.
Results Between February and July 2020, during the first wave of the US pandemic, there were marked downward trends in the 30-day incidence of hospitalization (44.2% to 15.8%), ICU admission (20.3% to 5.3%), mechanical ventilation (12.7% to 2.2%), and death (12.5% to 4.4%), with subsequent stabilization between July and September 2020. These trends persisted after adjustment for sociodemographic characteristics, comorbid conditions, and documented symptoms and after additional adjustment for laboratory test results among hospitalized patients, including among subgroups admitted to the ICU and treated with mechanical ventilation. Among hospitalized patients, use of hydroxychloroquine (56.5% to 0%), azithromycin (48.3% to 16.6%) vasopressors (20.6% to 8.7%), and dialysis (11.6% to 3.8%) decreased while use of dexamethasone (3.4% to 53.1%), other corticosteroids (4.9% to 29.0%) and remdesivir (1.7% to 45.4%) increased from February to September.
Conclusions Among patients who tested positive for SARS-CoV-2 in a large national US healthcare system, risk for a range of adverse outcomes decreased markedly between February and July, with subsequent stabilization from July to September. These trends were not explained by changes in measured baseline patient characteristics.
Introduction
Infection with SARS-CoV-2 results in a broad spectrum of clinical severity ranging from asymptomatic infection to life-threatening illness.1, 2 Among individuals who test positive for SARS-CoV-2, there is substantial variation across different populations in rates of hospital admission (from 8% to 80.7%),3-5 mechanical ventilation (from 2.3% of the Chinese population to 93.2% of critically ill admitted to New York area hospitals6-18), and mortality (from 2.8% of all patients in the US to 10.2-67% of hospitalized and high-risk patients).3, 7, 8, 12, 17-23 To date, only a few studies have described trends in these measures over time within populations.24-26
Since the first case of SARS-CoV-2 was diagnosed in the US on January 19, 2020, the pandemic has spread across the country from early epicenters in the Pacific Northwest, New York, New Orleans and Detroit. During this time, there have been marked changes in the availability and approach to testing, in population screening strategies, in mitigation strategies both in the community and within health systems, in how symptomatic patients are managed medically and most recently in the availability of vaccination. However, even prior to the availability of a vaccine and as the number deaths due to COVID-19 across the world continued to climb27, several studies described improved outcomes among infected patients after the initial months of the pandemic. Improvements in mortality were reported for patients treated in an 8-hospital health system in Houston, Texas between mid-March and mid-July, in a 3-hospital system in New York between March and August, 202026, and in a single hospital in Milano, Italy between late February and mid-May, 2020.24, 28 Similar improvements in outcomes over time have also been described for the most severe cases, including among patients admitted to the ICU in England, Wales and Northern Ireland between February and May, 2020 and among patients admitted to the ICU in England from March through June, 2020.25, 29 To our knowledge, no prior studies have described temporal trends in the range of severe complications of COVID-19 including hospital admission, intensive care unit (ICU) admission, mechanical ventilation, dialysis and death in a national US health care system and extending beyond the first wave of the pandemic i.e. after June 2020. The Department of Veterans Affairs (VA) healthcare system supports the largest, integrated, comprehensive national healthcare system in the US, providing a unique view of national trends in outcomes among those infected with SARS-CoV-2. We used national data from the VA to measure trends in hospitalization, ICU admission, mechanical ventilation, dialysis, medication use and death among individuals who tested positive for SARS-CoV-2 between February 28 and September 30, 2020, with follow-up extending through November 19, 2020.
Methods
Data source and study population
The VA provides care for more than six million Veterans annually and uses a single, national comprehensive electronic healthcare information network. We used data from the VA’s Corporate Data Warehouse (CDW), a relational database of VA enrollees’ electronic health records (EHR), developed by the VA Informatics and Computing Infrastructure (VINCI) to support research and clinical operations, including datasets and analytic variables specifically related to SARS-CoV-2 developed by VINCI as part of the “COVID-19 Shared Data Resource.”30.
This study was approved by the Institutional Review Board of the Veterans Affairs Puget Sound Healthcare System, which granted a waiver of informed consent.
Study Population, definition of SARS-CoV-2 infection and index date
We identified all VA enrollees (n=559,616) with available results of testing for SARS-CoV-2 using approved polymerase chain reaction (PCR) tests from 02/28/2020 (the date of the first SARS-CoV-2 test in the VA system) to 9/30/2020, including 55,952 (10.0%) who tested positive at least once. Final adjudication of SARS-CoV-2 status was performed by the VA National Surveillance Tool, which is intended to be the single, authoritative data source for determination of positive and negative cases within VA. The index date was the earliest date of a documented positive test or, in rare instances, if the patient was already hospitalized prior to testing positive, the date of hospitalization.
Adverse outcomes
We ascertained the following four outcomes within 30 days of the index date: 1) hospitalization, 2) ICU admission, 3) mechanical ventilation, and 4) all-cause mortality. Deaths that occurred both within and outside the VA are comprehensively captured in CDW through a variety of sources including VA inpatient files, VA Beneficiary Identification and Records Locator System (BIRLS), Social Security Administration (SSA) death files, and the Department of Defense.31 However, episodes of mechanical ventilation, ICU admission and hospitalization that occurred outside the VA are captured only if these were paid for by the VA and thus do not include episodes of care covered by Medicare, Medicaid, the Department of Defense or private health insurance. Patients were followed through 11/19/2020 for study outcomes allowing for a minimum follow-up time of 50 days from the most recent index (09/30/21), which was adequate for outcomes to be electronically recorded in CDW.
Baseline characteristics
Baseline sociodemographic characteristics, previously associated with adverse SARS-CoV-2 outcomes,32 included age, sex, race, ethnicity, body mass index (BMI), urban vs. rural residence (based on zip codes), and geographic location (divided into 10 US Federal Regions33) (Table 2).
Comorbid conditions were based on ICD-10 codes recorded in VA electronic health records during the 2-year period before the index date.30 We used the Deyo modification34 of the Charlson Comorbidity Index (CCI)35 to estimate overall burden of comorbidity (Table 2). We also extracted the following comorbid conditions: diabetes, chronic kidney disease, hypertension, obstructive sleep apnea, and obesity hypoventilation, which were associated with adverse outcomes among SARS-CoV-2 infected Veterans.32
We identified symptoms potentially associated with SARS-CoV-2 documented on or within 30 days prior to the index date based on a combination of natural language processing to capture relevant documentation in the EHR and searching for relevant ICD-10 codes,30 occurring on or within 30 days prior to the index date. In our multivariable models, we used only two symptoms (dyspnea and fever).
Finally, we extracted 5 routine laboratory blood tests (albumin, aspartate aminotransferase, creatinine, white blood cell count, neutrophil-to-lymphocyte ratio), which we recently showed to be associated with adverse outcomes in SARS-CoV-232. For each test, we extracted the value closest to the index date, on or within 10 days before the index date, or if absent, within 5 days after the index date (93% of tests were performed within 2 days of the index date).
Pharmacotherapy, Noninvasive Ventilation and Dialysis
We identified medications prescribed for SARS-CoV-2 on or within 60 days after the index date. We also report use of non-invasive ventilation (continuous (CPAP) or bilevel (BiPAP) positive airway pressure) and dialysis during the same time frame. Use of high-flow oxygen supplementation could not be reliably captured in available data sources for this project.
Statistical Analysis
We divided the observation period into monthly intervals based on the index date. Using the Kaplan-Meier method, we calculated for the persons infected each month the 30-day cumulative incidence of hospitalization, ICU admission, mechanical ventilation, and mortality from the index date. Patients who did not experience the relevant outcome of interest were censored at the time of death or after 30 days of follow-up. We used Cox proportional hazards regression to compare the risk of each outcome among patients infected during each subsequent month relative to the earliest time period after adjusting for the following baseline characteristics, which we found in a prior study to be associated with adverse outcomes among SARS-CoV-2-infected VA enrollees32: age, sex, race, ethnicity, BMI, regional COVID-19 disease burden, CCI, diabetes, chronic kidney disease, hypertension, obstructive sleep apnea, obesity hypoventilation, fever, and dyspnea. In analyses limited to hospitalized patients, we additionally adjusted for the 5 laboratory tests listed above that have also been associated with adverse outcomes32. Cox proportional hazards analyses were stratified by the VA medical center where patients were tested for SARS-CoV-2. The aim of these multivariable adjustments was to attempt to account for changes over time in baseline characteristics of SARS-CoV-2 positive patients (resulting, for example, from increased availability of testing) that could have resulted in changes in the risk of adverse outcomes. Interrupted time series analysis of trends methods were used to test monthly trends in outcomes over the testing period from February through September 202036. Time trends were modeled using ordinary least squares with adjustment for an auto-correlation lag of one month. Robust variances estimators were used to address heteroskedasticity.
Results
Trends in baseline characteristicsof patients with SARS-CoV-2 infection
Among all 55,952 individuals who tested positive for SARS-CoV-2 during the study period, the majority were male (89.3%), mean age was 60.0 years, 58.2% were White, 31.0% were Black, 12.3% were Hispanic, 12.4% had a CCI≥5, 25.4% had documentation of fever and 11% had documentation of dyspnea (Table 1).
Baseline characteristics of VA patients who tested positive for SARS-CoV-2 between 2/28/20 and 9/30/20, by month
The number of patients testing positive each month ranged from 13,967 in the earliest time period (2/28/20 to 03/31/20) to 124,771 in July. There were substantial trends in certain baseline socio-demographic characteristics of patients who tested positive for SARS-CoV-2 across successive months, including a decrease over time in the proportion who were Black (from 51.9% before March 31 to 21.3% in September), an increase in the proportion who were white (from 51.9% to 68.2%) and an increase in the proportion who resided in rural/highly rural (vs. urban) areas (from 18.8% to 51.9% in September). The proportion of patients originating from each U.S. geographical region varied dramatically over time consistent with well-described changes in disease epicenters (e.g 24.2% of positive patients from region 2 (NY, NJ, PR) in February/March versus only 2.8% in September).
The proportion of patients with severe comorbidity burden (CCI≥5) decreased over time (from 18.1% to 11%) as did the proportion with documented fever (from 57.7% to 20%) or dyspnea (from 26.6% to 8.9%). Among patients who were hospitalized, there were no clear temporal trends the proportion who had abnormalities in the five laboratory tests examined.
Trends in pharmacotherapy, dialysis and non-invasive ventilation in patients with SARS-CoV-2 infection
The proportion of hospitalized patients who received hydroxychloroquine on or after the index date decreased from 56.6% in February/March to 2.4% in May 2020 and ranged from 0-0.6% thereafter (Table 2, Figure 1a). The proportion of hospitalized patients prescribed azithromycin decreased from 48.3% in February/March to 16.6% in September and use of vasopressors decreased from 20.6% to 8.7%. Tocilizumab prescription rates in hospitalized patients ranged from a peak of 8.7% in April to 0.8% in September 2020. In contrast, remdesivir prescriptions increased from 1.7% to 45.4% of hospitalized patients from February/March to September. Dexamethasone prescriptions were low during the first three months of the ascertainment period (<4% for hospitalized patients) but from May to September increased from 3% to 53.1% among hospitalized patients and prescription of corticosteroids other than dexamethasone also increased during the same time from 7.8% to 29%. Prescription of anticoagulants among hospitalized patients did not change appreciably over time. Use of investigational therapies including eculizumab and androgen deprivation therapy such as degarelix was minimal throughout the ascertainment period. The percentage of hospitalized patients who received dialysis decreased from a high of 11.6% in the earliest time period to a low of 3.8% in August, and use of non-invasive ventilation increased from a low of 8.4% in February/March to a high of 12.5% in August.
Trends over time in pharmacotherapy and adverse outcomes of SARS-CoV-2 infection in a national cohort of VA patients
1a. Trends over time in pharmacotherapy among SARS-CoV-2-infected hospitalized patients
(“Vasopressors” include Norepinephrine, Dobutamine, Vasopressin)
1b. Trends over time in 30-day adverse outcomes among all SARS-CoV-2-infected patients
1c. Trends over time in 30-day adverse outcomes among SARS-CoV-2-infected hospitalized patients
Changes in pharmacotherapy, noninvasive ventilation and dialysis over time in patients with SARS-CoV-2 infection
Analogous decreases in prescription of hydroxychloroquine, azithromycin and vasopressors and increases in prescription of remdesivir, dexamethasone and other corticosteroids were observed when presented as a proportion of all infected patients (rather than the hospitalized subgroup) – Table 2.
Trends in risk of adverse outcomes: hospitalization, ICU admission, mechanical ventilation, death
Between February and July 2020, there were downward trends in the 30-day incidence of hospitalization (44.2% to 15.8%), ICU admission (20.3% to 5.3%), mechanical ventilation (12.7% to 2.2%), and death (12.5% to 4.4%) among patients who tested positive for SARS-CoV-2, with relative stabilization from July through September, 2020 (15.8% to 17.1% for hospitalization; 5.3% to 6.6% for ICU admission; 2.2% to 2.2% for mechanical ventilation; and 4.4% to 4.9% for mortality) (Table 3, Figure 1b).
Association between time period and adverse outcomes in patients with SARS-CoV-2 infection
Among hospitalized patients, the crude 30-day incidence of ICU admission decreased from 47.6% in February/March to 34.8% in July then increased slightly between July and September (34.8% to 39.5%) (Table 3, Figure 1c). The 30-day incidence of mechanical ventilation decreased from 29.2% in February/March to 11.4% in June and ranged between 9.7 and 11.9% from July to September. The 30-day incidence of mortality decreased from 24.1% in February/March to 12.0% in July and ranged between 12.5% to 12.9% thereafter. Among those admitted to the ICU, the crude 30-day incidence of mechanical ventilation ranged from 57.0% in February/March to 22.0% in August, increasing slightly to 23.9% in September. The crude 30-day mortality among those admitted to the ICU decreased from 39.2% in February/March to 22.7% in May, ranging from 20.2% to 23.2% thereafter. Among patients who received mechanical ventilation, the crude 30-day incidence of death decreased from 58.3% in February/March to 31.3% in May, increasing thereafter to between 39.8% to 46.7%. After adjustment for sociodemographic characteristics, comorbid conditions, symptoms and baseline laboratory characteristics, risk for all adverse outcomes decreased most markedly during the early months of the pandemic with stabilization during late summer and early fall (p<0.05 for interrupted time series test for difference in 30-day outcome rate slope before and after July, 2020 for all outcomes except for mortality among mechanically ventilated patients) (Table 3 and Figure 2).
Discussion
In our study of 55,952 VA enrollees who tested positive for SARS-CoV-2 between February 28 and September 30, 2020, there were marked decreases in the 30-day incidence of hospitalization, ICU admission, mechanical ventilation and death during the first wave of the US pandemic from February/March through July, 2020 with stabilization during the late summer and early fall. Similar trends were noted in analyses adjusted for baseline sociodemographic characteristics, comorbid conditions, symptoms and laboratory tests (among hospitalized patients) and in all subgroups examined, including those admitted to the hospital, those admitted to the ICU and those who received mechanical ventilation
While there have been marked changes over time in the characteristics of Veterans infected with SARS-CoV-2 that could perhaps have explained our findings, trends in improved outcomes during the first wave of the US pandemic with subsequent stabilization persisted after adjustment for a wide range of measured sociodemographic characteristics and comorbid conditions. While increased availability of testing over time would be expected to identify a greater number of asymptomatic or mildly symptomatic individuals infected with SARS-CoV-2, improvements in adverse outcomes during the initial months of the pandemic persisted after stratification by illness severity and among the sickest cohort members (i.e., those who were hospitalized, admitted to the ICU or received mechanical ventilation). It is also noteworthy that among hospitalized patients, results were similar after adjustment for a range of laboratory markers for illness acuity.
While it is possible that temporal trends in the frequency of adverse outcomes over time might reflect geographic variation in unmeasured factors (e.g., treatment practices) as the pandemic moved across the country, multivariable models were stratified by VA facility and our results were robust to adjustment for geographic region. Another possibility is that our results reflect changing thresholds for hospitalization, ICU admission and mechanical ventilation over the course of the pandemic. While relaxation in criteria for hospitalization, ICU admission and mechanical ventilation over the pandemic would be expected to be associated with reductions in mortality among patients who were hospitalized, admitted to the ICU or mechanically ventilated, we saw a decrease rather than increase in rates of each of these outcomes over time.
Our findings do raise the question of whether marked improvements in outcomes among patients infected with SARS-CoV-2 during the first wave of the pandemic might reflect changes in the clinical approach to caring for members of this population.37 We did observe very substantial increases in the use of potentially effective pharmacotherapy (e.g. remdesivir and dexamethasone) and decreases in the use of ineffective (e.g. azithromycin) or potentially harmful (e.g. hydroxychloroquine) therapy toward the end of the first wave of the pandemic, although inflection points in trends in the use of each of these agents did not exactly parallel those for adverse outcomes. Perhaps most importantly, there were probably other shifts in clinical practice that occurred within the first months of the pandemic that could not be reliably captured in the data sources available to us. Such shifts include high flow nasal canula [HFNC], pronation and noninvasive positive-pressure ventilation (NIPPV) which were increasingly used to avoid or delay intubation and may have resulted in improved outcomes38.
Another possible explanation for declining rates of adverse outcomes over time is that the severity of disease caused by SARS-CoV-2 infection may have declined between February and June 2020. This might occur, for example, through reduction in inoculum related to physical distancing or wearing face masks, which became more common during that time period.39 Also, as an RNA virus, SARS-CoV-2 does undergo rapid mutation and these mutations could impact the pathogenicity and clinical characteristics of infection40, 41. It is now clear that a SARS-CoV-2 variant carrying the Spike protein amino acid change D614G became progressively more common after March 2020 throughout the world, including the United States, and was the most prevalent variant by July 202042. This strain has enhanced replication and is more infectious42, 43 with but may be associated with decreased disease severity. It is important to emphasize that such a reduction in disease severity caused by SARS-CoV-2 has not been proven, but it could explain declining trends in adverse outcomes over time that persisted after adjustment for baseline characteristics. Since the end of our study period, many viral variants have been described, such as the “United Kingdom variant” (B.1.1.7), the “South Africa variant” (B.1.351) and the “Brazil variant” (P.1), all of which have now been detected in the US. All of these variants have greater transmissibility, but it is not clear whether they are associated with worse disease severity or higher risk of death.
There are several limitations to our study. First, while use of a Veteran cohort provided a useful window on trends in complications of SARS-CoV-2 over time within a racially diverse national health care system, our cohort includes relatively few women which may limit the generalizability of our findings. We did not capture episodes of hospitalization, ICU admission or mechanical ventilation of VA patients in non-VA facilities (unless they were paid for by the VA), but there is no evidence that such non-VA care increased substantially over the course of our study. There may also have been a delay in capturing episodes of care occurring in the community under VA purchased care. However, this would not explain the temporal trends in adverse outcomes among hospitalized patients reported here. Finally, although we adjusted for a number of sociodemographic characteristics, comorbid conditions, symptoms and laboratory tests, some residual confounding by disease severity and other unmeasured factors likely persisted.
In conclusion, we identified marked improvements in outcomes of SARS-CoV-2 infection among Veterans infected during the first wave of the US pandemic with stabilization in late summer and early fall. These trends may reflect changing treatment practices, public health measures and viral pathogenicity. Studies should continue to track trends in adverse outcomes of SARS-CoV-2 infection.
Data Availability
Data are not available due to restrictions imposed by the Department of Veterans Affairs on access to VA data.
Acknowledgements
We acknowledge the VA Informatics and Computing Infrastructure (VINCI) group who worked tirelessly to create the “COVID 19 Shared Data Resource”.
Footnotes
Funding: The study was supported in part by the Department of Veterans Affairs, CSR&D grant COVID19-8900-11 to GNI.
Declaration of Personal Interests: None of the authors has any conflicts of interest to disclose.
Disclaimer: The contents do not represent the views of the U.S. Department of Veterans Affairs or the US Government.
Authors’ Contributions and Authorship Statement George Ioannou is the guarantor of this paper. All authors approved the final version of the manuscript.
Ioannou: Study concept and design, analysis of data, interpretation of results, drafting of manuscript, critical revision of manuscript, obtaining funding.
Locke: Study concept and design
Green: Extraction of data, creation of analytic variables, analysis of data, interpretation of results.
Berry: Study concept and design, analysis of data, interpretation of results, drafting of manuscript, critical revision of manuscript
O’Hare: Study concept and design, interpretation of results, critical revision of manuscript
Shah: Study concept and design, critical revision of manuscript
Crothers: Study concept and design, interpretation of results, critical revision of manuscript
Eastment: Study concept and design, interpretation of results, drafting of manuscript, critical revision of manuscript
Fan: Study concept and design, interpretation of results, critical revision of manuscript
Dominitz: Study concept and design, interpretation of results, drafting of manuscript, critical revision of manuscript