A meta-analysis of the GWAS results of the 28 discovery samples was performed through fixed-effects meta-analysis in METAL, using SNP p values weighted by effective sample size. We meta-analyzed the BroadABC data (N = 50,252) with case-control data of the Psychiatric Genetics Consortium/iPSYCH consortium meta-analysis (N cases = 3802, N controls = 31,305, N effective = 6323) in METAL, leading to a total effective sample size of 56,575. Since some individuals participating in the eight QIMR studies may be related to some degree, we first meta-analyzed those samples with the sample overlap option in METAL before meta-analyzing those results with the rest of the samples. Although throughout the paper we report the total sample size (N = 85,359), we used the effective sample size (the METAL output) to calculate the PRS, SNP-based heritability (h2SNP) and genetic correlations. After combining all cleaned GWAS data files, meta-analysis results were filtered to exclude any variants with N < 30,000. Consequently, we removed 2,134,049 SNPs, resulting in 7,392,849 SNPs available for analysis. To investigate sex-specific genetic effects, we also ran the meta-analysis in the
