Heritable influences on behavioural problems from early childhood to mid-adolescence: evidence for genetic stability and innovation

Background Although behavioural problems (e.g. anxiety, conduct, hyperactivity, peer problems) are known to be heritable both in early childhood and in adolescence, limited work has examined prediction across these ages, and none using a genetically informative sample. Method We examined, first, whether parental ratings of behavioural problems (indexed by the Strengths and Difficulties questionnaire) at ages 4, 7, 9, 12, and 16 years were stable across these ages. Second, we examined the extent to which stability reflected genetic or environmental effects through multivariate quantitative genetic analysis on data from a large (n > 3000) population (UK) sample of monozygotic and dizygotic twins. Results Behavioural problems in early childhood (age 4 years) showed significant associations with the corresponding behavioural problem at all subsequent ages. Moreover, stable genetic influences were observed across ages, indicating that biological bases underlying behavioural problems in adolescence are underpinned by genetic influences expressed as early as age 4 years. However, genetic and environmental innovations were also observed at each age. Conclusion These observations indicate that genetic factors are important for understanding stable individual differences in behavioural problems across childhood and adolescence, although novel genetic influences also facilitate change in such behaviours.


Introduction
Twin and family studies have indicated that behavioural problems in childhood (Saudino et al. 2005) and adolescence (Scourfield et al. 2004) are underpinned by substantial genetic influences. However, less is known about the genetic stability of these influences. It is possible that genetic influences are broadly stable across these ages, reflective of an early maturing biological basis to behavioural problems that endures throughout childhood and adolescence. By contrast, the heritable effects at each of these ages may reflect genetic innovation, with novel heritable effects becoming apparent over the course of development, perhaps reflecting prolonged neural maturation (Gogtay et al. 2004;Blakemore & Choudhury, 2006). An intermediate position, positing that stable genetic effects are evident over childhood and adolescence accompanied by novel genetic effects emerging over time, is also plausible. Determining which of these perspectives best characterizes the genetic architecture of behavioural problems over childhood and adolescence will be important in order to better understand the manner in which gene effects exert their influences on behavioural problems. In addition, if genes are to be discovered for child and adolescent behavioural problems it is critical to know whether the same genetic influences are apparent across these stages in life. To address these issues, the current study used a longitudinal twin design in order to investigate the extent to which genetic influences contribute to the stability and innovation of genetic and environmental effects on a set of core behavioural problems.

Child and adolescent behavioural problems
Understanding the aetiology of adolescent behavioural problemssuch as problems with peer relationships, hyperactivity, conduct, and anxiety, and (low) prosocialityis of importance as such behaviours show links to broad-based life and psychosocial outcomes (Rutter, 1995). Research has suggested that around one in five adolescents satisfy diagnostic criteria for at least one psychiatric disorder (Merikangas et al. 2010), with specific behaviour problems such as conduct disorder, hyperactivity disorder, and anxiety/emotional problems all relatively common among adolescents (Ford et al. 2003;Merikangas et al. 2010). These behavioural problems are of clear and broad interest because of their links to mental health, school achievement, antisocial behaviour, and delinquency (Rutter, 1995;Aunola et al. 2000;Fergusson et al. 2005).
Several measurement instruments exist to screen for such behaviours (Achenbach, 1991a-c), includingas well as being the focus of the current studythe Strengths and Difficulties Questionnaire (SDQ; Goodman, 1997). The SDQ is a brief dimensional measure that covers the core domains of child and adolescent psychopathology (conduct problems, hyperactivity-inattention, emotional symptoms/ anxiety, and peer relationship problems) alongside personal strengths (prosociality).

Genetic origins of behavioural problems
Research examining the aetiology of such problem behaviours has indicated a role for genetics. For instance, children's and adolescents' scores (ages 5-17 years) on the conduct problems scale of the SDQ have been reported to be substantially heritable (35-77%) (Scourfield et al. 2004). Similarly high heritabilities were found for all of the SDQ scales at age 7 years (Saudino et al. 2005) and age 16 years (Lewis et al. 2014). Moderate to substantial heritabilities (50.40) have also been found for related measures in childhood and adolescence, such as antisocial behaviour and callous-unemotional traits (Viding et al. 2005) and internalizing and externalizing behaviours (Bartels et al. 2004). In sum, child and adolescent behaviour problems across multiple domains exhibit moderate-to-large heritable effects.

Temporal stability of genetic effects
Less clear at this stage is whether these heritable influences on behavioural problems represent stable and enduring effects or instead are more transient and time-specific. While a number of phenotypically oriented studies have noted that a range of temperaments and behaviours (both normal and disordered) in adolescence and adulthood can be predicted by measures in early childhood (e.g. Caspi & Silva, 1995;Caspi et al. 1996;Hampson & Goldberg, 2006;Slutske et al. 2012), thus illustrating the temporal stability of behavioural problems, less work has addressed this question using a genetically informative design. Of work addressing SDQ domains, evidence for genes influencing both stability and innovation has been reported for prosociality from ages 2-7 (Knafo & Plomin, 2006) and for anxiety from ages 7-9 (Trzaskowski et al. 2012). More broadly, a study of externalizing and internalizing behaviours across four time-points from ages 3-12 found that a transmission model fitted the data best, such that genetic influences were stable over time but that novel genetic influences were apparent at each measured age (Bartels et al. 2004). More recent work found broad support for this finding, noting genetic stability on externalizing (Wichers et al. 2013), antisociality (Tuvblad et al. 2011), and attention problems (Chang et al. 2013) from ages 8-19, although again novel genetic influences were also observed for each of these variables over this time period. Broadly consistent findings have been reported for common childhood fears (e.g. snakes, spiders; Kendler et al. 2008).
Some reported findings, however, have diverged from this pattern of genetic stability. For example, a study examining conduct and delinquency reported that the genetic influences evident in early childhood were dissociable from those acting in mid-adolescence (Van Hulle et al. 2009). This observation (as the authors note) may have emerged as an artifact of the study designparental reports were used in early childhood, self-report was used in adolescence; but as this study is unique in its examination the time-period from early childhood through to adolescence the findings raise important questions concerning the nature of genetic stability in behavioural problems across development.

The current study
While research in recent years has begun to establish that heritable influences on different aspects of behavioural problems are fairly stable across broad stretches of developmentmost notably from age 10 onwardsspecific questions remain unanswered. First, as noted above, the sole study to address stability from early childhood into adolescence reported that antisocial behaviour was underpinned by genetic influences at ages 4-9 which were wholly distinct to those genetic effects acting on antisocial behaviour at ages 14-17. Moreover, it is unclear if a similar pattern of stability and change exists for broader components of behavioural problems, such as those indexed by the SDQ. More generally, while important longitudinal genetically informative research has been focused on behavioural problems, no longitudinal genetically informative work has directly addressed prosociality and peer problems across from early childhood into adolescence, and work addressing anxiety and hyperactivity, while of value in its own right, has largely focused at the aggregated level of internalizing and externalizing (e.g. Bartels et al. 2004;Wichers et al. 2013).
A better understanding of the nature of genetic influences on behavioural problems over this important period of development is important for several reasons. First, knowledge of the dynamics of genetic influences on behavioural problems will serve to more accurately characterize the origins of behavioural problems. Second, environmental influences, and specifically non-shared-environment effects, are known to be large for behavioural problems, but this variance component also includes measurement error and so it is unclear whether this aetiology reflects a stable influence on behavioural problems or is better understood as a more transient factor. Third, if distinct sets of genes underpin behavioural problems at different ages, locating the causal variants will be significantly more challenging if researchers aggregate samples across ages. As such, knowledge of the genetic stability of behavioural problems across these ages will inform gene discovery strategies on whether aggregation across ages is a valid approach.
To rectify these gaps of knowledge in the literature, here we examined the extent to which genetic influences on childhood and adolescent behavioural problems are underpinned by stable genetic effects present at each measured age, or are more accurately described by genetic effects specific to each developmental age. We used the SDQ instrument to tap behavioural problems with parental scores taken at each age in order to maintain a consistent mode of rating across each wave of measurement. Ratings were available for ages 4, 7, 9, and 12 for all SDQ scales, and also for age 16 for conduct, hyperactivity, and prosociality.

Participants
Participants were drawn from the Twins Early Development Study (TEDS), which is an ongoing longitudinal study following monozygotic (MZ) and dizygotic (DZ) twins born in England and Wales between 1994 and 1996 (Oliver & Plomin, 2007). The TEDS sample is representative of the UK population (Kovas et al. 2007) and the project received approval from the Institute of Psychiatry Ethics Committee. Twin zygosity was determined using a parental rating measure of similarity and DNA genotyping (Price et al. 2000). The number of complete twin pairs at each age were as follows: MZ male pairs: n = 533-1200; MZ female pairs: n = 670-1368; DZ male pairs: n = 505-1198; DZ female pairs: n = 555-1255; and DZ opposite-sex pairs: n = 996-2355.

SDQ
The SDQ is a short but reliable instrument (25 items; Goodman, 2001;Stone et al. 2010) for measuring psychosocial problems in children (Goodman, 1997). The SDQ consists of five scales measuring anxiety, conduct problems, hyperactivity-inattention, peer problems, and prosocial behaviour. Scores for all subscales were acquired by parental rating when the child was 4, 7, 9, and 12 years old. Scores for prosociality, conduct, and hyperactivity were also acquired by parental rating when the individual was 16 years old. Cronbach's α was low (<0.60) for conduct (all ages) and for peer problems (ages 4 and 7), although in line with previously reported values (Goodman, 2001), but broadly acceptable for the rest of the SDQ measures.

Analysis
Correlations between twins differing in their degrees of genetic relatedness (e.g. MZ and DZ twins) are informative as a guiding heuristic to estimate relative magnitudes of genetic and environmental effects (Plomin et al. 2013). The presence of genetic effects is inferred if correlations between MZ twins are larger than correlations for DZ twins. The presence of shared-environment effects is inferred if correlations for DZ twins are larger than half the magnitude of the correlations for MZ twins. Finally, nonshared-environment effects are inferred if correlations for MZ twins are less than 1.0, and so this variance component also contains measurement error. These correlation analyses were extended using formal model-fitting of variance-covariance matrices for the twin data. This approach allows parameter estimates in univariate models to be formally tested for significance as well as allowing multivariate modelsthe core focus of the current studyto be analysed.
In the current study, longitudinal analyses were central to our tests. We sought to determine the extent to which genetic effects underlying SDQ measures across the ages reflected stable v. novel genetic influences. To perform this analysis, we compared four classes of model: the Cholesky decomposition, the common pathway model, the independent pathway model, and the simplex/transmission model, which are each detailed below (see Fig. 1 and online Supplementary  Figs S1-S3).
The Cholesky decomposition specifies as many factors as there are variables for each source of variance, with each subsequent factor having one fewer pathway than the preceding factor (see Fig. 1). In other words, for additive genetic effects (A) the first latent factor loads on all of the n measured variables. The subsequent latent factors load on n −1 , n −2 , . . ., n −i variables. In this way each factor accounts for as much of the remaining variance as possible, until the last factor accounts for just the remaining variance in the last measured variable. This is repeated for the sharedenvironment factors (C) and non-shared-environment factors (E). This design can also be viewed as allowing the genetic and environmental influences explaining SDQ measures at age 4 to also explain SDQ measures at the later ages, while also leaving specific genetic and environmental paths available to explain SDQ variation at the later ages that does not covary with SDQ measures at age 4. The common pathway model is a more restricted model, constraining all genetic and environmental variance though a single pathway (Kendler et al. 1987) (see online Supplementary Fig. S1). The independent pathway model also constrains genetic and environmental covariation, but instead of requiring all sources of this covariance to be channelled through a single pathway, this model allows for independent general effects of genes, shared environments, and non-shared environments (see online Supplementary  Fig. S2). Finally, the simplex/transmission model estimates genetic and environmental effects specific to each age, with these effects also allowed to 'carry-over' from one age to the subsequent age (see online Supplementary Fig. S3).

Results
Descriptive statistics for all measures are presented in Table 1. Assumption testing using all twins indicated that means and variances could be equated across twin order, zygosity, and sex for most variables, with the small number of significant differences observed consistent with the large number of tests performed.
Of note, however, was evidence for modestto-moderate mean sex differences, particularly for hyperactivity and prosociality (see Table 1 and text below for more detail). Sex-limitation modelling (testing for quantitative and qualitative genetic and environmental differences across sex) largely indicated that genetic and environmental influences could be equated across sex, with the significant differences that were observed being either small in magnitude or, again, consistent with the large number of tests performed. Following these observations we pooled our sample across sex, but used sex-residualized variables.

Phenotypic associations
Zero-order correlations across age for the SDQ subscales are shown in Table 2. In all cases significant associations were noted with correlations between ages 4-16 ranging from 0.25 to 0.31. These associations indicate a moderate degree of phenotypic stability for behavioural problems across childhood and into adolescence.

Multivariate twin analyses
We next assessed genetic and environmental links from SDQ measures across ages. We fitted each of our models individually for each SDQ scale using all available ages. In each case, the common pathway, independent pathway, and simplex models provided a significantly poorer fit to the data compared to the Cholesky (see online Supplementary Table S1) and so we retained the Cholesky as our final model for each of the SDQ measures (see Tables 3 and 4). Several key points are notable from these final models. First, heritable influences on all of the SDQ measures at age 4 are stable across all subsequent ages. Second, genetic innovation was also apparent, with significant novel heritable influences apparent at each age, and these novel heritable effects were in almost all cases stable over subsequent ages. While levels of innovation were broadly mirrored for non-shared-environment influences, the stability of non-shared-environment influences was considerably more modest for each measure. Finally, sharedenvironment effects were significant, albeit modest in magnitude, for anxiety, conduct, and prosociality, at age 4. These effects were broadly stable over subsequent ages. Some evidence for novel sharedenvironment effects was also observednamely, anxiety at age 7 and conduct at ages 7 and 9although these effects were typically modest.

Discussion
The current study provides several important insights into the aetiology of childhood and adolescent behavioural problems. At the phenotypic level we saw significant prediction across childhood and adolescence such that scores on the SDQ scales at age 16 (age 12 for anxiety and peer problems) were predicted by the corresponding SDQ score at age 4. These phenotypic overlaps were mostly accounted for by genetic factors, although evidence for modest shared environmental stability was observed (with the exception of hyperactivity). While research has suggested that genetic influences on antisociality shows distinct genetic underpinnings in early childhood v. early adolescence (Van Hulle et al. 2009), the current results using parental report across all measured ages (Van Hulle et al. relied on different raters in childhood v. adolescence) indicates that aspects of antisociality (i.e. conduct, hyperactivity), as well as broader elements of behavioural problems, are genetically stable from as early as age 4, in line with longitudinal work examining such phenotypes from later in childhood (e.g. Wichers et al. 2013). We also noted substantial genetic innovation at each age for each of the SDQ measures, which in turn were then seen to remain stable across the following ages. Environmental influences, and specifically non-shared-environment effects, largely acted to provide time-specific sources of variance. These findings converge with related work stressing that behavioural measures taken in early childhood predict later life psychopathology (Caspi et al. 1996;Slutske et al. 2012), although these effects tend to be modest-to-moderate when spanning substantial time ranges and the current results are consistent with such observations. The current results also converge with recent research demonstrating gene variants linked to risk for psychopathology in adults predict neuroanatomical variation in infants (Knickmeyer et al. 2014) and suggest that individual differences in measurable neurophysiology (such as regional greyor white-matter structure) in early childhood may also serve as a useful biomarker for subsequent behavioural problems. Our observation that behavioural problems are underpinned by novel heritable influences across development is also in line with recent work demonstrating heritable effects on cortical development operate in similar fashion (Schmitt et al. 2014). Assuming that cortical maturation mediates the impact of novel genetic influences on behavioural problems raises important questions regarding the brain bases of the early v. later emergent genetic bases of behavioural problems. Given that brain regions governing high-level cognitive control, such One randomly chosen member from each twin pair, n = 1886-5653.
All correlations are p < 0.001. as the dorsolateral prefrontal cortex, appear to mature comparatively late in adolescence (Gogtay et al. 2004;Blakemore & Choudhury, 2006), one possibility is that the novel genetic influences on behavioural problems, particularly those reflecting aspects of disinhibition or poor impulse control (e.g. hyperactivity), may reflect genetic influences on the latter-stages in development of these neural structures. Specific recommendations for future research are warranted. First, from a methodological standpoint, the degree of genetic innovation observed across the ages suggests that gene discovery should either focus on adult individuals, where genetic innovation is less evident (Johnson et al. 2005), or samples of children of equal age. Second, we relied entirely on parental report data. While this approach had the advantage of allowing us to use the same source of measurement across multiple ages, multiple raters would have potentially provided a more accurate assessment of the individual being rated. In particular, parental ratings may result in underestimates of heritability and overestimates of shared environment. While we cannot rule out this possibility, here we observed high heritabilities across all measures, with only modest shared-environment effects, suggesting that this possible source of bias had no major effect on our estimates. Third, the classical twin design is subject to a number of assumptions, such as the equal environment assumption (Neale & Cardon, 1992). Future studies that can capitalize on additional family structures in order to provide more assumption-free estimates would be valuable, although it is noteworthy that research testing whether violations of the equal environment assumption are apparent for psychopathology has found little evidence for this potential source of bias (Kendler et al. 1993;Derks et al. 2006). Finally, genetically informative latent growth modelling may help to further delineate the factors underlying developmental changes in behavioural problems and so would be a valuable direction to follow in future research.
In summary, here we report for the first time evidence that genetic influences on behavioural problems in early childhood are stable across childhood and into adolescence. Of importance, we also observed novel genetic influences at each age, with non-sharedenvironment effects also providing additional novel sources of variance across this important developmental phase. These findings support models of behavioural problems that posit early emerging and enduring individual differences in biology reflect such phenotypes, although also suggesting that the genetic influences underlying behavioural problems unfold across childhood and adolescence.

Supplementary material
For supplementary material accompanying this paper visit http://dx.doi.org/10.1017/S0033291715000173. Bold values indicate that the confidence intervals (95%) did not cross zero.