COVID-19 Outcomes and Sequencing of SARS-CoV-2 isolated from Veterans in New England

Introduction

The Coronavirus Disease 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has infected nearly 30 million people in the United States and caused over 500,000 deaths.¹ Disease severities and outcomes vary among individuals who become infected, and throughout the course of the pandemic, researchers have tried to identify patient characteristics and comorbidities that may predispose patients to having worse outcomes. It still remains unclear why some people have severe disease requiring intubation or death, while others have asymptomatic infections. Some factors that have been shown to impact disease severity and mortality include Charlson comorbidity index score, age, and body mass index (BMI).^2-6

The genomic sequence of the virus itself may also impact the severity of disease. Previous research has divided globally circulating SARS-CoV-2 viruses into 5 major clades based on an 11-nucleotide barcode.⁷ These clades and mutations within them could impact infectivity and virulence. Early during the pandemic course in the United States, viral isolates on the East Coast predominantly saw the Spike Protein D614G mutation.⁸ This mutation has been suggested to have higher transmissibility(Ref),⁹ supported by the higher upper airway viral loads and higher viral replication in lung epithelial cells.¹⁰At the same time, the D614G variant has also been found to be more susceptible to neutralization by antibodies generated by infection with wild-type virus.^11,12 With the increasing focus on possible mutations such as the “UK variant” (B.1.1.7) being more lethal,¹³ more research is required on how mutations affect patient outcomes.

Given the high variability and conflicting data in predicting who will have poor outcomes, it becomes important to assess specific populations to give providers the best clinical picture on their patients. Clinical outcomes of Veterans with COVID-19 in New England and respective genomic variants of SARS-CoV-2 have not been described. We aimed to (1) describe patient characteristics, comorbidities, and disease factors that impact patient outcomes and (2) present data on the genomes of the viruses in these patients who were infected. Factors impacting outcomes will inform triage and management algorithms in these patients.

Methods

Our cohort was composed of all Veterans receiving care at Veterans Administration centers who tested positive in six New England states (CT, MA, ME, NH, RI, VT) from April 8, 2020, to September 16, 2020. We manually reviewed charts and recorded demographics (age, gender, race, BMI, long-term care facility (LTC) status, state). We recorded history of comorbidities, including immunosuppression, dementia, diabetes mellitus, chronic kidney disease (CKD) stage 3 or higher, chronic liver disease, coronary artery disease (CAD), heart failure, atrial fibrillation (aFib), chronic obstructive pulmonary disease (COPD), asthma, and current or recent tobacco use (based on patients eligible for lung cancer screening).¹⁴ Protocols were approved by Veterans Administration (VA) Connecticut Institutional Review Board.

We recorded hospitalization status of patients, mortality, and whether patients required oxygen (O2) within 24 hours of admission. Lastly, we divided patients based on their peak severity status depending on oxygenation requirements, dividing them into 5 categories: 1) no O2 requirement, 2) 1-3 liters (L) by nasal cannula (NC), 3) 4-6 L NC, 4) >6 L O2 or non-invasive positive pressure ventilation (NIPPV), and 5) mechanical ventilation. We used STATA v16 for univariate and multivariate logistic regressions on finding predictors of the outcomes of interest.

For the genetic data, we sequenced 299 whole genomes (20x coverage, ≥70% of genome) using Illumina (n = 238) and Nanopore (n = 61 platforms). Using BWA-MEM version 0.7.15, we aligned reads to the Wuhan-Hu-1 reference genomes (GenBank MN908937.3). With iVar v1.2.1, we trimmed sequencing adaptors and primer sequences and called bases by simple majority (>50% frequency) at each site to generate consensus genomes. An ambiguous N was used when fewer than 10 reads were present at a site.

We aligned these consensus genomes with MAFFT,¹⁵ and masked problematic sites.^16,17 We then built a phylogenetic tree with IQTree¹⁷ using an HKY substitution model and 1000 bootstraps. We visualized this tree using the Python module baltic v0.1.5. We used Pangolin to assign lineages to our genomes.

Results

Of 274 Veterans, 92.7% were male, 83.2% white, and mean age was 63 years ± 17.6 years (Table 1). Over a third resided in a long-term care facility (LTC) (n=92), of which 76% were at the VA’s community living center. Over half the residents were from Massachusetts, and one third from Connecticut. Most common comorbidities were coronary artery disease (27%), diabetes (25%), and tobacco use (23%).

Overall, 20.8% of all patients required oxygen support above their baseline home oxygen dose throughout their hospitalization. Additionally, 11.7% of all patients required O2 within 24 hours of admission. In terms of peak severity, 78.5% required only room air, 10.9% required 1-3 L NC, 4.0% required 4-6 L NC, 3.6% required >6L or NIPPV, and 2.2% required mechanical intubation. The hospitalization rate was 28.8% (Figure 1), with the highest rates in people with COPD (50%), CAD (45.9%), and those over age 80 years (44.4%). The overall mortality rate was 10.6% (Figure 2), with the highest rates in people with dementia (33.3%), age >80 (28.9%), and those from LTC (23.9%).

• Download figure
• Open in new tab

Figure 1. Percent of patients hospitalized, based on patient demographics and comorbidities

• Download figure
• Open in new tab

Figure 2. Percent of patients who died, based on patient demographics and comorbidities

On univariate regression, significant predictors of hospitalization were age (Odds Ratio (OR): 1.06), non-white race (OR: 1.97), LTC residence (OR: 1.78), COPD (OR: 3.04), CAD (OR: 2.93), and aFib (OR: 2.05) (Table 2). Peak severity increased with age (OR: 1.06), LTC residence (OR: 2.8), CAD (OR: 2.1), aFib (OR: 2.97), and O2 requirement on admission (OR: 46.2). Higher mortality was seen with age (OR: 1.10), LTC residence (OR: 7.86), BMI<30 (OR: 5.80), dementia (OR: 7.23), aFib (OR: 3.11), hospitalization (OR: 2.58), and O2 requirement on admission (OR: 4.34).

On multivariate regression, significant predictors of hospitalization were age (OR: 1.05) and Non-white race (OR: 2.39) (Table 3). Peak severity also varied by age (OR: 1.07) and O2 requirement on admission (OR: 45.7). Mortality was predicted by age (OR: 1.06), dementia (OR:3.44), and O2 requirement on admission (OR: 6.74). This means that for every year increase in age, the odds of hospitalization increased by 5%, peak severity increased by 7%, and mortality increased by 6%.

For the genomes, we found that all but one is part of 2 main lineages: A (2.36%) and B (97.3%) (Figure 3). One genome is part of the N.1 lineage. Notably, the majority of genomes (97%) are part of the B1 lineage, which is defined by a D614G substitution.

• Download figure
• Open in new tab

Figure 3. Maximum Likelihood tree of genomes

Discussion

Our study found that in an older cohort of Veterans from the six New England states with a high comorbidity burden, age was the single strongest predictor of hospitalization, peak severity, and mortality. Our study is unique in that it is the first to look at several COVID-19 outcomes in this population of Veterans. The CDC currently has a list of chronic medical conditions that predispose individuals to severe illness from COVID-19 infection.¹⁸ However, the list is so extensive that one study found over 75% of United States adult would fall under a high-risk category.¹⁹ Thus, understanding clinical factors that impact outcomes in Veterans in the Northeast will help clinicians in risk stratification for patients with similar demographic profiles.

Our primary outcome is that age is associated with worse outcomes and has been supported by nearly all studies on COVID-19 outcomes.^20-24 In our study, the impact of age was so large that most comorbidities became insignificant after adjusting for it. Interestingly, long term care status predicted all three of our outcomes on univariate analysis, but became insignificant after adjusting for age on multivariate analysis. Early during the COVID-19 era, residents of nursing homes had higher rates of infection as well as severe illness and mortality.²⁵ Our study shows that for Veterans in a long term care facility, a majority of which was the VA’s community living center, their disease outcomes were not impacted by their residence status.

Concurrent work in our lab suggests that oxygen requirement on admission predicts poor outcomes in Veterans. Similarly, our study found that patients who required oxygen within 24 hours of admission had worse outcomes. This is an important finding because other ways of using oxygenation status, such as the arterial oxygen saturation/ fraction of inspired oxygen (SaO2/FiO2) ratio,²⁶ can frequently change and thus become difficult for clinicians to use in practice. Knowing that patients who require oxygen within 24 hours of admission can help providers more easily triage who will likely end up requiring a higher level of care and are at higher risk of death throughout the admission.

Our finding that non-White participants were more likely to be hospitalized but experienced similar peak severity and mortality outcomes was repeated by several studies across different hospital systems.^27,28 One study of 5902 patients with COVID-19 in New York found that non-Hispanic Black patients were more likely to test positive for coronavirus but had similar outcomes.²⁹ Similar to our study, Azar et al. found that in 1052 cases of COVID in California, non-Hispanic African American patients had 2.7 increase in odds of being hospitalized.³⁰ Cardemil et al. also found a 4.6 and 4.2-fold increase in rates of hospitalization for Hispanic and Black veterans, respectively.²⁸ Many studies have shown that minorities often have delays in care due to mistrust of the medical system because of providers’ implicit biases,^31,32 and our results may indicate that non-white patients had more advanced or severe illness at the time of presentation. It is critical to continue ongoing efforts to combat medical inequities and target prevention efforts and education to communities and racial groups most affected by COVID-19.

Even after adjusting for age and other comorbidities, we found that patients with dementia had a higher risk of death. This is similar to other studies on patients with COVID-19 and dementia.^24,33,34 In a meta-analysis of 46,391 patients, Hariyanto et al. found that dementia was significantly associated with a higher risk of severe infection (relative risk: 2.63) and death (relative risk: 2.62).³⁵ This may be explained by a host of biological factors,³⁵ but it emphasizes the importance of extra care and monitoring required when approached with a patient with dementia. In this vulnerable population that may not be able to voice their symptoms, providers and caregivers must ensure they receive sufficient support throughout the course of their illness.

Limitations of this work include the smaller sample size. We are working on a larger study based off of these results by increasing our date ranges and capturing more patients. Furthermore, our study is specific to Veterans, which is a largely male and older cohort, and may not be generalizable to women. Our study time period was prior to established medical therapies for COVID and our outcomes reported are likely worse than you would expect today. Since this was during the first wave of the pandemic when standard guidelines were not yet solidified, many patients not requiring oxygen were admitted which is not standard of practice today. Strengths of our study include its comprehensive scope, wide geographic range, manual chart review allowing for the capturing of all comorbidities and oxygenation parameters that may not be available otherwise in a database, and multivariate analysis of many potential risk factors.

Conclusion

Our study found that in an older cohort of Veterans from the six New England states with a high comorbidity burden, age was the single strongest predictor of hospitalization, peak severity, and mortality. Non-white Veterans were more likely to be hospitalized, and patients who required oxygen on admission were more likely to have severe disease and higher rates of mortality. Furthermore, patients with dementia were more likely to die. Multiple genomic variants of SARS-CoV-2 were distributed in patients in New England early in the COVID-19 era, mostly in the B1 clade with a D614G mutation.