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. Here, utilizing the Million Veteran Program cohort, we conducted a genome-wide association study in individuals of European and African ancestry. Adding other published data, we performed genome-wide association study meta-analysis for each of the five personality traits with sample sizes ranging from 237,390 to 682,688. We identified 208, 14, 3, 2 and 7 independent genome-wide significant loci associated with neuroticism, extraversion, agreeableness, conscientiousness and openness, respectively. These findings represent 62 novel loci for neuroticism, as well as the first genome-wide significant loci discovered for 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.

Genome-wide association study of major anxiety disorders in 122,341 European-ancestry cases identifies 58 loci and highlights GABAergic signaling.

The major anxiety disorders (ANX; including generalized anxiety disorder, panic disorder, and phobias) are highly prevalent, often onset early, persist throughout life, and cause substantial global disability. Although distinct in their clinical presentations, they likely represent differential expressions of a dysregulated threat-response system. Here we present a genome-wide association meta-analysis comprising 122,341 European ancestry ANX cases and 729,881 controls. We identified 58 independent genome-wide significant ANX risk variants and 66 genes with robust biological support. In an independent sample of 1,175,012 self-report ANX cases and 1,956,379 controls, 51 of the 58 associated variants were replicated. As predicted by twin studies, we found substantial genetic correlation between ANX and depression, neuroticism, and other internalizing phenotypes. Follow-up analyses demonstrated enrichment in all major brain regions and highlighted GABAergic signaling as one potential mechanism underlying ANX genetic risk. These results advance our understanding of the genetic architecture of ANX and prioritize genes for functional follow-up studies.

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; European ancestry N = 161,405; African ancestry N = 15,846), cannabis use disorder (CanUD; European ancestry N = 886,025; African ancestry N = 120,208), and ever-regular tobacco smoking (Smk; European ancestry N = 805,431; African ancestry N = 24,278) 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). Genetic instrumental variable analyses suggested the presence of shared heritable factors, but evidence for bidirectional causal relationships was also found between all three phenotypes even after correcting for these shared genetic factors. We identified 327 pleiotropic loci with 439 lead SNPs 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 shared genetic factors 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.

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.

Polygenic liability for anxiety in association with comorbid anxiety in multiple sclerosis.

OBJECTIVE: Comorbid anxiety occurs often in MS and is associated with disability progression. Polygenic scores offer a possible means of anxiety risk prediction but often have not been validated outside the original discovery population. We aimed to investigate the association between the Generalized Anxiety Disorder 2-item scale polygenic score with anxiety in MS. METHODS: Using a case-control design, participants from Canadian, UK Biobank, and United States cohorts were grouped into cases (MS/comorbid anxiety) or controls (MS/no anxiety, anxiety/no immune disease or healthy). We used multiple anxiety measures: current symptoms, lifetime interview-diagnosed, and lifetime self-report physician-diagnosed. The polygenic score was computed for current anxiety symptoms using summary statistics from a previous genome-wide association study and was tested using regression. RESULTS: A total of 71,343 individuals of European genetic ancestry were used: Canada (n = 334; 212 MS), UK Biobank (n = 70,431; 1,390 MS), and the USA (n = 578 MS). Meta-analyses identified that in MS, each 1-SD increase in the polygenic score was associated with ~50% increased odds of comorbid moderate anxious symptoms compared to those with less than moderate anxious symptoms (OR: 1.47, 95% CI: 1.09-1.99). We found a similar direction of effects in the other measures. MS had a similar anxiety genetic burden compared to people with anxiety as the index disease. INTERPRETATION: Higher genetic burden for anxiety was associated with significantly increased odds of moderate anxious symptoms in MS of European genetic ancestry which did not differ from those with anxiety and no comorbid immune disease. This study suggests a genetic basis for anxiety in MS.

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.

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.

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.

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.

Genetic contribution to the comorbidity between attention-deficit/hyperactivity disorder and substance use disorders.

We characterized the genetic architecture of the attention-deficit hyperactivity disorder-substance use disorder (ADHD-SUD) relationship by investigating genetic correlation, causality, pleiotropy, and common polygenic risk. Summary statistics from genome-wide association studies (GWAS) were used to investigate ADHD (Neff = 51,568), cannabis use disorder (CanUD, Neff = 161,053), opioid use disorder (OUD, Neff = 57,120), problematic alcohol use (PAU, Neff = 502,272), and problematic tobacco use (PTU, Neff = 97,836). ADHD, CanUD, and OUD GWAS meta-analyses included cohorts with case definitions based on different diagnostic criteria. PAU GWAS combined information related to alcohol use disorder, alcohol dependence, and the items related to alcohol problematic consequences assessed by the alcohol use disorders identification test. PTU GWAS was generated a multi-trait analysis including information regarding Fagerström Test for Nicotine Dependence and cigarettes per day. Linkage disequilibrium score regression analyses indicated positive genetic correlation with CanUD, OUD, PAU, and PTU. Genomic structural equation modeling showed that these genetic correlations were related to two latent factors: one including ADHD, CanUD, and PTU and the other with OUD and PAU. The evidence of a causal effect of PAU and PTU on ADHD was stronger than the reverse in the two-sample Mendelian randomization analysis. Conversely, similar strength of evidence was found between ADHD and CanUD. CADM2 rs62250713 was a pleiotropic SNP between ADHD and all SUDs. We found seven, one, and twenty-eight pleiotropic variants between ADHD and CanUD, PAU, and PTU, respectively. Finally, OUD, CanUD, and PAU PRS were associated with increased odds of ADHD. Our findings demonstrated the contribution of multiple pleiotropic mechanisms to the comorbidity between ADHD and SUDs.

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.

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.

Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference.

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals.

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Multi-ancestry genome-wide association study of cannabis use disorder yields insight into disease biology and public health implications.

As recreational use of cannabis is being decriminalized in many places and medical use widely sanctioned, there are growing concerns about increases in cannabis use disorder (CanUD), which is associated with numerous medical comorbidities. Here we performed a genome-wide association study of CanUD in the Million Veteran Program (MVP), followed by meta-analysis in 1,054,365 individuals (ncases = 64,314) from four broad ancestries designated by the reference panel used for assignment (European n = 886,025, African n = 123,208, admixed American n = 38,289 and East Asian n = 6,843). Population-specific methods were applied to calculate single nucleotide polymorphism-based heritability within each ancestry. Statistically significant single nucleotide polymorphism-based heritability for CanUD was observed in all but the smallest population (East Asian). We discovered genome-wide significant loci unique to each ancestry: 22 in European, 2 each in African and East Asian, and 1 in admixed American ancestries. A genetically informed causal relationship analysis indicated a possible effect of genetic liability for CanUD on lung cancer risk, suggesting potential unanticipated future medical and psychiatric public health consequences that require further study to disentangle from other known risk factors such as cigarette smoking.

Genetic Decomposition of the Heritable Component of Reported Childhood Maltreatment.

Background: Decades of research have shown that environmental exposures, including self-reports of trauma, are partly heritable. Heritable characteristics may influence exposure to and interpretations of environmental factors. Identifying heritable factors associated with self-reported trauma could improve our understanding of vulnerability to exposure and the interpretation of life events. Methods: We used genome-wide association study summary statistics of childhood maltreatment, defined as reporting of abuse (emotional, sexual, and physical) and neglect (emotional and physical) (N = 185,414 participants). We calculated genetic correlations (rg) between reported childhood maltreatment and 576 traits to identify phenotypes that might explain the heritability of reported childhood maltreatment, retaining those with |rg| > 0.25. We specified multiple regression models using genomic structural equation modeling to detect residual genetic variance in childhood maltreatment after accounting for genetically correlated traits. Results: In 2 separate models, the shared genetic component of 12 health and behavioral traits and 7 psychiatric disorders accounted for 59% and 56% of heritability due to common genetic variants (single nucleotide polymorphism-based heritability [h2SNP]) of childhood maltreatment, respectively. Genetic influences on h2SNP of childhood maltreatment were generally accounted for by a shared genetic component across traits. The exceptions to this were general risk tolerance, subjective well-being, posttraumatic stress disorder, and autism spectrum disorder, identified as independent contributors to h2SNP of childhood maltreatment. These 4 traits alone were sufficient to explain 58% of h2SNP of childhood maltreatment. Conclusions: We identified putative traits that reflect h2SNP of childhood maltreatment. Elucidating the mechanisms underlying these associations may improve trauma prevention and posttraumatic intervention strategies.

Genome-wide association studies and cross-population meta-analyses investigating short and long sleep duration.

Sleep duration has been linked to a wide range of negative health outcomes and to reduced life expectancy. We present genome-wide association studies of short ( ≤ 5 h) and long ( ≥ 10 h) sleep duration in adults of European (N = 445,966), African (N = 27,785), East Asian (N = 3141), and admixed-American (N = 16,250) ancestry from UK Biobank and the Million Veteran Programme. In a cross-population meta-analysis, we identify 84 independent loci for short sleep and 1 for long sleep. We estimate SNP-based heritability for both sleep traits in each ancestry based on population derived linkage disequilibrium (LD) scores using cov-LDSC. We identify positive genetic correlation between short and long sleep traits (rg = 0.16 ± 0.04; p = 0.0002), as well as similar patterns of genetic correlation with other psychiatric and cardiometabolic phenotypes. Mendelian randomisation reveals a directional causal relationship between short sleep and depression, and a bidirectional causal relationship between long sleep and depression.

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.