Genome-wide association meta-analysis of childhood and adolescent internalising symptoms

Eshim S Jami, Anke R Hammerschlag, Hill F Ip, Andrea G Allegrini, Beben Benyamin, Richard Border, Elizabeth W Diemer, Chang Jiang, Ville Karhunen, Yi Lu, Qing Lu, Travis T Mallard, Pashupati P Mishra, Ilja M Nolte, Teemu Palviainen, Roseann E Peterson, Hannah M Sallis, Andrey A Shabalin, Ashley E Tate, Elisabeth Thiering, Natàlia Vilor-Tejedor, Carol Wang, Ang Zhou, Daniel E Adkins, Silvia Alemany, Helga Ask, Qi Chen, Robin P Corley, Erik A Ehli, Luke M Evans, Alexandra Havdahl, Fiona A Hagenbeek, Christian Hakulinen, Anjali K Henders, Jouke Jan Hottenga, Tellervo Korhonen, Abdullah Mamun, Shelby Marrington, Alexander Neumann, Kaili Rimfeld, Fernando Rivadeneira, Judy L Silberg, Catharina E van Beijsterveldt, Eero Vuoksimaa, Alyce M Whipp, Tong Xiaoran, Ole A Andreassen, Dorret Boomsma, Sandra A Brown, S Alexandra Burt, William Copeland, Elizabeth J Costello, Danielle M Dick, Lindon J Eaves, K Paige Harden, Kathleen Mullan Harris, Catharina A Hartman, Joachim Heinrich, John K Hewitt, Christian Hopfer, Elina Hypponen, Marjo-Riitta Jarvelin, Jaakko Kaprio, Liisa Keltikangas-Järvinen, Kelly L Klump, Kenneth Krauter, Ralf Kuja-Halkola, Henrik Larsson, Terho Lehtimäki, Paul Lichtenstein, Sebastian Lundstrom, Hermine H Maes, Per Magnus, Marcus R Munafò, Jake M Najman, Pål R Njølstad, Albertine J Oldehinkel, Craig E Pennell, Robert Plomin, Ted Reichborn-Kjennerud, Chandra Reynolds, Richard J Rose, Andrew Smolen, Harold Snieder, Michael Stallings, Marie Standl, Jordi Sunyer, Henning Tiemeier, Sally Wadsworth, Tamara L Wall, Andrew J O Whitehouse, Gail M Williams, Eivind Ystrom, Michel G Nivard, Meike Bartels, Christel M Middeldorp
doi: https://doi.org/10.1101/2020.09.11.20175026

Abstract

Internalising symptoms in childhood and adolescence are as heritable as adult depression and anxiety, yet little is known of their molecular basis. This genome-wide association meta-analysis of internalising symptoms included repeated observations from 64,641 individuals, aged between 3 and 18. The N-weighted meta-analysis of overall internalising symptoms (INToverall) detected no genome-wide significant hits and showed low SNP heritability (1.66%, 95% confidence intervals 0.84-2.48%, Neffective=132,260). Stratified analyses indicated rater-based heterogeneity in genetic effects, with self-reported internalising symptoms showing the highest heritability (5.63%, 95% confidence intervals 3.08-8.18%). Additive genetic effects on internalising symptoms appeared stable over age, with overlapping estimates of SNP heritability from early-childhood to adolescence. Genetic correlations were observed with adult anxiety, depression, and the wellbeing spectrum (|rg|> 0.70), as well as with insomnia, loneliness, attention-deficit hyperactivity disorder, autism, and childhood aggression (range |rg|=0.42-0.60), whereas there were no robust associations with schizophrenia, bipolar disorder, obsessive-compulsive disorder, or anorexia nervosa. The pattern of genetic correlations suggests that childhood and adolescent internalising symptoms share substantial genetic vulnerabilities with adult internalising disorders and other childhood psychiatric traits, which could partially explain both the persistence of internalising symptoms over time and the high comorbidity amongst childhood psychiatric traits. Reducing phenotypic heterogeneity in childhood samples will be key in paving the way to future GWAS success.

Introduction

Internalising disorders, including anxiety and depression, are substantial contributors to the global burden of disease (1, 2). Whilst the estimated 12-month prevalence of depression and anxiety disorders in adults is 15% (3), internalising disorders are also present in early life, with an estimated prevalence of 2-3% of depression and 6-7% of anxiety in childhood and adolescence (4). Prior to the diagnosis of internalising disorders, as many as one in five children self-report internalising symptoms (5). These early symptoms of anxiety and depression appear to pose a long-term risk, as longitudinal studies show that internalising symptoms in childhood are associated with mood disorders, anxiety, and suicidality in adulthood (6-8). Findings from twin research show that internalising symptoms have a moderately strong genetic component. 40-50% of individual differences in internalising symptoms are explained by genetic factors (9-11). Moreover, research suggests that both stability and change in anxious and depressive symptoms from early childhood to adulthood are genetically influenced (10, 12-14). However, unlike adult anxiety and depression, investigation of the molecular genetic architecture of internalising symptoms in early life has received little attention thus far and to date, only two studies have applied a genome-wide approach (15, 16).

Published in 2013 and 2014, the first genome-wide association studies (GWASs) on childhood internalising symptoms did not identify any genome-wide significant hits for maternal-reported anxiety-related behaviours in children aged seven (N=2,810) (15), or internalising problems in children aged three (N=4,596) (16). Estimates of SNP-based heritability (the proportion of phenotypic variance captured by single nucleotide polymorphisms (SNPs) included in the GWAS), using genome-wide complex trait analysis (GCTA), were not robust in both studies (15, 16). Other GCTA studies similarly show mostly inconsistent and broad estimates of SNP heritability, mainly due to small sample sizes (17-22). Large-scale GWASs have led to significant discoveries in adult samples, with now 102 variants identified for depression (23) and 5 variants for anxiety (24). Given the comparable heritability estimates of adult and childhood internalising phenotypes, the next step in this line of research is to increase childhood sample sizes in order to generate sufficient power to capture the small effects of common variants that have been observed in adult studies.

Here, we present a genome-wide association meta-analysis which aims to identify common genetic variants associated with the development and course of internalising symptoms. The study combines repeated measurements of dimensional symptom scores from 22 independent cohorts of European ancestry, resulting in an overall sample of 64,641 individuals and 251,152 observations in children and adolescents aged between 3 and 18. All datasets were combined to produce a GWAS of overall internalising symptoms (INToverall), with an effective sample size of 132,260. Stratified analyses were used to investigate age, rater, and instrument-specific genetic effects. The overall GWAS of INToverall was followed up with gene-based analyses. Genetic overlap with external traits was examined by computing genetic correlations, with a focus on psychiatric phenotypes. Non-psychiatric traits were also investigated if they were previously found to be genetically correlated with adult anxiety and depression (23-25). Finally, polygenic scores were computed to test prediction of internalising symptoms in independent samples. With this study, we aim to gain insight into the genetic underpinnings of internalising symptoms throughout childhood and adolescence in order to improve our understanding of the development and progression of internalising disorders.

Methods

This project was pre-registered at the Open Science Framework (https://osf.io/edas6). Minor deviations from the pre-registration are explained in the Supplementary Note.

Sample and univariate analyses

The sample includes cohorts that are part of the EArly Genetics and Lifecourse Epidemiology (EAGLE) consortium behaviour and cognition working group (https://www.eagle-consortium.org/) (26) and additional cohorts with appropriate data. In total, 22 cohorts of European ancestry participated in the study. Ethical approval was provided by local committees at cohort level. Many cohorts were longitudinal birth or childhood cohort studies with long-term follow-up and multiple raters, e.g., mother, father, self and teacher. Repeated assessments of internalising symptoms within childhood and adolescence, from age 3 to age 18, were included. All cohorts performed univariate GWASs stratified by (i) age, (ii) rater, and (iii) instrument, with a minimum of 450 observations in each analysis. In the absence of diagnostic data, Internalising symptoms were dimensionally measured and positively scored on continuous scales, with higher scores indicating more internalising symptoms. Data was not dichotomized into a case-control design as this would have resulted in a reduction of statistical power (27). Detailed descriptions of the cohorts, phenotypic measures, and genotyping and imputation procedures can be found in Supplementary Tables 1-6 and the Supplementary Note.

Cohorts that included only unrelated subjects applied a linear regression model. Cohorts with a sample of related individuals corrected for non-independence of observations by either applying a mixed linear model or a sandwich correction of the standard errors. As there is no evidence for sex differences in genetic variation in psychiatric traits (28), analyses were not stratified by sex but gender was included as a covariate. Further details about the univariate GWASs are provided in the Supplementary Note. In total, 125 univariate GWASs were collated, with 251,152 observations based on 64,641 unique participants. The observations included ratings by mothers (40.7%), fathers (6.8%), teachers (18.3%), self (19.7%) and siblings (0.7%). An additional 13.8% of ratings were parental reports, where the informant was either the mother or the father. 15.1% of observations were in early-childhood (3 to 6 years), 36.0% in mid-childhood (7 to 10 years), 18.4% in late-childhood (11 to 12 years), and 30.0% in adolescence (13 to 18 years). Twelve instruments were used to measure internalising symptoms, of which the most commonly used were the Strengths and Difficulties Questionnaire (SDQ; 38.2%) (29), Achenbach System of Empirically Based Assessment (ASEBA; 36.7%) (30) and Rutter Children’s Questionnaires (8.2%) (31, 32).

Meta-analyses and the calculation of SNP heritabilities stratified by age, rater, and instrument

Quality control for each univariate GWAS was performed using EasyQC (Supplementary Text) (33). After QC, most cohorts retained between 3.4 and 7.1 million autosomal SNPs per GWAS (Supplementary Table 7). An exception was the Philadelphia Neurodevelopmental Cohort which retained fewer SNPs after merging data from different genotyping platforms. To account for dependency of repeated measurements of internalising symptoms within cohorts, the N-weighed meta-analysis approach was applied (34, 35). In short, two N × N matrices, representing sample overlap and phenotypic covariance within cohorts, were created, where N was the total number of univariate GWASs. As there was no overlap across cohorts, sample overlap and phenotypic covariance between cohorts were set to zero. Using the observed sample overlap within cohorts and their phenotypic covariance matrices, expected pairwise cross-trait intercept (CTI) values between GWASs were calculated. The pairwise CTI is approximately equal to the covariance between the test statistics from univariate GWASs. N-weighted meta-analyses were performed to obtain a multivariate test statistic per SNP, which represents a weighted sum of test statistics, adjusted by the CTI in order to account for sample overlap between GWASs. Formulas for the calculation of the multivariate test statistic for each SNP in the meta-analyses, the CTI between GWASs and estimation of effective sample size to account for repeated measurements (Neff) are provided in Ip et al. supplementary text (35).

A meta-analysis was performed based on the results of all available GWASs on internalising symptoms: INToverall SNP-based heritabilities (h2) were estimated using linkage disequilibrium score regression (LDSC) (36), first for INToverall, and next based on results of meta-analyses stratified according to rater, age, rater-by-age, and instrument (Supplementary Table 8). To ensure that the stratified analyses had sufficient power, a sample size threshold was set so that the total number of observations (Nobs) for each meta-analysis was at least 15,000. Rater-specific SNP heritabilities were estimated using assessments from parents (mother and/or father), mothers, fathers, teachers, and self, respectively. Age-specific SNP heritabilities focused on internalising symptoms during early childhood (3 to 6 years), mid-childhood (7 to 10 years), late-childhood (11 to 12 years), and adolescence (13 to 18 years). Rater-by-age SNP heritabilities assessed age effects within and between raters, provided that the univariate Nobs exceeded 15,000. Lastly, instrument-specific SNP heritabilities were calculated for SDQ, ASEBA, and Rutter for which the Nobs exceeded 15,000.

Genetic correlations across stratified GWAMAs were calculated using LDSC, but only if the z-score of the heritability estimate was ≥ 4, given that the heritability z-score is a good indicator of power and a score less than 4 is considered too noisy for meaningful estimates (37).

SNPs with minor allele frequency < 5% or Neff < 15,000 were removed from further analyses.

Gene-based analysis

Using summary statistics for INToverall, a MAGMA (38) gene-based test (v1.8, implemented in FUMA (39)) was performed to identify genes with a significant effect on internalising symptoms. The gene-based test applies a multiple regression model in which p-values from individual SNPs in a gene are combined into a test-statistic for each gene, while accounting for linkage disequilibrium between SNPs. European populations from the 1000 Genomes Phase III reference panel were used to estimate linkage disequilibrium. A total of 18,592 protein-coding genes were assessed for an association with internalising symptoms. A Bonferroni correction was applied to correct for multiple testing (α = 0.05 / 18,592; p < 2.69×10−06).

Tissue expression and gene-set analyses

Tissue enrichment and gene-set analyses were conducted in FUMA. The tissue enrichment analyses used two types of tissues from GTEx version 8: 30 general tissue types from multiple organs and 53 specific tissue types within these organs. A MAGMA gene-property test was performed to test one-sided relationships between cell type-specific gene expression and disease–gene associations. Bonferroni corrections were applied to correct for multiple testing for the general (α = 0.05 / 30; p < 1.7×10−04) and specific (α = 0.05 / 53; p < 9.4×10−04) tissue types.

The gene-set analysis was performed with default parameters in MAGMA v1.8. Gene-based P-values were converted to Z values and a between-gene correlation matrix was used as input to perform gene-set enrichment tests. Predefined gene sets from the molecular signature database MsigDB v7.0 were used. In total, 15,484 gene sets were tested. A Bonferroni correction was applied to correct for multiple testing (α = 0.05 / 15,484; p < 3.2×10−06).

Genetic correlations with external traits

Genetic correlations between internalising symptoms and other phenotypes were investigated using publicly available summary statistics for a curated set of traits (N=27). These primarily included adult psychiatric traits, in addition to other phenotypes selected based on previously identified correlations with adult anxiety and depression (23-25). Additionally, we obtained summary statistics from the GWA meta-analyses of overall and rater-specific childhood and adolescent aggression (35), that were based on overlapping cohorts and similar statistical methods, and calculated genetic correlations with these traits. The external traits and source studies are summarised in Supplementary Table 9. Summary statistics from INToverall and INTself (for which the z-score of the h2 was ≥ 4 (37)) were used. Genetic correlations were calculated using LDSC (36), which calculates genetic covariance between two traits based on all polygenic effects captured by included SNPs. Overlapping samples or population differences in GWAS summary statistics do not bias the computation of genetic correlations in LDSC. LDSC corrects for sample overlap by including a covariance matrix of the cross-trait LD score intercept, which is an estimate of sample overlap and phenotypic correlation. The genetic correlation estimate was based on the estimated slope from regressing the product of z-scores from two GWASs on the LD score. The LD scores used were computed using 1000 Genomes Phase III European data (37). Genetic correlations were considered significant at p < 9.26×10−04, after applying a Bonferroni correction for 54 independent tests.

Sensitivity analysis: polygenic score prediction

Polygenic score prediction of INToverall was tested as a sensitivity analysis. The Netherlands Twin Register (NTR) was used as the target sample to examine prediction of internalising symptoms in childhood and adolescence. We considered maternal-reported internalising symptoms at age 7 (N=3,845), and self-reported internalising symptoms during adolescence (age 13 to 18, N=2,679), using the ASEBA Child Behaviour Checklist and the Youth Self Report scales (30), respectively. A leave-one-cohort-out meta-analysis omitting data from NTR was performed for INToverall. The NTR target dataset was restricted to SNPs with minor allele frequency > 5% and imputation quality of R2 > 90%. Polygenic scores were constructed using LDpred (40), using a prior value of 0.5 to account for high polygenicity. Associations between polygenic scores of internalising symptoms and internalising problems were examined using Generalized Estimating Equations as implemented in the “gee” package in R (version 3.5.2). To account for relatedness in the target sample, the exchangeable working correlation matrix in gee was used, which applies a sandwich correction over the standard errors to account for clustering in the data. Age, sex, genotyping array, and the first 10 genetic principal components were included as covariates. Polygenic prediction was considered significant at p < 0.025, after applying a Bonferroni correction for 2 independent tests.

Results

Overall meta-analysis of childhood and adolescent internalising symptoms

The genome-wide association meta-analysis of INToverall found no genome-wide significant hits (Figure 1). Assuming a Neff of 132,260, SNP-based heritability of INToverall was estimated at 1.66% (95% confidence interval (CI) 0.84-2.48%). The mean chi-squared statistic was 1.086, with an LDSC-intercept of 1.043 (standard error (SE)= 0.0075), indicating that a small part of the inflation in test statistics might have been due to confounding biases, such as population stratification.

Figure 1
Figure 1

Manhattan plot of overall meta-analysis for childhood and adolescent internalising symptoms (INToverall). The solid line represents the significance threshold (p < 5×10−08) and the dotted line represents the suggestive threshold (p < 1×10−05).

Stratified SNP heritabilities and within-trait genetic correlations

Estimates of SNP heritability from stratified meta-analyses are shown in Figure 2 and Supplementary Table 8. In rater-specific meta-analyses, self-reported internalising symptoms showed the highest heritability (5.63%; 95% CI 3.08–8.18%), followed by teacher, maternal, and parental report, which were all significant. Although father-reported internalising symptoms had the highest SNP heritability in rater-specific analyses (8.98%), the wide confidence intervals overlapped zero (−0.06–18.02%). In age-specific meta-analyses, SNP h2 for internalising symptoms in adolescence was highest (1.97%, 95% CI 0.30– 3.64%), whereas estimates for early childhood, mid-childhood, and late childhood were similar, but not robust to significance testing. In rater-by-age meta-analyses, self-reported internalising symptoms during adolescence showed the highest SNP h2 (3.20%, 95% CI 0.34–6.06%). Instrument-specific meta-analyses showed that variance in internalising symptoms explained by ASEBA and SDQ scales were comparable, ∼3%. The estimate for Rutter was smaller (.3%), but the difference was not substantial, based on the overlapping confidence intervals.

Figure 2
Figure 2

SNP heritabilities based on N-weighted meta-analyses of internalising symptoms. Error bars represent 95% confidence intervals.

INToverall and self-reported internalising symptoms were highly genetically correlated (rg= 0.84, SE= 0.12, p=2.08×10−12). The other stratified meta-analyses were insufficiently powered to estimate genetic correlations (heritability z-score < 4).

Gene-based analysis, tissue expression and gene-set analyses

The genome-wide gene-based analysis did not reveal any genes significantly associated with internalising symptoms, but the top 10 genes are reported in Supplementary Table 10. MAGMA tissue expression analyses of 30 general and 53 specific tissue types did not show any statistically significant associations with internalising symptoms (Supplementary Table 11). The gene-set analysis did not show any significant associations (Supplementary Table 12).

Genetic correlations with external traits

Genetic correlations between INToverall and INTself (for which the z-score of the h2 was ≥ 4(37)), and a set of preselected external traits are shown in Figure 3 and Supplementary Table 13. INToverall held strong positive genetic correlations (rg > 0.7) with Major Depressive Disorder, anxiety, and neuroticism, and a strong negative correlation (rg < -0.7) with the well-being spectrum. High correlations (|rg| > 0.5) with other adult and childhood psychiatric and psychological traits, including attention-deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), depressive symptoms, loneliness, and overall and maternal-reported aggression were found. Moderate genetic correlations (|rg| > 0.3) with insomnia, age at first birth, cigarettes per day, educational attainment, and intelligence were also observed. INTself showed a similar pattern, but generally weaker genetic associations with external traits, with some exceptions. ASD, overall and maternal-reported aggression, age at first birth, and intelligence were correlated with INToverall, but showed weaker correlations with INTself, whereas self-reported aggression, smoking initiation and body-mass index (BMI) were correlated with INTself, but showed weaker or no correlation with INToverall.

Figure 3
Figure 3

Genetic correlations with external phenotypes. The left panel shows genetic correlations with the meta-analysis for overall int ernalising symptoms in childhood and adolescence (INToverall) and the right panel shows genetic correlations with self-reported internalising symptoms (INTself). Error bars represent 95% confidence intervals. Correlation points in red are statistically significant after correction for multiple testing.

Polygenic score prediction

Prediction of internalising symptoms in childhood and adolescence by polygenic scores based on INToverall are shown in Supplementary Table 14. After correction for multiple testing, polygenic scores for INToverall (Neff=132,260) were significantly associated with maternal-reported internalising problems in 7-year-olds, and explained up to 0.38% of the phenotypic variance. Polygenic scores for INToverall were not associated with self-reported internalising problems in adolescence.

Discussion

This genome-wide association meta-analysis of childhood and adolescent internalising symptoms included data from 64,641 individuals aged between 3 and 18. No genome-wide significant loci were detected at single SNP level. Genetic correlations with external traits suggested that childhood and adolescent internalising symptoms share substantial genetic vulnerabilities with adult internalising disorders and other childhood psychiatric traits, which could partially explain both the persistence of internalising symptoms over time and the high comorbidity amongst childhood psychiatric traits. The overall pattern of results pointed to rater-based heterogeneous effects (as discussed further on), which suggests that in addition to further increases of sample size, approaches that reduce heterogeneity will be essential in future GWAS investigations. As no biological pathways were identified, this study does not provide leads for new targets for treatment.

The most striking findings of this study are the direction and strength of genetic correlations with external traits, which point to an overlapping genetic architecture between internalising symptoms and other traits. This may initially be surprising given the low SNP heritability observed here, but while SNP heritability estimates the overall variance in a trait explained by genome-wide SNPs, a genetic correlation reflects the extent to which the same set of genetic factors are involved in two traits. As such, even traits with low SNP heritability can have high genetic correlations if the underlying set of genetic factors influencing the traits are overlapping. Strong genetic correlations (|rg| > 0.7) with adult depression, anxiety, neuroticism, and the wellbeing spectrum were of note, and suggest a substantial shared genetic etiology between childhood internalising symptoms and adult internalising disorders and related traits, that has also been observed in previous studies (41-43). Viewed in combination with the overlapping estimates of SNP heritability from early-childhood to adolescence in this study, these findings point to a stable set of genetic factors that partially explain the persistence of symptoms over time. Comparisons with other psychiatric disorders showed high genetic correlations (|rg| > 0.5) with childhood-onset disorders ADHD and ASD, but no robust associations with bipolar disorder, obsessive-compulsive disorder, or anorexia nervosa. A small genetic correlation with schizophrenia was observed (rg=0.2, p=.0025), which, albeit not significant due to the strict correction for multiple testing applied here, is in line with previous studies showing successful prediction of internalising symptoms in childhood using polygenic scores for schizophrenia (43-46). The overall pattern of genetic correlations with other psychiatric traits is comparable to adult cross-disorder genetic correlations, where depression shows stronger associations with ADHD and ASD than with schizophrenia or bipolar disorder (47). It appears that like adult depression, the broader (and perhaps also milder) symptomatology captured by dimensional measures of childhood internalising symptoms shares fewer genetic similarities with severe and less common disorders such as schizophrenia, bipolar disorder, OCD, and anorexia, but is more closely tied to childhood-onset disorders ADHD and ASD. This also resembles findings from the recent GWAS of total child psychiatric problems, which similarly found no robust genetic correlations with less common disorders (48). Correlations with other traits, including insomnia, loneliness, intelligence, educational attainment, cigarettes per day, and age at first birth were observed, as also seen in GWASs of adult depression and anxiety (23, 24), but unlike adult depression, no robust associations with coronary artery disease, BMI, smoking initiation, or age at menarche were found. On the other hand, both BMI and smoking initiation held robust associations with INTself, for which ratings were only available during adolescence. This could indicate that genetic factors during adolescence are particularly important in these associations. Age-specific genetic effects may also explain why coronary artery disease was not associated with INToverall, in contrast to the small but robust genetic correlation that coronary artery disease shares with both adult depression and anxiety (23, 24). This may point to genetic innovation (the involvement of novel genetic variants) in adulthood which could explain the genetic commonalities between adult internalising disorders and coronary artery disease. Alternatively, the lack of genetic correlation between INToverall and coronary artery disease, as well as age of menarche (which also genetically correlates with adult depression), could be due to a lack of power. This is also shown by the wide confidence intervals for some genetic correlations (Figure 3), which can be a consequence of low SNP heritability.

Focusing on childhood traits, as well as sharing high genetic correlations with childhood-onset disorders ADHD and ASD, internalising symptoms were also highly correlated with childhood aggression. The high correlations observed across childhood traits indicate the presence of specific genetic effects that are common between childhood disorders within the neurodevelopmental spectrum. These shared genetic effects could partially explain the high comorbidity between psychiatric traits in childhood (49-51). In further examining the association between childhood internalising symptoms and aggression, INToverall shared high genetic correlations with overall and maternal-reported aggression, but not with teacher or self-report. On the other hand, self-reported aggression and self-reported internalising symptoms were highly correlated, whereas INTself did not share robust associations with overall, teacher, or maternal reported aggression. These patterns of rater-stratified genetic correlations suggest that observed genetic effects on childhood phenotypes can vary substantially due to differences in the phenotype captured by different raters, with the same set of raters showing the highest correlation between traits.

The difficulty in identifying causal loci for early-life internalising symptoms is not novel and resembles the trajectory of GWAS investigations of adult internalising disorders. GWAS studies of adult depression also made slow progress due to limited sample sizes and heterogeneity (52-54). As depression has several potential sources of heterogeneity, including a diverse presentation of symptoms, large case-control sample sizes were required to achieve success in identifying specific genomic loci (23, 25). GWAS studies of anxiety similarly saw increased success as sample sizes grew (24, 55). However, here we found that a substantial increase in sample size in comparison to previous studies (15, 16) did not yield the expected increase in power. In addition to heterogeneity due to a broad symptomatology, our findings indicate that GWAS investigations of childhood internalising symptoms are further disadvantaged by rater-based heterogeneous effects. Unlike adult studies where measurements are typically self or clinician reports, childhood studies, particularly those focusing on early childhood, rely heavily on parent and teacher report, which act as an additional source of heterogeneity. Rater-based differences in genetic correlations with external traits have been discussed above. The current study also observed varying estimates of SNP-heritability in rater-stratified analyses (Figure 2). Although these estimates did not appear significantly different (likely due to sample size limitations), the partial genetic correlation between INToverall and INTself points to incomplete overlap in relevant SNPs, indicating the presence of rater-specific genetic effects. Additionally, polygenic scores based on INToverall did not predict self-reported internalising symptoms in the NTR cohort, which could also indicate heterogeneity between the target and discovery traits (56). Rater-specific genetic effects and rater disagreement on internalising symptoms are noted in previous research (57-60) and rater-based heterogeneity is also reported in the GWAS of childhood aggression (35). Heterogeneity in this study is further reflected in the low LDSC-estimate of SNP heritability in the INToverall meta-analysis. This can partly be explained by the methods; estimates of SNP heritability from summary statistics are typically lower than estimates from raw genotypic data and potential over-correction of biases in LDSC could have led to more conservative estimates. However, heritability estimates rose when analyses were stratified, indicating that heterogeneity also diluted overall SNP effects in the meta-analysis.

Heterogeneous effects underlying childhood internalising symptoms can be accounted for in multivariate GWAS approaches, but our study shows that current childhood samples seem unable to meet the power requirements of these types of analyses. Another way of reducing heterogeneity and helping signal detection is to focus on diagnoses. The case-control approach has proven to be more successful than dimensional measures in adult studies of depression and anxiety (23, 24) and overcomes the limitations of treating symptom scales as continuous traits. However, diagnostic data are currently not available for childhood phenotypes in large enough samples. Instead, we expect that reducing heterogeneity at phenotypic level will be key in paving the way to success in future GWAS investigations in childhood samples. This could be tackled by examining symptom level phenotypes or separating childhood anxiety and depression into two distinct phenotypes. Another promising approach would be to eliminate heterogenous effects entirely through factor analysis. Factor analysis can be used to derive a stable phenotype which captures the core behaviour that multiple measurements (e.g. from different informants or at different time points) have in common. This eliminates variability from rater, age, or situational effects. Evidence from both twin and molecular research shows that focusing on the common part of multiple assessments results in a more reliable phenotype which shows higher heritability than that captured by individual measurements separately (41, 59, 61, 62). This way of managing rater bias has broader applicability in genetic studies within child psychiatry, but is dependent on the availability of multiple informants on behaviour at one time point.

To conclude, in this large GWAS of childhood and adolescent internalising symptoms in population-based cohorts of European ancestry, no individual loci with strong associations with the outcome were detected. However, strong genetic correlations with adult internalising traits and childhood psychiatric traits indicate that there is signal buried in the noise. Future GWAS success is likely to lie in reducing heterogeneity in childhood samples by focusing on a more stable phenotype of internalising symptoms. Finally, an important goal for future GWASs is the funding and inclusion of multi-ancestry cohorts to allow better representation of diverse populations and ensure broader applicability of findings.

Data Availability

Summary statistics are available upon request.

https://osf.io/edas6

https://www.eagle-consortium.org/

https://www.ebi.ac.uk/gwas/

https://www.omim.org

Footnotes

  • Gene-based analysis updated, author details updated, supplementary files updated

References

  1. 1.↵
    Erskine H, Moffitt TE, Copeland W, Costello E, Ferrari A, Patton G, et al. A heavy burden on young minds: the global burden of mental and substance use disorders in children and youth. Psychological medicine. 2015;45(7):1551–63.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Organization WH. Depression and other common mental disorders: global health estimates. World Health Organization; 2017.
  3. 3.↵
    Steel Z, Marnane C, Iranpour C, Chey T, Jackson JW, Patel V, et al. The global prevalence of common mental disorders: a systematic review and meta-analysis 1980–2013. International journal of epidemiology. 2014;43(2):476–93.
    OpenUrlCrossRefPubMed
  4. 4.↵
    Polanczyk GV, Salum GA, Sugaya LS, Caye A, Rohde LA. Annual Research Review: A meta-analysis of the worldwide prevalence of mental disorders in children and adolescents. Journal of Child Psychology and Psychiatry. 2015;56(3):345–65.
    OpenUrlCrossRefPubMed
  5. 5.↵
    Husky MM, Boyd A, Bitfoi A, Carta MG, Chan-Chee C, Goelitz D, et al. Self-reported mental health in children ages 6–12 years across eight European countries. European child & adolescent psychiatry. 2017:1–11.
  6. 6.↵
    Weissman MM, Wolk S, Wickramaratne P, Goldstein RB, Adams P, Greenwald S, et al. Children with prepubertal-onset major depressive disorder and anxiety grown up. Archives of general psychiatry. 1999;56(9):794–801.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.
    Roza SJ, Hofstra MB, van der Ende J, Verhulst FC. Stable prediction of mood and anxiety disorders based on behavioral and emotional problems in childhood: A 14-year follow- up during childhood, adolescence, and young adulthood. American Journal of Psychiatry. 2003;160(12):2116–21.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    Gregory AM, Eley TC. Genetic Influences on Anxiety in Children: What we’ve Learned and Where we’re Heading. Clinical Child and Family Psychology Review. 2007;10(3):199–212.
    OpenUrlCrossRefPubMed
  9. 9.↵
    Franic S, Dolan CV, Borsboom D, van Beijsterveldt CE, Boomsma DI. Three-and-a-half- factor model? The genetic and environmental structure of the CBCL/6–18 internalizing grouping. Behavior genetics. 2014;44(3):254–68.
    OpenUrl
  10. 10.↵
    Nivard M, Dolan C, Kendler K, Kan K-J, Willemsen G, van Beijsterveldt C, et al. Stability in symptoms of anxiety and depression as a function of genotype and environment: a longitudinal twin study from ages 3 to 63 years. Psychological medicine. 2015;45(5):1039–49.
    OpenUrlCrossRefPubMed
  11. 11.↵
    Polderman TJ, Benyamin B, De Leeuw CA, Sullivan PF, Van Bochoven A, Visscher PM, et al. Meta-analysis of the heritability of human traits based on fifty years of twin studies. Nature genetics. 2015;47(7):702.
    OpenUrlCrossRefPubMed
  12. 12.↵
    Hannigan L, Walaker N, Waszczuk M, McAdams T, Eley T. Aetiological influences on stability and change in emotional and behavioural problems across development: a systematic review. Psychopathology review. 2017;4(1):52.
    OpenUrl
  13. 13.
    Franic S, Middeldorp CM, Dolan CV, Ligthart L, Boomsma DI. Childhood and Adolescent Anxiety and Depression: Beyond Heritability. Journal of the American Academy of Child & Adolescent Psychiatry. 2010;49(8):820–9.
    OpenUrl
  14. 14.↵
    Kendler KS, Gardner CO, Lichtenstein P. A developmental twin study of symptoms of anxiety and depression: evidence for genetic innovation and attenuation. Psychological Medicine. 2008;38(11):1567–75.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    Trzaskowski M, Eley TC, Davis OS, Doherty SJ, Hanscombe KB, Meaburn EL, et al. First genome-wide association study on anxiety-related behaviours in childhood. PloS one. 2013;8(4):e58676.
    OpenUrl
  16. 16.↵
    Benke KS, Nivard MG, Velders FP, Walters RK, Pappa I, Scheet PA, et al. A genome- wide association meta-analysis of preschool internalizing problems. Journal of the American Academy of Child & Adolescent Psychiatry. 2014;53(6):667-76. e7.
    OpenUrl
  17. 17.↵
    Trzaskowski M, Dale PS, Plomin R. No genetic influence for childhood behavior problems from DNA analysis. Journal of the American Academy of Child & Adolescent Psychiatry. 2013;52(10):1048-56. e3.
    OpenUrl
  18. 18.
    Pappa I, Fedko IO, Mileva-Seitz VR, Hottenga J-J, Bakermans-Kranenburg MJ, Bartels M, et al. Single nucleotide polymorphism heritability of behavior problems in childhood: genome-wide complex trait analysis. Journal of the American Academy of Child & Adolescent Psychiatry. 2015;54(9):737–44.
    OpenUrl
  19. 19.
    Sallis H, Evans J, Wootton R, Krapohl E, Oldehinkel AJ, Smith GD, et al. Genetics of depressive symptoms in adolescence. BMC psychiatry. 2017;17(1):321.
    OpenUrl
  20. 20.
    Cheesman R, Selzam S, Ronald A, Dale PS, McAdams TA, Eley TC, et al. Childhood behaviour problems show the greatest gap between DNA-based and twin heritability. Translational psychiatry. 2017;7(12):1284.
    OpenUrl
  21. 21.
    Jami ES, Eilertsen EM, Hammerschlag AR, Qiao Z, Evans DM, Ystrom E, et al. Maternal and paternal effects on offspring internalizing problems: Results from genetic and family-based analyses. American Journal of Medical Genetics Part B-Neuropsychiatric Genetics. 2020;183(5):258–67.
    OpenUrl
  22. 22.↵
    Cheesman R, Eilertsen EM, Ahmadzadeh YI, Gjerde LC, Hannigan LJ, Havdahl A, et al. How important are parents in the development of child anxiety and depression? A genomic analysis of parent-offspring trios in the Norwegian Mother Father and Child Cohort Study (MoBa). medRxiv. 2020:2020.04.14.20064782.
  23. 23.↵
    Howard DM, Adams MJ, Clarke T-K, Hafferty JD, Gibson J, Shirali M, et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nature neuroscience. 2019;22(3):343.
    OpenUrlPubMed
  24. 24.↵
    Purves KL, Coleman JR, Meier SM, Rayner C, Davis KA, Cheesman R, et al. A major role for common genetic variation in anxiety disorders. Molecular psychiatry. 2019:1–12.
  25. 25.↵
    Wray NR, Ripke S, Mattheisen M, Trzaskowski M, Byrne EM, Abdellaoui A, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nature genetics. 2018;50(5):668.
    OpenUrlCrossRefPubMed
  26. 26.↵
    Middeldorp CM, Felix JF, Mahajan A, McCarthy MI, consortium EGG. The Early Growth Genetics (EGG) and EArly Genetics and Lifecourse Epidemiology (EAGLE) consortia: design, results and future prospects. European journal of epidemiology. 2019;34(3):279–300.
    OpenUrlCrossRef
  27. 27.↵
    Van der Sluis S, Posthuma D, Nivard M, Verhage M, Dolan C. Power in GWAS: lifting the curse of the clinical cut-off. Molecular psychiatry. 2013;18(1):2–3.
    OpenUrlCrossRefPubMed
  28. 28.↵
    Martin J, Khramtsova EA, Goleva SB, Blokland GA, Traglia M, Walters RK, et al. Examining sex-differentiated genetic effects across neuropsychiatric and behavioral traits. Biological psychiatry. 2021.
  29. 29.↵
    Goodman R. The Strengths and Difficulties Questionnaire: a research note. Journal of child psychology and psychiatry. 1997;38(5):581–6.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    Achenbach TM. The Achenbach system of empirically based assessment (ASEBA): Development, findings, theory, and applications: University of Vermont, Research Center for Children, Youth, & Families; 2009.
  31. 31.↵
    Rutter M. A children’s behaviour questionnaire for completion by teachers: preliminary findings. Journal of child Psychology and Psychiatry. 1967;8(1):1–11.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    Rutter M, Tizard J, Whitmore K. Education, Health and Behaviour. London: Longman; 1970.
  33. 33.↵
    Winkler TW, Day FR, Croteau-Chonka DC, Wood AR, Locke AE, Mägi R, et al. Quality control and conduct of genome-wide association meta-analyses. Nature protocols. 2014;9(5):1192.
    OpenUrl
  34. 34.↵
    Baselmans BM, Jansen R, Ip HF, van Dongen J, Abdellaoui A, van de Weijer MP, et al. Multivariate genome-wide analyses of the well-being spectrum. Nature genetics. 2019;51(3):445–51.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Ip HF, van der Laan CM, Brikell I, Sánchez-Mora C, Nolte IM, St Pourcain B, et al. Genetic Association Study of Childhood Aggression across raters, instruments and age. bioRxiv. 2019:854927.
  36. 36.↵
    Bulik-Sullivan BK, Loh P-R, Finucane HK, Ripke S, Yang J, Patterson N, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nature genetics. 2015;47(3):291–5.
    OpenUrlCrossRefPubMed
  37. 37.↵
    Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh P-R, et al. An atlas of genetic correlations across human diseases and traits. Nature genetics. 2015;47(11):1236.
    OpenUrlCrossRefPubMed
  38. 38.↵
    de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: generalized gene-set analysis of GWAS data. PLoS computational biology. 2015;11(4).
  39. 39.↵
    Watanabe K, Taskesen E, Van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nature communications. 2017;8(1):1–11.
    OpenUrl
  40. 40.↵
    Vilhjálmsson Bjarni J, Yang J, Finucane Hilary K, Gusev A, Lindström S, Ripke S, et al. Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores. The American Journal of Human Genetics. 2015;97(4):576–92.
    OpenUrlCrossRefPubMed
  41. 41.↵
    Cheesman R, Purves KL, Pingault J-B, Breen G, Plomin R, Eley TC. Extracting stability increases the SNP heritability of emotional problems in young people. Translational psychiatry. 2018;8(1):1–9.
    OpenUrl
  42. 42.
    Akingbuwa WA, Hammerschlag AR, Jami ES, Allegrini AG, Karhunen V, Sallis H, et al. Genetic Associations Between Childhood Psychopathology and Adult Depression and Associated Traits in 42L998 Individuals: A Meta-Analysis. JAMA Psychiatry. 2020.
  43. 43.↵
    Musci RJ, Masyn KE, Benke K, Maher B, Uhl G, Ialongo NS. The effects of the interplay of genetics and early environmental risk on the course of internalizing symptoms from late childhood through adolescence. Development and psychopathology. 2016;28(1):225–37.
    OpenUrl
  44. 44.
    Nivard MG, Gage SH, Hottenga JJ, van Beijsterveldt CE, Abdellaoui A, Bartels M, et al. Genetic Overlap Between Schizophrenia and Developmental Psychopathology: Longitudinal and Multivariate Polygenic Risk Prediction of Common Psychiatric Traits During Development. Schizophrenia Bulletin. 2017: sbx031.
  45. 45.
    Riglin L, Collishaw S, Richards A, Thapar AK, Rice F, Maughan B, et al. The impact of schizophrenia and mood disorder risk alleles on emotional problems: investigating change from childhood to middle age. Psychological medicine. 2017:1–6.
  46. 46.↵
    Jansen PR, Polderman TJ, Bolhuis K, Ende J, Jaddoe VW, Verhulst FC, et al. Polygenic scores for schizophrenia and educational attainment are associated with behavioural problems in early childhood in the general population. Journal of Child Psychology and Psychiatry. 2018;59(1):39–47.
    OpenUrlCrossRef
  47. 47.↵
    Lee PH, Anttila V, Won H, Feng Y-CA, Rosenthal J, Zhu Z, et al. Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders. Cell. 2019;179(7):1469-82. e11.
    OpenUrlCrossRef
  48. 48.↵
    Neumann A, Nolte IM, Pappa I, Ahluwalia TS, Pettersson E, Rodriguez A, et al. A genome-wide association study of total child psychiatric problems scores. medRxiv. 2020:2020.06.04.20121061.
  49. 49.↵
    Rhee SH, Lahey BB, Waldman ID. Comorbidity Among Dimensions of Childhood Psychopathology: Converging Evidence From Behavior Genetics. Child Development Perspectives. 2015;9(1):26–31.
    OpenUrl
  50. 50.
    Allegrini AG, Cheesman R, Rimfeld K, Selzam S, Pingault J-B, Eley TC, et al. The p factor: genetic analyses support a general dimension of psychopathology in childhood and adolescence. Journal of Child Psychology and Psychiatry. 2020;61(1):30–9.
    OpenUrlCrossRef
  51. 51.↵
    Bartels M, Hendriks A, Mauri M, Krapohl E, Whipp A, Bolhuis K, et al. Childhood aggression and the co-occurrence of behavioural and emotional problems: results across ages 3– 16 years from multiple raters in six cohorts in the EU-ACTION project. European child & adolescent psychiatry. 2018;27(9):1105–21.
    OpenUrl
  52. 52.↵
    Levinson DF, Mostafavi S, Milaneschi Y, Rivera M, Ripke S, Wray NR, et al. Genetic studies of major depressive disorder: Why are there no GWAS findings, and what can we do about it? Biological psychiatry. 2014;76(7):510.
    OpenUrlCrossRefPubMed
  53. 53.
    Wray NR, Maier R. Genetic Basis of Complex Genetic Disease: The Contribution of Disease Heterogeneity to Missing Heritability. Current Epidemiology Reports. 2014;1(4):220–7.
    OpenUrl
  54. 54.↵
    Cai N, Choi KW, Fried EI. Reviewing the genetics of heterogeneity in depression: operationalizations, manifestations and etiologies. Human Molecular Genetics. 2020.
  55. 55.↵
    Otowa T, Hek K, Lee M, Byrne EM, Mirza SS, Nivard MG, et al. Meta-analysis of genome-wide association studies of anxiety disorders. Molecular psychiatry. 2016;21(10):1391.
    OpenUrlCrossRefPubMed
  56. 56.↵
    de Vlaming R, Okbay A, Rietveld CA, Johannesson M, Magnusson PKE, Uitterlinden AG, et al. Meta-GWAS Accuracy and Power (MetaGAP) Calculator Shows that Hiding Heritability Is Partially Due to Imperfect Genetic Correlations across Studies. PLoS Genet. 2017;13(1):e1006495–e.
    OpenUrl
  57. 57.↵
    Ask H, Torgersen S, Seglem KB, Waaktaar T. Genetic and environmental causes of variation in adolescent anxiety symptoms: A multiple-rater twin study. Journal of anxiety disorders. 2014;28(4):363–71.
    OpenUrlCrossRef
  58. 58.
    Fedko IO, Hottenga J-J, Medina-Gomez C, Pappa I, van Beijsterveldt CE, Ehli EA, et al. Estimation of genetic relationships between individuals across cohorts and platforms: application to childhood height. Behavior genetics. 2015;45(5):514–28.
    OpenUrl
  59. 59.↵
    Wesseldijk LW, Fedko IO, Bartels M, Nivard MG, van Beijsterveldt CE, Boomsma DI, et al. Psychopathology in 7-year-old children: Differences in maternal and paternal ratings and the genetic epidemiology. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics. 2017;174(3):251–60.
    OpenUrl
  60. 60.↵
    Bartels M, Boomsma DI, Hudziak JJ, van Beijsterveldt TC, van den Oord EJ. Twins and the study of rater (dis) agreement. Psychological methods. 2007;12(4):451.
    OpenUrlPubMed
  61. 61.↵
    Lubke GH, Miller PJ, Verhulst B, Bartels M, van Beijsterveldt T, Willemsen G, et al. A powerful phenotype for gene-finding studies derived from trajectory analyses of symptoms of anxiety and depression between age seven and 18. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics. 2016;171(7):948–57.
    OpenUrlCrossRef
  62. 62.↵
    Fedko IO, Wesseldijk LW, Nivard MG, Hottenga J-J, van Beijsterveldt CE, Middeldorp CM, et al. Heritability of Behavioral Problems in 7-Year Olds Based on Shared and Unique Aspects of Parental Views. Behavior genetics. 2017;47(2):152–63.
    OpenUrl
Back to top
Posted July 31, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Genome-wide association meta-analysis of childhood and adolescent internalising symptoms
Eshim S Jami, Anke R Hammerschlag, Hill F Ip, Andrea G Allegrini, Beben Benyamin, Richard Border, Elizabeth W Diemer, Chang Jiang, Ville Karhunen, Yi Lu, Qing Lu, Travis T Mallard, Pashupati P Mishra, Ilja M Nolte, Teemu Palviainen, Roseann E Peterson, Hannah M Sallis, Andrey A Shabalin, Ashley E Tate, Elisabeth Thiering, Natàlia Vilor-Tejedor, Carol Wang, Ang Zhou, Daniel E Adkins, Silvia Alemany, Helga Ask, Qi Chen, Robin P Corley, Erik A Ehli, Luke M Evans, Alexandra Havdahl, Fiona A Hagenbeek, Christian Hakulinen, Anjali K Henders, Jouke Jan Hottenga, Tellervo Korhonen, Abdullah Mamun, Shelby Marrington, Alexander Neumann, Kaili Rimfeld, Fernando Rivadeneira, Judy L Silberg, Catharina E van Beijsterveldt, Eero Vuoksimaa, Alyce M Whipp, Tong Xiaoran, Ole A Andreassen, Dorret Boomsma, Sandra A Brown, S Alexandra Burt, William Copeland, Elizabeth J Costello, Danielle M Dick, Lindon J Eaves, K Paige Harden, Kathleen Mullan Harris, Catharina A Hartman, Joachim Heinrich, John K Hewitt, Christian Hopfer, Elina Hypponen, Marjo-Riitta Jarvelin, Jaakko Kaprio, Liisa Keltikangas-Järvinen, Kelly L Klump, Kenneth Krauter, Ralf Kuja-Halkola, Henrik Larsson, Terho Lehtimäki, Paul Lichtenstein, Sebastian Lundstrom, Hermine H Maes, Per Magnus, Marcus R Munafò, Jake M Najman, Pål R Njølstad, Albertine J Oldehinkel, Craig E Pennell, Robert Plomin, Ted Reichborn-Kjennerud, Chandra Reynolds, Richard J Rose, Andrew Smolen, Harold Snieder, Michael Stallings, Marie Standl, Jordi Sunyer, Henning Tiemeier, Sally Wadsworth, Tamara L Wall, Andrew J O Whitehouse, Gail M Williams, Eivind Ystrom, Michel G Nivard, Meike Bartels, Christel M Middeldorp
medRxiv 2020.09.11.20175026; doi: https://doi.org/10.1101/2020.09.11.20175026
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area