Genetic correlations estimated from GWAS reveal pervasive pleiotropy across a wide variety of phenotypes. We introduce genomic structural equation modeling (Genomic SEM), a multivariate method for analyzing the joint genetic architecture of complex traits. Genomic SEM synthesizes genetic correlations and SNP-heritabilities inferred from GWAS summary statistics of individual traits from samples with varying and unknown degrees of overlap. Genomic SEM can be used to model multivariate genetic associations among phenotypes, identify variants with effects on general dimensions of cross-trait liability, calculate more predictive polygenic scores, and identify loci that cause divergence between traits. We demonstrate several applications of Genomic SEM, including a joint analysis of summary statistics from five psychiatric traits. We identify 27 independent SNPs not previously identified in the contributing univariate GWASs. Polygenic scores from Genomic SEM consistently outperform those from univariate GWAS. Genomic SEM is flexible, open ended, and allows for continuous innovation in multivariate genetic analysis.

Cannabis use is a heritable trait that has been associated with adverse mental health outcomes. In the largest genome-wide association study (GWAS) for lifetime cannabis use to date (N = 184,765), we identified eight genome-wide significant independent single nucleotide polymorphisms in six regions. All measured genetic variants combined explained 11% of the variance. Gene-based tests revealed 35 significant genes in 16 regions, and S-PrediXcan analyses showed that 21 genes had different expression levels for cannabis users versus nonusers. The strongest finding across the different analyses was CADM2, which has been associated with substance use and risk-taking. Significant genetic correlations were found with 14 of 25 tested substance use and mental health-related traits, including smoking, alcohol use, schizophrenia and risk-taking. Mendelian randomization analysis showed evidence for a causal positive influence of schizophrenia risk on cannabis use. Overall, our study provides new insights into the etiology of cannabis use and its relation with mental health.

We introduce two novel methods for multivariate genome-wide association meta-analysis (GWAMA) of related traits that correct for sample overlap. A broad range of simulation scenarios supports the added value of our multivariate methods relative to univariate GWAMA. We applied the novel methods to life satisfaction, positive affect, neuroticism, and depressive symptoms, collectively referred to as the well-being spectrum (N obs = 2,370,390), and found 304 significant independent signals. Our multivariate approaches resulted in a 26% increase in the number of independent signals relative to the four univariate GWAMA, and in a ~ 57% increase in the predictive power of polygenic risk scores. Supporting transcriptome -and methylome-wide analyses (TWAS/MWAS) uncovered an additional 17 and 75 independent loci, respectively. Bioinformatic analyses, based on gene expression in brain tissues and cells, showed that genes differentially expressed in the subiculum and GABAergic interneurons are enriched in their effect on the wellbeing spectrum.In the past decade, genome-wide association studies (GWAS) have provided insights into the genetic basis of quantitative variation in complex traits 1 . With summary statistics of these GWASs becoming public and the development of linkage disequilibrium score regression (LDSC) 2,3 , genetic correlations between traits can be systematically estimated (e.g. Brainstorm consortium 4 ). Levering this widely observed genetic overlap between traits, we introduce two novel methods for multivariate genome-wide association metaanalysis, where we define a multivariate model as a model where the effect of a single SNP is considered for multiple traits: 1) N-weighted multivariate GWAMA (N-GWAMA), with a unitary effect of the SNP on all traits, and 2) model averaging GWAMA (MA-GWAMA), where we relaxed the assumption of a unitary effect of the SNP on all traits. Both methods are well equipped to deal with (unknown) sample overlap. The dependence between effect sizes (error correlation) induced by possible sample overlap is estimated from the univariate GWAMA using LDSC 2,3 . Furthermore, the univariate LDSC intercept is used to correct for population stratification and cryptic relatedness. Both methods have advantages over existing methods. In contrast to MultiPhen 5 , CCA (mv-PLINK) 6 , Combined-PC 7 , and mv-BIMBAM 8 , both our methods can be applied without the need of individual-level genotypic data as only GWAS/GWAMA summary-statistics are required. Additionally, in contrast to S Hom 9 , N -and MA-GWAMA take a more precise estimate of the error correlation into account. In contrast to MTAG 10 , MA-GWAMA , similar to S Het 9 , generates trait specific estimates for each SNP allowing for a certain degree of heterogeneity (see online methods). Finally, in contrast to TATES 11 , both N-GWAMA and MA-GWAMA generate effect sizes for the multivariate effect where TATES only generates a P-value. The absence of a signed statistic in TATES complicates or even prohibits polygenic prediction.

Little is known about the genetic architecture of traits affecting educational attainment other than cognitive ability. We used Genomic Structural Equation Modeling and prior genome-wide association studies (GWAS) of educational attainment ( n = 1,131,881) and cognitive test performance ( n = 257,841) to estimate SNP associations with educational attainment variation that is independent of cognitive ability.We identified 157 genome-wide significant loci and a polygenic architecture accounting for 57% of genetic variance in educational attainment. Non-cognitive genetics were enriched in the same brain tissues and cell types as cognitive performance but showed different associations with gray-matter brain volumes. Non-cognitive genetics were further distinguished by associations with personality traits, less risky behavior,and increased risk for certain psychiatric disorders.For socioeconomic success and longevity, non-cognitive and cognitive-performance genetics demonstrated similar-magnitude associations. By conducting a GWAS of a phenotype that was not directly measured, we offer a first view of genetic architecture of non-cognitive skills influencing educational success.

General rightsCopyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.• Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal ? Take down policyIf you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

Childhood aggression and its resulting consequences inflict a huge burden on affected children, their relatives, teachers, peers and society as a whole. Aggression during childhood rarely occurs in isolation and is correlated with other symptoms of childhood psychopathology. In this paper, we aim to describe and improve the understanding of the co-occurrence of aggression with other forms of childhood psychopathology. We focus on the co-occurrence of aggression and other childhood behavioural and emotional problems, including other externalising problems, attention problems and anxiety–depression. The data were brought together within the EU-ACTION (Aggression in Children: unravelling gene-environment interplay to inform Treatment and InterventiON strategies) project. We analysed the co-occurrence of aggression and other childhood behavioural and emotional problems as a function of the child’s age (ages 3 through 16 years), gender, the person rating the behaviour (father, mother or self) and assessment instrument. The data came from six large population-based European cohort studies from the Netherlands (2x), the UK, Finland and Sweden (2x). Multiple assessment instruments, including the Child Behaviour Checklist (CBCL), the Strengths and Difficulties Questionnaire (SDQ) and Multidimensional Peer Nomination Inventory (MPNI), were used. There was a good representation of boys and girls in each age category, with data for 30,523 3- to 4-year-olds (49.5% boys), 20,958 5- to 6-year-olds (49.6% boys), 18,291 7- to 8-year-olds (49.0% boys), 27,218 9- to 10-year-olds (49.4% boys), 18,543 12- to 13-year-olds (48.9% boys) and 10,088 15- to 16-year-olds (46.6% boys). We replicated the well-established gender differences in average aggression scores at most ages for parental ratings. The gender differences decreased with age and were not present for self-reports. Aggression co-occurred with the majority of other behavioural and social problems, from both externalising and internalising domains. At each age, the co-occurrence was particularly prevalent for aggression and oppositional and ADHD-related problems, with correlations of around 0.5 in general. Aggression also showed substantial associations with anxiety–depression and other internalizing symptoms (correlations around 0.4). Co-occurrence for self-reported problems was somewhat higher than for parental reports, but we found neither rater differences, nor differences across assessment instruments in co-occurrence patterns. There were large similarities in co-occurrence patterns across the different European countries. Finally, co-occurrence was generally stable across age and sex, and if any change was observed, it indicated stronger correlations when children grew older. We present an online tool to visualise these associations as a function of rater, gender, instrument and cohort. In addition, we present a description of the full EU-ACTION projects, its first results and the future perspectives.

Childhood aggressive behavior (AGG) has a substantial heritability of around 50%. Here we present a genome-wide association meta-analysis (GWAMA) of childhood AGG, in which all phenotype measures across childhood ages from multiple assessors were included. We analyzed phenotype assessments for a total of 328 935 observations from 87 485 children aged between 1.5 and 18 years, while accounting for sample overlap. We also meta-analyzed within subsets of the data, i.e., within rater, instrument and age. SNP-heritability for the overall meta-analysis (AGGoverall) was 3.31% (SE = 0.0038). We found no genome-wide significant SNPs for AGGoverall. The gene-based analysis returned three significant genes: ST3GAL3 (P = 1.6E–06), PCDH7 (P = 2.0E–06), and IPO13 (P = 2.5E–06). All three genes have previously been associated with educational traits. Polygenic scores based on our GWAMA significantly predicted aggression in a holdout sample of children (variance explained = 0.44%) and in retrospectively assessed childhood aggression (variance explained = 0.20%). Genetic correlations (rg) among rater-specific assessment of AGG ranged from rg = 0.46 between self- and teacher-assessment to rg = 0.81 between mother- and teacher-assessment. We obtained moderate-to-strong rgs with selected phenotypes from multiple domains, but hardly with any of the classical biomarkers thought to be associated with AGG. Significant genetic correlations were observed with most psychiatric and psychological traits (range $$\left| {r_g} \right|$$ r g : 0.19–1.00), except for obsessive-compulsive disorder. Aggression had a negative genetic correlation (rg = ~−0.5) with cognitive traits and age at first birth. Aggression was strongly genetically correlated with smoking phenotypes (range $$\left| {r_g} \right|$$ r g : 0.46–0.60). The genetic correlations between aggression and psychiatric disorders were weaker for teacher-reported AGG than for mother- and self-reported AGG. The current GWAMA of childhood aggression provides a powerful tool to interrogate the rater-specific genetic etiology of AGG.

Abstract Humans are social animals that experience intense suffering when they perceive a lack of social connection. Modern societies are experiencing an epidemic of loneliness. Although the experience of loneliness is universally human, some people report experiencing greater loneliness than others. Loneliness is more strongly associated with mortality than obesity, emphasizing the need to understand the nature of the relationship between loneliness and health. Although it is intuitive that circumstantial factors such as marital status and age influence loneliness, there is also compelling evidence of a genetic predisposition toward loneliness. To better understand the genetic architecture of loneliness and its relationship with associated outcomes, we extended the genome-wide association study meta-analysis of loneliness to 511 280 subjects, and detect 19 significant genetic variants from 16 loci, including four novel loci, as well as 58 significantly associated genes. We investigated the genetic overlap with a wide range of physical and mental health traits by computing genetic correlations and by building loneliness polygenic scores in an independent sample of 18 498 individuals with EHR data to conduct a PheWAS with. A genetic predisposition toward loneliness was associated with cardiovascular, psychiatric, and metabolic disorders and triglycerides and high-density lipoproteins. Mendelian randomization analyses showed evidence of a causal, increasing, the effect of both BMI and body fat on loneliness. Our results provide a framework for future studies of the genetic basis of loneliness and its relationship to mental and physical health.