Unraveling shared susceptibility loci and Mendelian genetic associations linking educational attainment with multiple neuropsychiatric disorders

Background Empirical studies have demonstrated that educational attainment (EA) is associated with neuropsychiatric disorders (NPDs), suggesting a shared etiological basis between them. However, little is known about the shared genetic mechanisms and causality behind such associations. Methods This study explored the shared genetic basis and causal relationships between EA and NPDs using the high-definition likelihood (HDL) method, cross phenotype association study (CPASSOC), transcriptome-wide association study (TWAS), and bidirectional Mendelian randomization (MR) with summary-level data for EA (N = 293,723) and NPDs (N range = 9,725 to 455,258). Results Significant genetic correlations between EA and 12 NPDs (rg range − 0.49 to 0.35; all p < 3.85 × 10−3) were observed. CPASSOC identified 37 independent loci shared between EA and NPDs, one of which was novel (rs71351952, mapped gene: ARFGEF2). Functional analyses and TWAS found shared genes were enriched in brain tissue, especially in the cerebellum and highlighted the regulatory role of neuronal signaling, purine nucleotide metabolic process, and cAMP-mediated signaling pathways. CPASSOC and TWAS supported the role of three regions of 6q16.1, 3p21.31, and 17q21.31 might account for the shared causes between EA and NPDs. MR confirmed higher genetically predicted EA lower the risk of ADHD (ORIVW: 0.50; 95% CI: 0.39 to 0.63) and genetically predicted ADHD decreased the risk of EA (Causal effect: −2.8 months; 95% CI: −3.9 to −1.8). Conclusion These findings provided evidence of shared genetics and causation between EA and NPDs, advanced our understanding of EA, and implicated potential biological pathways that might underlie both EA and NPDs.


Highlights
• The strongest negative genetic correlation was observed between educational attainment and attention deficit/hyperactivity disorder; the strongest positive genetic correlation was observed between educational attainment and bipolar disorder.• Cross-trait meta-analysis identifies 37 independent genomic loci shared between educational attainment and neuropsychiatric disorders, one of which is novel.• Associated genetic variants may indicate potential therapeutic targets and implicated signaling pathways in low educational attainment phenotypes and neuropsychiatric disorders.• Mendelian randomization analysis provided strong evidence in support of low educational attainment causally increasing the risk of attention deficit/hyperactivity disorder, Alzheimer's dementia, and cannabis use disorder.• Cross-trait meta-analysis and transcriptome-wide association studies support the role of three particular region 6q16.1,3p21.31 and 17q21.31may account for the shared causes between educational attainment and neuropsychiatric disorders.
More importantly, the genetic underpinnings of EA and NPDs lack a comprehensive understanding.Large-scale genome-wide association studies (GWASs) have identified 74 susceptibility loci for EA (18), including specific loci at 3p21.31 and 17q21.31.These loci are also associated with the risk of neuropsychiatric disorders (NPDs) (19)(20)(21)(22)(23)(24)(25), indicating a potential shared genetic architecture between them.Genetic studies can provide valuable insights into the specific biological processes contributing to the co-occurrence of EA and NPDs, potentially serving as targets for early prevention or treatment of NPDs.Yet, our understanding of the shared genetics underlying possible genetic correlations, specific shared genetic variants, and the corresponding biological pathways remains poorly characterized.
In this study, we conducted a large-scale cross-trait genetic analysis to explore the genetic correlations and potential causality between EA and 13 NPDs, and attempted to delineate the precise shared genetic variants and biological pathways linking them.

Study design, data source, and study population
Our overall study design was shown in Figure 1.GWAS summary statistics for EA were obtained from the meta-analysis of 64 independent cohorts (N = 293,723) (18).Summary statistics for ADHD, AN, ASD, AUD, BIP, CUD, MDD, OCD, PTSD, SCZ, and TS were available at PGC (psychiatric genomics consortium).Summary statistics for ALS and AD were available at the Project MinE and CTGlab (complex trait genetics lab), respectively.All participants were of European ancestry to avoid confounding bias.Details of study designs and each of the GWASs were present in Supplementary Table S1.The potential sample overlap between the GWASs of EA and NPDs is detailed in the Supplementary methods "sample overlap" section and Supplementary Table S2.

Genetic correlation analysis
Genetic correlations (r g ) were estimated using the High-Definition Likelihood (HDL) 1 method, which yielded more precise genetic correlation estimates compared to linkage disequilibrium score regression (LDSC) 2 (26).We estimated r g through comparisons within specific traits, utilizing LD reference data derived from 335,265 individuals of British ancestry from the United Kingdom Biobank (UKB).We also estimated annotation-specific genetic correlations using 13 extensive annotations through partitioned LD score regression (27) (Supplementary methods).

Cross trait meta-analysis
Cross-trait association analysis integrates effect estimates and standard errors from GWAS summary statistics to assess the association between an SNP and two traits (28).We employed the heterogenous variant of the cross-phenotype association test statistic (SHet), calculated using a fixed-effects model, offering enhanced statistical power in the presence of heterogeneity among studies.A significance level of p < 3.8 × 10 −9 (5 × 10 −8 /13 for 13 meta-analysis tests) was employed in the genome-wide cross-trait meta-analysis (Supplementary methods).

Mendelian randomization
Bidirectional Mendelian randomization (MR) analysis was performed to estimate the causal relationship between EA and NPDs (29).Additionally, given the potential impact of horizontal pleiotropy as a confounding factor that may bias the estimates and inflate test statistics in MR analysis, we employed MR-Egger regression and the MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO)3 methods (30) to identify and address potential horizontal pleiotropy.We chose the inverse variance-weighted (IVW) method as the primary approach.Bonferroni correction was used for the p value (p < 0.05/13).A detailed explanation of the causal effect estimates can be found in the Supplementary material.

Transcriptome-wide association studies
To uncover common genes that elucidate the shared mechanisms behind the genetic correlations observed between EA and NPDs, we further delved into the gene level to examine tissue-specific expression-trait associations, as well as shared expression-trait associations between EA and each genetically correlated NPD, using transcriptome-wide association studies (TWAS) 4 (31).We applied the Benjamini-Hochberg correction to identify significant expression-trait associations, accounting for multiple comparisons across all gene-tissue pairs tested for each trait (totaling ~230,000 gene-tissue pairs in total; significant expression-trait associations were defined as P Benjamini-Hochberg < 0.05) (Supplementary methods).

Secondary analyses 2.3.1 Fine-mapping credible set analysis
For each shared locus between EA and NPDs that met the significance criteria in the cross-trait meta-analysis, variants within 500 kb of the lead SNP were extracted, and then a 99% credible set of potentially causal SNPs was identified utilizing the Bayesian likelihood fine-mapping algorithm (FM-summary) 5 (32).If no variants with a linkage disequilibrium r 2 greater than or equal to 0.9 with the lead SNP were found, the FM-summary cannot be performed.

Co-localization analysis
Summary statistics for variants within 500 kb of the lead SNP at shared loci between EA and NPDs were extracted.We employed the 'coloc'6 package in R to conduct genetic colocalization analysis, calculating the probability of shared common genetic causal variants between the two traits.In this study, loci with a probability (H4) equal or greater than 0.5 were considered to be colocalized (33).Workflow of the study.

Tissue-specific enrichment analysis
To assess whether the shared gene sets identified through crosstrait meta-analysis exhibited significant enrichment or tissue-specific expression, we performed tissue-specific enrichment analysis (TSEA) using the 'deTS'7 package in R (34) (Supplementary methods).

Functional enrichment analysis
For gaining biological insights into the shared genes identified (P meta < 3.8 × 10 −9 ) in the cross-trait meta-analysis, we employed the ClueGO8 plug-in within the Cytoscape tool to assess the enrichment of gene sets in the Gene Ontology (GO) biological processes and visualized the relationships between genes and GO terms (35,36).We applied the Benjamini-Hochberg procedure to correct for multiple testing (False Discovery Rate < 0.05).
To further investigate whether specific functional regions played a disproportionate role in the association between EA and related NPDs, we performed a partitioned genetic correlation analysis across 13 major functional categories (27).In general, annotation-specific LDSC yielded results consistent with the HDL estimate, but indicated a notably stronger genetic correlation at conserved regions and DNase I hypersensitivity sites (DHS) (Supplementary Figure S1; Supplementary Table S3).These findings provided robust evidence of substantial genetic correlations between EA and related NPDs.

Cross-trait meta-analysis between EA and NPDs
The QQ plot and corresponding inflation lambda for the crosstrait GWAS could be found in Supplementary Figure S2.In total, we identified 37 independent loci shared between EA and NPDs (P meta < 3.8 × 10 −9 and single-trait p < 1 × 10 −5 ) (Table 2).We also determined the credible set of SNPs for each of these shared loci (Supplementary Tables S4-S12).We discovered a previously unreported locus at 20q13.13 (index SNP: rs71351952) that is shared between EA and PTSD.This locus is associated with the proteincoding gene ARFGEF2.ARFGEF2 encodes the large brefeldin A (BFA)-inhibited GEF2 protein (BIG2), which is involved in vesicle and membrane trafficking originating from the trans-Golgi network (37).Vesicle trafficking plays a crucial role in regulating cell proliferation and migration during the development of the human cerebral cortex (38).
We identified seven, seven, three, six, three, five, two, and three independent loci shared between EA and ADHD, AN, ASD, AUD, BIP, CUD, MDD, or TS, respectively (Table 2).We did not find a shared locus between EA and AD, ALS, or OCD.Our findings revealed both commonalities and distinctions in the shared genetic factors between EA and various NPDs.One significant region, located at 6q16.1(near MIR2113), was shared between EA with    .Interestingly, this region was also the strongest shared signal between EA with ASD, AUD, BIP, or TS.MIR2113 was found to be associated with regulatory elements, including open chromatin, histone modifications, DNase hypersensitive sites, and position weight matrix sites, indicating that the associated SNPs are located within regions of active transcription and may have a regulatory role in transcription (39).

Colocalization analysis
Out of the 37 independent loci shared between EA and NPDs, seventeen of these loci exhibited colocalization at the same candidate causal variant (PPH 4 > 0.5), sixteen loci exhibited colocalization with different candidate causal variants (PPH 3 > 0.5), while four loci did not exhibit colocalization (Supplementary Table S13).

Tissue-specific enrichment analysis
We detected ten distinct tissues that exhibited significantly enriched expression of genes associated with cross-traits (The query gene list was extracted from the shared genes between EA with ADHD, ASD, AUD, BIP, and TS) (Figure 2).The three most significantly enriched tissues were all part of the nervous system, including the brain-anterior cingulate cortex, brain-frontal cortex, and brain-hypothalamus.

Functional enrichment analysis
Our finding showed that the 117 genes shared between EA and NPDs were enriched in axon development, positive regulation of nervous system development, synapse organization, learning or memory, purine nucleotide metabolic process, cAMP (cyclic adenosine monophosphate)-mediated signaling, cellular response to abiotic stimulus, protein deubiquitinating, monocarboxylic acid transport, and mRNA splicing, via spliceosome (Figure 3; Supplementary Table S14).

Causal inference
We then tested the causality between EA and related NPDs using Two sample MR. 9 MR-Egger regression and MR-PRESSO method indicated little significant horizontal pleiotropy (most P MR_Egger > 0.05) (Supplementary Tables S15, S16).Table 3 indicates a causal relationship between EA and ADHD.The findings suggest that a genetic predisposition to higher EA is associated with a reduced risk of ADHD [OR IVW (odds ratio): 0.50 per SD increment; 95% CI: 0.39-0.63;P IVW = 1.62 × 10 −8 ], and ADHD is estimated to have a causal effect of −2.8 months on EA, corresponding to a doubling in ADHD prevalence [Causal effect = −2.8;95%CI: −3.9 to −1.8; P IVW = 8.46 × 10 −8 ].In addition, there was evidence suggesting a negative causal effect of genetic predisposition to higher EA on a reduced risk of AD [OR IVW : 0.94 per SD increment; 95%CI: 0.91 to 0.98; P IVW = 0.0028] and CUD [OR IVW : 0.50 per SD increment; 95%CI: 0.35 to 0.73; P IVW = 0.0003].Moreover, we observed that higher EA had a positive causal effect on ASD [OR IVW : 1.45 per SD increment; 95%CI: 1.15 to 1.83; P IVW = 0.0019].Other related NPDs (including ALS, AN, AUD, BIP, MDD, OCD, PTSD, SCZ, and TS) were not related to higher EA and vice versa after correction.Sensitivity analysis confirmed the reliability of our results, indicating that they were not influenced by horizontal pleiotropy or outliers (Supplementary Table S17).

Single trait TWAS and shared genetics between EA and NPDs from TWAS
We then conducted a gene-level analysis to investigate the common genes expressed in both EA and related NPDs.In total, we identified 5,874 gene-tissue pairs that exhibited a significant association with EA across 48 GTEx tissues after applying the Benjamini-Hochberg correction.Additionally, we found 495, 2,108, 962, 487, 1,227, 1707, 128, and 207 gene-tissue pairs associated with ADHD, AD, AN, ASD, AUD, BIP, CUD, and MDD, respectively (Supplementary Figures S3-S5).
Interestingly, in strong alignment with the findings from CPASSOC,10 certain loci were common to both EA and multiple NPDs.For instance, locus 3p21.31,harboring multiple genes, was commonly TWAS significant for EA, ADHD, AN, BIP, CUD, and MDD and was simultaneously the most prominent shared locus among EA and AN, BIP, or CUD.Notably, we found NAT6 (3p21.31)was commonly significant for EA, ADHD, AN, BIP, and CUD from TWAS.NAT6 encodes a member of the N-acetyltransferase family.N-acetyltransferases modify proteins by transferring acetyl groups from acetyl-CoA to the N-termini of protein substrates.The encoded protein is a cytoplasmic N-acetyltransferase with a substrate specificity for N-termini that is enriched for acidic residues (47).This gene was significantly related to CUD through genetically regulated gene expression (23).Locus 17q21.31 was co-significant for EA, ASD, and AUD and was also the most significant shared locus between them (Supplementary results).

Discussion
Our analysis yielded four main findings.First and foremost, we presented compelling evidence demonstrating a significant genetic correlation between EA and 12 NPDs using the HDL method and partitioned LDSC.Second, cross-trait meta-analysis identified independent shared loci between EA and NPDs.Functional analysis revealed that these shared loci were predominantly enriched in brain tissues, associated with neurogenesis, synapse formation, and the purine nucleotide metabolic and cAMP-mediated signaling pathways.Third, TWAS uncovered shared gene expressions between EA and NPDs in tissues from the nervous system, endocrine system, and digestive system.Fourth, MR analyses revealed bidirectional causality between EA and ADHD, i.e., higher genetically predicted EA was associated with lower risk of ADHD, and genetically predicted ADHD was associated with lower EA.Our results advanced our understanding of the shared genetic contributions to EA and NPDs, elucidated the causal connections between these traits, provided insight into the potential regulatory roles of shared genetics, and shed light on the etiologies and mechanisms underlying the coexistence of low EA and NPDs.Low educational attainment is indicative of the neurobiological abnormalities that underlie NPDs and has been proposed as a potential premorbid marker for these disorders (48).Based on our results, we found significantly negative genetic correlations of EA with ADHD, ALS, AD, CUD, MDD, PTSD, and TS, and positive genetic associations of EA with AN, ASD, AUD, BIP, and OCD.The positive association of EA with AN, ASD, AUD, BIP, and OCD may be driven by perfectionism (49).Meta-analysis suggested that perfectionism contributes to higher academic achievement, and several studies have confirmed that higher levels of perfectionism are associated with a greater risk of AN (50), ASD (51), AUD (52), BIP (53), and OCD (54).We observed that EA exhibited a positive association with SCZ (r g = 0.07, p = 1.01 × 10 −3 ), although this significance was no longer present after Bonferroni correction.However, recent studies have demonstrated a negative association between higher EA and SCZ using the proxy phenotype method (55).The mechanisms underlying this counterintuitive finding of an inverse relationship between EA and SCZ remained unknown.It may be partially attributed to the genetic heterogeneity of EA, which can be deconstructed into sub-phenotypes (including higher IQ, greater openness, and increased conscientiousness) that exhibit imperfect genetic correlations among themselves.If the various symptoms of SCZ also possess non-identical genetic architectures, this could lead to conflicting phenotypic and genetic correlation outcomes (56).
According to a previous study, the hallmark biology process of EA was neural development (18).It was well known that excitatory synapse, developing brain, neuronal morphogenesis, and central nervous system morphogenesis were significant biological pathways for ADHD (19), ASD (24), MDD (22), and CUD (23), respectively.Functional enrichment for the consensus set.Functional enrichment analysis results using the ClueGO method in Cytoscape (34).Each dot represented a gene or a GO term.Dots in the same color were considered from the same functional group by ClueGO annotation.Gene names were highlighted in red.Each edge indicated the gene was a component gene of the linked GO term.Several of these pathways were associated with EA regulation, consistent with the results of our functional enrichment analysis, suggesting a common biological mechanism for these disorders.Shared genes identified from CPASSOC and TWAS also supported the shared biological pathways between EA and NPDs.Noteworthy, our findings emphasized the potentially intriguing roles of novel associations involving ARFGEF2 in the context of the relationship between EA and PTSD.Additionally, our attention was directed towards genes located within the 6q16.1,3p21.31, and 17q21.31loci.Genes shared between EA and NPDs exhibited enrichment in various brain tissues, notably the cerebellum, suggesting that these disorders may be attributed to nervous system dysfunction.For instance, mutations in the ARFGEF2 gene, responsible for encoding the BIG2 protein and identified in patients with autosomal recessive periventricular heterotopia with microcephaly, are closely associated with neuronal migration disorders (57).Hong et al. (58) reported that BIG2 is a RhoA interactor that regulates dendrite morphogenesis and Golgi polarization during neuronal cell development.BIG2 also contains A-kinase anchoring protein (AKAP) sequences that can act as scaffolds for multimolecular assemblies that facilitate and limit cAMP signaling temporally and spatially (59) indicating the axon development and cAMP-mediated signaling pathways shared between EA and PTSD.Notably, MIR2113 (6q16.1)has been associated with accelerated decline in episodic memory (60) and neurodegenerative diseases such as AD and ALS (61).It is also biologically and genetically implicated in EA (62), ASD (24), AUD (63), BIP (64), and TS (65).The genealogical study demonstrated that the 6q16.1 deletion will cause developmental delay and intellectual disability (66).
Additionally, a microdeletion located at 3p21.31 was previously associated with mild cognitive disability (67), CNS abnormalities, and developmental delay (43).A comprehensive transcriptomic study identified CAMKV within the synaptic neuropil (68), and its protein synthesis is governed by neuronal activity.Moreover, the knockdown of CAMKV in mouse hippocampal CA1 pyramidal neurons adversely  (25) and ASD (24).Recurrent microdeletions of 17q21.31,which include the CHRH1 gene, have been identified as a relatively common cause of intellectual disability.This underscores the potential involvement of CHRH1 in neurodevelopmental disorders (71).WNT3 encoded the proto-oncogene protein Wnt-3, which played a crucial role in brain development and patterning.It was also implicated in the pathogenesis of diseases such as Parkinson's, Alzheimer's, and schizophrenia (72).More importantly, 17q21.31 was also co-significant for EA, ASD, and AUD from TWAS.CPASSOC, functional analysis, and TWAS underlined the molecular pathways in the regulation of axon development, nervous system development, synapse organization, learning or memory, purine nucleotide metabolic process, cAMP-mediated signaling, and protein deubiquitinating, further underscores the potential significance of neurological dysregulation throughout the entire body in predisposition to both low EA and NPDs.
In addition to the nervous system, what particularly interested us is the enrichment of tissues in the cardiovascular system, endocrine system, and musculoskeletal system.This suggested that the shared pathway between EA and NPDs may have significant functions beyond brain tissue.For example, AN was associated with cardiac abnormalities including hypotension, and dysregulation in peripheral vascular contractility (73).Medication treatment in ADHD patients was associated with an increased risk of arrhythmias and myocardial infarctions (74).EA may influence the risk of AN and ADHD through alterations in blood pressure or high-density lipoprotein cholesterol in the cardiovascular system (75).Consistent with the result from Shadrin (76), musculoskeletal tissue was TWAS significant tissue between EA and ADHD.Taken together, TWAS provided additional evidence that the expressed enriched genes may not be specifically expressed in brain tissue; rather, it appeared to be prevalent throughout the body in brain, metabolic, musculoskeletal, and endocrine-related organs.
In our bidirectional MR analysis, we observed a causal relationship between higher EA and a reduced risk of ADHD (OR IVW = 0.33).This finding aligns with results obtained from Mendelian randomization analysis, which indicated a 0.3-fold decrease in ADHD risk per SD increase in EA (15).Conversely, the onset of ADHD was associated with a subsequent reduction in EA (Causal effect = −2.8months), suggesting that targeting EA could be a promising therapeutic strategy for ADHD.The bidirectional causal association between EA and ADHD is highly similar to the results of a recent MR study conducted by Artigas et al. (77).Research suggested that attentional strengths may be a potential phenotype that simultaneously drives the association between lower EA and ADHD (51,78).Our results also provided evidence of a causal relationship between higher EA and a decreased risk of AD and CUD, while increasing the risk of ASD.These findings were in line with discoveries from recent studies (14,15,79,80).MR analysis demonstrated that lower EA in individuals may exacerbate the risk of NPDs (especially for ADHD, AD, and CUD).Therefore, to improve this relationship, increased efforts to reduce dropout rates and create more educational opportunities are critical to ensure optimal long-term EA for this vulnerable group (81).The mechanisms underlying the causality between EA and NPDs remain to be elucidated, but the shared CPASSOC and TWAS genes could offer new insights into revealing shared etiologies and mechanisms between them.
Our results have provided clear insights into the shared genetic architecture between EA and NPDs.We emphasized the importance of implementing NPD prevention programs for individuals with lower levels of EA.However, our study has several limitations.First, while TWAS has enhanced the power to detect significant expression trait associations, the relatively smaller sample size for metabolic traits and GTEx reference panels in certain tissues may not be sufficient to detect signals with small to moderate effects.Second, although we observed higher instrumental estimates from EA to NPDs than in the reverse path, we cannot rule out the possibility that EA mediates the effects of certain upstream latent factors on NPDs.Third, given that the traits we investigated involved EA, it cannot be ruled out that the detected causal effects were caused by dynastic effects or assortative mating (82).Fourth, our study was confined to individuals of European ancestry, and the shared genetics in other ethnic groups remain uncertain.Therefore, future research should encompass other ethnic groups, and exploration into how mediating variables, such as family background, personality characteristics, and intelligence, influence the relationship between EA and NPDs is encouraged.
In conclusion, our study provided evidence of significant genetic correlations and instrumental causality between EA and NPDs.The identification of genetic loci associated with both EA and the risk of NPDs presents therapeutic opportunities to enhance EA and reduce the risk of NPDs.Furthermore, our study advanced our understanding of EA and revealed a common etiology for the coexistence of low EA and NPDs.

FIGURE 2 Tissue
FIGURE 2Tissue-specific enrichment result for query gene list.X-axis: groups of genes defined according to different thresholds based on Pascal value of p. Y-axis: 47 GTEx tissues were used as the reference panel.The top three significant tissues (adjusted p < 0.1 from Fisher's Exact Test after Benjamini-Hochberg correction) were marked in numbers.

FIGURE 4
FIGURE 4Numbers of significant genes related to ADHD, AN, ASD, AUD, CUD, AD, BIP, and MDD, and the number of shared genes with EA.Significant genes were identified by P BH < 0.05.GTEx = genotype-tissue expression project; GWAS = genome-wide association studies; TWAS = transcriptome-wide association study; NSTSGR = Number of Shared TWAS Significant Genes with EA; EA: Educational attainment; ADHD: Attention-deficit/hyperactivity disorder; AD: Alzheimer's dementia; AN: Anorexia nervosa; ASD: Autism spectrum disorders; AUD: Alcohol use disorders; BIP: Bipolar disorder; CUD: Cannabis use disorder; MDD: Major depressive disorder.(A): No. of TWAS Significant Genes for ADHD and No. of Overlapped Genes with EA. (B): No. of TWAS Significant Genes for AN and No. of Overlapped Genes with EA. (C): No. of TWAS Significant Genes for ASD and No. of Overlapped Genes with EA. (D): No. of TWAS Significant Genes for AUD and No. of Overlapped Genes with EA. (E): No. of TWAS Significant Genes for CUD and No. of Overlapped Genes with EA. (F): No. of TWAS Significant Genes for AD and No. of Overlapped Genes with EA. (G): No. of TWAS Significant Genes for BIP and No. of Overlapped Genes with EA. (H): No. of TWAS Significant Genes for MDD and No. of Overlapped Genes with EA.

TABLE 1
Genetic correlation of educational attainment with related neuropsychiatric traits (NPDs) using high-definition likelihood method.

TABLE 3
Causal inference between educational attainment and neuropsychiatric traits using two-sample MR.The causal effect size was the beta coefficient from linear or logistic regression models for corresponding outcome; MR: Mendelian randomization; Direction: Forward means the causal effect size of EA on neuropsychiatric traits, Reverse means the causal effect size of neuropsychiatric traits on EA; N_snp: number of instrumental variables; The threshold of significance was at the Bonferroni-adjusted level p-value < 0.0038(0.05/13).AD, Alzheimer's dementia; ADHD, Attention deficit/hyperactivity disorder; ALS, Amyotrophic lateral sclerosis; AN, Anorexia nervosa; ASD:,Autism spectrum disorders; AUD, Alcohol use disorders; BIP, Bipolar disorder; CUD, Cannabis use disorder; MDD, Major depressive disorder; OCD, Obsessive compulsive disorder; PTSD, Posttraumatic stress disorder; TS, Tourette's syndrome. *,