ORIGINAL RESEARCH article

Front. Psychiatry, 04 January 2024

Sec. Mood Disorders

Volume 14 - 2023 | https://doi.org/10.3389/fpsyt.2023.1303430

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

  • 1. Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Genetics, Peking University Cancer Hospital & Institute, Beijing, China

  • 2. Shenzhen Health Development Research and Data Management Center, Shenzhen, China

  • 3. Department of Occupational and Environmental Health Sciences, School of Public Health, Peking University, Beijing, China

Abstract

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.31 may account for the shared causes between educational attainment and neuropsychiatric disorders.

1 Introduction

Educational attainment (EA) is moderately heritable (heritability = 40.0%) from behavior-genetic studies and an important correlate of socioeconomic status, and health outcomes (1, 2). Neuropsychiatric disorders (NPDs) are characterized by affective, cognitive, and behavioral abnormalities stemming from underlying cerebral dysfunction or secondary effects of systemic disease (3). Previous observational studies have found negative associations between EA and risk of attention-deficit/hyperactivity disorder (ADHD) (4), Alzheimer’s dementia (AD) (5), amyotrophic lateral sclerosis (ALS) (6), alcohol use disorders(AUD) (7), cannabis use disorder (CUD) (7), major depression disorder (MDD) (8), obsessive-compulsive disorder (OCD) (9), posttraumatic stress disorder (PTSD) (10), or Tourette’s syndrome (TS) (11). Other observational studies suggest a positive relation between EA and bipolar disorder(BIP) (8), autism spectrum disorder (ASD) (12), or anorexia nervosa (AN) (13). Therefore, not all the correlation results between EA and NPDs exhibit isotropic behavior in nature. ADHD is a prevalent psychiatric condition manifesting as persistent symptoms of inattention, hyperactivity, and impulsivity, potentially leading to worse education outcomes, suggesting a potential bidirectional association between educational attainment (EA) and ADHD. Given that genetic factors are less susceptible to reverse causation nor external confounding factors, it is informative to investigate the bidirectional causality behind these associations using Mendelian randomization. To date, several Mendelian randomization studies have reported causal associations between EA and different NPDs [e.g., AD (14), ADHD (15), ASD (15), ALS (16), and CUD (17)]. The present study will provide further evidence of potential causal associations between EA and NPDs.

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–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.

2 Methods

2.1 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.

Figure 1

2.2 Main analyses

2.2.1 Genetic correlation analysis

Genetic correlations (rg) 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 rg 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).

2.2.2 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).

2.2.3 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.

2.2.4 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 PBenjamini-Hochberg < 0.05) (Supplementary methods).

2.3 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 r2 greater than or equal to 0.9 with the lead SNP were found, the FM-summary cannot be performed.

2.3.2 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).

2.3.3 Tissue-specific enrichment analysis

To assess whether the shared gene sets identified through cross-trait 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).

2.3.4 Functional enrichment analysis

For gaining biological insights into the shared genes identified (Pmeta < 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).

3 Results

3.1 Genetic correlations of EA with NPDs

EA was negatively genetically correlated with ADHD (rg = −0.49; p = 1.36 × 10−111), ALS (rg = −0.29; p = 2.77 × 10−04), AD (rg = −0.25; p = 1.84 × 10−07), CUD (rg = −0.41; p = 5.44 × 10−5), MDD (rg = −0.19; p = 5.52 × 10−15), PTSD (rg = −0.29; p = 3.37 × 10−10), and TS (rg = −0.17; p = 6.87 × 10−15), and positively genetically associated with AN (rg = 0.23; p = 8.15 × 10−20), ASD (rg = 0.17; p = 3.08 × 10−13), AUD (rg = 0.20; p = 3.97 × 10−19), BIP (rg = 0.20; p = 2.27 × 10−22), and OCD (rg = 0.35; p = 1.45 × 10−7) (Table 1). No significant genetic correlation was observed for EA and SCZ.

Table 1

NPDs traitrgrg, SErg, 95% CIp value*
ADHD−0.490.022−0.53 to −0.451.36E-111
ALS−0.290.079−0.44 to −0.142.77E-04
AD−0.250.048−0.34 to −0.161.84E-07
AN0.230.0250.18 to 0.288.15E-20
ASD0.170.0230.12 to 0.223.08E-13
AUD0.200.0220.16 to 0.243.97E-19
BIP0.200.0210.16 to 0.242.27E-22
CUD−0.410.043−0.49 to −0.335.44E-22
MDD−0.190.024−0.24 to −0.145.52E-15
OCD0.350.0660.22 to 0.481.45E-07
PTSD−0.290.046−0.38 to −0.23.37E-10
SCZ0.070.0210.03 to 0.111.01E-03
TS−0.170.043−0.25 to −0.096.87E-05

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

Summary statistics for each trait were combined with Hapmap3 SNPs, excluding the HLA region, to estimate rg. *p value < 0.05/13. ADHD, Attention-deficit/hyperactivity disorder; ALS, Amyotrophic lateral sclerosis; 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; OCD, Obsessive-compulsive disorder; PTSD, Posttraumatic stress disorder; SCZ, Schizophrenia; TS, Tourette’s syndrome.

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.

3.2 Cross-trait meta-analysis between EA and NPDs

The QQ plot and corresponding inflation lambda for the cross-trait GWAS could be found in Supplementary Figure S2. In total, we identified 37 independent loci shared between EA and NPDs (Pmeta < 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 protein-coding 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).

Table 2

TraitIndex SNPPOSEANEATrait 1Trait 2PCPASSOCGenes
BETAPBETAP
ADHDrs673253chr1p34.2TC0.0158.28E-10−0.0921.53E-111.17E-14[KDM4A, PTPRF]
rs3791101chr1p34.1AG0.0127.14E-06−0.0975.3E-117.98E-11[MIR6079, ST3GAL3]
rs549845chr1p34.2AG0.0155.78E-08−0.0931.38E-101.832E-10[KDM4A, PTPRF, ST3GAL3]
rs62260755chr3p21.31CG0.0213.37E-12−0.0806.74E-071.026E-12[CAMKV, TRAIP, CDHR4, MST1R, UBA7]
rs304137chr5q14.3AG−0.0161.78E-100.0641.79E-063.228E-16[MEF2C]
rs12653396chr5q14.3AT−0.0132.14E-070.0731.13E-073.725E-11[LINC00461]
rs4839923chr6q16.1AG−0.0131.1E-070.0651.9E-064.286E-11[MIR2113]
ANrs13093385chr3p21.31AT−0.0283.43E-22−0.0711E-064.377E-25[APEH, BSN, MST1, RNF123]
rs13096760chr3p21.31TC0.0155.18E-090.0651.45E-062.463E-17[APEH, BSN, MST1, RNF123]
rs9821797chr3p21.31AT0.022.4E-07−0.1576.99E-157.55E-15[IP6K2, NCKIPSD, MIR4793]
rs73073015chr3p21.31AG0.0239.32E-08−0.1733.24E-133.887E-13[DAG1, NICN1, AMT, RHOA, GPX1]
rs9832454chr3p21.31TC0.0232.27E-07−0.1251.68E-071.428E-11[GNAI2, GNAT1, MIR566,
SEMA3F, RBM5, RBM6]
rs6967776chr7q33AG−0.0161.54E-07−0.0759.43E-067.207E-10[EXOC4]
rs705696chr12q13.2AG0.0182.5E-11−0.0647.67E-069.125E-12[ERBB3]
ASDrs9320913chr6q16.1AC0.0242.05E-210.0671.51E-061.196E-36[MIR2113]
rs1106761chr8q24.3AG−0.0174.08E-11−0.0657.16E-061.514E-09[MROH5]
rs62057121chr17q21.31AG−0.0182.8E-080.0785.94E-061.066E-11[CRHR1, KANSL1-AS1, MAPT,
MAPT-AS1, SPPL2C, STH]
AUDrs2101975chr4q24AG0.0132.75E-07−0.0065.5E-061.655E-14[TET2, PPA2]
rs1338549chr6q16.1TG0.0218.26E-17−0.0056.73E-061.4E-30[MIR2113]
rs13266268chr8q24.3TC−0.0174.1E-11−0.0063.95E-066.634E-12[MROH5]
rs113925422chr17q21.31AT−0.0181.02E-07−0.0082.6E-083.904E-14[MAPT, MAPT-AS1, SPPL2C, STH]
rs538628chr17q21.31CG−0.0182.16E-08−0.0083.27E-084.408E-14[NSF, WNT3]
rs1879581chr17q21.31TC0.0162.1E-06−0.0075.16E-065.599E-11[PLEKHM1]
BIPrs12754946chr1p31.1TC0.0121.77E-060.0662.39E-061.19E-10[AK5]
rs9320913chr6q16.1AC0.0242.05E-210.0701.27E-075.566E-37[MIR2113]
rs10429537chr9p21.3CG−0.0251.02E-21−0.0671.21E-066.978E-19[ELAVL2]
CUDrs12122664chr1p22.2AC0.0136.22E-07−0.0831.92E-062.768E-09[BARHL2]
rs11711407chr3p21.31AG−0.0203.73E-150.0776.41E-062.34E-23[GNAT1, SEMA3F, SLC38A3]
rs35926495chr3p21.31TC−0.0171.31E-110.0882.73E-071.608E-12[GNAT1, GNAI2, LSMEM2, MIR566,
NAT6, SEMA3B, SLC38A3]
rs9467773chr6p22.2AG0.0123E-06−0.0749.13E-061.868E-11[BTN1A1, HCG11]
rs7783012chr7q31.1AG−0.0134.69E-070.1011.84E-094.17E-12[FOXP2]
MDDrs7531118chr1p31.1TC−0.0121.84E-06−0.0452.15E-084.392E-10[RPL31P12]
rs76025409chr5q21.2CG−0.0125.14E-060.0572.32E-115.972E-11[RP11-6 N13.1]
PTSDrs71351952chr20q13.13TC0.0139.6E-060.0784.94E-061.772E-09[ARFGEF2, CSE1L,
LOC102723483, STAU1]
TSrs9401593chr6q16.1AC−0.0241.71E-210.1287.87E-074.122E-35[MIR2113]
rs12154193chr6q16.1CG−0.0171.38E-110.1169.62E-063.961E-18[MIR2113]
rs1619561chr12q24.31CG−0.0231.44E-140.1445.06E-066.059E-15[MPHOSPH9]

Cross-trait meta-analysis results between educational attainment and neuropsychiatric traits (PCPASSOC < 3.8 × 10−9 and Psingle-trait < 1 × 10−5).

EA, Effect allele; NEA, non-effect allele; Trait 1 refers to the EA trait; Trait 2 refers to the corresponding neuropsychiatric traits. PCPASSOC means the cross-trait meta-analysis p value. EA, educational attainment.; ADHD, Attention deficit/hyperactivity disorder; AN, Anorexia nervosa; ASD, Autism spectrum disorders; AUD, Alcohol use disorders; BIP, Bipolar disorder; CUD, Cannabis use disorder; MDD, Major depressive disorder; PTSD, Posttraumatic stress disorder; TS, Tourette’s syndrome.

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 ADHD (index-SNP: rs4839923, Pmeta = 4.29 × 10−11), ASD (index-SNP: rs9320913, Pmeta = 1.20 × 10−36), AUD (index-SNP: rs1338549, Pmeta = 1.40 × 10−30), BIP (index-SNP: rs9320913, Pmeta = 5.57 × 10−37), and TS (index-SNP: rs9401593, Pmeta = 4.12 × 10−35; index-SNP: rs12154193, Pmeta = 3.96 × 10−18). 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).

Another prominent multigenic locus, lied in 3p21.31 (mapped genes: CAMKV, and BSN) was common for EA and ADHD (index-SNP: rs62260755, Pmeta = 1.03 × 10−12), AN (index-SNP: rs13093385, Pmeta = 4.38 × 10−25; index-SNP: rs13096760, Pmeta = 2.46 × 10−17; index-SNP: rs9821797, Pmeta = 7.55 × 10−15; index-SNP: rs73073015, Pmeta = 3.89 × 10−13; index-SNP: rs9832454, Pmeta = 1.43 × 10−11), CUD (index-SNP: rs11711407, Pmeta = 2.34 × 10−23; index-SNP: rs35926495, Pmeta = 1.61 × 10−12). The genes mapped to these loci are all protein-coding and have a role in brain neurology and brain function. For example, CAMKV acts as a convergence point for transducing Ca2+ signals to the neuronal cytoskeleton through RhoA (Ras homolog gene family, member A), playing a crucial role in normal synaptic transmission that underlies brain function (40). BSN (Bassoon) is known to participate in presynaptic transmission in excitatory neurons; its loss results in reduced transmission of the excitatory neurotransmitter glutamate in both the brain and retina (41, 42). The loss of 3p21.31 may contribute to intellectual disability and developmental delays, potentially playing a functional role in the observed intellectual disability and developmental delays in individuals with ADHD and low academic attainment (43). Additional results are available in the Supplementary materials.

Additionally, two loci at chromosome 8q24.3 (MROH5, index-SNP: rs1106761, Pmeta = 1.51 × 10−9 for ASD; index-SNP: rs13266268, Pmeta = 6.63 × 10−12 for AUD) and 17q21.31 (CRHR1 and WNT3, index-SNP: rs62057121, Pmeta = 1.07 × 10−11 for ASD; index-SNP: rs113925422, Pmeta = 3.90 × 10−14; index-SNP: rs538628, Pmeta = 4.41 × 10−14; index-SNP: rs1879581, Pmeta = 1.19 × 10−10 for AUD) were common for EA, ASD, and AUD. MROH5 is a risk factor for ASD (24), but its function is not yet clear. Chromosome 17q21.31 region spanned several genes. For example, CRHR1 (corticotropin receptor gene) has been linked to alcohol use in both animals and humans (44). WNT3 has demonstrated its role not only as an axon guidance molecule but also notably as a gradient for retinotopic mapping along the medial-lateral axis (45). Furthermore, we identified additional shared loci between EA and ADHD, AN, AUD, BIP, CUD, or MDD, consistent with previous research (19, 21–23, 25, 46).

3.3 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 (PPH4 > 0.5), sixteen loci exhibited colocalization with different candidate causal variants (PPH3 > 0.5), while four loci did not exhibit colocalization (Supplementary Table S13).

3.4 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.

Figure 2

3.5 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).

Figure 3

3.6 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 PMR_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 [ORIVW (odds ratio): 0.50 per SD increment; 95% CI: 0.39–0.63; PIVW = 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; PIVW = 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 [ORIVW: 0.94 per SD increment; 95%CI: 0.91 to 0.98; PIVW = 0.0028] and CUD [ORIVW: 0.50 per SD increment; 95%CI: 0.35 to 0.73; PIVW = 0.0003]. Moreover, we observed that higher EA had a positive causal effect on ASD [ORIVW: 1.45 per SD increment; 95%CI: 1.15 to 1.83; PIVW = 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).

Table 3

OutcomeDirectionMethodN_snpCausal_Effect_Size95% CIp_value
ADForwardInverse variance weighted63−0.059(−0.098, −0.020)0.0028
ReverseInverse variance weighted240.030(−0.029, 0.088)0.3191
ADHDForwardInverse variance weighted54−0.699(−0.941, −0.456)1.62E-08
ReverseInverse variance weighted8−0.095(−0.130, −0.060)8.46E-08
ALSForwardInverse variance weighted58−0.264(−0.548, 0.019)0.0679
ReverseInverse variance weighted4−0.015(−0.037, 0.007)0.1799
ANForwardInverse variance weighted440.333(0.078, 0.588)0.01055
ReverseInverse variance weighted110.012(−0.025, 0.049)0.5185
ASDForwardInverse variance weighted530.369(0.136, 0.602)0.0019
ReverseInverse variance weighted90.004(−0.051, 0.059)0.8804
AUDForwardInverse variance weighted610.040(0.009, 0.071)0.01048
ReverseInverse variance weighted6−0.178(−0.539, 0.183)0.3344
BIPForwardInverse variance weighted470.245(−0.016, 0.506)0.0662
ReverseInverse variance weighted110.0360(0.000, 0.071)0.0481
CUDForwardInverse variance weighted63−0.688(−1.058, −0.318)0.0003
ReverseInverse variance weighted11−0.039(−0.065, −0.014)0.0024
MDDForwardInverse variance weighted39−0.075(−0.290, 0.140)0.4956
ReverseInverse variance weighted40.010(−0.133, 0.153)0.8921
OCDForwardInverse variance weighted46−0.072(−0.678, 0.534)0.8150
ReverseInverse variance weighted11−0.004(−0.015, 0.006)0.3963
PTSDForwardInverse variance weighted610.0160(−0.222, 0.254)0.8954
ReverseInverse variance weighted4−0.034(−0.083, 0.014)0.1662
SCZForwardInverse variance weighted430.152(−0.063, 0.368)0.1661
ReverseInverse variance weighted650.003(−0.017, 0.023)0.7721
TSForwardInverse variance weighted420.009(−0.488, 0.505)0.9731
ReverseInverse variance weighted16−0.015(−0.043,0.014)0.3224

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.

3.7 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).

We identified 55, 86, 34, 49, 26, 39, 51, and 16 TWAS-significant genes shared between EA and ADHD, AN, ASD, AUD, CUD, AD, BIP, or MDD, respectively (Supplementary Table S18). Most of these genes were observed in tissues from the nervous system, cardiovascular system, endocrine system, digestive system, and musculoskeletal system (Figure 4). We did not observe TWAS-significant genes shared between EA and ALS, OCD, PTSD, or TS.

Figure 4

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).

4 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 (rg = 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. 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.31 loci.

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 affects synaptic transmission and plasticity in vivo, leading to hyperactivity and spatial memory impairment (40). Okerlund et al. (69) discovered that BSN plays a pivotal role in the local regulation of presynaptic autophagy, a process dependent on poly-ubiquitination but not on the E3 ubiquitin ligase Siah1. A recent study found that neuroinflammation leads to the induction and toxic accumulation of the synaptic protein BSN in the neuronal somata of mice and patients with multiple sclerosis, implying that BSN may contribute to neurodegeneration (70). In particular, locus 3p21.31 was also confirmed to be shared for EA, ADHD, AN, BIP, CUD, and MDD using TWAS.

Moreover, CHRH1 and WNT3 have been shown to have strong associations with AUD (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 (ORIVW = 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.8 months), 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.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

DC: Conceptualization, Data curation, Formal analysis, Investigation, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing. YiZ: Conceptualization, Data curation, Investigation, Resources, Software, Supervision, Writing – original draft, Writing – review & editing. YaZ: Formal analysis, Investigation, Methodology, Project administration, Writing – review & editing. HZ: Conceptualization, Data curation, Investigation, Writing – review & editing. LW: Conceptualization, Data curation, Investigation, Supervision, Writing – review & editing. YL: Conceptualization, Data curation, Investigation, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by Sanming Project of Medicine in Shenzhen [No. SZSM202111001].

Acknowledgments

The authors thank all the participants and researchers who contributed and collected data.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2023.1303430/full#supplementary-material

    Glossary

  • EA

    Educational attainment

  • NPDs

    Neuropsychiatric disorders

  • ADHD

    Attention-deficit/hyperactivity disorder

  • ALS

    Amyotrophic lateral sclerosis

  • AD

    Alzheimer’s dementia

  • AN

    Anorexia nervosa

  • ASD

    Autism spectrum disorder

  • AUD

    Alcohol use disorder

  • BIP

    Bipolar disorder

  • CUD

    Cannabis use disorder

  • MDD

    Major depressive disorder

  • OCD

    Obsessive-compulsive disorder

  • PTSD

    Posttraumatic stress disorder

  • SCZ

    Schizophrenia

  • TS

    Tourette’s syndrome

  • HDL

    High-definition likelihood

  • GWAS

    Genome-wide association study

  • LDSC

    Linkage disequilibrium score regression

  • CPASSOC

    Cross phenotype association study

  • TSEA

    Tissue-specific enrichment analysis

  • MR

    Mendelian randomization

  • TWAS

    Transcriptome-wide association study

  • SNP

    Single nucleotide polymorphism

  • GTEx

    Genotype-tissue expression portal

References

  • 1.

    BonaccioMdi CastelnuovoACostanzoSPersichilloMDonatiMBde GaetanoGet al. Interaction between education and income on the risk of all-cause mortality: prospective results from the MOLI-SANI study. Int J Public Health. (2016) 61:76576. doi: 10.1007/s00038-016-0822-z

  • 2.

    BraniganARMcCallumKJFreeseJ. Variation in the heritability of educational attainment: an international Meta-analysis. Soc Forces. (2013) 92:10940. doi: 10.1093/sf/sot076

  • 3.

    HymanSE. A glimmer of light for neuropsychiatric disorders. Nature. (2008) 455:8903. doi: 10.1038/nature07454

  • 4.

    VoigtRGKatusicSKColliganRCKillianJMWeaverALBarbaresiWJ. Academic achievement in adults with a history of childhood attention-deficit/hyperactivity disorder: a population-based prospective Study. J Dev Behav Pediatr. (2017) 38:111. doi: 10.1097/DBP.0000000000000358

  • 5.

    XuWTanLWangHFJiangTTanMSTanLet al. Meta-analysis of modifiable risk factors for Alzheimer's disease. J Neurol Neurosurg Psychiatry. (2015) 86:1299306. doi: 10.1136/jnnp-2015-310548

  • 6.

    SutedjaNAVeldinkJHFischerKKromhoutHWokkeJHJHuismanMHBet al. Lifetime occupation, education, smoking, and risk of ALS. Neurology. (2007) 69:150814. doi: 10.1212/01.wnl.0000277463.87361.8c

  • 7.

    GrantJDScherrerJFLynskeyMTAgrawalADuncanAEHaberJRet al. Associations of alcohol, nicotine, cannabis, and drug use/dependence with educational attainment: evidence from cotwin-control analyses. Alcohol Clin Exp Res. (2012) 36:141220. doi: 10.1111/j.1530-0277.2012.01752.x

  • 8.

    BreslauJLaneMSampsonNKesslerRC. Mental disorders and subsequent educational attainment in a US national sample. J Psychiatr Res. (2008) 42:70816. doi: 10.1016/j.jpsychires.2008.01.016

  • 9.

    Pérez-VigilAFernández de la CruzLBranderGIsomuraKJangmoAFeldmanIet al. Association of Obsessive-Compulsive Disorder with Objective Indicators of educational attainment: a Nationwide register-based sibling control Study. JAMA Psychiatry. (2018) 75:4755. doi: 10.1001/jamapsychiatry.2017.3523

  • 10.

    FrankenbergESikokiBSumantriCSuriastiniWThomasD. Education, vulnerability, and resilience after a natural disaster. Ecol Soc. (2013) 18:16. doi: 10.5751/ES-05377-180216

  • 11.

    Pérez-VigilAFernández de la CruzLBranderGIsomuraKJangmoAKuja-HalkolaRet al. Association of Tourette Syndrome and Chronic tic Disorders with Objective Indicators of educational attainment: a population-based sibling comparison Study. JAMA Neurol. (2018) 75:1098105. doi: 10.1001/jamaneurol.2018.1194

  • 12.

    KeenDWebsterARidleyG. How well are children with autism spectrum disorder doing academically at school? An overview of the literature. Autism. (2016) 20:27694. doi: 10.1177/1362361315580962

  • 13.

    SundquistJOhlssonHWinklebyMASundquistKCrumpC. School achievement and risk of eating disorders in a Swedish National Cohort. J Am Acad Child Adolesc Psychiatry. (2016) 55:4146.e1. doi: 10.1016/j.jaac.2015.09.021

  • 14.

    AndersonELHoweLDWadeKHBen-ShlomoYHillWDDearyIJet al. Education, intelligence and Alzheimer's disease: evidence from a multivariable two-sample Mendelian randomization study. Int J Epidemiol. (2020) 49:116372. doi: 10.1093/ije/dyz280

  • 15.

    DardaniCRiglinLLeppertBSandersonERaiDHoweLDet al. Is genetic liability to ADHD and ASD causally linked to educational attainment?Int J Epidemiol. (2022) 50:201123. doi: 10.1093/ije/dyab107

  • 16.

    ZhangLTangLXiaKHuangTFanD. Education, intelligence, and amyotrophic lateral sclerosis: a Mendelian randomization study. Ann Clin Transl Neurol. (2020) 7:16427. doi: 10.1002/acn3.51156

  • 17.

    ChenDWangXHuangTJiaJ. Genetic support of a causal relationship between cannabis use and educational attainment: a two-sample Mendelian randomization study of European ancestry. Addiction. (2023) 118:698710. doi: 10.1111/add.16090

  • 18.

    OkbayABeauchampJPFontanaMALeeJJPersTHRietveldCAet al. Genome-wide association study identifies 74 loci associated with educational attainment. Nature. (2016) 533:53942. doi: 10.1038/nature17671

  • 19.

    DemontisDWaltersRKMartinJMattheisenMAlsTDAgerboEet al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet. (2019) 51:6375. doi: 10.1038/s41588-018-0269-7

  • 20.

    WatsonHJYilmazZThorntonLMHübelCColemanJRGasparHAet al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. Nat Genet. (2019) 51:120714. doi: 10.1038/s41588-019-0439-2

  • 21.

    StahlEABreenGForstnerAJMcQuillinARipkeSTrubetskoyVet al. Genome-wide association study identifies 30 loci associated with bipolar disorder. Nat Genet. (2019) 51:793803. doi: 10.1038/s41588-019-0397-8

  • 22.

    WrayNRRipkeSMattheisenMTrzaskowskiMByrneEMAbdellaouiAet al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. (2018) 50:66881. doi: 10.1038/s41588-018-0090-3

  • 23.

    JohnsonECDemontisDThorgeirssonTEWaltersRKPolimantiRHatoumASet al. A large-scale genome-wide association study meta-analysis of cannabis use disorder. Lancet Psychiatry. (2020) 7:103245. doi: 10.1016/S2215-0366(20)30339-4

  • 24.

    GroveJRipkeSAlsTDMattheisenMWaltersRKWonHet al. Identification of common genetic risk variants for autism spectrum disorder. Nat Genet. (2019) 51:43144. doi: 10.1038/s41588-019-0344-8

  • 25.

    Sanchez-RoigeSPalmerAAFontanillasPElsonSLthe 23andMe Research Team, the Substance Use Disorder Working Group of the Psychiatric Genomics ConsortiumAdamsMJet al. Genome-wide association Study Meta-analysis of the alcohol use disorders identification test (AUDIT) in two population-based cohorts. Am J Psychiatry. (2019) 176:10718. doi: 10.1176/appi.ajp.2018.18040369

  • 26.

    NingZPawitanYShenX. High-definition likelihood inference of genetic correlations across human complex traits. Nat Genet. (2020) 52:85964. doi: 10.1038/s41588-020-0653-y

  • 27.

    FinucaneHKBulik-SullivanBGusevATrynkaGReshefYLohPRet al. Partitioning heritability by functional annotation using genome-wide association summary statistics. Nat Genet. (2015) 47:122835. doi: 10.1038/ng.3404

  • 28.

    ZhuXFengTTayoBOLiangJYoungJHFranceschiniNet al. Meta-analysis of correlated traits via summary statistics from GWASs with an application in hypertension. Am J Hum Genet. (2015) 96:2136. doi: 10.1016/j.ajhg.2014.11.011

  • 29.

    HemaniGTillingKDavey SmithG. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. (2017) 13:e1007081. doi: 10.1371/journal.pgen.1007081

  • 30.

    VerbanckMChenCYNealeBDoR. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:6938. doi: 10.1038/s41588-018-0099-7

  • 31.

    GusevAKoAShiHBhatiaGChungWPenninxBWJHet al. Integrative approaches for large-scale transcriptome-wide association studies. Nat Genet. (2016) 48:24552. doi: 10.1038/ng.3506

  • 32.

    FarhKKMarsonAZhuJKleinewietfeldMHousleyWJBeikSet al. Genetic and epigenetic fine mapping of causal autoimmune disease variants. Nature. (2015) 518:33743. doi: 10.1038/nature13835

  • 33.

    GiambartolomeiCVukcevicDSchadtEEFrankeLHingoraniADWallaceCet al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. (2014) 10:e1004383. doi: 10.1371/journal.pgen.1004383

  • 34.

    PeiGDaiYZhaoZJiaP. deTS: tissue-specific enrichment analysis to decode tissue specificity. Bioinformatics. (2019) 35:38425. doi: 10.1093/bioinformatics/btz138

  • 35.

    BindeaGMlecnikBHacklHCharoentongPTosoliniMKirilovskyAet al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. (2009) 25:10913. doi: 10.1093/bioinformatics/btp101

  • 36.

    ShannonPMarkielAOzierOBaligaNSWangJTRamageDet al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498504. doi: 10.1101/gr.1239303

  • 37.

    SheenVLGaneshVSTopcuMSebireGBodellAHillRSet al. Mutations in ARFGEF2 implicate vesicle trafficking in neural progenitor proliferation and migration in the human cerebral cortex. Nat Genet. (2004) 36:6976. doi: 10.1038/ng1276

  • 38.

    de WitMCde CooIFHalleyDJLequinMHManciniGM. Movement disorder and neuronal migration disorder due to ARFGEF2 mutation. Neurogenetics. (2009) 10:3336. doi: 10.1007/s10048-009-0192-2

  • 39.

    DaviesGArmstrongNBisJCBresslerJChourakiVGiddaluruSet al. Genetic contributions to variation in general cognitive function: a meta-analysis of genome-wide association studies in the CHARGE consortium (N=53949). Mol Psychiatry. (2015) 20:18392. doi: 10.1038/mp.2014.188

  • 40.

    LiangZZhanYShenYWongCCLYatesJRIIIPlattnerFet al. The pseudokinase CaMKv is required for the activity-dependent maintenance of dendritic spines. Nat Commun. (2016) 7:13282. doi: 10.1038/ncomms13282

  • 41.

    AltrockWDtom DieckSSokolovMMeyerACSiglerABrakebuschCet al. Functional inactivation of a fraction of excitatory synapses in mice deficient for the active zone protein bassoon. Neuron. (2003) 37:787800. doi: 10.1016/S0896-6273(03)00088-6

  • 42.

    DickOtom DieckSAltrockWDAmmermüllerJWeilerRGarnerCCet al. The presynaptic active zone protein bassoon is essential for photoreceptor ribbon synapse formation in the retina. Neuron. (2003) 37:77586. doi: 10.1016/S0896-6273(03)00086-2

  • 43.

    Haldeman-EnglertCRGaiXPerinJCCianoMHalbachSSGeigerEAet al. A 3.1-Mb microdeletion of 3p21.31 associated with cortical blindness, cleft lip, CNS abnormalities, and developmental delay. Eur J Med Genet. (2009) 52:2658. doi: 10.1016/j.ejmg.2008.11.005

  • 44.

    ZorrillaEPLogripMLKoobGF. Corticotropin releasing factor: a key role in the neurobiology of addiction. Front Neuroendocrinol. (2014) 35:23444. doi: 10.1016/j.yfrne.2014.01.001

  • 45.

    SchmittAMShiJWolfAMLuCCKingLAZouY. Wnt-Ryk signalling mediates medial-lateral retinotectal topographic mapping. Nature. (2006) 439:317. doi: 10.1038/nature04334

  • 46.

    DuncanLYilmazZGasparHWaltersRGoldsteinJAnttilaVet al. Significant locus and metabolic genetic correlations revealed in genome-wide association Study of anorexia nervosa. Am J Psychiatry. (2017) 174:8508. doi: 10.1176/appi.ajp.2017.16121402

  • 47.

    WiameETahayGTytecaDVertommenDStroobantVBommerGTet al. NAT6 acetylates the N-terminus of different forms of actin. FEBS J. (2018) 285:3299316. doi: 10.1111/febs.14605

  • 48.

    ChongSASubramaniamMLeeIMPekECheokCVermaSet al. Academic attainment: a predictor of psychiatric disorders?Soc Psychiatry Psychiatr Epidemiol. (2009) 44:9991004. doi: 10.1007/s00127-009-0027-3

  • 49.

    MadiganDJJEPR. A Meta-analysis of perfectionism and academic achievement. Educ Psychol Rev. (2019) 31:96789. doi: 10.1007/s10648-019-09484-2

  • 50.

    NilssonKSundbomEHägglöfB. A longitudinal study of perfectionism in adolescent onset anorexia nervosa-restricting type. Eur Eat Disord Rev. (2008) 16:38694. doi: 10.1002/erv.850

  • 51.

    DupuisAMudiyanselagePBurtonCLArnoldPDCrosbieJSchacharRJ. Hyperfocus or flow? Attentional strengths in autism spectrum disorder. Front Psych. (2022) 13:886692. doi: 10.3389/fpsyt.2022.886692

  • 52.

    MackinnonSPRayCMFirthSMO'ConnorRM. Perfectionism, negative motives for drinking, and alcohol-related problems: a 21-day diary study. J Res Pers. (2019) 78:17788. doi: 10.1016/j.jrp.2018.12.003

  • 53.

    CorryJGreenMRobertsGFranklandAWrightALauPet al. Anxiety, stress and perfectionism in bipolar disorder. J Affect Disord. (2013) 151:101624. doi: 10.1016/j.jad.2013.08.029

  • 54.

    WuKDCortesiGT. Relations between perfectionism and obsessive-compulsive symptoms: examination of specificity among the dimensions. J Anxiety Disord. (2009) 23:393400. doi: 10.1016/j.janxdis.2008.11.006

  • 55.

    Escott-PriceVBracher-SmithMMenziesGWaltersJKirovGOwenMJet al. Genetic liability to schizophrenia is negatively associated with educational attainment in UK biobank. Mol Psychiatry. (2020) 25:7035. doi: 10.1038/s41380-018-0328-6

  • 56.

    BansalVMitjansMBurikCAPLinnérRKOkbayARietveldCAet al. Genome-wide association study results for educational attainment aid in identifying genetic heterogeneity of schizophrenia. Nat Commun. (2018) 9:3078. doi: 10.1038/s41467-018-05510-z

  • 57.

    SheenVLJansenAChenMHParriniEMorganTRavenscroftRet al. Filamin a mutations cause periventricular heterotopia with Ehlers-Danlos syndrome. Neurology. (2005) 64:25462. doi: 10.1212/01.WNL.0000149512.79621.DF

  • 58.

    HongEHKimJYKimJHLimDSKimMKimJY. BIG2-ARF1-RhoA-mDia1 signaling regulates dendritic Golgi polarization in hippocampal neurons. Mol Neurobiol. (2018) 55:770116. doi: 10.1007/s12035-018-0954-7

  • 59.

    LiCCLeKKatoJMossJVaughanM. Enhancement of β-catenin activity by BIG1 plus BIG2 via Arf activation and cAMP signals. Proc Natl Acad Sci U S A. (2016) 113:594651. doi: 10.1073/pnas.1601918113

  • 60.

    AndrewsSJDasDAnsteyKJEastealS. Association of AKAP6 and MIR2113 with cognitive performance in a population-based sample of older adults. Genes Brain Behav. (2017) 16:4728. doi: 10.1111/gbb.12368

  • 61.

    LauPBossersKJankyRSaltaEFrigerioCSBarbashSet al. Alteration of the microRNA network during the progression of Alzheimer's disease. EMBO Mol Med. (2013) 5:161334. doi: 10.1002/emmm.201201974

  • 62.

    DaviesGMarioniRELiewaldDCHillWDHagenaarsSPHarrisSEet al. Genome-wide association study of cognitive functions and educational attainment in UK biobank (N=112 151). Mol Psychiatry. (2016) 21:75867. doi: 10.1038/mp.2016.45

  • 63.

    KranzlerHRZhouHKemberRLVickers SmithRJusticeACDamrauerSet al. Genome-wide association study of alcohol consumption and use disorder in 274,424 individuals from multiple populations. Nat Commun. (2019) 10:1499. doi: 10.1038/s41467-019-09480-8

  • 64.

    HouLBergenSEAkulaNSongJHultmanCMLandénMet al. Genome-wide association study of 40,000 individuals identifies two novel loci associated with bipolar disorder. Hum Mol Genet. (2016) 25:338394. doi: 10.1093/hmg/ddw181

  • 65.

    YuDSulJHTsetsosFNawazMSHuangAYZelayaIet al. Interrogating the genetic determinants of Tourette's syndrome and other tic disorders through genome-wide association studies. Am J Psychiatry. (2019) 176:21727. doi: 10.1176/appi.ajp.2018.18070857

  • 66.

    KasherPRSchertzKEThomasMJacksonAAnnunziataSBallesta-MartinezMJet al. Small 6q16.1 deletions encompassing POU3F2 cause susceptibility to obesity and variable developmental delay with intellectual disability. Am J Hum Genet. (2016) 98:36372. doi: 10.1016/j.ajhg.2015.12.014

  • 67.

    LamMTrampushJWYuJKnowlesEDaviesGLiewaldDCet al. Large-scale cognitive GWAS Meta-analysis reveals tissue-specific neural expression and potential nootropic drug targets. Cell Rep. (2017) 21:2597613. doi: 10.1016/j.celrep.2017.11.028

  • 68.

    CajigasIJTushevGWillTJtom DieckSFuerstNSchumanEM. The local transcriptome in the synaptic neuropil revealed by deep sequencing and high-resolution imaging. Neuron. (2012) 74:45366. doi: 10.1016/j.neuron.2012.02.036

  • 69.

    OkerlundNDSchneiderKLeal-OrtizSMontenegro-VenegasCKimSAGarnerLCet al. Bassoon controls presynaptic autophagy through Atg5. Neuron. (2017) 93:897913.e7. doi: 10.1016/j.neuron.2017.01.026

  • 70.

    SchattlingBEnglerJBVolkmannCRothammerNWooMSPetersenMet al. Bassoon proteinopathy drives neurodegeneration in multiple sclerosis. Nat Neurosci. (2019) 22:88796. doi: 10.1038/s41593-019-0385-4

  • 71.

    GregorAKrumbiegelMKrausCReisAZweierC. De novo triplication of the MAPT gene from the recurrent 17q21.31 microdeletion region in a patient with moderate intellectual disability and various minor anomalies. Am J Med Genet A. (2012) 158A:176570. doi: 10.1002/ajmg.a.35427

  • 72.

    WibergANgMal OmranYAlfaro-AlmagroFMcCarthyPMarchiniJet al. Handedness, language areas and neuropsychiatric diseases: insights from brain imaging and genetics. Brain. (2019) 142:293847. doi: 10.1093/brain/awz257

  • 73.

    SachsKVHarnkeBMehlerPSKrantzMJ. Cardiovascular complications of anorexia nervosa: a systematic review. Int J Eat Disord. (2016) 49:23848. doi: 10.1002/eat.22481

  • 74.

    EliaJVetterVL. Cardiovascular effects of medications for the treatment of attention-deficit hyperactivity disorder: what is known and how should it influence prescribing in children?Paediatr Drugs. (2010) 12:16575. doi: 10.2165/11532570-000000000-00000

  • 75.

    CireraLTormoMJChirlaqueMDNavarroC. Cardiovascular risk factors and educational attainment in southern Spain: a study of a random sample of 3091 adults. Eur J Epidemiol. (1998) 14:75563. doi: 10.1023/A:1007596222217

  • 76.

    ShadrinAASmelandOBZayatsTSchorkAJFreiOBettellaFet al. Novel loci associated with attention-deficit/hyperactivity disorder are revealed by leveraging polygenic overlap with educational attainment. J Am Acad Child Adolesc Psychiatry. (2018) 57:8695. doi: 10.1016/j.jaac.2017.11.013

  • 77.

    Soler ArtigasMSánchez-MoraCRoviraPVilar-RibóLRamos-QuirogaJARibasésM. Mendelian randomization analysis for attention deficit/hyperactivity disorder: studying a broad range of exposures and outcomes. Int J Epidemiol. (2023) 52:386402. doi: 10.1093/ije/dyac128

  • 78.

    McDougalERibyDMHanleyM. Profiles of academic achievement and attention in children with and without autism Spectrum disorder. Res Dev Disabil. (2020) 106:103749. doi: 10.1016/j.ridd.2020.103749

  • 79.

    DalsgaardSMcGrathJØstergaardSDWrayNRPedersenCBMortensenPBet al. Association of Mental Disorder in childhood and adolescence with subsequent educational achievement. JAMA Psychiatry. (2020) 77:797805. doi: 10.1001/jamapsychiatry.2020.0217

  • 80.

    VerweijKJHuizinkACAgrawalAMartinNGLynskeyMT. Is the relationship between early-onset cannabis use and educational attainment causal or due to common liability?Drug Alcohol Depend. (2013) 133:5806. doi: 10.1016/j.drugalcdep.2013.07.034

  • 81.

    EricksonJel-GabalawyRPalitskyDPattenSMackenzieCSSteinMBet al. Educational attainment as a protective factor for psychiatric disorders: findings from a nationally representative longitudinal study. Depress Anxiety. (2016) 33:101322. doi: 10.1002/da.22515

  • 82.

    BrumptonBSandersonEHeilbronKHartwigFPHarrisonSVieet al. Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses. Nat Commun. (2020) 11:3519. doi: 10.1038/s41467-020-17117-4

Summary

Keywords

educational attainment, neuropsychiatric disorders, genetic correlation, shared genetic architecture, Mendelian randomization

Citation

Chen D, Zhou Y, Zhang Y, Zeng H, Wu L and Liu Y (2024) Unraveling shared susceptibility loci and Mendelian genetic associations linking educational attainment with multiple neuropsychiatric disorders. Front. Psychiatry 14:1303430. doi: 10.3389/fpsyt.2023.1303430

Received

28 September 2023

Accepted

11 December 2023

Published

04 January 2024

Volume

14 - 2023

Edited by

Panagiotis Ferentinos, National and Kapodistrian University of Athens, Greece

Reviewed by

Zhongshan Cheng, St. Jude Children's Research Hospital, United States; Vasilios Karageorgiou, University of Exeter, United Kingdom

Updates

Copyright

*Correspondence: Yi Zhou, Dongze Chen,

†These authors have contributed equally to this work and share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics