Design and Pilot Results from Million Veteran Program Return Of Actionable Genetic Results (MVP-ROAR) Study

BACKGROUND

As a mega-biobank linked to a national healthcare system, the Million Veteran Program (MVP) can directly improve the health care of participants. To determine the feasibility and outcomes of returning medically actionable genetic results to MVP participants, the program launched the MVP Return Of Actionable Results (MVP-ROAR) Study, with familial hypercholesterolemia (FH) as an exemplar actionable condition.

METHODS

The MVP-ROAR Study consists of a completed single-arm pilot phase and an ongoing randomized clinical trial (RCT), in which MVP participants are recontacted and invited to receive clinical confirmatory gene sequencing testing and a telegenetic counseling intervention. The primary outcome of the RCT is 6-month change in low-density lipoprotein cholesterol (LDL-C) between participants receiving results at baseline and those receiving results after 6 months.

RESULTS

The pilot developed processes to identify and recontact participants nationally with probable pathogenic variants in low-density lipoprotein receptor (LDLR) on the MVP genotype array, invite them to clinical confirmatory gene sequencing, and deliver a telegenetic counseling intervention. Among participants in the pilot phase, 8 (100%) had active statin prescriptions after 6 months. Results were shared with 16 first-degree family members. Six-month ΔLDL-C (low-density lipoprotein cholesterol) after the genetic counseling intervention was -37 mg/dL (95% CI: -12 to -61; p=0.03). The ongoing RCT will determine between-arm differences in this primary outcome.

CONCLUSION

While underscoring the importance of clinical confirmation of research results, the pilot phase of the MVP-ROAR Study marks a turning point in MVP and demonstrates the feasibility of returning genetic results to participants and their providers. The ongoing RCT will contribute to understanding how such a program might improve patient health care and outcomes.

Characterization of Post-COVID-19 Definitions and Clinical Coding Practices: Longitudinal Study

BACKGROUND

Post-COVID-19 condition (colloquially known as "long COVID-19") characterized as postacute sequelae of SARS-CoV-2 has no universal clinical case definition. Recent efforts have focused on understanding long COVID-19 symptoms, and electronic health record (EHR) data provide a unique resource for understanding this condition. The introduction of the International Classification of Diseases, Tenth Revision (ICD-10) code U09.9 for "Post COVID-19 condition, unspecified" to identify patients with long COVID-19 has provided a method of evaluating this condition in EHRs; however, the accuracy of this code is unclear.

OBJECTIVE

This study aimed to characterize the utility and accuracy of the U09.9 code across 3 health care systems-the Veterans Health Administration, the Beth Israel Deaconess Medical Center, and the University of Pittsburgh Medical Center-against patients identified with long COVID-19 via a chart review by operationalizing the World Health Organization (WHO) and Centers for Disease Control and Prevention (CDC) definitions.

METHODS

Patients who were COVID-19 positive with either a U07.1 ICD-10 code or positive polymerase chain reaction test within these health care systems were identified for chart review. Among this cohort, we sampled patients based on two approaches: (1) with a U09.9 code and (2) without a U09.9 code but with a new onset long COVID-19-related ICD-10 code, which allows us to assess the sensitivity of the U09.9 code. To operationalize the long COVID-19 definition based on health agency guidelines, symptoms were grouped into a "core" cluster of 11 commonly reported symptoms among patients with long COVID-19 and an extended cluster that captured all other symptoms by disease domain. Patients having ≥2 symptoms persisting for ≥60 days that were new onset after their COVID-19 infection, with ≥1 symptom in the core cluster, were labeled as having long COVID-19 per chart review. The code's performance was compared across 3 health care systems and across different time periods of the pandemic.

RESULTS

Overall, 900 patient charts were reviewed across 3 health care systems. The prevalence of long COVID-19 among the cohort with the U09.9 ICD-10 code based on the operationalized WHO definition was between 23.2% and 62.4% across these health care systems. We also evaluated a less stringent version of the WHO definition and the CDC definition and observed an increase in the prevalence of long COVID-19 at all 3 health care systems.

CONCLUSIONS

This is one of the first studies to evaluate the U09.9 code against a clinical case definition for long COVID-19, as well as the first to apply this definition to EHR data using a chart review approach on a nationwide cohort across multiple health care systems. This chart review approach can be implemented at other EHR systems to further evaluate the utility and performance of the U09.9 code.

The Role of X Chromosome in Alzheimer's Disease Genetics

Importance

The X chromosome has remained enigmatic in Alzheimer's disease (AD), yet it makes up 5% of the genome and carries a high proportion of genes expressed in the brain, making it particularly appealing as a potential source of unexplored genetic variation in AD.

Objectives

Perform the first large-scale X chromosome-wide association study (XWAS) of AD. Primary analyses are non-stratified, while secondary analyses evaluate sex-stratified effects.

Design

Meta-analysis of genetic association studies in case-control, family-based, population-based, and longitudinal AD-related cohorts from the US Alzheimer's Disease Genetics Consortium (ADGC) and Alzheimer's Disease Sequencing Project (ADSP), the UK Biobank (UKB), the Finnish health registry (FinnGen), and the US Million Veterans Program (MVP). Risk for AD evaluated through case-control logistic regression analyses. Data were analyzed between January 2023 and March 2024.

Setting

Genetic data available from high-density single-nucleotide polymorphism (SNP) microarrays and whole-genome sequencing (WGS). Summary statistics for multi-tissue expression and protein quantitative trait loci (QTL) available from published studies, enabling follow-up genetic colocalization analyses.

Participants

1,629,863 eligible participants were selected from referred and volunteer samples, of which 477,596 were excluded for analysis exclusion criteria. Number of participants who declined to participate in original studies was not available.

Main Outcome and Measures

Risk for AD (odds ratio; OR) with 95% confidence intervals (CI). Associations were considered at X-chromosome-wide (P-value<1e-5) and genome-wide (P-value<5e-8) significance.

Results

Analyses included 1,152,284 non-Hispanic White European ancestry subjects (57.3% females), including 138,558 cases. 6 independent genetic loci passed X-chromosome-wide significance, with 4 showing support for causal links between the genetic signal for AD and expression of nearby genes in brain and non-brain tissues. One of these 4 loci passed conservative genome-wide significance, with its lead variant centered on an intron of SLC9A7 (OR=1.054, 95%-CI=[1.035, 1.075]) and colocalization analyses prioritizing both the SLC9A7 and nearby CHST7 genes.

Conclusion and Relevance

We performed the first large-scale XWAS of AD and identified the novel SLC9A7 locus. SLC9A7 regulates pH homeostasis in Golgi secretory compartments and is anticipated to have downstream effects on amyloid beta accumulation. Overall, this study significantly advances our knowledge of AD genetics and may provide novel biological drug targets.

Centralized Interactive Phenomics Resource: an integrated online phenomics knowledgebase for health data users

OBJECTIVE

Development of clinical phenotypes from electronic health records (EHRs) can be resource intensive. Several phenotype libraries have been created to facilitate reuse of definitions. However, these platforms vary in target audience and utility. We describe the development of the Centralized Interactive Phenomics Resource (CIPHER) knowledgebase, a comprehensive public-facing phenotype library, which aims to facilitate clinical and health services research.

MATERIALS AND METHODS

The platform was designed to collect and catalog EHR-based computable phenotype algorithms from any healthcare system, scale metadata management, facilitate phenotype discovery, and allow for integration of tools and user workflows. Phenomics experts were engaged in the development and testing of the site.

RESULTS

The knowledgebase stores phenotype metadata using the CIPHER standard, and definitions are accessible through complex searching. Phenotypes are contributed to the knowledgebase via webform, allowing metadata validation. Data visualization tools linking to the knowledgebase enhance user interaction with content and accelerate phenotype development.

DISCUSSION

The CIPHER knowledgebase was developed in the largest healthcare system in the United States and piloted with external partners. The design of the CIPHER website supports a variety of front-end tools and features to facilitate phenotype development and reuse. Health data users are encouraged to contribute their algorithms to the knowledgebase for wider dissemination to the research community, and to use the platform as a springboard for phenotyping.

CONCLUSION

CIPHER is a public resource for all health data users available at https://phenomics.va.ornl.gov/ which facilitates phenotype reuse, development, and dissemination of phenotyping knowledge.

A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Heterogeneous associations between interleukin-6 receptor variants and phenotypes across ancestries and implications for therapy

The Phenome-Wide Association Study (PheWAS) is increasingly used to broadly screen for potential treatment effects, e.g., IL6R variant as a proxy for IL6R antagonists. This approach offers an opportunity to address the limited power in clinical trials to study differential treatment effects across patient subgroups. However, limited methods exist to efficiently test for differences across subgroups in the thousands of multiple comparisons generated as part of a PheWAS. In this study, we developed an approach that maximizes the power to test for heterogeneous genotype-phenotype associations and applied this approach to an IL6R PheWAS among individuals of African (AFR) and European (EUR) ancestries. We identified 29 traits with differences in IL6R variant-phenotype associations, including a lower risk of type 2 diabetes in AFR (OR 0.96) vs EUR (OR 1.0, p-value for heterogeneity = 8.5 × 10-3), and higher white blood cell count (p-value for heterogeneity = 8.5 × 10-131). These data suggest a more salutary effect of IL6R blockade for T2D among individuals of AFR vs EUR ancestry and provide data to inform ongoing clinical trials targeting IL6 for an expanding number of conditions. Moreover, the method to test for heterogeneity of associations can be applied broadly to other large-scale genotype-phenotype screens in diverse populations.

Development and utility of a clinical research informatics application for participant recruitment and workflow management for a return of results pilot trial in familial hypercholesterolemia in the Million Veteran Program

Objective

The development of clinical research informatics tools and workflow processes associated with re-engaging biobank participants has become necessary as genomic repositories increasingly consider the return of actionable research results.

Materials and Methods

Here we describe the development and utility of an informatics application for participant recruitment and enrollment management for the Veterans Affairs Million Veteran Program Return Of Actionable Results Study, a randomized controlled pilot trial returning individual genetic results associated with familial hypercholesterolemia.

Results

The application is developed in Python-Flask and was placed into production in November 2021. The application includes modules for chart review, medication reconciliation, participant contact and biospecimen logging, survey recording, randomization, and documentation of genetic counseling and result disclosure. Three primary users, a genetic counselor and two research coordinators, and 326 Veteran participants have been integrated into the system as of February 23, 2023. The application has successfully handled 3367 task requests involving greater than 95 000 structured data points. Specifically, application users have recorded 326 chart reviews, 867 recruitment telephone calls, 158 telephone-based surveys, and 61 return of results genetic counseling sessions, among other available study tasks.

Conclusion

The development of usable, customizable, and secure informatics tools will become increasingly important as large genomic repositories begin to return research results at scale. Our work provides a proof-of-concept for developing and using such tools to aid in managing the return of results process within a national biobank.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed health problems, prior pharmacological treatments, and polygenic scores (PGS) has potential to inform risk stratification. We examined self-reported SB and ideation using the Columbia Suicide Severity Rating Scale (C-SSRS) among 3,942 SCZ and 5,414 BPI patients receiving care within the Veterans Health Administration (VHA). These cross-sectional data were integrated with electronic health records (EHRs), and compared across lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. PGS were constructed using available genomic data for related traits. Genome-wide association studies were performed to identify and prioritize specific loci. Only 20% of the veterans who reported SB had a corroborating ICD-9/10 EHR code. Among those without prior SB, more than 20% reported new-onset SB at follow-up. SB were associated with a range of additional clinical diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking initiation, suicide attempt, and major depressive disorder were associated with SB. The GWAS for SB yielded no significant loci. Among individuals with a diagnosed mental illness, self-reported SB were strongly associated with clinical variables across several EHR domains. Analyses point to sequelae of substance-related and psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in health records, underscoring the value of regular screening with direct, in-person assessments, especially among high-risk individuals.

Long-Term Body Mass Index Variability and Adverse Cardiovascular Outcomes

Importance

Body mass index (BMI; calculated as weight in kilograms divided by height in meters squared) is a commonly used estimate of obesity, which is a complex trait affected by genetic and lifestyle factors. Marked weight gain and loss could be associated with adverse biological processes.

Objective

To evaluate the association between BMI variability and incident cardiovascular disease (CVD) events in 2 distinct cohorts.

Design, Setting, and Participants

This cohort study used data from the Million Veteran Program (MVP) between 2011 and 2018 and participants in the UK Biobank (UKB) enrolled between 2006 and 2010. Participants were followed up for a median of 3.8 (5th-95th percentile, 3.5) years. Participants with baseline CVD or cancer were excluded. Data were analyzed from September 2022 and September 2023.

Exposure

BMI variability was calculated by the retrospective SD and coefficient of variation (CV) using multiple clinical BMI measurements up to the baseline.

Main Outcomes and Measures

The main outcome was incident composite CVD events (incident nonfatal myocardial infarction, acute ischemic stroke, and cardiovascular death), assessed using Cox proportional hazards modeling after adjustment for CVD risk factors, including age, sex, mean BMI, systolic blood pressure, total cholesterol, high-density lipoprotein cholesterol, smoking status, diabetes status, and statin use. Secondary analysis assessed whether associations were dependent on the polygenic score of BMI.

Results

Among 92 363 US veterans in the MVP cohort (81 675 [88%] male; mean [SD] age, 56.7 [14.1] years), there were 9695 Hispanic participants, 22 488 non-Hispanic Black participants, and 60 180 non-Hispanic White participants. A total of 4811 composite CVD events were observed from 2011 to 2018. The CV of BMI was associated with 16% higher risk for composite CVD across all groups (hazard ratio [HR], 1.16; 95% CI, 1.13-1.19). These associations were unchanged among subgroups and after adjustment for the polygenic score of BMI. The UKB cohort included 65 047 individuals (mean [SD] age, 57.30 (7.77) years; 38 065 [59%] female) and had 6934 composite CVD events. Each 1-SD increase in BMI variability in the UKB cohort was associated with 8% increased risk of cardiovascular death (HR, 1.08; 95% CI, 1.04-1.11).

Conclusions and Relevance

This cohort study found that among US veterans, higher BMI variability was a significant risk marker associated with adverse cardiovascular events independent of mean BMI across major racial and ethnic groups. Results were consistent in the UKB for the cardiovascular death end point. Further studies should investigate the phenotype of high BMI variability.

Red Meat Intake and the Risk of Cardiovascular Diseases: A Prospective Cohort Study in the Million Veteran Program

BACKGROUND

Red meat consumption was associated with an increased risk of cardiovascular disease (CVD) in prospective cohort studies and a profile of biomarkers favoring high CVD risk in short-term controlled trials. However, several recent systematic reviews and meta-analyses concluded with no or weak evidence for limiting red meat intake.

OBJECTIVES

To prospectively examine the associations between red meat intake and incident CVD in an ongoing cohort study with diverse socioeconomic and racial or ethnic backgrounds.

METHODS

Our study included 148,506 participants [17,804 female (12.0%)] who were free of cancer, diabetes, and CVD at baseline from the Million Veteran Program. A food frequency questionnaire measured red meat intakes at baseline. Nonfatal myocardial infarction and acute ischemic stroke were identified through a high-throughput phenotyping algorithm, and fatal CVD events were identified by searching the National Death Index.

RESULTS

Comparing the extreme categories of intake, the multivariate-adjusted relative risks of CVD was 1.18 (95% CI: 1.01, 1.38; P-trend < 0.0001) for total red meat, 1.14 (95% CI: 0.96, 1.36; P-trend = 0.01) for unprocessed red meat, and 1.29 (95% CI: 1.04, 1.60; P-trend = 0.003) for processed red meat. We observed a more pronounced positive association between red meat intake and CVD in African American participants than in White participants (P-interaction = 0.01). Replacing 0.5 servings/d of red meat with 0.5 servings/d of nuts, whole grains, and skimmed milk was associated with 14% (RR: 0.86; 95% CI: 0.83, 0.90), 7% (RR: 0.93; 95% CI: 0.89, 0.96), and 4% (RR: 0.96; 95% CI: 0.94, 0.99) lower risks of CVD, respectively.

CONCLUSIONS

Red meat consumption is associated with an increased risk of CVD. Our findings support lowering red meat intake and replacing red meat with plant-based protein sources or low-fat dairy foods as a key dietary recommendation for the prevention of CVD.

New statin use, mortality, and first cardiovascular events in older US Veterans by frailty status

BACKGROUND

Statins are part of long-term medical regimens for many older adults. Whether frailty modifies the protective relationship between statins, mortality, and major adverse cardiovascular events (MACE) is unknown.

METHODS

This was a retrospective study of US Veterans ≥65, without CVD or prior statin use seen in 2002-2012, followed through 2017. A 31-item frailty index was used. The co-primary endpoint was all-cause mortality or MACE (MI, stroke/TIA, revascularization, or cardiovascular death). Cox proportional hazards models were developed to evaluate the association of statin use with outcomes; propensity score overlap weighting accounted for confounding by indication.

RESULTS

We identified 710,313 Veterans (mean age (SD) 75.3(6.5), 98% male, 89% white); 86,327 (12.1%) were frail. Over mean follow-up of 8 (5) years, there were 48.6 and 72.6 deaths per 1000 person-years (PY) among non-frail statin-users vs nonusers (weighted Incidence Rate Difference (wIRD)/1000 person years (PY), -24.0[95% CI, -24.5 to -23.6]), and 90.4 and 130.4 deaths per 1000PY among frail statin-users vs nonusers (wIRD/1000PY, -40.0[95% CI, -41.8 to -38.2]). There were 51.7 and 60.8 MACE per 1000PY among non-frail statin-users vs nonusers (wIRD/1000PY, -9.1[95% CI, -9.7 to -8.5]), and 88.2 and 102.0 MACE per 1000PY among frail statin-users vs nonusers (wIRD/1000PY, -13.8[95% CI, -16.2 to -11.4]). There were no significant interactions by frailty for statin users vs non-users by either mortality or MACE outcomes, p-interaction 0.770 and 0.319, respectively. Statin use was associated with lower risk of all-cause mortality (HR, 0.61 (0.60-0.61)) and MACE (HR 0.86 (0.85-0.87)).

CONCLUSIONS

New statin use is associated with a lower risk of mortality and MACE, independent of frailty. These findings should be confirmed in a randomized clinical trial.

Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

Clinical predictors of right ventricular dysfunction and association with adverse outcomes in peripartum cardiomyopathy

AIMS

We sought to identify factors associated with right ventricular (RV) dysfunction and elevated pulmonary artery systolic pressure (PASP) and association with adverse outcomes in peripartum cardiomyopathy (PPCM).

METHODS AND RESULTS

We conducted a multi-centre cohort study to identify subjects with PPCM with the following criteria: left ventricular ejection fraction (LVEF) < 40%, development of heart failure within the last month of pregnancy or 5 months of delivery, and no other identifiable cause of heart failure with reduced ejection fraction. Outcomes included a composite of (i) major adverse events (need for extracorporeal membrane oxygenation, ventricular assist device, orthotopic heart transplantation, or death) or (ii) recurrent heart failure hospitalization. RV function was obtained from echocardiogram reports. In total, 229 women (1993-2017) met criteria for PPCM. Mean age was 32.4 ± 6.8 years, 28% were of African descent, 50 (22%) had RV dysfunction, and 38 (17%) had PASP ≥ 30 mmHg. After a median follow-up of 3.4 years (interquartile range 1.0-8.8), 58 (25%) experienced the composite outcome of adverse events. African descent, family history of cardiomyopathy, LVEF, and PASP were significant predictors of RV dysfunction. Using Cox proportional hazards models, we found that women with RV dysfunction were three times more likely to experience the adverse composite outcome: hazard ratio 3.21 (95% confidence interval: 1.11-9.28), P = 0.03, in a multivariable model adjusting for age, race, body mass index, preeclampsia, hypertension, diabetes, kidney disease, and LVEF. Women with PASP ≥ 30 mmHg had a lower probability of survival free from adverse events (log-rank P = 0.04).

CONCLUSIONS

African descent and family history of cardiomyopathy were significant predictors of RV dysfunction. RV dysfunction and elevated PASP were significantly associated with a composite of major adverse cardiac events. This at-risk group may prompt closer monitoring or early referral for advanced therapies.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

LATTE: Label-efficient incident phenotyping from longitudinal electronic health records

Electronic health record (EHR) data are increasingly used to support real-world evidence studies but are limited by the lack of precise timings of clinical events. Here, we propose a label-efficient incident phenotyping (LATTE) algorithm to accurately annotate the timing of clinical events from longitudinal EHR data. By leveraging the pre-trained semantic embeddings, LATTE selects predictive features and compresses their information into longitudinal visit embeddings through visit attention learning. LATTE models the sequential dependency between the target event and visit embeddings to derive the timings. To improve label efficiency, LATTE constructs longitudinal silver-standard labels from unlabeled patients to perform semi-supervised training. LATTE is evaluated on the onset of type 2 diabetes, heart failure, and relapses of multiple sclerosis. LATTE consistently achieves substantial improvements over benchmark methods while providing high prediction interpretability. The event timings are shown to help discover risk factors of heart failure among patients with rheumatoid arthritis.

The Million Veteran Program 1990-1991 Gulf War Era Survey: An Evaluation of Veteran Response, Characteristics, and Representativeness of the Gulf War Era Veteran Population

To address gaps in understanding the pathophysiology of Gulf War Illness (GWI), the VA Million Veteran Program (MVP) developed and implemented a survey to MVP enrollees who served in the U.S. military during the 1990-1991 Persian Gulf War (GW). Eligible Veterans were invited via mail to complete a survey assessing health conditions as well as GW-specific deployment characteristics and exposures. We evaluated the representativeness of this GW-era cohort relative to the broader population by comparing demographic, military, and health characteristics between respondents and non-respondents, as well as with all GW-era Veterans who have used Veterans Health Administration (VHA) services and the full population of U.S. GW-deployed Veterans. A total of 109,976 MVP GW-era Veterans were invited to participate and 45,270 (41%) returned a completed survey. Respondents were 84% male, 72% White, 8% Hispanic, with a mean age of 61.6 years (SD = 8.5). Respondents were more likely to be older, White, married, better educated, slightly healthier, and have higher socioeconomic status than non-respondents, but reported similar medical conditions and comparable health status. Although generally similar to all GW-era Veterans using VHA services and the full population of U.S. GW Veterans, respondents included higher proportions of women and military officers, and were slightly older. In conclusion, sample characteristics of the MVP GW-era cohort can be considered generally representative of the broader GW-era Veteran population. The sample represents the largest research cohort of GW-era Veterans established to date and provides a uniquely valuable resource for conducting in-depth studies to evaluate health conditions affecting 1990-1991 GW-era Veterans.

Adaption and National Validation of a Tool for Predicting Mortality from Other Causes Among Men with Nonmetastatic Prostate Cancer

BACKGROUND

An electronic health record-based tool could improve accuracy and eliminate bias in provider estimation of the risk of death from other causes among men with nonmetastatic cancer.

OBJECTIVE

To recalibrate and validate the Veterans Aging Cohort Study Charlson Comorbidity Index (VACS-CCI) to predict non-prostate cancer mortality (non-PCM) and to compare it with a tool predicting prostate cancer mortality (PCM).

DESIGN, SETTING, AND PARTICIPANTS

An observational cohort of men with biopsy-confirmed nonmetastatic prostate cancer, enrolled from 2001 to 2018 in the national US Veterans Health Administration (VA), was divided by the year of diagnosis into the development (2001-2006 and 2008-2018) and validation (2007) sets.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Mortality (all cause, non-PCM, and PCM) was evaluated. Accuracy was assessed using calibration curves and C statistic in the development, validation, and combined sets; overall; and by age (<65 and 65+ yr), race (White and Black), Hispanic ethnicity, and treatment groups.

RESULTS AND LIMITATIONS

Among 107 370 individuals, we observed 24 977 deaths (86% non-PCM). The median age was 65 yr, 4947 were Black, and 5010 were Hispanic. Compared with CCI and age alone (C statistic 0.67, 95% confidence interval [CI] 0.67-0.68), VACS-CCI demonstrated improved validated discrimination (C statistic 0.75, 95% CI 0.74-0.75 for non-PCM). The prostate cancer mortality tool also discriminated well in validation (C statistic 0.81, 95% CI 0.78-0.83). Both were well calibrated overall and within subgroups. Owing to missing data, 18 009/125 379 (14%) were excluded, and VACS-CCI should be validated outside the VA prior to outside application.

CONCLUSIONS

VACS-CCI is ready for implementation within the VA. Electronic health record-assisted calculation is feasible, improves accuracy over age and CCI alone, and could mitigate inaccuracy and bias in provider estimation.

PATIENT SUMMARY

Veterans Aging Cohort Study Charlson Comorbidity Index is ready for application within the Veterans Health Administration. Electronic health record-assisted calculation is feasible, improves accuracy over age and Charlson Comorbidity Index alone, and might help mitigate inaccuracy and bias in provider estimation of the risk of non-prostate cancer mortality.

Moderate alcohol consumption on the risk of stroke in the Million Veteran Program

BACKGROUND

There is inconsistent evidence on the association of moderate alcohol consumption and stroke risk in the general population and is not well studied among U.S. Veterans. Furthermore, it is unclear whether primarily drinking beer, wine, or liquor is associated with a difference in stroke risk.

METHODS

The study included 185,323 Million Veteran Program participants who self-reported alcohol consumption on the Lifestyle Survey. Moderate consumption was defined as 1-2 drinks/day and beverage preference of beer, wine or liquor was defined if ≥ 50% of total drinks consumed were from a single type of beverage. Strokes were defined using ICD-9 and ICD-10 codes from the participants' electronic health record.

RESULTS

The mean (sd) age of the sample was 64 (13) years and 11% were women. We observed 4,339 (94% ischemic; 6% hemorrhagic) strokes over a median follow-up of 5.2 years. In Cox models adjusted for age, sex, race, education, income, body mass index, smoking, exercise, diet, cholesterol, prevalent diabetes, prevalent hypertension, lipid-lowering medication, antihypertensive medication, and diabetes medication, moderate alcohol consumption (1-2 drinks/day) was associated with a 22% lower risk of total stroke compared with never drinking [Hazards ratio (HR) 95% confidence interval (CI): 0.78 (0.67, 0.92)]. When stratifying by stroke type, we observed a similar protective association with moderate consumption and ischemic stroke [HR (95% CI): 0.76 (0.65, 0.90)], but a non-statistically significant higher risk of hemorrhagic stroke [HR (95% CI): 1.29 (0.64, 2.61)]. We did not observe a difference in ischemic or hemorrhagic stroke risk among those who preferred beer, liquor or wine vs. no beverage preference. When stratifying by prior number of hospital visits (≤ 15, 16-33, 34-64, ≥ 65) as a proxy for health status, we observed attenuation of the protective association with greater number of visits [HR (95% CI): 0.87 (0.63, 1.19) for ≥ 65 visits vs. 0.80 (0.59, 1.08) for ≤ 15 visits].

CONCLUSIONS

We observed a lower risk of ischemic stroke, but not hemorrhagic stroke with moderate alcohol consumption and did not observe substantial differences in risk by beverage preference among a sample of U.S. Veterans. Healthy user bias of moderate alcohol consumption may be driving some of the observed protective association.

Statins, Mortality, and Major Adverse Cardiovascular Events Among US Veterans With Chronic Kidney Disease

Importance

There are limited data for the utility of statins for primary prevention of atherosclerotic cardiovascular disease (ASCVD) and death in adults with chronic kidney disease (CKD).

Objective

To evaluate the association of statin use with all-cause mortality and major adverse cardiovascular events (MACE) among US veterans older than 65 years with CKD stages 3 to 4.

Design, Setting, and Participants

This cohort study used a target trial emulation design for statin initiation among veterans with moderate CKD (stages 3 or 4) using nested trials with a propensity weighting approach. Linked Veterans Affairs (VA) Healthcare System, Medicare, and Medicaid data were used. This study considered veterans newly diagnosed with moderate CKD between 2005 and 2015 in the VA, with follow-up through December 31, 2017. Veterans were older than 65 years, within 5 years of CKD diagnosis, had no prior ASCVD or statin use, and had at least 1 clinical visit in the year prior to trial baseline. Eligibility criteria were assessed for each nested trial, and Cox proportional hazards models with bootstrapping were run. Analysis was conducted from July 2021 to October 2023.

Exposure

Statin initiation vs none.

Main Outcomes and Measures

Primary outcome was all-cause mortality; secondary outcome was time to first MACE (myocardial infarction, transient ischemic attack, stroke, revascularization, or mortality).

Results

Included in the analysis were 14 828 veterans. Mean (SD) age at CKD diagnosis was 76.9 (8.2) years, 14 616 (99%) were men, 10 539 (72%) White, and 2568 (17%) Black. After expanding to person-trials and assessing eligibility at each baseline, there were 151 243 person-trials (14 685 individuals) of nonstatin initiators and 2924 person-trials (2924 individuals) of statin initiators included. Propensity score adjustment via overlap weighting with nonparametric bootstrapping resulted in covariate balance, with mean (SD) follow-up of 3.6 (2.7) years. The hazard ratio for all-cause mortality was 0.91 (95% CI, 0.85-0.97) comparing statin initiators to noninitiators. The hazard ratio for MACE was 0.96 (95% CI, 0.91-1.02). Results remained consistent in prespecified subgroup analyses.

Conclusions and Relevance

In this target trial emulation of statin initiation in US veterans older than 65 years with CKD stages 3 to 4 and no prior ASCVD, statin initiation was significantly associated with a lower risk of all-cause mortality but not MACE. Results should be confirmed in a randomized clinical trial.

Multi-ancestry genome-wide association study of cannabis use disorder yields insight into disease biology and public health implications

As recreational use of cannabis is being decriminalized in many places and medical use widely sanctioned, there are growing concerns about increases in cannabis use disorder (CanUD), which is associated with numerous medical comorbidities. Here we performed a genome-wide association study of CanUD in the Million Veteran Program (MVP), followed by meta-analysis in 1,054,365 individuals (ncases = 64,314) from four broad ancestries designated by the reference panel used for assignment (European n = 886,025, African n = 123,208, admixed American n = 38,289 and East Asian n = 6,843). Population-specific methods were applied to calculate single nucleotide polymorphism-based heritability within each ancestry. Statistically significant single nucleotide polymorphism-based heritability for CanUD was observed in all but the smallest population (East Asian). We discovered genome-wide significant loci unique to each ancestry: 22 in European, 2 each in African and East Asian, and 1 in admixed American ancestries. A genetically informed causal relationship analysis indicated a possible effect of genetic liability for CanUD on lung cancer risk, suggesting potential unanticipated future medical and psychiatric public health consequences that require further study to disentangle from other known risk factors such as cigarette smoking.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Association of potato consumption with calcified atherosclerotic plaques in the coronary arteries: The NHLBI Family Heart Study

BACKGROUND AND AIMS

While the association of potato consumption with risk factors for coronary artery disease has been inconsistent, no data are available in the literature on the influence of potato consumption on subclinical disease. Thus, we sought to examine whether baked/mashed potato consumption is associated with calcified atherosclerotic plaques in the coronary arteries.

METHODS AND RESULTS

In a cross-sectional design, we studied 2208 participants of the NHLBI Family Heart Study. These subjects were selected based on their elevated cardiovascular disease risk compared to the general population. Potato consumption was assessed by a semi-quantitative food frequency questionnaire. We defined prevalent CAC using an Agatston score of at least 100 and fitted generalized estimating equations to calculate prevalence odds ratios of CAC. Mean age at initial clinic visit was 58.2 years and 55% were female. Median consumption of potatoes was 2-4/week. There was no statistically significant association between frequency of potato consumption and prevalent CAC: odds ratios (95% CI) for CAC were 1.0 (reference), 0.85 (0.56-1.30), 0.85 (0.58-1.26), and 0.95 (0.60-1.53) among subjects reporting potato consumption of <1/week, 1/week, 2-4/week, and 5+/week, respectively (p for linear trend 0.83), adjusting for age, sex, BMI, smoking, exercise, diabetes, hypertension, total calories, prevalent coronary heart disease, income, education, and daily red meat intake.

CONCLUSIONS

We found no significant association between baked/mashed potato consumption and CAC in older adults.

STUDY REGISTRATION NUMBER

NCT00005136. Study registration date: 5/25/2000.

Serum Cholesterol and Impact of Age on Coronary Heart Disease Death in More Than 4 Million Veterans

Background The lipid hypothesis postulates that lower blood cholesterol is associated with reduced coronary heart disease (CHD) risk, which has been challenged by reports of a U-shaped relation between cholesterol and death in recent studies. We sought to examine whether the U-shaped relationship is true and to assess the impact of age on this association. Method and Results We conducted a prospective cohort study of 4 467 942 veterans aged >18 years, with baseline outpatient visits from 2002 to 2007 and follow-up to December 30, 2018, in the Veterans Health Administration electronic health record system. We observed a J-shaped relation between total cholesterol (TC) and CHD mortality after a comprehensive adjustment of confounding factors: flat for TC <180 mg/dL, and greater risk was present at higher cholesterol levels. Compared with veterans with TC between 180 and 199 mg/dL, the multiadjusted hazard ratios (HRs) for CHD death were 1.03 (95% CI, 1.02-1.04), 1.07 (95% CI, 1.06-1.09), 1.15 (95% CI, 1.13-1.18), 1.25 (95% CI, 1.22-1.28), and 1.45 (95% CI, 1.42-1.49) times greater among veterans with TC (mg/dL) of 200 to 219, 220 to 239, 140 to 259, 260 to 279 and ≥280, respectively. Similar J-shaped TC-CHD mortality patterns were observed among veterans with and without statin use at or before baseline. Conclusions The cholesterol paradox, for example, higher CHD death in patients with a low cholesterol level, was a reflection of reverse causality, especially among older participants. Our results support the lipid hypothesis that lower blood cholesterol is associated with reduced CHD. Furthermore, the hypothesis remained true when TC was low due to use of statins or other lipid-lowering medication.

Genetic Inhibition of APOL1 Pore-Forming Function Prevents APOL1-Mediated Kidney Disease

SIGNIFICANCE STATEMENT

African Americans are at increased risk of CKD in part due to high-risk (HR) variants in the apolipoprotein L1 ( APOL1 ) gene, termed G1/G2. A different APOL1 variant, p.N264K , reduced the risk of CKD and ESKD among carriers of APOL1 HR variants to levels comparable with individuals with APOL1 low-risk variants in an analysis of 121,492 participants of African ancestry from the Million Veteran Program (MVP). Functional genetic studies in cell models showed that APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR mutations. Pharmacologic inhibitors that mimic this mutation blocking APOL1 -mediated pore formation may be able to prevent and/or treat APOL1 -associated kidney disease.

BACKGROUND

African Americans are at increased risk for nondiabetic CKD in part due to HR variants in the APOL1 gene.

METHODS

We tested whether a different APOL1 variant, p.N264K , modified the association between APOL1 HR genotypes (two copies of G1/G2) and CKD in a cross-sectional analysis of 121,492 participants of African ancestry from the MVP. We replicated our findings in the Vanderbilt University Biobank ( n =14,386) and National Institutes of Health All of Us ( n =14,704). Primary outcome was CKD and secondary outcome was ESKD among nondiabetic patients. Primary analysis compared APOL1 HR genotypes with and without p.N264K . Secondary analyses included APOL1 low-risk genotypes and tested for interaction. In MVP, we performed sequential logistic regression models adjusting for demographics, comorbidities, medications, and ten principal components of ancestry. Functional genomic studies expressed APOL1 HR variants with and without APOL1 p.N264K in cell models.

RESULTS

In the MVP cohort, 15,604 (12.8%) had two APOL1 HR variants, of which 582 (0.5%) also had APOL1 p.N264K . In MVP, 18,831 (15%) had CKD, 4177 (3%) had ESKD, and 34% had diabetes. MVP APOL1 HR, without p.N264K , was associated with increased odds of CKD (odds ratio [OR], 1.72; 95% confidence interval [CI], 1.60 to 1.85) and ESKD (OR, 3.94; 95% CI, 3.52 to 4.41). In MVP, APOL1 p.N264K mitigated the renal risk of APOL1 HR, in CKD (OR, 0.43; 95% CI, 0.28 to 0.65) and ESKD (OR, 0.19; CI 0.07 to 0.51). In the replication cohorts meta-analysis, APOL1 p.N264K mitigated the renal risk of APOL1 HR in CKD (OR, 0.40; 95% CI, 0.18 to 0.92) and ESKD (OR, 0.19; 95% CI, 0.05 to 0.79). In the mechanistic studies, APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR variants.

CONCLUSIONS

APOL1 p.N264K is associated with reduced risk of CKD and ESKD among carriers of APOL1 HR to levels comparable with individuals with APOL1 low-risk genotypes.

Leveraging the Million Veteran Program Infrastructure and Data for a Rapid Research Response to COVID-19

Background

The Veterans Health Administration Office of Research and Development (ORD) played a key role in the federal government's response to the COVID-19 pandemic. The ORD effectively leveraged existing resources to answer questions related to the SARS-CoV-2 virus and COVID-19.

Observations

When the COVID-19 pandemic hit in 2020, the Million Veteran Program (MVP), one of the largest genomic cohorts in the world, extended the centralized recruitment and enrollment infrastructure to develop a COVID-19 research volunteer registry to assist enrollment in the vaccine and treatment trials in which the US Department of Veterans Affairs (VA) participated. In addition, the MVP allowed for new data collection and a large genomic cohort to understand host contributions to COVID-19. This article describes ways the MVP contributed to the VA's rapid research response to COVID-19. Several host genetic factors believed to play a role in the development and severity of COVID-19 were identified. Furthermore, existing MVP partnerships with other federal agencies, particularly with the Department of Energy, were leveraged to improve understanding and management of COVID-19.

Conclusions

A previously established enterprise approach and research infrastructure were essential to the VA's successful and timely COVID-19 research response. This infrastructure not only supported rapid recruitment in vaccine and treatment trials, but also leveraged the unique MVP and VA electronic health record data to drive rapid scientific discovery and inform clinical operations. Extending the models that VA research applied to the federal government at large and establishing centralized resources for shared or federated data analyses across federal agencies will better equip the nation to respond to future public health crises.

Comparative effectiveness and cost-effectiveness of cardioprotective glucose-lowering therapies for type 2 diabetes in Brazil: a Bayesian network model

BACKGROUND

The escalating prevalence of type 2 diabetes (T2DM) poses an unparalleled economic catastrophe to developing countries. Cardiovascular diseases remain the primary source of costs among individuals with T2DM, incurring expenses for medications, hospitalizations, and surgical interventions. Compelling evidence suggests that the risk of cardiovascular outcomes can be reduced by three classes of glucose-lowering therapies (GLT), including SGLT2i, GLP-1A, and pioglitazone. However, an evidence-based and cost-effective protocol is still unavailable for many countries. The objective of the current study is to compare the effectiveness and cost-effectiveness of GLT in individuals with T2DM in Brazil.

METHODS

We employed Bayesian Networks to calculate the incremental cost-effectiveness ratios (ICER), expressed in international dollars (Int$) per disease-adjusted life years [DALYs] averted. To determine the effectiveness of GLT, we conducted a systematic review with network meta-analysis (NMA) to provide insights for our model. Additionally, we obtained cardiovascular outcome incidence data from two real-world cohorts comprising 851 and 1337 patients in primary and secondary prevention, respectively. Our cost analysis took into account the perspective of the Brazilian public health system, and all values were converted to Int$.

RESULTS

In the NMA, SGLT2i [HR: 0.81 (95% CI 0.69-0.96)], GLP-1A [HR: 0.79 (95% CI 0.67-0.94)], and pioglitazone [HR: 0.73 (95% CI 0.59-0.91)] demonstrated reduced relative risks of non-fatal cardiovascular events. In the context of primary prevention, pioglitazone yielded 0.2339 DALYs averted, with an ICER of Int$7,082 (95% CI 4,521-10,770) per DALY averted when compared to standard care. SGLT2i and GLP-1A also increased effectiveness, resulting in 0.261 and 0.259 DALYs averted, respectively, but with higher ICERs of Int$12,061 (95% CI: 7,227-18,121) and Int$29,119 (95% CI: 23,811-35,367) per DALY averted. In the secondary prevention scenario, all three classes of treatments were deemed cost-effective at a maximum willingness-to-pay threshold of Int$26,700. Notably, pioglitazone consistently exhibited the highest probability of being cost-effective in both scenarios.

CONCLUSIONS

In Brazil, pioglitazone presented a higher probability of being cost-effective both in primary and secondary prevention, followed by SGLT2i and GLP-1A. Our findings support the use of cost-effectiveness models to build optimized and hierarchical therapeutic strategy in the management of T2DM.

TRIAL REGISTRATION

CRD42020194415.

More intensive therapy as more effective treatment for frail patients with multiple myeloma [corrected]

Although randomized controlled trial data suggest that the more intensive triplet bortezomib-lenalidomide-dexamethasone (VRd) is superior to the less intensive doublet lenalidomide-dexamethasone (Rd) in patients newly diagnosed with multiple myeloma (MM), guidelines have historically recommended Rd over VRd for patients who are frail and may not tolerate a triplet. We identified 2573 patients (median age, 69.7 years) newly diagnosed with MM who were initiated on VRd (990) or Rd (1583) in the national US Veterans Affairs health care System from 2004 to 2020. We measured frailty using the Veterans Affairs Frailty Index. To reduce imbalance in confounding, we matched patients for MM stage and 1:1 based on a propensity score. Patients who were moderate-severely frail had a higher prevalence of stage III MM and myeloma-related frailty deficits than patients who were not frail. VRd vs Rd was associated with lower mortality (hazard ratio [HR], 0.81; 95% confidence interval [CI], 0.70-0.94) in the overall matched population. Patients who were moderate-severely frail demonstrated the strongest association (HR 0.74; 95% CI, 0.56-0.97), whereas the association weakened in those who were mildly frail (HR, 0.80; 95% CI, 0.61-1.05) and nonfrail (HR, 0.86; 95% CI, 0.67-1.10). VRd vs Rd was associated with a modestly higher incidence of hospitalizations in the overall population, but this association weakened in patients who were moderate-severely frail. Our findings confirm the benefit of VRd over Rd in US veterans and further suggest that this benefit is strongest in patients with the highest levels of frailty, arguing that more intensive treatment of myeloma may be more effective treatment of frailty itself.

Knowledge-Driven Online Multimodal Automated Phenotyping System

Though electronic health record (EHR) systems are a rich repository of clinical information with large potential, the use of EHR-based phenotyping algorithms is often hindered by inaccurate diagnostic records, the presence of many irrelevant features, and the requirement for a human-labeled training set. In this paper, we describe a knowledge-driven online multimodal automated phenotyping (KOMAP) system that i) generates a list of informative features by an online narrative and codified feature search engine (ONCE) and ii) enables the training of a multimodal phenotyping algorithm based on summary data. Powered by composite knowledge from multiple EHR sources, online article corpora, and a large language model, features selected by ONCE show high concordance with the state-of-the-art AI models (GPT4 and ChatGPT) and encourage large-scale phenotyping by providing a smaller but highly relevant feature set. Validation of the KOMAP system across four healthcare centers suggests that it can generate efficient phenotyping algorithms with robust performance. Compared to other methods requiring patient-level inputs and gold-standard labels, the fully online KOMAP provides a significant opportunity to enable multi-center collaboration.

Genome-wide association studies and cross-population meta-analyses investigating short and long sleep duration

Sleep duration has been linked to a wide range of negative health outcomes and to reduced life expectancy. We present genome-wide association studies of short ( ≤ 5 h) and long ( ≥ 10 h) sleep duration in adults of European (N = 445,966), African (N = 27,785), East Asian (N = 3141), and admixed-American (N = 16,250) ancestry from UK Biobank and the Million Veteran Programme. In a cross-population meta-analysis, we identify 84 independent loci for short sleep and 1 for long sleep. We estimate SNP-based heritability for both sleep traits in each ancestry based on population derived linkage disequilibrium (LD) scores using cov-LDSC. We identify positive genetic correlation between short and long sleep traits (rg = 0.16 ± 0.04; p = 0.0002), as well as similar patterns of genetic correlation with other psychiatric and cardiometabolic phenotypes. Mendelian randomisation reveals a directional causal relationship between short sleep and depression, and a bidirectional causal relationship between long sleep and depression.

Correction for Collider Bias in the Genome-wide Association Study of Diabetes-Related Heart Failure due to Bidirectional Relationship between Heart Failure and Type 2 Diabetes

Aims

Type 2 diabetes (T2D) is a major risk factor for heart failure (HF) across demographic groups. On the other hand, metabolic impairment, including elevated T2D incidence is a hallmark of HF pathophysiology. We investigated the bidirectional relationship between T2D and HF, and identified genetic associations with diabetes-related HF after correction for potential collider bias.

Methods

We performed a genome-wide association study (GWAS) of HF to identify genetic instrumental variables (GIVs) for HF, and to enable bidirectional Mendelian Randomization (MR) analysis between T2D and HF. Since genetics and HF can independently influence T2D, collider bias may occur when T2D (i.e., collider) is controlled for by design or analysis. Thus, we conducted GWAS of diabetes-related HF with correction for collider bias.

Results

We first identified 61 genomic loci, including 24 novel loci, significantly associated with all-cause HF in 114,275 HF cases and over 1.5 million controls of European ancestry. Combined with the summary statistics of a T2D GWAS, we obtained 59 and 82 GIVs for HF and T2D, respectively. Using a two-sample bidirectional MR approach, we estimated that T2D increased HF risk (OR 1.07, 95% CI 1.04-1.10), while HF also increased T2D risk (OR 1.60, 95% CI 1.36-1.88). Then we performed a GWAS of diabetes-related HF corrected for collider bias due to prevalent HF affecting incidence of T2D. After removing the spurious association of TCF7L2 locus due to collider bias, we identified two genome-wide significant loci close to PITX2 (chromosome 4) and CDKN2B-AS1 (chromosome 9) associated with diabetes-related HF in the Million Veteran Program, and replicated the associations in the UK Biobank study.

Conclusion

We identified novel HF-associated loci to enable bidirectional MR study of T2D and HF. Our MR findings support T2D as a HF risk factor and provide strong evidence that HF increases T2D risk. As a result, collider bias leads to spurious genetic associations of diabetes-related HF, which can be effectively corrected to identify true positive loci. Evaluation of collider bias should be a critical component when conducting GWAS of complex disease phenotypes such as diabetes-related cardiovascular complications.

Predictive value of ASCVD risk score for mortality and major adverse cardiovascular events in the year following a COVID-19 infection among US Veterans

BACKGROUND

Morbidity and mortality following COVID-19 infection may be influenced by baseline atherosclerotic cardiovascular disease (ASCVD) risk, yet limited data are available to identify those at highest risk. We examined the association between baseline ASCVD risk with mortality and major adverse cardiovascular events (MACE) in the year following COVID-19 infection.

METHODS

We evaluated a nationwide retrospective cohort of US Veterans free of ASCVD who were tested for COVID-19. The primary outcome was absolute risk of all-cause mortality in the year following a COVID-19 test among those hospitalized vs. not stratified by baseline VA-ASCVD risk scores. Secondarily, risk of MACE was examined.

RESULTS

There were 393,683 Veterans tested for COVID-19 and 72,840 tested positive. Mean age was 57 years, 86% were male, and 68% were white. Within 30 days following infection, hospitalized Veterans with VA-ASCVD scores >20% had an absolute risk of death of 24.6% vs. 9.7% (P ≤0.0001) for those who tested positive and negative for COVID-19 respectively. In the year following infection, risk of mortality attenuated with no difference in risk after 60 days. The absolute risk of MACE was similar for Veterans who tested positive or negative for COVID-19.

CONCLUSIONS

Veterans without clinical ASCVD experienced an increased absolute risk of death within 30 days of a COVID-19 infection compared to Veterans with the same VA-ASCVD risk score who tested negative, but this risk attenuated after 60 days. Whether cardiovascular preventive medications can lower the risk of mortality and MACE in the acute period following COVID-19 infection should be evaluated.

Genetic Susceptibility to Nonalcoholic Fatty Liver Disease and Risk for Pancreatic Cancer: Mendelian Randomization

BACKGROUND

There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer.

METHODS

Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes.

RESULTS

No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88-1.22; MR-Egger OR, 0.89; 95% CI, 0.65-1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90-1.27; MR-Egger OR, 0.93; 95% CI, 0.67-1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample.

CONCLUSIONS

Genetic predisposition to NAFLD is not associated with pancreatic cancer risk.

IMPACT

Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer.

Genome-Wide Association Study of CKD Progression

SIGNIFICANCE STATEMENT

Rapid progression of CKD is associated with poor clinical outcomes. Most previous studies looking for genetic factors associated with low eGFR have used cross-sectional data. The authors conducted a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD, focusing on longitudinal data. They identified three loci (two of them novel) associated with longitudinal eGFR decline. In addition to the known UMOD/PDILT locus, variants within BICC1 were associated with significant differences in longitudinal eGFR slope. Variants within HEATR4 also were associated with differences in eGFR decline, but only among Black/African American individuals without diabetes. These findings help characterize molecular mechanisms of eGFR decline in CKD and may inform new therapeutic approaches for progressive kidney disease.

BACKGROUND

Rapid progression of CKD is associated with poor clinical outcomes. Despite extensive study of the genetics of cross-sectional eGFR, only a few loci associated with eGFR decline over time have been identified.

METHODS

We performed a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD-defined by two outpatient eGFR measurements of <60 ml/min per 1.73 m 2 , obtained 90-365 days apart-from the Million Veteran Program and Vanderbilt University Medical Center's DNA biobank. The primary outcome was the annualized relative slope in outpatient eGFR. Analyses were stratified by ethnicity and diabetes status and meta-analyzed thereafter.

RESULTS

In cross-ancestry meta-analysis, the strongest association was rs77924615, near UMOD / PDILT ; each copy of the G allele was associated with a 0.30%/yr faster eGFR decline ( P = 4.9×10 -27 ). We also observed an association within BICC1 (rs11592748), where every additional minor allele was associated with a 0.13%/yr slower eGFR decline ( P = 5.6×10 -9 ). Among participants without diabetes, the strongest association was the UMOD/PDILT variant rs36060036, associated with a 0.27%/yr faster eGFR decline per copy of the C allele ( P = 1.9×10 -17 ). Among Black participants, a significantly faster eGFR decline was associated with variant rs16996674 near APOL1 (R 2 =0.29 with the G1 high-risk genotype); among Black participants with diabetes, lead variant rs11624911 near HEATR4 also was associated with a significantly faster eGFR decline. We also nominally replicated loci with known associations with eGFR decline, near PRKAG2, FGF5, and C15ORF54.

CONCLUSIONS

Three loci were significantly associated with longitudinal eGFR change at genome-wide significance. These findings help characterize molecular mechanisms of eGFR decline and may contribute to the development of new therapeutic approaches for progressive CKD.

Venous thromboembolism risk in cancer patients receiving first-line immune checkpoint inhibitor versus chemotherapy

It remains unclear if immune checkpoint inhibitor (ICI) therapy is associated with higher rate of venous thromboembolism (VTE) compared with cytotoxic chemotherapy (chemo) in patients with comparable cancer type, staging, and comorbidities. Using the national Veterans Affairs healthcare system database from 2016 to 2021, we performed a propensity score (PS)-weighted retrospective cohort study to compare the incidence of VTE in patients with selected stage III/IV cancer receiving first-line ICI versus chemo. The PS model utilized overlap weights to balance age, sex, race, treatment year, VTE history, paralysis/immobilization, prolonged hospitalization, cancer type, staging, time between diagnosis and treatment, and National Cancer Institute comorbidity index. Weighted Cox regressions with robust standard error were used to assess the hazard ratio (HR) and 95% confidence interval (CI). We found that among comparable advanced cancers, first-line ICI (n = 1823) and first-line chemo (n = 6345) had similar rates of VTE (8.49% for ICI and 8.36% for chemo at 6 months). The weighted HR was 1.06 (95% CI 0.88-1.26) for ICI versus chemo. In a subgroup analysis restricted to lung cancers, first-line ICI/chemo (n = 828), ICI monotherapy (n = 428), and chemo monotherapy (n = 4371) had similar rates of VTE (9.60% for ICI/chemo, 10.04% for ICI, and 8.91% for chemo at 6 months). The weighted HR was 1.05 (95% CI 0.77-1.42) for ICI versus chemo, and 1.08 (95% CI 0.83-1.42) for ICI/chemo versus chemo. In conclusion, ICI as a systemic therapy has a similarly elevated risk as cytotoxic chemo for VTE occurrence in cancer patients. This finding can inform future prospective studies exploring thromboprophylaxis strategies.

Dietary ω-3 fatty acids and the incidence of atrial fibrillation in the Million Veteran Program

BACKGROUND

Although recent large randomized clinical trials have reported an increased risk of atrial fibrillation (AF) with marine ω-3 fatty acid supplements, it is unclear whether dietary marine ω-3 fatty acids assessed through food frequency questionnaires are associated with AF risk.

OBJECTIVES

We sought to test the hypothesis that dietary eicosapentaenoic acid/docosahexaenoic acid/docosapentaecnoic acid (EPA/DHA/DPA) is associated with a higher risk of AF in a large prospective cohort of US Veterans.

METHODS

We analyzed data from Million Veteran Program participants who completed self-reported food frequency questionnaires. We used multivariable Cox regression to estimate the HRs of AF across quintiles of ω-3 fatty acid consumption and a cubic spline analysis to assess the dose-response relations between ω-3 fatty acids and AF.

RESULTS

Of the 301,294 veterans studied, the median intake of ω-3 fatty acids (EPA/DHA/DPA) was 219 mg/d (IQR: 144-575), and the mean age was 64.9 y (SD: 12.0); 91% were men, and 84% were White. Consumption of EPA/DHA/DPA exhibited a nonlinear inverse relation with incident AF characterized by an initial decline to 11% at 750 mg/d of marine ω-3 fatty acid intake followed by a plateau.

CONCLUSIONS

Contrary to our hypothesis, dietary EPA/DHA/DPA was not associated with a higher risk of AF but was inversely related to AF risk in a nonlinear manner.

Endogenous sex steroid hormones and risk of liver cancer among US men: Results from the Liver Cancer Pooling Project

Background & Aims

Incidence rates of liver cancer in most populations are two to three times higher among men than women. The higher rates among men have led to the suggestion that androgens are related to increased risk whereas oestrogens are related to decreased risk. This hypothesis was investigated in the present study via a nested case-control analysis of pre-diagnostic sex steroid hormone levels among men in five US cohorts.

Methods

Concentrations of sex steroid hormones and sex hormone-binding globulin were quantitated using gas chromatography-mass spectrometry and a competitive electrochemiluminescence immunoassay, respectively. Multivariable conditional logistic regression was used to calculate odds ratios (ORs) and 95% CIs for associations between hormones and liver cancer among 275 men who subsequently developed liver cancer and 768 comparison men.

Results

Higher concentrations of total testosterone (OR per one-unit increase in log₂ = 1.77, 95% CI = 1.38-2.29), dihydrotestosterone (OR = 1.76, 95% CI = 1.21-2.57), oestrone (OR = 1.74, 95% CI = 1.08-2.79), total oestradiol (OR = 1.58, 95% CI=1.22-20.05), and sex hormone-binding globulin (OR = 1.63, 95% CI = 1.27-2.11) were associated with increased risk. Higher concentrations of dehydroepiandrosterone (DHEA), however, were associated with a 53% decreased risk (OR = 0.47, 95% CI = 0.33-0.68).

Conclusions

Higher concentrations of both androgens (testosterone, dihydrotestosterone) and their aromatised oestrogenic metabolites (oestrone, oestradiol) were observed among men who subsequently developed liver cancer compared with men who did not. As DHEA is an adrenal precursor of both androgens and oestrogens, these results may suggest that a lower capacity to convert DHEA to androgens, and their subsequent conversion to oestrogens, confers a lower risk of liver cancer, whereas a greater capacity to convert DHEA confers a greater risk.

Impact and implications

This study does not fully support the current hormone hypothesis as both androgen and oestrogen levels were associated with increased risk of liver cancer among men. The study also found that higher DHEA levels were associated with lower risk, thus suggesting the hypothesis that greater capacity to convert DHEA could be associated with increased liver cancer risk among men.

Early Outcomes of SARS-CoV-2 Infection in a Multisite Prospective Cohort of Inpatient Veterans

Background

Over 870 000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have occurred among Veterans Health Administration users, and 24 000 have resulted in death. We examined early outcomes of SARS-CoV-2 infection in hospitalized veterans.

Methods

In an ongoing, prospective cohort study, we enrolled veterans age ≥18 tested for SARS-CoV-2 and hospitalized at 15 Department of Veterans Affairs medical centers between February 2021 and June 2022. We estimated adjusted odds ratios (aORs), adjusted incidence rate ratios (aIRRs), and adjusted hazard ratios (aHRs) for maximum illness severity within 30 days of study entry (defined using the 4-category VA Severity Index for coronavirus disease 2019 [COVID-19]), as well as length of hospitalization and rehospitalization within 60 days, in relationship with demographic characteristics, Charlson comorbidity index (CCI), COVID-19 vaccination, and calendar period of enrollment.

Results

The 542 participants included 329 (61%) who completed a primary vaccine series (with or without booster; "vaccinated"), 292 (54%) enrolled as SARS-CoV-2-positive, and 503 (93%) men, with a mean age of 64.4 years. High CCI scores (≥5) occurred in 61 (44%) vaccinated and 29 (19%) unvaccinated SARS-CoV-2-positive participants. Severe illness or death occurred in 29 (21%; 6% died) vaccinated and 31 (20%; 2% died) unvaccinated SARS-CoV-2-positive participants. SARS-CoV-2-positive inpatients per unit increase in CCI had greater multivariable-adjusted odds of severe illness (aOR, 1.21; 95% CI, 1.01-1.45), more hospitalization days (aIRR, 1.06; 95% CI, 1.03-1.10), and rehospitalization (aHR, 1.07; 95% CI, 1.01-1.12).

Conclusions

In a cohort of hospitalized US veterans with SARS-CoV-2 infection, those with a higher CCI had more severe COVID-19 illness, more hospital days, and rehospitalization, after adjusting for vaccination status, age, sex, and calendar period.

Evidence of Novel Susceptibility Variants for Prostate Cancer and a Multiancestry Polygenic Risk Score Associated with Aggressive Disease in Men of African Ancestry

BACKGROUND

Genetic factors play an important role in prostate cancer (PCa) susceptibility.

OBJECTIVE

To discover common genetic variants contributing to the risk of PCa in men of African ancestry.

DESIGN, SETTING, AND PARTICIPANTS

We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness.

RESULTS AND LIMITATIONS

Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4).

CONCLUSIONS

This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry.

PATIENT SUMMARY

In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.

Agent Orange exposure and prostate cancer risk in the Million Veteran Program

Purpose

Exposure to Agent Orange, a known carcinogen, might increase risk of prostate cancer (PCa). We sought to investigate the association of Agent Orange exposure and PCa risk when accounting for race/ethnicity, family history, and genetic risk in a diverse population of US Vietnam War veterans.

Methods & Materials

This study utilized the Million Veteran Program (MVP), a national, population-based cohort study of United States military veterans conducted 2011-2021 with 590,750 male participants available for analysis. Agent Orange exposure was obtained using records from the Department of Veterans Affairs (VA) using the US government definition of Agent Orange exposure: active service in Vietnam while Agent Orange was in use. Only veterans who were on active duty (anywhere in the world) during the Vietnam War were included in this analysis (211,180 participants). Genetic risk was assessed via a previously validated polygenic hazard score calculated from genotype data. Age at diagnosis of any PCa, diagnosis of metastatic PCa, and death from PCa were assessed via Cox proportional hazards models.

Results

Exposure to Agent Orange was associated with increased PCa diagnosis (HR 1.04, 95% CI 1.01-1.06, p=0.003), primarily among Non-Hispanic White men (HR 1.09, 95% CI 1.06- 1.12, p<0.001). When accounting for race/ethnicity and family history, Agent Orange exposure remained an independent risk factor for PCa diagnosis (HR 1.06, 95% CI 1.04-1.09, p<0.05). Univariable associations of Agent Orange exposure with PCa metastasis (HR 1.08, 95% CI 0.99-1.17) and PCa death (HR 1.02, 95% CI 0.84-1.22) did not reach significance on multivariable analysis. Similar results were found when accounting for polygenic hazard score.

Conclusions

Among US Vietnam War veterans, Agent Orange exposure is an independent risk factor for PCa diagnosis, though associations with PCa metastasis or death are unclear when accounting for race/ethnicity, family history, and/or polygenic risk.

Derivation and Validation of a Clinical Risk Assessment Model for Cancer-Associated Thrombosis in Two Unique US Health Care Systems

PURPOSE

Venous thromboembolism (VTE), especially pulmonary embolism (PE) and lower extremity deep vein thrombosis (LE-DVT), is a serious and potentially preventable complication for patients with cancer undergoing systemic therapy.

METHODS

Using retrospective data from patients diagnosed with incident cancer from 2011-2020, we derived a parsimonious risk assessment model (RAM) using least absolute shrinkage and selection operator regression from the Harris Health System (HHS, n = 9,769) and externally validated it using the Veterans Affairs (VA) health care system (n = 79,517). Bootstrapped c statistics and calibration curves were used to assess external model discrimination and fit. Dichotomized risk strata using integer scores were created and compared against the Khorana score (KS).

RESULTS

Incident VTE and PE/LE-DVT at 6 months occurred in 590 (6.2%) and 437 (4.6%) patients in HHS and 4,027 (5.1%) and 3,331 (4.2%) patients in the VA health care system. Assessed at the time of systemic therapy initiation, the new RAM included components of the KS with the modified cancer subtype, cancer staging, systemic therapy class, history of VTE, history of paralysis/immobility, recent hospitalization, and Asian/Pacific Islander race. The c statistic was 0.71 in HHS and 0.68 in the VA health care system (compared with 0.65 and 0.60, respectively, for KS). Furthermore, the new RAM appropriately reclassified 28% of patients and increased the proportion of VTEs in the high-risk group from 37% to 68% in the validation data set.

CONCLUSION

The novel RAM stratified patients with cancer into a high-risk group with 8%-10% cumulative incidence of VTE and 7% PE/LE-DVT at 6 months (v 3% and 2%, respectively, in the low-risk group). The model had improved performance over the original KS and doubled the number of VTE events in the high-risk stratum. We encourage additional external validation from prospective studies.[Media: see text].

Cardiovascular Disease Risk Assessment Using Traditional Risk Factors and Polygenic Risk Scores in the Million Veteran Program

Importance

Primary prevention of atherosclerotic cardiovascular disease (ASCVD) relies on risk stratification. Genome-wide polygenic risk scores (PRSs) are proposed to improve ASCVD risk estimation.

Objective

To 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 Participants

This 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.

Exposures

PRSs 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 Measures

Incident nonfatal myocardial infarction (MI), ischemic stroke, ASCVD death, and composite ASCVD events.

Results

A 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 Relevance

Study 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.

Epidemiology of Cancer-Associated Venous Thromboembolism in Patients With Solid and Hematologic Neoplasms in the Veterans Affairs Health Care System

Importance

Identifying changes in epidemiologic patterns of the incidence and risk of cancer-associated thrombosis (CAT), particularly with evolving cancer-directed therapy, is essential for risk stratification.

Objective

To assess the incidence of CAT over time and to determine pertinent patient-specific, cancer-specific, and treatment-specific factors associated with its risk.

Design, Setting, and Participants

This longitudinal, retrospective cohort study was conducted from 2006 to 2021. Duration of follow-up was from the date of diagnosis until first venous thromboembolism (VTE) event, death, loss of follow-up (defined as a 90-day gap without clinical encounters), or administrative censoring on April 1, 2022. The study took place within the US Department of Veterans Affairs national health care system. Patients with newly diagnosed invasive solid tumors and hematologic neoplasms were included in the study. Data were analyzed from December 2022 to February 2023.

Exposure

Newly diagnosed invasive solid tumors and hematologic neoplasms.

Main Outcomes

Incidence of VTE was assessed using a combination of International Classification of Diseases, Ninth Revision, Clinical Modification and International Statistical Classification of Diseases, Tenth Revision, Clinical Modification and natural language processing confirmed outcomes. Cumulative incidence competing risk functions were used to estimate incidence of CAT. Multivariable Cox regression models were built to assess the association of baseline variables with CAT. Pertinent patient variables included demographics, region, rurality, area deprivation index, National Cancer Institute comorbidity index, cancer type, staging, first-line systemic treatment within 3 months (time-varying covariate), and other factors that could be associated with the risk of VTE.

Results

A total of 434 203 patients (420 244 men [96.8%]; median [IQR] age, 67 [62-74] years; 7414 Asian or Pacific Islander patients [1.7%]; 20 193 Hispanic patients [4.7%]; 89 371 non-Hispanic Black patients [20.6%]; 313 157 non-Hispanic White patients [72.1%]) met the inclusion criteria. Overall incidence of CAT at 12 months was 4.5%, with yearly trends ranging stably from 4.2% to 4.7%. The risk of VTE was associated with cancer type and stage. In addition to confirming well-known risk distribution among patients with solid tumors, a higher risk of VTE was observed among patients with aggressive lymphoid neoplasms compared with patients with indolent lymphoid or myeloid hematologic neoplasms. Compared with no treatment, patients receiving first-line chemotherapy (hazard ratio [HR], 1.44; 95% CI, 1.40-1.49) and immune checkpoint inhibitors (HR, 1.49; 95% CI, 1.22-1.82) had a higher adjusted relative risk than patients receiving targeted therapy (HR, 1.21; 95% CI, 1.13-1.30) or endocrine therapy (HR, 1.20; 95% CI, 1.12-1.28). Finally, adjusted VTE risk was significantly higher among Non-Hispanic Black patients (HR, 1.23; 95% CI, 1.19-1.27) and significantly lower in Asian or Pacific Islander patients (HR, 0.84; 95% CI, 0.76-0.93) compared with Non-Hispanic White patients.

Conclusions and Relevance

In this cohort study of patients with cancer, a high incidence of VTE was observed, with yearly trends that remained stable over the 16-year study period. Both novel and known factors associated with the risk of CAT were identified, providing valuable and applicable insights in this current treatment landscape.

Overview of Efforts to Increase Women Enrollment in the Veterans Affairs Million Veteran Program

Background

Ensuring enhanced delivery of care to women Veterans is a top Veterans Affairs (VA) priority; however, women are historically underrepresented in research that informs evidence-based health care. A primary barrier to women's participation is the inability to engage with research in person due to a number of documented challenges. The VA Million Veteran Program (MVP) is committed to increasing access for women Veterans to participate in research, thereby better understanding conditions specific to this population and how disease manifests differently in women compared to men. The goal of this work is to describe the results of the MVP Women's Campaign, an effort designed to increase outreach to and awareness of remote enrollment options for women Veterans.

Materials and Methods

The MVP Women's Campaign launched two phases between March 2021 and April 2022: the Multimedia Phase leveraged a variety of strategic multichannel communication tactics and the Email Phase focused on direct email communication to women Veterans. The effect of the Multimedia Phase was determined using t-tests and chi-square tests, as well as logistic regression models to compare demographic subgroups. The Email Phase was evaluated using comparisons of the enrollment rate across demographic groups through a multivariate adjusted logistic regression model.

Results

Overall, 4694 women Veterans enrolled during the MVP Women's Campaign (54% during the Multimedia Phase and 46% during the Email Phase). For the Multimedia Phase, the percentage of older women online enrollees increased, along with women from the southwest and western regions of the United States. Differences for women Veteran online enrollment across different ethnicity and race groups were not observed. During the Email Phase, the enrollment rate increased with age. Compared to White women Veterans, Blacks, Asians, and Native Americans were significantly less likely to enroll while Veterans with multiple races were more likely to enroll.

Conclusion

The MVP Women's Campaign is the first large-scale outreach effort focusing on recruitment of women Veterans into MVP. The combination of print and digital outreach tactics and direct email recruitment resulted in over a fivefold increase in women Veteran enrollees during a 7-month period. Attention to messaging and communication channels, combined with a better understanding of effective recruitment methods for certain Veteran populations, allows MVP the opportunity to advance health and health care not only for women Veterans, but beyond. Lessons learned will be applied to increase other populations in MVP such as Blacks, Hispanics, Asians, Native Americans, younger Veterans, and Veterans with certain health conditions.

ARCH: Large-scale Knowledge Graph via Aggregated Narrative Codified Health Records Analysis

Objective

Electronic health record (EHR) systems contain a wealth of clinical data stored as both codified data and free-text narrative notes, covering hundreds of thousands of clinical concepts available for research and clinical care. The complex, massive, heterogeneous, and noisy nature of EHR data imposes significant challenges for feature representation, information extraction, and uncertainty quantification. To address these challenges, we proposed an efficient Aggregated naRrative Codified Health (ARCH) records analysis to generate a large-scale knowledge graph (KG) for a comprehensive set of EHR codified and narrative features.

Methods

The ARCH algorithm first derives embedding vectors from a co-occurrence matrix of all EHR concepts and then generates cosine similarities along with associated p -values to measure the strength of relatedness between clinical features with statistical certainty quantification. In the final step, ARCH performs a sparse embedding regression to remove indirect linkage between entity pairs. We validated the clinical utility of the ARCH knowledge graph, generated from 12.5 million patients in the Veterans Affairs (VA) healthcare system, through downstream tasks including detecting known relationships between entity pairs, predicting drug side effects, disease phenotyping, as well as sub-typing Alzheimer's disease patients.

Results

ARCH produces high-quality clinical embeddings and KG for over 60,000 EHR concepts, as visualized in the R-shiny powered web-API (https://celehs.hms.harvard.edu/ARCH/). The ARCH embeddings attained an average area under the ROC curve (AUC) of 0.926 and 0.861 for detecting pairs of similar EHR concepts when the concepts are mapped to codified data and to NLP data; and 0.810 (codified) and 0.843 (NLP) for detecting related pairs. Based on the p -values computed by ARCH, the sensitivity of detecting similar and related entity pairs are 0.906 and 0.888 under false discovery rate (FDR) control of 5%. For detecting drug side effects, the cosine similarity based on the ARCH semantic representations achieved an AUC of 0.723 while the AUC improved to 0.826 after few-shot training via minimizing the loss function on the training data set. Incorporating NLP data substantially improved the ability to detect side effects in the EHR. For example, based on unsupervised ARCH embeddings, the power of detecting drug-side effects pairs when using codified data only was 0.15, much lower than the power of 0.51 when using both codified and NLP concepts. Compared to existing large-scale representation learning methods including PubmedBERT, BioBERT and SAPBERT, ARCH attains the most robust performance and substantially higher accuracy in detecting these relationships. Incorporating ARCH selected features in weakly supervised phenotyping algorithms can improve the robustness of algorithm performance, especially for diseases that benefit from NLP features as supporting evidence. For example, the phenotyping algorithm for depression attained an AUC of 0.927 when using ARCH selected features but only 0.857 when using codified features selected via the KESER network[1]. In addition, embeddings and knowledge graphs generated from the ARCH network were able to cluster AD patients into two subgroups, where the fast progression subgroup had a much higher mortality rate.

Conclusions

The proposed ARCH algorithm generates large-scale high-quality semantic representations and knowledge graph for both codified and NLP EHR features, useful for a wide range of predictive modeling tasks.

Multi-ancestry GWAS of Fuchs corneal dystrophy highlights roles of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular pathophysiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program (MVP) and meta-analyzed with the previous largest FECD GWAS, finding twelve significant loci (eight novel). We further confirmed the TCF4 locus in admixed African and Hispanic/Latino ancestries, and found an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and co-localization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Validation of algorithms to select patients with multiple myeloma and patients initiating myeloma treatment in the national Veterans Affairs Healthcare System

BACKGROUND

We aimed to evaluate and compare the performance of multiple myeloma (MM) selection algorithms for use in Veterans Affairs (VA) research.

METHODS

Using the VA Corporate Data Warehouse (CDW), the VA Cancer Registry (VACR), and VA pharmacy data, we randomly selected 500 patients from 01/01/1999 to 06/01/2021 who had (1) either one MM diagnostic code OR were listed in the VACR as having MM AND (2) at least one MM treatment code. A team reviewed oncology notes for each veteran to annotate details regarding MM diagnosis and initial treatment within VA. We evaluated inter-annotator agreement and compared the performance of four published algorithms (two developed and validated external to VA data and two used in VA data).

RESULTS

A total of 859 patients were reviewed to obtain 500 patients who were annotated as having MM and initiating MM treatment in VA. Agreement was high among annotators for all variables: MM diagnosis (98.3% agreement, Kappa = 0.93); initial treatment in VA (91.8% agreement; Kappa = 0.77); and initial treatment classification (87.6% agreement; Kappa = 0.86). VA Algorithms were more specific and had higher PPVs than non-VA algorithms for both MM diagnosis and initial treatment in VA. We developed the "VA Recommended Algorithm," which had the highest PPV among all algorithms in identifying patients diagnosed with MM (PPV = 0.98, 95% CI = 0.95-0.99) and in identifying patients who initiated their MM treatment in VA (PPV = 0.93, 95% CI = 0.90-0.96).

CONCLUSION

Our VA Recommended Algorithm optimizes sensitivity and PPV for cohort selection and treatment classification.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder

Objective

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed mental and physical health problems, prior pharmacological treatments, and aggregate genetic factors has potential to inform risk stratification and mitigation strategies.

Methods

In this study of 3,942 SCZ and 5,414 BPI patients receiving VA care, self-reported SB and ideation were assessed using the Columbia Suicide Severity Rating Scale (C-SSRS). These cross-sectional data were integrated with electronic health records (EHR), and compared by lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. Polygenic scores (PGS) for traits related to psychiatric disorders, substance use, and cognition were constructed using available genomic data, and exploratory genome-wide association studies were performed to identify and prioritize specific loci.

Results

Only 20% of veterans who self-reported SB had a corroborating ICD-9/10 code in their EHR; and among those who denied prior behaviors, more than 20% reported new-onset SB at follow-up. SB were associated with a range of psychiatric and non-psychiatric diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking, suicide attempt, and major depressive disorder were also associated with attempt and ideation.

Conclusions

Among individuals with a diagnosed mental illness, a GWAS for SB did not yield any significant loci. Self-reported SB were strongly associated with clinical variables across several EHR domains. Overall, clinical and polygenic analyses point to sequelae of substance-use related behaviors and other psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in clinical settings, underscoring the value of regular screening based on direct, in-person assessments, especially among high-risk individuals.

Framework of the Centralized Interactive Phenomics Resource (CIPHER) standard for electronic health data-based phenomics knowledgebase

The development of phenotypes using electronic health records is a resource-intensive process. Therefore, the cataloging of phenotype algorithm metadata for reuse is critical to accelerate clinical research. The Department of Veterans Affairs (VA) has developed a standard for phenotype metadata collection which is currently used in the VA phenomics knowledgebase library, CIPHER (Centralized Interactive Phenomics Resource), to capture over 5000 phenotypes. The CIPHER standard improves upon existing phenotype library metadata collection by capturing the context of algorithm development, phenotyping method used, and approach to validation. While the standard was iteratively developed with VA phenomics experts, it is applicable to the capture of phenotypes across healthcare systems. We describe the framework of the CIPHER standard for phenotype metadata collection, the rationale for its development, and its current application to the largest healthcare system in the United States.