Educational attainment and mental health conditions: a within-sibship Mendelian randomization study
===================================================================================================

* María Fernanda Vinueza Veloz
* Laxmi Bhatta
* Paul Remy Jones
* Martin Tesli
* George Davey Smith
* Neil Martin Davies
* Ben M. Brumpton
* Øyvind Erik Næss

## Abstract

**Importance** Observational studies have demonstrated consistent protective effects of higher educational attainment (EA) on the risk of suffering mental health conditions (MHC). Determining whether these beneficial effects are causal is challenging given the potential role of dynastic effects and demographic factors (assortative mating and population structure) in this association.

**Objective** To evaluate to what extent the relationship between EA and various MHC is independent from dynastic effects and demographic factors.

**Design** Within-sibship Mendelian randomization (MR) study.

**Setting** One-sample MR analyses included participants’ data from the Trøndelag Health Study (HUNT, Norway) and UK Biobank (United Kingdom). For two-sample MR analyses we used summary statistics from publicly available genome-wide-association-studies.

**Participants** 61 880 siblings (27 507 sibships).

**Exposure** Years of education.

**Main outcomes** Scores for symptoms of anxiety, depression and neuroticism using the Hospital Anxiety Depression Scale (HADS), the 7-item Generalized Anxiety Disorder Scale (GAD-7), the 9-item Patient Health Questionnaire (PHQ-9), and the Eysenck Personality Questionnaire, as well as self-reported consumption of psychotropic medication.

**Results** One standard deviation (SD) unit increase in years of education was associated with a lower symptom score of anxiety (−0.20 SD [95%CI: −0.26, −0.14]), depression (−0.11 SD [−0.43, 0.22]), neuroticism (−0.30 SD [−0.53, −0.06]), and lower odds of psychotropic medication consumption (OR: 0.60 [0.52, 0.69]). Estimates from the within-sibship MR analyses showed some attenuation, which however were suggestive of a causal association (anxiety: −0.17 SD [−0.33, −0.00]; depression: −0.18 SD [-1.26, 0.89]; neuroticism: −0.29 SD [−0.43, −0.15]); psychotropic medication consumption: OR, 0.52 [0.34, 0.82]).

**Conclusions and Relevance** Associations between EA and MHC in adulthood, although to some extend explained by dynastic effects and demographic factors, overall remain robust, indicative of a causal effect. However, larger studies are warranted to improve statistical power and further validate our conclusions.

**Question** Are the protective effects of higher educational attainment on the risk of suffering mental health conditions independent from familial (i.e.,dynastic effects) and demographic factors (i.e., assortative mating and population structure)?

**Findings** Although to some extend explained by familial and demographic factors, overall the protective effects of higher educational attainment remain robust.

**Meaning** Our findings are consistent with those of previous observational studies, which suggests that higher educational attainment may causally reduce the risk of suffering a mental health condition.

Key words
*   educational attainment
*   Mendelian randomization
*   depression
*   anxiety
*   neuroticism
*   psychotropic medication

## Introduction

Mental health conditions (MHC) are mental disorders, psycho-social disabilities as well as mental states associated with significant distress, impaired functioning, or risk of self-harm 1. Common MHC including anxiety and mood disorders (e.g., depression) are among the top causes of years lived with disability across all age groups 2. Moreover, people with MHC have a mortality rate that is at least twice higher than that of the general population or people without such conditions 3. Socioeconomic circumstances including those related to level of educational attainment (EA) have been related to the risk of suffering MHC 4.

Higher EA is associated with a lower risk of suffering MHC including, anxiety, depression, personality disorders, and substance abuse 5–7. This inverse relationship might be explained by the wider availability of mental/physical resources in adulthood to cope with daily hassles among educated individuals, which could assist preventing MHC 8. Nevertheless, it has been shown that family-related factors influence both educational level and mental health status, which suggests that alternative causal pathways in adulthood may have a role 9–12.

Mendelian randomization (MR) is a method based on instrumental variable analysis used to investigate relationships between exposures and outcomes by using genetic variants or polygenic scores (PGS) as instruments for the modifiable factor 13. MR relies on the premise that there is a causal pathway from an individual’s genotype to the individual’s phenotype (e.g., a genotype-phenotype association). If its three core assumptions hold, the approach can overcome both exposure-outcome confounding and reverse causation.

To provide evidence of causality i) the genetic instrument must be associated to the exposure of interest and this relationship is required to be reasonably strong (relevance); ii) there mustn’t be common factors associated to the genetic instrument or the outcome of interest (independence); iii) the genetic instrument must influence the outcome only through the exposure of interest, which implies that the effect of the genetic instrument on an outcome is fully mediated by the exposure (exclusion) 13.

Two recent population-based MR studies have provided evidence that supports the role of EA in the aetiology of MHC like anxiety and depression in adulthood 14,15. However, since the phenotypic variation of complex traits cannot be fully attributed to direct genetic effects, such MR estimates might be biased 16,17. In fact, violations of the independence assumption resulting from indirect genetic effects (i.e., dynastic effects) and confounding due demographic factors (i.e., assortative mating and population stratification) are now well documented for various traits including education 16–18.

Dynastic effects are the result of genetic nurturing, namely the effect of parents’ or other relatives’ genotypes on an individual’s phenotype through their contributions to the environment 17,18. Assortative mating is when partners are not selected at random but based on particular characteristics 19. Population structure is the presence of subpopulation differences in allele frequencies that correlate with both the phenotype and the genotype 16. In the context of MR all three mechanisms have the potential of amplifying the effect of a given genotype-phenotype association 17,20.

The objective of the present study is to evaluate if the association between EA and symptoms of anxiety, depression, and neuroticism as well as the consumption of psychotropic medication are explained by indirect genetic effects and/or are confounded by demographic factors. We applied a within-sibship design coupled with genomic data to perform MR to account for dynastic effects, assortative mating, and population structure.

### Methodology

#### Study design and data sources

The present is the report of a within-sibship MR study. One-sample and two-sample MR methods were applied using individual-level data, and summary statistics from genome wide association studies (GWAS). Individual-level data came from the Trøndelag Health Study (HUNT) and UK Biobank (UKB) 21,22. A brief description of contributing GWAS can be found in eTable 1.

#### Setting and participants

HUNT is a population-based cohort study that is held in the Trøndelag County in Norway and started in 1984 21. We used data from all participants of the second (HUNT2) and third HUNT wave (HUNT3) who had been genotyped. From those, we selected all individuals who were > 30 years of age when they participated in the survey, and had at least one sibling. The final sample from HUNT2 and HUNT3 included 26 770 (10 428 sibships) and 16 718 siblings (7010 sibships), respectively (eTable 2).

The UKB is a prospective cohort study that began in 2006. UKB is following nearly 500 000, 40 to 69 years old participants from across the UK, who volunteered to be part of the study and provided consent for follow-up through linkage to their health records 22. We included all participants who had been genotyped. After restricting the sample to sibships with two or more individuals, our analysis sample included 35 118 participants from 17 079 sibships (eTable 2).

#### Genetic variants

For HUNT participants we used a weighted PGS as an instrumental variable for EA (PGS-edu). The calculation of the PGS-edu was based on the genetic variants reported as significantly associated to years of education at the genome-wide level (*p* < 5 × 10-8) in a recent GWAS 23. From the 3952 independent genetic variants reported by Okbay et al., we included those that were well imputed in the target population (eTable 3) 24. For further information see eMethods.

#### Exposure and outcomes

EA was the exposure of interest and was assessed trough the question: What is your highest level of education? for HUNT participants, or Which of the following qualification do you have? for UKB participants. A number of years of attained education was assigned for each of the answers based on The International Standard Classification of Education (ISCED) mapping 1997 (eTable 4).

Symptoms of anxiety, depression and neuroticism as well as consumption of psychotropic medication were the outcomes of interest. In HUNT symptoms of anxiety and depression were assessed by the Hospital Anxiety and Depression Scale (HADS), and in UKB by the 7-item Generalized Anxiety Disorder Scale (GAD-7) and 9-item Patient Health Questionnaire (PHQ-9) 25–27. In both HUNT and UKB, neuroticism was measured using the Eysenck Personality Questionnaire and consumption of psychotropic medication was self-reported 28. A detailed description on how the exposure and outcomes were processed can be found in eMethods.

#### Ethics approval and informed consent

See eMethods.

#### Statistical analysis

Before running the analyses using individual-level data we standardized all numerical variables so that they had a mean of 0 and standard deviation (SD) of 1. EA, as well as symptoms of anxiety, depression and neuroticism were analysed as continuos, while consumption of psychotropic medication as categorical (yes/no). The association between EA and the outcomes was analysed by applying one-sample MR (two-stage least squares regression), and ordinary least squares (OLS) or logistic (LOG) regression for comparison.

All models were adjusted by sex and age. However, when the PGS-edu was included as predictor, the model was also adjusted by the first 10 principal components of ancestry (PCA) and genotyping batch to account for population structure. In all cases we assumed that standard errors were correlated within sibships and therefore clustered standard errors were computed. OLS, LOG and one-sample MR analyses were performed using the “feols” (for continuos) and “felgm” (for categorical) functions of the “fixest” package in R 29,30.

Any difference between families due to indirect genetic effects (i.e. dynastic effects) or bias due to assortative mating was accounted for by using family fixed effects (anxiety, depression and neuroticism) or the sibling difference method (consumption of psychotropic medication) 16. All analyses were performed using the “fixest” package in R (see eMethods for a detailed description)30.

Individual-level data for each of the two cohorts were analysed separately using the same model specifications and R packages. Then, results were meta-analysed using the “rma” function from the “metafor” R package 31. We applied fixed effect models, except when heterogeneity between HUNT and UK Biobank estimates was detected, that is Cochran’s Q Chi2 < 0.05 and I2 > 50% (eTable 5).

Two-sample MR analyses were performed using summary statistics from GWAS described in eTable 1 and the R package “TwoSampleMR” 32. The inverse-variance estimator weighted (IVW) estimator and its 95% confidence intervals (CI) are reported in the present work. Robust estimators including weighted median, weighted mode, and MR-Egger coefficients were calculated to investigate pleiotropy and reported in supplementary tables. The directionality of the effect was evaluated using Steiger test of directionality. MR-Egger intercept test was performed to assess pleiotropy. Within-sibship two-sample MR was performed for depressive symptoms and neuroticism, as only for them we had access to summary statistics from within-sibship GWAS (eTable 1) 33.

#### Handling of missing information

We imputed missing data for various questions of the HADS score as well as education and consumption of psychotropic medication in HUNT (see eTable 6 for details on missing information). We applied multivariate imputation to each HUNT survey, using fully conditional specification implemented by the MICE algorithm 34. For each imputed data set, we then calculated the corresponding score. For further information see eMethods.

## Results

### Descriptive statistics

In both cohorts participation rate among females was higher than that among males. Participants from HUNT were younger, but had similar level of education than those from the UKB (Table 1). Symptoms of anxiety were more often reported by HUNT than from UKB participants (8% vs. 3%). About 5% and 9% of HUNT and UKB participants reported symptoms of depression and consumption of psychotropic medication, respectively (Table 1). However, the consumption of psychotropic medication was slightly higher among HUNT than UKB participants. Differences between females and males can be seen in Table 1. For each of the cohorts, the total number of siblings varied depending on the outcome (eTable 2).

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/08/11/2024.08.10.24311789/T1)

Table 1. General characteristics of HUNT and UKB samples
. In HUNT symptoms of anxiety and depression were assessed by the Hospital Anxiety and Depression Scale (HADS), and in UK Biobank by the 7-item Generalized Anxiety Disorder Scale (GAD-7) and 9-item Patient Health Questionnaire (PHQ-9), respectively. Neuroticism was assessed using a six-item and a 12-item Eysenck Personality Questionnaire in HUNT and UK Biobank, respectively (see eMethods).

The PGS-edu was associated with years of education, conditional from age, sex, first 10 PCA, and batch in both HUNT and UK Biobank. In HUNT, each SD unit increase of the PGS-edu was associated with a 0.19 SD increase in years of education (95%CI, 0.18: 0.20, *p =* 4.40×10-187, *F-test stat.* = 162.99, *r2* = 0.03). In UK Biobank, each SD unit increase of the PGS-edu was associated with a 0.24 SD increase in years of education (95%CI, 0.23: 0.24, *p* = 2.20×10-16, *F-test stat.* = 103.90, *r2 =* 0.06). This association was attenuated after including a family fixed effect (HUNT: 0.13, 95%CI: 0.11 – 0.15, *p =* 3.47×10-44, *F-test stat.* = 1349.62, *r2* = 0.64; UK Biobank: 0.13, 95%CI, 0.11: 0.14, *p* = 3.38×10-43, *F-test stat.* = 2362.40, *r2* = 0.33). The associations between the PGS ant the outcomes are depicted in eTable 7.

### Main analyses

The direction of the regression and MR estimates was consistent among both cohorts and across all analyses. There were however some differences in the strength of the associations between HUNT and UKB (Table 2). These differences were more evident for the depression and neuroticism MR estimates.

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/08/11/2024.08.10.24311789/T2)

Table 2. Association between educational attainment and mental health outcomes
. HUNT and UK Biobank data were analysed separately and results then meta-analysed (see Methodology).

The results of the MR analyses suggest that EA has a beneficial effect on anxiety, depression and neuroticism symptomatology (Figure 1, 3 and 4). A higher EA was also associated with lower consumption of psychotropic medication (Table 2). A SD unit increase in years of education was associated with lower scores of anxiety (−0.20 SD [−0.26, −0.14]), and neuroticism (−0.30 SD [−0.53, −0.06]) as well as lower odds of psychotropic medication consumption (OR: 0.60 [0.52, 0.69]). There was no association between EA and depression (−0.11 SD [−0.43, 0.22]) (Figure 2). Estimates from one-sample MR were larger than those of the regression analyses, except from depression, but were consistent in direction between each other (Table 2).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/11/2024.08.10.24311789/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/08/11/2024.08.10.24311789/F1)

Figure 1. Educational attainment and symptoms of anxiety.
Standard deviation (SD) changes in the anxiety score and its 95% confidence interval per SD increase in years of education are shown. Estimated associations are displayed for ordinary least square regression (OLS) and Mendelian randomization models. *Abbreviations: SD, standard deviation unit; OLS EA, ordinary least square regression model with educational attainment as exposure; OLS PGS-edu, ordinary least square regression model with the educational attainment polygenic score as exposure; 1SMR, one-sample Mendelian randomization*.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/11/2024.08.10.24311789/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/08/11/2024.08.10.24311789/F2)

Figure 2. Educational attainment and symptoms of depression.
Standard deviation (SD) changes in the depression score and its 95% confidence interval per SD increase in years of education are shown. Estimated associations are displayed for ordinary least square regression (OLS) and Mendelian randomization models. *Abbreviations: SD, standard deviation unit; OLS EA, ordinary least square regression model with educational attainment as exposure; OLS PGS-edu, ordinary least square regression model with the educational attainment polygenic score as exposure; 1SMR, one-sample Mendelian randomization; 2SMR, two-sample Mendelian randomization*.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/11/2024.08.10.24311789/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/08/11/2024.08.10.24311789/F3)

Figure 3. Educational attainment and neuroticism.
Standard deviation (SD) changes in the neuroticism score and its 95% confidence interval per SD increase in years of education are shown. Estimated associations are displayed for ordinary least square regression (OLS) and Mendelian randomization models. *Abbreviations: SD, standard deviation unit; OLS EA, ordinary least square regression model with educational attainment as exposure; OLS PGS-edu, ordinary least square regression model with the educational attainment polygenic score as exposure; 1SMR, one-sample Mendelian randomization; 2SMR, two-sample Mendelian randomization*.

Most of the estimates of the within-sibling MR analyses showed some attenuation in comparison to the MR estimates, the exception was depression (Table 2 and Figure 1 to 4). The effect of EA attenuated from −0.20 to −0.17 [−0.33, −0.00] for anxiety, and from −0.30 to −0.29 [−0.43, −0.15] for neuroticism. The odds ration of psychotropic medication decreased from −0.60 to −0.52 [0.33, 0.82]. Some attenuation was also observed among the estimates of the within-sibling regression analyses, except from consumption of psychotropic medication (Table 2).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/11/2024.08.10.24311789/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/08/11/2024.08.10.24311789/F4)

Figure 4. Educational attainment and consumption of psychotropic medication.
Log odds changes in consumption of psychotropic medication and its 95% confidence interval per SD increase in years of education are shown. Estimated associations are displayed for logistic regression (LOG) and Mendelian randomization models. *Abbreviations: SD, standard deviation unit; LOG EA, logistic regression model with educational attainment as exposure; LOG PGS-edu, logistic regression model with the educational attainment polygenic score as exposure; 1SMR, one-sample Mendelian randomization*.

Two-sample MR analyses were performed for depression and neuroticism and in both cases, results were consistent with those from one-sample MR and likewise suggest a protective role of education (eTable9). Each SD increase in EA caused a −0.22 [−0.25, −0.18] and −0.21 [−0.25, −0.17] SD decrease in the scores of depression and neuroticism. Weighted median, weighted mode and MR-Egger estimators were consistent with IVW. In both cases, the directionality test indicated that the causal direction is likely to be correct, and that horizontal pleiotropy is unlikely (eTable 9).

Within-sibship two-sample MR estimates suggested an attenuation of both associations, which were however less precise, and overlapped with those of the population-based MR and the null.

## Discussion

The objective of the present study was to assess to what extent the relationships between EA and anxiety, depression, neuroticism and consumption of psychotropic medication are explained by family and/or demographic factors in adulthood. Regression and population-based MR analyses indicate that higher EA is associated with lower levels of anxiety, neuroticism, and consumption of psychotropic medication, but not depression. The protective effect of EA attenuated after accounting for family and/or demographic factors. The observed attenuation was however modest, which is indicative of a causal association between EA and the studied mental health outcomes.

Various potential mechanisms that could explain a causal relationship between EA and the propensity to develop mental health symptoms have been described, and include availability of resources, cognition, knowledge and social integration 35. Higher educated people are more likely to have higher income than less educated people, and this difference in earning can affect health by providing easier access to health care services, better nutrition and opportunities for being physical active. Similarly, EA can influence cognition, which can affect health outcomes by leading people to adopt healthier lifestyles and avoiding risk factors associated with mental health conditions such as stress.

However, as family is a driving force of adult health, the existence of other pathways than those in adulthood should be considered 9–12,36. Such pathways might involve environment-mediated genetic factors during early life (specially pertained to parents’ characteristics such as genetic-nurturing), similar but reasonably less strong genetic nurturing effect might be expected from the interaction of the individual and other relatives (dynastic-effects) 17.

EA might also associate with mental health outcomes because of parental cross-trait assortative-mating on the exposure and outcome or on variables genetically correlated with them (e.g., highly educated individuals are more likely to select partners with lower neuroticism symptomatology) 20,37. In some situations parental single-trait assortative-mating on a phenotype genetically correlated with both educational level and the mental health outcome might also explain their association (e.g., partner selection based on cognitive ability that is correlated with EA and the mental health outcome) 20,38. Nonetheless, with the present design we are unable to estimate the relative influence of dynastic effects isolated from that of assortative mating or that of population structure.

Since we observed that the level of change was dissimilar among the examined MHC, it is possible that the contribution of indirect genetic effects may vary depending on the outcome of interest and its clinical severity. Accordingly, MHC whose aetiology cannot be fully attributed to direct genetic effects may be influenced by genetic nurturing to a higher extent in comparison to others with a stronger genetic basis. For example, a less pronounced change may be expected for the association between EA and attention-deficit/hyperactivity disorder or schizophrenia, since both seem to rely strongly on direct genetic effects 39,40.

### Strengths and limitations

A major strength of our work is that we applied a robust method to evaluate the impact dynastic effects, assortative mating, and population structure may have on the association between EA and the risk of suffering a MHC 16,33. As genetic variants are randomly assigned at conception the design is generally likely to avoid bias due to confounding or reverse causality. Moreover, the application of a within-family design has allowed us to overcome the limitations of population-based MR such as non-compliance with the independence criterion 16. Importantly, our results were replicated in two different cohorts and across different analyses: Results of observational and MR models were consistent in both cohorts. Nevertheless, our approach has limitations that must be considered when interpreting our findings.

For highly polygenic phenotypes such as EA, limitations arise from the uncertainty regarding the compliance of core assumptions including relevance and exclusion criteria 36. Even though first-stage F-statistic indicated our analyses did not suffer from weak instrument bias, our results might still be underpowered to detect an effect of interest, as the PGS-edu only explained nearly 2% of the variation on EA. Similarly, we found little evidence of horizontal pleiotropy.

Importantly, the phenotypical expression of the genetic architecture of EA is contingent on social, historical, and cultural contexts, which might limit the informativeness of the genetic instruments in real scenarios 36. EA in the context of MR can be thought of as propensity to achieve a given level of education, which is environmental-dependent, and hence the impact of EA on MHC may vary depending on the environmental context 36. Moreover, MR may not be well suited to inform either novel or actionable policy recommendations, as it does not provide much insight on how or why EA is associated to any of the specific MHC we analysed 36.

Our within-family analyses had limited power, and could not rule out a causal effect EA on none of the analysed MHC in adulthood. Additionally, as MR relies on available GWAS data (e.g., to calculate PGS-edu) our study suffers from lack of representativeness of populations from other geographical contexts and ancestry which might limit the generalizability of our findings. Our results may be affected by selection bias, as participants of large studies such as UK Biobank were non-randomly sampled 41. However, simulations have shown that in the context of MR the impact of selection bias is likely to be less than that of other biases such as pleiotropy or population stratification 42.

## Conclusions

Our findings suggest that previously observed associations between EA and MHC in adulthood may to some extent originate and be explained by family-related factors (dynastic effects) and/or assortative mating, and population structure. However, overall all estimates remain robust, indicative of a causal protective effect of EA on the risk of suffering MHC. Future larger family based studies are needed to provide more precise evidence about these effects.

## Supporting information

eMethods [[supplements/311789_file06.doc]](pending:yes)

eTables [[supplements/311789_file07.doc]](pending:yes)

eTable2 [[supplements/311789_file08.xlsx]](pending:yes)

## Data Availability

Researchers associated with Norwegian research institutes can apply for the use of HUNT material: data and samples - given approval by a Regional Committee for Medical and Health Research Ethics. Researchers from other countries can also apply in cooperation with a Norwegian Principle Investigator. Information for data access can be found at [https://www.ntnu.edu/hunt/data](https://www.ntnu.edu/hunt/data). UK Biobank individual-level participant data are available via enquiry to access{at}ukbiobank.ac.uk. All GWAS summary statistics used in the present manuscript are publicly available and can be download from [https://gwas.mrcieu.ac.uk/](https://gwas.mrcieu.ac.uk/) and [https://thessgac.com/papers/](https://thessgac.com/papers/).

## Information on author access to data

MFVV and LB had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

## Authors’ contributions

Conception and design: MFVV, LB, BMB and ØEN. Analysis: MFVV and LB. Interpretation of data for the manuscript: MFVV, LB, PRJ, MT, GDS, NMD, BMB and ØEN. Drafting of the manuscript: MFVV and LB. Critical review of the manuscript for important intellectual content: MFVV, LB, PRJ, MT, GDS, NMD, BMB and ØEN. Final approval of the version to be published: MFVV, LB, PRJ, MT, GDS, NMD, BMB and ØEN. Agreement to be accountable for all aspects of the work: MFVV, LB, PRJ, MT, GDS, NMD, BMB and ØEN.

## Conflict of interest disclosures

ØEN reports having received funding from the Norwegian University of Science and Technology for a review assignment. GDS reports a Scientific Advisory Board Membership for Relation Therapeutics and Insitro. MFVV, LB, PRJ, MT, NMD and BMB declare no conflict of interest to disclose.

## Funding/Support

The K. G. Jebsen Center for Genetic Epidemiology is financed by Stiftelsen Kristian Gerhard Jebsen and the Faculty of Medicine and Health Sciences, NTNU, Norway. MFVV, PRJ and ØEN are funded by The Research Council of Norway, Project Inequalities in non-communicable diseases: Indirect selection or social causation (INDI-INEQ #287347). BMB is funded by the Research Council of Norway (FRIMEDBIO #287347). LB and BMB are supported by grants from the Liaison Committee for Education, Research and Innovation in Central Norway and the Joint Research Committee between St Olavs Hospital and the Faculty of Medicine and Health Sciences, NTNU. GDS works at the Medical Research Council (MRC) Integrative Epidemiology Unit at the University of Bristol, which is supported by MRC (MC_UU_00011/1). NMD is supported by The Research Council of Norway (#295989) and NIMH (MH130448).

## Role of the Funder

The funders had no role in the design, analysis, interpretation of the data, drafting of the manuscript, review, or final approval of the version to be published, nor in the decision to submit the manuscript for publication.

## Acknowledgments

The Trøndelag Health Study (HUNT) is a collaboration between HUNT Research Center (Faculty of Medicine and Health Sciences, NTNU, Norwegian University of Science and Technology), Trøndelag County Council, Central Norway Regional Health Authority, and the Norwegian Institute of Public Health. The genotyping in HUNT was financed by the National Institutes of Health; University of Michigan; the Research Council of Norway; the Liaison Committee for Education, Research and Innovation in Central Norway; and the Joint Research Committee between St Olavs hospital and the Faculty of Medicine and Health Sciences, NTNU. We also like to thank to Dr. David Carslake for his his valuable comments.

*   Received August 10, 2024.
*   Revision received August 10, 2024.
*   Accepted August 11, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.World Health Organization. World Mental Health Report: Transforming Mental Health for All. World Health Organization; 2022. Accessed July 8, 2022. [https://www.who.int/publications/i/item/9789240049338](https://www.who.int/publications/i/item/9789240049338)
    
    

2.  2.Abbafati C, Abbas KM, Abbasi-Kangevari M, et al. Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396(10258):1222. doi:10.1016/S0140-6736(20)30925-9
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30925-9&link_type=DOI) 

3.  3.Walker ER, McGee RE, Druss BG. Mortality in Mental Disorders and Global Disease Burden Implications: A Systematic Review and Meta-analysis. JAMA Psychiatry. 2015;72(4):334. doi:10.1001/JAMAPSYCHIATRY.2014.2502
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/JAMAPSYCHIATRY.2014.2502&link_type=DOI) 

4.  4.Kivimäki M, Batty GD, Pentti J, et al. Association between socioeconomic status and the development of mental and physical health conditions in adulthood: a multi-cohort study. The Lancet Public Health. 2020;5(3):e140–e149. doi:10.1016/S2468-2667(19)30248-8
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2468-2667(19)30248-8&link_type=DOI) 

5.  5.Bjelland I, Krokstad S, Mykletun A, Dahl AA, Tell GS, Tambs K. Does a higher educational level protect against anxiety and depression? The HUNT study. Social Science and Medicine. 2008;66(6):1334–1345. doi:10.1016/j.socscimed.2007.12.019
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.socscimed.2007.12.019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18234406&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254266900008&link_type=ISI) 

6.  6.Erickson J, El-Gabalawy R, Palitsky D, et al. Educational attainment as a protective factor for psychiatric disorders: Findings from a nationally representative longitudinal study. Depression and Anxiety. 2016;33(11):1013–1022. doi:10.1002/da.22515
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/da.22515&link_type=DOI) 

7.  7.Cohen AK, Nussbaum J, Weintraub MLR, Nichols CR, Yen IH. Association of adult depression with educational attainment, aspirations, and expectations. Preventing Chronic Disease. 2020;17:200098. doi:10.5888/pcd17.200098
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.5888/pcd17.200098&link_type=DOI) 

8.  8.Niemeyer H, Bieda A, Michalak J, Schneider S, Margraf J. Education and mental health: Do psychosocial resources matter? SSM - Population Health. 2019;7:100393. doi:10.1016/J.SSMPH.2019.100392
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.SSMPH.2019.100392&link_type=DOI) 

9.  9.Jami ES, Hammerschlag AR, Bartels M, Middeldorp CM. Parental characteristics and offspring mental health and related outcomes: a systematic review of genetically informative literature. Translational Psychiatry. 2021;11(1):1–38. doi:10.1038/S41398-021-01300-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41398-021-01694-z&link_type=DOI) 

10. 10.McLaughlin KA, Gadermann AM, Hwang I, et al. Parent psychopathology and offspring mental disorders: results from the WHO World Mental Health Surveys. The British Journal of Psychiatry. 2012;200(4):290–299. doi:10.1192/BJP.BP.111.101253
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImJqcHJjcHN5Y2giO3M6NToicmVzaWQiO3M6OToiMjAwLzQvMjkwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMTEvMjAyNC4wOC4xMC4yNDMxMTc4OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

11. 11.Davis-Kean PE, Tighe LA, Waters NE. The Role of Parent Educational Attainment in Parenting and Children’s Development. Current Directions in Psychological Science. 2021;30(2):186–192. doi:10.1177/0963721421993116/ASSET/IMAGES/LARGE/10.1177_0963721421993116-FIG1.JPEG
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0963721421993116&link_type=DOI) 

12. 12.Ayano G, Lin A, Dachew BA, Tait R, Betts K, Alati R. The impact of parental mental health problems on the educational outcomes of their offspring: Findings from the Raine Study. The Australian and New Zealand journal of psychiatry. 2022;56(5):510–524. doi:10.1177/00048674211025633
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/00048674211025633&link_type=DOI) 

13. 13.Evans DM, Davey Smith G. Mendelian randomization: New applications in the coming age of hypothesis-free causality. Annual Review of Genomics and Human Genetics. 2015;16:327–350. doi:10.1146/ANNUREV-GENOM-090314-050016
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-genom-090314-050016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25939054&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 

14. 14.Jones DP, Wootton RE, Gill D, et al. Mental health as a mediator of the association between educational inequality and cardiovascular disease: A mendelian randomization study. Journal of the American Heart Association. 2021;10(17):19340. doi:10.1161/JAHA.120.019340
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/JAHA.120.019340&link_type=DOI) 

15. 15.Yuan S, Xiong Y, Michaëlsson M, Michaëlsson K, Larsson SC. Genetically predicted education attainment in relation to somatic and mental health. Scientific Reports. 2021;11(1):4296. doi:10.1038/s41598-021-83801-0
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-021-83801-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33619316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 

16. 16.Brumpton B, Sanderson E, Heilbron K, et al. Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses. Nature Communications 2020 11:*1*. 2020;11(1):1-13. doi:10.1038/s41467-020-17117-4
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-17117-4&link_type=DOI) 

17. 17.Kong A, Thorleifsson G, Frigge ML, et al. The nature of nurture: Effects of parental genotypes. Science. 2018;359(6374):424-428. doi:10.1126/science.aan6877
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTkvNjM3NC80MjQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8xMS8yMDI0LjA4LjEwLjI0MzExNzg5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

18. 18.Morris TT, Davies NM, Hemani G, Smith GD. Population phenomena inflate genetic associations of complex social traits. Science Advances. 2020;6(16). doi:10.1126/SCIADV.AAY0328/SUPPL\_FILE/AAY0328\_SM.PDF
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access\_num=10.1126/SCIADV.AAY0328/SUPPL_FILE/AAY0328_SM.PDF&link_type=DOI) 

19. 19.Hartwig FP, Davies NM, Davey Smith G. Bias in Mendelian randomization due to assortative mating. Genetic Epidemiology. 2018;42(7):608–620. doi:10.1002/gepi.22138
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.22138&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29971821&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 

20. 20.Hartwig FP, Davies NM, Davey Smith G. Bias in Mendelian randomization due to assortative mating. Genetic Epidemiology. 2018;42(7):620. doi:10.1002/GEPI.22138
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/GEPI.22138&link_type=DOI) 

21. 21.Olav B, Svold Å, Langhammer A, et al. Cohort Profile Update: The HUNT Study, Norway. International Journal of Epidemiology. 2022;2022:1–12. doi:10.1093/IJE/DYAC095
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/IJE/DYAC095&link_type=DOI) 

22. 22.Bycroft C, Freeman C, Petkova D, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562(7726):203-209. doi:10.1038/s41586-018-0579-z
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 

23. 23.Okbay A, Wu Y, Wang N, et al. Polygenic prediction of educational attainment within and between families from genome-wide association analyses in 3 million individuals. Nature Genetics. 2022;54(4):437–449. doi:10.1038/S41588-022-01016-Z
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588--01016-z&link_type=DOI) 

24. 24.Choi SW, Mak TSH, O’Reilly PF. A guide to performing Polygenic Risk Score analyses. Nat Protoc. 2020;15(9):2759–2772. doi:10.1038/s41596-020-0353-1
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41596-020-0353-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32709988&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 

25. 25.Spitzer RL, Kroenke K, Williams JBW, Löwe B. A Brief Measure for Assessing Generalized Anxiety Disorder: The GAD-7. Archives of Internal Medicine. 2006;166(10):1092–1097. doi:10.1001/archinte.166.10.1092
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinte.166.10.1092&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16717171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237720500008&link_type=ISI) 

26. 26.Kroenke K, Spitzer RL, Williams JBW. The PHQ-9. J Gen Intern Med. 2001;16(9):606–613. doi:10.1046/j.1525-1497.2001.016009606.x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1525-1497.2001.016009606.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11556941&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171184700005&link_type=ISI) 

27. 27.Zigmond AS, Snaith RP. The hospital anxiety and depression scale. Acta psychiatrica Scandinavica. 1983;67(6):361–370. doi:10.1111/J.1600-0447.1983.TB09716.X
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1600-0447.1983.tb09716.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6880820&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1983QV25600001&link_type=ISI) 

28. 28.Eysenck SBG, Tambs K. Cross-cultural comparison of personality: Norway and England. Scand J Psychol. 1990;31(3):191–197. doi:10.1111/j.1467-9450.1990.tb00830.x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1467-9450.1990.tb00830.x&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990EC47900003&link_type=ISI) 

29. 29.R Core Team. R: A Language and Environment for Statistical Computing. Published online 2021. [https://www.r-project.org/](https://www.r-project.org/)
    
    

30. 30.Bergé L. Efficient estimation of maximum likelihood models with multiple fixed-effects: the R package FENmlm. DEM Discussion Paper Series. Published online 2018. Accessed May 30, 2023. [https://ideas.repec.org//p/luc/wpaper/18-13.html](https://ideas.repec.org//p/luc/wpaper/18-13.html)
    
    

31. 31.Viechtbauer W. Conducting Meta-Analyses in R with the metafor Package. Journal of Statistical Software. 2010;36:1–48. doi:10.18637/jss.v036.i03
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v036.i03&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21401589&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 

32. 32.Hemani G, Zheng J, Elsworth B, et al. The MR-base platform supports systematic causal inference across the human phenome. eLife. 2018;7:e34408. doi:10.7554/ELIFE.34408
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/ELIFE.34408&link_type=DOI) 

33. 33.Howe LJ, Nivard MG, Morris TT, et al. Within-sibship GWAS improve estimates of direct genetic effects. *bioRxiv*. Published online March 7, 2021:2021.03.05.433935. doi:10.1101/2021.03.05.433935
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMS4wMy4wNS40MzM5MzV2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzExLzIwMjQuMDguMTAuMjQzMTE3ODkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

34. 34.van Buuren S, Groothuis-Oudshoorn K. mice: Multivariate Imputation by Chained Equations in R. Journal of Statistical Software. 2011;45(3):1–67. doi:10.18637/JSS.V045.I03
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/JSS.V045.I03&link_type=DOI) 

35. 35.Jiang W, Lu Y, Xie H. Education and mental health: Evidence and mechanisms. Journal of Economic Behavior & Organization. 2020;180:407–437. doi:10.1016/j.jebo.2020.09.032
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jebo.2020.09.032&link_type=DOI) 

36. 36.Mcmartin A, Conley D. Commentary: Mendelian randomization and education–Challenges remain. International Journal of Epidemiology. 2020;49(4):1193–1206. doi:10.1093/IJE/DYAA160
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/IJE/DYAA160&link_type=DOI) 

37. 37.Border R, Athanasiadis G, Buil A, et al. Cross-trait assortative mating is widespread and inflates genetic correlation estimates. Science. 2022;378(6621):754–761. doi:10.1126/SCIENCE.ABO2059
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abo2059&link_type=DOI) 

38. 38.Polushina T, Banerjee N, Giddaluru S, et al. Identification of pleiotropy at the gene level between psychiatric disorders and related traits. Translational Psychiatry. 2021;11(1):410. doi:10.1038/S41398-021-01530-4
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41398-021-01530-4&link_type=DOI) 

39. 39.Hilker R, Helenius D, Fagerlund B, et al. Heritability of schizophrenia and schizophrenia spectrum based on the nationwide danish twin register. Biological Psychiatry. 2018;83(6):492–498. doi:10.1016/J.BIOPSYCH.2017.08.017
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2017.08.017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F11%2F2024.08.10.24311789.atom) 

40. 40.de Zeeuw EL, Hottenga JJ, Ouwens KG, et al. Intergenerational transmission of education and ADHD: Effects of parental genotypes. Behavior Genetics. 2020;50(4):232. doi:10.1007/S10519-020-09992-W
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/S10519-020-09992-W&link_type=DOI) 

41. 41.Tyrrell J, Zheng J, Beaumont R, et al. Genetic predictors of participation in optional components of UK Biobank. Nature Communications. 2021;12(1):886. doi:10.1038/S41467-021-21073-Y
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41467-021-21073-Y&link_type=DOI) 

42. 42.Gkatzionis A, Burgess S. Contextualizing selection bias in Mendelian randomization: how bad is it likely to be? International Journal of Epidemiology. 2019;48(3):701. doi:10.1093/IJE/DYY202
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/IJE/DYY202&link_type=DOI)