Glucocorticoids, state-dependent reproductive investment and success in the face of danger in a long-lived bird

Glucocorticoid hormones may mediate trade-offs between current and future reproduction. However, understanding their role is complicated by predation risk, which simultaneously affects the value of the current reproductive investment and elevates glucocorticoid levels. Here, we shed light on these issues in long-lived female Eiders (Somateria mollissima) by investigating how current reproductive investment (clutch size) and hatching success relate to faecal glucocorticoid metabolite [fGCM] level and residual reproductive value (minimum years of breeding experience, body condition, relative telomere length) under spatially variable predation risk. Our results showed a positive relationship between colony-specific predation risk and mean colony-specific fGCM levels. Clutch size and female fGCM were negatively correlated only under high nest predation and in females in good body condition, previously shown to have a longer life expectancy. We also found that younger females with longer telomeres had smaller clutches. The drop in hatching success with increasing fGCM levels was least pronounced under high nest predation risk, suggesting that elevated fGCM levels may allow females to ensure some reproductive success under such conditions. Hatching success was positively associated with female body condition, with relative telomere length, particularly in younger females, and with female minimum age, particularly under low predation risk, showing the utility of these metrics as indicators of individual quality. In line with a trade-off between current and future reproduction, our results show that high potential for future breeding prospects and increased predation risk shift the balance toward investment in future reproduction, with glucocorticoids playing a role in the resolution of this trade-off.


Introduction
The 'cost of reproduction' hypothesis (Williams 1966) assumes a trade-off between current reproductive investment (CRI) and residual reproductive value (RRV) (Stearns 1992;Clark 1994). In iteroparous animals, RRV is expected to be higher while CRI is expected to be lower in young individuals, as future reproductive opportunities decrease late in life (Velando et al. 2006). However, the relative value of CRI and RRV is also affected by extrinsic factors such as predation which shapes current reproductive success (Ghalambor and Martin 2000). Reduced reproductive effort under high predation risk may save resources for reproduction during safer periods and thus increase lifetime reproductive success. For example, in an elaborate experimental study eliminating direct predation but using playbacks of predator calls, perceived predation risk alone significantly reduced both clutch size and chick survival in songbirds .
Considerable interest has been given to glucocorticoid hormones as candidate physiological mediators of the tradeoff between CRI and RRV. For example, bird species with a higher value of the current brood relative to RRV mount weaker corticosterone responses to stressors, suggesting a putative role of glucocorticoids in mediating the cost of reproduction (Bokony et al. 2009). Besides their important role in regulating energy expenditure (Jimeno et al. 2018), glucocorticoids are also secreted upon exposure to the presence or cues of predators (Cockrem and Silverin 2002;Ibáñez-Álamo et al. 2011). The resulting elevated glucocorticoid levels should favour self-maintenance and survival processes at the expense of less immediately vital functions such as reproduction (Wingfield et al. 1998). Thus, recurring predator stimuli may lead to chronically elevated glucocorticoid levels (Monaghan 2014) and result in lower breeding propensity (Öst et al. 2018), increased nest desertion (Thierry et al. 2013) and, ultimately, reduced fitness (Bonier et al. 2009b). Nevertheless, the relationship between baseline glucocorticoid levels and fitness is closely dependent on context (Bonier et al. 2009a;Jaatinen et al. 2013). For example, moderate increases in glucocorticoids may positively influence parental performance, by enhancing foraging effort before and during egg-laying (Crossin et al. 2012). Furthermore, glucocorticoids are simultaneously affected by predation risk and involved in the regulation of the trade-off between CRI and RRV. This dual role complicates interpretations on the role of glucocorticoids in mediating reproductive trade-offs under variable risk of predation.
Estimating an organism's potential for future reproduction can be difficult, complicating efforts to resolve the trade-off between current and future reproduction. Recent advancements in molecular ecology suggest that telomeres may serve as a potential proxy for RRV (Heidinger et al. 2012). These nucleoprotein structures found at the ends of chromosomes shorten with every cell division (Vaziri et al. 1994) and in response to physiological (Hau et al. 2015) or reproductive stress (Heidinger et al. 2012). Telomeres are increasingly regarded as suitable proxies for biological state (Mizutani et al. 2013;Hau et al. 2015;Angelier et al. 2019), where individuals in a poor biological state face reduced life expectancy. Consequently, there tends to be a positive correlation between telomere length and survival (e.g., Heidinger et al. 2012), likely resulting in higher RRV. Although increased reproductive effort may accelerate telomere shortening (e.g., Sudyka et al. 2014), this relationship is far from universal (Beaulieu et al. 2011;Plot et al. 2012;Bauch et al. 2013), suggesting that some individuals may be able to allocate more resources to maintaining both telomeres and reproductive effort (Cerchiara et al. 2017).
Here, we investigate the trade-off between CRI and RRV under varying colony-specific predation risk and for individually variable faecal glucocorticoid metabolite [fGCM] levels. Elevated baseline GCs could promote the re-allocation of resources from reproduction to survival when the value of CRI is expected to be lower (e.g. Wingfield et al. 1998), such as under high predation risk. First, we tested the assumption that higher predation risk is associated with elevated baseline fGCM levels. Drawing on the cost of reproduction hypothesis assuming a trade-off between CRI and RRV, we predicted that (1) individuals exposed to higher predation risk will decrease their CRI, (2) this effect is pronounced in individuals with higher RRV (Clark 1994) and that (3) adjustment of initial reproductive investment to prevailing predation risk will positively influence reproductive success. Finally, we assessed the validity of female body condition, telomere length and breeding experience as individual quality attributes, expecting positive associations with hatching success (Clutton-Brock 1988).
To this end, we studied long-lived female Eider ducks (Somateria mollissima) breeding in the northern Baltic Sea, in a study population subject to intense predation pressure on adults, eggs and ducklings (Öst et al. 2018). We used clutch size as a proxy for CRI, because Eiders mainly utilize stored reserves to produce the yolk (Hobson et al. 2015) and because larger clutches are more costly to care for (Hanssen et al. 2005;Neggazi et al. 2016). As proxies for RRV, we used female minimum breeding experience (in years), body condition at hatching and relative telomere length. The body condition of first-time breeders in Eiders has been found to be positively related to their subsequent survival (Ekroos et al. 2012) and female body condition is positively correlated with reproductive success (Lehikoinen et al. 2010). Furthermore, female body condition exhibits limited plasticity between seasons Ekroos et al. 2012). Thus, body condition should be regarded as a suitable proxy for RRV. Telomere length may represent a complimentary measure of RRV. Thus, although the telomere length of adult female Eiders is not significantly associated with age or body condition, reproductive success may be achieved at the expense of telomere length (Noreikiene et al. 2017). Female reproductive success was measured as proportional hatching success. Faecal glucocorticoid metabolite [fGCM] levels provided an integrated measure of adrenocortical activity over hours to days (Nilsson et al. 2008;Legagneux et al. 2011).

Study area and species
This study was conducted in 2011 at Tvärminne Zoological Station, western Gulf of Finland (59°50′ N, 23°15′ E). A map of the study area can be found in Waldeck et al. (2008). Eiders breed on the majority of islands in this archipelago, comprising both forested and non-forested islands. Eiders arrive at the breeding grounds at the end of March-first half of April and egg-laying typically occurs at the end of Aprilfirst half of May (Lehikoinen et al. 2006). Re-nesting after failed breeding attempts is extremely rare (Öst et al. 2018). Females take regular recesses from incubation, mainly to drink, as fatty acid profiling of incubating females suggests that feeding during incubation is uncommon (Hellman 2017).

Female trapping and sampling
Nesting Eider females (N = 197) were captured with hand nets on 20 different islands (mean ± SD = 9.9 ± 11.7 individuals per island), weighed to the nearest 10 g, measured for structural size (radius-ulna length to the nearest 1 mm), and ringed with a standard metal ring. Trapping was timed to roughly coincide with the latter part of the incubation period because the risk of human-induced nest desertion is significantly lower if Eider nests are visited during the later phases of incubation (Bolduc and Guillemette 2003). Because females are not ringed as ducklings, we used a relative age indicator which estimates (minimum) years of breeding experience (mean ± SD = 2.6 ± 3.68 years; range [0-16]), given by the number of years since first capture (Öst and Steele 2010). This variable approximates female age (Jaatinen and Öst 2011) since females exhibit strong breeding philopatry (Öst et al. 2011) and therefore previously unringed females are younger than those already ringed. The majority of breeding females are trapped annually (2011: 61.3%), with a constant trapping effort for the past decades . Nevertheless, we acknowledge a lack of precision in the variable, because it scores all unringed females as inexperienced, and ignores variation in the age at first breeding (typically 3 years, range 2-5 years; Hario and Rintala 2009). Nevertheless, this imprecision renders the current analysis conservative in its ability to detect effects related to female age.
Female blood samples (< 1 ml) were extracted from a brachial vein. When alarmed, incubating female Eiders readily spray faeces on their eggs that deter predators. Faecal samples were immediately collected from the female or from the nest cup and placed into Whirl-Pak (Nasco) bags. Concentrations of fGCM were unaffected by our handling since all samples were collected before any capture-induced increases in fGCM levels are detectable (6-10 min after capture in geese, Legagneux et al. 2011). Blood and faecal samples were immediately stored on ice and transferred to − 20 °C within 4 h of sampling.

Reproductive investment and female body condition
Mean clutch size was 4.5 eggs (range 2-7). Conspecific brood parasitism is rare (ca 6% of eggs; Waldeck et al. 2004). We used egg flotation (Kilpi and Lindström 1997) to determine the incubation stage (mean ± SD = 16.49 ± 4.50 days of incubation), laying and hatching dates. The estimated and true incubation stages do not statistically differ from each other (Kilpi and Lindström 1997). Female body condition indices were determined if eggs had been incubated for at least 8 days (egg-laying may otherwise still be in progress; Öst et al. 2008). This index was given by the standardized residuals of a regression of log-transformed projected weight at hatching (response variable) on log-transformed radiusulna length (Öst and Steele 2010). Projected weight at hatching was obtained by subtracting an estimate of the expected body weight loss during the remaining incubation time from measured incubation body weight at capture. Because females abstain from feeding during incubation and incubation stage at capture is variable, we derived an estimate of mean weight loss rate during incubation as the slope of the regression of log-transformed body mass (response variable) on log-transformed incubation time and projected hatching date (Öst et al. 2008).

Hatching success
Hatching success was determined during a nest revisit at the expected hatching date based on egg flotation. Successful nests contained either live ducklings or eggshells with intact membranes (Öst and Steele 2010). Such eggshell remnants are clearly distinct from those associated with nest depredation; in the latter case eggshells are shattered and inner egg membranes are bloody and attached to the eggshell fragments. Since hatching status was based on systematic daily nest monitoring, any signs of successful hatching or egg predation were fresh and not appreciably affected by exposure to weather conditions. If the brood had hatched and left prior to our arrival, leaving only intact egg membranes (indicating successful hatching) with no remains of depredated eggs, the viable clutch size was deemed the same as the clutch size recorded at female capture . We recorded the number of successfully hatched eggs (hatched ducklings surviving to leave the nest) and unhatched eggs (eggs failing to hatch due to depredation, abandonment and occasionally inviability). We were able to unequivocally determine nest fate in 165 out of the 197 nests (83.8%).

Predation risk
The two most important predators of adult female Eiders are the White-tailed Sea Eagle (Haliaeetus albicilla) and the American Mink (Neovison vison) (Öst et al. 2018). Clutches are depredated by Hooded Crows (Corvus corone cornix), Ravens (Corvus corax) and large gulls (Larus spp.) . We quantified colony-specific predation risk (mean ± SD = 0.17 ± 0.11; range [0-0.67]) by calculating the proportion of depredated nests at the first encounter of all nesting attempts on a given island (Öst et al. 2011). The total number of nesting attempts was the sum of actively incubated, early-depredated and abandoned nests, and nests that had already hatched when first found. Basing the predation index on the number of initial nest predation events was necessary for two reasons. First, this yielded an index unaffected by our nest monitoring that may induce additional nest depredation. Second, nests of incubating and blood-sampled females that were later found depredated (see "Hatching success" above) were not included in the count of early-depredated nests, which was used for calculating the colony-specific predation risk. Thus, we avoided statistical circularity otherwise arising due to an association between proportional hatching success (response variable) and colony-specific predation risk (independent variable). Our predation risk index is likely unbiased by variation in the timing of initial nest monitoring relative to island-specific nesting phenology. Thus, the time difference between the island-specific median date of nest visits and median estimated hatching dates (mean ± SD = − 12.19 ± 3.35 day s) was unrelated to colony-specific predation risk (r = 0.02, p = 0.93, N = 21). The predation index accurately captures colony-specific danger because there is a tight correlation (r = 0.75) between predation on incubating females and eggs (Öst et al. 2018), based on the same nest predation index as here.
Inclusion of proportional nest cover as a covariate was necessary to account for a potential increase in clutch size with increasing nest concealment, as (often undetected) partial clutch predation may disproportionately affect nests offering less protective covering (Öst et al. 2008). Female Eiders with different glucocorticoid profiles may also differ in their nest cover preferences (D'Alba et al. 2011;Seltmann et al. 2014). To quantify the structural cover of the nest, digital photographs were taken by placing an upwardfacing Olympus C-740 camera equipped with a 42-mm Opteka semi-fisheye lens in the nest (Öst and Steele 2010). The software Image Tool 3.00 (University of Texas Health Science Center, San Antonio) was used to convert images to greyscale and the pixels were then converted to either black or white using a threshold so that sky was white and all rocks and vegetation were black (Öst et al. 2008). The proportion of black pixels was used as a quantitative measure of nest cover.

Faecal glucocorticoid metabolites
Faecal glucocorticoid metabolites (fGCM) were measured by radioimmunoassay (RIA) using a double-antibody kit (ImmuChemTM Double Antibody, Corticosterone, 1251 RIA Kit, MP Biomedicals, Orangeburg, NY). The RIA kit protocol for Eider fGCM measurements in this population is described by Jaatinen et al. (2013). The mean sensitivity of the fGCM assay was 10.6 ng/g and the mean (± SD) fGCM concentration in the samples was 222.25 ng/g (± 136.77) (range 43.0-961.1 ng/ml, N = 185). The intra-assay CV was 5% and inter-assay CV was 12.29%. As a biological validation of the assay, injections of a corticosterone-releasing hormone (adrenocorticotropic hormone, ACTH) induce increases in fGCM concentrations 1-4 h post-injection in two other waterfowl species (Nilsson et al. 2008;Legagneux et al. 2011). Although the determination of female fGCMs took place after clutch completion, this is not a critical limitation. First, female faeces start to accumulate immediately after incubation onset, typically occurring after the second egg is laid (e.g., Andersson and Waldeck 2006). Thus, fGCM concentrations are likely to provide an integrated measure of faecal hormone profiles accumulated over a period of several days. Second, fGCM levels do not systematically change with advancing incubation (Jaatinen et al. 2013).

Telomere length
We measured relative telomere length (TL) in female Eider blood using real-time quantitative PCR (qPCR ;Cawthon 2002;Criscuolo et al. 2009). This method has been validated for use in ecological studies (Criscuolo et al. 2009;Heidinger et al. 2012), and optimized for use in Eiders (Noreikiene et al. 2017;Öst et al. 2020). We, therefore, provide a detailed account of the telomere qPCR as supplementary information (Supplementary file 1).

Data analysis
All analyses were conducted using the statistical software R 3.4.3 (R Development Core Team 2017). Because predation risk was assessed at the colony level, we tested the association between colony-specific predation risk and mean colony-specific fGCM level, after log-transforming individual fGCM values to meet the normality assumption required of Pearson correlation.
To analyse variation in clutch size, we constructed a linear mixed model including measures of RRV (minimum breeding experience in years, body condition, TL), colonyspecific predation risk, and female fGCM level. The initial model included the two-way interactions between fGCM and each RRV variable and between fGCM and predation risk, respectively. These interactions were selected because rather than expecting a simple, linear association between fGCM and the response variables (clutch size and hatching success), we should consider both RRV and predation risk as potential moderators between fGCM levels and our response variables (see "Introduction"). An additional interaction between minimum breeding experience and predation risk was added to control for possible age-associated changes in investment decisions under predation risk. Proportional nest cover was included as a covariate. Colony identity was included as a random effect to account for the non-independence of observations from the same colony. This random effect may potentially exhibit redundancy with the colony-level variable predation risk. Nevertheless, inclusion of colony identity was warranted because of its nonzero variance (variance (± 95% CI) = 0.175 (0.021-1.47)). Furthermore, a model without a random effect did not offer a better fit (likelihood-ratio test (LRT): χ 2 1 = 0.21, p = 0.64), indicating the presence of variation among colonies not captured by the fixed effects.
To analyse variation in proportional hatching success, created by combining the number of viable and nonviable offspring using 'cbind' in R, we constructed a generalized linear mixed model with binomial error distribution (logit link function) and Laplace parameter estimation. The fixed effects were identical as in the clutch size model, except for the addition of timing of breeding, to control for a seasonal decrease in nest success (Jaatinen et al. 2013), and clutch size, to control for potential clutch-size dependent hatchability. Inclusion of colony identity controlled for pseudoreplication. Colony identity had a variance component exceeding zero [variance (± 95% CI) = 0.40 (0.07-1.75)], and its inclusion significantly improved model fit (LRT: χ 2 1 = 11.69, p < 0.001). Additionally, because the incubation stage at capture may influence nest desertion (Bolduc and Guillemette 2003), we tested whether the addition of incubation stage to the final model of proportional hatching success (containing only significant terms; see below) was necessary. Incubation stage at capture was not significant (p = 0.09), and the results remained quantitatively almost identical, with no loss or gain of significance for any term. We, therefore, excluded the incubation stage at capture from the final model presented here.
The significance of each fixed effect was assessed by comparing models with and without the term using likelihood ratio tests until only significant terms (α < 0.05) remained (Tables 1, 2). Backward stepwise elimination is a conservative model selection method that, despite being criticized, often performs similarly to alternative approaches (Murtaugh 2009).
The residuals of all models adhered to the assumption of independence and normality. We tested for multicollinearity by calculating the variance inflation factors (VIFs) for the independent variables in the two models. VIFs, calculated using the vif function in the car package in program R, are defined as 1/(1 − r 2 j ), where r 2 j is the coefficient of determination of the linear regression of regressor j on the other regressors in the model. All variance inflation factors were below 1.13, indicating that no multicollinearity was present, thus allowing the variables to be present in the same models.
Significant interaction terms between continuous variables were illustrated following Aiken and West (1991). Briefly, regression equations were restructured to reflect the regression of the criterion on one predictor and simple slope regressions were plotted to display the interaction at low (25% quartile), intermediate (median) and high (75% quartile) levels of the moderator variable. All other independent variables in the model are held constant, by centring them around their respective sample means. Although there is no consensus on how to estimate the effect size of interactions in mixed models, we reported the difference between the marginal coefficient of determination (R 2 ) (i.e., the proportion of variance explained by the fixed effects alone) of the full model and the marginal R 2 of the model lacking the focal significant interaction term, as proposed by Nakagawa and Schielzeth (2013). In addition, we reported the total amount of variation accounted for by all the significant interactions together, to assess the relative importance of main effects versus interactions in explaining the variation in clutch size and hatching success, respectively. The marginal R 2 of models were calculated with the MuMIn package (Barton 2018).

Association between fGCM levels and predation
There was a significant positive correlation between colony-specific predation risk and mean colony-specific logtransformed fGCM levels of females (Pearson correlation: r = 0.49, p = 0.025, N = 21).

Reproductive investment
Colony-specific predation risk and fGCM interactively explained variation in clutch size (colony-specific predation × fGCM, p < 0.001; Table 1). Thus, clutch size decreased with increasing fGCM levels when predation risk was high, but was nearly independent of fGCM levels when predation risk was low or intermediate (Fig. 1). The colonyspecific predation by fGCM interaction explained 4.2% of the total variance (marginal R 2 = 0.197 with interaction vs. marginal R 2 = 0.155 excluding the interaction). Similarly, female body condition and fGCM level had a significant interactive effect on clutch size (body condition × fGCM, p = 0.016; Table 1). The clutch size of good-condition females decreased with increasing fGCM, whereas clutch size showed no obvious relationship with fGCM levels in intermediate and poor-condition females (Fig. 2). The body condition by fGCM interaction explained 4.2% of the total variance (marginal R 2 = 0.197 with interaction vs. marginal R 2 = 0.155 excluding the interaction). Female minimum breeding experience and telomere length, in an interactive manner, explained variation in clutch size (breeding experience × telomere length, p = 0.049; Table 1). Younger females tended to show a slightly negative correlation between clutch size and telomere length whereas more experienced (i.e. older) females displayed a slightly positive correlation between clutch size and telomere length (Fig. 3). This third significant interaction explained 1.8% of the total variance . The values for clutch size are jittered to make overlapping points visible. Note that this figure is for illustrative purposes, and that in our analyses we treated the variables included in the interaction term as continuous predictors (marginal R 2 = 0.197 with interaction vs. marginal R 2 = 0.179 excluding the interaction). Clutch size also increased with nest cover (Table 1). The three aforementioned significant interactions together accounted for 13% of the total variance in clutch size (marginal R 2 = 0.197 with interactions vs. marginal R 2 = 0.0671 excluding the interactions), a value to be compared with 6.7% of the variation explained by the fixed main effects alone.

Hatching success
Colony-specific predation risk and fGCM interactively explained variation in hatching success (colony-specific predation × fGCM, p < 0.001; Table 2). This interaction revealed that the drop in hatching success with increasing fGCM levels was least pronounced under high nest predation risk (Fig. 4). The colony-specific predation by fGCM interaction explained 4.5% of the total variance (marginal R 2 = 0.374 with interaction vs. marginal R 2 = 0.33 excluding the interaction). Hatching success increased with female relative telomere length, but this increase was modulated by breeding experience (breeding experience × telomere length, p = 0.002; Table 2). Older females generally showed higher hatching success than younger females, but the increase in hatching success with relative telomere length was accentuated for younger breeders (Fig. 5). This second interaction explained 0.8% of the total variance (marginal R 2 = 0.374 with interaction vs. marginal R 2 = 0.366 excluding the interaction). There was also a significant interaction between breeding experience and colony-specific predation risk (breeding experience × colony-specific predation, p = 0.018; Hatching success increased with relative telomere length of female Eiders, but this increase was modulated by minimum years of maternal experience, a proxy for age. The increase in hatching success with telomere length was most pronounced for young females (25 percentile, dotted line, open circles) compared to medium-aged (50 percentile, dashed line, grey dots) and older females (75 percentile, solid line, black dots). Note that this figure is for illustrative purposes, and that in our analyses we treated the variables included in the interaction term as continuous predictors Table 2). Hatching success increased with female breeding experience, but most markedly so under low predation risk (Fig. 6). The variance explained by this third interaction was 2.1% (marginal R 2 = 0.374 with interaction vs. marginal R 2 = 0.353 excluding the interaction). Hatching success also increased with maternal body condition, but decreased with later laying date and increasing nest cover ( Table 2). The three significant interactions together accounted for 6.9% of the total variance in clutch size (marginal R 2 = 0.374 with interactions vs. marginal R 2 = 0.305 excluding the interactions), a value to be compared with 30.5% of the variation explained by the fixed main effects alone.

Discussion
Our results supported the expected positive relationship between predation risk and fGCM. Consistent with our first prediction of reduced CRI with increasing predation risk, clutch size decreased with increasing fGCM levels under high predation risk. Elevated perceived predation risk may promote physiological changes that aid survival but reduce investment in egg production (Bonier et al. 2009a;Clinchy et al. 2011;Zanette et al. 2011). In contrast, females breeding in areas with lower predation risk showed, if anything, a weak positive association between clutch size and maternal fGCM levels (Fig. 1). Particularly pertinent is the variation in fGCM levels under high and low threat from predators. Female fGCM levels under high predation risk were numerically, although not statistically significantly, higher (238.9 ± 138.3 (SD) ng/g) than under low predation risk (196.3 ± 92.1 ng/g) (t 78.71 = − 1.73, p = 0.088). Thus, moderately elevated fGCMs, potentially characterizing females nesting in safer colonies, may be associated with increased reproductive allocation, e.g. through more efficient acquisition of food resources for egg production ('cort-adaptation hypothesis'; Bonier et al. 2009a;Crossin et al. 2012;Love et al. 2014). Potentially, the apparent decline in clutch size under high predation risk in females with high fGCM levels may be an artefact of undetected partial clutch predation that is particularly common during early laying, when nest attendance is low (Öst et al. 2008). However, this possibility seems unlikely for two reasons. First, the proportion of successfully hatched offspring decreased less with increasing fGCM in more risky areas (Fig. 4). Second, well-covered nests had lower proportional hatching success (Table 2), whereas (early) partial clutch predation is expected to be lower in such nests (Öst et al. 2008).
It is important to recognize that clutch size is but one measure of reproductive investment. In Eiders, egg mass is closely positively correlated with duckling body condition (Öst et al. 2020). Considering the potential trade-off between clutch size and egg mass (Stearns 1992), the question arises as to whether the observed changes in clutch size may be mitigated by compensatory changes in egg mass. However, this possibility is unlikely in the present study. Thus, there is no correlation between clutch size and mean egg mass per clutch (r = − 0.094, p = 0.19, N = 190), standardized for incubation stage according to Öst et al. (2020). Thus, from a cross-sectional perspective, female Eiders producing larger clutches appear able to do so without compromising offspring quality.
In line with our second prediction that individuals with high RRV should reduce CRI to ensure future reproductive opportunities, females in good body condition produced smaller clutches with increasing fGCM levels (Fig. 2). Females in good body condition are likely to have high RRV, as evidenced by their higher survival (Ekroos et al. 2012), nest success (Lehikoinen et al. 2010) and post-hatching duckling survival (Bustnes and Erikstad 1991) than females in poor body condition. Higher RRV is expected to make animals risk averse when confronted with environmental challenges or reproductive stress (Clark 1994;Wolf et al. 2007), in this case, reflected in elevated fGCMs. Producing a smaller clutch may reduce metabolic demands and the time spent laying and incubating eggs, so minimizing exposure to predators (Martin 2002). One possible mechanism may relate to condition-dependent resource allocation to reproduction. Thus, female Eiders in good body condition breed earlier Laursen et al. 2019), and more extensively utilize stored reserves for producing their eggs  . Note that this figure is for illustrative purposes, and that in our analyses we treated the variables included in the interaction term as continuous predictors be optimal for females in good body condition to promptly initiate breeding after arrival, so as not to deplete stored energy reserves. In contrast, females in poor condition may maximize their reproductive output by intense foraging after arrival, which delays breeding. Such intense feeding activity may be facilitated by increased corticosterone levels (Crossin et al. 2012). Indeed, late-nesting females show a pronounced increase in baseline plasma corticosterone levels with increased reproductive effort (indexed by total clutch mass) (Jaatinen et al. 2013). Differential breeding resource allocation strategies may, therefore, potentially explain the differential relationship between clutch size and fGCM levels (Fig. 2) when comparing females in good and poor body condition, respectively. Reproductive effort has been shown to cause telomere attrition in mammals (e.g., Kotrschal et al. 2007) and birds (e.g., Heidinger et al. 2012). We showed that younger females with longer telomeres tended to produce smaller clutches (Fig. 3), although we advise caution in interpreting this result given the low effect size of the interaction between breeding experience and telomere length. It is plausible that such females have high RRV, providing an incentive for them not to put excessive effort into their current reproduction which may shorten their telomeres and reduce their future reproductive prospects (Kotrschal et al. 2007). The clutch sizes of older females, on the other hand, were weakly positively correlated with telomere length (Fig. 3). Telomere length may predict remaining lifespan (Heidinger et al. 2012), and reproductive capacity tends to increase with age (Figs. 5, 6;Clutton-Brock 1988). If so, a positive correlation between lifespan, potentially associated with longer telomeres, and reproductive performance may be expected in older, more experienced females. It is noteworthy that the relative telomere length of adult female eiders shows no significant trend with age (Noreikiene et al. 2017). However, because both this study and our current findings are based on cross-sectional observational data, longitudinal data on telomere length will now be needed to infer any causal connections between telomere length, breeding experience and fitness.
We observed a negative association between fGCM and hatching success (Fig. 4) in support of the corticosteronefitness hypothesis (Bonier et al. 2009a). However, we find that females in high predation areas sustained relatively better performance with increasing fGCM, concomitant with a reduction in clutch size (Fig. 1). Hence, reduced clutch size in response to increased predation risk may ultimately provide a relative reproductive advantage, consistent with our third prediction. The higher hatching success of these females' clutches was not a by-product of clutch sizedependent hatching success (Table 2).
Increased glucocorticoid reactivity upon immediate threat (plasma CORT: Jaatinen et al. 2014) or during prolonged exposure (fGCM; this study) may confer benefits under predation risk by preparing individuals for potential threats (Romero 2002). One possible mechanism is through a more efficient escape from predation. Thus, higher handlinginduced plasma corticosterone reactivity is associated with an earlier escape response from the nest (Seltmann et al. 2012). Female eiders may successfully hatch a clutch after managing to escape a predator attack (pers. obs. using wildlife cameras), indicating a potential association between increased glucocorticoid reactivity and hatching success, mediated through more efficient avoidance of being captured once attacked. Second, higher handling-induced glucocorticoid levels are associated with a shorter incubation period (Seltmann et al. 2012), which should reduce exposure time to predators. To conclude, females with elevated maternal fGCM levels may suffer a relatively smaller loss of reproductive success in areas of high nest predation risk, where predation on adult females is also intense (Öst et al. 2018).
Our results indicate the utility of telomere length and age as indicators of individual quality: hatching success increased with maternal telomere length (Fig. 5) and years of breeding experience (Figs. 5,6). Likewise, telomere length is correlated with individual quality in other longlived birds (Pauliny et al. 2006;Vaillant et al. 2015;Angelier et al. 2019), and breeding experience is a key determinant of avian breeding performance (e.g., Forslund and Pärt 1995). A potential mechanism for age-related improvements in breeding competence may involve the pituitary hormone prolactin. Prolactin is essential for regulating parental care and older parents are typically more adept at maintaining adequate levels of this hormone (Angelier et al. 2007). Alternatively, the observed age-associated improvement in reproductive performance may reflect selective disappearance (Curio 1983). Longitudinal studies would be needed to separate within-individual (individual improvement) from between-individual (selective disappearance) effects of age on reproductive performance.
To conclude, increased predation risk was coupled with elevated fGCM levels of breeding females. We further showed that high potential for future breeding prospects and increased predation risk may restrain current reproductive investment, with glucocorticoids playing a role in the resolution of this trade-off. Thus, elevated fGCM levels were negatively correlated with CRI under high predation threat and for females in good body condition. Furthermore, elevated maternal glucocorticoid levels may serve as a buffer against the loss of reproductive success under severe predation pressure. Lastly, we showed that body condition, telomere length and breeding experience (age) can serve as phenotypic quality indicators in Eiders, the latter two having interactive effects on CRI and hatching success. However, two final cautionary notes are in order. First, it should be noted that the interactive effects involving hatching success were generally relatively weak compared to those involving clutch size. Second, because our correlational data preclude causal interpretation, the intricate relationships indicated by our study await further elucidation by longitudinal and experimental work involving not only predation but also competition and food limitation.