Search

ImportancePrimary prevention of atherosclerotic cardiovascular disease (ASCVD) relies on risk stratification. Genome-wide polygenic risk scores (PRSs) are proposed to improve ASCVD risk estimation.ObjectiveTo determine whether genome-wide PRSs for coronary artery disease (CAD) and acute ischemic stroke improve ASCVD risk estimation with traditional clinical risk factors in an ancestrally diverse midlife population.Design, Setting, and ParticipantsThis was a prognostic analysis of incident events in a retrospectively defined longitudinal cohort conducted from January 1, 2011, to December 31, 2018. Included in the study were adults free of ASCVD and statin naive at baseline from the Million Veteran Program (MVP), a mega biobank with genetic, survey, and electronic health record data from a large US health care system. Data were analyzed from March 15, 2021, to January 5, 2023.ExposuresPRSs for CAD and ischemic stroke derived from cohorts of largely European descent and risk factors, including age, sex, systolic blood pressure, total cholesterol, high-density lipoprotein (HDL) cholesterol, smoking, and diabetes status.Main Outcomes and MeasuresIncident nonfatal myocardial infarction (MI), ischemic stroke, ASCVD death, and composite ASCVD events.ResultsA total of 79 151 participants (mean [SD] age, 57.8 [13.7] years; 68 503 male [86.5%]) were included in the study. The cohort included participants from the following harmonized genetic ancestry and race and ethnicity categories: 18 505 non-Hispanic Black (23.4%), 6785 Hispanic (8.6%), and 53 861 non-Hispanic White (68.0%) with a median (5th-95th percentile) follow-up of 4.3 (0.7-6.9) years. From 2011 to 2018, 3186 MIs (4.0%), 1933 ischemic strokes (2.4%), 867 ASCVD deaths (1.1%), and 5485 composite ASCVD events (6.9%) were observed. CAD PRS was associated with incident MI in non-Hispanic Black (hazard ratio [HR], 1.10; 95% CI, 1.02-1.19), Hispanic (HR, 1.26; 95% CI, 1.09-1.46), and non-Hispanic White (HR, 1.23; 95% CI, 1.18-1.29) participants. Stroke PRS was associated with incident stroke in non-Hispanic White participants (HR, 1.15; 95% CI, 1.08-1.21). A combined CAD plus stroke PRS was associated with ASCVD deaths among non-Hispanic Black (HR, 1.19; 95% CI, 1.03-1.17) and non-Hispanic (HR, 1.11; 95% CI, 1.03-1.21) participants. The combined PRS was also associated with composite ASCVD across all ancestry groups but greater among non-Hispanic White (HR, 1.20; 95% CI, 1.16-1.24) than non-Hispanic Black (HR, 1.11; 95% CI, 1.05-1.17) and Hispanic (HR, 1.12; 95% CI, 1.00-1.25) participants. Net reclassification improvement from adding PRS to a traditional risk model was modest for the intermediate risk group for composite CVD among men (5-year risk >3.75%, 0.38%; 95% CI, 0.07%-0.68%), among women, (6.79%; 95% CI, 3.01%-10.58%), for age older than 55 years (0.25%; 95% CI, 0.03%-0.47%), and for ages 40 to 55 years (1.61%; 95% CI, −0.07% to 3.30%).Conclusions and RelevanceStudy results suggest that PRSs derived predominantly in European samples were statistically significantly associated with ASCVD in the multiancestry midlife and older-age MVP cohort. Overall, modest improvement in discrimination metrics were observed with addition of PRSs to traditional risk factors with greater magnitude in women and younger age groups.

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.

High dimensional predictive models of Major Adverse Cardiac Events (MACE), which includes heart attack (AMI), stroke, and death caused by cardiovascular disease (CVD), were built using four longitudinal cohorts of Veterans Administration (VA) patients created from VA medical records. We considered 247 variables / risk factors measured across 7.5 years for millions of patients in order to compare predictions for the first reported MACE event using six distinct modelling methodologies. The best-performing methodology varied across the four cohorts. Model coefficients related to disease pathophysiology and treatment were relatively constant across cohorts, while coefficients dependent upon the confounding variables of age and healthcare utilization varied considerably across cohorts. In particular, models trained on a retrospective case-control (Rcc) cohort (where controls are matched to cases by date of birth cohort and overall level of healthcare utilization) emphasize variables describing pathophysiology and treatment, while predictions based on the cohort of all active patients at the start of 2017 (C-17) rely much more on age and variables reflecting healthcare utilization. In consequence, directly using an Rcc-trained model to evaluate the C-17 cohort resulted in poor performance (C-statistic = 0.65). However, a simple reoptimization of model dependence on age, demographics, and five other variables improved the C-statistic to 0.74, nearly matching the 0.76 obtained on C-17 by a C-17-trained model. Dependence of MACE risk on biomarkers for hypertension, cholesterol, diabetes, body mass index, and renal function in our models was consistent with the literature. At the same time, including medications and procedures provided important indications of both disease severity and the level of treatment. More detailed study designs will be required to disentangle these effects.

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 Genetic studies may help identify causal pathways; therefore, we sought to identify genetic determinants of ideal CVH and their association with CVD outcomes in the multi-population Veteran Administration Million Veteran Program. Methods An ideal health score (IHS) was calculated from 3 clinical factors (blood pressure, total cholesterol, and blood glucose levels) and 3 behavioral factors (smoking status, physical activity, and BMI), ascertained at baseline. Multi-population genome-wide association study (GWAS) was performed on IHS and binary ideal health using linear and logistic regression, respectively. Using the genome-wide significant SNPs from the IHS GWAS, we created a weighted IHS polygenic risk score (PRSIHS) which was used (i) to conduct a phenome-wide association study (PheWAS) of associations between PRSIHS and ICD-9 phenotypes and (ii) to further test for associations with mortality and selected CVD outcomes using logistic and Cox regression and, as an instrumental variable, in Mendelian Randomization. Results The discovery and replication cohorts consisted of 142,404 (119,129 European American (EUR); 16,495 African American (AFR)), and 45,766 (37,646 EUR; 5,366 AFR) participants, respectively. The mean age was 65.8 years (SD = 11.2) and 92.7% were male. Overall, 4.2% exhibited ideal CVH based on the clinical and behavioral factors. In the multi-population meta-analysis, variants at 17 loci were associated with IHS and each had known GWAS associations with multiple components of the IHS. PheWAS analysis in 456,026 participants showed that increased PRSIHS was associated with a lower odds ratio for many CVD outcomes and risk factors. Both IHS and PRSIHS measures of ideal CVH were associated with significantly less CVD outcomes and CVD mortality. Conclusion A set of high interest genetic variants contribute to the presence of ideal CVH in a multi-ethnic cohort of US Veterans. Genetically influenced ideal CVH is associated with lower odds of CVD outcomes and mortality.

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.

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.