Search

To effectively reduce vision loss due to age-related macular generation (AMD) on a global scale, knowledge of its genetic architecture in diverse populations is necessary. A critical element, AMD risk profiles in African and Hispanic/Latino ancestries, remains largely unknown. We combined data in the Million Veteran Program with five other cohorts to conduct the first multi-ancestry genome-wide association study of AMD and discovered 63 loci (30 novel). We observe marked cross-ancestry heterogeneity at major risk loci, especially in African-ancestry populations which demonstrate a primary signal in a major histocompatibility complex class II haplotype and reduced risk at the established CFH and ARMS2/HTRA1 loci. Dissecting local ancestry in admixed individuals, we find significantly smaller marginal effect sizes for CFH risk alleles in African ancestry haplotypes. Broadening efforts to include ancestrally distinct populations helped uncover genes and pathways that boost risk in an ancestry-dependent manner and are potential targets for corrective therapies., Competing Interests: Competing interests: A.S. and B.L.Y. are employees of Genentech/Roche and hold stock and stock options in Roche. E.J. is an employee and a stockholder of Regeneron Genetics Center. The other authors declare no competing interests., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Lung cancer remains the leading cause of cancer mortality, despite declining smoking rates. Previous lung cancer GWAS have identified numerous loci, but separating the genetic risks of lung cancer and smoking behavioral susceptibility remains challenging. Here, we perform multi-ancestry GWAS meta-analyses of lung cancer using the Million Veteran Program cohort (approximately 95% male cases) and a previous study of European-ancestry individuals, jointly comprising 42,102 cases and 181,270 controls, followed by replication in an independent cohort of 19,404 cases and 17,378 controls. We then carry out conditional meta-analyses on cigarettes per day and identify two novel, replicated loci, including the 19p13.11 pleiotropic cancer locus in squamous cell lung carcinoma. Overall, we report twelve novel risk loci for overall lung cancer, lung adenocarcinoma, and squamous cell lung carcinoma, nine of which are externally replicated. Finally, we perform PheWAS on polygenic risk scores for lung cancer, with and without conditioning on smoking. The unconditioned lung cancer polygenic risk score is associated with smoking status in controls, illustrating a reduced predictive utility in non-smokers. Additionally, our polygenic risk score demonstrates smoking-independent pleiotropy of lung cancer risk across neoplasms and metabolic traits., (© 2024. The Author(s).)

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 = .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., Clinical Trial Registration: ClinicalTrials.gov ID NCT04178122., (Published by Elsevier Inc.)

Large-scale genome-wide association studies of schizophrenia have uncovered hundreds of associated loci but with extremely limited representation of African diaspora populations. We surveyed electronic health records of 200,000 individuals of African ancestry in the Million Veteran and All of Us Research Programs, and, coupled with genotype-level data from four case-control studies, realized a combined sample size of 13,012 affected and 54,266 unaffected persons. Three genome-wide significant signals - near PLXNA4 , PMAIP1 , and TRPA1 - are the first to be independently identified in populations of predominantly African ancestry. Joint analyses of African, European, and East Asian ancestries across 86,981 cases and 303,771 controls, yielded 376 distinct autosomal loci, which were refined to 708 putatively causal variants via multi-ancestry fine-mapping. Utilizing single-cell functional genomic data from human brain tissue and two complementary approaches, transcriptome-wide association studies and enhancer-promoter contact mapping, we identified a consensus set of 94 genes across ancestries and pinpointed the specific cell types in which they act. We identified reproducible associations of schizophrenia polygenic risk scores with schizophrenia diagnoses and a range of other mental and physical health problems. Our study addresses a longstanding gap in the generalizability of research findings for schizophrenia across ancestral populations, underlining shared biological underpinnings of schizophrenia across global populations in the presence of broadly divergent risk allele frequencies.

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., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

One of the justifiable criticisms of human genetic studies is the underrepresentation of participants from diverse populations. Lack of inclusion must be addressed at-scale to identify causal disease factors and understand the genetic causes of health disparities. We present genome-wide associations for 2068 traits from 635,969 participants in the Department of Veterans Affairs Million Veteran Program, a longitudinal study of diverse United States Veterans. Systematic analysis revealed 13,672 genomic risk loci; 1608 were only significant after including non-European populations. Fine-mapping identified causal variants at 6318 signals across 613 traits. One-third ( n = 2069) were identified in participants from non-European populations. This reveals a broadly similar genetic architecture across populations, highlights genetic insights gained from underrepresented groups, and presents an extensive atlas of genetic associations.

To efficiently distribute blood flow to cardiac muscle, the coronary artery tree must follow a specific branching pattern over the heart. How this pattern arises in humans is unknown due to the limitations of studying human heart development. Here, we leveraged a natural variation of coronary artery anatomy, known as coronary dominance, in genetic association studies to identify the first known driver of human coronary developmental patterning. Coronary dominance refers to whether the right, left, or both coronary arteries branch over the posterior left ventricle, but whether this variability is heritable and how it would be genetically regulated was completely unknown. By conducting the first large-scale, multi-ancestry genome-wide association study (GWAS) of coronary dominance in 61,043 participants of the VA Million Veteran Program, we observed moderate heritability (27.7%) with ten loci reaching genome wide significance. An exceptionally strong association mapped DNA variants to a non-coding region near the chemokine CXCL12 in both European and African ancestries, which overlapped with variants associated with coronary artery disease. Genomic analyses predicted these variants to impact CXCL12 levels, and imaging revealed dominance to develop during fetal life coincident with CXCL12 expression. Reducing Cxcl12 in mice to model the human genetics altered septal artery dominance patterns and caused coronary branches to develop away from Cxcl12 expression domains. Cxcl12 heterozygosity did not compromise overall artery coverage as seen with full deletion, but instead changed artery patterning, reminiscent of the human scenario. Together, our data support CXCL12 as a critical determinant of human coronary artery growth and patterning and lay a foundation for the utilization of developmental pathways to guide future precision 'medical revascularization' therapeutics., Competing Interests: Competing interests: A.K. is on the scientific advisory board of SerImmune, TensorBio and OpenTargets, and a consultant with Arcadia Science and Inari Agriculture. A.K. was a scientific co-founder of RavelBio, a paid consultant with Illumina, was on the SAB of PatchBio and owns shares in DeepGenomics, Immunai, Freenome, and Illumina. All other authors declare they have no competing interests.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals 1,2 , but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Central serous chorioretinopathy (CSC) is a fluid maculopathy whose etiology is not well understood. Abnormal choroidal veins in CSC patients have been shown to have similarities with varicose veins. To identify potential mechanisms, we analyzed genotype data from 1,477 CSC patients and 455,449 controls in FinnGen. We identified an association for a low-frequency (AF=0.5%) missense variant (rs113791087) in the gene encoding vascular endothelial protein tyrosine phosphatase (VE-PTP) (OR=2.85, P=4.5×10 -9 ). This was confirmed in a meta-analysis of 2,452 CSC patients and 865,767 controls from 4 studies (OR=3.06, P=7.4×10 -15 ). Rs113791087 was associated with a 56% higher prevalence of retinal abnormalities (35.3% vs 22.6%, P=8.0×10 -4 ) in 708 UK Biobank participants and, surprisingly, with varicose veins (OR=1.31, P=2.3×10 -11 ) and glaucoma (OR=0.82, P=6.9×10 -9 ). Predicted loss-of-function variants in VEPTP, though rare in number, were associated with CSC in All of Us (OR=17.10, P=0.018). These findings highlight the significance of VE-PTP in diverse ocular and systemic vascular diseases., Competing Interests: Declaration of Competing Interests Dr. Rossin and Dr. Rämö are named inventors on a provisional patent application that describes the secondary use of intravitreal Anti-Ang2 medications for use in the treatment of central serous chorioretinopathy. Dr. Ellinor receives sponsored research support from Bayer AG, IBM Research, Bristol Myers Squibb, Pfizer and Novo Nordisk; he has also served on advisory boards or consulted for MyoKardia and Bayer AG. The remaining authors declare no competing interests.

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

Objective: Determine the incidence of vestibular disorders in patients with SARS-CoV-2 compared to the control population., Study Design: Retrospective., Setting: Clinical data in the National COVID Cohort Collaborative database (N3C)., Methods: Deidentified patient data from the National COVID Cohort Collaborative database (N3C) were queried based on variant peak prevalence (untyped, alpha, delta, omicron 21K, and omicron 23A) from covariants.org to retrospectively analyze the incidence of vestibular disorders in patients with SARS-CoV-2 compared to control population, consisting of patients without documented evidence of COVID infection during the same period., Results: Patients testing positive for COVID-19 were significantly more likely to have a vestibular disorder compared to the control population. Compared to control patients, the odds ratio of vestibular disorders was significantly elevated in patients with untyped (odds ratio [OR], 2.39; confidence intervals [CI], 2.29-2.50; P < 0.001), alpha (OR, 3.63; CI, 3.48-3.78; P < 0.001), delta (OR, 3.03; CI, 2.94-3.12; P < 0.001), omicron 21K variant (OR, 2.97; CI, 2.90-3.04; P < 0.001), and omicron 23A variant (OR, 8.80; CI, 8.35-9.27; P < 0.001)., Conclusions: The incidence of vestibular disorders differed between COVID-19 variants and was significantly elevated in COVID-19-positive patients compared to the control population. These findings have implications for patient counseling and further research is needed to discern the long-term effects of these findings., Competing Interests: None declared., (Copyright © 2024 The Authors. Published by Wolters Kluwer Health, Inc. on behalf of Otology & Neurotology, Inc.)

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., Competing Interests: J.G. reports a relationship with Novartis that includes grant funding. All other authors report no competing interests., (Published by Oxford University Press on behalf of the American Medical Informatics Association 2024.)

Importance: Non-Hispanic Black (hereafter, Black) individuals experience worse prostate cancer outcomes due to socioeconomic and racial inequities of access to care. Few studies have empirically evaluated these disparities across different health care systems., Objective: To describe the racial and ethnic and neighborhood socioeconomic status (nSES) disparities among residents of the same communities who receive prostate cancer care in the US Department of Veterans Affairs (VA) health care system vs other settings., Design, Setting, and Participants: This cohort study obtained data from the VA Central Cancer Registry for veterans with prostate cancer who received care within the VA Greater Los Angeles Healthcare System (VA cohort) and from the California Cancer Registry (CCR) for nonveterans who received care outside the VA setting (CCR cohort). The cohorts consisted of all males with incident prostate cancer who were living within the same US Census tracts. These individuals received care between 2000 and 2018 and were followed up until death from any cause or censoring on December 31, 2018. Data analyses were conducted between September 2022 and December 2023., Exposures: Health care setting, self-identified race and ethnicity (SIRE), and nSES., Main Outcomes and Measures: The primary outcome was all-cause mortality (ACM). Cox proportional hazards regression models were used to estimate hazard ratios for associations of SIRE and nSES with prostate cancer outcomes in the VA and CCR cohorts., Results: Included in the analysis were 49 461 males with prostate cancer. Of these, 1881 males were in the VA cohort (mean [SD] age, 65.3 [7.7] years; 833 Black individuals [44.3%], 694 non-Hispanic White [hereafter, White] individuals [36.9%], and 354 individuals [18.8%] of other or unknown race). A total of 47 580 individuals were in the CCR cohort (mean [SD] age, 67.0 [9.6] years; 8183 Black individuals [17.2%], 26 206 White individuals [55.1%], and 13 191 individuals [27.8%] of other or unknown race). In the VA cohort, there were no racial disparities observed for metastasis, ACM, or prostate cancer-specific mortality (PCSM). However, in the CCR cohort, the racial disparities were observed for metastasis (adjusted odds ratio [AOR], 1.36; 95% CI, 1.22-1.52), ACM (adjusted hazard ratio [AHR], 1.13; 95% CI, 1.04-1.24), and PCSM (AHR, 1.15; 95% CI, 1.05-1.25). Heterogeneity was observed for the racial disparity in ACM in the VA vs CCR cohorts (AHR, 0.90 [95% CI, 0.76-1.06] vs 1.13 [95% CI, 1.04-1.24]; P = .01). No evidence of nSES disparities was observed for any prostate cancer outcomes in the VA cohort. However, in the CCR cohort, heterogeneity was observed for nSES disparities with ACM (AHR, 0.82; 95% CI, 0.80-0.84; P = .002) and PCSM (AHR, 0.86; 95% CI, 0.82-0.89; P = .007)., Conclusions and Relevance: Results of this study suggest that racial and nSES disparities were wider among patients seeking care outside of the VA health care system. Health systems-related interventions that address access barriers may mitigate racial and socioeconomic disparities in prostate cancer.

Purpose: Chronic age-related imbalance is a common cause of falls and subsequent death in the elderly and can arise from dysfunction of the vestibular system, an elegant neuroanatomical group of pathways that mediates human perception of acceleration, gravity, and angular head motion. Studies indicate that 27-46% of the risk of age-related chronic imbalance is genetic; nevertheless, the underlying genes remain unknown., Methods: The cohort consisted of 50,339 cases and 366,900 controls in the Million Veteran Program. The phenotype comprised cases with two ICD diagnoses of vertigo or dizziness at least 6 months apart, excluding acute or recurrent vertiginous syndromes and other non-vestibular disorders. Genome-wide association studies were performed as individual logistic regressions on European, African American, and Hispanic ancestries followed by trans-ancestry meta-analysis. Downstream analysis included case-case-GWAS, fine mapping, probabilistic colocalization of significant variants and genes with eQTLs, and functional analysis of significant hits., Results: Two significant loci were identified in Europeans, another in the Hispanic population, and two additional in trans-ancestry meta-analysis, including three novel loci. Fine mapping revealed credible sets of intronic single nucleotide polymorphisms (SNPs) in MLLT10 - a histone methyl transferase cofactor, BPTF - a subunit of a nucleosome remodeling complex implicated in neurodevelopment, and LINC01224 - a proto-oncogene receptor tyrosine kinase., Conclusion: Despite the difficulties of phenotyping the nature of chronic imbalance, we replicated two loci from previous vertigo GWAS studies and identified three novel loci. Findings suggest candidates for further study and ultimate treatment of this common elderly disorder., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

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., Competing Interests: Author disclosures: The authors report no actual or potential conflicts of interest or outside sources of funding with regard to this article., (Copyright © 2023 Frontline Medical Communications Inc., Parsippany, NJ, USA.)

Background: To address the COVID-19 pandemic and future threats, VA leadership assembled research and clinical teams to coordinate a unified response, which included creating the VA Science and Health Initiative to Combat Infectious and Emerging Life-Threatening Diseases (VA SHIELD)., Observations: VA SHIELD is a comprehensive specimen and data repository. It links specific types of biospecimens with data regarding genetics, exposure, and disease risk by connecting data sources and the collections of biospecimens across clinical and research environments. Researchers can test novel diagnostic platforms and therapeutics for new and existing diseases, allowing for an expedited, more robust, and informed response. The existing longitudinal disease risk-factor information, records of causal processes, and outcomes data present an unparalleled opportunity to optimize prevention, diagnosis, and treatment of many acute and chronic diseases., Conclusions: VA SHIELD will expand to become an enterprise resource for investigators and public health officials. The alignment of basic science, clinical, and translational research goals under one governance is a significant advancement. VA SHIELD has the opportunity to transform the VA research enterprise by creating an entirely new biorepository., Competing Interests: Author disclosures: The authors report no actual or potential conflicts of interest or outside sources of funding with regard to this article., (Copyright © 2023 Frontline Medical Communications Inc., Parsippany, NJ, USA.)

Background: The US Department of Veterans Affairs (VA) has dedicated significant resources toward countering the COVID-19 pandemic. Sequencing for Research Clinical and Epidemiology (SeqFORCE) and Sequencing Collaborations United for Research and Epidemiology (SeqCURE) were developed as clinical and research consortiums, respectively, focused on the genetic COVID-19 surveillance., Observations: Through genetic sequencing, VA SeqFORCE and SeqCURE collaborations contributed to the COVID-19 pandemic response and scientific understanding. Future directions for each program include the assessment of the unique impact of COVID-19 on the veteran population, as well as the adaptation of these programs to future infectious disease threats. We foresee the use of these established platforms beyond infectious diseases., Conclusions: VA SeqFORCE and SeqCURE were established as clinical and research programs dedicated to sequencing COVID-19 as part of ongoing clinical and surveillance efforts. In the future, we anticipate that having these programs embedded within the largest integrated health care system in the US will enable the study of pathogens and pandemics beyond COVID-19 and at an unprecedented scale. The investment in these programs will form an integral part of our nation's response to emerging infectious diseases, with future applications to precision medicine and beyond., Competing Interests: Author disclosures: VCM has received support from the Emory CFAR (P30 AI050409) and received investigator-initiated research grants (to the institution) and consultation fees (both unrelated to the current work) from Eli Lilly, Bayer, Gilead Sciences, and ViiV. CWW has a consulting relationship with Biomeme, Bavarian-Nordic, Pfizer, and Regeneron. CWW has also received research grants from Pfizer and Sanofi. All other authors report no actual or potential conflicts of interest or outside sources of funding with regard to this article., (Copyright © 2023 Frontline Medical Communications Inc., Parsippany, NJ, USA.)

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have a significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors; however, recent work suggests prevalence may differ between sub-groups. Here, we examined a large cohort of diverse adults in the Bio Me biobank in New York City. We observed the prevalence of PAD at 1.7% in EAs vs. 8.5% and 9.4% in AAs and HLs, respectively, and among HL sub-groups, the prevalence was found at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs [OR = 3.15 (95% CI 2.33-4.25), p < 6.44 × 10 -14 ]. To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry and PAD in Dominican Bio Me participants (N = 1,813) separately from European, African, and Native American (NAT) continental ancestry tracts. The top association with PAD was an NAT ancestry tract at chromosome 2q35 [OR = 1.96 (SE = 0.16), p < 2.75 × 10 -05 ) with 22.6% vs. 12.9% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607 , a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. Efforts to reproduce the signal in other Hispanic cohorts were unsuccessful. In summary, we showed how leveraging health system data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a putative risk locus in a Dominican population., Competing Interests: EK received personal fees from Illumina, 23andMe, Allelica, and Regeneron Pharmaceuticals, received research funding from Allelica, and serves as a scientific advisory board member for Encompass Bio, Foresite Labs, and Galateo Bio. NA-H is an employee and equity holder of 23andMe; serves as a scientific advisory board member for Allelica; received personal fees from Genentech, Allelica, and 23andMe; received research funding from Akcea; and was previously employed by Regeneron Pharmaceuticals. ES is an employee of Calico Life Sciences LLC. SA received consulting fees from Variant Bio and Biogen Inc. CG owns stock in 23andMe and serves as an advisory board member for Encompass Bio. DK serve as scientific advisor to Bitterroot Bio, Inc. SMD receives research support from RenalytixAI and Novo Nordisk, outside the scope of the current research. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Cullina, Wojcik, Shemirani, Klarin, Gorman, Sorokin, Gignoux, Belbin, Pyarajan, Asgari, Tsao, Damrauer, Abul-Husn and Kenny.)

Importance: The reported phenotypes of men with 47,XXY and 47,XYY syndromes include tall stature, multisystem comorbidities, and poor health-related quality of life (HRQoL). However, knowledge about these sex chromosome aneuploidy (SCA) conditions has been derived from studies in the <15% of patients who are clinically diagnosed and also lack diversity in age and genetic ancestry., Objectives: Determine the prevalence of clinically diagnosed and undiagnosed X or Y chromosome aneuploidy among men enrolled in the Million Veteran Program (MVP); describe military service metrics of men with SCAs; compare morbidity and mortality outcomes between men with SCA with and without a clinical diagnosis to matched controls., Design: Cross-sectional, case-control., Setting: United States Veterans Administration Healthcare System., Participants: Biologic males enrolled in the MVP biobank with genomic identification of an additional X or Y chromosome (cases); controls matched 1:5 on sex, age, and genetic ancestry., Main Outcomes and Measures: Prevalence of men with SCAs from genomic analysis; clinical SCA diagnosis; Charlson Comorbidity Index (CCI); rates of outpatient, inpatient, and emergency encounters per year; self-reported health outcomes; standardized mortality ratio (SMR)., Results: An additional X or Y chromosome was present in 145 and 125 per 100,000 males in the MVP, respectively, with the highest prevalence among men with European and East Asian ancestry. At a mean age of 61±12 years, 74% of male veterans with 47,XXY and >99% with 47,XYY remained undiagnosed. Individuals with 47,XXY (n=862) and 47,XYY (n=747) had similar military service history, all-cause SMR, and age of death compared to matched controls. CCI and healthcare utilization were higher among individuals with SCA, while several measures of HRQoL were lower. Men with a clinical diagnosis of 47,XXY had higher healthcare utilization but lower comorbidity score compared to those undiagnosed., Conclusion and Relevance: One in 370 males in the MVP cohort have SCA, a prevalence comparable to estimates in the general population. While these men have successfully served in the military, they have higher morbidity and report poorer HRQoL with aging. Longer longitudinal follow-up of this sample will be informative for clinical and patient-reported outcomes, the role of ancestry, and mortality statistics., Key Points: Comparable to the general population, approximately 1 in 370 male veterans have a sex chromosome aneuploidy, but most are undiagnosed.Men with X or Y chromosome aneuploidy successfully complete US miliary duty with similar service history compared to their 46,XY peers.Medical comorbidities and healthcare utilization metrics are higher in male veterans with 47,XXY and 47,XYY during aging, however life expectancy is similar to matched controls.

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8 + T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8 + T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs., Competing Interests: Declaration of interests G.P.M. receives research support from Thermo Fisher, CareDx, and PIRCHE. M.Z. is a board member of Invectys Inc., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

The current understanding of the genetic determinants of thoracic aortic aneurysms and dissections (TAAD) has largely been informed through studies of rare, Mendelian forms of disease. Here, we conducted a genome-wide association study (GWAS) of TAAD, testing ~25 million DNA sequence variants in 8,626 participants with and 453,043 participants without TAAD in the Million Veteran Program, with replication in an independent sample of 4,459 individuals with and 512,463 without TAAD from six cohorts. We identified 21 TAAD risk loci, 17 of which have not been previously reported. We leverage multiple downstream analytic methods to identify causal TAAD risk genes and cell types and provide human genetic evidence that TAAD is a non-atherosclerotic aortic disorder distinct from other forms of vascular disease. Our results demonstrate that the genetic architecture of TAAD mirrors that of other complex traits and that it is not solely inherited through protein-altering variants of large effect size., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Genome-wide association studies (GWAS) have underrepresented individuals from non-European populations, impeding progress in characterizing the genetic architecture and consequences of health and disease traits. To address this, we present a population-stratified phenome-wide GWAS followed by a multi-population meta-analysis for 2,068 traits derived from electronic health records of 635,969 participants in the Million Veteran Program (MVP), a longitudinal cohort study of diverse U.S. Veterans genetically similar to the respective African (121,177), Admixed American (59,048), East Asian (6,702), and European (449,042) superpopulations defined by the 1000 Genomes Project. We identified 38,270 independent variants associating with one or more traits at experiment-wide P < 4.6 × 10 - 11 significance; fine-mapping 6,318 signals identified from 613 traits to single-variant resolution. Among these, a third (2,069) of the associations were found only among participants genetically similar to non-European reference populations, demonstrating the importance of expanding diversity in genetic studies. Our work provides a comprehensive atlas of phenome-wide genetic associations for future studies dissecting the architecture of complex traits in diverse populations., Competing Interests: Competing interests CJD and JPC are employed full-time by the Novartis Institute of Biomedical Interest (their major contributions to this project were while employed at VA Boston Healthcare System). H.K. is a member of advisory boards for Dicerna Pharmaceuticals, Sophrosyne Pharmaceuticals; Enthion Pharmaceuticals; and Clearmind Medicine; a consultant to Sobrera Pharmaceuticals; the recipient of research funding and medication supplies for an investigator-initiated study from Alkermes member of the American Society of Clinical Psychopharmacology’s Alcohol Clinical Trials Initiative; which was supported in the last three years by Alkermes, Dicerna, Ethypharm, Lundbeck, Mitsubishi, Otsuka, and Pear Therapeutics; and holder of U.S. patent 10,900,082 titled: “Genotype-guided dosing of opioid agonists,” issued 26 January 2021. JG and RP are paid for their editorial work in the journal Complex Psychiatry. RP reports a research grant from Alkermes. SD reports grants from Alnylam Pharmaceuticals, Inc, grants from Astellas Pharma, Inc; grants from AstraZeneca Pharmaceuticals LP; grants from Biodesix, grants from Celgene Corporation; grants from Cerner Enviza; grants from GlaxoSmithKline PLC, grants from Janssen Pharmaceuticals, Inc.; grants from Kantar Health; grants from Myriad Genetic Laboratories, Inc.; grants from Novartis International AG; grants from Parexel International Corporation through the University of Utah or Western Institute for Veteran Research outside the submitted work. SMD receives research support from RenalytixAI and Novo Nordisk, outside the scope of the current research, and is named as a co-inventor on a Government-owned US Patent application related to the use of genetic risk prediction for venous thromboembolic disease and for the use of PDE3B inhibition for preventing cardiovascular disease, both filed by the US Department of Veterans Affairs in accordance with Federal regulatory requirements.

Background: The Million Veteran Program (MVP) participants represent 100 years of US history, including significant social and demographic changes over time. Our study assessed two aspects of the MVP: (i) longitudinal changes in population diversity and (ii) how these changes can be accounted for in genome-wide association studies (GWAS). To investigate these aspects, we divided MVP participants into five birth cohorts (N-range = 123,888 [born from 1943 to 1947] to 136,699 [born from 1948 to 1953])., Results: Ancestry groups were defined by (i) HARE (harmonized ancestry and race/ethnicity) and (ii) a random-forest clustering approach using the 1000 Genomes Project and the Human Genome Diversity Project (1kGP + HGDP) reference panels (77 world populations representing six continental groups). In these groups, we performed GWASs of height, a trait potentially affected by population stratification. Birth cohorts demonstrate important trends in ancestry diversity over time. More recent HARE-assigned Europeans, Africans, and Hispanics had lower European ancestry proportions than older birth cohorts (0.010 < Cohen's d < 0.259, p < 7.80 × 10 -4 ). Conversely, HARE-assigned East Asians showed an increase in European ancestry proportion over time. In GWAS of height using HARE assignments, genomic inflation due to population stratification was prevalent across all birth cohorts (linkage disequilibrium score regression intercept = 1.08 ± 0.042). The 1kGP + HGDP-based ancestry assignment significantly reduced the population stratification (mean intercept reduction = 0.045 ± 0.007, p < 0.05) confounding in the GWAS statistics., Conclusions: This study provides a characterization of ancestry diversity of the MVP cohort over time and compares two strategies to infer genetically defined ancestry groups by assessing differences in controlling population stratification in genome-wide association studies., (© 2023. The Author(s).)

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors, however recent work suggests prevalence may differ between sub-groups. Here we examined a large cohort of diverse adults in the Bio Me biobank in New York City (NYC). We observed the prevalence of PAD at 1.7% in EAs vs 8.5% and 9.4% in AAs and HLs, respectively; and among HL sub-groups, at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs (OR=3.15 (95% CI 2.33-4.25), P <6.44×10 -14 ). To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry (LA) and PAD in Dominican Bio Me participants (N=1,940) separately for European (EUR), African (AFR) and Native American (NAT) continental ancestry tracts. We identified a NAT ancestry tract at chromosome 2q35 that was significantly associated with PAD (OR=2.05 (95% CI 1.51-2.78), P <4.06×10 -6 ) with 22.5% vs 12.5% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607 , a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. In summary, we showed how leveraging health systems data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a novel risk locus in a Dominican population.

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.

We present allele frequencies of pharmacogenomics relevant variants across multiple ancestry in a sample representative of the US population. We analyzed 658,582 individuals with genotype data and extracted pharmacogenomics relevant single nucleotide variant (SNV) alleles, human leukocyte antigens (HLA) 4-digit alleles and an important copy number variant (CNV), the full deletion/duplication of CYP2D6. We compiled distinct allele frequency tables for European, African American, Hispanic, and Asian ancestry individuals. In addition, we compiled allele frequencies based on local ancestry reconstruction in the African-American (2-way deconvolution) and Hispanic (3-way deconvolution) cohorts., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

Increasing rates of human papillomavirus (HPV)-driven oropharyngeal cancer (OPC) have largely offset declines in tobacco-associated head and neck squamous cell carcinoma (HNSCC) at non-OPC sites. Host immunity is an important modulator of HPV infection, persistence, and clearance, and also of immune evasion in both virally- and nonvirally-driven cancers. However, the association between collective known cancer-related immune gene variants and HNSCC susceptibility has not been fully characterized. Here, we conducted a genetic association study in the multiethnic Veterans Affairs Million Veteran Program cohort, evaluating 16,050 variants in 1,576 immune genes in 4,012 HNSCC cases (OPC = 1,823; non-OPC = 2,189) and 16,048 matched controls. Significant polymorphisms were further examined in a non-Hispanic white (NHW) validation cohort (OPC = 1,206; non-OPC = 955; controls = 4,507). For overall HNSCC susceptibility in NHWs, we discovered and validated a novel 9q31.1 SMC2 association and replicated the known 6p21.32 HLA-DQ-DR association. Six loci/genes for overall HNSCC susceptibility were selectively enriched in African-Americans (6p21.32 HLA-G, 9q21.33 GAS1, 11q12.2 CD6, 11q23.2 NCAM1/CD56, 17p13.1 CD68, 18q22.2 SOCS6); all 6 genes function in antigen-presenting regulation and T-cell activation. Two additional loci (10q26 DMBT1, 15q22.2 TPM1) were uncovered for non-OPC susceptibility, and three loci (11q24 CRTAM, 16q21 CDH5, 18q12.1 CDH2) were identified for HPV-positive OPC susceptibility. This study underscores the role of immune gene variants in modulating susceptibility for both HPV-driven and non-HPV-driven HNSCC. Additional large studies, particularly in racially diverse populations, are needed to further validate the associations and to help elucidate other potential immune factors and mechanisms that may underlie HNSCC risk., Significance: Several inherited variations in immune system genes are significantly associated with susceptibility to head and neck cancer, which could help improve personalized cancer risk estimates., (©2022 The Authors; Published by the American Association for Cancer Research.)

Mosaic loss of the X chromosome (mLOX) is the most commonly occurring clonal somatic alteration detected in the leukocytes of women, yet little is known about its genetic determinants or phenotypic consequences. To address this, we estimated mLOX in >900,000 women across eight biobanks, identifying 10% of women with detectable X loss in approximately 2% of their leukocytes. Out of 1,253 diseases examined, women with mLOX had an elevated risk of myeloid and lymphoid leukemias and pneumonia. Genetic analyses identified 49 common variants influencing mLOX, implicating genes with established roles in chromosomal missegregation, cancer predisposition, and autoimmune diseases. Complementary exome-sequence analyses identified rare missense variants in FBXO10 which confer a two-fold increased risk of mLOX. A small fraction of these associations were shared with mosaic Y chromosome loss in men, suggesting different biological processes drive the formation and clonal expansion of sex chromosome missegregation events. Allelic shift analyses identified alleles on the X chromosome which are preferentially retained, demonstrating that variation at many loci across the X chromosome is under cellular selection. A novel polygenic score including 44 independent X chromosome allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Collectively our results support a model where germline variants predispose women to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of subsequent clonal expansion.

Background: The coronavirus disease 2019 (COVID-19) pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has demonstrated the need to share data and biospecimens broadly to optimize clinical outcomes for US military Veterans., Methods: In response, the Veterans Health Administration established VA SHIELD (Science and Health Initiative to Combat Infectious and Emerging Life-threatening Diseases), a comprehensive biorepository of specimens and clinical data from affected Veterans to advance research and public health surveillance and to improve diagnostic and therapeutic capabilities., Results: VA SHIELD now comprises 12 sites collecting de-identified biospecimens from US Veterans affected by SARS-CoV-2. In addition, 2 biorepository sites, a data processing center, and a coordinating center have been established under the direction of the Veterans Affairs Office of Research and Development. Phase 1 of VA SHIELD comprises 34 157 samples. Of these, 83.8% had positive tests for SARS-CoV-2, with the remainder serving as contemporaneous controls. The samples include nasopharyngeal swabs (57.9%), plasma (27.9%), and sera (12.5%). The associated clinical and demographic information available permits the evaluation of biological data in the context of patient demographics, clinical experience and management, vaccinations, and comorbidities., Conclusions: VA SHIELD is representative of US national diversity with a significant potential to impact national healthcare. VA SHIELD will support future projects designed to better understand SARS-CoV-2 and other emergent healthcare crises. To the extent possible, VA SHIELD will facilitate the discovery of diagnostics and therapeutics intended to diminish COVID-19 morbidity and mortality and to reduce the impact of new emerging threats to the health of US Veterans and populations worldwide., Competing Interests: Potential conflicts of interest. The authors have noo reported conflicts of interest., (© The Author(s) 2022. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Osteoarthritis is a common progressive joint disease. As no effective medical interventions are available, osteoarthritis often progresses to the end stage, in which only surgical options such as total joint replacement are available. A more thorough understanding of genetic influences of osteoarthritis is essential to develop targeted personalized approaches to treatment, ideally long before the end stage is reached. To date, there have been no large multiancestry genetic studies of osteoarthritis. Here, we leveraged the unique resources of 484,374 participants in the Million Veteran Program and UK Biobank to address this gap. Analyses included participants of European, African, Asian and Hispanic descent. We discovered osteoarthritis-associated genetic variation at 10 loci and replicated findings from previous osteoarthritis studies. We also present evidence that some osteoarthritis-associated regions are robust to population ancestry. Drug repurposing analyses revealed enrichment of targets of several medication classes and provide potential insight into the etiology of beneficial effects of antiepileptics on osteoarthritis pain., (© 2022. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

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

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (N cases  = 4,325; N controls  = 507,640; OR = 0.89 [0.82-0.97]; P  = 6.86 × 10 -3 ) and joint meta-analyses with the HGI (N cases  = 13,320; N controls  = 1,508,841; OR, 0.90 [0.86-0.95]; P  = 8.99 × 10 -5 ). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (N cases  = 19,168/N controls  = 492,854; OR, 0.98 [0.95-1.01]; P  = 0.06) but was nominally significant ( P  < 0.05) in the joint meta-analysis with the HGI (N cases  = 44,820; N controls  = 1,775,827; OR, 0.97 [0.95-1.00]; P  = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P  = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes 1 . Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel 2 ) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries., (© 2022. The Author(s).)

Importance: Serious mental illnesses, including schizophrenia, bipolar disorder, and depression, are heritable, highly multifactorial disorders and major causes of disability worldwide., Objective: To benchmark the penetrance of current neuropsychiatric polygenic risk scores (PRSs) in the Veterans Health Administration health care system and to explore associations between PRS and broad categories of human disease via phenome-wide association studies., Design, Setting, and Participants: Extensive Veterans Health Administration's electronic health records were assessed from October 1999 to January 2021, and an embedded cohort of 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder were found. The performance of schizophrenia, bipolar disorder, and major depression PRSs were compared in participants of African or European ancestry in the Million Veteran Program (approximately 400 000 individuals), and associations between PRSs and 1650 disease categories based on ICD-9/10 billing codes were explored. Last, genomic structural equation modeling was applied to derive novel PRSs indexing common and disorder-specific genetic factors. Analysis took place from January 2021 to January 2022., Main Outcomes and Measures: Diagnoses based on in-person structured clinical interviews were compared with ICD-9/10 billing codes. PRSs were constructed using summary statistics from genome-wide association studies of schizophrenia, bipolar disorder, and major depression., Results: Of 707 299 enrolled study participants, 459 667 were genotyped at the time of writing; 84 806 were of broadly African ancestry (mean [SD] age, 58 [12.1] years) and 314 909 were of broadly European ancestry (mean [SD] age, 66.4 [13.5] years). Among 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder, 8962 (95.6%) were correctly identified using ICD-9/10 codes (2 or more). Among those of European ancestry, PRSs were robustly associated with having received a diagnosis of schizophrenia (odds ratio [OR], 1.81 [95% CI, 1.76-1.87]; P < 10-257) or bipolar disorder (OR, 1.42 [95% CI, 1.39-1.44]; P < 10-295). Corresponding effect sizes in participants of African ancestry were considerably smaller for schizophrenia (OR, 1.35 [95% CI, 1.29-1.42]; P < 10-38) and bipolar disorder (OR, 1.16 [95% CI, 1.11-1.12]; P < 10-10). Neuropsychiatric PRSs were associated with increased risk for a range of psychiatric and physical health problems., Conclusions and Relevance: Using diagnoses confirmed by in-person structured clinical interviews and current neuropsychiatric PRSs, the validity of an electronic health records-based phenotyping approach in US veterans was demonstrated, highlighting the potential of PRSs for disentangling biological and mediated pleiotropy.

We report a genome-wide association study (GWAS) of coronary artery disease (CAD) incorporating nearly a quarter of a million cases, in which existing studies are integrated with data from cohorts of white, Black and Hispanic individuals from the Million Veteran Program. We document near equivalent heritability of CAD across multiple ancestral groups, identify 95 novel loci, including nine on the X chromosome, detect eight loci of genome-wide significance in Black and Hispanic individuals, and demonstrate that two common haplotypes at the 9p21 locus are responsible for risk stratification in all populations except those of African origin, in which these haplotypes are virtually absent. Moreover, in the largest GWAS for angiographically derived coronary atherosclerosis performed to date, we find 15 loci of genome-wide significance that robustly overlap with established loci for clinical CAD. Phenome-wide association analyses of novel loci and polygenic risk scores (PRSs) augment signals related to insulin resistance, extend pleiotropic associations of these loci to include smoking and family history, and precisely document the markedly reduced transferability of existing PRSs to Black individuals. Downstream integrative analyses reinforce the critical roles of vascular endothelial, fibroblast, and smooth muscle cells in CAD susceptibility, but also point to a shared biology between atherosclerosis and oncogenesis. This study highlights the value of diverse populations in further characterizing the genetic architecture of CAD., (© 2022. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Importance: Sickle cell trait (SCT), defined as the presence of 1 hemoglobin beta sickle allele (rs334-T) and 1 normal beta allele, is prevalent in millions of people in the US, particularly in individuals of African and Hispanic ancestry. However, the association of SCT with COVID-19 is unclear., Objective: To assess the association of SCT with the prepandemic health conditions in participants of the Million Veteran Program (MVP) and to assess the severity and sequelae of COVID-19., Design, Setting, and Participants: COVID-19 clinical data include 2729 persons with SCT, of whom 353 had COVID-19, and 129 848 SCT-negative individuals, of whom 13 488 had COVID-19. Associations between SCT and COVID-19 outcomes were examined using firth regression. Analyses were performed by ancestry and adjusted for sex, age, age squared, and ancestral principal components to account for population stratification. Data for the study were collected between March 2020 and February 2021., Exposures: The hemoglobin beta S (HbS) allele (rs334-T)., Main Outcomes and Measures: This study evaluated 4 COVID-19 outcomes derived from the World Health Organization severity scale and phenotypes derived from International Classification of Diseases codes in the electronic health records., Results: Of the 132 577 MVP participants with COVID-19 data, mean (SD) age at the index date was 64.8 (13.1) years. Sickle cell trait was present in 7.8% of individuals of African ancestry and associated with a history of chronic kidney disease, diabetic kidney disease, hypertensive kidney disease, pulmonary embolism, and cerebrovascular disease. Among the 4 clinical outcomes of COVID-19, SCT was associated with an increased COVID-19 mortality in individuals of African ancestry (n = 3749; odds ratio, 1.77; 95% CI, 1.13 to 2.77; P = .01). In the 60 days following COVID-19, SCT was associated with an increased incidence of acute kidney failure. A counterfactual mediation framework estimated that on average, 20.7% (95% CI, -3.8% to 56.0%) of the total effect of SCT on COVID-19 fatalities was due to acute kidney failure., Conclusions and Relevance: In this genetic association study, SCT was associated with preexisting kidney comorbidities, increased COVID-19 mortality, and kidney morbidity.

Background: Height has been associated with many clinical traits but whether such associations are causal versus secondary to confounding remains unclear in many cases. To systematically examine this question, we performed a Mendelian Randomization-Phenome-wide association study (MR-PheWAS) using clinical and genetic data from a national healthcare system biobank., Methods and Findings: Analyses were performed using data from the US Veterans Affairs (VA) Million Veteran Program in non-Hispanic White (EA, n = 222,300) and non-Hispanic Black (AA, n = 58,151) adults in the US. We estimated height genetic risk based on 3290 height-associated variants from a recent European-ancestry genome-wide meta-analysis. We compared associations of measured and genetically-predicted height with phenome-wide traits derived from the VA electronic health record, adjusting for age, sex, and genetic principal components. We found 345 clinical traits associated with measured height in EA and an additional 17 in AA. Of these, 127 were associated with genetically-predicted height at phenome-wide significance in EA and 2 in AA. These associations were largely independent from body mass index. We confirmed several previously described MR associations between height and cardiovascular disease traits such as hypertension, hyperlipidemia, coronary heart disease (CHD), and atrial fibrillation, and further uncovered MR associations with venous circulatory disorders and peripheral neuropathy in the presence and absence of diabetes. As a number of traits associated with genetically-predicted height frequently co-occur with CHD, we evaluated effect modification by CHD status of genetically-predicted height associations with risk factors for and complications of CHD. We found modification of effects of MR associations by CHD status for atrial fibrillation/flutter but not for hypertension, hyperlipidemia, or venous circulatory disorders., Conclusions: We conclude that height may be an unrecognized but biologically plausible risk factor for several common conditions in adults. However, more studies are needed to reliably exclude horizontal pleiotropy as a driving force behind at least some of the MR associations observed in this study., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: CJO is a full-time employee of Novartis Institutes of Biomedical Research. The remaining authors have declared that no competing interests exist.

Nonalcoholic fatty liver disease (NAFLD) is a growing cause of chronic liver disease. Using a proxy NAFLD definition of chronic elevation of alanine aminotransferase (cALT) levels without other liver diseases, we performed a multiancestry genome-wide association study (GWAS) in the Million Veteran Program (MVP) including 90,408 cALT cases and 128,187 controls. Seventy-seven loci exceeded genome-wide significance, including 25 without prior NAFLD or alanine aminotransferase associations, with one additional locus identified in European American-only and two in African American-only analyses (P < 5 × 10 -8 ). External replication in histology-defined NAFLD cohorts (7,397 cases and 56,785 controls) or radiologic imaging cohorts (n = 44,289) replicated 17 single-nucleotide polymorphisms (SNPs) (P < 6.5 × 10 -4 ), of which 9 were new (TRIB1, PPARG, MTTP, SERPINA1, FTO, IL1RN, COBLL1, APOH and IFI30). Pleiotropy analysis showed that 61 of 77 multiancestry and all 17 replicated SNPs were jointly associated with metabolic and/or inflammatory traits, revealing a complex model of genetic architecture. Our approach integrating cALT, histology and imaging reveals new insights into genetic liability to NAFLD., (© 2022. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)