Personality over ontogeny in zebra finches: long-term repeatable traits but unstable behavioural syndromes

Repeatability ('differential consistency')

If a single trait is stable over time, it has a high repeatability. Repeatability is the proportion of variance in a given trait that can be attributed to the difference between individuals (and is thus always a population level measure [7]). That is, when multiple measurements of the same individual are similar (low within-individual variance), while at the same time individuals differ substantially from each other (high between-individual variance), the behavioural variable is repeatable.

Average repeatability of behavioural traits has been shown to be around 0.4 [8] and 0.48 [9] in two meta-analyses. It can be influenced by several factors, for instance the type of behaviour, the time interval between tests, age or ontogenetic stage and the sex of the individuals [8]. Unfortunately, the vast majority of studies merely use two measurements per individual often quantified over short time scales. This was indicated in a meta-analysis [8], with a mode of 2 and a mean of 4.4 measures per individual over 759 studies [10] and only 9% of these studies covering a 'long' time interval, defined as over one year. This seems rather insufficient, as especially for traits with low repeatability, we expect the precision of the estimate to increase with higher sample sizes (number of subjects and/or test rounds) [9].

Behavioural syndrome structure ('structural consistency')

Not only can single traits show consistency, but also the between-individual correlations of two or more behavioural traits, forming a behavioural syndrome. Behavioural syndromes represent a functional coupling of traits [11] which has been argued to be either adaptive [12–14] or maladaptive [1, 15, 16]. They can emerge quickly and may change over time, or with critical life experiences [17–19], though a meta-analysis on behavioural syndromes found on average weak phenotypic correlations between behavioural traits (mean effect size r = 0.264 with 95% CI 0.210-0.316) [20].

Investigating personality over ontogeny

Although consistency is widely acknowledged as an important prerequisite for animal personality, studies often identify and describe animal personality and behavioural syndromes at a single point in time (and possibly their implications e.g. on fitness, or correlations to physiological parameters [21]), leading to a gap in our understanding how animal personality develops. Even when personality is tested over time, the test-retest intervals are mostly short [8] and not representative of the species' life span and different life history stages [6]. So far, most studies only dealt with adult animals, providing an incomplete picture of animal personality [but see e.g. [22–25]] and resulting in misinterpretations of the consistency of personality. This lack of studies across key stages of ontogeny probably arises from the intuitive contradiction between the two topics, with personality focusing on stability, while ontogeny deals with how organisms change over time [26]. However, there have been recent calls for a broader incorporation of developmental aspects (measurements over different ontogenetic stages) into animal personality research [27–29], arguing that this could reveal the proximate mechanisms underlying personalities. As recent findings indicate mixed evidence with regard to the consistency of personality over ontogeny ((partly) inconsistent: [15, 17, 30], consistent: [23, 31]), this challenges the assumption that personality is fixed across a lifetime.

Even though the literature focusing on ontogenetic stages currently increases significantly, few studies assess more than two (or three) personality traits (but see e.g. [32]). Furthermore, only little is known about how behavioural syndrome structure changes over time. Juveniles are facing different challenges (ecological or social niches) compared with adults (e.g. their behaviour should focus on survival rather than reproduction), hence adaptations needed for early life stages might disappear later [29], resulting in behavioural syndromes either forming or decoupling during adolescence. If the behaviours forming a syndrome are regulated by a common underlying physiological mechanism, they might be linked particularly tightly, potentially imposing a developmental constraint and thus limiting optimal behaviour, because single traits cannot change independently [15]. Indeed, behavioural syndrome structure was shown to be consistent over time (two weeks) [32] or similar across two cohorts [33].

So far, the existence of animal personality and behavioural syndromes in general is well documented, leading to a current shift to investigate the causes and consequences of personality [34, 35]. One of the most extensively studied model organisms is the zebra finch, Taeniopygia guttata. In this species, several personality traits have been identified, such as activity [22, 36, 37], aggression [37, 38], fearfulness [36], exploratory behaviour and struggling rate [37, 39] and boldness/neophobia [37, 40] including their heritability [40]. Furthermore, the existence of behavioural syndromes has been documented in zebra finches [37], as well as implications of personality type for fitness correlates e.g. reproductive success [41, 42].

Aims of this study

With our study, we intended to address the following widespread shortcomings regarding consistency of animal personality: (a) testing over short intervals or within single developmental stages, neglecting (b) the aspect of multidimensionality (assessing only one or two traits) and (c) the development of behavioural syndrome structure.

Maximum life spans of 1.3 to 5 years have been reported for different zebra finch populations in the wild [43], while under laboratory conditions, life spans may be at least 5 to 9 years [44, 45].

We repeated our tests after approximately one year ('mature adult' stage), which is a long time span in zebra finches, representing about 15 - 20% of the lifespan of captive animals. So far, the longest test-retest intervals in zebra finches we are aware of is three to seven months for exploration behaviour [39, 42], seven months for struggling rate [39], 1.5 years for neophilia [40] and even two years for activity levels [22].

We therefore provide the first extensive longitudinal study on a multitude of personality tests, following the same individuals over a substantial amount of their lifespan. We specifically asked the following questions: (1) How stable are the single personality traits within and across key ontogenetic stages; (2) which factors influence repeatability (trait, interval, life stage, sex); (3) how stable are the correlations between the traits (behavioural syndromes)?

Repeatabilities ('differential consistency')

Repeatability estimates are presented in Figure 1. Fearlessness (tonic immobility test, TI) and exploration (novel environment test, NE) were repeatable within all life stages and across the whole measurement period including all six test rounds. TI, however, was the only trait with a significant repeatability if the estimate was derived from only two data points with an interval of about one year ('long interval' using test rounds 1 and 6).

General activity (GA) was repeatable within the subadult stage and across all life stages, but not in young or mature adults.

In general, repeatability estimates for boldness (novel object test, NO) were very low (R ranging between 0.03 and 0.28), which might be due to the differences in the novel objects used in the test rounds (for details see table 1). NO in subadults and aggression (AG) in young adults showed significant p-values but confidence intervals overlapping with zero indicated non-significance. AG repeatability estimates were significant across life stages, but not within life stages.

Sex differences in repeatability

GA was highly repeatable for both males and females at the subadult life stage, but at the young adult stage only for females, and in mature adults in neither of the sexes. Aggression (AG) repeatability estimates of males and females showed confidence intervals overlapping with zero within all life stages. Across life stages, only female aggression was repeatable. NO repeatability estimates were similar for the sexes and non-significant in all life stages.

Behavioural syndrome structure ('structural consistency')

We additionally calculated behavioural syndrome structures separately for the sexes. Comparison with the dataset for both sexes revealed that the significant correlations in young adults seem to be highly driven by the female sex. While the correlations for females were similar, but with higher correlation coefficients (AG-NO: r = -0.69; CI -0.85 - (-0.41); p < 0.001; GA-NO: r = -0.61; CI -0.81 - (-0.28); p < 0.01 and GA-AG: r = 0.53; CI 0.17 - 0.76; p < 0.05, for more details see part 3 in additional file 1), these correlations were completely absent in the male subset. There were no significant correlations for males, except for a positive trend between exploration and aggression (r = 0.57; CI 0.18 - 0.80; p = 0.07) at the mature life stage.

Generalized linear mixed effects models

Results of the GLMMs are given in detail in part 4 (complete, 'inflated' datasets caused by zero- and maximum-inflated data) and part 5 ('reduced' datasets excluding these floor or ceiling effects) of additional file 1. Only in a few cases did we find evidence for significant effects of test round and sex. Fixed effects are given in brackets as estimates ± SE.

For fearlessness (TI), the interaction between sex and test round was (marginally) significant at the mature life stage in the inflated dataset (0.61 ± 0.30, t₁₃ = 2.01, p = 0.05), but not the reduced dataset (0.07 ± 0.21, t₁₁ = 0.33, p = 0.75). Conversely, at the young adult life stage, the interaction was significant in the reduced dataset (0.44 ± 0.21, t₁₃ = 2.15, p = 0.03), but not in the inflated dataset (0.48 ± 0.34, t₁₁ = 1.39, p = 0.17).

In young adults, aggression (AG) seemed to be lower in round 4 (-0.22 ± 0.11, t₁₁ = -2.05, p = 0.04), but this effect disappeared in the reduced dataset (-0.13 ± 0.10, t₁₁ = -1.27, p = 0.21).

Young adult males were more active than females (GA reduced dataset: 0.36 ± 0.15, t₁₂ = 2.36, p = 0.02; inflated dataset: 0.51 ± 0.27, t₁₁ = 1.84, p = 0.07).

Boldness (NO) was influenced by round, but this effect only became apparent in the reduced datasets, but not in the zero-inflated dataset. Thus, after removal of the birds that did not approach the novel object at all from the analysis, subadults were less bold in round 2 where a white object was used (estimate from unweighted model: -0.96 ± 0.21, t₁₁ = -4.54, p < 0.001). At both the young and mature adult stages, there was an increase in boldness between test rounds, when blue objects were used (estimates from unweighted models; young adult: 2.12 ± 0.26, t₁₂ = 8.16, p < 0.001; mature adult: 0.91 ± 0.41, t₁₂ = 2.24, p = 0.03). Furthermore, males were bolder than females at the young and mature adult life stages (young adult: 0.55 ± 0.24, t₁₂ = 2.29, p = 0.03; mature adult: 1.06 ± 0.39, t₁₂ = 2.72, p = 0.01). At the subadult stage, sex was only marginally significant (0.41 ± 0.22, t₁₁ = 1.85, p = 0.07). All these effects were not detected in the inflated datasets.

Exploration behaviour (NE) was influenced by sex at the mature life stage, with males being more explorative than females. This effect was significant both for the inflated dataset (0.84 ± 0.31, t₁₂ = 2.69, p = 0.01) and reduced dataset (0.49 ± 0.11, t₁₂ = 4.29, p < 0.001).

Effect of selection line origin

The birds chosen for this study originated from selection line experiments on animal personality. However, as they were from the F1-generation, we did not expect pronounced differences between lines yet, which is supported by the fact that line was not significant in most models except for NO at the subadult stage. Birds from the lines selected on low NE and high TI were less bold (estimates from unweighted model: NE low: -1.10 ± 0.52, t₁₁ = -2.13; TI high: -1.47 ± 0.66, t₁₁ = -2.22; p = 0.03). This effect of line did not occur in the GLMM using the reduced dataset without floor and ceiling effects.

Repeatabilities ('differential consistency')

We found evidence for considerable differences in repeatability of the five behavioural traits within and across life stages. We showed that in zebra finches, fearlessness (TI) and exploration (NE) were highly repeatable both within and across life stages. General activity (GA) and aggression (AG) were only repeatable at specific life stages and also across life stages, while boldness (NO) was not repeatable at all. Thus our results are within the range of repeatabilities described for behavioural traits in general [8, 9].

The factors possibly influencing repeatability include the type of behaviour, the interval between observations, age and sex of the individuals, which we will discuss in the following.

Behavioural traits differ in their repeatability

A categorization of behaviours into functional classes showed that exploratory behaviour and aggression were among the most repeatable traits, while activity was one of the least repeatable [8]. Our results partly confirm this as exploratory behaviour (NE) was highly repeatable throughout, and activity (GA) was only repeatable in subadults. In contrast to previous findings, we found aggression to be of only low repeatability across life stages, and non-repeatable within life stages. This might be explained by the fact that zebra finches are a highly gregarious species with comparatively low levels of aggression, hence there might be little need for consistent aggressive behaviour e.g. to assess or stabilize a dominance hierarchy through aggressive confrontations.

Tonic immobility (TI) has been described as the endpoint of a reflexive defence cascade in response to predators occurring in a wide range of animals and also in humans [46]. It represents an unlearned response reducing the likelihood of being attacked or increasing the survival after a predator attack [47] which seems to be evolutionarily conserved. Because TI is so closely linked to the underlying mechanisms [48, 49], this might be a reason why it is so highly repeatable. Also in previous analyses with a larger sample size [36], we have shown repeatabilities of adult zebra finches to be 0.44 for TI latency after 48 ± 17 days (CI 0.21-0.56, p < 0.001, N = 144).

NO was not repeatable at all, probably due to different amounts of excitement induced by the types of objects. Birds reacted particularly strong with avoidance to the blue coloured objects used in rounds 4 and 6. About half of the birds did not approach the blue objects at all, but the ones that did reacted even more boldly, approaching it quicker or more often than the light coloured objects. This notion is also supported by the subadult stage, when only wooden and white coloured objects were used, where we found the highest repeatability estimate for NO among the three life stages (R = 0.28, but CI still overlapping with zero). We should have used objects that elicit similar states of excitement and did not consider beforehand that a bright blue object would be too different from the other ones. Yet our aim in using different objects was to avoid habituation which might have occurred after using the same object or too similar types of objects repeatedly.

The time interval between tests influences repeatability

Longer intervals between tests are commonly expected to lead to decreased repeatabilities, for several reasons: firstly, because of environmental effects, as the animal is more likely to be in a similar state after a short interval. Secondly, the same is true for genetic influences (expression levels or epigenetic effects) on the phenotypic trait under study: after a longer interval the genetic changes might have been more pronounced. Thirdly, if the interval between tests even covers a major developmental re-organization such as sexual maturation (or a niche shift after metamorphosis e.g. [31]), it will probably also influence repeatability.

For instance, a study on male bush crickets (Sciarasaga quadrata) which were tested multiple times over their lifespan revealed significant repeatability in three parameters of calling activity throughout their lifespan, but a comparison of short-term and long-term intervals showed repeatability estimates decreasing over time [50].

Our study suggests differing results depending on the trait in question: in some cases, traits showed similar levels of short- and long-term repeatability (TI; AG in females) or were not repeatable at all (NO). In other cases, we found that traits were repeatable in the short-term (but also not in all life stages), but showed lower long-term repeatability (GA). Thus we found that repeatability estimates are not higher in general for short-term than long-term intervals, but it also has to be considered that our findings are based on a rather small sample size. The importance of taking multiple samples as suggested by Biro [51] becomes evident when comparing the repeatability estimates obtained for all six test rounds ('across lifestages') opposed to only rounds 1 and 6 ('long interval'). The power to detect significant repeatabilities if individuals are sampled over one year with only two measurements per individual seems to be too low in most cases (except TI).

Development of repeatability across life stages

With a test-retest interval representing about 15-20% of zebra finches' life span, we were able to show that TI and NE were repeatable throughout.

In AG and GA, the repeatability estimates across life stages (rounds 1 to 6) were significant, although within life stages this was not always the case. This might be due to the bigger sample size across life stages compared with the within life stage subsets, allowing for more precise estimation of repeatabilities. This point is also supported by the generally lower repeatability estimates for the 'long interval' using only test rounds 1 and 6. The resulting large confidence intervals overlapping with zero may be an artefact, because with a lower number of measurements per individual the confidence interval tends to be overestimated [9]. In a previous study [36], we found position diversity index, one behavioural variable of GA, to be repeatable (R = 0.44, CI 0.11 - 0.64, p = 0.013, N = 77) in adult zebra finches tested with an interval of 148 ± 67 days.

The notion that personality traits undergo a gradual development at early life stages and become stable only thereafter has been challenged recently [29] and indeed we found no evidence for this to be the case. Our results rather suggest the opposite, as we did not find increasing or stabilizing repeatability estimates, but rather that traits were either in the same range of repeatability estimates in all life stages, or as in GA, birds were only showing repeatable behaviour at the earliest life stage but not anymore later on.

Other studies also found differences in repeatability across life stages and depending on the trait. For instance in yellow-bellied marmots (Marmota flaviventris), boldness was repeatable in only one age class (yearlings) and docility was repeatable in all three age classes (juveniles, yearlings, adults) [17].

One argument which contradicts the expectation of life-long stability is that proximate mechanisms, such as brain structure development, finishing only in adolescence, lead to rather gradual adjustments taking place during sensitive periods throughout ontogeny [29].

Repeatability can differ between the sexes

Bell and co-workers [8] expected males to be more repeatable than females, which was confirmed in their overall meta-analysis. However, after exclusion of mate preference behaviours, females were the more repeatable sex.

We found that the sexes differed in repeatabilities, with males showing less repeatable behaviour than females. This resulted in overall repeatabilities being mainly driven by females' high estimates in the combined dataset, for instance in exploratory behaviour (NE). Females also became less repeatable in NE, while males became more repeatable (at the mature adult stage). Showing consistent exploratory behaviour might become important for zebra finch males at the adult stage after pair formation, as they are generally leading their females in feeding flocks [43].

In a study on domesticated female zebra finches [39], exploratory behaviour was repeatable over short-term (3 d and 7 d) and long-term (7 months) intervals. In contrast, another study [52] showed exploration was not repeatable in female zebra finches but in males. This suggests that either these differences result from the studies being undertaken with animals at different life stages or because there might be population differences in repeatability [53].

There is varying evidence regarding repeatability differences between the sexes in the literature. Supporting our finding that there was a tendency for lower repeatabilities in males, Burtka and Grindstaff [54] studied nest defence behaviour in Eastern bluebirds (Siala sialis) and found that females were more repeatable between years (within years also for one year, but not for the other). A study on field crickets (Gryllus integer) [55] reports that females were consistent in boldness (measured as latency to emerge from a refuge) tested across metamorphosis, while this was not the case for males. In common voles (Microtus arvalis), there were no sex differences in repeatabilities of boldness, exploration or activity [25].

The expectation to find sex differences in repeatability probably depends on the type of behaviour in general, but also on the species and their biology. Especially if traits are under sexual selection and related to mate choice, it might be advantageous to be predictable [56]. For instance, Nakagawa et al. [57] found that male house sparrows (Passer domesticus) were more repeatable than females (within and between years) in parental care behaviour (feeding rate).

Development of behavioural syndrome structure over ontogeny

Three of the investigated traits formed a behavioural syndrome, but only in young adults. Thus, we found that the behavioural syndrome was not stable across life stages but specific for a certain age. Furthermore, it consisted only of non-repeatable traits and was mostly driven by the female sex.

At the young adult stage, the zebra finches were about 100 d old which is the time they are sexually mature and able to find a partner and reproduce, given favourable environmental conditions. This suggests that during the time of sexual maturation, traits are being restructured, so that previously independent traits become shortly linked together in young females when they reach reproductive age, but correlations then disappear again in mature animals. Young females that were more active also interacted more aggressively with the mirror, but were less bold towards a novel object. Possibly during the time when females are usually choosing a partner they become more active and more aggressive against female competitors. As the traits stayed uncoupled in males, there might be different selection regimes on the sexes.

It is expected that selection pressures on juveniles and adults differ significantly [27], leading to differences in personality axes as well, which in turn may influence the evolution of personality (e.g. if selection mainly acts on the juvenile stage). Differences in syndrome structure over ontogeny have been shown to occur in a number of study species. Bell and Stamps [15], in their study on two stickleback populations, found a stable boldness-aggression syndrome in one population, while the single behaviours constituting the syndrome were not repeatable. They also found (in the same population) correlations that were present during the juvenile stage, then disappeared in subadults but re-emerged in adults. In cavies, boldness and exploration were correlated in mature animals, but not in juveniles [24]. In wild brown trout, a behavioural syndrome including activity, aggression and exploration developed after an interval of two months covering a survival bottleneck [19]. In contrast, a stable syndrome structure comprising activity, exploration and boldness was detected in firebugs measured across life stage transitions [23]. In a meta-analysis on behavioural syndromes [20] the mean effect size was rather small (r = 0.264; 95% CI 0.210-0.316).

Thus it seems that behavioural syndromes are long-term consistent only in few cases, but can also quite readily change in others, with the correlations being broken or new ones forming over development. Changes in correlation structure can be expected to occur during development for several reasons, for instance because during major life stage transitions the most advantageous suite of behaviours changes. This is especially the case if the environment changes between age classes or if there is a prominent niche shift, as in the most extreme case of metamorphosis [1]. Another possibility is that when traits are affected by the same hormones, a shift in hormonal levels or a general hormonal reorganization such as during sexual maturation may influence the links between behavioural traits in such a way that they influence syndrome structure leading to instability [15]. On the other hand, it has been proposed that behavioural syndromes will be particularly stable if they are caused by genetic correlations, such as pleiotropy, or correlational selection [16, 33]. However, the mechanisms generating, maintaining or disrupting behavioural syndromes are far from well understood and deserve further investigations.

We conclude that the stability of behavioural traits as well as their correlations cannot be assumed but need to be tested. This leads to the prediction that if single traits are regulated (more or less) independently from one another, animal personality might be much more flexible than expected so far, which is also supported by endocrinological findings [29]. Of course we are the first to admit that care needs to be taken as the low sample size, especially per sex, is a limitation in our study, meaning that we possibly could not detect further correlations due to this restriction [20]. Nevertheless we believe that more attention is needed when animal personality and behavioural syndromes are simply assumed to exist, rather than investigated.

Animals and housing conditions

The zebra finch (Taeniopygia guttata castanotis) is a well-studied model organism [60] which can be bred and kept in captivity readily and has a short generation time, as juveniles are independent after 35 days and sexually mature after 65 days.

We conducted experiments on animal personality in zebra finches in a project on divergent bi-directional selection lines. The monitoring of the birds included five different personality tests, three of which served as selection parameters for the differing selection lines. The long-term design of the selection line study allowed us to assess repeatability and the development of behavioural syndromes over a longer time scale than is possible in most studies on long-lived organisms.

We chose 52 wild-type zebra finches from the F1-generation, 22 males and 30 females, to study long-term consistency of personality traits. They were bred from 31 different pairs, with between 5 and 14 individuals coming from the same selection line and maximum 3 full-sibs per family. They showed no behavioural trait values at the extreme end of the personality axis in the direction their parents were selected on.

The zebra finches were reared in outdoor aviaries (6 × 2 × 3 m) at Bielefeld University in the different selection line groups, each consisting of five breeding pairs plus their offspring.

After independence from their parents at ca. 40 days of age (median 43 d, min 30 d, max 62 d, IQR 8 d), subadult birds were translocated from their natal aviary groups into indoor mixed-sex tutor groups (cages: 81 × 60 × 50 cm). Tutor groups consisted of seven to ten juveniles with a pair of unrelated, unfamiliar adults. Before sexual maturation, at about 60 days of age (median 60 d, min 47 d, max 81 d, IQR 9 d), birds were then transferred to double cages (82 × 40 × 30 cm) in groups of three to four same-sex individuals. Tutor groups and same-sex cage groups were both arranged so they comprised animals originating from several different selection lines.

All animals were maintained in the same room with auditory and visual contact between cages. At all times, birds had ad libitum access to commercial zebra finch seed mix (Elles, Mischfutter für Exoten, L. Stroetmann Saat, 48163 Münster, Germany) and fresh water. Additionally, a mixture of germinated seeds and egg food (Cédé N.V., 9940 Evergem, Belgium) was provided daily in outdoor aviaries and tutor groups and thrice weekly after the transfer to double cages. The diet was weekly supplemented with fresh greens. In the housing room, birds were kept on a 14:10 light-dark-cycle, additional to the natural light conditions.

Test schedule and description of behavioural tests

Behavioural tests started while birds were still kept in the tutor groups. All individuals were tested in five behavioural tests

(a) twice at the 'subadult' stage

(b) twice shortly after reaching adulthood ('young adult')

(c) twice at approximately one year of age ('mature adult')

leading to a total of six test rounds, with approximately 14 days apart within a given life stage. The interval between life stages was about one month for the transition to adulthood, and seven to ten months for the transition between young adult and mature adult. For all rounds, mean age and intervals between tests are given in table 1.

Experiments were conducted always in the same order: (TI - AG - GA - NO - NE) between 8:30 and 20:00 hours, and we took care that between tests, birds had at least one day break.

Tonic immobility (TI): An empty wire cage (40 × 31.5 × 21 cm) located on a table in a sound proof chamber (to exclude noise disturbances during the test) was used to conduct the tonic immobility test. The bird was placed on its back with wings pressed to its body on a metal holding cradle padded with foam rubber. An experimenter locked the bird in this position for five seconds by gently pressing index and middle fingers on its chest and then retracting the hands. Movements by the observer were minimized to prevent disturbances. Birds were considered successfully immobilized if they stayed in this position for at least five seconds. If a bird was not immobilized at once the procedure was repeated up to a maximum of 10 inductions. The test was terminated after a maximum of 20 minutes immobility (ceiling value of latency = 1200 s). Tests were performed by YW and student assistants Simon Tiersch and Nele Heitland (rounds 1-4) and Ivonne Kienast (rounds 5 and 6). Tonic immobility latencies and number of inductions are considered a measure of fearfulness [61] or boldness in predator contexts [62].

Aggression (AG): The aggression test cage (after [42]) was a small wire cage (40 × 31.5 × 21 cm) with one perch (39 cm). At the far end of the perch, a mirror (15 × 15 cm) was fixed and a food dispenser placed underneath.

The focal individual was transferred to the test cage with the mirror being covered with a piece of cardboard. After an habituation phase of five minutes, the mirror was uncovered and the recording started for five minutes. Birds were tested with the mirror fixed on either the right (= R) or left (= L) side of the cage in a given sequence (R-L-L-R-R-L). During the test, another cage of birds was present in the test room for auditory contact (always the same birds).

The first five seconds after removal of the division were discarded. Then, the frequency of the following interactions was determined within five minutes of video recording:

- Pecking/beak contact: each contact of the bill with the mirror.

- Flying against the mirror: body contact with the mirror while feet are not touching the perch.

- Breast contact: touching mirror with breast, usually after straightening up to full body height.

- Head contact: Usually occurring when sitting quite close to the mirror, it is touched with the head while the feet still remain on the perch.

We used the sum of aggressive interactions with the mirror as measure for aggressive behaviour.

General activity (GA): The general activity test was conducted in the birds' home cage. All birds were removed from the cage prior to testing and relocated to a waiting cage (without visual contact to the focal individual). The food dispenser remained in the home cage during testing. After replacing the focal individual to its home cage and a two minute recovery time, the activities of the bird were video-recorded for ten minutes.

In the video analysis, the first 20 s were discarded (settling period). Then, the number of flights in the cage was determined. A flight was defined as a movement between perches or between perch and cage floor. We also determined the events of visits to all seven possible positions in the cage (four perches, left, middle and right bottom part of the cage) and calculated a position diversity index (PDI, equivalent to the calculation of the Exploration Diversity Index (EDI) as described in [53]). A higher position diversity index means that an animal visited more different places in the cage and thus indicates a higher activity.

Novel object (NO): Directly following the activity test, the divider was placed in the cage and a novel object inserted in the empty cage compartment. The position of the novel object (far left or far right perch) was randomly assigned beforehand. The novel object was fixed to a perch and placed inside the cage, while the focal individual was in the other compartment. The divider was then removed after the recording started, and the experimenter left the room. Behaviour was video-recorded for ten minutes.

This procedure was repeated with the remaining cage mates, while the already tested birds were transferred to the waiting cage.

As novel objects we used wooden blocks, styrofoam balls and chipboard discs of different colours. Details regarding the type of novel object used in each test round are given in part 6 in additional file 1.

Data of the activity test as well as the novel object data were analysed using CowLog [63]. The events and durations of visits to each of the seven positions in the cage were determined, as well as the latency to land on the novel object perch. We used the novel object latency, the frequency to visit the NO perch and the percentage of time spent there as variables for boldness. GA and NO tests were both conducted in the animal housing room.

Novel environment (NE): As a novel environment, we used standard double cages whose interiors were covered with adhesive foil. Otherwise, they were equipped as the usual home cages (sand on floor, four wooden perches, food dish with seed mixture, water bowl). Descriptions of the cage interiors used as novel environments are given in part 6 in additional file 1. Up to four individuals were tested simultaneously (in separate cages); in general all cage mates at a time. The focal individuals were transferred from their home cage into separated deprivation cages (30 × 40 × 40 cm) without food or water about 3.5 h (207 ± 10 min) prior to testing. Each focal individual was then put into a start box (11 × 17 × 11 cm) which was attached on either the left or right side of the experimental cage. A cardboard division between cages ensured that there was no visual contact between focal individuals in the start boxes. In the cages a food dish was placed on the same side as the start box, a water tray in the other compartment. The division between start box and novel environment was lifted and the behaviour video-recorded for one hour. A cage with flock mates was present in the experimental room during testing for auditory contact with the focal individuals.

As explorative behaviour, we determined how many of the seven possible positions in the cage were visited within one hour and we measured the latency to enter the novel environment cage and the latency to visit all seven positions if applicable. If the start box was not left or not all positions were explored, birds received a maximum value of 3601 seconds.

Statistical analysis

All statistical analyses were run in statistical software R [64] within the RStudio environment [65], using the packages lme4, rptR, lattice, psych and psychometric [66–70].

Principal components analysis

For each round per test we ran a principal component analysis (PCA) calculated with correlation matrix to extract the first axis for each personality trait. PCA were performed using the maximum available dataset (max. N = 52, for sample sizes per test and round see part 2, additional file 1). The included behavioural variables and resulting PC loadings are given in part 2 in additional file 1. Because seven data points from six individuals were missing, the final dataset for subsequent analyses was reduced to 46 individuals (n = 276). For general activity (GA) round 1, the derived scores were inverted by multiplication with (-1), so that higher scores represent bolder behaviours in all tests and rounds (GA = more active; NO = bolder towards object; NE = faster exploration; TI = lower latency and more inductions; AG = more aggressive). As the loadings differed between test rounds for aggression (AG) and the derived PC scores were not normally distributed, we instead used the decimal logarithm of the sum of aggressive interactions with the mirror as a score for later analysis, which showed a distribution closer to normality.

Generalized Linear Mixed Models for extraction of best linear unbiased predictors

For each personality axis in each life stage, we ran a univariate Generalized Linear Mixed Model (GLMM) with Gaussian error structure fit with Maximum Likelihood. This was done to extract the best linear unbiased predictors (BLUPs) for behavioural syndrome structure calculation and to check for effects of fixed factors such as sex differences in mean level behaviour.

Models were run in R using the function lmer from the R package lme4 [66]. In the full model we included the three fixed effects sex, test round and selection line and the two-way interaction of sex and test round as well as the random effects batch (to control for possible season effects, as all birds from the F1-generation were transferred from outdoor aviaries to indoor cages in a total of 10 batches between mid-June and mid-December 2012) and individual ID nested within mother ID. We conducted a stepwise deletion of fixed effects if model comparison using a likelihood-ratio analysis between the full and reduced model was non-significant. The random effects as well as test round and selection line as a fixed effect to control for were always retained in the final model.

Because Fligner-Kileen tests indicated heterogeneity of variance in NO scores for the interaction of the fixed factors sex and round, we weighted this effect (by the ratio between the variances of the factor levels) to improve the fit of the model. NO models were calculated with and without weights and the model results for both are presented in parts 4 and 5 in additional file 1. Although likelihood-ratio tests showed that the weights significantly improved the fit, we decided to use the unweighted models for subsequent analyses after visual inspection of the residuals plotted against the fitted values.

Because of ceiling or floor effects of the data distribution, we created subsets of data containing only valid data points (without zero- or maximum-inflated data) as follows: TI was maximum-inflated as individuals that did not enter the immobile state after ten inductions received a latency of one second (leading to high PC scores). The reduced dataset included data points with latencies larger than one. AG, NO, NE and GA all were zero-inflated, as some individuals did not interact with the mirror (interactions = 0), not approach the novel object (events on perch = 0), not explore all seven positions the novel environment (received a maximum latency of 3601 s resulting in lowest PC scores) or not move in the home cage (flights = 0) at all. In the respective reduced datasets these data points were removed. Each model (one per life stage and personality trait) was run for the inflated dataset and reduced dataset separately and the model outputs are given in parts 4 and 5 in additional file 1. Although results of the reduced dataset deviated from results of the inflated dataset in some cases, we continued calculations with the complete (inflated) dataset, as using the reduced dataset would have decreased our overall sample size (see part 5 additional file 1 for reduced dataset sample sizes) immensely. We discuss these cases carefully.

Repeatabilities ('differential consistency')

We used the R package rptR [67] to calculate repeatabilities within and between life stages. All response variables (PC axes and logarithmized sum of aggressive interactions) were approximately normally distributed and thus analysed with Gaussian fit and REML estimation. We report repeatability estimates with their significance levels and 95% confidence intervals, calculated with 1000 permutations and 1000 bootstrappings. We calculated repeatabilities within each life stage (comprising two subsequent test rounds), for the complete dataset including all six test rounds (referred to as 'across life stages') and also including only test rounds 1 and 6, to get an estimate across life stages using only two data points (referred to as 'long interval'). We compared the overlap of confidence intervals, showing non-significance at the 5% threshold if applicable [71]. As the mixed models had indicated a significant effect of sex in some cases, we calculated repeatabilities for males and females separately as well.

Behavioural syndrome structure ('structural consistency')

We extracted the best linear unbiased predictors (BLUPs) of the GLMMs. BLUPs represent scores that are "controlling for" the effects included in the model. Behavioural syndrome structure at the three life stages was assessed by computing Spearman rank correlation matrices with Holm's correction using the life stage mean of the fitted values of the minimum adequate models of the five personality traits for each life stage (R package psych [69]). The lower and upper bounds of 95% confidence intervals were estimated with the R package psychometric [70]. We also calculated behavioural syndrome structure separately for males and females.