Bodyweight
There were no obvious effects of ELS on the bodyweight gain from PND21 to PND50. From PND21 to PND64 and PND245, the siblings in the control group with no separation (SH) gained more weight than the other groups (Table S2; Separation group × PND). Consequently, the overall bodyweight of the SH group was higher than the other three groups on PND245 at the end of the experiment (multiple comparisons: NS < SH; estimate ± SE = 20.64 ± 5.85, t = 3.53, p = 0.02, CS < SH; estimate ± SE = 18.58 ± 5.45, t = 3.41, p = 0.030, IS < SH; estimate ± SE = 24.34 ± 5.25, t = 4.64, p < 0.001). In males, from PND21 to PND245, the SH group gained more weight than other groups (Figure 3a, Table S3; Separation group × PND). The bodyweight of the SH group was heavier than that of the IS group on PND245 (Figure 3a, multiple comparisons: IS < SH; estimate ± SE = 25.95 ± 7.22, t = 3.60, p = 0.02). In females, from PND21 to PND64 and PND245, the SH group gained more weight than the other groups (Figure 3b, Table S4; Separation group × PND). The bodyweight of the SH group tended to be heavier than that of the NS group on PND64 (Figure 3b, multiple comparisons: NS < SH; estimate ± SE = 30.88 ± 9.35, t = 3.30, p = 0.08).
Behavior
Effects of separation and PND on freezing, rearing, and grooming behaviors in all individuals
Overall, the time spent freezing for all individuals was not different between groups (Table S5a; Separation group). By contrast, the time spent freezing was shorter in this order: PND50 < PND21 < PND245 (Table S5a; PND). From PND21 to PND50, the time spent freezing decreased less in the NS group than in the SH group, whereas it decreased more in the CS group than in the SH group (Table S5a; Separation group × PND). From PND21 to PND245, the time spent freezing decreased in the NS and CS groups while it increased in the SH group (Table S5a; Separation group × PND). On PND21, the time spent freezing tended to be shorter in the IS group than in the NS group (Figure 4a, multiple comparisons: IS < NS; estimate ± SE = 0.96 ± 0.33, z = 2.92, p = 0.06). On PND50, the time spent freezing was shorter in the CS and IS groups than in the NS group (Figure 4a, multiple comparisons: CS < NS; estimate ± SE = 1.70 ± 0.36, z = 4.74, p < 0.001, IS < NS; estimate ± SE = 1.09 ± 0.34, z = 3.25, p = 0.02), and the time spent freezing in the CS group tended to be shorter than that in the SH group (Figure 4a, CS < SH; estimate ± SE = 0.97 ± 0.34, z = 2.82, p = 0.09). On PND245, the time spent freezing was shorter in the CS group than in the SH and NS groups (Figure 4a, multiple comparisons: CS < SH; estimate ± SE = 1.30 ± 0.33, z = 3.90, p = 0.002, CS < NS; estimate ± SE = 1.16 ± 0.35, z = 3.28, p = 0.02). Especially, the time spent freezing in the SH group was long, with a maximum value of 117 s for PND245.
The total number of rearing behaviors was higher in the CS and IS groups than in the SH group (Table S5b; Separation group). The total number of rearing behaviors was higher on PND50 than on PND21 (Table S5b; PND). From PND21 to PND50, the total number of rearing increased less in the CS and IS groups than in the SH group, whereas from PND21 to PND245, it increased larger in the NS and CS groups than in the SH group (Table S5b; Separation group × PND). On PND21, the total number of rearing behaviors was higher in the IS group than in the SH and NS groups (Figure 4b, multiple comparisons: IS > SH; estimate ± SE = -0.69 ± 0.22, z = -3.11, p = 0.03, IS > NS; estimate ± SE = -0.96 ± 0.25, z = -3.84, p = 0.002), and tended to be higher in the CS group than in the NS group (Figure 4b, multiple comparisons: CS > NS ; estimate ± SE = -0.75 ± 0.26, z = -2.84, p = 0.08). On PND50, the total number of rearing behaviors did not differ between these separation groups (Figure 4b, multiple comparisons: NS < SH; estimate ± SE = 0.49 ± 0.23, z = 2.16, p =0.55, CS > SH; estimate ± SE = -0.07 ± 0.20, z = -0.33, p = 1.0, IS < SH; estimate ± SE = 0.08 ± 0.19, z = 0.41, p = 1.0, CS > NS; estimate ± SE = -0.56 ± 0.24, z = -2.39, p = 0.31, IS > NS; estimate ± SE = -0.41 ± 0.23, z = -1.81, p = 1.0, IS < CS; estimate ± SE = 0.15 ± 0.20, z = 0.73, p = 1.0). On PND245, it was higher in the CS group than in the SH and NS groups (Figure 4b, multiple comparisons: CS > SH; estimate ± SE = -0.91 ± 0.23, z = -3.92, p = 0.002, CS > NS; estimate ± SE = -0.70 ± 0.25, z = -2.80, p = 0.09). Moreover, the total number of rearing behaviors tended to be higher in IS group than in SH group (Figure 4b, multiple comparisons: IS > SH; estimate ± SE = -0.67 ± 0.23, z = -2.94, p = 0.06).
Overall, the total number of grooming behaviors tended to be higher in the IS group than in the SH group and most evidently on PND50 than on PND21 (Table S5c; Separation group, PND). By contrast, the total number of grooming behaviors were not different between separation groups on each observation day (Figure 4c).
Sex differences in the effects of separation and PND on behaviors
Regarding fixed effects, males showed similar patterns to the case of all individuals pooled as they showed no differences in the time spent freezing between the separation groups, whereas females did not. Here the time females spent freezing tended to be longer in the NS group than in the SH group (Table S7a; Separation group). Similar to the case of all individuals pooled, the time males spent freezing showed the order: PND50 < PND21 < PND245, whereas that of females did not differ between PND21 and PND245 (Table S6a, Table S7a; PND). Conversely, the effects of separation group × PND on the time spent freezing males and females followed a similar pattern to the case of all individuals. On PND21, the time spent freezing of males did not differ between the groups. On PND50, the time spent freezing was shorter in the CS group than in the SH and NS groups and tended to be shorter in IS group than in NS group (Figure 5a, multiple comparisons: CS < SH; estimate ± SE = 1.48 ± 0.43, z = 3.44, p = 0.01, CS < NS; estimate ± SE = 1.99 ± 0.46, z = 4.37, p < 0.001). On PND245, the time spent freezing was shorter in the CS group than in the SH and NS groups (Figure 5a, multiple comparisons: CS < SH; estimate ± SE = 1.60 ± 0.41, z = 3.91, p = 0.002, CS < NS; estimate ± SE = 1.41 ± 0.44, z = 3.21, p = 0.02). By contrast, the time females spent freezing was not different between separation groups on each observation day (Figure 5b).
The total number of rearing of males was higher in the CS and IS groups than in the SH group, like the case of all individuals pooled, whereas those of females were lower in the NS group than in the SH group (Table S6b, Table S7b; Separation group). In males and females, the total number of rearing was higher on PND50 than on PND21 (Table S6b, Table S7b; PND). The effects of separation group × PND on the total number of rearing of males and females had a similar pattern to the case of all individuals pooled. The differences between groups according to age were more distinct in males than in females. On PND21, the total number of rearing of males was not different between the groups (Figure 5c), whereas that of females was higher in the IS group than in the NS group (Figure 5d, multiple comparisons: IS > NS; estimate ± SE = -1.49 ± 0.40, z = -3.72, p =0.004). On PND50, the total number of rearing of both males and females was not different between the groups. On PND245, the total number of rearing in males was higher in the CS and IS groups than in the SH group (Figure 5c, multiple comparisons: CS > SH; estimate ± SE = -1.50 ± 0.34, z = -4.47, p < 0.001, IS > SH; estimate ± SE = -1.27 ± 0.33, z = -3.81, p = 0.003).
The total number of grooming in males tended to be higher on PND50 than on PND21 like the case of all individuals pooled, whereas that of females did not (Table S6c). From PND21 to PND 50, the total number of grooming of males increased in the SH group but decreased in the IS group, whereas such a tendency was not observed in females (Table S6c). Also, in males and females, there were no differences between groups on each observation day (Figure 5e, Figure 5f).