Biomedical subjects

VA Million Veteran Program

Publications and source records attributed to VA Million Veteran Program.

13 recordsLinked to original sources

PubMed2026

Cross-Ancestry and Phenome-Wide Associations of Cancer-Specific Polygenic Risk Scores.

PURPOSE: Genome-wide association studies have identified many common variants associated at low effect sizes with various cancers. Summing the effects of these variants into polygenic risk scores (PRS) can improve cancer risk prediction. However, cross-cancer and cross-phenotype pleiotropic associations of cancer-specific PRS are limited. METHODS: Using logistic regression models, we tested the association of 13 cancer-specific PRS with curated phenotypes representing the same 13 cancers and 340 cancer and cardiometabolic phecodes in 560,287 individuals (114,255 African ancestry [AFR] and 446,032 European ancestry [EUR]) from the Million Veteran Program. Models were stratified by ancestry and used age, principal components, and cancer-specific PRS per standard deviation as independent variables, and correction was applied for multiple comparisons. RESULTS: All 13 cancer-specific PRS were significantly associated with their respective cancers among EUR individuals with odds ratios per standard deviation of PRS (odds ratio [OR]) 1.05-1.70. Among AFR individuals, the effect sizes of the cancer PRS were lower, with OR 1.01-1.48, and cancer-specific PRS were significantly associated with their respective cancers for five of 13 cancers (bladder, breast in female patients, colorectal, prostate, and thyroid). In cancer-cancer pleiotropy studies, only the renal cancer-specific PRS was significantly associated with skin cancer (OR = 1.04, P = 4.5 × 10-06) among EUR individuals. PheWAS demonstrated five positive pleotropic associations with cardiometabolic conditions (thyroid cancer PRS with thyroid goiter, oral cancer PRS with diabetes phenotypes, and hypothyroidism) and two negative associations (oral and lung cancer PRS separately with coronary artery disease). CONCLUSION: Cancer PRS have stronger associations per cancer among EUR versus AFR individuals. In contrast to PRS of other chronic diseases, the majority of cancer-related PRS are highly specific and pleiotropic associations with other cancers and cardiometabolic traits are uncommon.

Female↗

PubMed2026

Phenotyping strategies for chronic overlapping pain conditions and internalizing disorders in Veterans: Prevalence, comorbidity, and latent structure as evidence for construct validity.

Chronic overlapping pain conditions (COPCs), internalizing (INT) disorders, and opioid use disorder (OUD) are common, comorbid, and difficult to phenotype at scale. Electronic health record (EHR) studies commonly define cases using Any Code (AC; ≥1 ICD-9/10 code) and Multiple Code (MC; ≥1 inpatient or ≥2 outpatient codes) phenotyping strategies, but it is unclear whether these thresholds change only case numbers or also the clinical relationships among conditions. This cross-sectional study included approximately 950,000 Million Veteran Program participants with ≥2 visits. AC and MC phenotypes for 17 conditions spanning COPCs, INT, and OUD were compared in prevalence, case characteristics, comorbidity, and latent structure. Random-thinning analysis compared AC-MC differences to case reduction alone. Construct validity was evaluated through correspondence with expected patterns of association and latent organization. Back pain (AC=58.1%; MC=48.1%), major depressive disorder (41.5%; 36.2%), and post-traumatic stress disorder (32.6%; 29.1%) were most prevalent. MC excluded 33.5% of AC cases on average, and MC cases had greater healthcare utilization, diagnostic burden, opioid exposure, and psychiatric medication use than AC-only cases. The observed mean absolute correlation change (mean |Δr|=0.014) was smaller than in all 1000 random-thinning replicates. Both strategies supported a correlated, four-factor model consistent with "Anxious Misery," "Fear," "Diffuse Pain," and "Head Pain" (AC: CFI=0.987, RMSEA=0.015; MC: CFI=0.987, RMSEA=0.014). The MC strategy reduced prevalence and altered case composition but maintained the expected comorbidity and latent organization patterns among conditions. Findings provide evidence of phenotype construct validity and inform selection of EHR phenotyping strategies for epidemiological and genomic research. PERSPECTIVE: Commonly used EHR phenotyping strategies tested in nearly one million Veterans produce broadly similar latent organization across comorbid and prevalent chronic overlapping pain conditions, internalizing disorders, and opioid use disorder. Findings support construct validity and clarify trade-offs involving case inclusion, recorded burden, and healthcare observation in large-scale research.

Chronic overlapping pain conditions↗

PubMed2026

Multi-ancestry genetic architecture of heart failure subtypes.

Heart failure (HF) affects 6.7 million people in the US and includes two major subtypes, HF with reduced ejection fraction (HFrEF) and HF with preserved ejection fraction (HFpEF), with distinct genetic architectures. We meta-analyze genome-wide association studies (GWAS) of 38,781 HFrEF cases, 38,163 HFpEF cases, and 526,135 controls across European, African, Hispanic, and Asian ancestries using the Million Veteran Program and Vanderbilt University DNA Databank (BioVU). We identify 46 genome-wide significant loci for HFrEF (9 novel) and 3 loci for HFpEF (1 novel). Four HFrEF loci are detected in African ancestry participants near CD36, SPI1, TRIM48, and SPNS3, with lead SNPs showing low risk-allele frequencies in European populations. In the all-cause HF meta-analysis (200,070 cases, 2,076,466 controls), we identify 136 loci (12 novel). Gene-based tests, tissue enrichment, transcriptome-wide association, and fine-mapping implicate vascular, metabolic, and TGF-β/Smad signaling pathways and nominate candidate causal genes, clarifying shared and subtype-specific risk across ancestries.

Humans↗

PubMed2026

A Common CD36 Variant and the Genetic Landscape of Dilated Cardiomyopathy in Individuals of African Ancestry.

IMPORTANCE: Dilated cardiomyopathy (DCM) is a major cause of heart failure that disproportionately affects individuals of African genetic ancestry (AFR), among whom familial clustering of disease is also more pronounced relative to those of European ancestry (EUR). However, established monogenic DCM genes, identified primarily in EUR populations, explain a smaller proportion of DCM cases in AFR populations. A recent study identified a common AFR-specific nonsense variant in CD36 that accounts for a substantial burden of DCM in AFR. How the risk and population impact of this variant compare with those of established genetic causes of DCM is unknown. OBJECTIVE: To compare the contribution of a CD36 nonsense variant to DCM risk with that of truncating variants in TTN and pathogenic or likely pathogenic (P/LP) variants in other established DCM genes. DESIGN SETTING AND PARTICIPANTS: Multicohort genetic association study including AFR and EUR participants with exome or genome sequence and DCM case status from four datasets: All of Us, Million Veteran Program, Penn Medicine Biobank, and the DCM Precision Medicine Study. EXPOSURE: Carrier status for TTN truncating variants, P/LP variants in 11 high confidence DCM genes, and the CD36 nonsense variant (Y325*; 0, 1, or 2 copies). MAIN OUTCOMES AND MEASURES: Odds of DCM; prevalence of risk-variant carriers among DCM cases; and population attributable fraction (PAF) for DCM. RESULTS: Among 82,623 AFR individuals across four studies, the mean age was 53.4 years and 1,625 had DCM. CD36 Y325* risk-allele homozygotes had 4.8-fold (95% CI, 3.1-7.3) increased odds of DCM, and CD36 Y325* heterozygotes had 1.4-fold (95% CI, 1.2-1.7) increased odds. TTN truncating variants also conferred elevated risk of DCM in AFR participants (OR, 8.46; 95% CI, 5.3-12.3). Among AFR DCM cases, 2.5% were CD36 homozygotes, second only to TTN truncating variants (4.3%) and exceeding all other high-confidence DCM genes combined (1.5%). In population-level analyses incorporating both heterozygous and homozygous CD36 Y325* carriers, the population-attributable fraction for CD36 (9.0%) surpassed that of TTN truncating variants (3.6%). CONCLUSIONS AND RELEVANCE: An ancestry-specific CD36 variant contributes more to DCM burden in AFR ancestry than established DCM genes, including TTN truncating variants, typically considered the most common genetic cause of DCM. These findings reshape the known genetic architecture of DCM in individuals of African ancestry and highlight the importance of representation in genomic research.

Journal Article↗

PubMed2026

SAIGE-GPU: accelerating genome- and phenome-wide association studies using GPUs.

MOTIVATION: Genome-wide association studies (GWAS) at biobank scale are computationally intensive, especially for admixed populations requiring robust statistical models. SAIGE is a widely used method for generalized linear mixed-model GWAS but is limited by its CPU-based implementation, making phenome-wide association studies impractical for many research groups. RESULTS: We developed SAIGE-GPU, a GPU-accelerated version of SAIGE that replaces CPU-intensive matrix operations with GPU-optimized kernels. The core innovation is distributing genetic relationship matrix calculations across GPUs and communication layers. Applied to 2068 phenotypes from 635 969 participants in the Million Veteran Program, including diverse and admixed populations, SAIGE-GPU achieved a 5-fold speedup in mixed model fitting on supercomputing infrastructure and cloud platforms. We further optimized the variant association testing step through multi-core and multi-trait parallelization. Deployed on Google Cloud Platform and Azure, the method provided substantial cost and time savings. AVAILABILITY AND IMPLEMENTATION: Source code and binaries are available for download at https://github.com/saigegit/SAIGE/tree/SAIGE-GPU-1.3.3. A code snapshot is archived at Zenodo for reproducibility (DOI: [10.5281/zenodo.17642591]). SAIGE-GPU is available in a containerized format for use across HPC and cloud environments and is implemented in R/C++ and runs on Linux systems.

Genome-Wide Association Study↗

PubMed2025

A multi-ancestry polygenic risk score for Alzheimer disease is associated with cognitive decline, hippocampal atrophy and neuropathological hallmarks in diverse populations.

Alzheimer disease (AD) has a strong genetic basis, yet previously derived polygenic risk scores (PRS) are heavily weighted by the APOE locus and perform inconsistently across diverse ancestries. We developed an APOE-independent multi-ancestry AD PRS using genome-wide association study summary statistics from cohorts in the United States, Europe and East Asia that were applied to European ancestry (EA), African American (AA), Caribbean Hispanic (CH), and East Asian cohorts from the Alzheimer's Disease Genetics Consortium. PRS performance was evaluated in the multi-ancestry Alzheimer's Disease Sequencing Project (ADSP) dataset and validated in several additional multi-ancestry cohorts. The PRS was significantly associated with AD in the ADSP EA, AA, CH, and Native American Hispanic groups with adjusted odds ratios (ORs) between 1.14 and 1.52 per standard deviation of the PRS. PRS performance was validated in the replication cohorts (ORs 1.21-1.65). The PRS was also associated with poorer memory, executive function, and language performance; greater AD-related neuropathological burden (including CERAD, Braak stage, and Thal phase scores); reduced hippocampal volume; lower CSF Aβ42; and elevated total tau and phosphorylated tau (p-tau), with stronger p-tau associations observed in women. Longitudinal analyses revealed that individuals in the highest PRS decile exhibited the steepest cognitive decline, particularly among those who progressed to AD. Our findings demonstrate the utility of an ancestry-aware and APOE-independent PRS for advancing understanding of the genetic basis of AD across diverse populations. Associations observed with early biological and cognitive changes and potential sex-specific differences support the incorporation of a PRS in clinical trials and personalized intervention and prevention strategies.

Journal Article↗

PubMed2025

Germline Variants Influence Chronic Liver Disease Progression through Distinct Pathways.

Cirrhosis and hepatocellular carcinoma (HCC) are long-term complications of chronic liver disease (CLD). In this large multi-ancestry genome-wide association study of all-cause cirrhosis (35,481 cases, 2.36M controls) and HCC (6,680 cases, 1.76M controls), we identified 27 loci associated with cirrhosis (10 novel) and 11 with HCC (three novel). Three novel cirrhosis loci were replicated in independent cohorts (e.g. FGF21, RPTOR, and IFNL3/4). Fifteen cirrhosis loci exhibited differential effects on cirrhosis risk via underlying etiologies, and six HCC loci influenced HCC risk indirectly via cirrhosis. In a gene-burden analysis of rare variants from whole-genome sequencing data in the VA Million Veteran Program (n=102,677), we identified GSTA5 as a novel cirrhosis-associated gene, while APOB and ATP9B were associated with and replicated for HCC. A high genetic risk score for cirrhosis was associated with a nearly doubled risk of CLD progressing to cirrhosis (HR=1.94, P=2×10-68) and of cirrhosis progressing to HCC (HR=1.65, P=7×10-08). Finally, among individuals with chronic hepatitis C who underwent antiviral therapy, cirrhosis risk was modified by variants in PNPLA3, IFNL3/4, and CD81 following pegylated interferon-α therapy, and by APOE lead variant following direct-acting antiviral therapy. These findings provide new insights into the complex genetic architecture of CLD progression with potential clinical and therapeutic implications.

Journal Article↗

PubMed2025

Genomics-informed drug-repurposing strategy identifies two therapeutic targets for preventing liver disease associated with metabolic dysfunction.

Identification of drug-repurposing targets with genetic and biological support is an economically and temporally efficient strategy for improving the treatment of diseases. We employed a cross-disciplinary approach to identify potential therapeutics for the prevention of metabolic-dysfunction-associated steatotic liver disease (MASLD) in at-risk individuals by using humans as a model organism. We identified 212 putative candidate genes associated with MASLD by using data from a large multi-ancestry genetic association study, of which 158 (74.5%) were previously unreported. From this set, we identified 57 genes that encode for druggable protein targets and for which the effects of increasing genetically predicted gene expression on MASLD risk align with the function of that drug on the protein target. We then used We then evaluated these potential targets for evidence of efficacy by using Mendelian randomization, pathway analysis, and protein structural modeling. Through these approaches, we present compelling evidence to suggest that the activation of FADS1 by icosapent ethyl, as well as S1PR2 by fingolimod, could be a promising therapeutic strategy for MASLD prevention.

Humans↗

PubMed2025

Polygenic scores for obstructive sleep apnoea reveal pathways contributing to cardiovascular disease.

BACKGROUND: Obstructive sleep apnoea (OSA) is a common chronic condition, with obesity its strongest risk factor. Polygenic scores (PGSs) summarise the genetic liability to phenotype and can provide insights into relationships between phenotypes. Recently, large datasets that include genetic data and OSA status became available, providing an opportunity to utilise PGS approaches to study the genetic relationship between OSA and other phenotypes, while differentiating OSA-specific from obesity-specific genetic factors. METHODS: Using race/ethnic diverse samples from over 1.2 million individuals from the Million Veteran Program, FinnGen, TOPMed, All of Us (AoU), Geisinger's MyCode, MGB Biobank, and the Human Phenotype Project, we developed and assessed PGSs for OSA, both without (BMIunadjOSA-PGS) and with adjustment for the genetic contributions of BMI (BMIadjOSA-PGS). FINDINGS: Adjusted odds ratios (ORs) for OSA per 1 standard deviation of the PGSs ranged from 1.38 to 2.75. The associations of BMIadjOSA- and BMIunadjOSA-PGSs with CVD outcomes in AoU shared both common and distinct patterns. Only BMIunadjOSA-PGS was associated with type 2 diabetes, heart failure, and coronary artery disease, while both BMIadjOSA- and BMIunadjOSA-PGSs were associated with hypertension and stroke. Sex stratified analyses revealed that BMIadjOSA-PGS association with hypertension was driven by females (OR = 1.1, p-value = 0.002, OR = 1.01 p-value = 0.2 in males). OSA PGSs were also associated with body fat measures with some sex-specific associations. INTERPRETATION: Distinct components of OSA genetic risk are related and independent of obesity. Sex-specific associations with body fat distribution measures may explain differing OSA risks and associations with cardiometabolic morbidities between sexes. FUNDING: R01AG080598.

Humans↗

PubMed2025

Genome-wide association study of long COVID.

Infections can lead to persistent symptoms and diseases such as shingles after varicella zoster or rheumatic fever after streptococcal infections. Similarly, severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) infection can result in long coronavirus disease (COVID), typically manifesting as fatigue, pulmonary symptoms and cognitive dysfunction. The biological mechanisms behind long COVID remain unclear. We performed a genome-wide association study for long COVID including up to 6,450 long COVID cases and 1,093,995 population controls from 24 studies across 16 countries. We discovered an association of FOXP4 with long COVID, independent of its previously identified association with severe COVID-19. The signal was replicated in 9,500 long COVID cases and 798,835 population controls. Given the transcription factor FOXP4's role in lung physiology and pathology, our findings highlight the importance of lung function in the pathophysiology of long COVID.

Humans↗

PubMed2025

Genetic study of von Willebrand factor antigen levels ≤ 50 IU/dL identifies variants associated with increased risk of von Willebrand disease and bleeding.

BACKGROUND: von Willebrand disease (VWD) is a common inherited bleeding disorder caused by low levels or activity of circulating von Willebrand factor (VWF). Genetic susceptibility to VWF antigen (VWF:Ag) below normal (≤ 50 IU/dL) in the general population is underexplored. OBJECTIVES: To identify genetic variants influencing VWF:Ag levels ≤ 50 IU/dL. METHODS: We performed a genome-wide association study in 926 cases with VWF:Ag levels ≤ 50 IU/dL and 12 846 controls from 7 studies from the Trans-Omics for Precision Medicine program. We then examined whether significant genome-wide findings were also associated with clinical diagnosis of VWD in 5 biobanks with 708 VWD cases and 1 286 069 controls, and with 6 bleeding and thrombotic disorders in FinnGen. RESULTS: Variants at 2 loci were associated (P < 5 × 10-9) with VWF:Ag levels ≤ 50 IU/dL: ABO and VWF. The VWF index variant, p.Tyr1584Cys, is a rare (0.22%) missense variant with odds ratio (OR) of 78.58, while the ABO index variant is a common intronic variant with a smaller effect (OR = 2.52). Notably, both VWF (OR = 7.16) and ABO (OR = 1.57) variants were also associated (P < .025) with diagnosed VWD. Among p.Tyr1584Cys heterozygotes, the penetrance of VWF:Ag levels ≤ 50 IU/dL was 24.2% and the penetrance of diagnosed VWD was 0.3%. p.Tyr1584Cys was associated (P < .0042) with increased odds of heavy menstrual bleeding (OR = 1.27), iron deficiency anemia (OR = 1.55), and intrapartum hemorrhage (OR = 2.20), but decreased odds of deep vein thrombosis (OR = 0.54). CONCLUSIONS: Although there are currently conflicting interpretations of pathogenicity p.Tyr1584Cys, our results suggest that it is a low penetrance pathogenic variant that contributes to VWF:Ag levels ≤ 50 IU/dL, bleeding, and VWD.

Humans↗

PubMed2025

Primary care providers' perspectives on receiving opportunistic genomic results from a national study: The Million Veteran Program Return Of Actionable Results (MVP-ROAR) Study.

PURPOSE: Patients are increasingly obtaining genetic health information and integrating it into their care with the help of their primary care provider (PCP). However, PCPs may not be adequately prepared to effectively utilize genetic results. Across the Veterans Health Administration health system, the Million Veteran Program Return Of Actionable Results-Familial Hypercholesterolemia (MVP-ROAR-FH) Study clinically confirms and returns genetic results associated with familial hypercholesterolemia (FH), identified in a national biobank program. METHODS: PCPs who received their patient's genetic results through the MVP-ROAR-FH study were invited to participate in semistructured interviews, which explored PCPs' familiarity with FH, how the results affected medical management, and suggestions for process improvement. Interviews were transcribed and analyzed using directed content analysis and constant comparison methods to identify key themes. RESULTS: Interviews with 9 PCPs revealed varied levels of familiarity with genetic testing and FH. Most PCPs did not distinguish FH from common high cholesterol issues and already used similar treatment approaches. Many PCPs did not recall receiving results from the MVP-ROAR-FH study. Alerts in medical records were deemed effective for communicating results. PCPs valued genetics in informing patient care and identifying at-risk family members but noted several implementation barriers, such as additional workload and unclear medical management benefits. Recommendations for improving results disclosure included simplifying the genetic testing report and associated support documents. CONCLUSION: The study represents the first investigation into PCPs' experiences with receiving genetic test results from a biobank linked to a national healthcare system. Results suggest that PCPs generally view genetic testing as beneficial, although they may not significantly alter medical management. PCPs expressed that integrating genetics into routine care may be burdensome and require additional training, which may not be practical. The study underscores the need for accessible genetic information, which could be aided by specialized support roles or different clinical specialties assisting with incorporating genetic results into patient care.

Humans↗

PubMed2024

Genome-wide association study of the common retinal disorder epiretinal membrane: Significant risk loci in each of three American populations.

Epiretinal membrane (ERM) is a common retinal condition characterized by the presence of fibrocellular tissue on the retinal surface, often with visual distortion and loss of visual acuity. We studied European American (EUR), African American (AFR), and Latino (admixed American, AMR) ERM participants in the Million Veteran Program (MVP) for genome-wide association analysis-a total of 38,232 case individuals and 557,988 control individuals. We completed a genome-wide association study (GWAS) in each population separately, and then results were meta-analyzed. Genome-wide significant (GWS) associations were observed in all three populations studied: 31 risk loci in EUR subjects, 3 in AFR, and 2 in AMR, with 48 in trans-ancestry meta-analysis. Many results replicated in the FinnGen sample. Several GWS variants associate to alterations in gene expression in the macula. ERM showed significant genetic correlation to multiple traits. Pathway enrichment analyses implicated collagen and collagen-adjacent mechanisms, among others. This well-powered ERM GWAS identified novel genetic associations that point to biological mechanisms for ERM.

Humans↗