Skip to main content

Natural variation in yolk fatty acids, but not androgens, predicts offspring fitness in a wild bird



In egg-laying animals, mothers can influence the developmental environment and thus the phenotype of their offspring by secreting various substances into the egg yolk. In birds, recent studies have demonstrated that different yolk substances can interactively affect offspring phenotype, but the implications of such effects for offspring fitness and phenotype in natural populations have remained unclear. We measured natural variation in the content of 31 yolk components known to shape offspring phenotypes including steroid hormones, antioxidants and fatty acids in eggs of free-living great tits (Parus major) during two breeding seasons. We tested for relationships between yolk component groupings and offspring fitness and phenotypes.


Variation in hatchling and fledgling numbers was primarily explained by yolk fatty acids (including saturated, mono- and polyunsaturated fatty acids) - but not by androgen hormones and carotenoids, components previously considered to be major determinants of offspring phenotype. Fatty acids were also better predictors of variation in nestling oxidative status and size than androgens and carotenoids.


Our results suggest that fatty acids are important yolk substances that contribute to shaping offspring fitness and phenotype in free-living populations. Since polyunsaturated fatty acids cannot be produced de novo by the mother, but have to be obtained from the diet, these findings highlight potential mechanisms (e.g., weather, habitat quality, foraging ability) through which environmental variation may shape maternal effects and consequences for offspring. Our study represents an important first step towards unraveling interactive effects of multiple yolk substances on offspring fitness and phenotypes in free-living populations. It provides the basis for future experiments that will establish the pathways by which yolk components, singly and/or interactively, mediate maternal effects in natural populations.


Understanding the causes of phenotypic variation is a major goal in evolutionary biology. Maternal effects are now recognized as important factors that contribute to phenotypic variation (meta-analysis by [1]), although the magnitude and direction of maternal effects is still under debate [2, 3]. Mothers shape offspring phenotype not only through the genes they pass on to their offspring but also by influencing the environment their offspring experience during development [4,5,6]. In egg-laying species, mothers can influence the embryonic environment by allocating different resources to the yolk, including hormones, antioxidants and fatty acids, among other components [7,8,9,10,11,12,13,14]. The importance of some of these maternally transmitted compounds for phenotypic traits has been documented in various taxa (e.g., [9, 15, 16]). In birds, yolk steroid hormones, such as androgens and glucocorticoids, can influence offspring growth, competitive ability and survival (e.g., Ficedula albicollis; [17, 18] and reviewed by [19]), but they can also increase chick susceptibility to oxidative stress by increasing the production of reactive oxygen species or by impairing antioxidant defenses (e.g., Gallus gallus; [20, 21]). Maternally derived antioxidants, such as carotenoids or vitamin E, can promote growth (reviewed by [22]) and limit the negative consequences of increased oxidative stress by scavenging the reactive oxygen species produced during growth (e.g., Gallus gallus, Larus michahellis, Parus major; [23,24,25]). Fatty acids are an important source of energy; they can enhance the general viability and proper development of both embryos and nestlings (reviewed by [26, 27]). In particular, polyunsaturated fatty acids (PUFAs) are essential components for the formation of cell membranes, heart function and brain development (e.g., [28]). At the same time, PUFAs are susceptible to damage from lipid peroxidation by reactive oxygen species that are generated as by-products of offspring metabolism [29]. Thus far, the evidence supporting a relevance of fatty acids for the phenotypic development of embryos and nestlings has come mainly from studies on poultry or captive birds (reviewed by [30]), in which environmental conditions are typically designed to be benign. Yolk fatty acids have been shown to vary with environmental conditions in free-living birds (Cyanistes caeruleus, Parus major [31, 32]), but it remains unclear whether they play a role in shaping offspring fitness and phenotype in natural populations that experience frequent fluctuations in food availability and weather conditions.

Maternal effects through egg deposition are often studied by measuring the effect of single yolk components (in particular androgens) or groups of related components (e.g., steroid hormones) on offspring fitness and phenotypic traits (reviewed by [33, 34]). This approach has been pivotal for advancing our knowledge of how mothers affect offspring phenotype. However, maternal effects are by nature multivariate (reviewed by [33, 34]), with different groups of yolk components influencing similar nestling traits, i.e., by having interactive effects on offspring phenotype [35,36,37]. For instance, hatchling mass was reduced and oxidative stress increased in Japanese quails (Coturnix japonica) hatching from eggs injected with either testosterone or carotenoids [35]. But when both components were administered together, neither hatchling mass nor oxidative stress were affected. Similar compensatory effects were found in yellow-legged gulls (Larus michahellis) after eggs were simultaneously injected with corticosterone and vitamin E [36, 37]. Such studies suggest that the consequences of maternal effects on offspring phenotype known from studies that focus on single components might actually be absent, weakened or potentiated if the presence and actions of other components are also considered.

Here we studied the relationship of multi-substance yolk composition with nestling fitness and phenotype in a wild population of great tits over two years. This observational study is powerful in assessing natural variation in yolk composition and offspring fitness and phenotypes, and provides an important basis for subsequent experimental work. We measured the concentrations of 31 yolk components in the fourth egg from 69 clutches, including 4 steroid hormones, 3 antioxidants and 24 fatty acids. We grouped yolk components using a principal component analysis, although a discussion of the associations is beyond the scope of this manuscript. We then analysed the relationships between these groupings of yolk components and offspring fitness proxies, such as hatchling and fledgling number. We also determined links between yolk components and phenotypic traits like growth (i.e., body mass and tarsus length) and oxidative status of the nestlings from a given nest. The oxidative status of an individual is defined as the concentrations of pro-oxidants (i.e., reactive oxygen species) and antioxidants (i.e., non-enzymatic and enzymatic compounds) present in cells and tissues (reviewed by [38]). A change in any of the molecular components of the oxidative system in favor of pro-oxidants can damage crucial molecules, such as lipids, proteins and DNA (reviewed by [39]), with potential fitness consequences (e.g., [40] and meta-analysis by [41]). Nestling oxidative status has therefore been proposed to mediate their survival and health (e.g., Fregata magnificens [42]). Because the strength of maternal effects declines throughout offspring ontogeny (meta-analysis by [1]) and because the nestling oxidative status can change during the growth period (e.g., antioxidant defences; Taeniopygia guttata; [43]), we measured chick phenotypic traits at two time points: before and after nestlings reached exponential growth (i.e., 6 and 12 days post-hatching, respectively).


Covariation of yolk components

Using PCA, we identified the relationships among all 31 yolk components (Additional file 1). PC1 was negatively associated with vitamin E, one specific monounsaturated fatty acid (MUFA; 20:1n − 9) and all ω − 6 PUFAs (16:2n − 6, 18:2n − 6, two 18:3n − 6, 20:2n − 6, 20:3n − 6, 20:4n − 6, 22:4n − 6). PC2 was positively associated with four saturated fatty acids (SFAs; 15:0, 16:0, 17:0 and 18:0), all but one MUFA (18:1n − 9, 16:1n − 9, 16:1n − 7, 18:1n − 7) and all ω − 3 PUFAs (18:3n − 3, 20:5n − 3, 22:5n − 3, 22:6n − 3). Androgens (androstenedione, 5α-dihydrotestosterone and testosterone) and carotenoids (lutein and zeaxanthin) loaded positively onto PC3 whereas two MUFAs (18:1n − 9, 16:1n − 9) loaded negatively onto PC3. To simplify the description of the results and the visualization of the figures, hereafter all PCs are discussed in terms of positive loadings (i.e., we transformed PC1 into positive values). Also, as almost all MUFAs were included in PC2, whereas only one or two were present in PC1 and PC3, and the effects of different groups fatty acids are more relevant to our research question rather than specific fatty acids, the effects of MUFAs on nestling phenotype are discussed primarily based on results for PC2. Corticosterone was the only yolk component for which only low values loaded on the PCs. Corticosterone concentration is the least repeatable trait in yolk (R ~ 0.18; [10]), thus yolk corticosterone concentrations in the fourth egg do not reliably predict concentrations in the other eggs of the same clutch. Therefore, we decided to forego further analysis of corticosterone in this study.

Egg components and fitness proxies

Nests with eggs that had higher SFA, MUFA and ω − 3 PUFA concentrations had higher numbers of hatchlings and fledglings (Fig. 1, Additional file 2) than nests containing eggs with lower concentrations of these fatty acids. Furthermore, fledgling number was lowest in nests with eggs with increased concentrations of vitamin E and ω − 6 PUFAs. Neither androgens nor carotenoids explained fitness proxies. Fledgling mass and tarsus length were not related to any yolk component.

Fig. 1
figure 1

Relationships between egg yolk components and fitness proxies in great tit nestlings. Yolk components are grouped based on a principal component analysis. PC1 was mainly represented by high concentrations of vitamin E (α-tocopherol) and ω-6 polyunsaturated fatty acids (PUFAs), PC2 by high concentrations of saturated fatty acids (SFAs), mono-unsaturated fatty acids (MUFAs) and ω-3 PUFAs, and PC3 by high concentrations of androgens (androstenedione, 5α-dihydrotestosterone and testosterone) and carotenoids (lutein and zeaxanthin; Additional file 1). Shown are the standardized effect sizes with their corresponding 95% credible intervals (CrIs). Filled circles indicate statistically meaningful support (i.e., if the mean difference between compared estimates is higher than 0.95) for the effect of a set of yolk components on fitness proxies

Egg components and phenotypic traits

On day 6, nestlings that hatched from clutches with eggs containing high concentrations of SFAs, MUFAs and ω − 3 PUFAs had longer tarsi than nestlings that hatched from nests with eggs containing low concentrations of these fatty acids (Fig. 2a, Additional file 3). Yolk androgen and carotenoid concentrations were negatively related to nestling non-enzymatic antioxidant concentrations (i.e., OXY), whereas mass was not associated with any group of components. On day 12, the oxidative status of great tit nestlings was explained by all yolk component groups, but in different ways (Fig. 2b, Additional file 4). In particular, nestlings hatching from nests with eggs that had high concentrations of vitamin E and ω − 6 PUFAs had high concentrations of reactive oxygen metabolites (i.e., ROMs) in plasma and low levels of enzymatic antioxidant glutathione peroxidase (i.e., GPX) in red blood cells. Plasma OXY concentrations were negatively related to high concentrations of SFAs, MUFAs and ω − 3 PUFAs, and, as on day 6, also to high concentrations of androgens and carotenoids. Body mass and tarsus length were not related to any group of yolk components.

Fig. 2
figure 2

Relationships between egg yolk components and physiological traits in great tit nestlings on days (a) 6 and (b) 12. Yolk components are grouped based on a principal component analysis. PC1 was mainly represented by high concentrations of vitamin E (α-tocopherol) and ω-6 polyunsaturated fatty acids (PUFAs), PC2 by high concentrations of saturated fatty acids (SFAs), mono-unsaturated fatty acids (MUFAs) and ω-3 PUFAs, and PC3 by high concentrations of androgens (androstenedione, 5α-dihydrotestosterone and testosterone) and carotenoids (lutein and zeaxanthin; Additional file 1). Nestling oxidative status and body condition are indicated with yellow and green shading, respectively. Shown are the standardized effect sizes with their corresponding 95% credible intervals (CrIs). Filled circles indicate statistically meaningful support (i.e., if the mean difference between compared estimates is higher than 0.95) for the effect of a set of yolk components on physiological traits. Oxidative status parameters: OXY, non-enzymatic antioxidants present in plasma; GPX, enzymatic antioxidant present in red blood cells; ROMs, reactive oxygen metabolites present in plasma


Natural variation in fatty acid content was the only group of yolk components that explained variation in nestling fitness proxies in a free-living population of great tits. In contrast, androgens and carotenoids—yolk components that were previously shown to affect chick development and fitness (e.g., [17, 18] and reviews by [19, 22])—did not account for variation in fitness. Furthermore, the two yolk component groupings that included fatty acids were related to more aspects of the physiological phenotype of the nestlings (size and oxidative status) than androgens and carotenoids (only one measure of oxidative status). Hence, natural variation in the fatty acid content of yolk may play a major role in determining the fitness and phenotype of offspring in a wild bird population.

Egg components and fitness proxies

In our study, yolk fatty acids were related to both hatchling and fledgling numbers. Yolk fatty acids have been reported to promote fitness in captive invertebrates (e.g., Penaeus chinensis, Acartia erythraea; [44, 45]) and vertebrates (e.g., birds: Gallus gallus [46]), but so far ours is the first study to indicate a relationship with fitness-related traits in a free-living population exposed to natural variation in environmental conditions.

In general, the importance of yolk fatty acid composition for hatchling and fledgling numbers might depend on the female’s diet during egg laying. During the laying period, great tits change from feeding on seeds to predominantly feeding on caterpillars [47, 48]. Caterpillars are rich in SFAs, ω − 3 PUFAs (i.e., α-linolenic acid) and antioxidants [49,50,51], i.e., in substances that are thought to be beneficial for offspring fitness. However, if caterpillars are scarce because of low ambient temperatures or poor habitat quality, adult great tits rely on other invertebrates like arachnids [52], which are richer in ω − 6 PUFAs [50, 51]. Whilst this dietary flexibility allows parents to maintain food intake, it may compromise the optimal balance of nutrients that females can deposit into their eggs. Females can store and mobilize all fatty acids from internal stores to plasma (reviewed by [53]) and then to the yolk. However, since the PUFA and antioxidant composition of females is strongly linked to their diet (e.g., Parus major; [54]), variation in maternal food resources may directly affect hatchling and fledgling numbers.

Hatchling and fledgling numbers were highest in the nests of great tit mothers that laid eggs whose overall concentrations of SFAs, MUFAs and ω − 3 PUFAs were high (Fig. 1). SFAs and MUFAs, which are the most abundant fatty acids in the yolk (e.g., [10, 32]), can be endogenously produced (i.e., they are not essential fatty acids; reviewed by [27]), and represent an important source of energy. ω − 3 PUFAs are particularly important in early life stages when nestlings are developing and growing fast, and can enhance bone formation and immune competence (in birds e.g., [55,56,57]). Hence, eggs containing high concentrations of SFAs, MUFAs and ω − 3 PUFAs might have provided the embryo with both the energy and essential components that enhance development, and thereby may have boosted the number of chicks that hatched in this study. These fatty acids may have promoted early growth, as nestlings in our study that emerged from nests that had such eggs had longer tarsi on day 6 (Fig. 2a, Additional file 3). In our population, long tarsi are a good proxy for survival during the nestling phase (see Additional file 5), thus probably also contributing to increased fledgling numbers. Alternatively, eggs containing high concentrations of SFAs, MUFAs and ω − 3 PUFAs could have also had an overall high lipid content (reviewed by [26]). Therefore, in the future it will be important to test whether hatchling and fledgling numbers are higher in eggs that contain high concentrations of specific fatty acid groups or have an overall high lipid content.

By contrast, high concentrations of one class of fatty acids, ω − 6 PUFAs, and that of the antioxidant vitamin E were negatively related to fledgling number (Fig. 1). We offer two non-mutually exclusive explanations for this finding. Firstly, in great tits one would expect to observe an inverse relationship between the two types of ω-PUFAs in the yolk given that ω − 3 PUFAs and ω − 6 PUFAs are obtained from different food sources (see above; [32]). Hence, nestlings developing from nests with eggs containing high concentrations of ω − 6 PUFAs and vitamin E might have had lower concentrations of other resources that might have boosted fledging success. Secondly, poultry studies indicate that a ratio of high ω − 6 PUFA concentrations relative to those of ω − 3 PUFAs is associated with the formation of pro-inflammatory eicosanoids (i.e., mediators of inflammation), whilst the inverse ratio is associated with the formation of anti-inflammatory eicosanoids (e.g., [56, 57]). If chicks needed to mount inflammatory processes during the nestling phase, those chicks developing from clutches with a ratio of high ω − 6 versus ω − 3 PUFAs in their yolk might have shown strong inflammatory responses that could have lowered their fledging success.

In the current study, offspring fitness proxies were not related to yolk androgen and carotenoid concentrations, in contrast to several publications showing an impact of these yolk components on early development and survival in birds (e.g., [17, 18] and reviews by [19, 22]). Here we offer three non-mutually exclusive explanations for this result, which could be tested in subsequent experiments. On the one hand, and unlike many previous studies that focused on one group of yolk components (e.g., androgens, e.g., [17]), we simultaneously measured three different groups of yolk components, and some, like fatty acids, explained more variation in offspring phenotype and fitness than androgens and carotenoids. Thus, a potentially small effect of androgens and/or carotenoids might have been masked by a stronger influence of fatty acids. On the other hand, our study design imposed two limitations. Firstly, we analyzed the relationship between fitness and androgens/carotenoids by using PC3. However, this PC3 also included two MUFAs, which might have obscured the link with androgens and carotenoids—although we consider this possibility unlikely. Secondly, to be able to assess offspring fitness and phenotypes, we collected the middle egg as a representative for each clutch and allowed its siblings to develop normally. This technique is commonly used to study the adaptive value of maternal yolk deposition (in this species see [58, 59]; see also e.g., Ficedula albicollis [18, 60]) because alternative methods, like taking a biopsy of the yolk from each egg, can result in a substantial reduction of offspring fitness (e.g., hatching success; Troglodytes aedon, [61]) while effects on offspring phenotypes are unclear. In our population of free-living great tits, the concentrations of many yolk components exhibit a considerable repeatability (androgens R = 0.30–0.64 and carotenoids R = 0.38–0.39; [10]), suggesting that females are relatively consistent in their deposition of yolk components within a clutch (see also [62]). Nevertheless, despite the repeatability there also remains a considerable amount of residual variance (~ 36–70%; (1 − R)), which likely includes plasticity along the laying sequence in the concentrations of yolk components that mothers deposit [10]. For those reasons, the lack of an association between yolk androgen and carotenoid concentrations with fitness traits requires further investigation in a multi-component study.

Egg components and phenotypic traits

Our study indicates that maternal yolk deposition also shapes offspring phenotypic variation, but that its influence varies across traits and developmental stages (Fig. 2a, b). In particular, yolk substances showed a stronger link with physiological traits considered to be indicators of nestling health (oxidative status) than with morphological traits. In fact, only one morphological trait (tarsus length) was explained by yolk content and only on day 6, before nestlings reached exponential growth.

Fatty acids, antioxidants and androgens were all associated with the oxidative status of nestlings, albeit via different components of the redox system (Fig. 2a, b). High yolk concentrations of ω − 6 PUFA and vitamin E, which are components negatively related to fledgling number, were associated with a low oxidative status of nestlings on day 12 (i.e., nestlings had low enzymatic antioxidant concentrations in red blood cells and high concentrations of oxidative damage in plasma). Conversely, high concentrations of SFAs, MUFAs and ω − 3 PUFAs, yolk components positively associated with hatchling and fledgling numbers, were linked to low concentrations of non-enzymatic antioxidants in plasma of nestlings on day 12 (but were not associated with nestling oxidative damage). Lastly, high yolk concentrations of androgens and carotenoids were also related to low concentrations of non-enzymatic antioxidants in nestlings on days 6 and 12. These results lead us to speculate that the concentrations of yolk fatty acids influence the oxidative status of the nestlings on day 12, thus also explaining the differences observed in fledgling numbers. However, this seems not to be the case in our population given that the majority of the nestlings recorded on day 12 (N = 160) survived to day 15 (N = 153).

In our study, maternal yolk deposition predicted fewer aspects of the nestlings’ oxidative status on day 6 compared to day 12 (Fig. 2). These results are in contradiction with the hypothesis that the strength of maternal effects declines throughout offspring development (e.g., [63]). However, our results could be explained if the development of the great tit nestling antioxidant system occurs during similar ontogenetic times as in poultry chicks. In poultry, carotenoids and vitamin E are transferred from the egg yolk to the chick’s tissues during the last week of incubation. Soon after hatching, nestlings are not yet able to effectively assimilate antioxidants from the diet and their antioxidant system is mainly determined by the concentration of yolk antioxidants deposited by the mother (reviewed by [64]). One to two weeks after hatching, a chick’s plasma concentrations of carotenoids and vitamin E decrease dramatically (reviewed by [64, 65]), and subsequently the antioxidant protection of a nestling is increasingly determined by its own production of enzymatic antioxidants like GPX (reviewed by [65]). In free-living populations, the temporal dynamics of chick antioxidant protection provided by maternal yolk resources versus the chick’s own antioxidant production remains to be studied. But if a similar temporal pattern of antioxidant protection as in poultry chicks occurs in free-living great tit nestlings, both a decrease in maternally derived antioxidants and an increase in metabolic demands during the nestling period could explain the stronger relationship between maternal yolk hormones and the nestlings’ oxidative status observed on day 12 compared to day 6. This could especially be the case for nestlings developing from clutches where the fourth egg had high concentrations of ω − 6 PUFA and vitamin E, since the former component is susceptible to damage from lipid peroxidation by reactive oxygen species that are generated as by-products of offspring metabolism.


Our study suggests a novel link between natural variation in yolk fatty acid content and offspring fitness and phenotype in a wild bird population. It also indicates that groupings of yolk components can be related to fitness and phenotypic traits of great tit nestlings in different ways. This finding, together with previous studies showing that substances are co-secreted into the yolk (e.g., [66]), supports the idea that the consequences of maternal yolk deposition for the offspring should be addressed by simultaneously studying several components rather than single yolk components [34]. Our study represents the first step towards an integrative understanding of how steroid hormones, antioxidants and fatty acids may jointly shape offspring development. Carefully designed experimental studies, in which egg components are being manipulated both in isolation and simultaneously, are now needed to test for the causal effect that each component—as well as in interaction with other yolk substances—have on the offspring.

Variation in offspring fitness proxies was explained by some yolk components that mothers have to acquire from their diet (i.e., vitamin E and both PUFA classes). This suggests that offspring fitness also depends on the female’s access to, or ability to exploit, high-quality food sources [30]. A female’s diet can be influenced by both genetic and environmental factors and it will now be important to determine the relative contributions of these two variables to yolk composition. Further, if hatchling or fledgling numbers were mainly related to concentrations of essential (dietary) yolk components, offspring fitness would greatly rely on the ability of females to find food. In this scenario fluctuations in food abundance, either because of natural variation in weather and habitat quality or human modification of the environment, could have direct fitness consequences for the offspring. Alternatively, if fitness traits were more strongly influenced by non-essential yolk components (i.e., SFA or MUFAs), offspring fitness would depend on the physiological capacity of females to produce and transfer these components into the yolk. Also, it is worth noting that in this study we cannot fully exclude the possibility that nestling fitness proxies and phenotypes could result from differences in habitat quality among breeding pairs that influenced the quality of the food delivered by the parents to their young. Future experimental studies should test to what extent differences in nestling fitness proxies and phenotypes are related solely to maternal yolk components versus nestling food quality.

By collecting the fourth egg, we determined the relationship between the average deposition of yolk components and offspring fitness and phenotype. Future studies should investigate how other sources of variation in female yolk deposition (i.e., the plasticity of deposition along the laying sequence and the covariance between average deposition and plasticity; e.g., [10, 67]) shape offspring traits.


Study species and field site

We studied a nest-box population of great tits from April to July in both 2015 and 2016 in the Dellinger Buchet, in Southern Germany (Bavaria; 48° 03′ N, 11° 13′ E, 620 m above sea level), a forest with a mosaic of deciduous and coniferous trees. Female great tits generally lay one egg each day. In our population, mean clutch size is (mean ± SD) 8.45 ± 1.13 eggs, incubation starts after the last egg is laid and lasts approximately 14 days (mean ± SD = 14.07 ± 2.77), and nestlings spend around 20 days in the nest (mean ± SD = 20.60 ± 1.43; our own unpubl. data). In this study, we only considered first clutches.

Nest monitoring and egg collection

We visited nests every second day from the beginning of the breeding season onwards. Once egg laying started, we visited nests every day and marked eggs with a pencil to identify the laying order. Great tit females vary substantially in the mean allocation of yolk components and in their plasticity along the laying sequence (e.g., [10, 62]); yet, in our population, we previously reported that yolk components show medium repeatabilities for androgens (R = 0.30–0.64) and antioxidants (R = 0.36–0.39), and medium–high repeatabilities for fatty acids (R > 0.5 for SFAs and MUFAs and R > 0.9 for PUFAs; [10]). These repeatability estimates suggest that the concentrations of yolk components in eggs laid by the same female are more similar to each other than to eggs laid by other females. These medium–high repeatabilities also suggest that the middle egg of each clutch represents the average yolk content for a given clutch. In the present study, we therefore collected the fourth egg because it represents, on average, the middle egg of a great tit clutch (for a similar approach in this species see [58, 59]. We collected eggs between 8:00 and 13:00 h on the day they were laid and replaced them with a dummy egg. All collected eggs from a total of 69 clutches were laid by different females, except for two focal females that were included in both years. On the day of collection, we weighed and opened each freshly laid egg in the laboratory and separated the yolk from the albumen by rolling it on a piece of paper following [68]. We then homogenized the yolk by mixing it with an equal amount of distilled water (1 µl per mg of yolk) and stored it at − 80 °C until further analysis.

Nestling growth monitoring and blood sampling

We checked nests until clutch completion, as identified by the absence of freshly laid eggs and incubation behaviour of females. Two days prior to expected hatching, we monitored each nest closely to record the hatching date (day 0 = the day the first hatchling was observed). On day 1, we identified each nestling individually by clipping several down feathers, and we recorded its body mass (to the nearest 0.1 g). On day 6 or 7 (hereafter referred to as day 6) and on day 12 or 13 after hatching (hereafter referred to as day 12), we measured nestlings’ oxidative status. We selected these two periods because they are just before and after great tit nestlings reach exponential growth (i.e., on day 9–10 post-hatching; [69]), and because nestlings are large enough to collect the minimum blood volumes required to analyse their oxidative status. In particular, on day 6 we collected a small blood sample (~ 20 μl) from the brachial vein with heparinized capillaries to determine nestling oxidative status. We sampled all nestlings (N = 182) from the same nest (N = 51) within 15 min of disturbance (mean ± SD = 12.72 ± 6.93), again recording body mass and tarsus length (to the nearest 0.1 mm). On day 12, we collected the second blood sample (~ 80 μl) to determine each nestling’s oxidative status. We took samples within 3 min of disturbance (mean ± SD = 1.86 ± 0.73) from two or three randomly selected chicks (N = 96) per nest (N = 36), and recorded body mass and tarsus length. We then fitted each nestling with a numbered aluminum ring. On day 15, we recorded final body mass and tarsus length, and fledging was monitored until all young had left the nest. Finally, we also recorded the environmental conditions experienced by great tits during the incubation and nestling periods (see Additional file 6).

Yolk analyses

In total, we measured 31 yolk components (Additional file 2). The equipment, methods and further assay details used to measure each yolk component are described in the Additional file 6 (see also [10]). Briefly, we separated steroid hormones (androstenedione, 5α-dihydrotestosterone, testosterone, and corticosterone) via diatomaceous earth columns. For this, we followed the method described by [70], modified by [71] with additional adjustments for the measurement of egg yolk following [72]. We quantified androgen concentrations using radio-immunoassays, assaying each sample in duplicate and correcting hormone concentration of each sample for its’ individual extraction efficiency. All samples were analysed in two assays. Intra-assay coefficients of variation, as determined from the positive controls containing stripped chicken plasma with a known quantity of hormone added, were: androstenedione = 9.6%, 5α-dihydrotestosterone = 14.9% and testosterone = 23.8%. The inter-assay coefficients of variation, as determined from the first positive control of each assay, were: androstenedione = 13.3%, 5α-dihydrotestosterone = 13.6% and testosterone = 22.3%. We determined corticosterone concentrations using enzyme immunoassays (lot numbers: 12041402D and 04281702, Enzo Life Sciences, Germany). We assayed each sample in duplicate and distributed samples across five assays. The intra-assay coefficients of variation, also determined from stripped chicken plasma with a known quantity of corticosterone added, were 6.2, 11.4, 19, 7 and 12.3%, and the inter-assay coefficient of variation was 10.8%.

We extracted and then quantified antioxidants (lutein, zeaxanthin and vitamin E) by high-performance liquid chromatography (HPLC) following [10] (for details see Additional file 6). We calculated antioxidant concentrations from standard curves made for lutein, zeaxanthin and vitamin E (α-tocopherol) along with corrections for their respective internal standards. Of the 69 eggs analysed, data on the concentration of lutein and zeaxanthin were missing for one and two eggs, respectively. We assigned the value of the average population of each antioxidant to those eggs with missing values [73].

We extracted and subsequently analysed fatty acids (saturated fatty acids: SFAs; mono-unsaturated fatty acids: MUFAs; ω − 3 and ω − 6 polyunsaturated fatty acids: PUFAs) using gas chromatography-mass spectrometry (GC–MS), according to previously established methods [10, 74] (for details see also Additional file 6). We quantified fatty acid concentrations based on the peak area of the internal standard (methyl cis-10-heptadecenoate; 16.65 µg added to each sample).

Oxidative biomarkers measurements

On days 6 and 12, we determined the concentrations of antioxidants and oxidative damage in nestling plasma (details on the equipment and methods used are described in Additional file 6). Briefly, we measured OXY (expressed as mM HOCl neutralized) using the OXY-Adsorbent test (Diacron International SRL, Grosseto, Italy) following [75], and oxidative damage (ROMs: produced by the oxidation of lipids, proteins and nucleic acids, and are expressed as mM H2O2 equivalents) using the d-ROMs test kit (Diacron International SRL, Grosseto, Italy) following [75]. We assayed each sample in duplicate and samples were analysed across 21 plates. Inter-assay coefficients of variation, as determined from the calibrators (OXY) or from known standards (ROMs), were 5.1% and 2.8% for OXY and ROMs respectively. On day 12, we additionally measured the enzymatic activity of GPX (expressed as U/ml) in red blood cells using the Ransel assay (Randox Laboratories, Germany) following [76].

Statistical analysis

To group yolk components, we ran a principal component analysis (PCA), including all 31 yolk components analysed (for further details on this analysis see Additional file 6). The three principal components (PC) identified accounted for 58% of the cumulative variance in our data (Additional file 1). The mean values of the principal components differed between years (Additional file 7). To account for year differences, we standardized the loadings because ‘year’ could not be included as a random (we only had two levels) or a fixed factor (limited sample size). Thus, for each nest we first subtracted and then divided the loading of each principal component by the mean value of that component in the corresponding year. This method allowed us to assess the importance of each yolk component irrespective of year differences [77].

To study the relationship between yolk components in the fourth egg of clutches and fitness proxies for these clutches, we first ran two generalized linear models (we initially ran two generalized linear mixed effect models by including ‘female identity’ as a random factor, but since this factor did not explain any variance, we excluded it from all final models). Hatchling and fledgling number were fitted as response variables, PC1, PC2, PC3, date of egg collection and clutch size were included as covariates (two clutches were excluded from analyses; see Additional file 6). Body mass at fledging is a key predictor of post-fledging survival in great tits [78]. As other proxies of fitness, we therefore assessed the relationship between yolk components and nestling mass and tarsus length on day 15. These morphological variables were studied in separate models because yolk components can influence structural body size (i.e., tarsus length) independently of mass (e.g., [55]). Because of sample size constraints and to reduce the number of explanatory variables, we fitted the residuals of a linear regression between these two estimates of body condition and clutch size as response variables in two separate linear mixed-effect models. PC1, PC2, PC3 and capture date were included as covariates in the model, and nest ID was fitted as a random factor. Environmental variables were initially included in the analyses, but since they did not influence the response variable we excluded them from all final models (for details see Additional file 6).

To study the relationship between egg components and the physiological condition of individual nestlings from a given brood, we ran linear mixed-effect models. We fitted OXY, GPX (only for day 12), ROMs, nestling mass and tarsus length (corrected for brood size) as response variables. We included PC1, PC2, PC3 and date of egg collection as covariates, and nest identification as a random factor. We initially included other covariates in the model based on their biological relevance to the study question (see Additional file 6), but because of sample size constraints, only variables that had an effect on the response variables were retained in the final model. In particular, clutch size was retained as a covariate in models analysing OXY concentrations, and on day 12, total sampling time was fitted as a covariate in the model for ROM concentrations.

All statistical analyses were performed in R statistical freeware R-3.3.3 [79]. Principal component analyses were performed using the ‘prcomp’ package. Statistical models were conducted using the ‘lme4’ and ‘arm’ packages in a Bayesian framework with non-informative priors. We assumed a Poisson error distribution for the generalized linear models and a Gaussian error distribution for the linear mixed-effect models. In all cases, residuals were checked visually for the model fit. Whenever necessary, response variables were transformed (details on transformations are provided in the tables). We mean-centered all covariates (i.e., mean value = 0, standard deviation = 1) because covariates differed in their scales of magnitude. We based model structure on the study question and the biology of the species rather than on model selection [80]. We subsequently used the ‘sim’ function to simulate values from the posterior distributions of model parameters. From 10,000 simulations, we extracted the 95% Bayesian credible interval (CrI) around the mean [77] and assessed statistical support by obtaining the posterior distribution of each parameter. CrIs provide more valuable information than p-values, such as the uncertainty around the estimates. We use the term ‘statistically meaningful’ when the estimated effect differed from zero with a posterior probability higher than 0.95. A threshold of 5% is equivalent to the significance level in a frequentist framework (for further details on statistical inference, see [80].

Availability of data and materials

The datasets and code supporting this article are available from the Dryad Digital Repository at


  1. 1.

    Moore MP, Whiteman HH, Martin RA. A mother’s legacy: the strength of maternal effects in animal populations. Ecol Lett. 2019;22:1620–8.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Sánchez-Tójar A, Lagisz M, Moran NP, Nakagawa S, Noble DWA, Reinhold K. The jury is still out regarding the generality of adaptive ‘transgenerational’ effects. Ecol Lett. 2020;23:1715–8.

    PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Yin J, Zhou M, Lin Z, Li QQ, Zhang Y. Transgenerational effects benefit offspring across diverse environments: a meta-analysis in plants and animals. Ecol Lett. 2019;22:1976–86.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Mousseau TA, Fox CW. The adaptive significance of maternal effects. Trends Ecol Evol. 1998;10:403–7.

    Article  Google Scholar 

  5. 5.

    Uller T. Developmental plasticity and the evolution of parental effects. Trends Ecol Evol. 2008;23:432–8.

    PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    Uller T. The evolution of parental care. In: Royle NJ, Smiseth PTKM, editors. Evol parent care. New York: Oxford University Press; 2012. p. 247–67.

    Google Scholar 

  7. 7.

    Fuiman LA. Egg boon fatty acids reveal effects of a climatic event on a marine food web. Ecol Monogr. 2018;88:585–99.

    Article  Google Scholar 

  8. 8.

    Lovern MB, Wade J. Maternal plasma and egg yolk testosterone concentrations during embryonic development in green anoles (Anolis carolinensis). Gen Comp Endocrinol. 2001;124:226–35.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  9. 9.

    McCormick MI. Experimental test of the effect of maternal hormones on larval quality of a coral reef fish. Oecologia. 1999;118:412–22.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Mentesana L, Isaksson C, Goymann W, Andersson MN, Trappschuh M, Hau M. Female variation in allocation of steroid hormones, antioxidants and fatty acids: a multilevel analysis in a wild passerine bird. J Avian Biol. 2019.

    Article  Google Scholar 

  11. 11.

    Noble RC, Deeming DC, Ferguson MWJ, McCartney R. Changes in the lipid and fatty acid composition of the yolk during embryonic development of the alligator (Alligator mississipiensis). Comp Biochem Physiol Part B Comp Biochem. 1990;96:183–7.

    Article  Google Scholar 

  12. 12.

    Saino N, Bertacche V, Ferrari RP, Martinelli R, Møller AP, Stradi R. Carotenoid concentration in barn swallow eggs is influenced by laying order, maternal infection and paternal ornamentation. Proc Biol Sci. 2002;269:1729–33.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Valcu CM, Scheltema RA, Schweiggert RM, Valcu M, Teltscher K, Walther DM, et al. Life history shapes variation in egg composition in the blue tit Cyanistes caeruleus. Commun Biol. 2019;2:1–14.

    Article  Google Scholar 

  14. 14.

    Weiss SL, Kennedy EA, Safran RJ, Mcgraw KJ. Pterin-based ornamental coloration predicts yolk antioxidant levels in female striped plateau lizards (Sceloporus virgatus). J Anim Ecol. 2011;80:519–27.

    PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Brown AC, Leonard HM, McGraw KJ, Clotfelter ED. Maternal effects of carotenoid supplementation in an ornamented cichlid fish. Funct Ecol. 2014;28:612–20.

    Article  Google Scholar 

  16. 16.

    Leone VA, Worzalla SP, Cook ME. Body compositional changes and growth alteration in chicks from hens fed conjugated linoleic acid. Lipids. 2009;44:437–47.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Hegyi G, Herényi M, Szöllősi E, Rosivall B, Török J, Groothuis TGG. Yolk androstenedione, but not testosterone, predicts offspring fate and reflects parental quality. Behav Ecol. 2011;22:29–38.

    Article  Google Scholar 

  18. 18.

    Tschirren B, Postma E, Gustafsson L, Groothuis TGG, Doligez B. Natural selection acts in opposite ways on correlated hormonal mediators of prenatal maternal effects in a wild bird population. Ecol Lett. 2014;17:1310–5.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    von Engelhardt N, Groothuis TGG. Maternal hormones in avian eggs. In: Norris DO, Lopez KH, editors. Hormones and reproduction in vertebrates, vol 4: Birds. 2011; p. 91–128.

  20. 20.

    Haussmann MF, Longenecker AS, Marchetto NM, Juliano SA, Bowden RM. Embryonic exposure to corticosterone modifies the juvenile stress response, oxidative stress and telomere length. Proc R Soc Lond B Biol Sci. 2012;279:1447–56.

    CAS  Google Scholar 

  21. 21.

    Treidel LA, Whitley BN, Benowitz-Fredericks ZM, Haussmann MF. Prenatal exposure to testosterone impairs oxidative damage repair efficiency in the domestic chicken (Gallus gallus). Biol Lett. 2013;9:20130684.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Deeming DC, Pike TW. Embryonic growth and antioxidant provision in avian eggs. Biol Lett. 2013;9:20130757.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Surai PF. Effect of selenium and vitamin E content of the maternal diet on the antioxidant system of the yolk and the developing chick. Br Poult Sci. 2000;41:235–43.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Parolini M, Khoriauli L, Possenti CD, Colombo G, Caprioli M, Santagostino M, et al. Yolk vitamin E prevents oxidative damage in gull hatchlings. R Soc Open Sci. 2017;4:170098.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  25. 25.

    Watson H, Salmón P, Isaksson C. Maternally derived yolk antioxidants buffer the developing avian embryo against oxidative stress induced by hyperoxia. J Exp Biol. 2018;221:jeb179465.

    PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Noble RC, Cocchi M. Lipid metabolism and the neonatal chicken. Prog Lipid Res. 1990;29:107–40.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    Surai PF, Speake BK. The natural fatty acid compositions of eggs of wild birds and the consequences of domestication. In: De Meester F, Watson RR, editors. Wild-type food in health promotion and disease prevention. Totowa: Humana Press; 2008. p. 121–37.

    Chapter  Google Scholar 

  28. 28.

    Hulbert AJ, Abbott SK. Nutritional ecology of essential fatty acids: an evolutionary perspective. Aust J Zool. 2011;59:369.

    Article  Google Scholar 

  29. 29.

    Bielski BHJ, Arudi RL, Sutherland MW. A study of the reactivity of HO2/O2- with unsaturated fatty acids. J Biol Chem. 1983;258:4759–61.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Twining CW, Brenna JT, Hairston NG, Flecker AS. Highly unsaturated fatty acids in nature: what we know and what we need to learn. Oikos. 2016;125:749–60.

    CAS  Article  Google Scholar 

  31. 31.

    Bourgault P, Thomas DW, Blondel J, Perret P, Lambrechts MM. Between-population differences in egg composition in Blue Tits (Cyanistes caeruleus). Can J Zool. 2007;85:71–80.

    CAS  Article  Google Scholar 

  32. 32.

    Toledo A, Andersson MN, Wang H-L, Salmón P, Watson H, Burdge GC, et al. Fatty acid profiles of great tit (Parus major) eggs differ between urban and rural habitats, but not between coniferous and deciduous forests. Sci Nat. 2016;103:55.

    CAS  Article  Google Scholar 

  33. 33.

    Groothuis TGG, Hsu B-Y, Kumar N, Tschirren B. Revisiting mechanisms and functions of prenatal hormone-mediated maternal effects using avian species as a model. Philos Trans R Soc B Biol Sci. 2019;374:20180115.

    CAS  Article  Google Scholar 

  34. 34.

    Williams TD, Groothuis T. Egg quality, embryonic development, and post-hatching phenotype: an integrated perspective. In: Deeming DC, James S, editors. Nests, eggs, and incubation: new ideas about avian reproduction. Oxford: Oxford University Press; 2015.

    Google Scholar 

  35. 35.

    Giraudeau M, Ziegler A-K, Pick JL, Ducatez S, Canale CI, Tschirren B. Interactive effects of yolk testosterone and carotenoid on prenatal growth and offspring physiology in a precocial bird. Behav Ecol. 2016.

    Article  Google Scholar 

  36. 36.

    Possenti CD, Secomandi S, Schiavon A, Caprioli M, Rubolini D, Romano A, et al. Independent and combined effects of egg pro- and anti-oxidants on gull chick phenotype. J Exp Biol. 2018;221:jeb174300.

    PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Parolini M, Possenti CD, Secomandi S, Carboni S, Caprioli M, Rubolini D, et al. Prenatal independent and combined effects of yolk vitamin E and corticosterone on embryo growth and oxidative status in the yellow-legged gull. J Exp Biol. 2019;222:jeb199265.

    PubMed  Article  Google Scholar 

  38. 38.

    Costantini D. Understanding diversity in oxidative status and oxidative stress: the opportunities and challenges ahead. J Exp Biol. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Bergamini C, Gambetti S, Dondi A, Cervellati C. Oxygen, reactive oxygen species and tissue damage. Curr Pharm Des. 2005;10:1611–26.

    Article  Google Scholar 

  40. 40.

    Vágási CI, Vincze O, Pătraș L, Osváth G, Pénzes J, Haussmann MF, et al. Longevity and life history coevolve with oxidative stress in birds. Funct Ecol. 2019;33:152–61.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Blount JD, Vitikainen EIK, Stott I, Cant MA. Oxidative shielding and the cost of reproduction. Biol Rev. 2016;91:483–97.

    Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Sebastiano M, Eens M, Abd Elgawad H, De TB, Lacoste V, Pineau K, et al. Oxidative stress biomarkers are associated with visible clinical signs of a disease in frigatebird nestlings. Sci Rep. 2017;7:1–13.

    CAS  Article  Google Scholar 

  43. 43.

    Noguera JC, Monaghan P, Metcalfe NB. Interactive effects of early and later nutritional conditions on the adult antioxidant defence system in zebra finches. J Exp Biol. 2015;218:2211–7.

    PubMed  PubMed Central  Google Scholar 

  44. 44.

    Xu XL, Ji WJ, Castell JD, O’Dor RK. Influence of dietary lipid sources on fecundity, egg hatchability and fatty acid composition of Chinese prawn (Penaeus chinensis) broodstock. Aquaculture. 1994;119:359–70.

    CAS  Article  Google Scholar 

  45. 45.

    Chen M, Liu H, Chen B. Effects of dietary essential fatty acids on reproduction rates of a subtropical calanoid copepod, Acartia erythraea. Mar Ecol Prog Ser. 2012;455:95–110.

    CAS  Article  Google Scholar 

  46. 46.

    Muma E, Palander S, Näsi M, Pfeiffer AM, Keller T, Griinari JM. Modulation of conjugated linoleic acid-induced loss of chicken egg hatchability by dietary soybean oil. Poult Sci. 2006;85:712–20.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. 47.

    Perrins CM. Tits and their caterpillar food supply. Ibis (Lond 1859). 1991;133:49–54.

    Article  Google Scholar 

  48. 48.

    Visser ME, Van Noordwijk AJ, Tinbergen JM, Lessells CM. Warmer springs lead to mistimed reproduction in great tits (Parus major). Proc R Soc B Biol Sci. 1998;265:1867–70.

    Article  Google Scholar 

  49. 49.

    Isaksson C, Andersson S. Carotenoid diet and nestling provisioning in urban and rural great tits Parus major. J Avian Biol. 2007;38:564–72.

    Article  Google Scholar 

  50. 50.

    Andersson MN, Wang H-L, Nord A, Salmón P, Isaksson C. Composition of physiologically important fatty acids in great tits differs between urban and rural populations on a seasonal basis. Front Ecol Evol. 2015;3:93.

    Article  Google Scholar 

  51. 51.

    Isaksson C, Hanson MA, Burdge GC. The effects of spatial and temporal ecological variation on fatty acid compositions of wild great tits Parus major. J Avian Biol. 2015;46:245–53.

    Article  Google Scholar 

  52. 52.

    Naef-Daenzer L, Naef-Daenzer B, Nager RG. Prey selection and foraging performance of breeding Great Tits Parus major in relation to food availability. J Avian Biol. 2000;31:206–14.

    Article  Google Scholar 

  53. 53.

    Raclot T. Selective mobilization of fatty acids from adipose tissue triacylglycerols. Prog Lipid Res. 2003;42:257–88.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Andersson MN, Nilsson J, Nilsson JÅ, Isaksson C. Diet and ambient temperature interact to shape plasma fatty acid composition, basal metabolic rate and oxidative stress in great tits. J Exp Biol. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Twining CW, Brenna JT, Lawrence P, Shipley JR, Tollefson TN, Winkler DW. Omega-3 long-chain polyunsaturated fatty acids support aerial insectivore performance more than food quantity. Proc Natl Acad Sci USA. 2016;113:10920–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  56. 56.

    Hall JA, Jha S, Skinner MM, Cherian G. Maternal dietary n-3 fatty acids alter immune cell fatty acid composition and leukotriene production in growing chicks. Prostaglandins Leukot Essent Fat Acids. 2007;76:19–28.

    CAS  Article  Google Scholar 

  57. 57.

    Cherian G, Bautista-Ortega J, Goeger DE. Maternal dietary n-3 fatty acids alter cardiac ventricle fatty acid composition, prostaglandin and thromboxane production in growing chicks. Prostaglandins Leukot Essent Fat Acids. 2009;80:297–303.

    CAS  Article  Google Scholar 

  58. 58.

    Giordano M, Groothuis TGG, Tschirren B. Interactions between prenatal maternal effects and posthatching conditions in a wild bird population. Behav Ecol. 2014;25:1459–66.

    Article  Google Scholar 

  59. 59.

    Remeš V, Matysioková B, Klejdus B. Egg yolk antioxidant deposition as a function of parental ornamentation, age, and environment in great tits Parus major. J Avian Biol. 2011;42:387–96.

    Article  Google Scholar 

  60. 60.

    Tschirren B, Sendecka J, Groothuis TGG, Gustafsson L, Doligez B. Heritable variation in maternal yolk hormone transfer in a wild bird population. Am Nat. 2009;174:557–64.

    Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Ellis LA, Borst DW, Thompson CF. Hatching asynchrony and maternal androgens in egg yolks of House Wrens. J Avian Biol. 2001;32:26–30.

    Article  Google Scholar 

  62. 62.

    Lessells CM, Ruuskanen S, Schwabl H. Yolk steroids in great tit Parus major eggs: variation and covariation between hormones and with environmental and parental factors. Behav Ecol Sociobiol. 2016;70:843–56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  63. 63.

    Bernardo J. Maternai effects in animal ecology. Am Zool. 1996;36:83–105.

    Article  Google Scholar 

  64. 64.

    Surai PF, Fisinin VI, Karadas F. Antioxidant systems in chick embryo development. Part 1. Vitamin E, carotenoids and selenium. Anim Nutr. 2016;2:1–11.

    Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Surai PF. Vitamin E in avian reproduction. Avian Poult Biol Rev. 1999;10:1–60.

    Google Scholar 

  66. 66.

    Giraudeau M, Ducatez S. Co-adjustment of yolk antioxidants and androgens in birds. Biol Lett. 2016;12:20160676.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  67. 67.

    Hsu B-Y, Verhagen I, Gienapp P, Darras VM, Visser ME, Ruuskanen S. Between- and within-individual variation of maternal thyroid hormone deposition in Wild Great Tits (Parus major). Am Nat. 2019;194:E96-108.

    Article  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Lessells CM, Dingemanse NJ, Both C. Egg weights, egg components weights, and laying gaps in great tits (Parus major) in relation to ambient temperature. Auk. 2002;119:1091.

    Article  Google Scholar 

  69. 69.

    Gibb J. The Breeding Biology of the Great and Blue Titmice. Ibis (Lond 1859). 1950;111.

  70. 70.

    Wingfield JC, Farner DS. Avian endocrinology: field investigations and methods. Condor. 1976;78:570.

    Article  Google Scholar 

  71. 71.

    Goymann W, Wittenzellner A, Schwabl I, Makomba M. Progesterone modulates aggression in sex-role reversed female African black coucals. Proc R Soc Lond B Biol Sci. 2008;275:1053–60.

    CAS  Google Scholar 

  72. 72.

    Schwabl H. Yolk is a source of maternal testosterone for developing birds. Proc Natl Acad Sci. 1993;90:11446–50.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  73. 73.

    Nakagawa S, Freckleton RP. Missing inaction: the dangers of ignoring missing data. Trends Ecol Evol. 2008;23:592–6.

    PubMed  Article  PubMed Central  Google Scholar 

  74. 74.

    Eikenaar C, Källstig E, Andersson MN, Herrera-Dueñas A, Isaksson C. Oxidative challenges of avian migration: a comparative field study on a partial migrant. Physiol Biochem Zool. 2016;90:223–9.

    PubMed  Article  PubMed Central  Google Scholar 

  75. 75.

    Costantini D, Casagrande S, De Filippis S, Brambilla G, Fanfani A, Tagliavini J, et al. Correlates of oxidative stress in wild kestrel nestlings (Falco tinnunculus). J Comp Physiol B. 2006;176:329–37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  76. 76.

    Costantini D, Monaghan P, Metcalfe NB. Biochemical integration of blood redox state in captive zebra finches (Taeniopygia guttata). J Exp Biol. 2011;214:1148–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  77. 77.

    Gelman A, Hill J. Data analysis using regression and multilevel/hierarchical models. Cambridge: Cambridge University Press; 2007.

    Google Scholar 

  78. 78.

    Tinbergen JM, Boerlijst MC. Nestling weight and survival in individual Great Tits (Parus major). J Anim Ecol. 1990;59:1113.

    Article  Google Scholar 

  79. 79.

    R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.

    Google Scholar 

  80. 80.

    Korner-Nievergelt F, Roth T, von Felten S, Guélat J, Almasi B, Korner-Nievergelt P. Bayesian data analysis in ecology using linear models with R, BUGS, and Stan. Amsterdam: Elsevier; 2015.

    Google Scholar 

Download references


We are grateful to Nicolás Adreani, Sabine Jörg, and Natalia Pérez-Ruiz for their invaluable help in the field, Caroline Deimel for logistical help, and Esteban Botero-Delgadillo for his statistical advice. Thanks also to Monika Trappschuh for her help during egg and steroid analysis, and to Hong-Lei Wang and Jürgen Kuhn for their contributions to fatty acid analyses and antioxidant quantification, respectively. Finally, LM wants to thank Nico and Inti for their support during the writing process.


Open Access funding enabled and organized by Projekt DEAL. The study was funded by the Max Planck Society (to MH). LM was supported by the International Max Planck Research School (IMPRS) for Organismal Biology.

Author information




LM and MH designed the study. LM conducted field work, analysed the data and drafted the manuscript. LM conducted egg and steroid analysis under the guidance of WG. LM extracted and analysed egg antioxidants and fatty acids under supervision from MNA and CI. LM and SC conducted the oxidative stress analyses. All authors contributed to manuscript preparation and approved the final version of the manuscript.

Corresponding author

Correspondence to Lucia Mentesana.

Ethics declarations

Ethics approval and consent to participate

All experimental procedures were conducted according to the legal requirements in Germany and were approved by the governmental authorities of Oberbayern, Germany.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Martin N. Andersson, Stefania Casagrande, Wolfgang Goymann, Caroline Isaksson ordered alphabetically.

Supplementary Information

Additional file 1.

Principal component analysis of 31 yolk components measured in the fourth egg of 69 wild great tit clutches.

Additional file 2.

Results of linear models and linear mixed-effects models to test for relationships between yolk composition and fitness traits.

Additional file 3.

 Results of linear mixed-effects models to test for relationships between yolk composition and phenotypic traits on nestling day 6.

Additional file 4.

Results of linear mixed-effects models to test for relationships between yolk composition and phenotypic traits on day 12.

Additional file 5.

Results of linear mixed-effects models used to test for relationships between nestling tarsus length on days 6 and 12 and survival in the nest.

Additional file 6.


Additional file 7.

Linear models testing for differences in yolk composition between years.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mentesana, L., Andersson, M.N., Casagrande, S. et al. Natural variation in yolk fatty acids, but not androgens, predicts offspring fitness in a wild bird. Front Zool 18, 38 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Maternal effects
  • Fitness
  • Phenotypic variance
  • Steroid hormones
  • Antioxidants
  • Fatty acids