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.

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.

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.

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.

Associations between moderate alcohol consumption, brain iron, and cognition in UK Biobank participants: Observational and mendelian randomization analyses.

BACKGROUND: Brain iron deposition has been linked to several neurodegenerative conditions and reported in alcohol dependence. Whether iron accumulation occurs in moderate drinkers is unknown. Our objectives were to investigate evidence in support of causal relationships between alcohol consumption and brain iron levels and to examine whether higher brain iron represents a potential pathway to alcohol-related cognitive deficits. METHODS AND FINDINGS: Observational associations between brain iron markers and alcohol consumption (n = 20,729 UK Biobank participants) were compared with associations with genetically predicted alcohol intake and alcohol use disorder from 2-sample mendelian randomization (MR). Alcohol intake was self-reported via a touchscreen questionnaire at baseline (2006 to 2010). Participants with complete data were included. Multiorgan susceptibility-weighted magnetic resonance imaging (9.60 ± 1.10 years after baseline) was used to ascertain iron content of each brain region (quantitative susceptibility mapping (QSM) and T2*) and liver tissues (T2*), a marker of systemic iron. Main outcomes were susceptibility (χ) and T2*, measures used as indices of iron deposition. Brain regions of interest included putamen, caudate, hippocampi, thalami, and substantia nigra. Potential pathways to alcohol-related iron brain accumulation through elevated systemic iron stores (liver) were explored in causal mediation analysis. Cognition was assessed at the scan and in online follow-up (5.82 ± 0.86 years after baseline). Executive function was assessed with the trail-making test, fluid intelligence with puzzle tasks, and reaction time by a task based on the "Snap" card game. Mean age was 54.8 ± 7.4 years and 48.6% were female. Weekly alcohol consumption was 17.7 ± 15.9 units and never drinkers comprised 2.7% of the sample. Alcohol consumption was associated with markers of higher iron (χ) in putamen (ß = 0.08 standard deviation (SD) [95% confidence interval (CI) 0.06 to 0.09], p < 0.001), caudate (ß = 0.05 [0.04 to 0.07], p < 0.001), and substantia nigra (ß = 0.03 [0.02 to 0.05], p < 0.001) and lower iron in the thalami (ß = -0.06 [-0.07 to -0.04], p < 0.001). Quintile-based analyses found these associations in those consuming >7 units (56 g) alcohol weekly. MR analyses provided weak evidence these relationships are causal. Genetically predicted alcoholic drinks weekly positively associated with putamen and hippocampus susceptibility; however, these associations did not survive multiple testing corrections. Weak evidence for a causal relationship between genetically predicted alcohol use disorder and higher putamen susceptibility was observed; however, this was not robust to multiple comparisons correction. Genetically predicted alcohol use disorder was associated with serum iron and transferrin saturation. Elevated liver iron was observed at just >11 units (88 g) alcohol weekly c.f. <7 units (56 g). Systemic iron levels partially mediated associations of alcohol intake with brain iron. Markers of higher basal ganglia iron associated with slower executive function, lower fluid intelligence, and slower reaction times. The main limitations of the study include that χ and T2* can reflect changes in myelin as well as iron, alcohol use was self-reported, and MR estimates can be influenced by genetic pleiotropy. CONCLUSIONS: To the best of our knowledge, this study represents the largest investigation of moderate alcohol consumption and iron homeostasis to date. Alcohol consumption above 7 units weekly associated with higher brain iron. Iron accumulation represents a potential mechanism for alcohol-related cognitive decline.

Genetic contributions to anxiety disorders: where we are and where we are heading.

Anxiety disorders are among the most common psychiatric disorders worldwide. They often onset early in life, with symptoms and consequences that can persist for decades. This makes anxiety disorders some of the most debilitating and costly disorders of our time. Although much is known about the synaptic and circuit mechanisms of fear and anxiety, research on the underlying genetics has lagged behind that of other psychiatric disorders. However, alongside the formation of the Psychiatric Genomic Consortium Anxiety workgroup, progress is rapidly advancing, offering opportunities for future research.Here we review current knowledge about the genetics of anxiety across the lifespan from genetically informative designs (i.e. twin studies and molecular genetics). We include studies of specific anxiety disorders (e.g. panic disorder, generalised anxiety disorder) as well as those using dimensional measures of trait anxiety. We particularly address findings from large-scale genome-wide association studies and show how such discoveries may provide opportunities for translation into improved or new therapeutics for affected individuals. Finally, we describe how discoveries in anxiety genetics open the door to numerous new research possibilities, such as the investigation of specific gene-environment interactions and the disentangling of causal associations with related traits and disorders.We discuss how the field of anxiety genetics is expected to move forward. In addition to the obvious need for larger sample sizes in genome-wide studies, we highlight the need for studies among young people, focusing on specific underlying dimensional traits or components of anxiety.

Dissecting the heterogeneity of posttraumatic stress disorder: differences in polygenic risk, stress exposures, and course of PTSD subtypes.

BACKGROUND: Definition of disorder subtypes may facilitate precision treatment for posttraumatic stress disorder (PTSD). We aimed to identify PTSD subtypes and evaluate their associations with genetic risk factors, types of stress exposures, comorbidity, and course of PTSD. METHODS: Data came from a prospective study of three U.S. Army Brigade Combat Teams that deployed to Afghanistan in 2012. Soldiers with probable PTSD (PTSD Checklist for Diagnostic and Statistical Manual of Mental Disorders-Fifth Edition ≥31) at three months postdeployment comprised the sample (N = 423) for latent profile analysis using Gaussian mixture modeling and PTSD symptom ratings as indicators. PTSD profiles were compared on polygenic risk scores (derived from external genomewide association study summary statistics), experiences during deployment, comorbidity at three months postdeployment, and persistence of PTSD at nine months postdeployment. RESULTS: Latent profile analysis revealed profiles characterized by prominent intrusions, avoidance, and hyperarousal (threat-reactivity profile; n = 129), anhedonia and negative affect (dysphoric profile; n = 195), and high levels of all PTSD symptoms (high-symptom profile; n = 99). The threat-reactivity profile had the most combat exposure and the least comorbidity. The dysphoric profile had the highest polygenic risk for major depression, and more personal life stress and co-occurring major depression than the threat-reactivity profile. The high-symptom profile had the highest rates of concurrent mental disorders and persistence of PTSD. CONCLUSIONS: Genetic and trauma-related factors likely contribute to PTSD heterogeneity, which can be parsed into subtypes that differ in symptom expression, comorbidity, and course. Future studies should evaluate whether PTSD typology modifies treatment response and should clarify distinctions between the dysphoric profile and depressive disorders.

Apolipoprotein E gene polymorphism, posttraumatic stress disorder, and cognitive function in older U.S. veterans: Results from the National Health and Resilience in Veterans Study.

BACKGROUND: Although the ε4 allele of the apolipoprotein E (APOE) gene and posttraumatic stress disorder (PTSD) have been linked to cognitive dysfunction and dementia risk, it is unknown whether they interact to predict cognitive dysfunction. METHODS: We analyzed data from European-American (EA) veterans who participated in the National Health and Resilience in Veterans Study (NHRVS): main sample (n = 1,386) and primary replication sample (n = 509). EAs from the Yale-Penn Study cohort (n = 948) served as a second replication sample. Multivariable analyses were conducted to evaluate the predictive effects of ε4 carrier status and PTSD on cognitive functioning, with a focus on whether PTSD moderates the effect of ε4 carrier status. RESULTS: APOE ε4 allele carrier status (d = 0.15 and 0.17 in the main and primary replication NHRVS samples, respectively) and PTSD (d = 0.31 and 0.17, respectively) were independently associated with lower cognitive functioning. ε4 carriers with PTSD scored lower than those without PTSD (d = 0.68 and 1.29, respectively) with the most pronounced differences in executive function (d's = 0.75-1.50) and attention/concentration (d's = 0.62-1.33). A significant interaction was also observed in the Yale-Penn sample, with ε4 carriers with PTSD making more perseverative errors on a measure of executive function than those without PTSD (24.7% vs. 17.6%; d = 0.59). CONCLUSIONS: APOE ε4 allele carriers with PTSD have substantially greater cognitive difficulties than ε4 carriers without PTSD. These results underscore the importance of assessing, monitoring, and treating PTSD in trauma-affected individuals who are at genetic risk for cognitive decline and dementia.

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

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

A Multivariate Genome-Wide Association Study Reveals Neural Correlates and Common Biological Mechanisms of Psychopathology Spectra.

There is considerable comorbidity across externalizing and internalizing behavior dimensions of psychopathology. We applied genomic structural equation modeling (gSEM) to genome-wide association study (GWAS) summary statistics to evaluate the factor structure of externalizing and internalizing psychopathology across 16 traits and disorders among European-ancestry individuals (n's = 16,400 to 1,074,629). We conducted GWAS on factors derived from well-fitting models. Downstream analyses served to identify biological mechanisms, explore drug repurposing targets, estimate genetic overlap between the externalizing and internalizing spectra, and evaluate causal effects of psychopathology liability on physical health. Both a correlated factors model, comprising two factors of externalizing and internalizing risk, and a higher-order single-factor model of genetic effects contributing to both spectra demonstrated acceptable fit. GWAS identified 409 lead single nucleotide polymorphisms (SNPs) associated with externalizing and 85 lead SNPs associated with internalizing, while the second-order GWAS identified 256 lead SNPs contributing to broad psychopathology risk. In bivariate causal mixture models, nearly all externalizing and internalizing causal variants overlapped, despite a genetic correlation of only 0.37 (SE = 0.02) between them. Externalizing genes showed cell-type specific expression in GABAergic, cortical, and hippocampal neurons, and internalizing genes were associated with reduced subcallosal cortical volume, providing insight into the neurobiological underpinnings of psychopathology. Genetic liability for externalizing, internalizing, and broad psychopathology exerted causal effects on pain, general health, cardiovascular diseases, and chronic illnesses. These findings underscore the complex genetic architecture of psychopathology, identify potential biological pathways for the externalizing and internalizing spectra, and highlight the physical health burden of psychiatric comorbidity.

A genome-wide investigation into the underlying genetic architecture of personality traits and overlap with psychopathology.

Personality is influenced by both genetic and environmental factors and is associated with other psychiatric traits such as anxiety and depression. The "Big Five" personality traits, which include neuroticism, extraversion, agreeableness, conscientiousness, and openness, are a widely accepted and influential framework for understanding and describing human personality. Of the big five personality traits, neuroticism has most often been the focus of genetic studies and is linked to various mental illnesses including depression, anxiety, and schizophrenia. Our knowledge of the genetic architecture of the other four personality traits is more limited. Utilizing the Million Veteran Program (MVP) cohort we conducted a genome-wide association study (GWAS) in individuals of European and African ancestry. Adding other published data, we performed GWAS meta-analysis for each of the five personality traits with sample sizes ranging from 237,390 to 682,688. We identified 158, 14, 3, 2, and 7 independent genome-wide significant (GWS) loci associated with neuroticism, extraversion, agreeableness, conscientiousness, and openness, respectively. These findings represent 55 novel loci for neuroticism, as well as the first GWS loci discovered for extraversion and agreeableness. Gene-based association testing revealed 254 genes showing significant association with at least one of the five personality traits. Transcriptome-wide and proteome-wide analysis identified altered expression of genes and proteins such as CRHR1, SLC12A5, MAPT, and STX4. Pathway enrichment and drug perturbation analyses identified complex biology underlying human personality traits. We also studied the inter-relationship of personality traits with 1,437 other traits in a phenome-wide genetic correlation analysis, identifying new associations. Mendelian randomization showed positive bidirectional effects between neuroticism and depression and anxiety while a negative bidirectional effect was observed for agreeableness and these psychiatric traits. This study improves our comprehensive understanding of the genetic architecture underlying personality traits and their relationship to other complex human traits.

Cross-ancestry genetic investigation of schizophrenia, cannabis use disorder, and tobacco smoking.

Individuals with schizophrenia frequently experience co-occurring substance use, including tobacco smoking and heavy cannabis use, and substance use disorders. There is interest in understanding the extent to which these relationships are causal, and to what extent shared genetic factors play a role. We explored the relationships between schizophrenia (Scz), cannabis use disorder (CanUD), and ever-regular tobacco smoking (Smk) using the largest available genome-wide studies of these phenotypes in individuals of African and European ancestries. All three phenotypes were positively genetically correlated (rgs = 0.17 - 0.62). Causal inference analyses suggested the presence of horizontal pleiotropy, but evidence for bidirectional causal relationships was also found between all three phenotypes even after correcting for horizontal pleiotropy. We identified 439 pleiotropic loci in the European ancestry data, 150 of which were novel (i.e., not genome-wide significant in the original studies). Of these pleiotropic loci, 202 had lead variants which showed convergent effects (i.e., same direction of effect) on Scz, CanUD, and Smk. Genetic variants convergent across all three phenotypes showed strong genetic correlations with risk-taking, executive function, and several mental health conditions. Our results suggest that both horizontal pleiotropy and causal mechanisms may play a role in the relationship between CanUD, Smk, and Scz, but longitudinal, prospective studies are needed to confirm a causal relationship.

Whole-exome sequencing in UK Biobank reveals rare genetic architecture for depression.

Nearly two hundred common-variant depression risk loci have been identified by genome-wide association studies (GWAS). However, the impact of rare coding variants on depression remains poorly understood. Here, we present whole-exome sequencing analyses of depression with seven different definitions based on survey, questionnaire, and electronic health records in 320,356 UK Biobank participants. We showed that the burden of rare damaging coding variants in loss-of-function intolerant genes is significantly associated with risk of depression with various definitions. We compared the rare and common genetic architecture across depression definitions by genetic correlation and showed different genetic relationships between definitions across common and rare variants. In addition, we demonstrated that the effects of rare damaging coding variant burden and polygenic risk score on depression risk are additive. The gene set burden analyses revealed overlapping rare genetic variant components with developmental disorder, autism, and schizophrenia. Our study provides insights into the contribution of rare coding variants, separately and in conjunction with common variants, on depression with various definitions and their genetic relationships with neurodevelopmental disorders.