Seasonal change in sexual differences in nestling size and survival: a framework to evaluate sex-dependent environmental sensitivity in the wild

Abstract

Background

Sexual size dimorphism (SSD) occurs in a wide range of species in birds and other animals, but the magnitude of SSD often varies with environmental conditions. In general, in the developmental stages, the larger sex is more vulnerable to adverse environmental conditions because the larger sex requires more energy than the smaller sex. However, this may not hold true for birds with large brood sizes; the larger sex can acquire more food by suppressing the smaller sex. In addition, most previous studies have been experimental, such as by manipulating clutch size and ectoparasites, which may not reflect natural conditions.

Methods

In the present study, we propose a general framework to assess sexual differences in environmental sensitivity in natural populations. Because environmental conditions change throughout the breeding season, seasonal changes of nestling SSD and sex ratio should reflect sexual differences in environmental sensitivity. We applied this approach to a large dataset (1555 nestlings over 5 years) of Japanese Tits (Parus minor). In this population, the male nestling is generally larger than the female (5% SSD in body weight).

Results

We found that the magnitude of SSD (weight, tarsus, wing) and fledgling sex ratio increased both in the beginning and the end of the breeding season.

Conclusion

Our study suggested that female nestlings are more valuable to poor environmental conditions in the relatively fecund species. This study underscores the importance of brood size on sexual differences in environmental stochasticity and our framework encourages comparative analysis among different bird species.

Background

Environmental conditions during early developmental stages can have long-term effects on individual life history (Lindström 1999; Monaghan 2008). For sexually dimorphic species, neonatal conditions are especially important because either one of the sexes could be more significantly influenced by environmental stochasticity. For example, female nestlings in raptors and male juveniles in mammals, the larger sexes, tend to be strongly influenced by poor environmental conditions due to their higher energy requirement (costly sex hypothesis, Clutton-Brock et al. 1985; Teather and Weatherhead 1989; Kalmbach and Benito 2007; Jones et al. 2009). Conversely, the larger sex will experience greater fitness benefits under favorable conditions (Jones et al. 2009). Therefore, sexual differences in environmental sensitivity can significantly affect individual fitness and even allocation of offspring sex ratio (Trivers and Willard 1973).

The costly sex hypothesis has been widely supported for species with small litter/brood sizes and large sexual size dimorphism (SSD) (Clutton-Brock et al. 1985; Teather and Weatherhead 1989; Kalmbach and Benito 2007; Jones et al. 2009). However, some recent studies of birds with large brood sizes and small SSD suggested that the smaller sex is more sensitive to poor conditions (Oddie 2000; Potti et al. 2002; Råberg et al. 2005; Dubiec et al. 2006; Nicolaus et al. 2009), although other studies have shown the opposite pattern (Tschirren et al. 2003; Chin et al. 2005; Dietrich-Bischoff et al. 2008). In the case of strong sibling competition for food delivered by parents, a larger body size may be advantageous (competitive advantage hypothesis, Nilsson and Gårdmark 2001; Nicolaus et al. 2009). Therefore, the mixed results among previous studies necessitate further research on species with large brood sizes. Moreover, most previous studies on birds with a large brood size and small SSD have been conducted experimentally, such as by adding/removing nestling to clutches or ectoparasites (Råberg et al. 2005; Kalmbach and Benito 2007), which may not reflect natural conditions. Thus, more observational studies are needed to confirm patterns under natural conditions.

We propose a general framework to assess the sexual differences in environmental sensitivity under natural conditions by monitoring seasonal change in the magnitude of SSD and fledging sex ratio (Fig. 1). First, seasonal change in nestling growth should be curvilinear (peak in mid-season) rather than a simple increasing or decreasing trend because in general, some pairs in a population start breeding early or late from the peak of food availability (i.e. abundance and/or quality of prey item, Naef-Daenzer and Keller 1999; Verboven et al. 2001; García-Navas and Sanz 2011). Therefore, if male and female nestlings respond differently to the environment, seasonal change in SSD should reflect this environmental vulnerability (Fig. 1a, b). When males are larger, and if males are more sensitive, the magnitude of SSD would be smaller both in the beginning and end of the breeding season (Fig. 1a). Conversely, if females are more sensitive, the magnitude of SSD would be larger both at the beginning and the end of the breeding season (Fig. 1b). A similar prediction applies to sexual differences in nestling mortality. If males are more sensitive, the proportion of males at fledging would be lower both at the beginning and the end of the breeding season (Fig. 1c). Alternatively, if females are more sensitive, the proportion of males would be higher both at the beginning and the end of the breeding season (Fig. 1d). In this study, we applied this approach to a large data set (1555 nestlings over 5 years) of Japanese Tits (Parus minor) in a temperate forest in northern Japan.

Fig. 1

Conceptual models of seasonal change of male and female nestling body size and fledgling sex ratio when males are more sensitive (a, c) and females are more sensitive (b, d). Dashed arrows show the difference between the male and female body size

Methods

Study species, study sites, and field procedures

The Japanese Tit is a small (ca. 14 g) hole-nesting passerine that readily accepts nest boxes for breeding. This species is sexually dimorphic and sex differences appear in the nestling period (5% of male biased SSD in body weight; Nomi et al. 2015). In our study area, the Japanese Tit has a large brood size in the first and later clutches. However, clutch size and fledging success were lower in later clutches (10.2 for the first, range 7–13, and 8.2 for the second clutches, range 5–10; Nomi et al. 2015), most likely because food quality and abundance were unfavorable later in the season (Murakami 2002).

Fieldwork was conducted over 5 years (2009–2010, 2012–2014) in the Tomakomai Experimental Forest (42.40°N, 141.36°E, 5–90 m a.s.l.), Hokkaido, northern Japan. The study area belongs to the cool temperate climate zone and is covered with mixed deciduous forests. We established two study sites in the forest (ca. 30 ha each) in 2008 with approximately 150 nest boxes at each site (Yuta and Koizumi 2012; Nomi et al. 2015). During the breeding season (late April to late August), we checked all the nest boxes weekly and recorded basic information, such as laying date and clutch size. Hatching date and fledging date of broods were confirmed by daily observation (date of one or more nestlings hatched or fledged). Parents were caught and ringed for individual identification when nestlings were 5–7 days old and measured at the same time. Body weight was measured using an electronic balance (accuracy: 0.1 g). Tarsus length was measured using a digital caliper to the nearest 0.01 mm. Wing length was measured using a ruler to the nearest 0.5 mm. Nestlings were caught and measured on day 10–15 post-hatching. After measuring, blood samples (10–50 μL) from the brachial vein were collected and stored in 99% ethanol for molecular sexing.

Sex determination

Sex determination was conducted following the protocol in Nomi et al. (2015). After extracting DNA from samples, we used the primers P8 and P2 to amplify sex specific regions of the CHD gene, which is often used for sex determination in birds (Griffiths et al. 1998). Birds were sexed according to the presence of the PCR products of CHD-Z (both sexes) and CHD-W (females only), which were separated by electrophoresis. We confirmed the accuracy of molecular sexing using the parents’ samples with known sexes.

Statistical analyses

Seasonal change in the magnitude of SSD

For the analysis of seasonal change in the magnitude of SSD, we used linear mixed models (LMM). The response variable was body weight, tarsus length, or wing length of nestlings. The explanatory variables in the full model included hatching date, square of hatching date, sex, interaction between the hatching date and sex, and interaction between the square of hatching date and sex as main effects. We included nest ID, study site, and year as random effects. Number of nestlings (brood size when nestlings were caught), nestling age, and maternal traits (weight, tarsus length, and wing length in each analysis) were also included as covariates. We did not include paternal traits because some nestlings were sired by extra-pair males in this population and their genetic fathers were not identified (Yuta and Koizumi 2016). Since the mean timing and length of the breeding season varied among years, we standardized the hatching date of each nest by compressing the hatching date to a range from 0 to 1 (standardized hatching date = (hatching date − earliest hatching date)/hatching date range). Standardized hatching date (hereafter termed hatching date) and sex were centered by subtracting the mean prior to the analyses to avoid multicollinearity between the variables and interactions (Robinson and Schumacker 2009).

To examine the seasonal change of the magnitude of SSD, seasonal change of male and female nestling body traits estimated by the best models were assigned to the Storer’s index of the magnitude of SSD (difference between the values of male and female body traits divided by the mean values of male and female traits, Benito and Gonzalez-Solis 2007).

Seasonal change in fledgling sex ratio

We defined the fledgling sex ratio as the proportion of male nestlings on the capturing dates (10–15 days old) because almost all nestlings successfully fledged if they were alive on the capturing dates. However, if dead nestlings were found after the capturing date, they were excluded from the analysis. We constructed generalized linear mixed models (GLMM) with a logit-link function with a binomial error distribution for the analysis. The response variable was nestling sex (male: 1, female: 0), and the explanatory variables were hatching date, square of hatching date, and brood size. Random effects included nest ID, female ID, site, and year. We excluded depredated nests from the analysis because nest predation potentially biases the fledgling sex ratio. To confirm that the fledgling sex ratio was independent of sex allocation at the time of laying (primary sex ratio), we conducted the same analysis using the data of nests in which the number of eggs and fledglings were equal.

We performed all the analyses using the statistical software R 2.15.3 (R Development Core Team 2013) with package “lme4” for LMMs and GLMMs. Multicollinearity of explanatory variables was tested using the variance inflation factor (VIF) with the package “car”. The VIF of all variables was under 3.5 after centering and were under the threshold value of 5.0 recommended by O’Brien (2007). To confirm whether curvilinear models were better than simple linear models, and whether the models with an interaction between hatching date and sex were better than the models without interaction, we compared all candidate models (i.e. different combinations of the explanatory variables) using AIC (Burnham and Anderson 2002). The model with the lowest AIC value in each analysis was determined to be the best model. The significance of the partial coefficient of each variable in the best model was examined by the likelihood ratio test comparing the models with and without the terms of interest.

Results

Seasonal change in the magnitude of SSD

In the 5 years of this study, we collected data on 1555 nestlings from 194 broods. We examined seasonal change in the magnitude of SSD in three nestling body traits (body weight, tarsus length, and wing length). We compared models with linear hatching date effects and the quadratic term of the hatching date, as well as the models with and without their interactions with sex, using AIC. In all three analyses, the models with the square of hatching date and the interaction of the square of hatching date and sex had the lowest AIC values (Table 1, Additional file 1: Table S1, Additional file 2: Table S2, Additional file 3: Table S3), meaning that the sex specific curvilinear model was the best model for predicting seasonal change in nestling body size. For body weight and tarsus length, the square of hatching date, sex, and the interaction between the square of hatching date and sex were all significant (Table 2). For wing length, although the square of hatching date was not significant, interaction between the square of hatching date and sex was significantly correlated, indicating that seasonal change in growth of wing length was sex specific. In all nestling body traits, the magnitude of SSD was higher both at the beginning and the end of the breeding season (Fig. 2).

Table 1 Model ranking of the GLMMs according to the AIC

Table 2 Parameter estimates and p values of the best models of GLMMs

Fig. 2

The relationship between hatching date and nestling body traits (left figures) and the magnitude of SSD (right figures). Parameters were estimated using the linear mixed model shown in Table 2. Lines are predicted values (male: blue, female: red). Filled circles represent males and open circles represent females. Seasonal changes of the magnitude of SSD were calculated using the best models shown in Table 2, where predicted lines of male and female body size were assigned to the Storer’s index of the magnitude of SSD (difference between the values of male and female body size divided by the mean values of male and female size, Benito and Gonzalez-Solis 2007)

Seasonal change in fledgling sex ratio

We also examined seasonal change in fledgling sex ratio to elucidate which sex was more vulnerable in poor environmental conditions. Comparisons of models by AIC showed that the curvilinear model was better than the single linear model (Table 1, Additional file 4: Table S4). In the analysis, the square of hatching date was significant (Table 2). The proportion of males was higher both in the beginning and end of the breeding season (Fig. 3). To examine whether the seasonal change in proportion of males was a result of primary sex allocation (sex ratio at eggs), we conducted the same analysis using nests with no loss of nestlings (the number of fledglings was equal to the number of eggs). Primary sex ratio (proportion of males in eggs) did not have a significant curvilinear relationship with hatching date (n = 936, partial coefficient of the square of hatching date; χ² = 2.82, p = 0.09), indicating that the seasonal change in fledgling sex ratio was not the result of primary sex allocation of female parents, but rather the result of female-biased mortality.

Fig. 3

The relationship between hatching date and fledgling sex ratio (proportion of males). Plots and bars represent mean values and standard errors, which are calculated by dividing the data into 5 rank categories starting from the lowest value of standardized hatching date (first 4 plots, 298 samples; last plot, 299 samples). Parameters were estimated using the generalized linear mixed model shown in Table 2. The line is the predicted value

Discussion

Our results consistently suggested that female Japanese Tit nestlings are more vulnerable to adverse environments. This agrees with most experimental studies on species with large brood sizes (Oddie 2000; Potti et al. 2002; Råberg et al. 2005; Dubiec et al. 2006; Nicolaus et al. 2009, but see Tschirren et al. 2003; Chin et al. 2005) and a few observational studies (Dhondt 1970; Eeva et al. 2012, but see Dietrich-Bischoff et al. 2008). Our approach was observational and, therefore, a clear demonstration of the causal link was not possible. In fact, the increase of SSD in the end of the breeding season can also be explained by factors other than poor environmental conditions. For example, multiple brooding is energy-consuming for parents, which could result in poorer parental body conditions, and consequently food delivery rate might have reduced. However, the main question in this topic is which sex is more valuable when the parental provisioning is not sufficient: for nestlings, parental condition and behavior can also be considered as environmental factors. Poor external environmental conditions and poor parental conditions should similarly affect the nestlings. In addition, past experimental studies have not controlled the external environment either; these studies examined the sexual differences of environmental stochasticity by changing the food intake per nestling by adding or removing clutches (Råberg et al. 2005; Dubiec et al. 2006; Nicolaus et al. 2009). In this sense, we need not change the main conclusion that female nestlings are more valuable under poor environmental condition.

On the other hand, observational study has some advantages. For example, while many experimental studies were conducted under extreme conditions, our approach and results reflect natural settings. Moreover, seasonal changes in nestling or fledgling characteristics are commonly monitored; thus, this approach can easily apply to many study systems, which in turn encourages comparative studies or meta-analysis. Based on past studies, it seems that male nestlings are more sensitive in species with small brood sizes, whereas the opposite is true in species with large brood sizes (Jones et al. 2009).

As suggested by Råberg et al. (2005), parents of large broods might not able to control distributing food and the larger sex has an advantage in sibling competition. For example, in Great Tits, the feeding position of individual parents is consistent through the breeding period (Poelman et al. 2006). Therefore, the larger nestlings could monopolize food delivered by the parents by moving toward the best position. Moreover, most species with large brood size, like Japanese Tits or Great Tits, are cavity nesters and therefore, it may be difficult for parents to distinguish between hungry and full nestlings in dark environment.

Alternatively, it is also possible that parents favored larger nestlings in poor condition and delivered more food to the larger nestlings (Jones et al. 2009; Caro et al. 2016). To disentangle whether it is a result of parental favoritism to larger nestling or of sibling competition, more observations are needed inside the nest boxes comparing good and poor environmental condition.

In addition to this study’s main results, we found patterns in the magnitude of SSD. First, the SSD of nestling body traits were larger at the end of the breeding season than in the beginning (Fig. 2). This may be because food abundance in this study site is lower in the end of the breeding season than in the beginning. Although caterpillar abundance and diversity in the studied forest had two peaks in the season (Yoshida 1985), clutch size and fledging success was lower in the later clutches (Yuta and Koizumi 2012), indicating that environmental condition including food availability were less favorable towards the end of the breeding season. Moreover, asynchronous hatching timing in a brood might also have increased sex difference later in the season. Hatching asynchrony, which is considered as a maternal adjustment under poor food availability, becomes more prominent later in the season (Theofanellis et al. 2008) and more severely affects the growth of late-hatched female nestling compared to male siblings (Oddie 2000). Some nestlings might have hatched 1–3 days later than the first hatched nestling in a brood, but we did not confirm hatching dates of all nestlings. Thus, larger SSD later in the breeding season may be the result of increased hatching asynchrony as well as decreased food availability.

Second, we found that the SSD of nestlings was smaller for tarsus length and wing length compared with body weight. This may be the result of differences in growth speed among body traits. For example, the tarsus grows faster than body weight or wing length and reaches a maximum around 10 days old (Royama 1966). In terms of wing length, sex differences in growth strategies have been reported in many studies (Oddie 2000; Råberg et al. 2005; Dubiec et al. 2006; Nicolaus et al. 2009). Since the female is the more dispersing sex, they may invest more in wing length than in body weight under poor conditions (Greenwood 1980). Therefore, while significant differences were detected in the Japanese Tit, the sexual difference was subtle in wing length. In contrast, males may allocate more into body weight because body weight plays an important role in competition among individuals and, thus, in the process of territory acquisition (Garnett 1981; Sandell and Smith 1991).

As discussed in some recent studies, differential effects of adverse condition on male and female traits may be a result of allocation of resources rather than the difference in sensitivity between the sexes (Tschirren et al. 2003). Similar to wing length, selective pressure may work differently on body traits between the sexes and future studies should focus more on this point.

Conclusion

Our study showed that female Japanese Tit nestlings are more valuable to poor environmental conditions, which is consistent with most experimental studies on species with large brood sizes. However, whether it is a result of parental favoritism to larger nestling or of sibling competition is still unclear. Our study underscores the importance of brood size on sexual differences in environmental stochasticity and our framework encourages comparative analysis among different bird species.