Disease severity and outcomes of COVID-19 caused by SARS-CoV-2 vary among individuals who become infected, with several factors that have been suggested as predictors of mortality, including Charlson comorbidity index score, age, and BMI [1-6]. Other comorbidities such as cardiovascular disease, diabetes, and dementia are prevalent in patients who are hospitalized and can predict worse outcomes and complications following infection, but the current literature shows variable impacts of comorbidities in different populations [7-15].

Virologic characteristics have been suggested to impact the severity of disease, and concern has been raised about several variants being more transmissible or more lethal [16-18]. A dynamic nomenclature system (Pangolin) was developed to classify SARS-CoV-2 and identify lineages and mutations that could impact infectivity and virulence [19]. In the early era of COVID-19, there was significant concern about the D614G variant being more infectious, and SARS-CoV-2 with this mutation had been found to cause infections in New York State [17,20]. Evaluation of prevalence of this virus in the nearby states therefore became of substantial interest. Viral epidemiology and regional evaluation of viral variants, with their clinical correlations, are important to provide a full understanding of the disease.

Given the high variability and conflicting data in predicting who will have poor outcomes, assessment of specific populations is necessary to give providers the best clinical picture on their patients. The US veteran population is predominantly male (89%) with an average age of 58 years, compared to the general US population that has a median age of 38.5 years and a gender distribution of 98.2 males to every 100 females [21,22]. Both older age and male sex have been found to be independent predictors of greater disease severity and mortality in COVID-19 infections [8,13,23-26]. Evaluation of clinical outcomes among US veterans, a population predisposed to poorer outcomes because of these demographic variables, therefore, is important. Clinical outcomes of veterans with COVID-19 in New England and SARS-CoV-2 genomics have not been described. Outcome assessment based on unique demographic and clinical variables in people infected with novel viruses cannot be overemphasized. Given the high mutability of SARS-CoV-2, changing epidemiological trends over time, and known impact of virological factors on clinical outcomes, this study becomes even more important as it will provide an insight into this unique regional population with COVID-19 at a relevant time period in the evolution of the virus.

We aimed to describe patient characteristics, comorbidities, and disease factors that impact patient outcomes and present data on the genomic composition of the SARS-CoV-2 infecting these patients.

This study aims to determine the clinical and virologic factors impacting outcomes in veterans with COVID-19.

We conducted a retrospective chart review to gather demographic and clinical variables as well as clinical outcomes. The Veterans Affairs (VA) health care system has electronic medical records that can be accessed to extract this information. We have conducted such chart reviews before to help inform management decisions based on predictors of outcomes [27].

In 2020, the VA health care system in West Haven, Connecticut had been entrusted with testing for SARS-CoV-2 for all six VA health care centers in New England states (Connecticut, Massachusetts, Maine, New Hampshire, Rhode Island, Vermont). The virus, isolated from samples testing positive, was then sent for whole genome sequencing (WGS) under an agreement funded by the Centers for Disease Control and Prevention (CDC).

This study included all veterans who tested positive for COVID-19 from April 8, 2020, to September 16, 2020, at any of the six New England VA hospitals. Inclusion criteria included patients with accessible chart records and a diagnosis of COVID-19 based on one of three polymerase chain reaction–based tests: Xpert Xpress SARS-CoV-2 (Cepheid), Simplexa COVID-19 Direct Kit (DiaSorin), and Roche cobas 6800 system.

We manually reviewed charts and recorded demographics (age, gender, race, BMI, long-term care [LTC] facility status, and state of residence when diagnosed with COVID-19). Comorbidities recorded were immunosuppression, dementia, diabetes mellitus, chronic kidney disease stage 3, chronic liver disease, coronary artery disease (CAD), heart failure, atrial fibrillation, chronic obstructive pulmonary disease, asthma, and active tobacco use. All data collection was retrospective after a diagnosis of COVID-19 had been confirmed. If chart review occurred while a veteran was hospitalized, the chart was again reviewed retrospectively after discharge from the hospital.

Handling of nasopharyngeal specimens or isolated virus was carried out by the VA Connecticut Healthcare System (VACHS) clinical laboratory as part of clinical care, following standardized Clinical Laboratory Improvement Amendments guidelines [28]. Our viral repository was populated by the positive test results of all New England veterans. The VACHS laboratory handled specimens, isolated the SARS-CoV-2 RNA, and shipped it for WGS to a non-VA laboratory. We obtained the details of the platform used to diagnose, the cycle threshold (Ct), and the date of test from the laboratory. Sequencing of viral genomes was conducted at the non-VA laboratory by our coauthors as follows.

Whole virus genomes were sequenced (≥20x coverage depth across ≥70% of the genome) using the Illumina (n=238) and Nanopore (n=61) platforms. WGS was conducted on SARS-CoV-2 isolates with a Ct value <36 and provided near-complete or complete genome results where the Ct value was <30. Using BWA-MEM version 0.7.15 (GNU Project), we aligned reads to the Wuhan-Hu-1 reference genomes (GenBank MN908937.3). With iVar v1.2.1 (Andersen Lab), 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 <10 reads were present at a site. We aligned consensus genomes with MAFFT (GNU Project) [29] and masked problematic sites [30]. We built a phylogenetic tree with IQTree (IQ-Tree) [31] using an HKY substitution model and 1000 bootstraps, visualized it with Python module baltic v0.1.5 (Python Software Foundation), and assigned lineages with Pangolin [19].

Our categorical outcomes, also derived from manual chart review, were hospitalization status, mortality, and oxygen (O₂) requirement within 24 hours of admission from manual chart review. We divided patients based on their peak disease severity (ordinal) into five categories depending on oxygenation requirements: (1) no O₂ requirement, (2) 1 to 3 liters by nasal cannula (NC), (3) 4 to 6 liters by NC, (4) >6 liters O₂ or noninvasive positive pressure ventilation, and (5) mechanical ventilation. Mortality was defined as death within 60 days from the date of diagnosis.

We used STATA v16 (StataCorp) for logistic regressions to predict our hospitalization and mortality, and ordinal logistic regression to predict peak disease severity. We first conducted a univariate analysis, then used significant variables from the univariate analysis (P<.05) to use in a multivariate model for each of our outcomes to assess the impact of several variables at once, which has been frequently used in COVID-19 literature [9,10,32-34]. Assumptions for logistic regressions (binary outcome, linearity, no outliers, and multicollinearity) were tested and met, with maximum variance inflation factors of 2. Genomic characteristics were reported descriptively.

The VACHS Institutional Review Board approved the creation and maintenance of a data repository of all veterans in New England diagnosed with COVID-19 and a viral repository of the SARS-CoV-2 RNA received from all six New England facilities. This study was conducted in accordance with the Declaration of Helsinki, keeping all private health information secure in approved secure folders behind a VA firewall. The RECORD (Reporting of Studies Conducted Using Observational Routinely-Collected Data) statement guidelines were used to maintain transparency in the reporting of this study [35].

Of 476 veterans in six New England states with confirmed SARS-CoV-2 during the study period, 274 had complete and accessible charts. Of 274 veterans, 93% (n=254) were men, 83% (n=228) were White, the mean age was 63 (SD 17.6) years, and over one-third resided were in LTC (n=92; Table 1). The most common comorbidities were CAD (n=74), diabetes (n=68), and tobacco use (n=62).

Notably, 12% (n=33) of patients required O₂ above their baseline home O₂ requirement within 24 hours of admission, and 21% (n=58) of all patients required O₂ support at some point during hospitalization. In terms of peak severity, 79% (n=216) required only room air, 11% (n=30) required 1 to 3 liters O₂, 4.0% (n=11) required 4 to 6 liters O₂, 3.6% (n=10) required >6 liters O₂ or noninvasive positive-pressure ventilation, and 2.2% (n=6) were intubated. The hospitalization rate was 29% (n=79; Figure 1), and the overall mortality rate was 11% (n=30; Figure 2).

Univariate regression analysis results are reported in Table 2. In the multivariate regression (Table 2), significant predictors of hospitalization (C statistic 0.75) were age (odds ratio [OR] 1.05, 95% CI 1.03-1.08) and non-White race (OR 2.39, 95% CI 1.13-5.01; Table 2). Peak severity (C statistic 0.70) also varied by age (OR 1.07, 95% CI 1.03-1.11) and O₂ requirement on admission (OR 45.7, 95% CI 18.79-111). Mortality (C statistic 0.87) was predicted by age (OR 1.06, 95% CI 1.01-1.11), dementia (OR 3.44, 95% CI 1.07-11.1), and O₂ requirement on admission (OR 6.74, 95% CI 1.74-26.1). In other words, 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 genomic characteristics, among the 476 patients, 299 patients’ genomes had adequate coverage for analysis. We found the majority of our specimens (154/299, 51.5%) were from SARS-CoV-2 lineage B.1 or a sublineage of B.1 (eg, B.1.302, B.1.303, B.1.356; 137/299, 46%), all of which are defined by D614G substitution (Figure 3, Multimedia Appendix 1). Only 2.4% (7/299) were from the lineage A that lack the D614G mutation. There were 41 different SARS-CoV-2 lineages detected in our cohort, and we did not have the power to test for clinical correlates. Our sequencing data does inform us that the outcomes presented in this VA cohort are dominated by the impacts of the B lineage D614G mutation.

Our study found that in a cohort of veterans with an average age of 63 years and a high comorbidity burden, age was significantly associated with risk of hospitalization, peak disease severity, and mortality. O₂ requirement upon admission correlated with peak disease severity and mortality, while dementia was an additional factor associated with higher mortality. The CDC provides a list of chronic medical conditions (May 2021) that predispose individuals to severe illness from SARS-CoV-2 infection [6]. Based on this list, >75% of US adults fall under a high-risk category [36], therefore making it important to have select populations evaluated for uniquely applicable risk factors. Veterans are a unique cohort because of advanced age on average [22] and more comorbidities. Understanding clinical factors that impact outcomes in veterans will help health care providers risk-stratify patients with similar demographic profiles, and future research should explore the impact of new treatments and vaccination on outcomes. The predominance of B lineage D614G in our study specimens provided valuable insight into the pace of epidemiological trend and the evolution of the virus early in the COVID-19 era through the New England region.

Many COVID-19 studies have found age to be a predictor of worse outcomes [3,37-40]. In our study, age was a significant predictor for all of the studied outcomes and was a confounder for other variables. Accordingly, LTC status predicted all three of our outcomes on univariate analysis but not on multivariate analyses, possibly because LTC units tend to have older residents. Earlier in the COVID-19 pandemic, residents of nursing homes had higher rates of infection and severe illness and mortality [41]. Our study shows disease outcomes were not impacted by their residence status after adjusting for age.

Concurrent work from our group suggests that O₂ requirement within 24 hours of admission predicts poor outcomes in veterans, which has helped inform the triage guidelines at our health care system. This is an important finding because other ways of determining oxygenation status can frequently change and thus become difficult for clinicians to use in practice [22].

Our study supports data from previous reports that non-White patients in the United States are at increased risk of hospitalization but have similar peak severity and mortality outcomes [42,43]. Many studies have shown that minorities often have delays in seeking care, causing higher risk of hospitalization when they do seek care [44,45]. This may explain the outcomes in our study. 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.

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 [37,46,47]. This may be explained by a host of biological factors but also may be a result of the inability to self-report symptoms. This finding emphasizes the importance of extra care and monitoring required when approaching a patient with dementia.

Limitations of this study include the smaller sample size. Furthermore, our study is specific to veterans, which is a largely male and older cohort, and results may not therefore be generalizable. The time period of this study was prior to established medical therapies for COVID-19, and our reported outcomes are likely worse than expected today. Strengths of our study include its comprehensive scope, wide geographic range, manual chart review process 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.

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 from a B.1 sublineage with the spike D614G mutation.