Search

AimsAnxiety disorders are prevalent and anxiety symptoms (ANX) co-occur with many psychiatric disorders. We aimed to identify genomic loci associated with ANX, characterize its genetic architecture, and genetic overlap with psychiatric disorders.MethodsWe included a genome-wide association study of ANX (meta-analysis of UK Biobank and Million Veterans Program, n = 301,732), schizophrenia (SCZ), bipolar disorder (BIP), major depression (MD), attention-deficit/hyperactivity disorder (ADHD), and autism spectrum disorder (ASD), and validated the findings in the Norwegian Mother, Father, and Child Cohort (n = 95,841). We employed the bivariate causal mixture model and local analysis of covariant association to characterize the genetic architecture including overlap between the phenotypes. Conditional and conjunctional false discovery rate analyses were performed to boost the identification of loci associated with anxiety and shared with psychiatric disorders.ResultsAnxiety was polygenic with 12.9k genetic variants and overlapped extensively with psychiatric disorders (4.1k-11.4k variants) with predominantly positive genetic correlations between anxiety and psychiatric disorders. We identified 119 novel loci for anxiety by conditioning on the psychiatric disorders, and loci shared between anxiety and MD n=47 , BIP n=33 , SCZ n=71 , ADHD n=20 , and ASD n=5 . Genes annotated to anxiety loci exhibit enrichment for a broader range of biological pathways including cell adhesion and neurofibrillary tangle compared with genes annotated to the shared loci.ConclusionsAnxiety is highly polygenic phenotype with extensive genetic overlap with psychiatric disorders, and we identified novel loci for anxiety implicating new molecular pathways. The shared genetic architecture may underlie the extensive cross-disorder comorbidity of anxiety, and the identified molecular underpinnings may lead to potential drug targets.

PurposeCognitive impairment is prevalent among individuals with epilepsy, and increasing evidence indicates that genetic factors can underlie this relationship. However, the extent to which epilepsy subtypes differ in their genetic relationship with cognitive function, and information about the specific genetic variants involved remain largely unknown.MethodsWe investigated the genetic relationship between epilepsies and general cognitive ability (COG) using complementary statistical tools, including linkage disequilibrium score (LDSC) regression, MiXeR and conjunctional false discovery rate (conjFDR). We analyzed genome-wide association study data on COG (n = 269,867) and common epilepsies (n = 27,559 cases, 42,436 controls), including the broad phenotypes 'all epilepsy', focal epilepsies and genetic generalized epilepsies (GGE), as well as specific subtypes. We functionally annotated the identified loci using several biological resources and validated the results in independent samples.ResultsUsing MiXeR, COG (11.2k variants) was estimated to be almost four times more polygenic than 'all epilepsy', GGE, juvenile myoclonic epilepsy (JME), and childhood absence epilepsy (CAE) (2.5k - 2.9k variants). The other epilepsy phenotypes were insufficiently powered for MiXeR analysis. We quantified extensive genetic overlap between COG and epilepsy types, but with varying negative genetic correlations (-0.23 to -0.04). COG was estimated to share 2.9k variants with both GGE and 'all epilepsy', and 2.3k variants with both JME and CAE. Using conjFDR, we identified 66 distinct loci shared between COG and epilepsies, including novel associations for GGE (27), 'all epilepsy' (5), JME (5) and CAE (5). The implicated genes were significantly expressed in multiple brain regions. The results were validated in independent samples (COG: p = 3.62 × 10-7; 'all epilepsy': p = 2.58 × 10-3).ConclusionOur study further dissects the substantial genetic basis shared between epilepsies and COG and identifies novel shared loci. An improved understanding of the genetic relationship between epilepsies and COG may lead to the development of novel comorbidity-targeted epilepsy treatments.

The basal ganglia are subcortical brain structures involved in motor control, cognition, and emotion regulation. We conducted univariate and multivariate genome-wide association analyses (GWAS) to explore the genetic architecture of basal ganglia volumes using brain scans obtained from 34,794 Europeans with replication in 4,808 white and generalization in 5,220 non-white Europeans. Our multivariate GWAS identified 72 genetic loci associated with basal ganglia volumes with a replication rate of 55.6% at P

Genome-wide association studies (GWAS) implicate broad genomic loci containing clusters of highly correlated genetic variants. Finemapping techniques can select and prioritize variants within each GWAS locus which are more likely to have a functional influence on the trait. Here, we present a novel method, Finemap-MiXeR, for finemapping causal variants from GWAS summary statistics, controlling for correlation among variants due to linkage disequilibrium. Our method is based on a variational Bayesian approach and direct optimization of the Evidence Lower Bound (ELBO) of the likelihood function derived from the MiXeR model. After obtaining the analytical expression for ELBOs gradient, we apply Adaptive Moment Estimation (ADAM) algorithm for optimization, allowing us to obtain the posterior causal probability of each variant. Using these posterior causal probabilities, we validated Finemap-MiXeR across a wide range of scenarios using both synthetic data, and real data on height from the UK Biobank. Comparison of Finemap-MiXeR with two existing methods, FINEMAP and SuSiE RSS, demonstrated similar or improved accuracy. Furthermore, our method is computationally efficient in several aspects. For example, unlike many other methods in the literature, its computational complexity does not increase with the number of true causal variants in a locus and it does not require any matrix inversion operation. The mathematical framework of Finemap-MiXeR is flexible and may also be applied to other problems including cross-trait and cross-ancestry finemapping.

BACKGROUND AND OBJECTIVES: Epilepsies are associated with differences in cortical thickness (TH) and surface area (SA). However, the mechanisms underlying these relationships remain elusive. We investigated the extent to which these phenotypes share genetic influences. METHODS: We analyzed genome-wide association study data on common epilepsies (n = 69,995) and TH and SA (n = 32,877) using Gaussian mixture modeling MiXeR and conjunctional false discovery rate (conjFDR) analysis to quantify their shared genetic architecture and identify overlapping loci. We biologically interrogated the loci using a variety of resources and validated in independent samples. RESULTS: The epilepsies (2.4 k-2.9 k variants) were more polygenic than both SA (1.8 k variants) and TH (1.3 k variants). Despite absent genome-wide genetic correlations, there was a substantial genetic overlap between SA and genetic generalized epilepsy (GGE) (1.1 k), all epilepsies (1.1 k), and juvenile myoclonic epilepsy (JME) (0.7 k), as well as between TH and GGE (0.8 k), all epilepsies (0.7 k), and JME (0.8 k), estimated with MiXeR. Furthermore, conjFDR analysis identified 15 GGE loci jointly associated with SA and 15 with TH, 3 loci shared between SA and childhood absence epilepsy, and 6 loci overlapping between SA and JME. 23 loci were novel for epilepsies and 11 for cortical morphology. We observed a high degree of sign concordance in the independent samples. DISCUSSION: Our findings show extensive genetic overlap between generalized epilepsies and cortical morphology, indicating a complex genetic relationship with mixed-effect directions. The results suggest that shared genetic influences may contribute to cortical abnormalities in epilepsies.

BackgroundFor many brain disorders, a subset of patients jointly exhibit alterations in cortical brain structure and elevated levels of circulating immune markers. This may be driven in part by shared genetic architecture. Therefore, we investigated the phenotypic and genetic associations linking global cortical surface area and thickness with blood immune markers (i.e., white blood cell counts and plasma C-reactive protein levels).MethodsLinear regression was used to assess phenotypic associations in 30,823 UK Biobank participants. Genome-wide and local genetic correlations were assessed using linkage disequilibrium score regression and local analysis of covariance annotation. The number of shared trait-influencing genetic variants was estimated using MiXeR. Shared genetic architecture was assessed using a conjunctional false discovery rate framework, and mapped genes were included in gene-set enrichment analyses.ResultsCortical structure and blood immune markers exhibited predominantly inverse phenotypic associations. There were modest genome-wide genetic correlations, the strongest of which were for C-reactive protein levels (rg_surface_area = -0.13, false discovery rate-corrected p = 4.17 × 10-3; rg_thickness = -0.13, false discovery rate-corrected p = 4.00 × 10-2). Meanwhile, local genetic correlations showed a mosaic of positive and negative associations. White blood cells shared on average 46.24% and 38.64% of trait-influencing genetic variants with surface area and thickness, respectively. Additionally, surface area shared 55 unique loci with the blood immune markers while thickness shared 15. Overall, monocyte count exhibited the largest genetic overlap with cortical brain structure. A series of gene enrichment analyses implicated neuronal-, astrocytic-, and schizophrenia-associated genes.ConclusionsThe findings indicate shared genetic underpinnings for cortical brain structure and blood immune markers, with implications for neurodevelopment and understanding the etiology of brain-related disorders.

Low vitamin D (vitD) levels have been consistently reported in schizophrenia (SCZ) suggesting a role in the etiopathology. However, little is known about the role of underlying shared genetic mechanisms. We applied a conditional/conjunctional false discovery rate approach (FDR) on large, nonoverlapping genome-wide association studies for SCZ (N cases = 53 386, N controls = 77 258) and vitD serum concentration (N = 417 580) to evaluate shared common genetic variants. The identified genomic loci were characterized using functional analyses and biological repositories. We observed cross-trait SNP enrichment in SCZ conditioned on vitD and vice versa, demonstrating shared genetic architecture. Applying the conjunctional FDR approach, we identified 72 loci jointly associated with SCZ and vitD at conjunctional FDR

C-reactive protein (CRP) tends to be elevated in individuals with psychiatric disorders. Recent findings have suggested a protective effect of the genetic liability to elevated CRP on schizophrenia risk and a causative effect on depression despite weak genetic correlations, while causal relationships with bipolar disorder were inconclusive. We investigated the shared genetic underpinnings of psychiatric disorders and variation in CRP levels. Genome-wide association studies for CRP (n = 575,531), bipolar disorder (n = 413,466), depression (n = 480,359), and schizophrenia (n = 130,644) were used in causal mixture models to compare CRP with psychiatric disorders based on polygenicity, discoverability, and genome-wide genetic overlap. The conjunctional false discovery rate method was used to identify specific shared genetic loci. Shared variants were mapped to putative causal genes, which were tested for overrepresentation among gene ontology gene-sets. CRP was six to ten times less polygenic (n = 1400 vs 8600-14,500 variants) and had a discoverability one to two orders of magnitude higher than psychiatric disorders. Most CRP-associated variants were overlapping with psychiatric disorders. We identified 401 genetic loci jointly associated with CRP and psychiatric disorders with mixed effect directions. Gene-set enrichment analyses identified predominantly CNS-related gene sets for CRP and each of depression and schizophrenia, and basic cellular processes for CRP and bipolar disorder. In conclusion, CRP has a markedly different genetic architecture to psychiatric disorders, but the majority of CRP associated variants are also implicated in psychiatric disorders. Shared genetic loci implicated CNS-related processes to a greater extent than immune processes, which may have implications for how we conceptualise causal relationships between CRP and psychiatric disorders.

IntroductionThere is a pressing need for non-invasive, cost-effective tools for early detection of Alzheimer's disease (AD).MethodsUsing data from the Alzheimer's Disease Neuroimaging Initiative (ADNI), Cox proportional models were conducted to develop a multimodal hazard score (MHS) combining age, a polygenic hazard score (PHS), brain atrophy, and memory to predict conversion from mild cognitive impairment (MCI) to dementia. Power calculations estimated required clinical trial sample sizes after hypothetical enrichment using the MHS. Cox regression determined predicted age of onset for AD pathology from the PHS.ResultsThe MHS predicted conversion from MCI to dementia (hazard ratio for 80th versus 20th percentile: 27.03). Models suggest that application of the MHS could reduce clinical trial sample sizes by 67%. The PHS alone predicted age of onset of amyloid and tau.DiscussionThe MHS may improve early detection of AD for use in memory clinics or for clinical trial enrichment.HighlightsA multimodal hazard score (MHS) combined age, genetics, brain atrophy, and memory. The MHS predicted time to conversion from mild cognitive impairment to dementia. MHS reduced hypothetical Alzheimer's disease (AD) clinical trial sample sizes by 67%. A polygenic hazard score predicted age of onset of AD neuropathology.

Autism spectrum disorder (ASD) is a highly heritable condition with a large variation in cognitive function. Here we investigated the shared genetic architecture between cognitive traits (intelligence (INT) and educational attainment (EDU)), and risk loci jointly associated with ASD and the cognitive traits. We analyzed data from genome-wide association studies (GWAS) of INT (n = 269,867), EDU (n = 766,345) and ASD (cases n = 18,381, controls n = 27,969). We used the bivariate causal mixture model (MiXeR) to estimate the total number of shared genetic variants, local analysis of co-variant annotation (LAVA) to estimate local genetic correlations, conditional false discovery rate (cond/conjFDR) to identify specific overlapping loci. The MiXeR analyses showed that 12.7k genetic variants are associated with ASD, of which 12.0k variants are shared with EDU, and 11.1k are shared with INT with both positive and negative relationships within overlapping variants. The majority (59-68%) of estimated shared loci have concordant effect directions, with a positive, albeit modest, genetic correlation between ASD and EDU (rg = 0.21, p = 2e-13) and INT (rg = 0.22, p = 4e-12). We discovered 43 loci jointly associated with ASD and cognitive traits (conjFDR

BACKGROUND: Immune mechanisms are indicated in schizophrenia (SCZ). Recent genome-wide association studies (GWAS) have identified genetic variants associated with SCZ and immune-related phenotypes. Here, we use cutting edge statistical tools to identify shared genetic variants between SCZ and white blood cell (WBC) counts and further understand the role of the immune system in SCZ. STUDY DESIGN: GWAS results from SCZ (patients, n = 53 386; controls, n = 77 258) and WBC counts (n = 56 3085) were analyzed. We applied linkage disequilibrium score regression, the conditional false discovery rate method and the bivariate causal mixture model for analyses of genetic associations and overlap, and 2 sample Mendelian randomization to estimate causal effects. STUDY RESULTS: The polygenicity for SCZ was 7.5 times higher than for WBC count and constituted 32%-59% of WBC count genetic loci. While there was a significant but weak positive genetic correlation between SCZ and lymphocytes (rg = 0.05), the conditional false discovery rate method identified 383 shared genetic loci (53% concordant effect directions), with shared variants encompassing all investigated WBC subtypes: lymphocytes, n = 215 (56% concordant); neutrophils, n = 158 (49% concordant); monocytes, n = 146 (47% concordant); eosinophils, n = 135 (56% concordant); and basophils, n = 64 (53% concordant). A few causal effects were suggested, but consensus was lacking across different Mendelian randomization methods. Functional analyses indicated cellular functioning and regulation of translation as overlapping mechanisms. CONCLUSIONS: Our results suggest that genetic factors involved in WBC counts are associated with the risk of SCZ, indicating a role of immune mechanisms in subgroups of SCZ with potential for stratification of patients for immune targeted treatment.

Personality and cognitive function are heritable mental traits whose genetic foundations may be distributed across interconnected brain functions. Previous studies have typically treated these complex mental traits as distinct constructs. We applied the pleiotropy-informed multivariate omnibus statistical test to genome-wide association studies of 35 measures of neuroticism and cognitive function from the UK Biobank (n = 336,993). We identified 431 significantly associated genetic loci with evidence of abundant shared genetic associations, across personality and cognitive function domains. Functional characterization implicated genes with significant tissue-specific expression in all tested brain tissues and brain-specific gene sets. We conditioned independent genome-wide association studies of the Big 5 personality traits and cognitive function on our multivariate findings, boosting genetic discovery in other personality traits and improving polygenic prediction. These findings advance our understanding of the polygenic architecture of these complex mental traits, indicating a prominence of pleiotropic genetic effects across higher order domains of mental function such as personality and cognitive function.

Anorexia nervosa (AN) is a heritable eating disorder (50-60%) with an array of commonly comorbid psychiatric disorders and related traits. Although significant genetic correlations between AN and psychiatric disorders and related traits have been reported, their shared genetic architecture is largely understudied. We investigated the shared genetic architecture of AN and schizophrenia (SCZ), bipolar disorder (BIP), major depression (MD), mood instability (Mood), neuroticism (NEUR), and intelligence (INT). We applied the conditional false discovery rate (FDR) method to identify novel risk loci for AN, and conjunctional FDR to identify loci shared between AN and related phenotypes, to summarize statistics from relevant genome-wide association studies (GWAS). Individual GWAS samples varied from 72,517 to 420,879 participants. Using conditional FDR we identified 58 novel AN loci. Furthermore, we identified 38 unique loci shared between AN and major psychiatric disorders (SCZ, BIP, and MD) and 45 between AN and psychological traits (Mood, NEUR, and INT). In line with genetic correlations, the majority of shared loci showed concordant effect directions. Functional analyses revealed that the shared loci are involved in 65 unique pathways, several of which overlapped across analyses, including the signal by MST1 pathway involved in Hippo signaling. In conclusion, we demonstrated genetic overlap between AN and major psychiatric disorders and related traits, and identified novel risk loci for AN by leveraging this overlap. Our results indicate that some shared characteristics between AN and related disorders and traits may have genetic underpinnings.

Psychiatric disorders and common epilepsies are heritable disorders with a high comorbidity and overlapping symptoms. However, the causative mechanisms underlying this relationship are poorly understood. Here we aimed to identify overlapping genetic loci between epilepsy and psychiatric disorders to gain a better understanding of their comorbidity and shared clinical features. We analysed genome-wide association study data for all epilepsies (n = 44 889), genetic generalized epilepsy (n = 33 446), focal epilepsy (n = 39 348), schizophrenia (n = 77 096), bipolar disorder (n = 406 405), depression (n = 500 199), attention deficit hyperactivity disorder (n = 53 293) and autism spectrum disorder (n = 46 350). First, we applied the MiXeR tool to estimate the total number of causal variants influencing the disorders. Next, we used the conjunctional false discovery rate statistical framework to improve power to discover shared genomic loci. Additionally, we assessed the validity of the findings in independent cohorts, and functionally characterized the identified loci. The epilepsy phenotypes were considerably less polygenic (1.0 K to 3.4 K causal variants) than the psychiatric disorders (5.6 K to 13.9 K causal variants), with focal epilepsy being the least polygenic (1.0 K variants), and depression having the highest polygenicity (13.9 K variants). We observed cross-trait genetic enrichment between genetic generalized epilepsy and all psychiatric disorders and between all epilepsies and schizophrenia and depression. Using conjunctional false discovery rate analysis, we identified 40 distinct loci jointly associated with epilepsies and psychiatric disorders at conjunctional false discovery rate

BACKGROUND: Irritable bowel syndrome (IBS) often co-occurs with psychiatric and gastrointestinal disorders. A recent genome-wide association study (GWAS) identified several genetic risk variants for IBS. However, most of the heritability remains unidentified, and the genetic overlap with psychiatric and somatic disorders is not quantified beyond genome-wide genetic correlations. Here, we characterize the genetic architecture of IBS, further, investigate its genetic overlap with psychiatric and gastrointestinal phenotypes, and identify novel genomic risk loci. METHODS: Using GWAS summary statistics of IBS (53,400 cases and 433,201 controls), and psychiatric and gastrointestinal phenotypes, we performed bivariate casual mixture model analysis to characterize the genetic architecture and genetic overlap between these phenotypes. We leveraged identified genetic overlap to boost the discovery of genomic loci associated with IBS, and to identify specific shared loci associated with both IBS and psychiatric and gastrointestinal phenotypes, using the conditional/conjunctional false discovery rate (condFDR/conjFDR) framework. We used functional mapping and gene annotation (FUMA) for functional analyses. RESULTS: IBS was highly polygenic with 12k trait-influencing variants. We found extensive polygenic overlap between IBS and psychiatric disorders and to a lesser extent with gastrointestinal diseases. We identified 132 independent IBS-associated loci (condFDR

Treatment resistant schizophrenia (TRS) is characterized by repeated treatment failure with antipsychotics. A recent genome-wide association study (GWAS) of TRS showed a polygenic architecture, but no significant loci were identified. Clozapine is shown to be the superior drug in terms of clinical effect in TRS; at the same time it has a serious side effect profile, including weight gain. Here, we sought to increase power for genetic discovery and improve polygenic prediction of TRS, by leveraging genetic overlap with Body Mass Index (BMI). We analysed GWAS summary statistics for TRS and BMI applying the conditional false discovery rate (cFDR) framework. We observed cross-trait polygenic enrichment for TRS conditioned on associations with BMI. Leveraging this cross-trait enrichment, we identified 2 novel loci for TRS at cFDR

BACKGROUND: The relationship between psychotic disorders and cannabis use is heavily debated. Shared underlying genetic risk is one potential explanation. We investigated the genetic association between psychotic disorders (schizophrenia and bipolar disorder) and cannabis phenotypes (lifetime cannabis use and cannabis use disorder). METHODS: We used genome-wide association summary statistics from individuals with European ancestry from the Psychiatric Genomics Consortium, UK Biobank, and International Cannabis Consortium. We estimated heritability, polygenicity, and discoverability of each phenotype. We performed genome-wide and local genetic correlations. Shared loci were identified and mapped to genes, which were tested for functional enrichment. Shared genetic liabilities to psychotic disorders and cannabis phenotypes were explored using causal analyses and polygenic scores, using the Norwegian Thematically Organized Psychosis cohort. FINDINGS: Psychotic disorders were more heritable than cannabis phenotypes and more polygenic than cannabis use disorder. We observed positive genome-wide genetic correlations between psychotic disorders and cannabis phenotypes (range 0·22-0·35) with a mixture of positive and negative local genetic correlations. Three to 27 shared loci were identified for the psychotic disorder and cannabis phenotype pairs. Enrichment of mapped genes implicated neuronal and olfactory cells as well as drug-gene targets for nicotine, alcohol, and duloxetine. Psychotic disorders showed a causal effect on cannabis phenotypes, and lifetime cannabis use had a causal effect on bipolar disorder. Of 2181 European participants from the Norwegian Thematically Organized Psychosis cohort applied in polygenic risk score analyses, 1060 (48·6%) were females and 1121 (51·4%) were males (mean age 33·1 years [SD 11·8]). 400 participants had bipolar disorder, 697 had schizophrenia, and 1044 were healthy controls. Within this sample, polygenic scores for cannabis phenotypes predicted psychotic disorders independently and improved prediction beyond the polygenic score for the psychotic disorders. INTERPRETATION: A subgroup of individuals might have a high genetic risk of developing a psychotic disorder and using cannabis. This finding supports public health efforts to reduce cannabis use, particularly in individuals at high risk or patients with psychotic disorders. Identified shared loci and their functional implications could facilitate development of novel treatments. FUNDING: US National Institutes of Health, the Research Council Norway, the South-East Regional Health Authority, Stiftelsen Kristian Gerhard Jebsen, EEA-RO-NO-2018-0535, European Unions Horizon 2020 Research and Innovation Programme, the Marie Skłodowska-Curie Actions, and University of Oslo Life Science.

Obesity and associated morbidities, metabolic associated fatty liver disease (MAFLD) included, constitute some of the largest public health threats worldwide. Body composition and related risk factors are known to be heritable and identification of their genetic determinants may aid in the development of better prevention and treatment strategies. Recently, large-scale whole-body MRI data has become available, providing more specific measures of body composition than anthropometrics such as body mass index. Here, we aimed to elucidate the genetic architecture of body composition, by conducting genome-wide association studies (GWAS) of these MRI-derived measures. We ran both univariate and multivariate GWAS on fourteen MRI-derived measurements of adipose and muscle tissue distribution, derived from scans from 33,588 White European UK Biobank participants (mean age of 64.5 years, 51.4% female). Through multivariate analysis, we discovered 100 loci with distributed effects across the body composition measures and 241 significant genes primarily involved in immune system functioning. Liver fat stood out, with a highly discoverable and oligogenic architecture and the strongest genetic associations. Comparison with 21 common cardiometabolic traits revealed both shared and specific genetic influences, with higher mean heritability for the MRI measures (h2 = .25 vs. .13, p = 1.8x10-7). We found substantial genetic correlations between the body composition measures and a range of cardiometabolic diseases, with the strongest correlation between liver fat and type 2 diabetes (rg = .49, p = 2.7x10-22). These findings show that MRI-derived body composition measures complement conventional body anthropometrics and other biomarkers of cardiometabolic health, highlighting the central role of liver fat, and improving our knowledge of the genetic architecture of body composition and related diseases.

OBJECTIVE: Mental disorders are heritable and polygenic, and genome-wide genetic correlations (rg) have indicated widespread shared genetic risk across multiple disorders and related traits, mirroring their overlapping clinical characteristics. However, rg may underestimate the shared genetic underpinnings of mental disorders and related traits because it does not differentiate mixtures of concordant and discordant genetic effects from an absence of genetic overlap. Using novel statistical genetics tools, the authors aimed to evaluate the genetic overlap between mental disorders and related traits when accounting for mixed effect directions. METHODS: The authors applied the bivariate causal mixture model (MiXeR) to summary statistics for four mental disorders, four related mental traits, and height from genome-wide association studies (Ns ranged from 53,293 to 766,345). MiXeR estimated the number of causal variants for a given trait (polygenicity), the number of variants shared between traits, and the genetic correlation of shared variants (rgs). Local rg was investigated using LAVA. RESULTS: Among mental disorders, ADHD was the least polygenic (5.6K causal variants), followed by bipolar disorder (8.6K), schizophrenia (9.6K), and depression (14.5K). Most variants were shared across mental disorders (4.4K-9.3K) and between mental disorders and related traits (5.2K-12.8K), but with disorder-specific variations in rg and rgs. Overlap with height was small (0.7K-1.1K). MiXeR estimates correlated with LAVA local rg (r=0.88, p

Genome-Wide Association studies have typically been limited to univariate analysis in which a single outcome measure is tested against millions of variants. Recent work demonstrates that a Multivariate Omnibus Statistic Test (MOSTest) is well powered to discover genomic effects distributed across multiple phenotypes. Applied to cortical brain MRI morphology measures, MOSTest has resulted in a drastic improvement in power to discover loci when compared to established approaches (min-P). One question that arises is how well these discovered loci replicate in independent data. Here we perform 10 times cross validation within 34,973 individuals from UK Biobank for imaging measures of cortical area, thickness and sulcal depth (>1,000 dimensionality for each). By deploying a replication method that aggregates discovered effects distributed across multiple phenotypes, termed PolyVertex Score (MOSTest-PVS), we demonstrate a higher replication yield and comparable replication rate of discovered loci for MOSTest (# replicated loci: 242-496, replication rate: 96-97%) in independent data when compared with the established min-P approach (# replicated loci: 26-55, replication rate: 91-93%). An out-of-sample replication of discovered loci was conducted with a sample of 4,069 individuals from the Adolescent Brain Cognitive Development® (ABCD) study, who are on average 50 years younger than UK Biobank individuals. We observe a higher replication yield and comparable replication rate of MOSTest-PVS compared to min-P. This finding underscores the importance of using well-powered multivariate techniques for both discovery and replication of high dimensional phenotypes in Genome-Wide Association studies.

Alcohol use disorder (AUD) is a pervasive and devastating mental illness with high comorbidity rates with other mental disorders. Understanding the genetic architecture of this comorbidity could be improved by focusing on intermediate traits that show positive genetic correlation with the disorders. Thus, we aimed to characterize the shared vs. unique polygenicity of AUD, alcohol consumption (AC) and mood instability (MOOD) -beyond genetic correlation, and boost discovery for jointly-associated loci. Summary statistics for MOOD (a binary measure of the tendency to report frequent mood swings), AC (number of standard drinks over a typical consumption week) and AUD GWASs (Ns > 200,000) were analyzed to characterize the cross-phenotype associations between MOOD and AC, MOOD and AUD and AC and AUD. To do so, we used a newly established pipeline that combines (i) the bivariate causal mixture model (MiXeR) to quantify polygenic overlap and (ii) the conjunctional false discovery rate (conjFDR) to discover specific jointly associated genomic loci, which were mapped to genes and biological functions. MOOD was highly polygenic (10.4k single nucleotide polymorphisms, SNPs, SD = 2k) compared to AC (4.9k SNPs, SD = 0.6k) and AUD (4.3k SNPs, SD = 2k). The polygenic overlap of MOOD and AC was twice that of MOOD and AUD (98% vs. 49%), with opposite genetic correlation (-0.2 vs. 0.23), as confirmed in independent samples. MOOD&AUD associated SNPs were significantly enriched for brain genes, conversely to MOOD&AC. Among 38 jointly associated loci, fifteen were novel for MOOD, AC and AUD. MOOD, AC and AUD were also strongly associated at the phenotypic level. Overall, using multilevel polygenic quantification, joint loci discovery and functional annotation methods, we evidenced that the polygenic overlap between MOOD and AC/AUD implicated partly shared biological underpinnings, yet, clearly distinct functional patterns between MOOD&AC and MOOD&AUD, suggesting new mechanisms for the comorbidity of AUD with mood disorders.

Recent genome-wide association studies of mood instability (MOOD) have found significant positive genetic correlation with major depression (DEP) and weak correlations with other psychiatric disorders. We investigated the polygenic overlap between MOOD and psychiatric disorders beyond genetic correlation to better characterize putative shared genetic determinants. GWAS summary statistics for schizophrenia (SCZ, n = 105,318), bipolar disorder (BIP, n = 413,466), DEP (n = 450,619), attention-deficit hyperactivity disorder (ADHD, n = 53,293), and MOOD (n = 363,705) were analyzed using the bivariate causal mixture model and conjunctional false discovery rate methods. MOOD correlated positively with all psychiatric disorders, but with wide variation in strength (rg  = 0.10-0.62). Of 10.4 K genomic variants influencing MOOD, 4 K-9.4 K influenced psychiatric disorders. Furthermore, MOOD was jointly associated with DEP at 163 loci, SCZ at 110, BIP at 60 and ADHD at 25. Fifty-three jointly associated loci were overlapping across two or more disorders, seven of which had discordant effect directions on psychiatric disorders. Genes mapped to loci associated with MOOD and all four disorders were enriched in a single gene-set, synapse organization. The extensive polygenic overlap indicates shared molecular underpinnings across MOOD and psychiatric disorders. However, distinct patterns of genetic correlation and effect directions may relate to differences in the core clinical features of each disorder.

Despite its major role in complex human functions across the lifespan, most notably navigation, learning and memory, much of the genetic architecture of the hippocampal formation is currently unexplored. Here, through multivariate genome-wide association analysis in volumetric data from 35,411 white British individuals, we reveal 177 unique genetic loci with distributed associations across the hippocampal formation. We identify genetic overlap with eight brain disorders with typical onset at different stages of life, where common genes suggest partly age- and disorder-independent mechanisms underlying hippocampal pathology.

Migraine is three times more prevalent in people with bipolar disorder or depression. The relationship between schizophrenia and migraine is less certain although glutamatergic and serotonergic neurotransmission are implicated in both. A shared genetic basis to migraine and mental disorders has been suggested but previous studies have reported weak or non-significant genetic correlations and five shared risk loci. Using the largest samples to date and novel statistical tools, we aimed to determine the extent to which migraines polygenic architecture overlaps with bipolar disorder, depression and schizophrenia beyond genetic correlation, and to identify shared genetic loci. Summary statistics from genome-wide association studies were acquired from large-scale consortia for migraine (n cases = 59 674; n controls = 316 078), bipolar disorder (n cases = 20 352; n controls = 31 358), depression (n cases = 170 756; n controls = 328 443) and schizophrenia (n cases = 40 675, n controls = 64 643). We applied the bivariate causal mixture model to estimate the number of disorder-influencing variants shared between migraine and each mental disorder, and the conditional/conjunctional false discovery rate method to identify shared loci. Loci were functionally characterized to provide biological insights. Univariate MiXeR analysis revealed that migraine was substantially less polygenic (2.8 K disorder-influencing variants) compared to mental disorders (8100-12 300 disorder-influencing variants). Bivariate analysis estimated that 800 (SD = 300), 2100 (SD = 100) and 2300 (SD = 300) variants were shared between bipolar disorder, depression and schizophrenia, respectively. There was also extensive overlap with intelligence (1800, SD = 300) and educational attainment (2100, SD = 300) but not height (1000, SD = 100). We next identified 14 loci jointly associated with migraine and depression and 36 loci jointly associated with migraine and schizophrenia, with evidence of consistent genetic effects in independent samples. No loci were associated with migraine and bipolar disorder. Functional annotation mapped 37 and 298 genes to migraine and each of depression and schizophrenia, respectively, including several novel putative migraine genes such as L3MBTL2, CACNB2 and SLC9B1. Gene-set analysis identified several putative gene sets enriched with mapped genes including transmembrane transport in migraine and schizophrenia. Most migraine-influencing variants were predicted to influence depression and schizophrenia, although a minority of mental disorder-influencing variants were shared with migraine due to the difference in polygenicity. Similar overlap with other brain-related phenotypes suggests this represents a pool of pleiotropic variants that influence vulnerability to diverse brain-related disorders and traits. We also identified specific loci shared between migraine and each of depression and schizophrenia, implicating shared molecular mechanisms and highlighting candidate migraine genes for experimental validation.

BackgroundPolygenic hazard scores (PHS) estimate age-dependent genetic risk of late-onset Alzheimer's disease (AD), but there is limited information about the performance of PHS on real-world data where the population of interest differs from the model development population and part of the model genotypes are missing or need to be imputed.ObjectiveThe aim of this study was to estimate age-dependent risk of late-onset AD using polygenic predictors in Nordic populations.MethodsWe used Desikan PHS model, based on Cox proportional hazards assumption, to obtain age-dependent hazard scores for AD from individual genotypes in the Norwegian DemGene cohort (n = 2,772). We assessed the risk discrimination and calibration of Desikan model and extended it by adding new genotype markers (the Desikan Nordic model). Finally, we evaluated both Desikan and Desikan Nordic models in two independent Danish cohorts: The Copenhagen City Heart Study (CCHS) cohort (n = 7,643) and The Copenhagen General Population Study (CGPS) cohort (n = 10,886).ResultsWe showed a robust prediction efficiency of Desikan model in stratifying AD risk groups in Nordic populations, even when some of the model SNPs were missing or imputed. We attempted to improve Desikan PHS model by adding new SNPs to it, but we still achieved similar risk discrimination and calibration with the extended model.ConclusionPHS modeling has the potential to guide the timing of treatment initiation based on individual risk profiles and can help enrich clinical trials with people at high risk to AD in Nordic populations.

Polygenic risk scores (PRSs) can boost risk prediction in late-onset Alzheimer's disease (LOAD) beyond apolipoprotein E (APOE) but have not been leveraged to identify genetic resilience factors. Here, we sought to identify resilience-conferring common genetic variants in (1) unaffected individuals having high PRSs for LOAD, and (2) unaffected APOE-ε4 carriers also having high PRSs for LOAD. We used genome-wide association study (GWAS) to contrast "resilient" unaffected individuals at the highest genetic risk for LOAD with LOAD cases at comparable risk. From GWAS results, we constructed polygenic resilience scores to aggregate the addictive contributions of risk-orthogonal common variants that promote resilience to LOAD. Replication of resilience scores was undertaken in eight independent studies. We successfully replicated two polygenic resilience scores that reduce genetic risk penetrance for LOAD. We also showed that polygenic resilience scores positively correlate with polygenic risk scores in unaffected individuals, perhaps aiding in staving off disease. Our findings align with the hypothesis that a combination of risk-independent common variants mediates resilience to LOAD by moderating genetic disease risk.

The folding of the human cerebral cortex is a highly genetically regulated process that allows for a much larger surface area to fit into the cranial vault and optimizes functional organization. Sulcal depth is a robust yet understudied measure of localized folding, previously associated with multiple neurodevelopmental disorders. Here, we report the first genome-wide association study of sulcal depth. Through the multivariate omnibus statistical test (MOSTest) applied to vertex-wise measures from 33,748 U.K. Biobank participants (mean age, 64.3 years; 52.0% female), we identified 856 genome-wide significant loci (P < 5 × 10−8). Comparisons with cortical thickness and surface area indicated that sulcal depth has higher locus yield, heritability, and effective sample size. There was a large amount of genetic overlap between these traits, with gene-based analyses indicating strong associations with neurodevelopmental processes. Our findings demonstrate sulcal depth is a promising neuroimaging phenotype that may enhance our understanding of cortical morphology.

Brain morphology has been shown to be highly heritable, yet only a small portion of the heritability is explained by the genetic variants discovered so far. Here we extended the Multivariate Omnibus Statistical Test (MOSTest) and applied it to genome-wide association studies (GWAS) of vertex-wise structural magnetic resonance imaging (MRI) cortical measures from N=35,657 participants in the UK Biobank. We identified 695 loci for cortical surface area and 539 for cortical thickness, in total 780 unique genetic loci associated with cortical morphology robustly replicated in 8,060 children of mixed ethnicity from the Adolescent Brain Cognitive Development (ABCD) Study®. This reflects more than 8-fold increase in genetic discovery at no cost to generalizability compared to the commonly used univariate GWAS methods applied to region of interest (ROI) data. Functional follow up including gene-based analyses implicated 10% of all protein-coding genes and pointed towards pathways involved in neurogenesis and cell differentiation. Power analysis indicated that applying the MOSTest to vertex-wise structural MRI data triples the effective sample size compared to conventional univariate GWAS approaches. The large boost in power obtained with the vertex-wise MOSTest together with pronounced replication rates and highlighted biologically meaningful pathways underscores the advantage of multivariate approaches in the context of highly distributed polygenic architecture of the human brain.

Increased risk-taking is a central component of bipolar disorder (BIP) and is implicated in schizophrenia (SCZ). Risky behaviours, including smoking and alcohol use, are overrepresented in both disorders and associated with poor health outcomes. Positive genetic correlations are reported but an improved understanding of the shared genetic architecture between risk phenotypes and psychiatric disorders may provide insights into underlying neurobiological mechanisms. We aimed to characterise the genetic overlap between risk phenotypes and SCZ, and BIP by estimating the total number of shared variants using the bivariate causal mixture model and identifying shared genomic loci using the conjunctional false discovery rate method. Summary statistics from genome wide association studies of SCZ, BIP, risk-taking and risky behaviours were acquired (n = 82,315-466,751). Genomic loci were functionally annotated using FUMA. Of 8.6-8.7 K variants predicted to influence BIP, 6.6 K and 7.4 K were predicted to influence risk-taking and risky behaviours, respectively. Similarly, of 10.2-10.3 K variants influencing SCZ, 9.6 and 8.8 K were predicted to influence risk-taking and risky behaviours, respectively. We identified 192 loci jointly associated with SCZ and risk phenotypes and 206 associated with BIP and risk phenotypes, of which 68 were common to both risk-taking and risky behaviours and 124 were novel to SCZ or BIP. Functional annotation implicated differential expression in multiple cortical and sub-cortical regions. In conclusion, we report extensive polygenic overlap between risk phenotypes and BIP and SCZ, identify specific loci contributing to this shared risk and highlight biologically plausible mechanisms that may underlie risk-taking in severe psychiatric disorders.

Patients with bipolar disorder (BIP) have a high risk of cardiovascular disease (CVD), despite considerable individual variation. The mechanisms underlying comorbid CVD in BIP remain largely unknown. We investigated polygenic overlap between BIP and CVD phenotypes, including CVD risk factors and coronary artery disease (CAD). We analyzed large genome-wide association studies of BIP (n = 51,710) and CVD phenotypes (n = 159,208-795,640), using bivariate causal mixture model (MiXeR), which estimates the total amount of shared genetic variants, and conjunctional false discovery rate (FDR), which identifies specific overlapping loci. MiXeR revealed polygenic overlap between BIP and body mass index (BMI) (82%), diastolic and systolic blood pressure (20-22%) and CAD (11%) despite insignificant genetic correlations. Using conjunctional FDR

Genome-wide association studies (GWAS) have identified several common genetic variants influencing major depression and general cognitive abilities, but little is known about whether the two share any of their genetic aetiology. Here we investigate shared genomic architectures between major depression (MD) and general intelligence (INT) with the MiXeR statistical tool and their overlapping susceptibility loci with conjunctional false discovery rate (conjFDR), which evaluate the level of overlap in genetic variants and improve the power for gene discovery between two phenotypes. We analysed GWAS data on MD (n = 480,359) and INT (n = 269,867) to characterize polygenic architecture and identify genetic loci shared between these phenotypes. Despite non-significant genetic correlation (rg = -0.0148, P = 0.50), we observed large polygenic overlap and identified 92 loci shared between MD and INT at conjFDR

We propose an extended Gaussian mixture model for the distribution of causal effects of common single nucleotide polymorphisms (SNPs) for human complex phenotypes that depends on linkage disequilibrium (LD) and heterozygosity (H), while also allowing for independent components for small and large effects. Using a precise methodology showing how genome-wide association studies (GWASs) summary statistics (z-scores) arise through LD with underlying causal SNPs, we applied the model to GWAS of multiple human phenotypes. Our findings indicated that causal effects are distributed with dependence on total LD and H, whereby SNPs with lower total LD and H are more likely to be causal with larger effects; this dependence is consistent with models of the influence of negative pressure from natural selection. Compared with the basic Gaussian mixture model it is built on, the extended model-primarily through quantification of selection pressure-reproduces with greater accuracy the empirical distributions of z-scores, thus providing better estimates of genetic quantities, such as polygenicity and heritability, that arise from the distribution of causal effects.

BackgroundParkinson's disease (PD) and schizophrenia (SCZ) are heritable brain disorders that involve dysregulation of the dopaminergic system. Epidemiological studies have reported potential comorbidity between the disorders, and movement disturbances are common in patients with SCZ before treatment with antipsychotic drugs. Despite this, little is known about shared genetic etiology between the disorders.MethodsWe analyzed recent large genome-wide association studies on patients with SCZ (N = 77,096) and PD (N = 417,508) using a conditional/conjunctional false discovery rate (FDR) approach to evaluate overlap in common genetic variants and improve statistical power for genetic discovery. Using a variety of biological resources, we functionally characterized the identified genomic loci.ResultsWe observed genetic enrichment in PD conditional on associations with SCZ and vice versa, indicating polygenic overlap. We then leveraged this cross-trait enrichment using conditional FDR analysis and identified 9 novel PD risk loci and 1 novel SCZ locus at conditional FDR < .01. Furthermore, we identified 9 genomic loci jointly associated with PD and SCZ at conjunctional FDR < .05. There was an even distribution of antagonistic and agonistic effect directions among the shared loci, in line with the insignificant genetic correlation between the disorders. Of 67 genes mapped to the shared loci, 65 are expressed in the human brain and show cell type-specific expression profiles.ConclusionsAltogether, the study increases understanding of the genetic architectures underlying SCZ and PD, indicating that common molecular genetic mechanisms may contribute to overlapping pathophysiological and clinical features between the disorders.

Clinical and epidemiological evidence suggest that loneliness is associated with severe mental disorders (SMDs) and increases the risk of cardiovascular disease (CVD). However, the mechanisms underlying the relationship between loneliness, SMDs, and CVD risk factors remain unknown. Here we explored overlapping genetic architecture and genetic loci shared between SMDs, loneliness, and CVD risk factors. We analyzed large independent genome-wide association study data on schizophrenia (SCZ), bipolar disorder (BD), major depression (MD), loneliness and CVD risk factors using bivariate causal mixture mode (MiXeR), which estimates the total amount of shared variants, and conditional false discovery rate to evaluate overlap in specific loci. We observed substantial genetic overlap between SMDs, loneliness and CVD risk factors, beyond genetic correlation. We identified 149 loci jointly associated with loneliness and SMDs (MD n = 67, SCZ n = 54, and BD n = 28), and 55 distinct loci jointly associated with loneliness and CVD risk factors. A total of 153 novel loneliness loci were found. Most of the shared loci possessed concordant effect directions, suggesting that genetic risk for loneliness may increase the risk of both SMDs and CVD. Functional analyses of the shared loci implicated biological processes related to the brain, metabolic processes, chromatin and immune system. Altogether, the study revealed polygenic overlap between loneliness, SMDs and CVD risk factors, providing new insights into their shared genetic architecture and common genetic mechanisms.

MotivationDetermining the relative contributions of functional genetic categories is fundamental to understanding the genetic etiology of complex human traits and diseases. Here, we present Annotation Informed-MiXeR, a likelihood-based method for estimating the number of variants influencing a phenotype and their effect sizes across different functional annotation categories of the genome using summary statistics from genome-wide association studies.ResultsExtensive simulations demonstrate that the model is valid for a broad range of genetic architectures. The model suggests that complex human phenotypes substantially differ in the number of causal variants, their localization in the genome and their effect sizes. Specifically, the exons of protein-coding genes harbor more than 90% of variants influencing type 2 diabetes and inflammatory bowel disease, making them good candidates for whole-exome studies. In contrast,