Age-related changes in sperm traits and evidence for aging costs of sperm production in a sexually promiscuous passerine

In many animal species, organismal performance declines with age in a process known as aging or senescence. Senescence typically leads to a deterioration of physiological functionality and can impact the development of primary sexual phenotypes. Sperm production is a complex and costly process that is sensitive to changes in individual physiological state, yet remarkably little is known about age-related changes in sperm performance and aging costs of sperm production. Here we use a non-linear generalized additive mixed models (GAMM) modelling to evaluate age-related changes in postcopulatory sexual traits in the European barn swallow (Hirundo rustica rustica), a relatively short lived sexually promiscuous passerine species, where male extra-pair fertilization success has been shown to increase with age. We confirmed a positive relationship between sperm midpiece length and sperm velocity in this species. Within-male changes in sperm morphology and sperm velocity were in general absent, with only sperm length decreasing linearly with increasing age, although this change was negligible compared to the overall variation in sperm size among males. In contrast, the cloacal protuberance (CP) size changed nonlinearly with age, with an initial increase between the first and third year of life followed by a plateau. The results further indicate the existence of a trade-off between investments in sperm production and survival as males with large CP tended to have a reduced lifespan. This seems consistent with the idea of expensive sperm production and survival aging costs associated with investments in post-copulatory traits in this sexually promiscuous species.


Introduction
Aging leads to a deterioration of physiological functionality and reduced individual fitness and reproductive performance (Partridge and Barton, 1993). A central idea of the contemporary theory of aging is the declining force of natural selection with increasing age (Medawar, 1952; see also the antagonistic pleiotropy theory; Williams, 1957). Because individuals typically have limited resources they invest in growth, reproduction and somatic maintenance (Kirkwood and Rose, 1991), this potentially leads to an evolutionary trade-off between reproductive effort early in adulthood and aging, i.e., there is a trade-off between allocation of resources to reproduction versus somatic Frontiers in Ecology and Evolution 02 frontiersin.org maintenance (disposable soma theory of aging, Kirkwood, 1977). So far, evidence for these allocation trade-offs is mostly available for females (Lemaître et al., 2015;Lemaître and Gaillard, 2017). In males, the existing studies evaluating potential trade-offs between reproduction and aging mostly focused on the investment in secondary sexual traits (e.g., Evans, 2003;Robinson et al., 2006), while less attention is currently paid to postcopulatory traits, such as sperm production, performance and morphology reviewed in Lemaître and Gaillard (2017) and Lemaître et al. (2020). This is despite the fact that in sexually promiscuous species, sperm quality and numbers typically determine male fitness reviewed in Snook (2005) and sperm production is presumably an energetically demanding process (Hayward and Gillooly, 2011). Thus, resource allocation to postcopulatory traits may play an important role in early vs. late life investment trade-offs, thereby determining male life-history strategies (Lemaître et al., 2020). There is considerable variability in the rate of aging among unrelated phylogenetic lineages (Jones et al., 2014), among closely related species (Fletcher and Selman, 2015), or within species among individuals (Lecomte et al., 2010). Aging is typically associated with changes in physiological and behavioral traits, such as foraging efficiency, immune system efficiency, health or individual condition (Reimers et al., 2005;Catry et al., 2006;Angelier et al., 2007). Similarly, individual reproductive performance may deteriorate with increasing age (reproductive senescence; e.g., Reid et al., 2003;Radwan et al., 2005; but see Bouwman et al., 2007). However, given the assumed trade-off between reproduction and somatic maintenance, it has also been hypothesized that individual strategies to optimize lifetime fitness may involve a change of life history allocation based on life expectancy, such that early in life individuals would invest primarily in survival, whereas towards the end of life they would invest more in reproduction at the cost of self-maintenance (Jones et al., 2014). This "terminal investment" hypothesis (Clutton-Brock, 1984) does not predict a deterioration of sperm traits with age but rather an increase of sperm quality and sperm production at terminal age.
Sperm phenotype is an important predictor of male reproductive success in sexually promiscuous species. The length of sperm (and its components) has been shown to play a role in determining sperm fertilization capacity across vertebrate taxa (LaMunyon and Ward, 1998;Knief et al., 2017). Sperm morphology, especially flagellum or midpiece length, can affect sperm velocity, another relevant factor affecting male fertilization success (Firman and Simmons, 2010;Knief et al., 2017;reviewed in Gomendio and Roldan, 2008). Sperm morphometry, and to a lesser extent sperm velocity, is genetically determined (Birkhead et al., 2005;Mossman et al., 2009), but there is evidence of plasticity caused by various factors, such as dietary quality (Støstad et al., 2019), redox state (Tomášek et al., 2017), stage of the breeding season (Edme et al., 2019), social environment (Immler et al., 2010), or immune reaction and health status (Losdat et al., 2011;Svobodová et al., 2018). These observations imply that, similar to precopulatory sexual traits (Andersson 1994), sperm traits may be condition dependent. Moreover, increased sperm production (and large ejaculate volumes) may allow males to transfer more sperm or increase the copulation frequency compared to other males, thus outcompeting males with low sperm production (Schulte- Hostedde and Millar, 2004;Boschetto et al., 2011); however, sperm production seems to be costly, as evidenced by the fact that males adjust ejaculate quality and quantity to the level of sperm competition (e.g., Ramm and Stockley, 2009) and that sperm production is often restricted to short periods of the year when sexually receptive females are present in the population (Calhim and Birkhead, 2007;Lüpold et al., 2012).
If sperm production is costly, it could be assumed that sperm quantity and performance change with age (Lemaître et al., 2020). However, the results of the few studies analyzing the correlation between age and sperm traits are inconsistent. There could be no age-related pattern in sperm sizes (Girndt et al., 2019) or older males produce longer sperm (Gasparini et al., 2010;Langen et al., 2017) or sperm with reduced velocity (Pasqualotto et al., 2005;Gasparini et al., 2010;Meunier et al., 2022). Similarly, old males could produce either higher (Schiavone et al., 2012;Jin et al., 2016) or lower (Pasqualotto et al., 2005;Langen et al., 2017;Meunier et al., 2022) sperm numbers compared to young males. Unfortunately, most available studies are cross-sectional and distinguish only between two age categories (young and old males). Only a few studies, mostly on farm animals, evaluated within-male changes in sperm traits (Dean et al., 2010;Waheed et al., 2015) but focused on a restricted part of individual life cycle (see also Vega-Trejo et al., 2019 for a similar approach in fish). Exceptions are studies on domestic fowl (Gallus gallus; Cornwallis et al., 2014) and houbara bustards (Chlamydotis undulata; Preston et al., 2011), which showed that the association between age and sperm phenotypes may not always be so straightforward and sperm quality may change nonlinearly with male age. Both studies showed a decline in sperm quality traits with progressing age. In contrast, sperm traits did not decrease with male age in a longitudinal study on ants (Metzler et al., 2018).
In sexually promiscuous songbird species (Passeriformes), old males tend to be more successful than young males in both protecting the within-pair paternity (Bowers et al., 2015;Hsu et al., 2017) and obtaining extra-pair mates (Lifjeld et al., 2011). Although age-related changes in sperm quality and sperm production have been hypothesized to partly account for this phenomenon (e.g., Laskemoen et al., 2010;Lifjeld et al., 2022), data on age-related changes in sperm traits are typically not available for species studied with regard to extrapair fertilization success of young and old males (but see DuVal, 2012;Sardell and DuVal, 2014), and the evidence for age-related changes in sperm morphology and velocity is inconsistent (Laskemoen et al., 2010;Saetre et al., 2018;Edme et al., 2019;Lifjeld et al., 2022). On the other hand, old males may have larger testes than young males, which might indicate intensified sperm production (Laskemoen et al., 2008;Lifjeld et al., 2022). Along with the seasonal increase in testis mass, male passerines produce a pronounced swelling of the cloaca, known as cloacal protuberance (hereafter CP; Wolfson, 1952). The cloacal protuberance is a sperm storage organ the size of which has been shown to positively correlate with testis mass (Birkhead et al., 1993;Laskemoen et al., 2008), the level of sperm competition (Birkhead and Møller, 1998), reflects ejaculate volumes at within-species level (Laskemoen et al., 2010;Girndt et al., 2019) and tends to be positively associated with male fertilization success (Laskemoen et al., 2008(Laskemoen et al., , 2010. As with sperm traits, there is contradictory evidence for an association between CP sizes and male age, with some cross-sectional studies reporting no effect (Girndt et al., 2019) or a positive association between CP and male age (Laskemoen et al., 2008).
To date, most available evidence for age-related changes in sperm traits and sperm production available is based on cross-sectional studies (see above) but longitudinal data are urgently needed to assess withinmale age-related patterns in postcopulatory traits, as well as to reveal hidden survival and aging costs of potentially expensive and energetically demanding sperm production reviewed in Lemaître et al. (2020). Here we assess age related changes in sperm morphology traits, velocity, and CP sizes in longitudinally observed European barn swallow Frontiers in Ecology and Evolution 03 frontiersin.org (Hirundo rustica rustica) males (Michálková et al., 2019;Kauzálová et al., 2022) using generalized additive mixed models (Wood, 2017;Cooper et al., 2021) to construct aging trajectories of each trait across individual lifespans. The studied population of barn swallows is characterized by moderate to high levels of sperm competition, with older males more successful in gaining extra-pair fertilizations than younger males (Michálková et al., 2019). Such pattern could reflect no age-related deterioration of sperm traits in barn swallows, increased investments in sperm production in old males, their increased attractiveness to females (Hsu et al., 2015) but also selective disappearance of low-quality individuals from the population (van de Pol and Verhulst, 2006). In addition, longitudinal data allowed, for the first time in a free-living bird species, to assess the potential survival and aging costs of sperm production.

Study area and general field procedure
Free-living barn swallows were captured with mist nest approximately every 3 weeks from beginning of May to the end of July at four localities in South Bohemia, in the protected landscape area of Třeboňsko -Hamr in Lužnice (49°3′24.217"N, 14°46′9.361″E), Šaloun in Lomnice nad Lužnicí (49°04′07.6"N, 14°42′37.7″E), Břilice (49°01′13.4"N 14°44′17.5″E), and Stará Hlína (49°02′21.2"N 14°49′06.5″E) in the breeding seasons of 2010-2021. All birds were marked with an aluminum National Museum of Prague ring and an individual combination of color rings. The age of birds that were ringed as a nestling at our study area (intensive ringing of the population started in 2008) was known. Birds that were captured un-ringed were assumed as 1-year old that hatched outside our observed populations. This approach has been used in previous studies (Costanzo et al., 2017;Kauzálová et al., 2022) and is possible because of high breeding philopatry and fidelity of barn swallows. Adult birds that did not return to the breeding locality were presumed dead (Costanzo et al., 2017;Kauzálová et al., 2022). Only resident birds (repeatedly occurring on breeding grounds) with a known year of birth and death were included in the analysis.
To obtain CP sizes, we measured the height of male cloaca and its maximum dimension along the two perpendicular axes (d 1 -transverse and d 2 -longitudinal, see also Laskemoen et al., 2010). The size of CP was calculated by the formula of an ellipsoidal cylinder: height × π × 0.5d 1 × 0.5d 2 (Laskemoen et al., 2008(Laskemoen et al., , 2010. CP sized were obtained immediately before sperm collection (see below) to avoid size changes after cloacal massage. Barn swallows usually produce sperm throughout the entire breeding season (April to August). The dataset only includes sexually active males that provided a good quality ejaculate sample and were sampled at the peak of the breeding season (early May to second half of July). Using this dataset, we tested the effect of sampling date (linear and second-order polynomial) on CP size and found no relationship (Supplementary Table S1). All measurements were performed blindly with respect to knowledge of individual age.

Sperm morphology
Sample of male ejaculate for sperm morphology measurement was taken from each male using a non-invasive method of cloacal massage (Wolfson, 1952;Albrecht et al., 2013). The sample was placed in a 5% formalin solution until the smears were prepared. Using an automatic pipette, 7 μL of each sample were transferred to a slide and allowed to air-dry before rinsing with distilled water to remove impurities. The slides were viewed at 400× magnification using a light microscope (BX51, Olympus, Japan), digital camera (DP71, Olympus, Japan) and imaging software (QuickPHOTO Industrial, Olympus Japan). Passerine sperms have a spiral conformation, the head is helical and the midpiece is distinctly elongated along the flagellum (Humphreys, 1972). Any sperm that did not conform to the characteristic helical conformation was marked as abnormal and was not included in the analysis. However, the number of abnormal sperm was very low (typically <1% sperm cells in the ejaculate). Using QuickPHOTO Industrial software, images of 10 spermatozoa with regular morphology from each ejaculate were taken. Sperm components (head length, midpiece length, and tail length) were measured to obtain an average value of sperm morphology traits for each ejaculate, which was used in further analyses. Total sperm length was calculated as the sum of these three sperm components, while flagellum length was calculated as the sum of the midpiece and tail length. All sperm morphological traits were significantly but rather weakly correlated with each other (N = 921 ejaculates, all r < 0.50; Supplementary Figure S1; Supplementary Table S2), except for midpiece and tail (r = −0.73) and total sperm length and flagellum (r = 0.97). Tail is the remaining part of the flagellum where the midpiece does not extend, and, in contrast to the midpiece (Knief et al., 2017), its biological relevance is unclear. Therefore, we did not use this trait in further analyses of age-related changes in sperm morphological traits. We also did not analyze age-related changes in flagellum length, as this was strongly correlated with total sperm length (above).

Sperm velocity
Samples for sperm velocity analysis were immediately transferred to Dulbecco's Modified Eagle Medium (Invitrogen) prewarmed at 40°C. Then we transferred a small amount of sample onto a prewarmed 4-chambered microscope Leja slide (20 μL deep, Leja, Netherlands). Sperm performance was recorded at 100× magnification using microscope CX41 (Olympus, Japan) fitted with the thermo plate (MATS-U55S Tokai Hit, Olympus), phase contrast, digital camera UI-1540-C (Olympus) and Olympus software QuickPHOTO Industrial. Sperm performance was recorded at multiple locations on the slide to record a sufficient number of sperms.
The recordings were analyzed using the CEROS computer-assisted sperm analysis system (Hamilton Thorne, Inc., United States). In the statistical analysis, the curvilinear velocity (VCL) value was used. This value characterized sperm speed over the entire trajectory by determining speed at each point of the pathway and it is an appropriate measure of swimming speed because of the absence of directional cues (Kleven et al., 2009;Laskemoen et al., 2010). Only progressive tracks were included in analysis (i.e., static and slow tracks were eliminated to remove the potential effect of drift in the chamber; Kleven et al., 2009). Moreover, cases where multiple sperm pathways merged were eliminated and samples containing less than 20 motile sperms were excluded from analyses. Sperm velocity measurements were performed blindly with respect to knowledge of the individual age.
Frontiers in Ecology and Evolution 04 frontiersin.org

Statistical analysis
Statistical analysis was performed using R 4.2.0 (R Core Team, 2020). Repeatability was estimated for males with more than one observation in life using the package rptR (the number of bootstraps was set to 1.000, Stoffel et al., 2017). Coefficient of variation (CV) was calculated as SD/mean × 100. Associations between sperm morphology and velocity were analyzed using lmer models in the lme4 (Bates et al., 2014) and lmerTest packages (Kuznetsova et al., 2017) with male identity and capture year as a random effect. To explain relationships between selected variables and age we used GAMM in the mgcv package with implemented nonparametric smoothing functions because of the possibility that aging trajectories may not follow parametric functions. GAMM are extended versions of generalized linear mixed models which include nonparametric terms and they are also not limited by parametric functions (Wood, 2017). The effect of age was modelled using penalized thin plate regression splines which estimated relationships between age and variables with penalized additive smoothing functions determined by restricted maximum likelihood. The penalized additive value of the number of smoothing functions is the effective number of degrees of freedom (EDF). An effective number of degrees of freedom value 1 indicates that the relationship is linear, a higher value determines the polynomial of the corresponding order (2 quadratic, 3 cubic, etc.). Only one measurement in each year for each male was used in the models. All models also included lifespan as one of the explanatory variables and male identity and capture year as random effects. The use of male identity as random affect allowed us to control for within-individual between-season variation in reproductive traits. Including both age and lifespan in the model allows us to identify within-male changes in sperm traits and sperm production and test selective disappearance of individuals in relation to postcopulatory traits investments. For example, in the presence of selective survival of males with better sperm performance, analyses would give false-positive age-related trends for these traits only because poorer-quality males disappeared (van de Pol and Verhulst, 2006). In case nonlinearity was only weakly supported for a given trait, we also fitted a linear mixed model using lmer in the lmerTest package. Since our dataset contained only one male older than 6 years (lifespan 8 years, i.e., one data point for age 7 and 8 years, respectively), the effect of age on sperm traits was also tested after removing this longest living potentially highly influential male and the results are presented simultaneously. To evaluate effect size of age in GAM models, we report the proportion of deviance explained by the variable. The value was calculated as the difference between the proportion of deviance explained by the model that did and did not include age among predictors (Wood, 2017). The significance level was 0.05 for all analyses. Means are reported with their associated standard error (SE).

Sperm traits, cloacal protuberance sizes, and their within-male repeatability between seasons
An overview of the phenotypic variation in sperm traits and CP sizes along with sample sizes is provided in Table 1. Repeated measurements (at least two per male in different years) of sperm morphology, sperm velocity and CP sizes were available for 225 males (592 observations), 179 males (453 observations) and 217 males (568 observations), respectively. The within-male between-season repeatability was significant in all sperm traits and ranged between 0.24 and 0.80 (Table 2). While sperm morphology traits exhibited a remarkably high within-male repeatability, sperm velocity was a more variable trait, and the CP size exhibited the highest variability between seasons.

Age-related changes in sperm traits and cloacal protuberance size
Spline and parametric effects are shown in Table 3. Detailed results for GAMM for each trait (including R 2 , relevant test statistics and the variance explained by random effects), and proportion of deviance explained by the variable are available in Supplementary Tables S4-S8. There was a tendency for total sperm length to decrease linearly with age (EDF = 1.488, F = 2.427, p = 0.063). Linear mixed model of the same data also revealed a decrease of sperm length with age with a linear trend best supported (estimate: −0.127 ± 0.058, p = 0.030; compared to the secondorder non-linear model: Chi 2 = 1.741, ∆df = 1, p = 0.187; Figure 1A). Results remained unchanged when the exceptionally long-lived male was removed from the analysis (Supplementary Table S9; Supplementary Figure S2). Proportion of deviance explained by age in both models was only 0.2%. Head (EDF = 1.001, F = 2.061, p = 0.152) and midpiece (EDF = 1.214, F = 0.060, p = 0.910) showed no age-related trend. Results remained unchanged when the exceptionally long-lived male was removed from the analysis (Supplementary Tables S10, S11).
Sperm velocity showed no change with male age (EDF = 1.003, F = 0.004, p = 0.967). Results were similar when the exceptionally longlived male was removed from the analysis (Supplementary Table S12).
The CP size changed nonlinearly with increasing age (EDF = 3.285, F = 10.416, p < 0.001 for the whole dataset; EDF = 2.851, F = 13.241, p < 0.001 when the exceptionally long-lived male was removed; Supplementary Table S13). In early life, the CP increases in size and remains stable after the third year of life ( Figure 1B; Supplementary Figure S3, respectively). Proportion of deviance explained by age in model was 5.9%, respectively 5.8% when the exceptionally longlived male was removed. At the same time, CP size tended to be negatively correlated with male lifespan (estimate: −1.724 ± 0.969, p = 0.076 for the whole dataset and − 2.138 ± 1.022, p = 0.037 when the exceptionally longlived male was removed from the analysis).

Discussion
In this study we used longitudinal data to analyze within-male age-related changes in selected sperm traits and in sperm production Frontiers in Ecology and Evolution 05 frontiersin.org (the CP sizes) in a small sexually promiscuous passerine. Previous studies on free-living vertebrates have typically used a cross-sectional approach to study age-related changes in sexual traits potentially involved in post-copulatory mate choice. Moreover, at least in promiscuous passerine species, age was usually assessed at a very coarse scale, distinguishing only two age categories, i.e., young males and old males, making it impossible to analyze changes in sexual traits with age or potential survival and aging costs of sperm production.
Our results indicate that sperm traits exhibit moderate to high within-individual repeatability between seasons. Specifically, morphological sperm traits showed higher repeatability compared to sperm velocity and CP sizes. Lower repeatability of velocity compared to sperm morphology was recorded also in other passerines (e.g., Birkhead et al., 2005;Mossman et al., 2009;Opatová et al., 2016;Saetre et al., 2018). This may indicate that while sperm morphology is under significant genetic control (Birkhead et al., 2005;Knief et al., 2017) and less prone to be affected by environment (e.g., reduced phenotypic plasticity, see also Tomášek et al., 2017), sperm velocity and sperm production may be susceptible to adverse environmental effects, may reflect individual condition (e.g., Simmons and Emlen, 2006) and be influenced by aging processes. Despite differences in withinindividual repeatability of sperm morphology traits and sperm velocity, our results suggest a positive association between sperm velocity and midpiece length (and a negative between tail length and sperm velocity). Sperm cells with a longer midpiece-mitochondria (and hence a shorter remaining tail part) are probably able to produce larger amounts of ATP and more energy (Vladić et al., 2002;Rowe et al., 2013) to swim faster (Firman and Simmons, 2010). The positive effect of midpiece length and relative midpiece length on sperm velocity seems to be a general phenomenon in passerine birds (e.g., Laskemoen et al., 2010;Knief et al., 2017;Tomášek et al., 2017;but see Cramer et al., 2015).
No age-related changes in the length of the head and midpiece of sperm were detected in barn swallows, consistent with the relatively high repeatability of morphological traits between seasons (see above). There was evidence for an age-related change in total sperm length, with sperm becoming shorter with progressing age. However, age explained only 0.2% of the deviance in sperm length (Supplementary Table S4), and the within-male sperm shortening (~ 0.13 μm/year) was negligible compared to the inter-male variability in sperm lengths (71 to 95 μm; Table 1). Our results based on a longitudinal approach cannot be directly compared with existing published crosssectional data. Probably the only available longitudinal study mapping changes in sperm size with age is on ants where age had no effect on sperm size (Metzler et al., 2018). However, cross-sectional studies mostly found no changes in sperm morphology with age of males in   Frontiers in Ecology and Evolution 06 frontiersin.org various passerine species (Laskemoen et al., 2008;Saetre et al., 2018;Edme et al., 2019;Girndt et al., 2019) including North American subspecies of barn swallows (Lifjeld et al., 2022). On the contrary, some studies, albeit using a limited dataset, found sperm of older males longer than in young males (Cramer et al., 2020). In general, sperm morphology seems to be a relatively stable trait within the life of the male, which is in agreement with the relatively high observed heritability of sperm morphology traits in passerine species (Birkhead et al., 2005;Mossman et al., 2009;Edme et al., 2019). Although some environmental stressors may influence sperm morphology, e.g., oxidative stress (Tomášek et al., 2017), heat stress (Armengol et al., 2015) or pesticides (Urióstegui-Acosta et al., 2014), these stressors may not be important in wild populations or the susceptibility of individuals to these stressors may not be age-related. Although sperm velocity is a trait exhibiting reduced year to year consistency compared to sperm morphological traits and may be influenced by individual condition or redox state (Tomášek et al., 2017), the lack of association between male age and sperm velocity may reflect an intimate association between sperm velocity and midpiece size, the trait that is also unaffected by age in our population of barn swallows. In cross-sectional studies, sperm velocity did not change with male age in other passerines (Laskemoen et al., 2010;Saetre et al., 2018) but sperm of old males tended to swim faster than sperm of young males in North American subspecies of barn swallows (Lifjeld et al., 2022). However, cross-sectional and longitudinal studies may provide contrasting results as shown in feral fowl (Gallus g. domesticus) where no within-male age related change in velocity was found in a longitudinal study, while a cross-sectional approach revealed that older males had reduced sperm velocity compared to young males (Dean et al., 2010).
We found evidence for a previously unknown age-related pattern in CP dynamics in barn swallow males -within individuals, CP sizes increased with age, reaching a plateau when at the age of 3 years, and then remaining stable in size until the male disappeared from the population. The proportion of deviance explained by the age in the model was 5.9% (Supplementary Table S8), suggesting that the effect of age on CP size is not strong, but neither it is negligible. Age related changes in CP sizes in passerines have been documented in several previous cross-sectional studies. In a typical scenario, CP size is larger in older males than young males (Bouwman et al., 2007;Laskemoen et al., 2008; but see Girndt et al., 2019). Age related changes in CP sizes may reflect age related dynamics in testes mass -in North American subspecies of barn swallows testes size increased with age in a crosssectional study (Lifjeld et al., 2022). CP sizes may reflect ejaculate sizes and investments in sperm production (Laskemoen et al., 2010;Girndt et al., 2019). Thus, an enlargement of CP size with age could explain high extrapair fertilization success of old males in our study population (Michálková et al., 2019).
Lifespan of males was included in all models to evaluate possible selective disappearance (van de Pol and Verhulst, 2006) in relation to postcopulatory traits. Overall, CP size tended to be negatively associated with lifespan and individuals with enlarged CPs were prone to disappear disproportionally from the population, indicating that the production of large ejaculates is costly. It was believed that sperm production is a relatively cheap process (in contrast to egg production; Parker, 1970Parker, , 1982. However, comparing energy expenditure on gamete production between the sexes is problematic (Hayward and Gillooly, 2011). Our results indicate an existence of a trade-off between reproduction and survival, and aging costs of sperm production (sensu Lemaître et al., 2020). Consistent with costly sperm production is the reduction in testicular size during nutritional stress (Perry and Rowe, 2010), outside the breeding season (Calhim and Birkhead, 2007) or a trade-off between ejaculate quality and either immunity (Simmons, 2012) or pre-copulatory sexual traits (Simmons and Emlen, 2006;Preston et al., 2011). Similar to our study, males in domestic fowl with higher sperm production lived shorter lives (Cornwallis et al., 2014). In contrast, we did not find a relationship between lifespan and sperm morphology traits or sperm velocity. This may indicate producing large sperm quantities, rather than maintaining stable sperm length or sperm velocity, is a costly life-history trait in birds, a taxon where sperm numbers may be a major factor determining male fertilization success (Immler et al., 2011).
To conclude, our results indicate that the production of highly motile sperm with optimal morphology is likely cheaper than the production of large sperm quantities, as only CP sizes were associated with lifespan and underwent a complex dynamic with male age. In general, old males may be more attractive to females as extra-pair partners (but see Girndt et al., 2018), may have easier access to extrapair copulations than young males (Hsu et al., 2015), and the production of large ejaculates (allowing more copulations) may be rewarding even at the expense of reduced survival. Further studies should aim to reveal the complex relationship between CP sizes, male survival, and extra-pair paternity component of male life-time fitness in sexually promiscuous free-living songbird populations. Effect of male age (years) on (A) total sperm length (μm) and (B) cloacal protuberance size (mm 3 ). Line is predicted value from GAMM with 95% confidence intervals. Points are partial residuals from model that included lifespan as a covariate, along with random effects of ring number (ID) and capture year.
Frontiers in Ecology and Evolution 07 frontiersin.org

Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement
The animal study was reviewed and approved by Animal Care and Use Committee, The Czech Academy of Science; Animal Care and Use Committee, Charles University in Prague.

Author contributions
TA conceived the study. KM, OT, and TA designed the study. KM, OT, VJ, MŠ, LP, JA, and TA collected the data. KM and JA analyzed sperm data. KM, OT, LP, and TA performed statistical analyses. KM and TA drafted the first version of the manuscript. All authors contributed to the final version of the manuscript.