Search

Introduction: Cardiovascular disease (CVD) is the leading cause of mortality and disability in the United States. Spatiotemporal modeling of disease incidence may guide prevention efforts toward regions most at risk in the future. Methods: We examined spatiotemporal trends in composite CVD events, defined as myocardial infarction, ischemic stroke, atrial fibrillation, heart failure, coronary artery disease, or cardiovascular death, using electronic health records from 9,718,107 U.S. Veterans who used VA healthcare facilities from 2003-2018. Age-standardized annual incidence by state was modeled using the Bernardinelli model, a Bayesian Poisson regression with linear temporal trends for the US and individual states. Results: There were 626,271 CVD events over 16 years of follow up and annual incidence fell gradually from 11.5% in 2003 to 6.6% in 2018. Incidence among all US veterans decreased an estimated 3.1% [95% CI: -3.3%, -2.9%] per year (RR year = 0.969). State-level incidences also decreased monotonically over the follow-up period (Figure panel A ), ranging from an absolute decrease 5.0% (RI) to 1.7% (OR). Estimated incidence for the five states most above (RI, FL, NC, MT, LA) and most below (OR, NV, CT, OK, MI) the national trend are plotted against the US average in Figure panel B . These outlying states differed significantly from the overall trend (posterior probabilities ≥ 99.7%), and the magnitude of trend was loosely correlated with initial incidence rate (in 2003). Notably, four states crossed the national average over the study period: CT and OK started below average in 2003 but ended above in 2018, whereas RI and NC started with high incidence and ended with below average incidence in 2018. Conclusions: Among US veterans, there was significant geographic variation in the rate of decline of incident CVD from 2003-2018. Future research will investigate the mechanisms underpinning the observed trends, particularly prevalence and control of risk factors.

Background: End-stage renal disease is associated with a high risk of cardiovascular events. It is unknown, however, whether mild-to-moderate kidney dysfunction is causally related to coronary heart disease (CHD) and stroke. Methods: Observational analyses were conducted using individual-level data from 4 population data sources (Emerging Risk Factors Collaboration, EPIC-CVD [European Prospective Investigation into Cancer and Nutrition–Cardiovascular Disease Study], Million Veteran Program, and UK Biobank), comprising 648 135 participants with no history of cardiovascular disease or diabetes at baseline, yielding 42 858 and 15 693 incident CHD and stroke events, respectively, during 6.8 million person-years of follow-up. Using a genetic risk score of 218 variants for estimated glomerular filtration rate (eGFR), we conducted Mendelian randomization analyses involving 413 718 participants (25 917 CHD and 8622 strokes) in EPIC-CVD, Million Veteran Program, and UK Biobank. Results: There were U-shaped observational associations of creatinine-based eGFR with CHD and stroke, with higher risk in participants with eGFR values 105 mL·min –1 ·1.73 m –2 , compared with those with eGFR between 60 and 105 mL·min –1 ·1.73 m –2 . Mendelian randomization analyses for CHD showed an association among participants with eGFR –1 ·1.73 m –2 , with a 14% (95% CI, 3%–27%) higher CHD risk per 5 mL·min –1 ·1.73 m –2 lower genetically predicted eGFR, but not for those with eGFR >105 mL·min –1 ·1.73 m –2 . Results were not materially different after adjustment for factors associated with the eGFR genetic risk score, such as lipoprotein(a), triglycerides, hemoglobin A1c, and blood pressure. Mendelian randomization results for stroke were nonsignificant but broadly similar to those for CHD. Conclusions: In people without manifest cardiovascular disease or diabetes, mild-to-moderate kidney dysfunction is causally related to risk of CHD, highlighting the potential value of preventive approaches that preserve and modulate kidney function.

We conducted genome-wide association analyses of over 250,000 participants of European (EUR) and African (AFR) ancestry from the Million Veteran Program using electronic health record-validated post-traumatic stress disorder (PTSD) diagnosis and quantitative symptom phenotypes. Applying genome-wide multiple testing correction, we identified three significant loci in European case-control analyses and 15 loci in quantitative symptom analyses. Genomic structural equation modeling indicated tight coherence of a PTSD symptom factor that shares genetic variance with a distinct internalizing (mood–anxiety–neuroticism) factor. Partitioned heritability indicated enrichment in several cortical and subcortical regions, and imputed genetically regulated gene expression in these regions was used to identify potential drug repositioning candidates. These results validate the biological coherence of the PTSD syndrome, inform its relationship to comorbid anxiety and depressive disorders and provide new considerations for treatment. Genome-wide association analyses of post-traumatic stress disorder and its symptom subdomains in individuals of European and African ancestry provide insights into its relationship with anxiety and depressive disorders and identify potential candidates for drug repositioning.

The prediction of myocardial infarction (MI), acute ischemic stroke (AIS), cardiovascular (CVD) death, and 3-point MACE (MI, AIS, or CVD death) in Million Veteran Program (MVP) participants over 8 years of follow up was evaluated using environmental/traditional risk factors (“E”), polygenic risk scores (“G”), and a combination of both approaches (“GxE”). Individuals free of atherosclerotic cardiovascular disease (ASCVD) at baseline were included. Risk factor levels for the participants close to the time of MVP enrollment for age, sex, systolic blood pressure, cholesterol, HDL cholesterol, smoking status, and diabetes status were analyzed. Outcomes were determined from VA electronic record data, Medicare/Medicaid data, and the National Death Index. Analyses were undertaken separately for non-Hispanic European (EUR), African American (AFR), and Hispanic (HIS) participants. There were 157,941 veterans at risk with 8,157 MI events, 2,024 AIS events, 1,778 CVD deaths, and 9,350 3-pt MACE events over 8 years of follow-up. The overall results showed good performance with the environmental model for all three racial-ethnic groups, with C-statistics ranging from 0.69 to 0.77. The G models showed very modest prediction capabilities and similarly modest improvement in the combined GxE models. In conclusion, traditional risk factor modeling has been shown to be highly effective and, in the MVP experience, the additional impact of genetic or genetic interaction information was small.

SARS-CoV-2 has caused symptomatic COVID-19 and widespread death across the globe. We sought to determine genetic variants contributing to COVID-19 susceptibility and hospitalization in a large biobank linked to a national United States health system. We identified 19,168 (3.7%) lab-confirmed COVID-19 cases among Million Veteran Program participants between March 1, 2020, and February 2, 2021, including 11,778 Whites, 4,893 Blacks, and 2,497 Hispanics. A multi-population genome-wide association study (GWAS) for COVID-19 outcomes identified four independent genetic variants (rs8176719, rs73062389, rs60870724, and rs73910904) contributing to COVID-19 positivity, including one novel locus found exclusively among Hispanics. We replicated eight of nine previously reported genetic associations at an alpha of 0.05 in at least one population-specific or the multi-population meta-analysis for one of the four MVP COVID-19 outcomes. We used rs8176719 and three additional variants to accurately infer ABO blood types. We found that A, AB, and B blood types were associated with testing positive for COVID-19 compared with O blood type with the highest risk for the A blood group. We did not observe any genome-wide significant associations for COVID-19 severity outcomes among those testing positive. Our study replicates prior GWAS findings associated with testing positive for COVID-19 among mostly White samples and extends findings at three loci to Black and Hispanic individuals. We also report a new locus among Hispanics requiring further investigation. These findings may aid in the identification of novel therapeutic agents to decrease the morbidity and mortality of COVID-19 across all major ancestral populations.

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.

Background The COVID-19 pandemic remains a significant global threat. However, despite urgent need, there remains uncertainty surrounding best practices for pharmaceutical interventions to treat COVID-19. In particular, conflicting evidence has emerged surrounding the use of hydroxychloroquine and azithromycin, alone or in combination, for COVID-19. The COVID-19 Evidence Accelerator convened by the Reagan-Udall Foundation for the FDA, in collaboration with Friends of Cancer Research, assembled experts from the health systems research, regulatory science, data science, and epidemiology to participate in a large parallel analysis of different data sets to further explore the effectiveness of these treatments. Methods Electronic health record (EHR) and claims data were extracted from seven separate databases. Parallel analyses were undertaken on data extracted from each source. Each analysis examined time to mortality in hospitalized patients treated with hydroxychloroquine, azithromycin, and the two in combination as compared to patients not treated with either drug. Cox proportional hazards models were used, and propensity score methods were undertaken to adjust for confounding. Frequencies of adverse events in each treatment group were also examined. Results Neither hydroxychloroquine nor azithromycin, alone or in combination, were significantly associated with time to mortality among hospitalized COVID-19 patients. No treatment groups appeared to have an elevated risk of adverse events. Conclusion Administration of hydroxychloroquine, azithromycin, and their combination appeared to have no effect on time to mortality in hospitalized COVID-19 patients. Continued research is needed to clarify best practices surrounding treatment of 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 (MR) 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=1.9×10−14) using cis-eQTL genetic instruments that also had strong evidence for colocalization with COVID-19 hospitalization. To disentangle the shared eQTL 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.

Individuals vary in their liability to develop Posttraumatic Stress Disorder (PTSD), the symptoms of which are highly heterogeneous, following exposure to life-threatening trauma. Understanding genetic factors that contribute to the biology of PTSD is critical for refining diagnosis and developing new treatments. Using genetic data from more than 250,000 participants in the Million Veteran Program, genomewide association analyses were conducted using a validated electronic health record-based algorithmically-defined PTSD diagnosis phenotype (48,221 cases and 217,223 controls), and PTSD quantitative symptom phenotypes (212,007 individuals). We identified several genome-wide significant loci in the case-control analyses, and numerous such loci in the quantitative trait analyses, including some (e.g., MAD1L1; TCF4; CRHR1) that were associated with multiple symptom sub-domains and total symptom score, and others that were more specific to certain symptom sub-domains (e.g., CAMKV to re-experiencing; SOX6 to hyperarousal). Genetic correlations between all pairs of symptom sub-domains and their total were very high (rg 0.93 – 0.98) supporting validity of the PTSD diagnostic construct. We also demonstrate strong shared heritability with a range of traits, show that heritability persists when conditioned on other major psychiatric disorders, present independent replication results, provide support for one of the implicated genes in postmortem brain of individuals with PTSD, and use this information to identify potential drug repositioning candidates. These results point to the utility of genetics to inform and validate the biological coherence of the PTSD syndrome despite considerable heterogeneity at the symptom level, and to provide new directions for treatment development.