Peer relationships play a key role in academic and adjustment outcomes of school-aged children (Kiuru et al. 2015; for a review, see also Parker et al. 2006), such as self-esteem, social and emotional skills, and mental health (for a review, see Parker et al. 2006; Shin et al. 2016). In turn, lack of peer acceptance in childhood has been shown to be a risk factor for developing various forms of internalizing problem behaviors, such as anxiety and depression (e.g., Healy and Sanders 2018; Hoglund and Chisholm 2014; Menting et al. 2015; Sakyi et al. 2015). Internalizing problem behaviors may have long-term consequences leading to interpersonal difficulties, social competence deficits, and academic problems (Barrett et al. 2005; see also Tennant et al. 2017). The impacts of low peer acceptance are not, however, necessarily the same for all: the consequences of low peer acceptance have been suggested to be weaker or non-existent for children who have supportive parents (Cohen and Wills 1985; Deci and Ryan 2000). Supportive parenting has, indeed, been shown to foster positive adjustment (e.g., Healy et al. 2015; Pettit et al. 1997), and buffer the negative impacts of problematic peer relationships in childhood (e.g., Bowes et al. 2010; Buckholdt et al. 2016; Hazel et al. 2014; Healy and Sanders 2018) and adolescence (e.g., Hazel et al. 2014; McDonald et al. 2010; Sentse and Laird 2010). The present study focusing on the critical transition to school aimed to examine to what extent supportive parenting can provide a buffer against the detrimental impact of the lack of peer acceptance at the beginning of school on development of internalizing problems during subsequent school years.

Peer Acceptance and Internalizing Problem Behaviors

During the period from middle to late childhood (6–12 years of age), the importance of peers is increased owing to the growth in children’s interpersonal skills and the development of their peer relationships (for a review, see Parker et al. 2006). Peer relationships have been suggested to be an absolute necessity for healthy social and cognitive development, and socialization (Johnson 1980). According to Bukowski (2001), friendship provides a sense of self-value and personal validation, serves a protective function, facilitates learning and development of new skills, and shapes development through shared cultures (Maunder and Monks 2018). Peer status has been found to be relatively stable in elementary school years and to predict later social-emotional development and school adjustment (for a review, see Jiang and Cillessen 2005). Peer acceptance, which refers to being liked or accepted by the members of one’s peer group (Bukowski and Hoza 1989), is an important developmental milestone in the first years of school. Peer acceptance has been linked, for example, to children’s high psychological well-being, high academic performance, and successful adaptation in the school context (Hymel et al. 2011; Parker et al. 2006; Shin et al. 2016).

However, some children fail to build positive relationships with their peers. Peer relationship problems in childhood have been linked to a variety of developmental and psychological difficulties in childhood, adolescence, and adulthood (e.g., Arseneault 2017; Forbes et al. 2018; Sakyi et al. 2015; Shin et al. 2016). In particular, low peer acceptance has been recognized as an important antecedent of internalizing problem behaviors, i.e. symptoms of withdrawal, fearfulness, anxiety and depression (Healy and Sanders 2018; Ladd 2006; McDougall et al. 2001; Menting et al. 2015; Rubin et al. 2006; Sakyi et al. 2015). For example, a study on 6–12-year-old Australian children, followed over a 9-month period, found that supportive relationships with peers, in the form of peer acceptance, predicted changes in internalizing problems over time, with higher levels of peer acceptance predicting lower levels of later internalizing problems (Healy and Sanders 2018). A study on a sample with a mean age of 7-year-old (first grade) Korean children with 6- to 8-year follow-ups found that children’s peer relationships had longitudinal effects on their mental health and adjustment. Children who had qualitative peer relation problems (less or no friends) were more likely to exhibit internalizing problems in adolescents (Shin et al. 2016). Further, a study which followed French children over 18 years (from the ages of 4–16 years to 22–35 years), found that individuals who had no friends at 4–16 years of age, were approximately twice more likely to experience internalizing symptoms compared to those who had at least one childhood friend (Sakyi et al. 2015). Furthermore, a study which followed children annually from the age of 6 (kindergarten) through to the age of 8 (second grade) in the Netherlands also found that emotional problems, such as internalizing symptoms, share poor peer acceptance as an underlying factor in their development (Menting et al. 2015).

The transition to first grade has been described as a critical period for children’s development (e.g., Bossaert et al. 2011; Denham 2006; Ladd and Price 1987). Transition to school involves spending more time with peers, which also increases the various peer-related demands such as forming new friendships and gaining acceptance in a new peer group (Coplan and Arbeau 2008). This transition can be both challenging and stressful for a child (Niesel and Griebel 2007), since children need to deal with increasing amounts of successes and failures not only in the academic area, but also in peer relations (Campbell and Stauffenberg 2007). The first grade, in particular, has significant prognostic value in terms of internalizing problems not only in the first grade, but also during late childhood (e.g., Bukowski et al. 2010; Ialongo et al. 1995; Shin et al. 2016) and adolescence (e.g., Grover et al. 2007; Shin et al. 2016). However, as far as we know less attention has been paid on possible cumulative and snowballing effects of peer problems, at the beginning of school career on subsequent development of internalizing problems (Bukowski et al. 2010; Shin et al. 2016).

The Buffering Effect of Supportive Parenting

There are some theoretical and empirical grounds suggesting that the negative impacts of problematic peer relationships in childhood and adolescence depends—at least to some extent—on the quality of the child–parent relationship. For example, the stress-buffering model (Cohen and Wills 1985) suggests that social support attenuates the relation between negative life events and the risk for internalizing problems. This model states that when faced with stressful life events, individuals with greater support are less likely to show internalizing symptoms than individuals with lower levels of support. According to the model (Cohen and Wills 1985), parents can act as buffers, protecting their children against the potentially pathogenic influence of stressful condition, such as low peer acceptance. When children are exposed to low peer acceptance, supportive parenting being evident, for example, in warmth, affection, monitoring, and positive engagement, may play an especially important role in their adjustment and behavioral outcomes. Supportive parenting may enhance a child’s self-esteem, self-efficacy, and confidence by increasing the child’s perception that he or she can cope effectively with low peer acceptance (Burton et al. 2004). In addition, supportive parents may help children in adapting problem-solving skills and coping strategies in the case of low peer acceptance, which then leads to lower internalizing problems (Hazel et al. 2014).

Previous cross-sectional studies have provided some support for the notion that supportive parenting may provide a buffer against the negative impacts of difficulties with peers on internalizing problems (e.g., Buckholdt et al. 2016; McDonald et al. 2010; McLachlan et al. 2010). For example, Buckholdt et al. (2016) found that positive maternal responses to a child’s sadness and anger gave protection against the negative association of poor peer relations (e.g., low peer acceptance) with feelings of loneliness in 10–12-year-old children (grades 4–6). In a study on 14-year-old adolescents (ninth grade), McDonald et al. (2010) found that supportive relationships with mothers and fathers provided a buffer against the association between poor peer relationships and adolescents’ internalizing difficulties. In another study, McLachlan et al. (2010) also reported that parental warmth played a protective role in the face of peer rejection, leading to a decrease in related anxious expectations among children and early adolescents (those aged 9 to 13).

There are also longitudinal studies demonstrating the buffering effect of supportive parenting against peer difficulties. For example, Patterson et al. (1989) found that high levels of maternal warmth before entering the first grade (at age 6), served as a protective factor against the psychosocial adjustment difficulties associated with peer rejection in the first grade. In another longitudinal study, Bowes et al. (2010) found that supportive parenting in terms of maternal warmth (measured when the child was 7 years old) buffered the negative consequences of peer victimization on children’s subsequent internalizing problems at ages 10 and 12. Furthermore, Hazel et al. (2014) found that maternal and paternal supportive parenting (evidenced as close and positive relationships with children) buffered against the negative effects of peer problems (e.g., feeling pressure from friends or fighting with them) on the development of depressive symptoms in middle childhood (third grade, 7–9 years old; sixth grade, 10–12 years old) to adolescence (ninth grade, 13–16 years old). Finally, in a study by Healy and Sanders (2018) on 6–12-year-old children who were followed over a nine-month period, supportive parenting (in the form of warm parenting), was found to protect children against ongoing victimization and predict lower levels of internalizing problems over time. There are, however, also studies that have found no evidence for the buffering role of supportive parenting. For example, Sentse and Laird (2010) followed early adolescents (mean age 12.9 years) over a 1-year period and found that although maternal support (i.e., acceptance, warmth and support) provide a buffer against the association between friendship conflict and the development of depression in early adolescence, it did not protect against the association between peer rejection and internalizing problems. Similarly, in a study by Bilsky et al. (2013) on 8–14-year-old students, over three time points 6 months apart, no evidence was found for interaction between supportive parenting and peer relationship problems (in the form of peer victimization), on depressive symptoms.

In the present study, we focused on the buffering role of supportive parenting against the long-term impacts of low peer acceptance after the critical transition to primary school. More specifically, we examined whether the impact of low peer acceptance on the level and subsequent development of internalizing problem behaviors moderated by supportive parenting. Middle childhood is an important period for a child’s development because it involves the transition to school, the development of independence, peer interaction and friendship, yet it continues to be strongly influenced by parents. Transition to school can be stressful and challenging for a child (Niesel and Griebel 2007), and as it can increase the risk of encountering failure in peer relations as well as academic areas (Campbell and Stauffenberg 2007), parental support at the start of the school year in particular can be assumed to play a highly significant role in a child ‘subsequent adjustment outcomes. In the present study, supportive parenting was approached from two viewpoints. First, in line with previous studies on the topic we operationalized supportive parenting as one of the most often investigated aspects of parenting style dimensions, i.e. affection or warmth. Second, another aspect of supportive parenting, namely lack of psychological control, was also investigated as a possible moderator of peer relations. Psychological control refers to a style of parenting that is characterized by attempts to control the child’s emotions and behavior through psychological means, such as withdrawal of love and induction of guilt (Barber 1996; Hart et al. 2003) and it can be seen to be an indicator of a lack of parental support (for a review, see McKee et al. 2008; see also Nanda et al. 2012; Zarra-Nezhad et al. 2014). Although previous research has mostly focused on the role of parental support in terms of warmth or positive affection as a buffer, i.e., the positive side of parenting, also lack of negative side of parenting can be considered to function as part of supportive parenting. For example, psychological control has been shown to be related to various forms of internalizing problems among children and adolescents (Barber 1996; Aunola and Nurmi 2005). Previous literature suggests that although parents show a high level of affection and are supportive, if they use psychological control in parenting, high affection and warmth may no longer provide protection problem behaviors (see, e.g., Aunola and Nurmi 2005; Aunola et al. 2018). Thus, when considering parents’ buffering role against peer problems and other stressors, it might be important to examine not only the role of the positive side of parenting, but also focus on the lack of the negative side of parenting (i.e., psychological control).

Besides examining the moderating role of parenting against the impacts of low peer acceptance on the level and subsequent development of internalizing behaviors, gender differences in this phenomenon were investigated. Thus far, little is known about whether supportive parenting serves a similar or different role as a moderator of peer acceptance for boys and girls. It has been suggested that, from middle childhood on, girls generally place more emphasis on friendships (Maccoby 1998), become more easily depressed (e.g., Hale et al. 2005), and are more likely to suffer from negative interpersonal experiences such as low peer acceptance. However, boys have been reported to spend more time in large peer groups than girls (for a review, see Rose and Rudolph 2006), which may make them more likely to suffer from stressors at the group level (e.g., low peer acceptance) than girls (Lewin et al. 1999). It has also been suggested that boys who are rejected by their peers are more likely to feel anxious and expect future rejection than girls with a similar experience (London et al. 2007). However, to our knowledge, gender differences in terms of the role of supportive parenting as a buffer against low peer acceptance have not yet been examined.

Current Study and Hypotheses

Three research objectives were addressed in this study. (1) To what extent does peer acceptance at first grade predict a child’s initial level and subsequent development of internalizing problem behaviors during their primary school years (i.e., across grades 1, 2, 3, 4, and 6)? Based on previous literature, we expected that experiences with peers at the beginning of a child’s school career (first grade) would not only be associated with the initial level of internalizing problems, but also predict the subsequent developmental trend of these (e.g., Bukowski et al. 2010; Shin et al. 2016). More specifically, low peer acceptance in the first grade was expected to be positively associated with both the level and the developmental trend (i.e., either larger increase or smaller decrease across time depending on the direction of the trend) of internalizing problems across grades 1 to 6. (2) To what extent is the impact of peer acceptance on the level and developmental trend of a child’s internalizing problem behaviors moderated by affection and psychological control exhibited by her/his mother and father? Based on the stress-buffering model (Cohen and Wills 1985), our assumption was that supportive parenting in terms of high parental affection and low psychological control in first grade would provide a buffer against the impacts of low peer acceptance on children’s subsequent internalizing problems across grades 1 to 6. Since childrearing approaches have been reported to be relatively stable in early childhood (e.g., Roberts et al. 1984) and over longer periods of time (i.e., from aged 7–8 years until children reach 15–26-years of old; McNally et al. 1991), parenting styles in first grade was expected to be predictive across primary school years. (3) To what extent are the buffering effects of high parental affection and low psychological control different for boys and girls? Because of the lack of previous research concerning the moderating role of the child’s gender, no hypotheses were set.

Method

Participants and Procedure

The present study is part of a longitudinal study (Lerkkanen et al. 2006–2016) that followed a larger sample of 1880 children from preschool up to ninth grade in different regions of Finland. In the present sub-study, a target sample of 608 children (264 girls, 344 boys) was drawn from this larger sample in order to decrease the data collection demands placed on teachers. The target sample was determined by randomly selecting a small number of students from each first-grade classroom. Typically, three children were selected from each classroom (M = 2.53, SD = 0.84), but this number varied from one to six depending on the size of the class.

Comparisons between the random target sample and the larger sample revealed no differences between the samples in terms of the type of family and parental level of education. The vast majority of children in the target sample (76%) came from nuclear families, 12% were from single-parent families, and 12% from blended families. A total of 25% of the children’s mothers had a master’s degree or higher, 37% had a bachelor’s or vocational college degree, 31% had completed secondary education, and 7% had no qualification beyond secondary school. A total of 22.80% of the children’s fathers had a master’s degree or higher, 33% had a bachelor’s or vocational college degree, 34.70% had completed secondary education, and 9.50% had no qualification beyond secondary school. The sample was fairly representative of the level of education among the general population in Finland (Statistics Finland 2007).

Children’s levels of internalizing problems were rated by their teachers in grade 1 (April 2008), grade 2 (April 2009), grade 3 (April 2010), grade 4 (April 2011), and grade 6 (April 2013). Mothers and fathers were asked to complete questionnaires at home concerning their parenting styles, independently and without conferring, when their children were in the first grade. Information on parenting styles was available to 432 of the mothers (71.05%) and 281 of the fathers (46.21%).

Measures

Peer Acceptance

Students’ peer acceptance was measured in the first grade using a sociometric nomination procedure. Accepted children were those who were liked or preferred by peers or were most often named as a friend. Students were asked to circle (from a list of all the students in the class) the names of those with whom they most liked to spend time with during breaks (see Kiuru et al. 2015). The number of permitted nominations varied by class size; if the number of students in the class was 12 or greater, then three nominations were permitted. If the number of students in the class was 8–11, two nominations were permitted. If the number of students in the class was fewer than eight, only one nomination was permitted. Peer acceptance represents the number of nominations each student received, standardized by class size. Sociometric nominations have been shown to provide valid, stable, and reliable assessments of peer acceptance during early childhood (Terry 2000).

Internalizing Problem Behaviors

During grades 1, 2, 3, 4, and 6, teachers were asked to rate each child on a 3-point rating scale (1 = does not apply, 2 = applies partly, 3 = certainly applies) using the Strengths and Difficulties Questionnaire (SDQ, for 4–17 years old; Goodman 1997). The Emotional Symptoms subscale (five items; e.g., “Is often unhappy, down-hearted or fearful”) of SDQ was used as an indicator of internalizing problem behaviors. Mean scores for this subscale were calculated at each of the time points. The Cronbach’s α of internalizing problems at different time points ranged from 0.74 to 0.76.

Parental Affection and Psychological Control

Mothers’ and fathers’ levels of affection and psychological control were measured using a revised Finnish version (Aunola and Nurmi 2004) of Block’s Child-Rearing Practices Report (CRPR; Roberts et al. 1984). Parents were asked to rate 14 items measuring affection and psychological control on a 5-point scale (1 = does not fit me at all; 5 = fits me very well). The affection dimension (10 items, e.g., “I often show my child that I love him/her”) measured parental warmth and responsiveness. The psychological control dimension (four items, e.g., “I believe my child should be aware of how much I have done for him/her”) measured parental attitudes stirring up guilt on the part of the child and expressing disappointment. Cronbach’s α regarding these dimensions was, respectively, 0.83 and 0.71 for mothers and 0.83 and 0.74 for fathers.

Attrition analyses were conducted to determine if children whose mothers did (n = 432, 71%) and did not (n = 176, 29%) return the questionnaires differed in the study variables. The results of independent samples t-tests showed that there were no significant differences between the two groups in peer acceptance (t (588) = 0.98 p = .32). Results further showed no significant differences between the two groups in internalizing problems except at T2 (t (482) = 2.52, p = .01) which internalizing problems were higher in children whose mothers did not return the questionnaire. Next, attrition analyses were conducted to determine if children whose fathers did (n = 281, 46%) and did not (n = 327, 54%) return the questionnaires differed in the study variables. The results of independent samples t tests showed that there were no significant differences between the two groups in peer acceptance (t (588) = − 1.27, p = .20). Results furthermore showed no significant differences between the two groups in internalizing problems except at T3 (t (451) = 3.12, p = .001) which internalizing problems were higher in children whose fathers did not return the questionnaire.

Analysis Strategy

In order to examine the impact of peer acceptance on the development of children’s internalizing problems and to determine the extent to which these associations vary depending on parental affection and psychological control, we employed latent growth curve modelling (LGCM) in the data analysis. In these models, the latent intercept variable representing the initial level of internalizing problems at first grade was estimated by setting factor loadings of observed internalizing problem variables at different time points to be 1. A latent slope variable representing a linear trend in internalizing problems was estimated by setting factor loadings of observed internalizing problem variables at T1–T5 to be 0, 1, 2, 3, and 5, respectively. Setting the first slope factor loading at T1 at zero meant that the intercept described the level of internalizing problems at the onset of the study. The residual variances of the observed variables were freely estimated. Next, we built conditional latent growth models by incorporating peer acceptance, parenting variables (separate models including mothers’ and fathers’ variables), and their interaction terms as predictors of inter-individual differences in the level and developmental trend of internalizing problems. Finally, we conducted multi-group analyses by dividing the sample into two groups, girls and boys, and then investigated whether the effects of predictors on the level of and changes in internalizing problems were different in these two groups.

To examine the interactions between parenting variables and peer acceptance in terms of children’s internalizing problems in greater detail, we used the Johnson–Neyman technique (Preacher et al. 2006). This procedure produces the precise regions of the continuum of parenting variable for which the regression slope of peer acceptance on the developmental trend (i.e., slope) of internalizing problems (X→Y) is estimated to be significantly different from zero. In these analyses, the regions of significance represent the range of moderator variable z (i.e., particular parenting variable), within which the association of independent variable x (i.e., peer acceptance) with dependent variable y (i.e., slope of internalizing problems) is significantly different from zero at the chosen α-level. As a result of the method, two values describing the upper and lower boundaries of the region of significance are calculated and plotted. The effect of peer acceptance on the slope of internalizing problems is significant at values of the parenting variable that are less than the lower boundary and greater than the upper boundary, while the effect of peer acceptance is non-significant at values of the parenting variable falling within the region.

The fit indices utilized included the Comparative Fit Index (CFI), the Tucker-Lewis Index (TLI), and the Root Mean Square Error of Approximation (RMSEA). For the CFI and the TLI, values between 0.90 and 1.00 indicate a good fit between the model and the data. An RMSEA value of 0.05 or less indicates a good fit, while values up to 0.08 indicate a reasonable error of approximation. The fits of the nested models were compared by examining the significant changes in the Chi square values and degrees of freedom. All the analyses were carried out using the Mplus statistical package (Version 7.0; L.K. Muthén and Muthén 1998–2015). The results of Little’s MCAR test for the study variables revealed that the data was completely missing at random: (χ2 (398) = 431.62, p = .12). The parameters of the models were estimated using the full information maximum likelihood (FIML) estimation with standard errors that are robust to non-normality (MLR estimator) (Muthén and Muthén 1998–2010). This method allowed the use of all observations in the data set to estimate the parameters in the models without the imputation of missing values. The means, standard deviations, internal consistencies, and correlations for the observed variables are presented in Table 1.

Table 1 Means, standard deviations (SD), and correlations for the study variables

Results

Developmental Trend of Children’s Internalizing Problem Behaviors

The results suggested a good model fit for the linear latent growth curve model for internalizing problems [χ2 (df = 10, N = 565) = 31.59, p < .001; RMSEA = 0.06; CFI = 0.95; TLI = 0.95]. The estimated mean values of the intercept and slope were 1.32 (SE = 0.02, p < .001) and − 0.02 (SE = 0.01, p = .05), respectively, suggesting that, on average, internalizing problems decreased statistically significantly across the school years. Also, the variance of intercept (variance = 0.10, SE = 0.01, p < .001) and slope (variance = 0.02, SE = 0.00, p < .001) were statistically significant, indicating individual differences in internalizing problems when children first entered school, as well as in the developmental trend of internalizing problems.

The Role of Supportive Parenting and Peer Acceptance

Next, conditional latent growth models were constructed by adding mothers’/fathers’ parenting variables, peer acceptance, and their interaction terms as predictors of the intercept and slope of internalizing problems for the LGCM. Separate models were carried out for each parenting variable. To test our hypothesis concerning the differential effects of parenting and peer acceptance and their interaction with children’s internalizing problems depending on the child’s gender, we started the analysis by conducted multiple-group LGCMs. By applying this method, we compared freely estimated models (all predictive paths estimated separately among boys and girls) to constrained models (all predictive paths constrained to be equal across boys’ and girls’ groups) via a Chi square difference test to reveal whether the relationships among parenting variables, separately, peer acceptance, and the interaction terms of parenting variables and peer acceptance, on children’s initial level and developmental trend of internalizing problems would vary between the boys and girls in our sample. If no gender differences were evident, final analyses were carried out for the whole sample.

Parental Affection

The results showed that in the case of mothers’ affection, the freely estimated model was not statistically better than the constrained model (Δχ2 (6) = 0.68, p  = .99), suggesting that the associations of maternal affection, peer acceptance, and their interactions with children’s internalizing problems were similar for boys and girls. Consequently, a conditional growth model was applied to the whole sample. The results of the models including mothers’ affection are shown in Table 2. The model for mothers’ affection and children’s peer acceptance and internalizing problems indicated an acceptable fit (Table 2). First, the results showed that mothers’ affection was not related to the initial level or developmental trend of internalizing problems in children. Peer acceptance was negatively related to the level of internalizing problems (β = − 0.26, p < .001) but not to the developmental trend of it: the lower the level of peer acceptance, the higher the level of internalizing problems in children. The results showed, however, that the interaction term between peer acceptance and maternal affection was statistically significant for the trend (β = 0.12, p = .04) of internalizing problems. In other words, the strength of the relationship between peer acceptance and the development of internalizing problems varied as a function of maternal affection. To interpret this interaction further, the Johnson–Neyman technique was utilized. The results are shown in Fig. 1.

Table 2 Main and interaction effects of parenting styles and peer acceptance predicting children’s internalizing problems
Fig. 1
figure 1

The effect of low peer acceptance on the developmental trend of children’s internalizing problem behaviors as a function of maternal affection. The X-axis depicts a continuous range of maternal affection (moderator variable) values. The Y-axis depicts a continuous range of values for the adjusted effect of peer acceptance on the trend of internalizing problem behaviors. The straight plot line represents values of the adjusted effect that correspond to the full range of all continuous values of maternal affection requested (i.e., − 2 SD to + 2 SD). The curved lines above and below the plot line represent 95% confidence bands around the adjusted effect of peer acceptance on the trend of internalizing problems. Gray shaded areas denote regions (the maternal affection values below − 1.3 SD) where the two lines statistically significantly differ

As shown in Fig. 1, the effect of peer acceptance on the developmental trend of internalizing problems was statistically significant and negative only when maternal affection was low (values below − 1.2 SD): the higher the peer acceptance in first grade, the greater the decrease in internalizing problems across grades 1 to 6—and the lower the peer acceptance in first grade, the smaller the decrease in internalizing problems from grade 1 to 6—among children whose mothers showed low levels of affection in first grade. However, when the mother showed moderate or high levels of affection, peer acceptance showed no effect on the trend of internalizing problems.

Next, a model including fathers’ affection was created. The results showed that for paternal affection the freely estimated model was not statistically better than the constrained models (Δχ2 (6) = 5.74, p = .55), suggesting that the associations of fathers’ affection, peer acceptance and their interactions with children’s internalizing problems was similar for boys and girls. Consequently, conditional growth model was conducted to the whole sample. The results for paternal affection and children’s peer acceptance and internalizing problems indicated an acceptable fit (Table 2). The results showed that paternal affection in first grade was not related to the initial level or developmental trend of internalizing problems in children across grades 1 to 6. No statistically significant interaction effects between peer acceptance and paternal affection in first grade were found either.

Parental Psychological Control

The model for maternal psychological control showed that the freely estimated model (χ2 (df = 39, N = 586) = 62.12, p < .001; RMSEA = 0.04; CFI = 0.96; TLI = 0.95) fitted the data better than the constrained model (χ2 (df = 45, N = 586) = 74.34, p < .001; RMSEA = 0.04; CFI = 0.95; TLI = 0.94; Δχ2 (6) = 12.22, p = .04), suggesting that the estimated paths are different among boys and girls. The results of the model estimating paths separately among boys and girls showed two differences between boys and girls. First, the results showed that mothers’ psychological control was related to the initial level of internalizing problems only among girls: the higher the maternal psychological control in first grade, the higher the level of internalizing problems at first grade among girls (β = 0.27, p = .001) but not among boys (β = 0.03, p = .71). Comparing the freely estimated model with the model constraining only this one path (path from maternal psychological control to the level of internalizing problem behaviors) to be same for boys and girls, confirmed statistically significant difference between boys and girls (Δχ2 (1) = 3.96, p = .04).

Second, the interaction between peer acceptance and mothers’ psychological control on the developmental trend of internalizing problems was statistically significant only among boys (estimate for boys: β = − 0.21, p = .04 estimate for girls: β = 0.11, p = .05). Comparing the freely estimated model with the model constraining only this one path (path from the interaction term to the developmental trend of internalizing problem behaviors) to be same for boys and girls, confirmed statistically significant difference between boys and girls (Δχ2 (5) = 11.45, p = .04). The results gained from using the Johnson–Neyman technique (Fig. 2) showed that when mothers showed a high level of psychological control (over + 2.5 SD), the effect of peer acceptance on the trend of internalizing problems among boys was negative: the lower the peer acceptance in first grade, the less the decrease in internalizing problems among boys across grades 1 to 6. When maternal psychological control was low (below − 2.3 SD), the effect of peer acceptance on the trend of internalizing problems among boys was reversed: the lower the peer acceptance in first grade, the greater the decrease in internalizing problems across grade 1 to 6 among boys.

Fig. 2
figure 2

The effect of low peer acceptance on the developmental trend of boys’ internalizing problem behaviors as a function of mothers’ psychological control. The X-axis depicts a continuous range of maternal psychological control (moderator variable) values. The Y-axis depicts a continuous range of values for the adjusted effect of peer acceptance on the trend of internalizing problem behaviors. The straight plot line represents values of the adjusted effect that correspond to the full range of all continuous values of maternal psychological control requested (i.e., − 2 SD to + 2 SD). The curved lines above and below the plot line represent 95% confidence bands around the adjusted effect of peer acceptance on the trend of internalizing problems. Gray shaded areas denote regions (the psychological control values above + 2.5 SD and below − 2.3 SD) where the two lines statistically significantly differ

When a similar model including psychological control was tested for fathers, the results showed that the freely estimated model was not statistically better than the constrained model (Δχ2 (6) = 3.60, p = .71), suggesting that the associations of fathers’ psychological control, peer acceptance and their interactions with children’s internalizing problems was similar for boys and girls. Consequently, the model for the whole sample was estimated. The fit of the model was good (see Table 2). The results showed that fathers’ psychological control was positively related to the developmental trend of children’s internalizing problems (β = 0.14, p = .05) but not to the level: the lower the level of paternal psychological control in first grade, the greater the decrease in internalizing problems from grades 1 to 6; and the higher the level of paternal psychological control in first grade, the smaller the decrease in internalizing problems from grades 1 to 6. No statistically significant interaction effects between paternal psychological control and peer acceptance were found.

Discussion

The principal aim of this study was to obtain a deeper understanding of the role of supportive parenting in terms of buffering the negative impacts of low peer acceptance on children’s internalizing problem behaviors from grades 1 to 6. We investigated whether peer acceptance at first grade after the critical transition to school predicts internalizing problems in subsequent primary school years. We also explored the moderating role of parental affection and psychological control on the relation between a child’s peer acceptance and internalizing problems, and possible gender differences in this.

Association Between Peer Acceptance and Internalizing Problem Behaviors

We hypothesized that low peer acceptance in the first grade is not only associated with a higher level of internalizing problems, but also positively predicts the developmental trend of these problems being evident either as greater increase of internalizing problems across time, or less decrease in these, depending on the general trend of internalizing problems. This hypothesis was partly supported. Low peer acceptance was associated with initial higher levels (but not with subsequent change) of internalizing problems at the beginning of primary school. This result was in line with previous studies showing that low peer acceptance is associated with various kinds of internalizing problems (e.g., Gifford-Smith et al. 2005; Klima and Repetti 2008; Ladd 2006; Menting et al. 2015). According to self-determination theory (Deci and Ryan 2000), relatedness (i.e., a feeling of belonging and being significant in the eyes of others) is a basic psychological need that is essential for healthy development and well-being. For example, lack of peer acceptance can negatively affect the children’s self-perception by making them believe that they are not worthy of being accepted (Ladd and Troop-Gordon 2003). Low peer acceptance may also undermine children’s sense of relatedness, which contributes to increased feelings of loneliness and other internalizing symptoms such as anxiety, depression and withdrawal (Klima and Repetti 2008; Ladd 2006; Parker et al. 2006). Because at the level of the whole sample low peer acceptance was associated only with the initial level, and not with the developmental trend of internalizing problem behaviors, this association could also be interpreted the other way around. It is possible that internalizing problem behaviors impact on peer acceptance rather than vice versa. For example, if a child has a high level of internalizing problems, it is possible that she/he is also more prone to withdraw from the peer group, which then affects how peers interact with him/her. Although low peer acceptance did not predict subsequent development of internalizing problems at the level of the whole sample, our results nevertheless demonstrated that the association of lack of peer acceptance on the developmental trend of internalizing problem behaviors was dependent on the level of maternal supportive parenting.

The Moderating Role of Supportive Parenting

We hypothesized that a high parental affection and low level of psychological control would moderate the relation between peer acceptance and internalizing problems by buffering the negative effect of low peer acceptance on the level and developmental trend of a child’s internalizing problems. The results concerning maternal affection supported our hypotheses by demonstrating that maternal affection protected children against the detrimental impacts of low peer acceptance: low peer acceptance predicted less positive developmental outcomes only when the level of maternal affection was low. When mothers exhibited moderate or high maternal affection, low peer acceptance showed no effects on the developmental trend of internalizing problems. This result is consistent with the stress-buffering model (Cohen and Wills 1985), and previous findings suggesting that parental affection and support can act as a buffer against the negative impacts of poor peer relationships on internalizing problems over time (Bowes et al. 2010; Hazel et al. 2014; Healy and Sanders 2018); Patterson et al. 1989). This result was, however, inconsistent with the findings of Bilsky et al. (2013) and Sentse and Laird (2010) that reported maternal support did not moderate the negative impact of peer relationship difficulties on internalizing problems. The discrepant results may be due to the fact that these two previous studies focused on adolescents, whereas the present study focused on primary school children. As children begin their transition into adolescence, supportive parenting in terms of parental warmth and affection may decrease (e.g., McGue et al. 2005). Thus, it is possible that parents are less influential as children age (Collins et al. 1997; see also, Hazel et al. 2014).

The results of the present study can be interpreted in the light of the stress-buffering model (Cohen and Wills 1985) according to which social support, or positive relationships with others (e.g., parents), protects against the potentially pathogenic influence of stressful life situation (such as poor peer relationships). According to this model, supportive parents help their children in adapting problem-solving skills and coping strategies and enhance their confidence, self-efficacy, and self-esteem, which may then lead to lower internalizing problems in the case of low peer acceptance. A second possible explanation for our results is that since warm and positive parenting is characterized by better parent–child communication and more social support (Ranson and Urichuk 2008), this kind of mothering may provide the emotional support that low peer-accepted children need. For example, children who have warm and positive mothers may be more willing to disclose and share their feelings of rejection, which makes mothers more aware of what is going on in their young children’s classrooms. Thus, warm and engaged mothers may invest more time and effort in helping their first-grade children get involved in extracurricular activities where they are able to find alternative opportunities to gain social skills outside of school. Third, it is possible that warm and supporting mothers are more attuned to their children’s difficulties in fitting into classroom social dynamics and do a better job of helping children cultivate personal strengths that do not depend so heavily on social feedback from peers. These mothers may be more likely to support and instruct their children in the use of coping strategies that provide them with the tools they need to deal with the pressures of negative feedback from peers and provide children with more positive and adaptive appraisals of stress (Hazel et al. 2014). This might then protect children from internalizing symptoms and provides them with alternative ways to build their self-esteem. It is nevertheless notable that we did not empirically test the above-discussed possible mechanisms underlying the buffering effect of high maternal affect. Future studies are needed to empirically test these alternative explanations.

Although our results showed that maternal affection buffered against the impacts of low peer acceptance on the developmental trend of internalizing problems, no such result was found concerning fathers’ affection. Thus, the findings suggest that it is particularly maternal rather than paternal affection that has buffering effect on children’s subsequent development of internalizing symptoms over the school years. This result is understandable when taking into account that mothers are often the primary caretakers, taking primary responsibility for children in families, and spend more time with their children (Lamb and Oppenheim 1989; Renk et al. 2003). Thus, their parenting styles may naturally play a more important role in a child’s emotional development than the father’s parenting does (Zarra-Nezhad et al. 2015). Alternatively, perhaps children tend to be more open to maternal than paternal influence (Aunola and Nurmi 2005) because interactions between mother and child are typically characterized by greater responsiveness, warmth and intimate exchanges than interactions between father and child (Collins and Russell 1991). However, because of the relatively small sample size of fathers in our sample, further research is needed to detect the buffering role of fathers’ affection.

The results concerning the moderating role of parental psychological control also partly supported our hypothesis in the case of mothers: when mothers showed a high level of psychological control, low peer acceptance was associated with relatively higher subsequent levels of internalizing problems among boys. On the contrary, when maternal psychological control was low, a lack of peer acceptance among boys was related to decreased rather than increased internalizing problems. For girls, this kind of result was not found suggesting some gender differences in the impact of low peer acceptance on internalizing problems across the school years depending on the level of maternal psychological control. The result for boys is consistent with the stress-buffering model (Cohen and Wills 1985), and our assumption, suggesting that low maternal psychological control as one form of supportive parenting can protect against the negative impacts of low peer acceptance. When interpreting the results in the light of the stress-buffering model (Cohen and Wills 1985), non-supportive mothering in terms of psychological control (i.e. controlling the child’s emotions and behavior through psychological means) might be assumed to decrease the child’s confidence, self-efficacy, and self-esteem in the face of stressful life situation such as low peer acceptance, and, by this way, foster rather than prevent internalizing problems.

Another possible explanation for this result is that high levels of maternal psychological control ignore the child’s fundamental need for autonomy (Soenens and Vansteenkiste 2010), consequently intensifying children’s symptoms of internalizing problems (Barber and Xia 2013; El-Sheikh et al. 2010). It is also possible that boys whose mothers use psychological control develop a negative sense of self as a consequence of maternal emotional unavailability, which then leads them to adopt maladaptive coping strategies in peer relationships. They may, for example, feel increased pressure to win social acceptance and get attention from peers, and may lack effective coping skills to deal with low levels of acceptance from peers. This then may put them at a higher risk of internalizing symptoms, such as low self-esteem and depression (for a review, see Soenens and Vansteenkiste 2010). On the other hand, boys whose mothers show very low levels of psychological control may help their sons to develop a stable, secure and positive sense of self, which then helps them to cope adaptively in the case of peer rejection. Although the combined effect of maternal psychological control and peer acceptance was only evident among boys, the main effect of mothers’ psychological control on the level of internalizing problems was found only among girls. Overall, these results suggest that although girls may be more prone to the negative impacts of maternal psychological control in general, boys showing low peer acceptance, in particular, may benefit from supportive mothering in terms of low psychological control. One explanation for this pattern of results is that because in middle childhood boys spend more time in large peer groups than girls (for a review, see Rose and Rudolph 2006), they suffer more by being rejected by their peers than girls (London et al. 2007) and, thus, also benefit more from compensating supportive mothering.

Another possible mechanism that might explain this result is that in middle childhood boys and girls use different coping strategies when dealing with peer difficulties (Visconti and Troop-Gordon 2010). For example, girls respond in a more prosocial manner to hypothetical peer conflict situations than boys in middle childhood (for a review, see Rose and Rudolph 2006). Similarly, in case of peer rejection, girls are more likely to use strategies to further the discussion of emotions within a situation, walk away from the conflict, or seek support if necessary, whereas boys engage less in prosocial behaviors, i.e., they show less control over negative emotions and less effective coping strategies (Fabes and Eisenberg 1992; Smith et al. 2001), which may put them at a higher risk of behavioral problems. Consequently, boys may benefit more from maternal supportive parenting, which provides them with an alternative source of emotional support and helps them to develop adaptive and effective coping skills to deal with experienced low levels of peer acceptance and decrease the risk of internalizing problems. Future studies on gender differences are needed, however, to provide support for these kinds of conclusions.

As was the finding in the case of parental affection, our results further revealed that low paternal psychological control did not protect children against the detrimental impacts of low peer acceptance either. Overall, the results suggest that mothers and fathers may play a different role when it comes to buffering against low peer acceptance. Although both mothers’ and fathers’ psychological control was found to be associated with high levels of internalizing problems among children, it was maternal rather than paternal affection and psychological control that showed a buffering effect on children’s subsequent development of internalizing symptoms over the school years.

Limitations

This study had some limitations that must be taken into account. First, peer acceptance and parenting styles was measured only once, at the beginning of the study. Peer status and parenting may, however, change over the school years. Consequently, further studies are needed to examine the developmental dynamics between peer acceptance and children’s internalizing problems as well as the buffering role of parenting in this development. Second, only parents’ self-report measures were used for measuring parenting styles and these self-reports may be impacted by social desirability. Third, using a teacher report for internalizing problem behaviors may limit assessment of this kind of behaviors in one context only (i.e., school). On the other hand, it is also a strength of our study that we were able to show rater-independent associations between the investigated constructs, that is, parenting styles (reported by parents), peer acceptance (rated by peers) and internalizing problems (rated by teachers). When common method bias could not affect our results, our longitudinal findings may be more robust. Fourth, the fathers’ sample (n = 281) was relatively small compared to the mothers’ sample (n = 432), which may have reduced the statistical power of some of the analyses concerning fathers. Fifth, although it is possible that parental affection and psychological control have also interactive effect on children’s developmental outcomes (e.g., Aunola and Nurmi 2004; Aunola et al. 2018), in the present study these effects were not empirically tested due to the relatively small sample and other moderator variables. In further studies with larger sample sizes the four-way interaction effects of Affection × Psychological control × Gender × Peer acceptance should also be tested in order to get more deeper and broader insight on the complex interplay of parenting and peer acceptance when predicting children’s internalizing problems. Finally, all the results were obtained on Finnish children and further research in other regions of the world is needed in order to generalize the findings.

Conclusions

Overall, the results of this study suggest that low peer acceptance at the very beginning of primary school is associated with higher levels of internalizing problem behaviors from first grade to sixth grade. The first grade can be challenging and stressful for a child, because after the transition to school, peer relationships are largely reformed, and children face a greater risk of encountering failure in peer relationships. Peer status has been found to be relatively stable in the elementary school years and to predict later social and emotional development, and school adjustment (for a review, see Jiang and Cillessen 2005). Therefore, low peer acceptance in this period can have a long-term negative impact on the child’s further social and emotional development, manifested in the form of increased internalizing problems, for example. However, our results showed that mothers’ supportive parenting style, characterized by high levels of affection and sensitive caregiving, can protect children from the negative effects of low peer acceptance throughout the primary school years.