Comparative Effectiveness of BNT162b2 and mRNA-1273 Vaccines in U.S. Veterans

BACKGROUND

The messenger RNA (mRNA)-based vaccines BNT162b2 and mRNA-1273 are more than 90% effective against coronavirus disease 2019 (Covid-19). However, their comparative effectiveness for a range of outcomes across diverse populations is unknown.

METHODS

We emulated a target trial using the electronic health records of U.S. veterans who received a first dose of the BNT162b2 or mRNA-1273 vaccine between January 4 and May 14, 2021, during a period marked by predominance of the SARS-CoV-2 B.1.1.7 (alpha) variant. We matched recipients of each vaccine in a 1:1 ratio according to their risk factors. Outcomes included documented severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, symptomatic Covid-19, hospitalization for Covid-19, admission to an intensive care unit (ICU) for Covid-19, and death from Covid-19. We estimated risks using the Kaplan-Meier estimator. To assess the influence of the B.1.617.2 (delta) variant, we emulated a second target trial that involved veterans vaccinated between July 1 and September 20, 2021.

RESULTS

Each vaccine group included 219,842 persons. Over 24 weeks of follow-up in a period marked by alpha-variant predominance, the estimated risk of documented infection was 5.75 events per 1000 persons (95% confidence interval [CI], 5.39 to 6.23) in the BNT162b2 group and 4.52 events per 1000 persons (95% CI, 4.17 to 4.84) in the mRNA-1273 group. The excess number of events per 1000 persons for BNT162b2 as compared with mRNA-1273 was 1.23 (95% CI, 0.72 to 1.81) for documented infection, 0.44 (95% CI, 0.25 to 0.70) for symptomatic Covid-19, 0.55 (95% CI, 0.36 to 0.83) for hospitalization for Covid-19, 0.10 (95% CI, 0.00 to 0.26) for ICU admission for Covid-19, and 0.02 (95% CI, -0.06 to 0.12) for death from Covid-19. The corresponding excess risk (BNT162b2 vs. mRNA-1273) of documented infection over 12 weeks of follow-up in a period marked by delta-variant predominance was 6.54 events per 1000 persons (95% CI, -2.58 to 11.82).

CONCLUSIONS

The 24-week risk of Covid-19 outcomes was low after vaccination with mRNA-1273 or BNT162b2, although risks were lower with mRNA-1273 than with BNT162b2. This pattern was consistent across periods marked by alpha- and delta-variant predominance. (Funded by the Department of Veterans Affairs and others.).

Risk factors and prediction models for incident heart failure with reduced and preserved ejection fraction

AIMS

This study aims to develop the first race-specific and sex-specific risk prediction models for heart failure with preserved (HFpEF) and reduced ejection fraction (HFrEF).

METHODS AND RESULTS

We created a cohort of 1.8 million individuals who had an outpatient clinic visit between 2002 and 2007 within the Veterans Affairs (VA) Healthcare System and obtained information on HFpEF, HFrEF, and several risk factors from electronic health records (EHR). Variables were selected for the risk prediction models in a 'derivation cohort' that consisted of individuals with baseline date in 2002, 2003, or 2004 using a forward stepwise selection based on a change in C-index threshold. Discrimination and calibration were assessed in the remaining participants (internal 'validation cohort'). A total of 66 831 individuals developed HFpEF, and 92 233 developed HFrEF (52 679 and 71 463 in the derivation cohort) over a median of 11.1 years of follow-up. The HFpEF risk prediction model included age, diabetes, BMI, COPD, previous MI, antihypertensive treatment, SBP, smoking status, atrial fibrillation, and estimated glomerular filtration rate (eGFR), while the HFrEF model additionally included previous CAD. For the HFpEF model, C-indices were 0.74 (SE = 0.002) for white men, 0.76 (0.005) for black men, 0.79 (0.015) for white women, and 0.77 (0.026) for black women, compared with 0.72 (0.002), 0.72 (0.004), 0.77 (0.017), and 0.75 (0.028), respectively, for the HFrEF model. These risk prediction models were generally well calibrated in each race-specific and sex-specific stratum of the validation cohort.

CONCLUSIONS

Our race-specific and sex-specific risk prediction models, which used easily obtainable clinical variables, can be a useful tool to implement preventive strategies or subtype-specific prevention trials in the nine million users of the VA healthcare system and the general population after external validation.

Pharmacoepidemiology, Machine Learning, and COVID-19: An Intent-to-Treat Analysis of Hydroxychloroquine, With or Without Azithromycin, and COVID-19 Outcomes Among Hospitalized US Veterans

Hydroxychloroquine (HCQ) was proposed as an early therapy for coronavirus disease 2019 (COVID-19) after in vitro studies indicated possible benefit. Previous in vivo observational studies have presented conflicting results, though recent randomized clinical trials have reported no benefit from HCQ amongst hospitalized COVID-19 patients. We examined the effects of HCQ alone, and in combination with azithromycin, in a hospitalized COVID-19 positive, United States (US) Veteran population using a propensity score adjusted survival analysis with imputation of missing data. From March 1, 2020 through April 30, 2020, 64,055 US Veterans were tested for COVID-19 based on Veteran Affairs Healthcare Administration electronic health record data. Of the 7,193 positive cases, 2,809 were hospitalized, and 657 individuals were prescribed HCQ within the first 48-hours of hospitalization for the treatment of COVID-19. There was no apparent benefit associated with HCQ receipt, alone or in combination with azithromycin, and an increased risk of intubation when used in combination with azithromycin [Hazard Ratio (95% Confidence Interval): 1.55 (1.07, 2.24)]. In conclusion, we assessed the effectiveness of HCQ with or without azithromycin in treating patients hospitalized with COVID-19 using a national sample of the US Veteran population. Using rigorous study design and analytic methods to reduce confounding and bias, we found no evidence of a survival benefit from the administration of HCQ.

Trajectories of Frailty in the 5 Years Prior to Death Among U.S. Veterans Born 1927-1934

BACKGROUND

Electronic frailty indices (eFIs) are increasingly used to identify patients at risk for morbidity and mortality. Whether eFIs capture the spectrum of frailty change, including decline, stability, and improvement, is unknown.

METHODS

In a nationwide retrospective birth cohort of U.S. Veterans, a validated eFI, including 31 health deficits, was calculated annually using medical record and insurance claims data (2002-2012). K-means clustering was used to assign patients into frailty trajectories measured 5 years prior to death.

RESULTS

There were 214 250 veterans born between 1927 and 1934 (mean [SD] age at death = 79.4 [2.8] years, 99.2% male, 90.3% White) with an annual eFI in the 5 years before death. Nine frailty trajectories were identified. Those starting at nonfrail or prefrail had 2 stable trajectories (nonfrail to prefrail, n = 29 786 and stable prefrail, n = 28 499) and 2 rapidly increasing trajectories (prefrail to moderately frail, n = 28 244 and prefrail to severely frail, n = 22 596). Those who were mildly frail at baseline included 1 gradually increasing trajectory (mildly to moderately frail, n = 33 806) and 1 rapidly increasing trajectory (mildly to severely frail, n = 15 253). Trajectories that started at moderately or severely frail included 2 gradually increasing trajectories (moderately to severely frail, n = 27 662 and progressing severely frail, n = 14 478) and 1 recovering trajectory (moderately frail to mildly frail, n = 13 926).

CONCLUSIONS

Nine frailty trajectories, including 1 recovering trajectory, were identified in this cohort of older U.S. Veterans. Future work is needed to understand whether prevention and treatment strategies can improve frailty trajectories and contribute to compression of morbidity toward the end of life.

Early Convalescent Plasma Therapy and Mortality Among US Veterans Hospitalized With Nonsevere COVID-19: An Observational Analysis Emulating a Target Trial

BACKGROUND

Early convalescent plasma transfusion may reduce mortality in patients with nonsevere coronavirus disease 2019 (COVID-19).

METHODS

This study emulates a (hypothetical) target trial using observational data from a cohort of US veterans admitted to a Department of Veterans Affairs (VA) facility between 1 May and 17 November 2020 with nonsevere COVID-19. The intervention was convalescent plasma initiated within 2 days of eligibility. Thirty-day mortality was compared using cumulative incidence curves, risk differences, and hazard ratios estimated from pooled logistic models with inverse probability weighting to adjust for confounding.

RESULTS

Of 11 269 eligible person-trials contributed by 4755 patients, 402 trials were assigned to the convalescent plasma group. Forty and 671 deaths occurred within the plasma and nonplasma groups, respectively. The estimated 30-day mortality risk was 6.5% (95% confidence interval [CI], 4.0%-9.7%) in the plasma group and 6.2% (95% CI, 5.6%-7.0%) in the nonplasma group. The associated risk difference was 0.30% (95% CI, -2.30% to 3.60%) and the hazard ratio was 1.04 (95% CI, .64-1.62).

CONCLUSIONS

Our target trial emulation estimated no meaningful differences in 30-day mortality between nonsevere COVID-19 patients treated and untreated with convalescent plasma. Clinical Trials Registration. NCT04545047.

A randomised trial of a web-based physical activity self-management intervention in COPD

Improving exercise capacity is a primary objective in COPD. Declines in exercise capacity result in reduced physical activity and health-related quality of life (HRQoL). Self-management interventions can teach patients skills and behaviours to manage their disease. Technology-mediated interventions have the potential to provide easily accessible support for disease self-management. We evaluated the effectiveness of a web-based self-management intervention, focused on physical activity promotion, on exercise capacity in COPD. This 6-month randomised controlled trial (NCT02099799) enrolled 153 persons with COPD at two US sites (VABoston, n=108; VABirmingham, n=45). Participants were allocated (1:1) to the web-based self-management intervention (physical activity promotion through personalised, progressive step-count goals, feedback, online COPD-related education and social support via an online community) or usual care. The primary outcome was exercise capacity (6-min walk distance (6 MWD)). Secondary outcomes included physical activity (daily steps per day), HRQoL (St. George's Respiratory Questionnaire Total Score), dyspnoea, COPD-related knowledge and social support. Change in step-count goals reflected intervention engagement. Participants' mean age was 69 (sd=7), and mean forced expiratory volume in 1 s % predicted was 61% (sd=21%). Change in 6MWD did not differ between groups. Intervention participants improved their mean daily step counts by 1312 more than those in the usual care group (p<0.001). Groups did not differ on other secondary outcomes. VABirmingham participants were significantly more engaged with the intervention, although site did not modify the effect of the intervention on 6MWD or secondary outcomes. The intervention did not improve exercise capacity but improved physical activity at 6 months. Additional intervention modifications are needed to optimise its COPD self-management capabilities.

A web-based self-management intervention improved physical activity but not exercise capacity. There is a need to develop and study accessible self-management interventions for COPD.https://bit.ly/3iT1yvU

Phenome-wide association of 1809 phenotypes and COVID-19 disease progression in the Veterans Health Administration Million Veteran Program

BACKGROUND

The risk factors associated with the stages of Coronavirus Disease-2019 (COVID-19) disease progression are not well known. We aim to identify risk factors specific to each state of COVID-19 progression from SARS-CoV-2 infection through death.

METHODS AND RESULTS

We included 648,202 participants from the Veteran Affairs Million Veteran Program (2011-). We identified characteristics and 1,809 ICD code-based phenotypes from the electronic health record. We used logistic regression to examine the association of age, sex, body mass index (BMI), race, and prevalent phenotypes to the stages of COVID-19 disease progression: infection, hospitalization, intensive care unit (ICU) admission, and 30-day mortality (separate models for each). Models were adjusted for age, sex, race, ethnicity, number of visit months and ICD codes, state infection rate and controlled for multiple testing using false discovery rate (≤0.1). As of August 10, 2020, 5,929 individuals were SARS-CoV-2 positive and among those, 1,463 (25%) were hospitalized, 579 (10%) were in ICU, and 398 (7%) died. We observed a lower risk in women vs. men for ICU and mortality (Odds Ratio (95% CI): 0.48 (0.30-0.76) and 0.59 (0.31-1.15), respectively) and a higher risk in Black vs. Other race patients for hospitalization and ICU (OR (95%CI): 1.53 (1.32-1.77) and 1.63 (1.32-2.02), respectively). We observed an increased risk of all COVID-19 disease states with older age and BMI ≥35 vs. 20-24 kg/m2. Renal failure, respiratory failure, morbid obesity, acid-base balance disorder, white blood cell diseases, hydronephrosis and bacterial infections were associated with an increased risk of ICU admissions; sepsis, chronic skin ulcers, acid-base balance disorder and acidosis were associated with mortality.

CONCLUSIONS

Older age, higher BMI, males and patients with a history of respiratory, kidney, bacterial or metabolic comorbidities experienced greater COVID-19 severity. Future studies to investigate the underlying mechanisms associated with these phenotype clusters and COVID-19 are warranted.

Egg consumption, overall diet quality, and risk of type 2 diabetes and coronary heart disease: A pooling project of US prospective cohorts

BACKGROUND AND AIMS

Data on the relation of egg consumption with risk of type 2 diabetes (T2D) and coronary heart disease (CHD) are limited and inconsistent. Few studies have controlled for overall dietary patterns in egg-T2D or egg-CHD analyses, and it is unclear whether any observed elevated risks of T2D and CHD with frequent egg consumption is real or due to confounding by dietary habits. We tested the hypothesis that frequent egg consumption is associated with a higher risk of T2D and CHD risk after adjustment for overall dietary patterns among adults.

DESIGN

We used prospective cohort design to complete time-to-event analyses.

METHODS

We pooled de novo, harmonized, individual-level analyses from nine US cohorts (n = 103,811). Cox regression was used to estimate hazard ratios separately in each cohort adjusting for age, ethnicity, body mass index (BMI), exercise, smoking, alcohol intake, and dietary patterns. We pooled cohort-specific results using an inverse-variance weighted method to estimate summary relative risks.

RESULTS

Median age ranged from 25 to 72 years. Median egg consumption was 1 egg per week in most of the cohorts. While egg consumption up to one per week was not associated with T2D risk, consumption of ≥2 eggs per week was associated with elevated risk [27% elevated risk of T2D comparing 7+ eggs/week with none (95% CI: 16%-37%)]. There was little evidence for heterogeneity across cohorts and we observed similar conclusions when stratified by BMI. Overall, egg consumption was not associated with the risk of CHD. However, in a sensitivity analysis, there was a 30% higher risk of CHD (95% CI: 3%-56%) restricted to older adults consuming 5-6 eggs/week.

CONCLUSIONS

Our data showed an elevated risk of T2D with egg consumption of ≥2 eggs per week but not with <2 eggs/week. While there was no overall association of egg consumption with CHD risk, the elevated CHD observed with consumption of 5-6 eggs/week in older cohorts merits further investigation.

Actionable druggable genome-wide Mendelian randomization identifies repurposing opportunities for COVID-19

Drug repurposing provides a rapid approach to meet the urgent need for therapeutics to address COVID-19. To identify therapeutic targets relevant to COVID-19, we conducted Mendelian randomization analyses, deriving genetic instruments based on transcriptomic and proteomic data for 1,263 actionable proteins that are targeted by approved drugs or in clinical phase of drug development. Using summary statistics from the Host Genetics Initiative and the Million Veteran Program, we studied 7,554 patients hospitalized with COVID-19 and >1 million controls. We found significant Mendelian randomization results for three proteins (ACE2, P = 1.6 × 10-6; IFNAR2, P = 9.8 × 10-11 and IL-10RB, P = 2.3 × 10-14) using cis-expression quantitative trait loci genetic instruments that also had strong evidence for colocalization with COVID-19 hospitalization. To disentangle the shared expression quantitative trait loci signal for IL10RB and IFNAR2, we conducted phenome-wide association scans and pathway enrichment analysis, which suggested that IFNAR2 is more likely to play a role in COVID-19 hospitalization. Our findings prioritize trials of drugs targeting IFNAR2 and ACE2 for early management of COVID-19.

Chocolate consumption and risk of coronary artery disease: the Million Veteran Program

Background

Although previous studies have suggested cocoa products may promote cardiovascular health in the general population, no public data are available from patients receiving care in a national integrated health care system.

Objectives

We tested the hypothesis that regular chocolate consumption is associated with a lower risk of coronary artery disease (CAD) events among participants of the Million Veteran Program (MVP). Secondary analysis examined if the main hypothesis was observed among participants with type 2 diabetes.

Methods

We analyzed data from MVP participants who completed the food frequency section of the MVP Lifestyle Survey and were free of CAD at the time of survey completion. CAD events during follow-up (International Statistical Classification of Diseases Ninth Revision codes 410-411 and 413-414, and Tenth Revision codes I20-I25 except I25.2) were assessed using electronic health records. We fitted a Cox proportional hazard model to estimate the RR of CAD.

Results

Of 188,447 MVP enrollees with survey data, mean ± SD age was 64 ± 12.0 y and 90% were men. For regular chocolate (28.3 g/serving) consumption of <1 serving/mo, 1-3 servings/mo, 1 serving/wk, 2-4 servings/wk, and ≥5 servings/wk, crude incidence rates (per 1000 person-years) for fatal and nonfatal CAD events or coronary procedures were 20.2, 17.5, 16.7, 17.1, and 16.9, respectively, during a mean follow-up of 3.2 y. After adjusting for age, sex, race, and lifestyle factors, the corresponding HRs (95% CIs) were 1.00 (ref), 0.92 (0.87, 0.96), 0.88 (0.83, 0.93), 0.89 (0.84, 0.95), and 0.89 (0.84, 0.96), respectively (P for linear trend < 0.0001). In a secondary analysis of 47,265 diabetics, we did not observe a decreasing trend in CAD mortality among those who consumed ≥1 serving of chocolate a month compared with those who consumed <1 serving/mo.

Conclusions

Regular chocolate consumption was associated with a lower risk of CAD among veterans, but was not associated with cardiovascular disease risk in veterans with type 2 diabetes.

Revisiting methods for modeling longitudinal and survival data: Framingham Heart Study

Background

Statistical methods for modeling longitudinal and time-to-event data has received much attention in medical research and is becoming increasingly useful. In clinical studies, such as cancer and AIDS, longitudinal biomarkers are used to monitor disease progression and to predict survival. These longitudinal measures are often missing at failure times and may be prone to measurement errors. More importantly, time-dependent survival models that include the raw longitudinal measurements may lead to biased results. In previous studies these two types of data are frequently analyzed separately where a mixed effects model is used for the longitudinal data and a survival model is applied to the event outcome.

Methods

In this paper we compare joint maximum likelihood methods, a two-step approach and a time dependent covariate method that link longitudinal data to survival data with emphasis on using longitudinal measures to predict survival. We apply a Bayesian semi-parametric joint method and maximum likelihood joint method that maximizes the joint likelihood of the time-to-event and longitudinal measures. We also implement the Two-Step approach, which estimates random effects separately, and a classic Time Dependent Covariate Model. We use simulation studies to assess bias, accuracy, and coverage probabilities for the estimates of the link parameter that connects the longitudinal measures to survival times.

Results

Simulation results demonstrate that the Two-Step approach performed best at estimating the link parameter when variability in the longitudinal measure is low but is somewhat biased downwards when the variability is high. Bayesian semi-parametric and maximum likelihood joint methods yield higher link parameter estimates with low and high variability in the longitudinal measure. The Time Dependent Covariate method resulted in consistent underestimation of the link parameter. We illustrate these methods using data from the Framingham Heart Study in which lipid measurements and Myocardial Infarction data were collected over a period of 26 years.

Conclusions

Traditional methods for modeling longitudinal and survival data, such as the time dependent covariate method, that use the observed longitudinal data, tend to provide downwardly biased estimates. The two-step approach and joint models provide better estimates, although a comparison of these methods may depend on the underlying residual variance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12874-021-01207-y.

Development and validation of a 30-day mortality index based on pre-existing medical administrative data from 13,323 COVID-19 patients: The Veterans Health Administration COVID-19 (VACO) Index

BACKGROUND

Available COVID-19 mortality indices are limited to acute inpatient data. Using nationwide medical administrative data available prior to SARS-CoV-2 infection from the US Veterans Health Administration (VA), we developed the VA COVID-19 (VACO) 30-day mortality index and validated the index in two independent, prospective samples.

METHODS AND FINDINGS

We reviewed SARS-CoV-2 testing results within the VA between February 8 and August 18, 2020. The sample was split into a development cohort (test positive between March 2 and April 15, 2020), an early validation cohort (test positive between April 16 and May 18, 2020), and a late validation cohort (test positive between May 19 and July 19, 2020). Our logistic regression model in the development cohort considered demographics (age, sex, race/ethnicity), and pre-existing medical conditions and the Charlson Comorbidity Index (CCI) derived from ICD-10 diagnosis codes. Weights were fixed to create the VACO Index that was then validated by comparing area under receiver operating characteristic curves (AUC) in the early and late validation cohorts and among important validation cohort subgroups defined by sex, race/ethnicity, and geographic region. We also evaluated calibration curves and the range of predictions generated within age categories. 13,323 individuals tested positive for SARS-CoV-2 (median age: 63 years; 91% male; 42% non-Hispanic Black). We observed 480/3,681 (13%) deaths in development, 253/2,151 (12%) deaths in the early validation cohort, and 403/7,491 (5%) deaths in the late validation cohort. Age, multimorbidity described with CCI, and a history of myocardial infarction or peripheral vascular disease were independently associated with mortality-no other individual comorbid diagnosis provided additional information. The VACO Index discriminated mortality in development (AUC = 0.79, 95% CI: 0.77-0.81), and in early (AUC = 0.81 95% CI: 0.78-0.83) and late (AUC = 0.84, 95% CI: 0.78-0.86) validation. The VACO Index allows personalized estimates of 30-day mortality after COVID-19 infection. For example, among those aged 60-64 years, overall mortality was estimated at 9% (95% CI: 6-11%). The Index further discriminated risk in this age stratum from 4% (95% CI: 3-7%) to 21% (95% CI: 12-31%), depending on sex and comorbid disease.

CONCLUSION

Prior to infection, demographics and comorbid conditions can discriminate COVID-19 mortality risk overall and within age strata. The VACO Index reproducibly identified individuals at substantial risk of COVID-19 mortality who might consider continuing social distancing, despite relaxed state and local guidelines.

Egg consumption and risk of coronary artery disease in the Million Veteran Program

BACKGROUND & AIMS

Limited and inconsistent data are available on the relation between egg consumption and risk of myocardial infarction (MI) and it is unclear if adiposity or type 2 diabetes modifies egg-MI relation. We tested the primary hypothesis that egg consumption is positively associated with incidence of MI among veterans. In secondary analyses, we examined potential effect modification of egg-MI relation by adiposity and type 2 diabetes.

METHODS

We analyzed data collected on 188,267 US veterans who were enrolled in the Million Veteran Program (MVP) from 2011 to 2018. Information on egg consumption was obtained via self-administered food frequency questionnaire and we used electronic health records to identify incident MI.

RESULTS

The mean age was 64.4 (SD = 12.0) years and 9.9% of the population were female. We ascertained 10,260 new cases of non-fatal MI during an average follow up of 3.24 years (range: 0.002 to 7.49 y). Hazard ratio (95% CI) for non-fatal MI were 1.00 (ref), 0.93 (0.85-0.1.02), 0.96 (0.87-1.05), 0.98 (0.89-1.07), 1.08 (0.98-1.19), 1.11 (1.00-1.24), and 1.13 (1.00-1.28) for egg consumption of <1/month, 1-3/month, 1/week, 2-4/week, 5-6/week, 1/d, and 2+/d, respectively, controlling for age, sex, race, body mass index, smoking, exercise, alcohol intake, and overall dietary pattern (p non-linear trend 0.019). In secondary analyses, we observed similar results with a composite endpoint including fatal MI, coronary angioplasty and revascularization.

CONCLUSIONS

Our data showed no association of infrequent consumption of eggs with non-fatal MI but a slightly elevated risk with intake of 1 or more eggs per day among US veterans.

Optimizing Atherosclerotic Cardiovascular Disease Risk Estimation for Veterans With Diabetes Mellitus

BACKGROUND

Estimated 10-year atherosclerotic cardiovascular disease (ASCVD) risk in diabetes mellitus patients is used to guide primary prevention, but the performance of risk estimators (2013 Pooled Cohort Equations [PCE] and Risk Equations for Complications of Diabetes [RECODe]) varies across populations. Data from electronic health records could be used to improve risk estimation for a health system's patients. We aimed to evaluate risk equations for initial ASCVD events in US veterans with diabetes mellitus and improve model performance in this population.

METHODS AND RESULTS

We studied 183 096 adults with diabetes mellitus and without prior ASCVD who received care in the Veterans Affairs Healthcare System (VA) from 2002 to 2016 with mean follow-up of 4.6 years. We evaluated model discrimination, using Harrell's C statistic, and calibration, using the reclassification χ2 test, of the PCE and RECODe equations to predict fatal or nonfatal myocardial infarction or stroke and cardiovascular mortality. We then tested whether model performance was affected by deriving VA-specific β-coefficients. Discrimination of ASCVD events by the PCE was improved by deriving VA-specific β-coefficients (C statistic increased from 0.560 to 0.597) and improved further by including measures of glycemia, renal function, and diabetes mellitus treatment (C statistic, 0.632). Discrimination by the RECODe equations was improved by substituting VA-specific coefficients (C statistic increased from 0.604 to 0.621). Absolute risk estimation by PCE and RECODe equations also improved with VA-specific coefficients; the calibration P increased from <0.001 to 0.08 for PCE and from <0.001 to 0.005 for RECODe, where higher P indicates better calibration. Approximately two-thirds of veterans would meet a guideline indication for high-intensity statin therapy based on the PCE versus only 10% to 15% using VA-fitted models.

CONCLUSIONS

Existing ASCVD risk equations overestimate risk in veterans with diabetes mellitus, potentially impacting guideline-indicated statin therapy. Prediction model performance can be improved for a health system's patients using readily available electronic health record data.

Association of Statin Use With All-Cause and Cardiovascular Mortality in US Veterans 75 Years and Older

IMPORTANCE

Data are limited regarding statin therapy for primary prevention of atherosclerotic cardiovascular disease (ASCVD) in adults 75 years and older.

OBJECTIVE

To evaluate the role of statin use for mortality and primary prevention of ASCVD in veterans 75 years and older.

DESIGN, SETTING, AND PARTICIPANTS

Retrospective cohort study that used Veterans Health Administration (VHA) data on adults 75 years and older, free of ASCVD, and with a clinical visit in 2002-2012. Follow-up continued through December 31, 2016. All data were linked to Medicare and Medicaid claims and pharmaceutical data. A new-user design was used, excluding those with any prior statin use. Cox proportional hazards models were fit to evaluate the association of statin use with outcomes. Analyses were conducted using propensity score overlap weighting to balance baseline characteristics.

EXPOSURES

Any new statin prescription.

MAIN OUTCOMES AND MEASURES

The primary outcomes were all-cause and cardiovascular mortality. Secondary outcomes included a composite of ASCVD events (myocardial infarction, ischemic stroke, and revascularization with coronary artery bypass graft surgery or percutaneous coronary intervention).

RESULTS

Of 326 981 eligible veterans (mean [SD] age, 81.1 [4.1] years; 97% men; 91% white), 57 178 (17.5%) newly initiated statins during the study period. During a mean follow-up of 6.8 (SD, 3.9) years, a total 206 902 deaths occurred including 53 296 cardiovascular deaths, with 78.7 and 98.2 total deaths/1000 person-years among statin users and nonusers, respectively (weighted incidence rate difference [IRD]/1000 person-years, -19.5 [95% CI, -20.4 to -18.5]). There were 22.6 and 25.7 cardiovascular deaths per 1000 person-years among statin users and nonusers, respectively (weighted IRD/1000 person-years, -3.1 [95 CI, -3.6 to -2.6]). For the composite ASCVD outcome there were 123 379 events, with 66.3 and 70.4 events/1000 person-years among statin users and nonusers, respectively (weighted IRD/1000 person-years, -4.1 [95% CI, -5.1 to -3.0]). After propensity score overlap weighting was applied, the hazard ratio was 0.75 (95% CI, 0.74-0.76) for all-cause mortality, 0.80 (95% CI, 0.78-0.81) for cardiovascular mortality, and 0.92 (95% CI, 0.91-0.94) for a composite of ASCVD events when comparing statin users with nonusers.

CONCLUSIONS AND RELEVANCE

Among US veterans 75 years and older and free of ASCVD at baseline, new statin use was significantly associated with a lower risk of all-cause and cardiovascular mortality. Further research, including from randomized clinical trials, is needed to more definitively determine the role of statin therapy in older adults for primary prevention of ASCVD.

Estimation of Atherosclerotic Cardiovascular Disease Risk Among Patients in the Veterans Affairs Health Care System

IMPORTANCE

Current guidelines recommend statin therapy for millions of US residents for the primary prevention of atherosclerotic cardiovascular disease (ASCVD). It is unclear whether traditional prediction models that do not account for current widespread statin use are sufficient for risk assessment.

OBJECTIVES

To examine the performance of the Pooled Cohort Equations (PCE) for 5-year ASCVD risk estimation in a contemporary cohort and to test the hypothesis that inclusion of statin therapy improves model performance.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study included adult patients in the Veterans Affairs health care system without baseline ASCVD. Using national electronic health record data, 3 Cox proportional hazards models were developed to estimate 5-year ASCVD risk, as follows: the variables and published β coefficients from the PCE (model 1), the PCE variables with cohort-derived β coefficients (model 2), and model 2 plus baseline statin use (model 3). Data were collected from January 2002 to December 2012 and analyzed from June 2016 to March 2020.

EXPOSURES

Traditional ASCVD risk factors from the PCE plus baseline statin use.

MAIN OUTCOMES AND MEASURES

Incident ASCVD and ASCVD mortality.

RESULTS

Of 1 672 336 patients in the cohort (mean [SD] baseline age 58.0 [13.8] years, 1 575 163 [94.2%] men, 1 383 993 [82.8%] white), 312 155 (18.7%) were receiving statin therapy at baseline. During 5 years of follow-up, 66 605 (4.0%) experienced an ASCVD event, and 31 878 (1.9%) experienced ASCVD death. Compared with the original PCE, the cohort-derived model did not improve model discrimination in any of the 4 age-sex strata but did improve model calibration. The PCE overestimated ASCVD risk compared with the cohort-derived model; 211 237 of 1 136 161 white men (18.6%), 29 634 of 218 463 black men (13.6%), 1741 of 44 399 white women (3.9%), and 836 of 16 034 black women (5.2%) would be potentially eligible for statin therapy under the PCE but not the cohort-derived model. When added to the cohort-derived model, baseline statin therapy was associated with a 7% (95% CI, 5%-9%) lower relative risk of ASCVD and a 25% (95% CI, 23%-28%) lower relative risk for ASCVD death.

CONCLUSIONS AND RELEVANCE

In this study, lower than expected rates of incident ASCVD events in a contemporary national cohort were observed. The PCE overestimated ASCVD risk, and more than 15% of patients would be potentially eligible for statin therapy based on the PCE but not on a cohort-derived model. In the statin era, health care professionals and systems should base ASCVD risk assessment on models calibrated to their patient populations.

Standardized Architecture for a Mega-Biobank Phenomic Library: The Million Veteran Program (MVP)

Electronic health records (EHRs) provide a wealth of data for phenotype development in population health studies, and researchers invest considerable time to curate data elements and validate disease definitions. The ability to reproduce well-defined phenotypes increases data quality, comparability of results and expedites research. In this paper, we present a standardized approach to organize and capture phenotype definitions, resulting in the creation of an open, online repository of phenotypes. This resource captures phenotype development, provenance and process from the Million Veteran Program, a national mega-biobank embedded in the Veterans Health Administration (VHA). To ensure that the repository is searchable, extendable, and sustainable, it is necessary to develop both a proper digital catalog architecture and underlying metadata infrastructure to enable effective management of the data fields required to define each phenotype. Our methods provide a resource for VHA investigators and a roadmap for researchers interested in standardizing their phenotype definitions to increase portability.

Diabetes Mellitus-Related All-Cause and Cardiovascular Mortality in a National Cohort of Adults

Background Diabetes mellitus is a risk factor for cardiovascular disease ( CVD ) and has been associated with 2- to 4-fold higher mortality. Diabetes mellitus-related mortality has not been reassessed in individuals receiving routine care in the United States in the contemporary era of CVD risk reduction. Methods and Results We retrospectively studied 963 648 adults receiving care in the US Veterans Affairs Healthcare System from 2002 to 2014; mean follow-up was 8 years. We estimated associations of diabetes mellitus status and hemoglobin A1c (HbA1c) with all-cause and CVD mortality using covariate-adjusted incidence rates and multivariable Cox proportional hazards regression. Of participants, 34% had diabetes mellitus. Compared with nondiabetic individuals, patients with diabetes mellitus had 7.0 (95% CI , 6.7-7.4) and 3.5 (95% CI, 3.3-3.7) deaths/1000-person-years higher all-cause and CVD mortality, respectively. The age-, sex-, race-, and ethnicity-adjusted hazard ratio for diabetes mellitus-related mortality was 1.29 (95% CI, 1.28-1.31), and declined with adjustment for CVD risk factors (hazard ratio, 1.18 [95% CI, 1.16-1.19]) and glycemia (hazard ratio, 1.03 [95% CI, 1.02-1.05]). Among individuals with diabetes mellitus, CVD mortality increased as HbA1c exceeded 7% (hazard ratios, 1.11 [95% CI, 1.08-1.14], 1.25 [95% CI, 1.22-1.29], and 1.52 [95% CI, 1.48-1.56] for HbA1c 7%-7.9%, 8%-8.9%, and ≥9%, respectively, relative to HbA1c 6%-6.9%). HbA1c 6% to 6.9% was associated with the lowest mortality risk irrespective of CVD history or age. Conclusions Diabetes mellitus remains significantly associated with all-cause and CVD mortality, although diabetes mellitus-related excess mortality is lower in the contemporary era than previously. We observed a gradient of mortality risk with increasing HbA1c >6% to 6.9%, suggesting HbA1c remains an informative predictor of outcomes even if causality cannot be inferred.

Association of pulse rate with outcomes in heart failure with reduced ejection fraction: a retrospective cohort study

BACKGROUND

In a real-world setting, the effect of pulse rate measured at the time of diagnosis and serially during follow-up and management, on outcomes in heart failure with reduced ejection fraction (HFrEF), has not been well-studied. Furthermore, how beta-blockade use in a real-world situation modifies this relation between pulse rate and outcomes in HFrEF is not well-known. Hence, we identified a large, national, real-world cohort of HFrEF to examine the association of pulse rate and outcomes.

METHODS

Using Veterans Affairs (VA) national electronic health records we identified incident HFrEF cases between 2006 and 2012. We examined the associations of both baseline and serially measured pulse rates, with mortality and days hospitalized per year for heart failure and for any cause, using crude and multivariable Cox proportional hazards and Poisson or negative binomial models, respectively. The exposure was examined as continuous, dichotomous, and categorical. Post-hoc analyses addressed the interaction of pulse rate and beta-blocker target dose.

RESULTS

We identified 51,194 incident HFrEF cases (67 ± 12 years, 98% male, 77% white. A significant positive, near linear relationship was observed for both baseline and serially measured pulse rates with all-cause mortality, all-cause hospitalization and heart failure hospitalization after adjusting for covariates including beta-blocker use. Patients who had a pulse rate ≥ 70 bpm in the past 6 months had 36% (95% CI: 31-42%), 25% (95% CI: 19-32%), and 51% (95% CI: 33-72%) increased rates of mortality, all-cause hospitalization, and heart failure hospitalization, respectively, compared to patients with pulse rates < 70 bpm. A minority of subjects (15%) were treated with guideline directed beta blockade ≥50% of recommended target dose, among whom better outcomes were seen compared to those who did not achieve target dose in patients with pulse rates both above and below 70 beats per minute.

CONCLUSIONS

High pulse rate, both at the time of diagnosis and during follow-up, is strongly associated with increased risk of adverse outcomes in HFrEF patients, independent of the use of beta-blockers. In a real-world setting, the majority of HFrEF patients do not achieve target dose of beta-blockade; greater use of strategies to reduce heart rate may improve outcomes in HFrEF.

Plasma vitamin D, past chest illness, and risk of future chest illness in chronic spinal cord injury (SCI): a longitudinal observational study

STUDY DESIGN

Observational study.

OBJECTIVE

Assess associations between vitamin D levels and other risk factors on future chest illness in a chronic spinal cord injury (SCI) cohort.

SETTING

Veterans Affairs Boston and the Boston, MA community.

METHODS

Between August 2009 and August 2017, 253 participants with chronic SCI were followed over a median of 3.2 years (up to 7.4 years) with two to four visits a median of 1.7 years apart. At each visit, plasma 25-hydroxyvitamin D level was obtained, spirometry performed, and a respiratory questionnaire assessing chest illnesses since last visit was completed. Repeated measures negative binomial regression was used to assess chest illness risk longitudinally.

RESULTS

At entry, 25% had deficient vitamin D levels (<20 nanograms/milliliter (ng/ml)), 52% were insufficient (20 to <30 ng/ml), and 23% were sufficient (≥30 ng/ml). Over 545 study visits, chest illnesses (n = 106) were reported by 60 participants. In multivariable models (including previous chest illness history), deficient vitamin D levels (compared with those with sufficient levels) were associated with future chest illness though with wide confidence limits (relative risk (RR) = 1.36, 95% confidence intervals (CI) = 0.74, 2.47). The strongest association with chest illness during the follow-up period was in persons who reported pneumonia/bronchitis after injury and a chest illness in the three years before study entry (RR = 7.62; 95% CI = 3.70, 15.71).

CONCLUSION

Assessed prospectively in chronic SCI, there was a suggestive association between deficient vitamin D levels and future chest illness. Past chest illness history was also strongly associated with future chest illness.

Long-term effects of web-based pedometer-mediated intervention on COPD exacerbations

BACKGROUND

Technology-based physical activity (PA) interventions have been shown to improve daily step counts and health-related quality of life, but their effect on long-term clinical outcomes like acute exacerbations (AEs) is unknown in persons with COPD.

METHODS

U.S. Veterans with stable COPD were randomized (1:1) to either pedometer alone (control) or pedometer plus a website with feedback, goal-setting, disease education, and a community forum (intervention) for 3 months. AEs were assessed every 3 months over a follow-up period of approximately 15 months. Pedometer-assessed daily step counts, health-related quality-of-life (HRQL), and self-efficacy were assessed at baseline, end-of-intervention at 3 months, and during follow-up approximately 6 and 12 months after enrollment. Zero-inflated Poisson models assessed the effect of the intervention on risk for AEs, compared to controls. Generalized linear mixed-effects models for repeated measures examined between-group and within-group changes in daily step count, HRQL, and self-efficacy.

RESULTS

There were no significant differences in age, FEV₁% predicted, baseline daily step count, AEs the year prior to enrollment, or duration of follow-up between the intervention (n = 57) and control (n = 52) groups. The intervention group had a significantly reduced risk of AEs (rate ratio = 0.51, [95%CI 0.31-0.85]), compared to the control group. There were no significant between-group differences in change in average daily step count, HRQL, or self-efficacy at 6 and 12 months after enrollment.

CONCLUSIONS

A 3-month internet-mediated, pedometer-based PA intervention was associated with reduced risk for AEs of COPD over 12-15 months of follow-up. ClinicalTrials.gov identifier: NCT01772082.

Association Between Early Hypertension Control and Cardiovascular Disease Incidence in Veterans With Diabetes

OBJECTIVE

Guidelines for hypertension treatment in patients with diabetes diverge regarding the systolic blood pressure (SBP) threshold at which treatment should be initiated and treatment goal. We examined associations of early SBP treatment with atherosclerotic cardiovascular disease (ASCVD) events in U.S. adults with diabetes.

RESEARCH DESIGN AND METHODS

We studied 43,986 patients with diabetes who newly initiated antihypertensive therapy between 2002 and 2007. Patients were classified into categories based on SBP at treatment initiation (130-139 or ≥140 mmHg) and after 2 years of treatment (100-119, 120-129, 130-139, 140-159, and ≥160 mmHg). The primary outcome was composite ASCVD events (fatal and nonfatal myocardial infarction and stroke), estimated using inverse probability of treatment-weighted Poisson regression and multivariable Cox proportional hazards regression.

RESULTS

Relative to individuals who initiated treatment when SBP was 130-139 mmHg, those with pretreatment SBP ≥140 mmHg had higher ASCVD risk (hazard ratio 1.10 [95% CI 1.02, 1.19]). Relative to those with pretreatment SBP of 130-139 mmHg and on-treatment SBP of 120-129 mmHg (reference group), ASCVD incidence was higher in those with pretreatment SBP ≥140 mmHg and on-treatment SBP 120-129 mmHg (adjusted incidence rate difference [IRD] 1.0 [-0.2 to 2.1] events/1,000 person-years) and in those who achieved on-treatment SBP 130-139 mmHg (IRD 1.9 [0.6, 3.2] and 1.1 [0.04, 2.2] events/1,000 person-years for those with pretreatment SBP 130-139 mmHg and ≥140 mmHg, respectively).

CONCLUSIONS

In this observational study, patients with diabetes initiating antihypertensive therapy when SBP was 130-139 mmHg and those achieving on-treatment SBP <130 mmHg had better outcomes than those with higher SBP levels when initiating or after 2 years on treatment.