Abstract
Disentangling direct and indirect genetic pathways underlying the intergenerational transmission of emotional problems could guide preventative strategies and further the understanding of the role of parental mental health in children’s outcomes. This study aimed to estimate the extent to which genetic pathways that are direct (via child genotype) and indirect (e.g., via parental phenotype) explain the well-established association between parent and child emotional problems. We leveraged data from seven European cohort studies with a combined population of Ntrios=15,475. Polygenic scores were calculated for parental and offspring neuroticism, as it represents a dispositional trait underlying emotional problems. Emotional problems in offspring were measured using validated scales across various developmental stages from early childhood to adulthood. We used neuroticism polygenic scores within a structural equation modelling framework to distinguish between direct genetic pathways from parental genotype to offspring outcome (acting through offspring genotype), and indirect genetic pathways (acting through parental phenotype and associated environment). Standard errors for direct genetic, indirect genetic and total effects were bootstrapped and meta-analyses pooled effect estimates at three developmental stages (childhood: 3-4 years, adolescence: 11-13 years, adulthood: 18+ years). We found evidence suggesting an indirect genetic pathway between mothers and child emotional problems during early childhood (pooled estimate, mean difference in standardised child emotional problems score per 1SD increase in maternal PGS for neuroticism=0.04, 95% CI: 0.01, 0.07). This association attenuated over child development, while direct genetic pathways strengthened. High attrition rates, measurement error and low variance explained by polygenic scores may have altered precision of the estimates, influencing the interpretation of the results. However, we provide the first multi-cohort study to provide evidence for an indirect genetic pathway from maternal neuroticism to early child emotional problems. This suggests that there are likely processes other than direct genetic pathways involved in the intergenerational transmission of emotional problems, highlighting the importance of timely support to prevent and reduce emotional issues in mothers as a preventative strategy for emotional difficulties.
Introduction
The prevalence of emotional problems, such as depression and anxiety(1), is increasing, particularly among young adults(2,3). Despite decades of research into ‘causal’ risk factors, policy changes, and intervention development, the burden of mental health disease has risen(4), with mental illness now being a leading cause of disability worldwide(5). Developing effective preventive interventions for emotional difficulties could have major societal and economic benefits by improving population health.
There is a well-established association between parent and child emotional problems. Family history of mental health problems for instance is considered one of the most important risk factors for the onset and maintenance of psychiatric disorders(6,7). In the UK, children with parents experiencing high levels of psychological distress are five times more likely to suffer from emotional disorders(8). This intergenerational risk is transmitted through both genetic and environmental mechanisms, as well as a likely interplay among the two. A key question is to identify potentially modifiable pathways to this intergenerational transmission of poor mental health, as different mechanisms of transmission will inform different preventative strategies.
Beyond an individuals’ genetic predisposition to a specific disorder (i.e., transmitted by the parents - direct genetic effect), non-transmitted parental genes can act as risk factors for the onset of emotional problems (also known as indirect genetic effect, “genetic nurture”, or familial effects(9)) by influencing aspects of the environment that are likely to contribute to child psychopathology and/or moderate the association between the child’s own genetic liability and mental health outcomes. Recent observational evidence and meta-analyses of genetically informed studies, including twin and adoption studies, have highlighted a clear role, albeit relatively small (12-20% of the variance), for the familiar environment in the development of psychopathology(10,11). A comprehensive review of genetically informative study designs, examining the association between parental and offspring emotional problems, reported strong evidence of both genetic and environmental transmission of parental depression to emotional problems across childhood, adolescence, and adulthood(12). Genetically informed studies can thus help to ascertain whether the associations between parental risk factors and offspring emotional outcomes are influenced by direct environmental effects or shared genetic factors. These designs offer improvements in accounting for genetic and shared environmental confounding factors compared to traditional observational studies(13,14).
Neuroticism is a personality trait indicating an individual’s tendency in experience negative emotions more frequently and intensely after exposure to stressors(15,16). It has also been conceptualised as a dispositional risk trait for emotional problems(16), including depression and anxiety, and its utility in research and clinical applications has been discussed elsewhere(17–19). Recent research shows that a genetic liability to neuroticism (proxied by a polygenic score – PGS) associates with emotional problems and depressive symptoms from infancy to young adulthood(20,21). Because neuroticism associates with both depression and anxiety, we decided to examine genetic liability to neuroticism as a common susceptibility to emotional problems. Using a genetic liability to a continuous personality trait as compared to a clinical disorder may be most useful when investigating symptom levels within general populations(22,23).
Attempting to isolate the role of parental phenotype is important because if parents’ emotional difficulties influence offspring via non-genetic paths, timely treatment of parental mental health symptoms alone could reduce the risk in the child. However, if there is only evidence for direct genetic pathways, this could suggest that focusing solely on parental mental health might not be as effective in prevention of child mental health risk as those including the child or with a specific focus on child’s mental health.
The aim of the current study is to estimate the extent to which the association between parent and child emotional problems is influenced by indirect genetic pathways and direct genetic pathways(9). Indirect genetic pathways, such as those mediated through parental phenotype, offer an opportunity to distinguish pathways operating through parental characteristics, distinct from pathways involving direct transmission of genetic information to the child (Figure 1).
To address this aim, we used polygenic scores (PGS) for neuroticism, obtained from mother-father-offspring trios, and a pathway approach disentangling direct genetic and indirect genetic pathways to emotional disorders in offspring in childhood, adolescence, and adulthood (Figure 2). While previous studies using similar methods have identified indirect genetic pathways in educational outcomes(24), this study aims to extend our current understanding to mental health outcomes(25) by aggregating data from seven cohorts to enhance statistical power and to balance different sources of bias (i.e., different cohorts have different attrition patterns and reporting of emotional problems)(26).
Understanding the potential contribution of direct and indirect pathways may help guiding future strategies to prevent intergenerational cycles of emotional difficulties and improve child mental health in those who have a parent with emotional difficulties.
Samples and methods
Participants and genotyped data
This study used data from seven cohort studies across five countries, constituting the recently established Genetic and Environmental Mental health INtergenerational Investigation (GEMINI) consortium. We included all cohorts with both genetic and phenotypic data on the family trio. Detailed information for each cohort is provided in the Supplementary Information. In brief, this study combined information from the Avon Longitudinal Study of Parents and Children (ALSPAC)(27–29), Millenium Cohort Study (MCS)(30), Generation Scotland (GS)(31), the Norwegian Mother, Father and Child Cohort Study (MoBa)(32), FinnBrain (FB)(33), CHILD-SLEEP (CS)(34), and the Pre-, Peri-, and Postnatal Stress: Epigenetic Impact on Depression study (POSEIDON)(35,36). Please note that the ALSPAC website contains details of all the data that is available through a fully searchable data dictionary and variable search tool (http://www.bristol.ac.uk/alspac/researchers/our-data/).
Measures
Offspring emotional problems and depressive symptoms were measured using different scales across development (validated for children’s ages) and reported by different assessors (i.e., parents and children – self-reported). Emotional problem scores were based on the following measures (STable 1) and standardised prior use in analyses to facilitate comparison and interpretability. Childhood emotional problems were measured using the Brief Infant-Toddler Social and Emotional Assessment (BITSEA)(37), the Child Behaviours Checklist (CBCL)(38), and Strengths and Difficulties Questionnaire (SDQ)(39). Adolescent emotional problems were measured using the Strengths and Difficulties Questionnaire (SDQ)(39) and the Short Moods and Feelings Questionnaire (SMFQ)(40). Finally, adult depressive symptoms were measured using The Clinical Interview Schedule-Revised (CIS- R)(41), General Health Questionnaire (GHQ)(42), and the Patient Health Questionnaire-2 (PHQ2)(43)
Polygenic scores indexing parental and offspring neuroticism
Information on genotyping in each cohort is available in Supplementary Information. SNPs were filtered on minor allele frequency (>1%) and imputation info score (>0.8). In each cohort, PGS for neuroticism were derived based on publicly available summary statistics from a genome-wide association study by Luciano et al. (44). This study analysed over 329,000 individuals of self-described white British ancestry in the UK Biobank Study (46.3% male). Neuroticism was measured using a 12-item scale from the Eysenck Personality Questionnaire-Revised (EPQ-R). Neuroticism scores were adjusted for age, sex, genotyping batch, genotyping array, assessment centre and 40 principal components of ancestry prior to analysis. PGS were derived in each of the cohorts using PRSice2 with default clumping options (clumping correlated SNPs within a 250kb window at a r2 threshold of 0.1) and based on a p-value threshold of p<0.05. This threshold was chosen to maximise the variance explained by the PGS while minimising noise. All scores were standardised as z- scores prior to use in analyses. Betas (β) indicate a unit difference increase or reduction in the standardised outcome score for each 1SD increase in the exposure (i.e., PGS for neuroticism of child, mother, or father).
Data analysis
Descriptive analyses and missing data
We calculated medians and interquartile ranges (IQR) or frequencies and percentages for basic demographic characteristics of the participants across cohort studies and reported these data for mothers, fathers, children and trios with genetic data available (that were included in analyses because had the outcome), and those with missing genetic and outcome data (STable 2).
Primary analyses
In each cohort we used neuroticism PGS from the mothers, fathers, and offspring within a SEM (Structural Equation Modelling) path analysis using Stata v16(45) (ALSPAC, MCS) or lavaan package in R (GS, MoBa, FinnBrain, ChildSleep, POSEIDON) (see supplementary file for exact coding) to separate out the direct genetic pathways of parental genotype (which acts via offspring genotype) on offspring emotional problems from the indirect genetic pathway that could act via maternal phenotype and environment (Figure 2). More specifically maternal and paternal PGS scores were regressed on child PGS scores and in turn child emotional outcomes and on child emotional outcomes directly. Because we included all pathways (both direct and indirect) from maternal and paternal PGS in a simultaneous model the path estimates allow us to deconstruct direct and indirect contributions. Among cohorts where population structure is more varied and thus more likely to introduce issues due to confounding by population stratification (e.g. GS, MCS), we accounted for 10 principal components of ancestry.
Standard errors were bootstrapped in each cohort and the effect estimates and 95% confidence intervals for the indirect and direct effects of maternal and paternal PRS were pooled in a random effects meta-analysis at three separate time points (childhood: age 2-4 years, adolescence: 11-13 years, adulthood: 18-99 years).
The measures were also combined according to rater (self-report or parent-report) across all time points in a random effects meta-analysis. This study follows STROBE reporting guidelines (STable 3).
Results
Samples descriptives
Sample demographic characteristics are summarised in Table 2. In ALSPAC, mothers with genetic data were more likely to be married or co-habiting, and were less likely to have ever smoked compared to mothers for whom genetic data was missing. In ALSPAC, MoBa, and GS, individuals with completed trio genetic data were more likely to have mothers with higher educational qualifications as compared to the full sample, while this is not the case for MCS(46,47). These descriptive results are summarised in STable 2 and reported in more details here(48).
Our pooled sample was 15,475 in the childhood analysis (aged 2-4 years), 4,507 in the adolescent analysis (11-13 years), and 4,244 in the adult analysis (18-99 years).
Direct and indirect genetic effects on childhood emotional problems in six cohorts
In childhood, data on parent-reported emotional problems in childhood was available from 6 of the cohorts (i.e., ALSPAC, Child Sleep, FinnBrain, MCS, MoBa, Poseidon). Evidence was strongest for an association between the maternal indirect genetic association (0.04SD difference in emotional problems per 1 standard deviation (SD) in maternal neuroticism PGS (95% CI: 0.01, 0.07), p=0.018) and child emotional problems. The effect estimates for direct genetic associations and paternal indirect genetic associations were in the same direction, but evidence was weaker for these effects (Figure 3, STable 4). However, the CIs for each of the estimates overlapped. We found some evidence of both maternal and paternal total effects.
Direct and indirect genetic effects on adolescent emotional problems and depressive symptoms in two cohorts
In adolescence, data on offspring emotional problems were available from ALSPAC and MCS cohorts (N=4,507). Here, we found some evidence of an association between both maternal and paternal direct genetic pathways and offspring emotional problems (Figure 4, STable 4). The direction of this association was consistent with that observed during childhood. For maternal direct genetic effects there was a 0.04SD difference in emotional problems per 1SD increase in maternal neuroticism PGS ((95% CI: 0.01, 0.06), p=0.010), while for paternal direct genetic effects there was a 0.03SD in difference in offspring emotional problems per 1SD increase in paternal neuroticism PGS (95% CI: 0.01, 0.06; p=0.006). While the direct genetic effect estimates are consistent between ALSPAC and MCS, the indirect genetic effect estimates appear to be in conflicting directions. In a similar way to early childhood, the CIs for all pathways overlap and the evidence for a total effect is less clear.
Direct and indirect genetic effects on adulthood depressive symptoms in three cohorts
In adulthood, self-reported depressive symptoms data were available from ALSPAC, GS, and MCS cohorts (N=4,210). We found strong evidence for an association between both maternal and paternal direct genetic effects and offspring emotional problems (Figure 5, STable 4). For maternal direct genetic effects there was a 0.05SD difference in emotional problems per 1SD increase in maternal neuroticism PGS (95% CI: 0.02, 0.08; p<0.001), while for paternal direct genetic effects there was a 0.05SD in difference in offspring depressive symptoms per 1SD increase in paternal neuroticism PGS (95% CI: 0.02, 0.07; p<0.001). We found no evidence of an association between indirect genetic effects and offspring depressive symptoms. Confidence intervals, in particular for the maternal effects, were overlapping across pathways, and evidence of total effects was unclear, particularly for paternal effects on offspring depressive symptoms in adulthood.
There were modest positive associations reflective of maternal and paternal potentially modifiable phenotypes and as well as both parent genetic inheritance. This was strongest for maternal phenotype (0.04SD difference in emotional problems per 1SD in maternal neuroticism PGS (95%CI: 0.01, 0.07)) and weakest in magnitude for maternal and paternal genetic inheritance (both 0.01 SD per 1SD PGS (-0.01, 0.03)). All pooled estimates were imprecise.
Discussion
This study leveraged data from seven cohort studies across five countries to elucidate the association between parental genetic predisposition to neuroticism and child emotional problems and depressive symptoms through both direct and indirect genetic pathways. Using a genetically informed approach and integrated polygenic scores within a structural equation model, we disentangled the influence of direct and indirect genetic pathways on this association. We also investigated whether this association persisted from early childhood to adulthood by performing a meta-analysis of emotional problems in early childhood (age 2-4 years), adolescence (11-13 years), and adulthood (18 years and over).
Overall, our pooled results report suggestive evidence of an association between indirect genetic maternal pathways and child emotional problems during early childhood (n=15,475). The indirect maternal genetic effects diminished over time while direct genetic effects intensified across the life course (in more than 4,000 individuals). The strength of the evidence supporting an indirect genetic pathway drastically reduced in adolescence and adulthood, with the magnitude of the effect size halving from childhood to adolescence and from adolescence to adulthood for the maternal indirect effect, and with inconsistent direction of the effect estimate for the paternal indirect effect. Studies with genetics data, including twins and children-of-twin studies, have suggested that familiar influences and indirect genetic effects may play a greater role on offspring phenotypes in early ages as compared to adolescence and adulthood where the effect of the own genetic liability is more likely to contribute to the phenotype(49–52). These hypotheses are also supported by evidence genetic heritability for several traits, including mental health related traits, increases over the lifespan(11,53).
Alternatively, it is possible that we observed stronger evidence for an indirect pathway in early childhood because of the use of parent vs self-reported measures at these time points, which may have induced differential biases in the reporting child emotional symptoms associated to maternal mental health characteristics. Nevertheless, some of these biases have been explored in other studies, even in ALSPAC, reporting often negligible effects(20,54) (see in limitations).
Unlike the direct genetic effects which show a consistent direction of association, the effect estimates for the indirect genetic pathway also appear to be in conflicting directions across ALSPAC and MCS in adolescence. The variability of the indirect pathway could suggest that the environment captured by the indirect genetic component is likely to differ across studies (e.g., ALSPAC participants were recruited a decade earlier and from a select region of the UK, while MCS was sampled from across the entire UK), thus being explained by differences specific to each cohort (e.g., measures, drop-out rates, selection processes).
The context within the literature
The estimates for direct genetic pathways become stronger in magnitude and strength of evidence during adolescence and adulthood. This may be potentially due to a diminishing role of the parental environment, as captured by the genetic component, in shaping the offspring’s environment during these developmental stages. An alternative explanation could be due to the use of polygenic scores derived from adult GWAS summary statistics in our analysis. These scores may therefore explain more of the variation in the phenotype as cohort members reach adulthood(55). However, previous work has shown that even at younger ages the neuroticism polygenic scores are still associated with the phenotypes(20).
Implications for intergenerational psychiatry
Our findings suggest that direct genetic pathways more strongly explain the associations between parental neuroticism and their children’s emotional difficulties in adolescence and adulthood. Contrarily, indirect genetic pathways, especially via maternal ‘genetic nurturing’, more strongly explain the association with child emotional problems in early childhood. These findings have separate potential implications for preventative interventions (Figure 6). In the case of direct genetic pathways, delivering preventative interventions exclusively targeting parental mental health is unlikely to have beneficial effects on child emotional problems in adolescence and adulthood. This, however, contrasts with trial evidence(56) that targeting parental mental health or parenting skills in parents with mental health disorders have modest effects in reducing internalising disorders in children regardless of their age (effects were observed in both childhood and adolescence). However, the magnitude of the effects is relatively small and interventions directly supporting child resilience, emotional regulation or parenting interventions with active involvement of the child may be more effective in reducing emotional problems. Such findings align though with meta-analyses of trials that found that parenting interventions with an active inclusion of the child and focus on the relationship are more effective in reducing internalising problems as interventions that involved parents only(57,58).
In the case of genetic nurture, the effects of interventions targeting parental mental health (especially maternal) would be expected to have beneficial downstream effects on child emotional problems as experienced in early childhood.
Limitations
This study may be subject to several limitations which should be considered when evaluating its findings. First, despite being the largest study of its kind, limitations persist due to sample size constraints. The need for genetic data from mothers, fathers, and offspring has restricted the available sample, as many studies only offer genetic data for parent- offspring pairs. However, combining data across multiple cohorts has maximised the sample size, especially in early childhood, where statistical power is greater as compared to the other timepoints. Nevertheless, missing genetic information data and severe attrition over time may have introduced bias into estimates, as those providing genetic information and remaining in the cohort were more likely to come from more socio-economically advantaged families, especially in ALSPAC, CHILD-SLEEP, MoBa, and GS (see STable 2). This pattern of missingness could have influenced our findings, most likely biasing them toward the null as it known that higher education moderates the effects of poor maternal mental health on offspring emotional and behavioural difficulties (47). Future studies should aim to employ methods such as parental genetic imputation (59), inverse probability weighting(61) or multiple imputation(62) to minimise these issues. In addition, we have used data from five European countries therefore limiting the generalisability to this study to non-European participants.
Second, despite the increase in sample size of GWAS, they are still limited in estimating all the phenotypic variance explained by common variants, often resulting in these scores explaining only a small portion of the phenotypic variance. Consequently, when the estimates are used as weights in the construction of the polygenic score, they introduce noise(9). This noise can lead to observed null associations between the resulting polygenic scores and outcomes of interest, potentially reflecting false negatives where a true association exists but remains undetected.
Third, depression or emotional problems were primarily collected by a single reporter within each study and at a single time point which is likely to increase measurement error. Different reporters of internalising problems lead, for example, to different heritability estimates(63). To address potential reporting biases, having multiple raters assessing the phenotype would have been beneficial. However, this is challenging in childhood due to limited validation of child or teacher reports. Moreover, subjective ratings of mental health difficulties present inherent limitations regardless of the reporter. According to the ‘distortion model’(64–66), the mental health characteristics of the reporters may systematically bias the ratings of mental health questionnaires, with implications for assessing causality in genetic studies (9). However, previous investigations into such biases have found limited support for dependent differential misclassification or measurement of outcomes (20,54).
Finally, we cannot rule out possible bias due to assortative mating or population stratification, even though we have likely accounted for the latter by including a large number of principal components of genetic ancestry.
Conclusions and direction for future research
In conclusion, our findings suggest that indirect genetic processes, such as the environment provided by the parent, are most strongly associated with offspring mental health in early life. Conversely, pathway estimates for directly transmitted genetic information become larger during adolescence and adulthood, although the individual reporting on emotional problems may play a role in this. This insight holds potential for refining intervention strategies aimed at mitigating the impact on offspring mental health and importantly the potential role of including both parent and child from as early as possible. The magnitude of the association was larger for the maternal indirect pathways as compared to the paternal indirect genetic pathways which may indicate the degree of parental involvement. Indeed, variation in parental involvement in children’s upbringing, such as time spent at home and involvement after separation, may also moderate the effects of indirect genetic effects on offspring phenotype with potentially more pronounced effects for the more involved parents. However, actual parental involvement was not explored in our study, and we recommend this for future investigations.
Data Availability
Data from each individual cohort can be obtained through individual data request and will vary depending on the agreement with each cohort study. For example, ALSPAC data access is through a system of managed open access. Information on how to request access to ALSPAC data can be found on the ALSPAC website (http://www.bristol.ac.uk/alspac/researchers/access/). Access to data sets requires approval from The Regional Committee for Medical and Health Research Ethics in Norway and an agreement with MoBa. For the MCS, data access is free of charge and it is detailed here https://cls.ucl.ac.uk/data-access-training/data-access/.
Fundings
The UK Medical Research Council and Wellcome (Grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. This work was supported by European Research Council (ERC) under the European Union’s Seventh Framework Programme (grant FP/2007-2013), European Research Council Grant Agreements (grants 758813; MHINT). HS, RW, DAL, NJT work in or are affiliated with a unit that receives funding from the University of Bristol and UK Medical Research Council (MC_UU_00032/05, MC_UU_00032/06 and MC_UU_00032/07). DAL is contribution to this study is supported by the European Research Council (Grant Agreement 101021566 ART-Health), GWAS data was generated by Sample Logistics and Genotyping Facilities at Wellcome Sanger Institute and LabCorp (Laboratory Corporation of America) using support from 23andMe. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgments.pdf). A. S. F. Kwong is supported by a Wellcome Early Career Award (Grant ref: 227063/Z/23/Z).
Generation Scotland: Scottish Family Health Study was supported by the Chief Scientist Office (CZD/16/6) and the Scottish Funding Council (HR03006) and by the Wellcome Trust (104036/Z/14/Z).
The POSEIDON study was supported by a grant of the Dietmar Hopp Foundation.
The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services, and the Ministry of Education and Research. A.H. is supported by the Research Council of Norway (#274611, #336085 and #288083), the South- Eastern Norway Regional Health Authority (#2020022), and her contribution to this study is supported by the European Union’s Horizon Europe Research and Innovation Programme (FAMILY, grant agreement No 101057529). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union, or the other funding agents. Neither the European Union nor the granting authorities can be held responsible for them.
Author contributions
This publication is the work of the authors and RMP, IC, and HS will serve as guarantors for the contents of this paper.
Acknowledgments
We are extremely grateful to all the families who took part in this study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists, and nurses. The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. We are grateful to all the participating families in Norway who take part in this on-going cohort study. We thank the Norwegian Institute of Public Health (NIPH) for generating high-quality genomic data. This research is part of the HARVEST collaboration, supported by the Research Council of Norway (NRC) (#229624). We also thank the NORMENT Centre for providing genotype data, funded by NRC (#223273), South East Norway Health Authority and KG Jebsen Stiftelsen. Further we thank the Center for Diabetes Research, the University of Bergen for providing genotype data and performing quality control and imputation of the data funded by the ERC AdG project SELECTionPREDISPOSED, Stiftelsen Kristian Gerhard Jebsen, Trond Mohn Foundation, NRC, the Novo Nordisk Foundation, the University of Bergen, and the Western Norway Health Authorities (Helse Vest). We thank all parents and children for taking part in the POSEIDON study. We thank our student employees and interns for their support with data acquisition and data entry. We are grateful to the GS:SHFS participants in their families. Analysis of GS:SHFS data made use of the resources provided by the Edinburgh Compute and Data Facility (ECDF) (http://www.ecdf.ed.ac.uk/).