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

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

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

Objective: The development of clinical research informatics tools and workflow processes associated with re-engaging biobank participants has become necessary as genomic repositories increasingly consider the return of actionable research results. Materials and Methods: Here we describe the development and utility of an informatics application for participant recruitment and enrollment management for the Veterans Affairs Million Veteran Program Return Of Actionable Results Study, a randomized controlled pilot trial returning individual genetic results associated with familial hypercholesterolemia. Results: The application is developed in Python-Flask and was placed into production in November 2021. The application includes modules for chart review, medication reconciliation, participant contact and biospecimen logging, survey recording, randomization, and documentation of genetic counseling and result disclosure. Three primary users, a genetic counselor and two research coordinators, and 326 Veteran participants have been integrated into the system as of February 23, 2023. The application has successfully handled 3367 task requests involving greater than 95 000 structured data points. Specifically, application users have recorded 326 chart reviews, 867 recruitment telephone calls, 158 telephone-based surveys, and 61 return of results genetic counseling sessions, among other available study tasks. Conclusion: The development of usable, customizable, and secure informatics tools will become increasingly important as large genomic repositories begin to return research results at scale. Our work provides a proof-of-concept for developing and using such tools to aid in managing the return of results process within a national biobank.

Diversity and Scale: Genetic Architecture of 2,068 Traits in the VA Million Veteran Program.

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.

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

Objective: Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed mental and physical health problems, prior pharmacological treatments, and aggregate genetic factors has potential to inform risk stratification and mitigation strategies. Methods: In this study of 3,942 SCZ and 5,414 BPI patients receiving VA care, self-reported SB and ideation were assessed using the Columbia Suicide Severity Rating Scale (C-SSRS). These cross-sectional data were integrated with electronic health records (EHR), and compared by lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. Polygenic scores (PGS) for traits related to psychiatric disorders, substance use, and cognition were constructed using available genomic data, and exploratory genome-wide association studies were performed to identify and prioritize specific loci. Results: Only 20% of veterans who self-reported SB had a corroborating ICD-9/10 code in their EHR; and among those who denied prior behaviors, more than 20% reported new-onset SB at follow-up. SB were associated with a range of psychiatric and non-psychiatric diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking, suicide attempt, and major depressive disorder were also associated with attempt and ideation. Conclusions: Among individuals with a diagnosed mental illness, a GWAS for SB did not yield any significant loci. Self-reported SB were strongly associated with clinical variables across several EHR domains. Overall, clinical and polygenic analyses point to sequelae of substance-use related behaviors and other psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in clinical settings, underscoring the value of regular screening based on direct, in-person assessments, especially among high-risk individuals.

Pharmacogenetic allele variant frequencies: An analysis of the VA's Million Veteran Program (MVP) as a representation of the diversity in US population.

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.

Immunogenetic Determinants of Susceptibility to Head and Neck Cancer in the Million Veteran Program Cohort.

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.

Penetrance and Pleiotropy of Polygenic Risk Scores for Schizophrenia, Bipolar Disorder, and Depression Among Adults in the US Veterans Affairs Health Care System.

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.

Actionable druggable genome-wide Mendelian randomization identifies repurposing opportunities for COVID-19.

Drug repurposing provides a rapid approach to meet the urgent need for therapeutics to address COVID-19. To identify therapeutic targets relevant to COVID-19, we conducted Mendelian randomization analyses, deriving genetic instruments based on transcriptomic and proteomic data for 1,263 actionable proteins that are targeted by approved drugs or in clinical phase of drug development. Using summary statistics from the Host Genetics Initiative and the Million Veteran Program, we studied 7,554 patients hospitalized with COVID-19 and >1 million controls. We found significant Mendelian randomization results for three proteins (ACE2, P = 1.6 × 10-6; IFNAR2, P = 9.8 × 10-11 and IL-10RB, P = 2.3 × 10-14) using cis-expression quantitative trait loci genetic instruments that also had strong evidence for colocalization with COVID-19 hospitalization. To disentangle the shared expression quantitative trait loci signal for IL10RB and IFNAR2, we conducted phenome-wide association scans and pathway enrichment analysis, which suggested that IFNAR2 is more likely to play a role in COVID-19 hospitalization. Our findings prioritize trials of drugs targeting IFNAR2 and ACE2 for early management of COVID-19.

Genome-Wide Association Studies of Schizophrenia and Bipolar Disorder in a Diverse Cohort of US Veterans.

BACKGROUND: Schizophrenia (SCZ) and bipolar disorder (BIP) are debilitating neuropsychiatric disorders, collectively affecting 2% of the world's population. Recognizing the major impact of these psychiatric disorders on the psychosocial function of more than 200 000 US Veterans, the Department of Veterans Affairs (VA) recently completed genotyping of more than 8000 veterans with SCZ and BIP in the Cooperative Studies Program (CSP) #572. METHODS: We performed genome-wide association studies (GWAS) in CSP #572 and benchmarked the predictive value of polygenic risk scores (PRS) constructed from published findings. We combined our results with available summary statistics from several recent GWAS, realizing the largest and most diverse studies of these disorders to date. RESULTS: Our primary GWAS uncovered new associations between CHD7 variants and SCZ, and novel BIP associations with variants in Sortilin Related VPS10 Domain Containing Receptor 3 (SORCS3) and downstream of PCDH11X. Combining our results with published summary statistics for SCZ yielded 39 novel susceptibility loci including CRHR1, and we identified 10 additional findings for BIP (28 326 cases and 90 570 controls). PRS trained on published GWAS were significantly associated with case-control status among European American (P < 10-30) and African American (P < .0005) participants in CSP #572. CONCLUSIONS: We have demonstrated that published findings for SCZ and BIP are robustly generalizable to a diverse cohort of US veterans. Leveraging available summary statistics from GWAS of global populations, we report 52 new susceptibility loci and improved fine-mapping resolution for dozens of previously reported associations.

Genotyping Array Design and Data Quality Control in the Million Veteran Program.

The Million Veteran Program (MVP), initiated by the Department of Veterans Affairs (VA), aims to collect biosamples with consent from at least one million veterans. Presently, blood samples have been collected from over 800,000 enrolled participants. The size and diversity of the MVP cohort, as well as the availability of extensive VA electronic health records, make it a promising resource for precision medicine. MVP is conducting array-based genotyping to provide a genome-wide scan of the entire cohort, in parallel with whole-genome sequencing, methylation, and other 'omics assays. Here, we present the design and performance of the MVP 1.0 custom Axiom array, which was designed and developed as a single assay to be used across the multi-ethnic MVP cohort. A unified genetic quality-control analysis was developed and conducted on an initial tranche of 485,856 individuals, leading to a high-quality dataset of 459,777 unique individuals. 668,418 genetic markers passed quality control and showed high-quality genotypes not only on common variants but also on rare variants. We confirmed that, with non-European individuals making up nearly 30%, MVP's substantial ancestral diversity surpasses that of other large biobanks. We also demonstrated the quality of the MVP dataset by replicating established genetic associations with height in European Americans and African Americans ancestries. This current dataset has been made available to approved MVP researchers for genome-wide association studies and other downstream analyses. Further data releases will be available for analysis as recruitment at the VA continues and the cohort expands both in size and diversity.

Genome-wide association study of post-traumatic stress disorder reexperiencing symptoms in >165,000 US veterans.

Post-traumatic stress disorder (PTSD) is a major problem among military veterans and civilians alike, yet its pathophysiology remains poorly understood. We performed a genome-wide association study and bioinformatic analyses, which included 146,660 European Americans and 19,983 African Americans in the US Million Veteran Program, to identify genetic risk factors relevant to intrusive reexperiencing of trauma, which is the most characteristic symptom cluster of PTSD. In European Americans, eight distinct significant regions were identified. Three regions had values of P < 5 × 10-10: CAMKV; chromosome 17 closest to KANSL1, but within a large high linkage disequilibrium region that also includes CRHR1; and TCF4. Associations were enriched with respect to the transcriptomic profiles of striatal medium spiny neurons. No significant associations were observed in the African American cohort of the sample. Results in European Americans were replicated in the UK Biobank data. These results provide new insights into the biology of PTSD in a well-powered genome-wide association study.

Genome-wide Association Study of Maximum Habitual Alcohol Intake in >140,000 U.S. European and African American Veterans Yields Novel Risk Loci.

BACKGROUND: Habitual alcohol use can be an indicator of alcohol dependence, which is associated with a wide range of serious health problems. METHODS: We completed a genome-wide association study in 126,936 European American and 17,029 African American subjects in the Veterans Affairs Million Veteran Program for a quantitative phenotype based on maximum habitual alcohol consumption. RESULTS: ADH1B, on chromosome 4, was the lead locus for both populations: for the European American sample, rs1229984 (p = 4.9 × 10-47); for African American, rs2066702 (p = 2.3 × 10-12). In the European American sample, we identified three additional genome-wide-significant maximum habitual alcohol consumption loci: on chromosome 17, rs77804065 (p = 1.5 × 10-12), at CRHR1 (corticotropin-releasing hormone receptor 1); the protein product of this gene is involved in stress and immune responses; and on chromosomes 8 and 10. European American and African American samples were then meta-analyzed; the associated region at CRHR1 increased in significance to 1.02 × 10-13, and we identified two additional genome-wide significant loci, FGF14 (p = 9.86 × 10-9) (chromosome 13) and a locus on chromosome 11. Besides ADH1B, none of the five loci have prior genome-wide significant support. Post-genome-wide association study analysis identified genetic correlation to other alcohol-related traits, smoking-related traits, and many others. Replications were observed in UK Biobank data. Genetic correlation between maximum habitual alcohol consumption and alcohol dependence was 0.87 (p = 4.78 × 10-9). Enrichment for cell types included dopaminergic and gamma-aminobutyric acidergic neurons in midbrain, and pancreatic delta cells. CONCLUSIONS: The present study supports five novel alcohol-use risk loci, with particularly strong statistical support for CRHR1. Additionally, we provide novel insight regarding the biology of harmful alcohol use.

Clinical and Molecular Characteristics of NF1-Mutant Lung Cancer.

PURPOSE: NF1 is a tumor suppressor that negatively regulates Ras signaling. NF1 mutations occur in lung cancer, but their clinical significance is unknown. We evaluated clinical and molecular characteristics of NF1 mutant lung cancers with comparison to tumors with KRAS mutations. EXPERIMENTAL DESIGN: Between July 2013 and October 2014, 591 non-small cell lung cancer (NSCLC) tumors underwent targeted next-generation sequencing in a 275 gene panel that evaluates gene mutations and genomic rearrangements. NF1 and KRAS cohorts were identified, with subsequent clinical and genomic analysis. RESULTS: Among 591 patients, 60 had NF1 mutations (10%) and 141 (24%) had KRAS mutations. 15 NF1 mutations (25%) occurred with other oncogenic mutations [BRAF (2); ERBB2 (2); KRAS (9); HRAS (1); NRAS (1)]. There were 72 unique NF1 variants. NF1 tumor pathology was diverse, including both adenocarcinoma (36, 60%) and squamous cell carcinoma (10, 17%). In contrast, KRAS mutations occurred predominantly in adenocarcinoma (136, 96%). Both mutations were common in former/current smokers. Among NF1 tumors without concurrent oncogenic alterations, TP53 mutations/2-copy deletions occurred more often (33, 65%) than with KRAS mutation (46, 35%; P < 0.001). No difference between cohorts was seen with other tumor suppressors. CONCLUSIONS: NF1 mutations define a unique population of NSCLC. NF1 and KRAS mutations present in similar patient populations, but NF1 mutations occur more often with other oncogenic alterations and TP53 mutations. Therapeutic strategies targeting KRAS activation, including inhibitors of MAP kinase signaling, may warrant investigation in NF1 mutant tumors. Tumor-suppressor inactivation patterns may help further define novel treatment strategies. Clin Cancer Res; 22(13); 3148-56. ©2016 AACR.

Genetic Basis for PD-L1 Expression in Squamous Cell Carcinomas of the Cervix and Vulva.

IMPORTANCE: Patients with squamous cell carcinoma (SCC) of the cervix or vulva have limited therapeutic options, and the potential for immunotherapy for this population has not been evaluated. Recent trials suggest that tumors with a genetic basis for PD-1 (programmed cell death protein 1) ligand expression are highly sensitive to therapeutic antibodies targeting PD-1. OBJECTIVE: To determine the genetic status of CD274 (encoding PD-L1 [programmed cell death 1 ligand 1]) and PDCD1LG2 (encoding PD-L2 [programmed cell death 1 ligand 2]) in SCCs of the cervix and vulva and to correlate the findings with PD-L1 protein expression. DESIGN, SETTING, AND PARTICIPANTS: We performed fluorescence in situ hybridization (FISH) using probes targeting CD274, PDCD1LG2, and the centromeric portion of chromosome 9, and immunohistochemistry (IHC) using an antibody recognizing PD-L1 on formalin-fixed, paraffin-embedded (FFPE) biopsy specimens from 48 cervical SCCs and 23 vulvar SCCs. MAIN OUTCOMES AND MEASURES: Tumors were categorized according to the genetic abnormality in CD274 and PDCD1LG2 (coamplification > cogain > polysomy > disomy) as detected by FISH, and evaluated on a semiquantitative scale (modified H score, the product of the percentage of tumor cells with positive staining and the maximum intensity of positive staining) for PD-L1 protein expression as detected by IHC. RESULTS: Overall, 71 samples of FFPE tissue from cases of cervical SCCs (n = 48) and vulvar SCCs (n = 23) were retrieved from the archives of Brigham and Women's Hospital and included in this study. We observed cogain or coamplification of CD274 and PDCD1LG2 in 32 of 48 cervical SCCs (67%) and 10 of 23 vulvar SCCs (43%). Median PD-L1 protein expression was highest among tumors with CD274 and PDCD1LG2 coamplification and lowest among tumors with disomy. CONCLUSIONS AND RELEVANCE: Recurrent copy number gain of the genes encoding the PD-1 ligands provides a genetic basis for PD-L1 expression in a subset of cervical and vulvar SCCs and identifies a class of patients that are rational candidates for therapies targeting PD-1.

OncDRS: An integrative clinical and genomic data platform for enabling translational research and precision medicine.

We live in the genomic era of medicine, where a patient's genomic/molecular data is becoming increasingly important for disease diagnosis, identification of targeted therapy, and risk assessment for adverse reactions. However, decoding the genomic test results and integrating it with clinical data for retrospective studies and cohort identification for prospective clinical trials is still a challenging task. In order to overcome these barriers, we developed an overarching enterprise informatics framework for translational research and personalized medicine called Synergistic Patient and Research Knowledge Systems (SPARKS) and a suite of tools called Oncology Data Retrieval Systems (OncDRS). OncDRS enables seamless data integration, secure and self-navigated query and extraction of clinical and genomic data from heterogeneous sources. Within a year of release, the system has facilitated more than 1500 research queries and has delivered data for more than 50 research studies.

Identification of miRNA changes in Alzheimer's disease brain and CSF yields putative biomarkers and insights into disease pathways.

MicroRNAs have essential functional roles in brain development and neuronal specification but their roles in neurodegenerative diseases such as Alzheimer's disease (AD) is unknown. Using a sensitive qRT-PCR platform we identified regional and stage-specific deregulation of miRNA expression in AD patient brains. We used experimental validation in addition to literature to reveal how the deregulated brain microRNAs are biomarkers for known and novel pathways in AD pathogenesis related to amyloid processing, neurogenesis, insulin resistance, and innate immunity. We additionally recovered miRNAs from cerebrospinal fluid and discovered AD-specific miRNA changes consistent with their role as potential biomarkers of disease.