Introduction

Prevention strategies and behavioural and pharmacological interventions for aggressive behaviour and conduct disorder are effective in some children, although a substantial number of children do not respond to prevention strategies, do not benefit from interventions or may even experience an escalation of symptom [9, 10]. One reason for this might be the heterogeneity of aggression. A second reason, which is related to the heterogeneous nature and occurrence of childhood aggressive problems, might be that children with aggressive problems often have co-occurring problems. Due to a multitude of problems, children may not respond to prevention or intervention targeting aggression, or the co-occurring problems may mask aggression, leaving it untreated. In 12 year olds, Bartels and colleagues [11] observed that at least half of the children who were deviant on aggressive behaviour (T score ≥ 67) also were deviant on rule-breaking behaviour, i.e. at least 50% of the children with clinical levels of aggression also showed a co-occurrence of clinically relevant rule-breaking behaviour. Strong links between aggression and attention-deficit/hyperactivity disorder (ADHD) [12] are often seen in the clinical presentation of ADHD [13], and it has been suggested that the strong association between ADHD and aggression may explain gender differences in clinical referral. For example, teachers rated boys with a DSM-based ADHD diagnosis as having higher levels of attention problems and aggression than girls with a similar ADHD diagnosis [14]. Aggression not only co-occurs with psychopathologies on the externalizing spectrum. Aggression also tends to co-occur with anxiety, and it has been proposed that anxiety needs to be given a central role in the treatment of aggression [15]. In more extreme cases, aggression was not found to co-occur solely with ADHD symptoms, such as attention problems, or anxiety but rather with both of these forms of psychopathology. This pattern of behavioural problems is referred to as the dysregulation profile [16,17,18], and has been described as a potential marker for severe childhood psychopathologies [19, 20].

To gain insight into the aetiology of individual differences in childhood aggression and in co-occurring behavioural and emotional problems, ACTION (Aggression in Children: unravelling gene-environment interplay to inform Treatment and InterventiON strategies; http://www.action-euproject.eu/) created a consortium with access to large childhood prospective twin, population-based and clinical cohorts. ACTION brings together multiple large cohort studies in genetically informative populations (see Table 1 and Appendix 1). The focus of ACTION is to inform on the aetiology of differences in aggression between children by unravelling its genetic architecture using univariate, multivariate and longitudinal genetic and epigenetic modelling in twin and genetic and epigenetic association studies. A strong focus of ACTION includes biomarker and metabolomics research [21].

Table 1 Sample sizes for different age groups of the ACTION cohort

In the current study, the aim is to describe and improve the understanding of the co-occurrence of aggression with other forms of childhood psychopathology by analysing data from the large ACTION phenotype databases in large samples of children. We analysed data on aggression and common emotional and behavioural problems in children aged 3–16 years. Multiple raters, i.e. fathers and mothers during childhood and also youngsters themselves during adolescence, provided information on different aggression measures. The two Dutch cohorts (The Netherlands Twin Register and Generation R) used the Achenbach System of Empirically Based Assessment (ASEBA [22]), which included the Child Behaviour Checklist (CBCL) and the Youth Self-Report (YSR). The UK-based Twins Early Development Study employed the Strengths and Difficulties Questionnaire (SDQ [23]). The Swedish Twin study of Child and Adolescent Development used the Autism–Tics, ADHD and other Comorbidities inventory (A-TAC [24]), and the Swedish Child and Adolescent Twin Study the ASEBA questionnaires. In Finland, the Multidimensional Peer Nomination Inventory (MPNI) was employed. For several age groups from different countries, aggression assessed with identical instruments was available. For example, parental ratings with the CBCL were available for 7- to 8-year-olds and 12- to 13-year-olds in the Netherlands (NTR) and Sweden (TCHAD). In addition to indicators of aggression, all instruments provided quantitative scores on other childhood psychopathologies from the externalising and internalising spectrum. We investigated patterns of co-occurrence across age, rater, instrument and gender.

Methods

Participants

Six large population-based cohorts (NTR and GenR from the Netherlands, TEDS from the UK, CATSS and TCHAD from Sweden and FinnTwin12 from Finland) analysed the co-occurrence of aggression measures with other psychopathologies. For a link to cohort-specific websites, see Table 1 and for a detailed description of the cohorts, please also see Appendix I. The twin cohorts were requested to randomly select one of the twins per pair, with an equal representation of first- and second-born children, to obtain parameter estimates that were not biased due to effects of family clustering. In our previous work [25], we have shown that children with an illness or disability that interfered with daily function tend to display more than twice as much problem behaviour across the entire age range compared to other twins, so they were excluded. Age-, gender- and rater- specific sample sizes are presented in Tables 2, 3, 4 and 5. Data were available for 30,523 3- to 4-year-olds (49.5% boys), 20,958 5- to 6-year-olds (49.6% boys), 18,291 7- to 8-year-olds 49% boys), 27,218 9- to 10-year-olds (49.4% boys), 18,543 11- to 12-year-olds (48.9% boys) and 10,088 15- to 16-year- olds (46.6% boys). Due to the longitudinal structure of most cohorts, these data points are not statistically independent observations, since overlapping groups of children were assessed at multiple ages. All data used in the current analyses were collected under protocols that have been approved by the appropriate ethics committees, and studies were performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.

Table 2 Means and standard deviations for the empirical scales of the Child Behaviour Checklist (CBCL) 1.5–5
Table 3 Means and standard deviations for the empirical scales of the ASEBA taxonomy (CBCL and YSR)
Table 4 Means and standard deviations for the scales of the Strengths and Difficulties Questionnaire (SDQ)
Table 5 Means and standard deviations for the scales of the Devereux Child Behaviour rating scale (DCB), Autism–Tics, ADHD and other Comorbidities inventory (A-TAC) and Multidimensional Peer Nomination Inventory (MPNI)

Measures

The Child Behaviour Checklist (CBCL) 1,5–5 [26] and 6–18 [22] were used by GenR (age 6 and 10), TCHAD (ages 8, 13 and 16) and NTR (ages 3, 7, 9 and 12). The Youth Self-Report (YSR) [22] was used by TCHAD (ages 13 and 16). The CBCL and YSR are part of the Achenbach System of Empirical-Based Assessment and designed to measure childhood and adolescent behavioural and emotional problems. The response format was on a three-point scale (with response options ‘not true’, ‘somewhat true or sometimes true’ and ‘very or often true’). With the CBCL 1,5–5 seven syndrome scales are obtained (emotionally reactive, anxious–depressed, somatic complaints, withdrawn, overactive behaviour, aggressive behaviour, sleep problems), while with the CBCL 6–18 eight syndrome scales are obtained (anxious–depressed, withdrawn, somatic complaints, social problems, thought problems, attention problems, rule-breaking behaviour, aggressive behaviour). With the YSR, eight syndrome scales are obtained (anxious–depressed, somatic complaints, withdrawn–depressed, social problems, thought problems, attention problems, rule-breaking behaviour and aggressive behaviour).

The Strengths and Difficulties Questionnaire (SDQ) [23] was used by TEDs (ages 4, 7, 9, 16) and CATSS (age 15). The SDQ is a 25-item questionnaire designed to measure common mental health problems during childhood and adolescence. Ratings were on a three-point scale (with response options ‘not true’, ‘somewhat true’ and ‘certainly true’). The 25 items form 5 scales, emotional symptoms, conduct problems, hyperactivity/inattention, peer relationship problem and prosocial behaviour. The conduct problem scale was used as a proxy for aggressive behaviour.

NTR used the short Devereux Child Behaviour (DCB) rating scale for 5 year olds. The DCB consists of questions about problem behaviour in children rated by the parents [27]. The short version includes 42 items that measure seven different aspects of problem behaviour in children. Parents were asked to indicate on a five-point scale whether the statements were applicable (0 = never, 1 = rarely, 2 = occasionally, 3 = frequently, 4 = very frequently). The items of the questionnaire cover the following aspects of problem behaviour: emotional liability (five items, e.g. “markedly impatient”), social isolation (three items, e.g. “quite timid or shy”), aggressive behaviour (seven items, e.g. “hits, bites and scratches other children”), attention problems (five items, e.g. “jumps from one activity to another”), dependency (five items, e.g. “does not want to do things for himself”), anxiety problems (six items, e.g. “concern about his physical health”) and physical coordination (five items, e.g. “gets dirty and untidy”).

In 9- and 12-year-old in the CATSS sample from Sweden, information on ODD/CD and other psychopathologies was gathered through a telephone interview with parents, using The AutismTics, ADHD and other Comorbidities inventory (A-TAC) [24]. A-TAC is a comprehensive screening interview for autism spectrum disorders (ASDs), attention-deficit/hyperactivity disorder (AD/HD), tic disorders (TD), developmental coordination disorder (DCD), learning disorders (LD) and other childhood mental disorders that have been associated with these neurodevelopmental disorders.

In the FinnTwin12 sample from Finland, aggressive behaviour was assessed at ages 12, 14 and 17 by versions of the Multidimensional Peer Nomination Inventory (MPNI). The MPNI includes 37 items comprising three subscales, the two subscales used here include: externalising behavioural problems (aggression, hyperactivity–impulsivity and inattention) and internalising emotional problems (anxiety and depression) [28]. For each question (e.g. ‘Does the child tease smaller or weaker children?’), the informant rated how well the description fit the twin in question on a scale from 0 (the characteristic does not fit the child at all) to 3 (the characteristic fits the child very well). Parents rated the children at age 12, and the child rated him or herself at ages 14 and 17 years.

Analyses

To ensure homogenous handling of data and identical analyses, all cohorts received a standard operating procedure that specified details of the comorbidity analyses. Following the SOP average scores and Pearson correlations for aggression with all other scales assessing psychopathology was obtained by a local analyst using their preferred statistical software. Average scores and correlations were computed by gender and age of children, separately for each rater and country. Results were uploaded to a shared server. Given the large datasets included in these analyses, leading to significance even if differences between average scores or between correlations being relatively small, we interpreted all results relative to each other and took the 95% confidence intervals into account. With the multi-instrument, multi-rater and multi-age assessments of aggression and of other emotional and behavioural problems, we established whether co-occurrence was stronger or weaker given different measurement instruments, raters and ages.

Results

Tables 2, 3, 4 and 5 provide an overview of the sample sizes and mean levels of aggression and all other traits. We replicated the well-established gender differences in average aggression scores at most ages for parental ratings. The gender difference was smaller or close to absent for self-reports. For example, while the difference between boys and girls is in general about 1.5–2 points on the CBCL and SDQ parental reports, the differences based on self-report ranged between 0.05 and 0.67.

Mean levels based on similar instruments across countries were almost identical. For example, the mean level of aggression based on maternal ratings of 7-year-old boys in the Netherlands was 5.74 (SD 5.29), while mean level of aggression based on parental ratings of 8-year-old boys in Sweden was 5.49 (SD 5.42).

We observed differences between raters in nearly every country in the same direction. Based on maternal ratings, higher levels of psychopathology were seen than when based on paternal ratings. These differences were observed both for boys and girls, at ages 3–12 for the CBCL and SDQ. The exception was an absence of differences in maternal and paternal ratings when using the Devereux Child Behaviour rating scale.

With respect to our main question of the co-occurrence of aggression and other behavioural and emotional problems, findings are presented in Tables 6, 7, 8 and 9. Strong correlations were found between aggression and other externalising traits, especially rule-breaking behaviour. Correlations of almost similar strength were also observed for aggression and attention problems and hyperactivity. However, correlations were lower between aggression and internalising behaviours including withdrawn–depression and somatic complaints. Correlations between aggression and all other emotional and behavioural problems and their 95% confidence intervals are also provided in an interactive application which can be found at http://www.action-euproject.eu/ComorbidityChildAggression.

Table 6 Phenotypic correlations between aggression and other empirical scales of the Child Behaviour Checklist (CBCL) 1.5-5
Table 7 Phenotypic correlations between aggression and other empirical scales of the achenbach system of empirically based assessment (ASEBA) 6–18
Table 8 Phenotypic correlations between aggression and other scales of the Strengths and Difficulties Questionnaire (SDQ)
Table 9 Phenotypic correlations between aggression and other scales of the Devereux Child Behaviour rating scale (DCB), Autism–Tics, ADHD and other Comorbidities inventory (A-TAC) and multidimensional peer nomination inventory (MPNI)

Some more remarkable findings included the relatively low correlation between aggression and obsessive–compulsive behaviour and the similarly relatively low correlation between aggression and social isolation and aggression and dependency. We, furthermore, observed a relatively low correlation between aggression and peer problems from the SDQ (ranging from 0.18 to 31). However, CBCL social problems showed stronger correlations with aggression (ranging from 0.34 to 0.66).

The overarching picture that emerged suggests that correlations are largely stable across rater and age. If any change is observed, it is indicative of stronger correlations when children grow older. The correlations patterns of boys are markedly similar to the correlational patterns of girls. The only exception was the ATAC-based correlation between ODD/CD and OCD based on parental ratings at age 12. Correlations were stronger when based on the CBCL in comparison to the other measures, especially for parental ratings, while the ATAC, which is a clinical interview rather than a survey, provided somewhat lower correlations. The Devereux Child Behaviour (DCB) rating scale provides the interesting finding of similar strength in correlations between aggressive behaviour and attention problems and anxiety problems, but also with physical coordination problems.

Discussion

One of the aims of ACTION is to describe and improve the understanding of the co-occurrence of aggression with other forms of childhood psychopathology. Here, we presented the correlations of aggression with other psychopathologies in large European samples of children between ages 3 and 16 years old. We showed that aggression co-occurred with almost all other behavioural and social problems. More specifically, aggression co-occurred with oppositional and ADHD-related problems, and at later ages with rule-breaking. In addition to the high correlations of aggression with externalising problems, we also observed substantial associations with anxiety–depression and other internalising symptoms. This co-occurrence of internalising and externalising problems has previously been shown to persist over childhood and adolescence [29]. Both for externalising and internalising problems, the patterns of co-occurrence were largely gender and rater independent, and were similar even when aggression and the other psychopathologies were assessed by different instruments, such as the CBCL and the SDQ. Also, there were large similarities in co-occurrence patterns across countries in the Northern part of Europe.

In ACTION, we compared co-occurrence patterns across different countries and cultures. These comparisons are somewhat hampered by the fact that in almost all cases more than one parameter varies between the different countries and cultures. For example, both NTR, a Dutch sample, and TEDs, a UK sample, have parental ratings at age 9, but NTR used the CBCL while TEDS used the SDQ. Any differences in correlations may thus be attributable to cultural differences or country differences between the Netherlands and the UK, instrument differences or any other protocol or unobserved difference. However, given all these sources of difference in this large co-occurrence study, it is even more striking that most correlations are so similar.

The large associations of aggression with other emotional and behavioural problems may form one of the obstacles for prevention and treatment of aggression. These findings indicate that an exclusive focus on aggression might not be the most feasible approach for the development of effective prevention and intervention programs. The complexity of psychopathology, partly due to the co-occurrence of behavioural and emotional problems, needs to be addressed and its aetiology explored through genetic, longitudinal and causal modelling: do the strong associations of aggression and other emotional and behavioural problems reflect a shared genetic vulnerability for multiple disorders, or do some disorders causally lead to other problems?

The absence of rater differences in co-occurrence patterns does not imply that rater’s views are interchangeable. Previous research suggested that, in general, mothers observe more behaviour problems in their children than fathers do [30]. We also see this pattern in the current paper, and consistently observe it across all counties. The differences in assessment between fathers and mothers in the levels of behavioural problems they observe may indicate that they both introduce their rater-specific view on the behaviour of the child [31], or that fathers and mothers interact with their offspring in different contexts.

The similarities across raters and countries indicate that large-scale gene-finding efforts of aggressive behaviour and its co-occurring psychopathologies across multiple cohorts will be feasible/successful. Such an effort is currently in progress within the ACTION consortium in collaboration with other cohorts and consortia that have collected measures of aggression in children as well as DNA samples for genotyping [32]. The results of this international genome-wide association meta-analysis (GWAMA) are expected to yield insight into the genetic variants that influence aggression across childhood and offer possibilities for the construction of polygenic scores which may be used in prediction models [33, 34] and gene-environment modelling [35]. Besides a GWAMA approach, which includes samples from multiple age groups, genome-wide epigenetic profiling will be done to compare methylation in several statistically well-powered contrasts (such as genetically identical twin pairs discordant for aggression) in children. Monozygotic (MZ) twins pairs who are longitudinally discordant of aggression, also offer a unique possibility to gain an understanding of the environmental risk factors associated with complex behaviour such as aggression [36].

Genetic and epigenetic effects do not act in isolation, so the results of these studies will need to be investigated in (epi)gene x environmental interplay models to understand the differences between children in aggression. Twin data may offer a first insight into the importance of gene-environment dependencies. Analyses of behavioural problems in 5-year-old twins showed strong evidence for larger environmental influences in children who were genetically more at risk for problem behaviour [37]. The available large-scale phenotypic, environmental and genotypic databases in ACTION will allow the development and application of these methods for gene-environment interaction and correlation.

Although it is known that co-occurrence is a risk factor for persisting symptoms (e.g. [38]), the implications for treatment are under-investigated. The current paper underlines that co-occurrence of behavioural and emotional problems with childhood aggression is highly prevalent. Instead of excluding children with multiple problems, specific trials should be undertaken to investigate the effectiveness of treatment and improve treatment for this group that requires our utmost attention. Of course, the question then arises what would be more effective, e.g. treatment targeting all psychopathologies at the same time or treatment at symptom level. It is essential to gain knowledge about the etiological and sequential effects of the comorbid disorders. If one disorder also is found to precede another disorder, treatment can be adjusted and specified. To be able to initiate such treatment specificity, we need to conduct cross-lag longitudinal analyses to examine whether aggression is driving the other psychopathologies, or if aggression is a result or outing of other problems. If one set of symptoms drives the rest, then intervention should focus on early detection and prevention.

We conclude that childhood aggression co-occurs with nearly all other behavioural, emotional and social problems, from both externalising and internalising domains, regardless of rater, gender, measurement instrument or country. These findings indicate that aggression during childhood and adolescence rarely occur in isolation, and that other behavioural and emotional problems are common in children with aggressive problems.

Future progress

The finding that aggression co-occurs with nearly all other behavioural, emotional and social problems during childhood puts aggression in the centre of scientific attention. If and when causes of differences in aggression during childhood are better understood, this information may aid in the development of prevention and intervention strategies. To this end, we designed the EU-ACTION project (see Fig. 1). The main objective of ACTION is to improve the understanding of the causes of individual differences in aggression among children to better inform the development of prevention and treatment strategies.

Fig. 1
figure 1

Work plan strategy of ACTION

ACTION has described current clinical practices in Europe with respect to childhood aggression and identified drawbacks in prevention and intervention of clinical aggression (also known as paediatric conduct disorders). An online semi-structured questionnaire investigating the status of national guidelines (N = 29 academic experts; 23 countries) and an online semi-structured questionnaire exploring clinical practices (N = 94 clinicians; 22 countries) on diagnosing and treating children with severe behavioural problems across Europe were developed. Several countries have official clinical guidelines, while others have at least some unofficial documents. In general, primary and secondary preventions were absent or poorly developed, whereas specific interventions for severe behavioural problems were very diverse across Europe. Improving parent–child interactions, parent/teacher interventions and collaborative approaches were most frequently identified as successful treatment elements. Several needs were listed by experts and clinicians, which will fuel further research within ACTION and beyond. The current findings on co-occurrence of aggression indicates that information on these current drawbacks could also be informative for other psychopathologies.

A challenge in combining large cohort studies carried out in different countries is the assessment of aggression. Within different countries and cohorts, different instruments are used. In a subsample of the Netherlands Twin Register, we have invited a group of parents of 9-year-old twins to complete multiple assessment instruments to have a ‘reference set’ or ‘backbone’ for phenotype imputation.

The first results with respect to genetic and environmental contributions to the variation and longitudinal stability in childhood aggressive behaviour [39] indicated high stability and heritability of aggressive behavioural problems. Heritability was on average around 60–80% without any large gender differences in the magnitude of genetic effects. In boys, shared environment explained around 20% of the variation in aggression across all ages, while in girls its influence was absent around age 7 and only came into play at later ages. Longitudinal genetic correlations explained most of the stability of aggressive behaviour. These results are encouraging for gene-finding studies. In earlier work, the first molecular genetic evidence for aggression in childhood was reported [40]. Using genomic relationship matrix restricted maximum likelihood (GREML) analyses significant influences of common SNPs were estimated for externalising problems (SNP h2 = 0.44), for attention problems (SNP h2 = 0.37–0.71) and total problems (SNP h2 = 0.18). A previous attempt to discover genomic locations of interest for childhood and adolescent aggression (N = 18,988) identified one region in chromosome 2 (2p12) at near genome-wide significance (top SNP rs11126630, P = 5.30 × 10− 8). The gene-based analysis indicated association with variation within AVPR1A with aggressive behaviour. It was concluded that common variants at 2p12 show suggestive evidence for association with childhood aggression [41]. To replicate this finding and to initiate new findings we will use newly developed multivariate genome-wide meta-analysis methodology, in which the power of sample overlap (e.g. due to having a paternal and maternal rating of the same child at the same age) is leveraged instead of omitted [42]. In line with the results, we include ADHD and ADHD-related problems, as well aggressive behaviour in this collaborative project. With this approach, we will be able to identify not only genomic regions of interest for aggression or ADHD, but also genomic regions that play a role in the co-occurrence of these psychopathologies.

In addition to existing genotype datasets, new DNA samples are collected for epigenetic research in clinical cases (children that are referred to child psychiatric clinics in the Netherlands) and in identical twins concordant and discordant for aggression. While DNA collection and epigenetic profiling in these children is in progress, we gained initial insight into the association between aggression and DNA methylation patterns by analysing available data on aggression available for adults [43]. DNA methylation was measured in whole blood by the Illumina HM450k array in more than 2000 adults for whom Adult Self-Report [44] data on aggression were available. No genome-wide significant methylation hits were identified, but gene-ontology (GO) analysis, in which categories of genes rather than single methylation sites were tested, highlighted that genes involved in developmental and central nervous system processes are enriched among the higher-ranking genes from the epigenome-wide meta-analysis (EWAS). This study is now followed by a meta-analysis EWAS (EWAMA) in children and adults across multiple cohorts. This EWAMA includes multiple cohorts with a sample size of over 10.000.

In addition to genetic and environmental factors acting additively to the development of childhood aggression, genes and environment may interact. Such interactions can be thought of as genes controlling sensitivity to the environment, or as the environment controlling the expression of genes. Genes and environment may also be correlated when genes alter the exposure to relevant environmental risk factors. We know that for traits such as aggression children are not randomly distributed over environments and describing environmental effects as “causal” may lead to wrong conclusions/interventions. Several mechanisms can be at play to explain the non-random distribution of genotypes over environments [45]: children who inherit genes that make them susceptible to exhibiting aggression are likely to grow up in aggressive homes (passive rGE), their genotypes may trigger aggression in others (reactive rGE) and they may seek out aggressive peer groups (active rGE). The analyses of rGE thus are closely related to issues of gene-environment independence and to questions of causality. The analyses of GxE interaction will employ several approaches that can make use of the large existing datasets. The first approach focuses on the estimation of the total contribution of genes when environmental exposures have been measured. In this approach genotypes and other, non-measured, influences are modelled as latent factors. Because of the presence of genome-wide marker data, a second approach is to estimate GxE interaction in a design with measured genotypes and environmental exposures (note that because of the twin design the remaining variance can still be attributed to latent G and E). The causal relation of environmental exposure and later outcome may be complex, but longitudinal twin data offer excellent opportunities to test models of causality versus other models of association between genes and environment [46, 47].

A final piece of the puzzle is sought in the assessment if biomarker and metabolomics profiles in clinical cases and MZ twins discordant for aggression. Given the fact that human aggressive behaviour is heterogeneous and that most effective therapeutic agents only work on the serotonergic system a comprehensive study of the role of the amino acid neurotransmitters (including both their precursors and degradation metabolites) and peptide-based neurotransmitters is warranted. In addition to this biomarker approach, ACTION will include a metabolomics approach and the platforms we are measuring include amines, organic acids and steroids.

Results from ACTION will be integrated into an empirical-based framework of aggression. The sample sizes of ACTION will allow us to examine the interplay between risk factors and test hypotheses to identify modifiable risk factors for childhood aggression. Thereby, our findings may inform prevention and treatment strategies, and assist in individual risk profiles based on combination of modifiable and non-modifiable risk indicators. Translation of results will be supported by several internet applications and disseminating the results via the ACTION website (http://www.action-euproject.eu/).