Introduction

The importance of nature (genetic vulnerability) and nurture (environment) in the development of behavioral traits has been widely accepted over the last 50 years [28]. In recent years, the interest of researchers in developmental psychopathology has shifted from assessing the relative contribution of genetic and environmental influences to investigating the role of gene–environment correlation and gene–environment interaction in the causation of developmental psychopathology.

The results of twin and adoption studies suggest that genetic factors contribute to the development of externalizing behaviors [13, 14, 31, 42], with heritability estimates of between 0.51 and 0.80 for broad externalizing behaviors, and of about 0.75 for attention problems and overactivity. However, many risk factors, including genetic factors, are not disorder-specific [14, 16]. High levels of comorbidity among ADHD, ODD and CD are found [8], which is likely to be due to a substantial degree of shared genetic liability, either operating directly, or indirectly through gene–environment correlations or interactions [24]. After controlling for the overlap between internalizing and externalizing symptoms, familial risk to externalizing behaviors (FR-EXT) is specifically associated with externalizing but not with internalizing psychopathology in the offspring [26]. In line with earlier studies [9, 25, 26, 43], we used a proxy for familial risk, which was based on data concerning life time parental externalizing psychopathology. Since the heritability of externalizing disorders is relatively high, and the etiologic contribution of common environmental risk factors to externalizing disorders is relatively modest, we may assume that mostly genetic factor drive the FR-EXT measure. Note that this familial risk might be a consequence of both genetic and environmental factors [39].

However, environmental influences (e.g. parenting styles) are also important to the causation of externalizing behaviors. A myriad of parenting styles and dimensions has been extensively studied since the 1930s [41]. As documented in several influential reviews [18], warm and accepting parenting styles are consistent predictors of favorable developmental outcomes in children, whereas hostility and rejection predict unfavorable outcomes. Adolescents who perceived a lack of parental warmth and high levels of rejection and overprotection exhibited more broad-band externalizing behaviors, aggressive behaviors, and delinquent behaviors [35]. Furthermore, perceived parental rejection was found to be the strongest predictor of the level of hostility in adolescents and adults [20], and of aggression/delinquency and attention problems [4].

Socialization research has mainly focused on finding connections between variation in child rearing and behavioral outcome in genetically related parent–child dyads. However, these associations between child-rearing styles and behavioral outcomes may be “spurious” (noncausal) [34] because genetic effects may cloud the interpretation of findings. First, the effect of genetic risk can be indirect, indicating that the effect of poor parenting may actually be the effect of susceptibility genes, or vice versa [27, 36]. The results of child-centered studies suggest that (retrospectively) perceived rearing styles are genetically influenced [29, 32, 33]. More specific, perceptions related to warmth in the family are genetically affected to a greater extent than perceptions of control. Second, children may differ—as a function of their personality and/or FR—in their susceptibility to rearing influences [6]; indicating heightened susceptibility to the negative effects of risky environments and to the beneficial effects of supportive environments [5, 6]. Therefore, research on gene–environment interaction effects of parenting and child characteristics may provide more insight into the complexity of child development.

In this study, we aimed to replicate and extend earlier research regarding the effects of FR-EXT and parenting styles on externalizing behaviors. We used the child’s perspective of parenting practices, because the child’s beliefs about parental behavior seem to have more influence on its social adjustment than the parents’ perspective of their own parenting behavior [12]. In order to avoid shared rater effects, we used parent’s and teacher’s evaluation of externalizing behaviors. By using measures for FR-EXT and perceived parenting styles in the same model, we could control for genetic influences on parenting and determine whether FR-EXT affects the way children respond to given environmental conditions [30]. Furthermore, many studies have focused on gene–environment interaction in antisocial behavior [21] and less in other domains of externalizing psychopathology. We attempted to answer the following questions:

  1. 1.

    does FR-EXT increase the risk of externalizing behaviors,

  2. 2.

    do perceived parenting styles [(lack of) emotional warmth, rejection, and overprotection] increase the risk of externalizing behaviors, and is a specific parenting style associated with specific externalizing behaviors, and

  3. 3.

    are children who have a FR-EXT differential susceptible to specific parenting styles with regard to externalizing behaviors?

Methods

Sample

The subjects were participants of the TRacking Adolescents’ Individual Lives Survey (TRAILS), a prospective cohort study of Dutch preadolescents who will be measured biennially until they are at least 25 years old. The key objective of TRAILS is to chart and explain the development of mental health from preadolescence into adulthood, in terms of underlying vulnerability and environmental risk. Participants were 10 to 12 years old and lived in the three largest cities and some rural areas in the north of the Netherlands. A detailed description of the sampling procedure and methods is provided elsewhere [11].

Briefly, the present study involves the first assessment wave of TRAILS, which ran from March 2001 to July 2002 [9, 11, 25, 43]. Of all children approached for enrollment in the study (i.e., children selected by the municipalities and attending a school that was willing to participate; N = 3,145 children from 122 schools, with 90.4% of the schools responding), 6.7% were excluded because of incapability or language problems. Of the remaining 2,935 children, 76.0% were enrolled in the study, yielding a sample size of 2,230. Both the child and the parent consented to participate. The mean age of the children was 11.09 years (SD = 0.55), 50.8% were girls, 10.3% were children who had at least one parent born in a non-Western country, and 32.6% of children had parents with a low educational level (i.e., a lower track of secondary education was the highest level attained). Responders and non-responders did not differ with respect to the prevalence rates of psychopathology and associations between sociodemographic variables and mental health outcomes [11].

Data collection

Well-trained interviewers visited one of the parents (preferably the mother, 95.6%) at home to administer an interview covering a wide range of topics, including the child’s developmental history and somatic health, parental psychopathology, and care utilization. Besides the interview, the parent was also asked to fill out a written questionnaire. Children were evaluated at school, where they filled out questionnaires in groups, under the supervision of TRAILS assistants, and were assessed individually. Teachers were asked to fill out a brief questionnaire for each TRAILS child in their class. Measures that were used in the present study are described below.

Familial risk to externalizing behaviors (FR-EXT)

Five dimensions of lifetime parental psychopathology were assessed (depressive disorders, anxiety disorders, substance dependence, antisocial behavior, and psychosis), using the TRAILS Family History Interview (FHI), which was administered at the parent interview [26]. Each dimension was introduced with a vignette describing the main DSM-IV [3] characteristics of the psychopathology, followed by a series of questions to assess lifetime occurrence, professional treatment, and medication use. Both biological parents were assessed during the interview, using a single informant, typically the mother. For each spectrum, the parents were assigned to one of the categories 0 = (probably) never had an episode, 1 = (probably) yes, or 2 = (probably) yes and treatment and/or medication were provided. For antisocial behavior, the last category was 2 = (probably) yes and picked up by the police. Prevalence rates in mother and fathers, respectively, were, for depression 27 and 15%, for anxiety 16 and 6%, for substance dependence 3 and 7%, and for antisocial behavior 3 and 7%. The FHI rates were by and large comparable to the CIDI-DSM-IV lifetime rates obtained by direct interviewing in NEMESIS [7]; the exception being fathers’ rates for anxiety disorder and substance dependence that were 40% too low [9, 26, 43]. We did not focus on parental psychosis, depressive and anxiety disorders.

The construction of FR-EXT was based on the reported path coefficients regarding substance abuse and antisocial behavior by Kendler et al. [16], who performed multivariate twin modeling to investigate the structure of genetic risk for common psychiatric and substance use disorders. Since twin studies provide compelling evidence that the familial transmission of alcohol and drug dependence and adult antisocial behavior is attributable to a highly heritable general vulnerability that contributes to a spectrum of externalizing behaviors (h 2 = 0.80) [10, 14, 16], we combined the coefficients reported by Kendler et al. [16] for alcohol dependence and drug abuse/dependence to create the variable substance abuse/dependence (SAD), and likewise we created the variable antisocial behavior (ASB) by combining the coefficients for antisocial behavior and conduct disorder. Subsequently, FR-EXT scores were computed by filling in the following regression equation: FR-EXT for externalizing behaviors = SAD mother + SAD father + ASB mother + ASB father, FR-EXT ranged from 0 to 8 (skewness 3.78, kurtosis, 17.63). The empirical justification for the construction of FR-EXT is provided elsewhere [26]. To enhance comparison with other TRAILS reports [9, 25, 43], missing values (N = 67) were replaced by the sample mean (0.18). Two groups were created on the basis of the distribution of FR-EXT (total N = 2,230); children with no FR-EXT (82.2%) and children with FR-EXT (17.8%). In the regression analyses, these two groups were used as a dummy variable.

Perceived parenting styles

The Egna Minnen Beträffande Uppfostran (My Memories of Upbringing) for Children [EMBU-C] [19] was developed to assess children’s and adolescents’ perception of parents’ rearing practices. The original EMBU-C contained 81 items. Markus et al. [19] developed a shorter version, which we used but without the favoring subject factor. We omitted this scale prior to administration, due to an internal consistency below 0.60. Each item, scored on a 4-point scale (1 = no, never, 2 = yes, sometimes, 3 = yes, often, 4 = yes, almost always), was presented for both the father and the mother. The EMBU-C contains the factors emotional warmth, rejection, and overprotection. The main concepts of emotional warmth are giving special attention, praising approved behavior, unconditional love, and being supportive/affectionately demonstrative. An example of an item is: ‘Do your parents make it obvious that they love you?’. Rejection is characterized by hostility, punishment (physical or not, abusive or not), and derogation and blaming of the child. An example of an item is: ‘Do your parents sometimes punish you even though you haven’t done anything wrong?’ Overprotection is characterized by fearfulness and anxiety for the child’s safety, guilt engendering, and intrusiveness. An example of an item is: ‘Are your parents concerned about what you do after school hours?’

Principal component analysis with three factors (emotional warmth, rejection and overprotection) as criterion, followed by VARIMAX rotation, mainly confirmed the results by Markus et al. [19]. All items loaded on the designated scale, with the exception of five items of the rejection scale (the items 8, 24, 35, 71, and 76 in the article of Markus et al. [19], which were rejected because they had loadings lower than 0.30 or had a loading that differed by less than 0.10 with the second highest loading. In the study by Markus et al. [19], the items had relatively low loadings, on average 0.36. The three factors explained 34.0 and 32.5% of the variance in the ratings on fathers and mothers.

The scale for emotional warmth contained 18 items with α = 0.91 for both fathers and mothers; the rejection scale contained 12 items with α = 0.84 for fathers and α = 0.83 for mothers, and the overprotection scale contained 12 items with α = 0.70 for fathers and α = 0.71 for mothers. The answers for both fathers and mothers were highly correlated (r = 0.79 for emotional warmth, 0.67 for rejection, and 0.81 for overprotection), so we felt it was justified to combine them. The test–retest stability of a shortened version of the EMBU-C (10-item scales) over a 2-month period is satisfactory, rs = 0.78 or higher [22]. Earlier studies have reported on the validity and cross-cultural equivalence of the EMBU-C [19].

Externalizing behaviors

Externalizing behaviors were assessed with the Child Behavior Checklist (CBCL), one of the most commonly used questionnaires in current child and adolescent psychiatric research [1, 44]. It contains a list of 112 behavioral and emotional problems which parents can rate as 0 = not true, 1 = somewhat or sometimes true, or 2 = very often true in the past 6 months. In addition to the CBCL, we administered the Teacher’s Checklist of Psychopathology (TCP). The TCP contains descriptions of problem behaviors corresponding to the syndromes scored with Achenbach’s Teacher Report [11]. Response options range from 0 (not applicable) to 4 (very clearly or frequently applicable). In this study, we focused on the CBCL syndromes attention problems (α = 0.81), aggressive behavior (α = 0.89), and delinquent behavior (α = 0.68). Consistent with other reports [2], the agreement between parent-reported and teacher-reported problems was only moderate (r = 0.47 for inattention, r = 0.37 for impulsivity/hyperactivity, r = 0.32 for aggression and r = 0.27 for delinquency). We feel that the two informants perceive different aspects of problem behavior and that differences between informants are meaningful.

Statistical analyses

To obtain comparable regression coefficients, z scores were used for all dependent variables and parenting style variables. All analyses were performed using the Statistical Package for the Social Sciences (SPSS for Windows, version 14.0).

Sex differences were examined by means of t tests; associations between variables by means of Pearson correlations. Multiple linear regression analyses were used to determine the relative contribution of FR-EXT, emotional warmth, rejection, overprotection and gender to parents’ and teachers’ ratings of inattention, hyperactivity/impulsivity, aggression, and delinquency. Interaction terms between FR-EXT and perceived parenting styles were also entered into the regression models.

Multicollinearity was not present in our data (greatest VIF value is 1.81 for the predictor parental rejection).

All analyses were repeated without implementing FR-EXT data. The results of the t tests, Pearson correlations, multiple regression analyses were largely comparable. This indicates that the likelihood of a bias due to implementing data is small [15].

Results

Descriptives

Table 1 contains means and standard deviation of FR-EXT, parent- and teacher-rated externalizing behaviors and perceived parenting styles, separately for boys and girls. Except for FR-EXT, all variables showed significant sex differences. Boys scored higher on all behavioral domains, indicating more problem behavior. Girls perceived more emotional warmth, and less rejection and less overprotection than boys.

Table 1 Sex differences in familial risk to externalizing behaviors, perceived parenting styles and externalizing behavior

Bivariate correlations between predictors and externalizing behaviors

Pearson correlations between the predictors (FR-EXT, emotional warmth, rejection, overprotection, and gender) and the dependent measures are summarized in Table 2. FR-EXT, rejection and overprotection were positively correlated with all dependent variables. Thus children with higher levels of FR-EXT showed more externalizing behaviors, as did children with higher levels of parental rejection and overprotection. Emotional warmth was negatively correlated with all dependent variables, indicating that children who perceived less parental emotional warmth showed more problem behavior, or conversely, children who perceived more parental warmth showed less problem behavior.

Table 2 Bivariate correlations between familial risk to externalizing behaviors, perceived parenting styles and externalizing behavior

FR-EXT was not related to emotional warmth or rejection, but was very weakly related to overprotection (r = 0.04), which indicates that there was no gene–environment correlation. Furthermore, overprotection was positively related with rejection (r = 0.36, p < 0.001) and emotional warmth (r = 0.18, p < 0.001), and rejection was negatively related with emotional warmth (r = −0.36, p < 0.001).

Multivariate models of inattention, hyperactivity, and impulsivity

Multiple linear regression analysis showed that FR-EXT, perceived parental styles and gender together explained 11% of the variance in parent-reported inattention (Table 3). There were main effects of FR-EXT (vs. no FR-EXT) (Fig. 1), emotional warmth, rejection, and overprotection, and gender, with boys having higher scores than girls. No significant interaction effects were found.

Table 3 Multiple regression analyses, for each separate (standardized) dependent variable and informant
Fig. 1
figure 1

Main effect of familial risk to externalizing behaviors (FR-EXT) as predictor of parent-reported inattention. Note The Y-axis represent the predicted means of the standardized CBCL score

The model explained 9–12% of the variance in teacher-rated inattention and hyperactivity/impulsivity. Main effects were found for FR-EXT, emotional warmth, overprotection, and gender. Two significant interactions were found for teacher-rated hyperactivity/impulsivity. A positive interaction effect was found between FR-EXT and rejection, indicating that children with FR-EXT were more vulnerable for parental rejection in exhibiting more hyperactive and impulsive behavior than children with no FR-EXT (Fig. 2). In contrast, a negative interaction effect was found between FR-EXT and overprotection. This means that the effect of parental overprotection was stronger for children with no FR-EXT, than for children with FR-EXT (Fig. 3).

Fig. 2
figure 2

Interaction effect between familial risk to externalizing behaviors (FR-EXT) and parental rejection as predictor of teacher-reported hyperactivity and impulsivity. Note The Y-axis represent the predicted means of the standardized TCP score, the X-axis the standardized score on the specific subscale of the EMBU

Fig. 3
figure 3

Interaction effect between familial risk to externalizing behaviors (FR-EXT) and parental overprotection as predictor for teacher-reported hyperactivity and impulsivity. Note The Y-axis represent the predicted means of the standardized TCP score, the X-axis the standardized score on the specific subscale of the EMBU

Multivariate models of aggressive behavior

The models explained 11 and 10% of the variance in aggressive behavior reported by parents and teachers, respectively (Table 2). Main predictors of parent-rated aggression were FR-EXT, emotional warmth, rejection, overprotection, and gender. Main predictors of teacher-rated aggression were FR-EXT, emotional warmth, overprotection, and gender.

Multivariate models of delinquent behavior

Parent-rated delinquent behavior could be explained by FR-EXT, emotional warmth, rejection, and gender (Table 2). These factors together explained 12% of the variance. The results for teacher-rated delinquent behaviors were comparable: FR-EXT, low levels of emotional warmth, high levels of overprotection, and male gender were main risk factors. One interaction effect was found between FR-EXT and rejection: children with FR-EXT were more vulnerable for parental rejection in showing more delinquent behavior than children with no FR-EXT.

Discussion

The study was set up on the premise that familial risk to externalizing behaviors (FR-EXT) and perceived parenting styles are main predictors of externalizing behaviors in preadolescents, and interact with each other. In line with earlier studies [9, 26, 43], we found FR-EXT to be a risk factor for externalizing behaviors in preadolescents, as evidenced by significant bivariate correlations between FR-EXT and all seven dependent measures of externalizing behaviors. In addition, in multivariate analyses, we found that FR-EXT was a main risk factor for the whole range of externalizing behaviors.

Beside the main effect of FR-EXT, three perceived parenting styles (emotional warmth, rejection, and overprotection) were found to be main predictors of externalizing behaviors, even after adjustment for other predictors in our multivariate models. In general, preadolescents who perceived a lack of parental emotional warmth, and high levels of parental rejection and overprotection were described as being more inattentive, aggressive, and delinquent by both parents and teachers. These findings are in line with other research [35]. In addition, perceived parental rejection had a greater effect on aggression and delinquent behavior than on inattention, impulsivity, and hyperactivity [23]. However, the effect of perceived parental rejection was present for parent-rated and not teacher-rated behavior. This finding may be explained by the fact parents and children are 50% genetically related. Therefore, the genetic influences of parent’s and children’s behavior may interact with each other and evoke negative parenting, which in turn, increases the likelihood of behavioral problems.

Our study design enabled us to determine the presence of gene–environment correlation and interaction on externalizing behaviors. We did not find correlations between FR-EXT and parental emotional warmth or rejection, and only a weak correlation between FR-EXT and parental overprotection. This indicates that FR-EXT and parenting styles seem to be two independent dimensions. This finding is in contrast with the idea that people modify and select their own environments, which results in (re)active gene–environment correlations [37, 38]. With regard to the moderating effects, we entered three interaction terms in seven multivariate analyses, and found three to be significant. Although this number of significant interactions is at change level (n = 21, test proportion 0.05, p < 0.10 by binomial test), these interactions seem to be meaningful. The interactions between FR-EXT and perceived parental rejection was significant for both teacher-rated hyperactivity and teacher-rated delinquency. This indicates that children with FR-EXT are more susceptible for parental rejection in showing more hyperactive/impulsive and delinquent behavior than children with no FR-EXT. Another interaction between FR-EXT and overprotection was significant for teacher-rated hyperactivity/impulsivity, indicating that children with no FR-EXT were more susceptible for parental overprotection. These findings are similar to results reported in studies of temperament–environment interactions, and problem behavior [25, 43]: perceived parental rejection and overprotection interacted both with temperament in predicting depressive and delinquent behaviors in preadolescence, respectively.

Our unexpected finding, that virtually no gene–environment correlation and only a few gene–environment interactions were present, may be explained in several ways. First, FR-EXT was defined at a phenotypic level and was based on family history as a proxy for genetic vulnerability. With the caveat in mind that FR-EXT may reflect both genetic and environmental influences, the availability of DNA analysis in the next future will allow us to refine our analyses by including genetic polymorphisms as risk factors. Second, the use of family history interviews, when compared with direct interviews of relatives, may have led to underreporting of lifetime parental psychopathology, and thus underestimation of associations is possible [26]. However, except for father’s rate for substance dependence, our prevalence rates were comparable to life time rates obtained by direct interviewing [7]. Third, we used a community sample, covering a wide range of environments and parenting styles; however, specific environments or parenting styles may have a more negative impact than other environments or styles [34, 40]. Therefore, future studies should include clinical samples and more extreme parenting styles. Fourth, we evaluated parenting style in terms of the child’s perception, which may have been influenced by the child’s own FR-EXT. Though the level of FR-EXT is reported by the parent, the genetic make-up of the child may be related with their own mental representation regarding parenting styles of their parents. In that way the genetic make-up of the child is discounted in the perception of the child, which may have resulted in and underestimation of gene–environment correlations. Lastly, we looked at the effect of perceived parenting in preadolescents but greater effects may be found in younger children [6, 17]. Moreover, it should be borne in mind that attention problems, aggression, and delinquent behavior are heterogeneous problem behaviors that may have different biological and etiological backgrounds.

Conclusion

FR-EXT and perceived parenting styles independently affected the expression of parent- and teacher-rated externalizing behaviors and can be seen as two relatively independent areas of causation in a population based sample of preadolescents. However, some environmental moderation of FR-EXT is present, albeit of a small effect. Further research is necessary to study whether our findings hold for more extreme forms of parenting styles and for longitudinal analyses of externalizing behaviors throughout adolescence.