Genetics of substance use disorders in the era of big data.

Substance use disorders (SUDs) are conditions in which the use of legal or illegal substances, such as nicotine, alcohol or opioids, results in clinical and functional impairment. SUDs and, more generally, substance use are genetically complex traits that are enormously costly on an individual and societal basis. The past few years have seen remarkable progress in our understanding of the genetics, and therefore the biology, of substance use and abuse. Various studies - including of well-defined phenotypes in deeply phenotyped samples, as well as broadly defined phenotypes in meta-analysis and biobank samples - have revealed multiple risk loci for these common traits. A key emerging insight from this work establishes a biological and genetic distinction between quantity and/or frequency measures of substance use (which may involve low levels of use without dependence), versus symptoms related to physical dependence.

Genetic influences and causal pathways shared between cannabis use disorder and other substance use traits.

Cannabis use disorder (CanUD) has increased with the legalization of the use of cannabis. Around 20% of individuals using cannabis develop CanUD, and the number of users has grown with increasing ease of access. CanUD and other substance use disorders (SUDs) are associated phenotypically and genetically. We leveraged new CanUD genomics data to undertake genetically-informed analyses with unprecedented power, to investigate the genetic architecture and causal relationships between CanUD and lifetime cannabis use with risk for developing SUDs and substance use traits. Analyses included calculating local and global genetic correlations, genomic structural equation modeling (genomicSEM), and Mendelian Randomization (MR). Results from the genetic correlation and genomicSEM analyses demonstrated that CanUD and cannabis use differ in their relationships with SUDs and substance use traits. We found significant causal effects of CanUD influencing all the analyzed traits: opioid use disorder (OUD) (Inverse variant weighted, IVW ß = 0.925 ± 0.082), problematic alcohol use (PAU) (IVW ß = 0.443 ± 0.030), drinks per week (DPW) (IVW ß = 0.182 ± 0.025), Fagerström Test for Nicotine Dependence (FTND) (IVW ß = 0.183 ± 0.052), cigarettes per day (IVW ß = 0.150 ± 0.045), current versus former smokers (IVW ß = 0.178 ± 0.052), and smoking initiation (IVW ß = 0.405 ± 0.042). We also found evidence of bidirectionality showing that OUD, PAU, smoking initiation, smoking cessation, and DPW all increase risk of developing CanUD. For cannabis use, bidirectional relationships were inferred with PAU, smoking initiation, and DPW; cannabis use was also associated with a higher risk of developing OUD (IVW ß = 0.785 ± 0.266). GenomicSEM confirmed that CanUD and cannabis use load onto different genetic factors. We conclude that CanUD and cannabis use can increase the risk of developing other SUDs. This has substantial public health implications; the move towards legalization of cannabis use may be expected to increase other kinds of problematic substance use. These harmful outcomes are in addition to the medical harms associated directly with CanUD.

Pleiotropy and genetically inferred causality linking multisite chronic pain to substance use disorders.

Individuals suffering from chronic pain develop substance use disorders (SUDs) more often than others. Understanding the shared genetic influences underlying the comorbidity between chronic pain and SUDs will lead to a greater understanding of their biology. Genome-wide association statistics were obtained from the UK Biobank for multisite chronic pain (MCP, Neffective = 387,649) and from the Million Veteran Program and the Psychiatric Genomics Consortium meta-analyses for alcohol use disorder (AUD, Neffective = 296,974), cannabis use disorder (CanUD, Neffective = 161,053), opioid use disorder (OUD, Neffective = 57,120), and problematic tobacco use (PTU, Neffective = 270,120). SNP-based heritability was estimated for each of the traits and genetic correlation (rg) analyses were performed to assess MCP-SUD pleiotropy. Bidirectional Mendelian Randomization analyses evaluated possible causal relationships. Finally, to identify and characterize individual loci, we performed a genome-wide pleiotropy analysis and a brain-wide analysis using imaging phenotypes available from the UK Biobank. MCP was positively genetically correlated with AUD (rg = 0.26, p = 7.55 × 10-18), CanUD (rg = 0.37, p = 8.21 × 10-37), OUD (rg = 0.20, p = 1.50 × 10-3), and PTU (rg = 0.29, p = 8.53 × 10-12). Although the MR analyses supported bi-directional relationships, MCP had larger effects on AUD (pain-exposure: beta = 0.18, p = 8.21 × 10-4; pain-outcome: beta = 0.07, p = 0.018), CanUD (pain-exposure: beta = 0.58, p = 2.70 × 10-6; pain-outcome: beta = 0.05, p = 0.014) and PTU (pain-exposure: beta = 0.43, p = 4.16 × 10-8; pain-outcome: beta = 0.09, p = 3.05 × 10-6) than the reverse. The genome-wide analysis identified two SNPs pleiotropic between MCP and all SUD investigated: IHO1 rs7652746 (ppleiotropy = 2.69 × 10-8), and CADM2 rs1248857 (ppleiotropy = 1.98 × 10-5). In the brain-wide analysis, rs7652746 was associated with multiple cerebellum and amygdala imaging phenotypes. When analyzing MCP pleiotropy with each SUD separately, we found 25, 22, and 4 pleiotropic variants for AUD, CanUD, and OUD, respectively. To our knowledge, this is the first large-scale study to provide evidence of potential causal relationships and shared genetic mechanisms underlying MCP-SUD comorbidity.

Modeling the longitudinal changes of ancestry diversity in the Million Veteran Program.

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.

Genetic and non-genetic predictors of risk for opioid dependence.

BACKGROUND: Elucidation of the interaction of biological and psychosocial/environmental factors on opioid dependence (OD) risk can inform our understanding of the etiology of OD. We examined the role of psychosocial/environmental factors in moderating polygenic risk for opioid use disorder (OUD). METHODS: Data from 1958 European ancestry adults who participated in the Yale-Penn 3 study were analyzed. Polygenic risk scores (PRS) were based on a large-scale multi-trait analysis of genome-wide association studies (MTAG) of OUD. RESULTS: A total of 420 (21.1%) individuals had a lifetime diagnosis of OD. OUD PRS were positively associated with OD (odds ratio [OR] 1.42, 95% confidence interval [CI] 1.21-1.66). Household income and education were the strongest correlates of OD. Among individuals with higher OUD PRS, those with higher education level had lower odds of OD (OR 0.92, 95% CI 0.85-0.98); and those with posttraumatic stress disorder (PTSD) were more likely to have OD relative to those without PTSD (OR 1.56, 95% CI 1.04-2.35). CONCLUSIONS: Results suggest an interplay between genetics and psychosocial environment in contributing to OD risk. While PRS alone do not yet have useful clinical predictive utility, psychosocial factors may help enhance prediction. These findings could inform more targeted clinical and policy interventions to help address this public health crisis.

Loneliness and depression: bidirectional mendelian randomization analyses using data from three large genome-wide association studies.

Major depression (MD) is a serious psychiatric illness afflicting nearly 5% of the world's population. A large correlational literature suggests that loneliness is a prospective risk factor for MD; correlational assocations of this nature may be confounded for a variety of reasons. This report uses Mendelian Randomization (MR) to examine potentially causal associations between loneliness and MD. We report on analyses using summary statistics from three large genome wide association studies (GWAS). MR analyses were conducted using three independent sources of GWAS summary statistics. In the first set of analyses, we used available summary statistics from an extant GWAS of loneliness to predict MD risk. We used two sources of outcome data: the Psychiatric Genomics Consortium (PGC) meta-analysis of MD (PGC-MD; N = 142,646) and the Million Veteran Program (MVP-MD; N = 250,215). Finally, we reversed analyses using data from the MVP and PGC samples to identify risk variants for MD and used loneliness outcome data from UK Biobank. We find robust evidence for a bidirectional causal relationship between loneliness and MD, including between loneliness, depression cases status, and a continuous measure of depressive symptoms. The estimates remained significant across several sensitivity analyses, including models that account for horizontal pleiotropy. This paper provides the first genetically-informed evidence that reducing loneliness may play a causal role in decreasing risk for depressive illness, and these findings support efforts to reduce loneliness in order to prevent or ameliorate MD. Discussion focuses on the public health significance of these findings, especially in light of the SARS-CoV-2 pandemic.

Identifying neurofunctional domains across substance use disorders.

Background: Substance use disorders (SUDs) are heterogeneous across multiple functional domains. Various frameworks posit that domains (e.g., executive function) contribute to the persistence of SUDs; however, the domains identified in different studies vary.Objectives: We used factor analysis to identify the underlying latent domains present in a large sample (N = 5,244, 55.8% male) with a variety of SUDs to yield findings more generalizable than studies with a narrower focus.Method: Participants (1,384 controls and 3,860 participants with one or more SUDs including alcohol, cocaine, cannabis, and/or opioid use disorders) completed the Semi-Structured Assessment for Drug Dependence and Alcoholism, the NEO Personality Inventory, and the Wisconsin Card Sorting Test. Exploratory factor analysis (EFA) and fit indices (root mean-squared error of approximation (RMSEA), Comparative Fit Index (CFI), and Tucker-Lewis Index (TLI)) were used to examine different latent variable models. A multiple indicators, multiple causes (MIMIC) approach-tested associations of the latent variables with sociodemographics, substance use, and a history of abuse/neglect.Results: A six-factor model (predominant alcohol, predominant cocaine, predominant opioid, externalizing, personality, and executive function) provided the best fit [RMSEA = 0.063 (90% CI 0.060, 0.066), CFI = 0.98, TLI = 0.96]. All factors were moderately correlated (coefficient = 0.25-0.55, p < .05) with the exception of executive function. MIMIC analysis revealed different patterns of associations (all p < .0001) with sociodemographics, substance use, and a history of abuse/neglect among the factors.Conclusions: The domains identified, particularly executive function, were parallel to those observed previously. These factors underscore the heterogeneous nature of SUDs and may be useful in developing more targeted clinical interventions.

Pleiotropic effects of telomere length loci with brain morphology and brain tissue expression.

Several studies have reported association between leukocyte telomere length (LTL) and neuropsychiatric disorders. Although telomere length is affected by environmental factors, genetic variants in certain loci are strongly associated with LTL. Thus, we aimed to identify the genomic relationship between genetic variants of LTL with brain-based regulatory changes and brain volume. We tested genetic colocalization of seven and nine LTL loci in two ancestry groups, European (EUR) and East-Asian (EAS), respectively, with brain morphology measures for 101 T1-magnetic resonance imaging-based region of interests (n = 21 821). The posterior probability (>90%) was observed for 'fourth ventricle', 'gray matter' and 'cerebellar vermal lobules I-IV' volumes. We then tested causal relationship using LTL loci for gene and methylation expression. We found causal pleiotropy for gene (EAS = four genes; EUR = five genes) and methylation expression (EUR = 17 probes; EAS = 4 probes) of brain tissues (P ≤ 2.47 × 10-6). Integrating chromatin profiles with LTL-single nucleotide polymorphisms identified 45 genes (EUR) and 79 genes (EAS) (P ≤ 9.78×10-7). We found additional 38 LTL-genes using chromatin-based gene mapping for EUR ancestry population. Gene variants in three LTL-genes-GPR37, OBFC1 and RTEL1/RTEL1-TNFRSF6B-show convergent evidence of pleiotropy with brain morphology, gene and methylation expression and chromatin association. Mapping gene functions to drug-gene interactions, we identified process 'transmission across chemical synapses' (P < 2.78 × 10-4). This study provides evidence that genetic variants of LTL have pleiotropic roles with brain-based effects that could explain the phenotypic association of LTL with several neuropsychiatric traits.

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.

Positive personality traits moderate persistent high alcohol consumption, determined by polygenic risk in U.S. military veterans: results from a 10-year, population-based, observational cohort study.

BACKGROUND: Understanding the interplay between psychosocial factors and polygenic risk scores (PRS) may help elucidate the biopsychosocial etiology of high alcohol consumption (HAC). This study examined the psychosocial moderators of HAC, determined by polygenic risk in a 10-year longitudinal study of US military veterans. We hypothesized that positive psychosocial traits (e.g. social support, personality traits, optimism, gratitude) may buffer risk of HAC in veterans with greater polygenic liability for alcohol consumption (AC). METHODS: Data were analyzed from 1323 European-American US veterans who participated in the National Health and Resilience in Veterans Study, a 10-year, nationally representative longitudinal study of US military veterans. PRS reflecting genome-wide risk for AC (AUDIT-C) was derived from a Million Veteran Program genome-wide association study (N = 200 680). RESULTS: Among the total sample, 328 (weighted 24.8%) had persistent HAC, 131 (weighted 9.9%) had new-onset HAC, 44 (weighted 3.3%) had remitted HAC, and 820 (weighted 62.0%) had no/low AC over the 10-year study period. AUDIT-C PRS was positively associated with persistent HAC relative to no/low AC [relative risk ratio (RRR) = 1.43, 95% confidence interval (CI) = 1.23-1.67] and remitted HAC (RRR = 1.63, 95% CI = 1.07-2.50). Among veterans with higher AUDIT-C PRS, greater baseline levels of agreeableness and greater dispositional gratitude were inversely associated with persistent HAC. CONCLUSIONS: AUDIT-C PRS was prospectively associated with persistent HAC over a 10-year period, and agreeableness and dispositional gratitude moderated this association. Clinical interventions designed to target these modifiable psychological traits may help mitigate risk of persistent HAC in veterans with greater polygenic liability for persistent HAC.

Genetic liability to suicidal thoughts and behaviors and risk of suicide attempt in US military veterans: moderating effects of cumulative trauma burden.

BACKGROUND: Little is known about environmental factors that may influence associations between genetic liability to suicidality and suicidal behavior. METHODS: This study examined whether a suicidality polygenic risk score (PRS) derived from a large genome-wide association study (N = 122,935) was associated with suicide attempts in a population-based sample of European-American US military veterans (N = 1664; 92.5% male), and whether cumulative lifetime trauma exposure moderated this association. RESULTS: Eighty-five veterans (weighted 6.3%) reported a history of suicide attempt. After adjusting for sociodemographic and psychiatric characteristics, suicidality PRS was associated with lifetime suicide attempt (odds ratio 2.65; 95% CI 1.37-5.11). A significant suicidality PRS-by-trauma exposure interaction emerged, such that veterans with higher levels of suicidality PRS and greater trauma burden had the highest probability of lifetime suicide attempt (16.6%), whereas the probability of attempts was substantially lower among those with high suicidality PRS and low trauma exposure (1.4%). The PRS-by-trauma interaction effect was enriched for genes implicated in cellular and developmental processes, and nervous system development, with variants annotated to the DAB2 and SPNS2 genes, which are implicated in inflammatory processes. Drug repurposing analyses revealed upregulation of suicide gene-sets in the context of medrysone, a drug targeting chronic inflammation, and clofibrate, a triacylglyceride level lowering agent. CONCLUSION: Results suggest that genetic liability to suicidality is associated with increased risk of suicide attempt among veterans, particularly in the presence of high levels of cumulative trauma exposure. Additional research is warranted to investigate whether incorporation of genomic information may improve suicide prediction models.

Investigation of convergent and divergent genetic influences underlying schizophrenia and alcohol use disorder.

BACKGROUND: Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders. METHODS: We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific. RESULTS: We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001). CONCLUSIONS: Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

Psychosocial moderators of polygenic risk for suicidal ideation: Results from a 7-year population-based, prospective cohort study of U.S. veterans.

Polygenic risk scores (PRS) may help inform the etiology of suicidal thoughts and behaviors. In this study, we evaluated whether a suicidality PRS derived from a large genome-wide association study (GWAS) of suicidality from the UK Biobank (N = 122,935) predicted suicidal ideation (SI) in a 7-year population-based, prospective cohort of European-American US veterans (N = 1326). Results revealed that 8.8% (n = 115) of veterans developed new-onset SI, 4.0% (n = 52) had chronic SI, 3.4% (n = 31) had remitted SI, and 83.8% (n = 1128) denied SI over the study period. Suicidality PRSstandardized was positively associated with chronic SI (relative risk ratio [RRR] = 4.54, 95% confidence interval [CI] = 1.01-20.48) and new-onset SI (RRR = 2.97, 95%CI = 1.22-7.23), and negatively associated with remitted SI (RRR = 0.12, 95% CI = 0.02-0.60). Among veterans with higher suicidality PRS, those with higher baseline dispositional optimism had a lower likelihood of chronic SI (RRR = 0.67, 95% CI = 0.49-0.91) and higher likelihood of remitted SI (RRR = 1.98, 95% CI = 1.18-3.31). Among veterans with higher suicidality PRS, those with higher baseline levels of social support were less likely to develop new-onset SI (RRR = 0.95, 95% CI = 0.92-0.99). These interaction effects were enriched for genes implicated in neuron recognition and development, while the PRS main effect was enriched for genes involved in mannosylation. Collectively, results of this study suggest that suicidality PRS is linked prospectively to symptomatic courses of SI, and that dispositional optimism and social support moderate these associations. Interventions targeting these modifiable psychosocial factors may help mitigate risk of SI in veterans with high polygenic risk for suicidality.

Integrating human brain proteomes with genome-wide association data implicates novel proteins in post-traumatic stress disorder.

Genome-wide association studies (GWAS) have identified several risk loci for post-traumatic stress disorder (PTSD); however, how they confer PTSD risk remains unclear. We aimed to identify genes that confer PTSD risk through their effects on brain protein abundance to provide new insights into PTSD pathogenesis. To that end, we integrated human brain proteomes with PTSD GWAS results to perform a proteome-wide association study (PWAS) of PTSD, followed by Mendelian randomization, using a discovery and confirmatory study design. Brain proteomes (N = 525) were profiled from the dorsolateral prefrontal cortex using mass spectrometry. The Million Veteran Program (MVP) PTSD GWAS (n = 186,689) was used for the discovery PWAS, and the Psychiatric Genomics Consortium PTSD GWAS (n = 174,659) was used for the confirmatory PWAS. To understand whether genes identified at the protein-level were also evident at the transcript-level, we performed a transcriptome-wide association study (TWAS) using human brain transcriptomes (N = 888) and the MVP PTSD GWAS results. We identified 11 genes that contribute to PTSD pathogenesis via their respective cis-regulated brain protein abundance. Seven of 11 genes (64%) replicated in the confirmatory PWAS and 4 of 11 also had their cis-regulated brain mRNA levels associated with PTSD. High confidence level was assigned to 9 of 11 genes after considering evidence from the confirmatory PWAS and TWAS. Most of the identified genes are expressed in other PTSD-relevant brain regions and several are preferentially expressed in excitatory neurons, astrocytes, and oligodendrocyte precursor cells. These genes are novel, promising targets for mechanistic and therapeutic studies to find new treatments for PTSD.

Genome-wide association study in individuals of European and African ancestry and multi-trait analysis of opioid use disorder identifies 19 independent genome-wide significant risk loci.

Despite the large toll of opioid use disorder (OUD), genome-wide association studies (GWAS) of OUD to date have yielded few susceptibility loci. We performed a large-scale GWAS of OUD in individuals of European (EUR) and African (AFR) ancestry, optimizing genetic informativeness by performing MTAG (Multi-trait analysis of GWAS) with genetically correlated substance use disorders (SUDs). Meta-analysis included seven cohorts: the Million Veteran Program, Psychiatric Genomics Consortium, iPSYCH, FinnGen, Partners Biobank, BioVU, and Yale-Penn 3, resulting in a total N = 639,063 (Ncases = 20,686;Neffective = 77,026) across ancestries. OUD cases were defined as having a lifetime OUD diagnosis, and controls as anyone not known to meet OUD criteria. We estimated SNP-heritability (h2SNP) and genetic correlations (rg). Based on genetic correlation, we performed MTAG on OUD, alcohol use disorder (AUD), and cannabis use disorder (CanUD). A leave-one-out polygenic risk score (PRS) analysis was performed to compare OUD and OUD-MTAG PRS as predictors of OUD case status in Yale-Penn 3. The EUR meta-analysis identified three genome-wide significant (GWS; p ≤ 5 × 10-8) lead SNPs-one at FURIN (rs11372849; p = 9.54 × 10-10) and two OPRM1 variants (rs1799971, p = 4.92 × 10-09; rs79704991, p = 1.11 × 10-08; r2 = 0.02). Rs1799971 (p = 4.91 × 10-08) and another OPRM1 variant (rs9478500; p = 1.95 × 10-08; r2 = 0.03) were identified in the cross-ancestry meta-analysis. Estimated h2SNP was 12.75%, with strong rg with CanUD (rg = 0.82; p = 1.14 × 10-47) and AUD (rg = 0.77; p = 6.36 × 10-78). The OUD-MTAG resulted in a GWAS Nequivalent = 128,748 and 18 independent GWS loci, some mapping to genes or gene regions that have previously been associated with psychiatric or addiction phenotypes. The OUD-MTAG PRS accounted for 3.81% of OUD variance (beta = 0.61;s.e. = 0.066; p = 2.00 × 10-16) compared to 2.41% (beta = 0.45; s.e. = 0.058; p = 2.90 × 10-13) explained by the OUD PRS. The current study identified OUD variant associations at OPRM1, single variant associations with FURIN, and 18 GWS associations in the OUD-MTAG. The genetic architecture of OUD is likely influenced by both OUD-specific loci and loci shared across SUDs.

Using phenotype risk scores to enhance gene discovery for generalized anxiety disorder and posttraumatic stress disorder.

UK Biobank (UKB) is a key contributor in mental health genome-wide association studies (GWAS) but only ~31% of participants completed the Mental Health Questionnaire ("MHQ responders"). We predicted generalized anxiety disorder (GAD), posttraumatic stress disorder (PTSD), and major depression symptoms using elastic net regression in the ~69% of UKB participants lacking MHQ data ("MHQ non-responders"; NTraining = 50%; NTest = 50%), maximizing the informative sample for these traits. MHQ responders were more likely to be female, from higher socioeconomic positions, and less anxious than non-responders. Genetic correlation of GAD and PTSD between MHQ responders and non-responders ranged from 0.636 to 1.08; both were predicted by polygenic scores generated from independent cohorts. In meta-analyses of GAD (N = 489,579) and PTSD (N = 497,803), we discovered many novel genomic risk loci (13 for GAD and 40 for PTSD). Transcriptomic analyses converged on altered regulation of prenatal dorsolateral prefrontal cortex in these disorders. Our results provide one roadmap by which sample size and statistical power may be improved for gene discovery of incompletely ascertained traits in the UKB and other biobanks with limited mental health assessment.

Understanding the comorbidity between posttraumatic stress severity and coronary artery disease using genome-wide information and electronic health records.

The association between coronary artery disease (CAD) and posttraumatic stress disorder (PTSD) contributes to the high morbidity and mortality observed for these conditions. To understand the dynamics underlying PTSD-CAD comorbidity, we investigated large-scale genome-wide association (GWA) statistics from the Million Veteran Program (MVP), the UK Biobank (UKB), the Psychiatric Genomics Consortium, and the CARDIoGRAMplusC4D Consortium. We observed a genetic correlation of CAD with PTSD case-control and quantitative outcomes, ranging from 0.18 to 0.32. To investigate possible cause-effect relationships underlying these genetic correlations, we performed a two-sample Mendelian randomization (MR) analysis, observing a significant bidirectional relationship between CAD and PTSD symptom severity. Genetically-determined PCL-17 (PTSD 17-item Checklist) total score was associated with increased CAD risk (odds ratio = 1.04; 95% confidence interval, 95% CI = 1.01-1.06). Conversely, CAD genetic liability was associated with reduced PCL-17 total score (beta = -0.42; 95% CI = -0.04 to -0.81). Because of these opposite-direction associations, we conducted a pleiotropic meta-analysis to investigate loci with concordant vs. discordant effects on PCL-17 and CAD, observing that concordant-effect loci were enriched for molecular pathways related to platelet amyloid precursor protein (beta = 1.53, p = 2.97 × 10-7) and astrocyte activation regulation (beta = 1.51, p = 2.48 × 10-6) while discordant-effect loci were enriched for biological processes related to lipid metabolism (e.g., triglyceride-rich lipoprotein particle clearance, beta = 2.32, p = 1.61 × 10-10). To follow up these results, we leveraged MVP and UKB electronic health records (EHR) to assess longitudinal changes in the association between CAD and posttraumatic stress severity. This EHR-based analysis highlighted that earlier CAD diagnosis is associated with increased PCL-total score later in life, while lower PCL total score was associated with increased risk of a later CAD diagnosis (Mann-Kendall trend test: MVP tau = 0.932, p < 2 × 10-16; UKB tau = 0.376, p = 0.005). In conclusion, both our genetically-informed analyses and our EHR-based follow-up investigation highlighted a bidirectional relationship between PTSD and CAD where multiple pleiotropic mechanisms are likely to be involved.

Genetically regulated multi-omics study for symptom clusters of posttraumatic stress disorder highlights pleiotropy with hematologic and cardio-metabolic traits.

Posttraumatic stress disorder (PTSD) is a psychiatric disorder that may arise in response to severe traumatic event and is diagnosed based on three main symptom clusters (reexperiencing, avoidance, and hyperarousal) per the Diagnostic Manual of Mental Disorders (version DSM-IV-TR). In this study, we characterized the biological heterogeneity of PTSD symptom clusters by performing a multi-omics investigation integrating genetically regulated gene, splicing, and protein expression in dorsolateral prefrontal cortex tissue within a sample of US veterans enrolled in the Million Veteran Program (N total = 186,689). We identified 30 genes in 19 regions across the three PTSD symptom clusters. We found nine genes to have cell-type specific expression, and over-representation of miRNA-families - miR-148, 30, and 8. Gene-drug target prioritization approach highlighted cyclooxygenase and acetylcholine compounds. Next, we tested molecular-profile based phenome-wide impact of identified genes with respect to 1678 phenotypes derived from the Electronic Health Records of the Vanderbilt University biorepository (N = 70,439). Lastly, we tested for local genetic correlation across PTSD symptom clusters which highlighted metabolic (e.g., obesity, diabetes, vascular health) and laboratory traits (e.g., neutrophil, eosinophil, tau protein, creatinine kinase). Overall, this study finds comprehensive genomic evidence including clinical and regulatory profiles between PTSD, hematologic and cardiometabolic traits, that support comorbidities observed in epidemiologic studies of PTSD.

Psychosocial Factors Associated With Accelerated GrimAge in Male U.S. Military Veterans.

OBJECTIVE: Veterans are at high risk for health morbidities linked to premature mortality. Recently developed "epigenetic clock" algorithms, which compute intra-individual differences between biological and chronological aging, can help inform prediction of accelerated biological aging and mortality risk. To date, however, scarce research has examined potentially modifiable correlates of GrimAge, a novel epigenetic clock comprised of DNA methylation surrogates of plasma proteins and smoking pack-years associated with various morbidities and time-to-death. The objective of the study was to examine psychosocial correlates of this novel epigenetic clock. DESIGN: Cross-sectional study. SETTING: U.S. veteran population. PARTICIPANTS: Participants were male, European American (EA), and derived from a nationally representative sample of U.S. veterans (N = 1,135, mean age = 63.3, standard deviation [SD] = 13.0). MEASUREMENTS: We examined the prevalence of accelerated GrimAge and its association with a broad range of health, lifestyle, and psychosocial variables. RESULTS: A total 18.3% of veterans had accelerated GrimAge (≥5 years greater GrimAge than chronological age; mean = 8.4 years acceleration, SD = 2.2). Fewer days of weekly physical exercise (relative variance explained [RVE] = 27%), history of lifetime substance use disorder (RVE = 21%), greater number of lifetime traumas (RVE = 19%), lower gratitude (RVE = 13%), reduced sleep quality (RVE = 7%), lower openness to experience (RVE = 7%), and unmarried/partnered status (RVE = 6%) were independently associated with increased odds of accelerated GrimAge. Increasing numbers of these risk factors were associated with greater odds of accelerated GrimAge, with greatest likelihood of acceleration for veterans with ≥3 risk factors (weighted 21.5%). CONCLUSIONS: These results suggest that nearly 1-of-5 EA male U.S. veterans have accelerated GrimAge, and highlight a broad range of health, lifestyle, and psychosocial variables associated with accelerated GrimAge. Given that many of these factors are modifiable, these findings provide promising leads for risk stratification models of accelerated biological aging and precision medicine-based targets for interventions to mitigate risk for premature mortality in this population.

Dissecting the epigenomic differences between smoking and nicotine dependence in a veteran cohort.

Smoking is a serious public health issue linked to more than 8 million deaths per year worldwide and may lead to nicotine dependence (ND). Although the epigenomic literature on smoking is well established, studies evaluating the role of epigenetics in ND are limited. In this study, we examined the epigenomic signatures of ND and how these differ from smoking exposure to identify biomarkers specific to ND. We investigated the peripheral epigenetic profile of smoking status (SS) and ND in a US male veteran cohort. DNA from saliva was collected from 1135 European American (EA) male US military veterans. DNAm was assessed using the Illumina Infinium Human MethylationEPIC BeadChip array. SS was evaluated as current smokers (n = 137; 12.1%) and non-current smokers (never and former; n = 998; 87.9%). NDFTND was assessed as a continuous variable using the Fagerström Test for ND (FTND; n = 1135; mean = 2.54 ± 2.29). Epigenome-wide association studies (EWAS) and co-methylation analyses were conducted for SS and NDFTND . A total of 450 and 22 genome-wide significant differentially methylated sites (DMS) were associated with SS and NDFTND , respectively (15 overlapped DMS). We identified 97 DMS (43 genes) in SS-EWAS previously reported in the literature, including AHRR and F2RL3 genes (p-value: 1.95 × 10-83 to 4.55 × 10-33 ). NDFTND novel DMS mapped to NEUROG1, ANPEP, and SLC29A1. Co-methylation analysis identified 386 modules (11 SS-related and 19 NDFTND -related). SS-related modules showed enrichment for alcoholism, while NDFTND -related modules were enriched for nicotine addiction. This study confirms previous findings associated with SS and identifies novel and-potentially specific-epigenetic biomarkers of ND that may inform prognosis and novel treatment strategies.