Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Lay Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (7)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Tinbergen, J. M.
Right arrow Articles by Sanz, J. J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Tinbergen, J. M.
Right arrow Articles by Sanz, J. J.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Behavioral Ecology Vol. 15 No. 4: 525-533
Behavioral Ecology vol. 15 no. 4 © International Society for Behavioral Ecology 2004; all rights reserved

Strong evidence for selection for larger brood size in a great tit population

Joost M. Tinbergena and Juan J. Sanzb

a Group Animal Ecology, Center for Ecological and Evolutionary Studies, University of Groningen, P.O. Box 14, 9750 AA, Haren, The Netherlands, and b Departemento de Ecología Evolutiva, Museo Nacional de Ciencias Naturales (CSIC), J Gutierrez Abascal 2, E-28006 Madrid, Spain

Address correspondence to J. M. Tinbergen. E-mail: j.m.tinbergen{at}biol.rug.nl

Received 7 August 2002; revised 1 May 2003; accepted 17 May 2003.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We measured the selection pressure on brood size in a recently established population of great tits (Parus major L.) in the northern Netherlands by manipulating brood size in three years (1995: n = 51, 1997: n = 66, 1998: n = 51), and we estimated fitness consequences in terms of local survival of both offspring and parents. Enlarged broods had highest fitness; the offspring fitness component was positively affected by manipulation and the parental fitness component was unaffected. Parental survival and the probability that parents produced a second clutch were not affected by the treatment. However, parents that had raised enlarged broods produced their second clutch later in the season. Clutch size, brood size, and laying date of birds recaptured in the next breeding season were largely independent of the treatment. We conclude that there is strong evidence for selection for larger brood size and reject the individual optimization hypothesis for this population because the number of young in the nest predicts fitness independently of the manipulation history.

Key words: brood size, clutch size, fitness, local adaptation, manipulation, Parus major.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
With the goal of understanding the adaptive significance of clutch size variation, many studies have manipulated brood size and measured subsequent fitness consequences. The logic is that experimental changes in brood size can be viewed as small mutations. By measuring the fitness effects we can examine whether such a mutation will increase in frequency in the population. The assumptions are: (1) parents react on the manipulation as if the resulting brood size was their own ‘decision’ (Lessells, 1993Go), and (2) selection pressure works primarily during the nestling phase and not during egg laying or incubation (Monaghan and Nager, 1997Go).

Results of brood size manipulation studies showed that both enlarging and reducing brood size may lower fitness, suggesting either that selection pressures stabilize brood size (collared flycatcher Ficedula albicollis: Gustafsson and Sutherland, 1988Go; great tit: Lindén, 1990Go; Tinbergen and Daan, 1990Go) or that birds are able to produce a brood size that maximizes their individual fitness (individual optimization of brood size, great tit: Perrins and Moss, 1975Go; Pettifor et al., 1988Go; kestrel (Falco tinnunculus ): Daan et al., 1990Go; blue tit (Parus caeruleus): Pettifor, 1993aGo,bGo; great tit: Pettifor et al., 2001Go; but see Tinbergen and Both, 1999Go). If individual birds adjust their clutch size to local circumstances, we need an explanation for the fact that different populations of great tits living under different ecological circumstances lay clutches of similar size (van Balen, 1973Go). Perhaps tits distribute themselves over the available habitat in such a way that resources per pair are relatively constant between habitats (ideal free distribution). Recently, Both et al. (2000)Go have presented evidence that optimal clutch size is negatively affected by population density and that birds adjust their clutch size accordingly, suggesting that clutch size is adjusted to both food resources and competition for these resources.

An alternative explanation for the small between-population variation in clutch size is that selection may not work at such a fine local scale, causing clutch size to be non-adaptive in some localities (Dhondt et al., 1990Go; Dias and Blondel, 1996Go). The observed clutch size is then viewed as a consequence of selection pressures acting over more generations and habitats; gene flow may prevent local adaptation (Kawecki, 1995Go; McNamara, 1997Go). In this case we would expect positive, negative, or stabilizing selection in different populations.

Actually, the first thorough brood size experiments led Nur (1984aGo,bGo, 1986Go; but see Pettifor 1993aGo,bGo) to conclude that blue tits do not adjust their clutch size to the local situation; larger broods than the population mean led to higher fitness. This idea has won impetus by the work of Verhulst (1995)Go, who showed experimentally that a reduction in brood size enhanced the lifetime reproductive success in an island population of great tits, thus revealing a selection pressure for smaller clutches. On the basis of brood size experiments, Blondel et al. (1998)Go also found indications for selection for smaller clutches in the Corsican population of blue tits. Both Verhulst (1995)Go and Blondel et al. (1998)Go suggested that in these populations gene flow might prevent local adaptation (see also Dias and Blondel, 1996Go).

In this paper we present the results of a brood size experiment in a recently established population of great tits. We find strong evidence that selection works towards larger broods and then review the brood size manipulation studies in great tits to place these contrasting findings in perspective. Consequences for microevolution are discussed briefly.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
General
The study area was a young (planted in 1974–1975), mixed deciduous forest in the Lauwersmeer area (The Netherlands, 53° 20' N, 06° 12' E). The study was started in 1993 in part of the area that since 1980 had a limited number of nest-boxes (80). From 1994 onwards, around 200 nest-boxes were available in seven plots of differing size (6–106 ha).

Nest-boxes were checked for occupation by great tits at least once a week and the number of eggs, occurrence of incubation, and the number of young were recorded. Laying date was estimated from the assumption that one egg per day was laid. We made daily nest visits around the expected day of hatching. Hatching date is the date on which we recorded the first egg in the clutch to hatch. Adults were captured with spring traps when the young were 7–10 days old (day of hatching = 0) and were ringed, weighed, and measured. At this age the young were ringed as well.

Brood size manipulation
During three years (1995, 1997, 1998), first broods with a similar hatching date (maximum difference of one day), clutch size (maximum difference one egg), and the same number of hatchlings were randomly assigned to one of the groups: reduced, control, or enlarged. Manipulation of brood size was carried out on day 2 after the nestlings hatched (see Sanz and Tinbergen, 1999Go). On this day three nestlings were transferred quickly between nests to create reduced and enlarged broods of the same age. For each experimental pair of nests, there was a control brood in which brood size was not altered but the nest was subjected to the same disturbance as the reduced and enlarged broods.

The probability of having a second clutch was determined by reading color rings of the females during incubation of the second clutch or by catching the parents with a spring trap in the nest-box during the nestling phase.

Before manipulation the birds did not differ in clutch or brood size between manipulation categories or years (ANOVA, all p >.1; Table 1). Laying date did differ between years (ANOVA, F2,163 = 106.0, p <.001) but not within years between manipulation categories (ANOVA, F2,163 = 0.71, p >.9). The group of birds used for the different manipulations can therefore be considered homogeneous within years.


View this table:
[in this window]
[in a new window]
 
Table 1 Basic parameters for the manipulated broods.

 
Offspring fitness
We analyzed first brood offspring fitness on a per nest basis in three components: (1) the probability that a brood was successful (at least one chick fledges), (2) the probability that a chick survived from manipulation to fledging from a successful brood, and (3) the probability to recruit when fledged. These components were integrated in the number of recruits per first brood and compared between manipulation categories.

Survival analysis adults
Local survival was estimated on the basis of recaptures of the breeding birds during the next breeding season in the study area. Recapture rate was estimated using the program MARK v 2.1 (White and Burnham, 1999Go). Experimental birds were included from their experimental year onwards. The animals that had their broods manipulated in more than one year (31 of the 326) were included in the analysis. They were coded as a new individual from the second experiment onwards. Both survival and recapture rates over the first and subsequent winters were estimated as a function of sex, manipulation (together six groups), and time (differences between years). Recapture rate was modeled before survival. Model selection was performed by backward deletion of higher order terms. Because we expected survival effects of manipulation to occur in the first year after manipulation, we specifically tested effects of manipulation on the survivals in these cells (three per group) and estimated the survival in the subsequent years with one sex-specific parameter. The most complicated model included sex, manipulation, and time, but we found no justification to include any of these variables in the model. We used Aikaike's Information Criterion (AIC) to select the most parsimonious model (White and Burnham, 1999Go).

Fitness calculations
Estimating overall fitness involves both the parental and the offspring components weighed for the age at which reproduction takes place (reproductive value; Fisher, 1930Go). As no effects of manipulation on reproductive parameters in the subsequent year could be detected, fitness was approximated by the sum of parents and their offspring identified as a breeding bird in the next year. Manipulated broods with unidentified males were excluded, reducing the sample size from 168 to 158. The resulting estimate is slightly biased, because in case of nest failure it is more likely that the identity of the males remains unknown. Because the identity of all females was known, we used the sum of the locally survived females identified in the next year (female parent + female offspring) as a second fitness estimate and compared it with an estimate for the males (male parent + male offspring) using Poisson regression.

In addition, we calculated a formal fitness estimate {lambda} (rate of increase) for the manipulation categories separately and the reproductive values of the ‘mutants’ reduced and enlarged in a population of control birds (V[x, x*]) to value the fitness differences between them. The reason to calculate {lambda} was to use a well defined fitness estimate that allows correction for recapture probability and enables calculation of V(x, x*).The estimate of {lambda} was based on P (average annual survival of parents per manipulation category) and F (annual fecundity in terms of recruits from first and second broods divided by two). We assumed age independent survival and fecundity, and no maximum age. In such a simple case, {lambda} can be calculated using the Euler-Lotka equation:


For P < l and n = {infty} (the maximum age) this reduces to:


where {lambda} is the rate of increase of the population, P is the annual survival, and F is the annual reproduction (see also Charnov and Krebs, 1974Go).

To see whether a ‘mutant’ (reduced or enlarged broods) would spread in the natural population (control group) we calculated the reproductive value of the mutant in the population of residents (control group) using:


For P(x) < {lambda}(x*) and n = {infty} and using Equation 2 this reduces to:


where x is the manipulation category, {lambda}(x*) is {lambda} of the resident population (control group), P(x*) and F(x*) are the annual survival and fecundity of the resident population and P(x) and F(x) are the annual survival and fecundity of the mutant population (reduced or enlarged), respectively. A mutant will spread when V(x, x*) > 1.

Detection of individual optimization
If individual birds maximize fitness by adjusting the clutch size to their local circumstances (individual optimization hypothesis) we expect fitness to peak at original brood or clutch size. Because, in this study population, there were no significant indications of diminishing returns with manipulated brood size, as expected under the individual optimization hypothesis, it is of interest to analyze the fitness effects of manipulation in more detail. This can be done in different ways. We did it by modeling total fitness as a function of the actual number of young after manipulation (Y) and, to allow for diminishing returns, its square (Y2). Such a model could explain fitness effects of brood size including diminishing returns if only the number of young determined fitness. When individual optimization occurs, one would expect that adding a factor for manipulation (M three levels) and perhaps the interaction M*Y (or M*Y2) to a Poisson model containing Y and Y2 should explain additional variance. If adding M would reduce deviance significantly, it means that fitness depends not only on the number of young actually in the nest, but also on the manipulation. This may be the number originally born or some correlate thereof, such as clutch size or local resources. One can also test this idea by including original brood size or clutch size in the analysis. Significant reductions in deviance again point in the direction that some pairs (or their territories) are ‘better’ than others. This result would be a first step in detecting individual optimization. More specific models are then needed to show that the original clutch size has highest fitness.

Statistics
In cases where the error structure was binomial, logistic regression was used, and Williams correction was used to correct for over-dispersion (Crawley, 1993Go). For Poisson regression we adjusted the scale parameter according to Crawley (1993Go, p. 262). Tests were based on likelihood ratios. Interactions were tested but not mentioned in the text when statistically non-significant.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Offspring fitness
Survival until next breeding season
The fraction of nests that fledged at least one chick did not differ between years or between manipulation categories (logistic regression controlled for year, manipulation {chi}2 = 4.977, df = 2, p >.08, full model deviance 46.699, df 163; Table 2). The control nests tended to be more successful. For the nests that had at least one fledgling (n = 162), the probability of fledging was not significantly dependent on manipulation (logistic regression controlled for year, after Williams correction: {chi}2 = 3.395, df = 2, p >.15, full model deviance 118.21, df 157). The probability of recruiting in the breeding population from fledging onwards was not affected by manipulation (logistic regression controlled for year, after Williams correction: {chi}2 = 2.545, df = 2, p >.3), nor did it differ between years (logistic regression controlled for manipulation, after Williams correction: {chi}2 = 2.390, df = 2, p >.3, full model deviance for both models 149.44, df 157).


View this table:
[in this window]
[in a new window]
 
Table 2 Effects of manipulation on offspring fitness components.

 
Poisson regression over the whole material (n = 168) showed that the number of recruits per first brood was affected by manipulation ({chi}2 = 13.44, df = 2, p =.001, scale parameter 1.491, full model deviance 152.15, df 165; Figure 1 and Table 3). The reduced group differed from the control group (p <.01) but the enlarged group did not differ significantly from the control group (p >.2) when tested pairwise. In a similar analysis using manipulation as a continuous variable ({chi}2 = 12.46, df = 1, p <.001, scale parameter 1.469, full model deviance 155.7, df 166), no effect of squared manipulation was detectable when added to the significant effect of manipulation ({chi}2 = 1.195, df = 1, p >.2), nor in interaction with year ({chi}2 = 1.119, df = 2, p >.5). There are thus no significant indications of diminishing returns with manipulated brood size.



View larger version (34K):
[in this window]
[in a new window]
 
Figure 1 The effect of brood size manipulations on different fitness components. Effects are shown per year separately (A-D: 1995 circles, 1997 triangles, and 1998 diamonds) and combined over the years (E-H). Given are: recruits per first brood (A and E), recruits per second brood (B and F), female survival (C and G), and male survival (D and H). Bars indicate SE. Manipulation effect could only be shown in the recruits per first brood

 

View this table:
[in this window]
[in a new window]
 
Table 3 Fitness components for the three manipulation categories.

 
Breeding in the next season
For 60 recruits, data on clutch size, laying date, and offspring born in their first breeding season were available (Table 4). Clutch size and the number of offspring born did not differ with manipulation history (ANOVA, controlled for year, F2,55 = 0.344, p >.7 and F2,55 = 0.474, p >.6, respectively), while differences between years were around significance (ANOVA, controlled for manipulation, F2,55 = 3.153, p =.051 and F2,55 = 3.371, p =.042, respectively). Laying date of the recruits, however, did differ both with manipulation history and year (ANOVA, manipulation effect controlled for year, F2,55 = 3.686, p =.031; ANOVA, year effect controlled for manipulation, F2,55 = 22.005, p <.001), the control group laying the latest and significantly later than the enlarged group (Tukey p <.05). We have no explanation for this effect. There were no exceptional differences in body mass, wing length, or tarsus in the nest between the survivors of the control group and the other manipulation categories.


View this table:
[in this window]
[in a new window]
 
Table 4 Effects of treatment on female offspring breeding performance in their first clutch in the next year.

 
Parental fitness
Second clutches
The probability of starting a second clutch did differ between years (logistic regression, {chi}2 = 30.0, df = 2, p <.001) but was not different between the manipulation categories ({chi}2 = 0.27, df = 2, p >.8, residual deviance 177.19, df 157). However, for those pairs that did produce a second clutch, the interval between the first and second clutch depended on manipulation (ANOVA, year effect, F2,49 = 5.814, p <.01, manipulation effect, F2,49 = 5.212, p <.01). Parents with enlarged broods took about one week longer to initiate a second clutch than parents with a reduced brood (Tukey p < 0.05). Recruitment of late clutches (second and repeat clutches) did not differ between manipulation categories (Figure 1b).

Parental survival and next year reproduction
The most adequate model that estimated survival and recapture rate of the adults simultaneously was the model with constant recapture and constant survival rate (Table 5). All more complicated models that were tested did not differ significantly from this model (likelihood ratio test all p >.25). We conclude that recapture rate was 0.897 (SE = 0.055) and the annual survival rate was 0.413 (SE = 0.0376), both not different between the sexes or with manipulation.


View this table:
[in this window]
[in a new window]
 
Table 5 Results of a MARK analysis of survival rate ({phi}) and recapture rate (p) as a constant (.) or a function of sex (s), manipulation (m), and time (t).

 
Also, the local survival of both males and females (when estimated simply as the fraction caught next year as a breeding bird) did not differ between years and were not related to manipulation (logistic regression males: n = 158, all p >.7; females: n = 168, all p >.09; Figure 1c,d). The large between-year variation in survival reduced the power of this analysis strongly. Clutch size, laying date, and the number of young born in the next year were not affected by the manipulation when controlled for year differences and analyzed for the females (n = 69) or for the males (n = 48).

Overall fitness
Fitness estimated as the sum of the parents and offspring surviving locally and identified in the next year was affected by manipulation (controlled for year, {chi}2 = 5.999, df = 2, p =.05, scale parameter 1.23, residual deviance 177.96, df 153). Mean values confirmed the positive selection pressure in this population (reduced = 1.245, control = 1.660, and enlarged = 1.891 locally survived birds per brood; Table 3).

Calculated for females separately, the same trend existed, but the manipulation effect was far from significant (Poisson regression: controlled for year, {chi}2 = 1.617, df = 2, p >.4, model deviance 122.30, df 167; mean values: reduced = 0.714, control = 0.893, and enlarged = 0.911 locally survived females per brood). This lack of significance relates to the high variation in adult female survival, especially in the control category (see Figure 1c). Recruitment of female offspring was strongly related to manipulation (Poisson regression: {chi}2 = 12.67, df = 2, p <.01, residual deviance 179.4, df = 165). In males the pattern was consistent and significant (controlled for year, {chi}2 = 7.574, df = 2, p <.05, scale parameter 1, residual deviance 177.1 df 153; mean values: reduced = 0.528, control = 0.800, and enlarged = 0.982 locally survived males per brood). This pattern is also generated by the significant manipulation effect on the number of male recruits (Poisson regression controlled for year, {chi}2 = 11.86, df = 2, p <.05, scale parameter 1, residual deviance 156.87, df 153). These results suggest different strength of selection on males and females, but this difference was not statistically significant (controlled for year and manipulation; sex: {chi}2 =0.376, df = 1, p >.5, interaction sex and manipulation: {chi}2 = 1.44, df = 2, p >.4, full model deviance 353.16, df 320).

We calculated {lambda} for each manipulation category on the basis of the data in Table 3. We assumed age-independent survival and reproductive rates and did not correct for dispersal. We used no maximum age but corrected for recapture probability. The fitness estimated in this way for the manipulation categories separately were reduced: {lambda} = 0.686, control: {lambda} = 0.935, and enlarged {lambda} = 1.048, indeed suggesting a positive selection pressure on clutch size. The reproductive values of the ‘mutants,’ reduced and enlarged in the resident population of control birds V(x, x*), were respectively 0.500 and 1.209, again showing directional selection pressure towards larger broods. This directional selection did not result in an increase in clutch size over the years (Figure 2).



View larger version (16K):
[in this window]
[in a new window]
 
Figure 2 Mean clutch size in the Lauwersmeer did not increase over the years. Bars are SD

 
Individual optimization of clutch size?
In this study there were no indications that fitness was maximized at the original clutch size. The total fitness was positively related to Y (the number of young in the nest) only. Y2 explained no additional variance (Table 6). The fact that M (manipulation), either alone or in interaction with Y or Y2, did not explain additional variance indicates that fitness depends on the number of young actually in the nest and not on traits of the birds or their resources (Table 6). This conclusion is consistent with the fact that original clutch size or brood size did not explain additional variance. Similar results were obtained in an analysis explaining variation in the number of recruits (Table 6). The fact that we cannot detect diminishing returns with manipulated brood size and that original clutch or brood size did not explain variation in fitness are strong arguments against the occurrence of individual optimization in this population.


View this table:
[in this window]
[in a new window]
 
Table 6 Results of Poisson regressions that do not support the individual optimization hypothesis.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Interpretation of brood size manipulations
In the Lauwersmeer population we have found strong evidence for positive selection on brood size but no response on this selection in terms of a changing clutch size over the years. There are a number of methodological and biological explanations possible for this lack of response that we will discuss here: (1) we estimated temporal variation in selection on clutch size rather than spatial variation, (2) we use incomplete fitness estimates because of dispersal, (3) there are potential clutch-size related costs in other phases of reproduction, and (4) the observed clutch size variation does not represent heritable variation.

(1) With regard to the temporal variation in selection on clutch size, the question is whether the detected differences in fitness are properties of the population under study or merely fluctuations over time, because most of the studies in the different populations have not been performed during the same years (see also Merilä et al., 2001Go). In our studies the Hoge Veluwe great tit population (1983–1987) and the Lauwersmeer population (1995, 1997, 1998) differed significantly in terms of the effect of manipulation on fitness when tested in one model (interaction population x manipulation, {chi}2 = 8.433, df = 2, p <.05, residual deviance 682.03, df 510). Only for the year 1998 can a direct comparison be made for the Hoge Veluwe (Visser and Lessells, 2001Go; recapture probabilities from Bauchau and van Noordwijk, 1995Go) and the Lauwersmeer (this study). In this year both populations showed a non-significant but positive selection pressure on clutch size (Hoge Veluwe: {lambda} control = 0.83 and {lambda} enlarged 1.12; Lauwersmeer: {lambda} control = 1.20 and {lambda} enlarged = 1.29). For the Hoge Veluwe, this is in contrast to the earlier results (1983–1987), in which selection was stabilizing (Tinbergen and Daan, 1990Go). Differences in selection between years are probably more a rule than an exception. Because we saw selection for larger brood size in all three years that we measured it in the Lauwersmeer, we believe that this is a population trait rather than a consequence of variation in time. Strictly speaking, longer synchronous data sets from brood size experiments in different areas are needed to come to a firmer conclusion as to what extent the differences found have a spatial or a temporal origin.

(2) Fitness estimates may be incomplete due to brood size dependent dispersal. Verhulst et al. (1997)Go showed that the dispersal pattern of young great tits can covary with their nestling mass, and thus it potentially biases fitness estimates that are based on local survival. The Lauwersmeer population differs from most other populations because it consists of a number of wood plots with a maximum distance of 10 km between the plots, surrounded by unsuitable breeding habitat for great tits. In this system, dispersal distance increased with experimental manipulation of brood size (Tinbergen, in preparation). Such dispersal may reduce the chance to find surviving offspring locally with increasing brood size. This could lead to underestimation of the positive selection pressure in the Lauwersmeer. If the same effect exists in other populations, it may seriously affect our judgment of local adaptation of brood size, and further analysis for other populations is also required.

(3) Selection may not act in the nestling phase alone; it has important components during other phases, such as egg laying or incubation. Recently, evidence has accumulated that the egg laying and incubation periods cannot be neglected when estimating selection pressures on clutch size (Monaghan and Nager, 1997Go). The question arises: what we can conclude about selection pressures on clutch size on the basis of brood size manipulation studies? Visser and Lessels (2001)Go estimated (in the Hoge Veluwe population) the fitness loss due to laying and incubation of extra eggs (comparison of their ‘full costs' treatment with their ‘free chicks’ treatment) to be around 0.3 female units for one to two eggs. If the cost of egg laying is indeed so high, it could also seriously change our view on selection acting on clutch size in tits.

(4) The observed clutch size variation did not represent heritable variation. When variation in clutch size is associated with variation in fitness, and this variation has a heritable component, we expect natural selection to occur (Endler, 1986Go). Using our experimental approach we showed that brood size variation was associated with fitness variation, and the general idea is that clutch size in great tits has heritable variation. When calculated for the whole Lauwersmeer population, the heritability of clutch size was indeed positive and significant (daughter-mother regression: h2 = 0.272, SE = 0.118, p =.02, n = 242). However, when controlling for plot for both daughter and mother, this value remained positive but became insignificant (h2 = 0.166, SE = 0.122, p >.15), suggesting that a considerable part of the variation was induced by the environment. The lack of response to the selection for larger clutches may thus be caused by a lack of (or a weak) heritable variation in clutch size in this population. Such an explanation for the lack of response to selection was formulated earlier by Price et al. (1988)Go and was recently illustrated by the work of Kruuk et al. (2002)Go on antler size in red deer. Variation in antler size was associated with lifetime breeding success and had a heritable component, but fitness differences seem to be largely associated with the environmentally caused variation in antler size.

Comparison to other studies
Only when we assume that the above-mentioned biases in the fitness measurements can be neglected is it useful to discuss the results of the current study in relation to others. We compiled the outcome of brood size manipulation studies in great tits over the last few decades to see what the conclusions of these experiments were regarding the selection on brood size (Table 7) within this well studied species. Fitness effects as judged by the authors are given for parents and brood separately. Their judgment of the overall effect is also given. For the great tit, there are nine studies based on data from five woods.


View this table:
[in this window]
[in a new window]
 
Table 7 Overview of the fitness effects of brood size manipulations in the great tit.

 
The great tit data of Wytham wood (UK) have been analyzed by different people using basically the same data set. In their brood size manipulation study, Perrins and Moss (1975)Go signaled that the most productive brood size estimated from experimental variation was close to the original brood size, in contrast to the larger most productive brood size calculated on the basis of natural variation in brood size. They suggested that this contrast might be a result of the fact that tits adjust their clutch size to produce their own most productive brood size in terms of fitness (later called the individual optimization hypothesis, Pettifor et al., 1988Go). Boyce and Perrins (1987)Go analyzed the effect of a fluctuating environment on the geometric mean fitness of the most productive brood size in the same data set extended with new years. This analysis led to the different view that also in the case of artificial variation in brood size, the most productive brood size was larger than the original brood size. Their point was that bad seasons affected larger broods relatively more, reducing optimal brood size in terms of geometric mean fitness below the average most productive brood size. No individual adjustment of clutch size was needed to explain the observed patterns. Pettifor et al. (1988Go, 2001Go) did a reanalysis of the same data set but used different statistical models. Their analysis, more consistent with the view in the paper of Perrins and Moss (1975)Go, led to the conclusion that stabilizing selection on individual brood size (individual optimization of clutch size) is the rule in the Wytham great tit population. Apparently, for Wytham the conclusion depends on the analysis and/or the fitness measure used. Because analyses are complicated and different fitness measures have been used, comparison between studies es hampered to some extent. Pettifor et al. (2001)Go clearly favor the individual optimization hypothesis as an explanation for the observed patterns, but they did not discuss why Boyce and Perrins (1987)Go reached a different conclusion using the same data.

Experimental studies on Gotland (Lindén, 1990Go) and on the Hoge Veluwe (Tinbergen and Daan, 1990Go) confirmed Pettifor's view. These authors also found evidence for stabilizing selection on brood size in line with the individual optimization hypothesis. However, Verhulst (1995)Go found a population with strong selection for smaller broods. He reduced brood size in an island population of great tits and showed that this reduction led to a fitness increase. Together with the results of the current study, we think that the claim of Pettifor et al. (2001)Go, that brood size experiments are generally in favor of the individual optimization hypothesis, does not hold.

The fact that there are populations with selection for larger (this study) and for smaller (Verhulst, 1995Go) broods is consistent with the gene flow hypothesis (Dhondt et al., 1990Go; Dias and Blondel, 1996Go). As a matter of fact, Perrins and Moss (1975)Go also suggested that clutch size could be a product of differing selection in different habitats. However, other explanations can be given for a lack of microevolution under directional selection in natural populations, as recently reviewed by Merilä et al. (2001)Go.

The outcome of our work urges us to formulate how ideas of individual optimization blend with ideas on microevolution. Suppose that in a population with perfect individual optimization of clutch size, a heritable component of clutch size exists. If larger clutches have higher fitness (Perrins and Moss, 1975Go; Tinbergen and Daan, 1990Go), we would expect directional selection towards larger clutches and need an explanation why this generally does not happen.

If individual optimization of clutch size were perfect, we would not expect a heritable component in the variation in clutch size. The flexible phenotype of the omniscient bird would produce exactly the number of eggs at a time that maximizes fitness. In this view, variation in clutch size would be environmentally determined and we expect competition for the better places, in concert with density dependent effects, to govern clutch size variation (Both et al., 2000Go).

In this paper we showed that flexible phenotypes are not always perfect. Great tits laid clutches that were not locally adapted. This is in line with the fact that the environment is not perfectly predictable. Even if individuals can adjust their reproduction on the basis of a learning process, adjustment will lag behind. On a longer time scale, adaptation may take place via microevolution. In this process, natural selection will shape the genetic architecture, and whether or not we will find locally adapted traits will depend on the interactions between selection, the genetic architecture, and spatial and temporal variation in the environment.


    ACKNOWLEDGEMENTS
 
We thank J. Willems, S. Krap, and H. Hut from ‘Staatsbosbeheer’ and J. H. Janssen, E. K. Hansen, and G. Altena from the ‘Koninklijke Landmacht’ for their permission to work in the Lauwersmeer area. We are grateful for the help of many students. We profited a lot from critical reading of the paper by C. Both, R. H. Drent, T. Piersma and J. Komdeur. J.J.S. was supported by a postdoctoral grant from the Ministerio Español de Educación y Cultura and by a contract from projects PB94-0070-C02-01 and PB97-1233-C02-01.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Bauchau V, Van Noordwijk AJ, 1995. Comparisons of survival estimates obtained from three different methods of recapture in the same population of the great tit. J Appl Stat 22:1031-1037.[CrossRef]

Blondel J, Maistre M, Perret P, Hurtrez-Bouses S, Lambrechts MM, 1998. Is the small clutch size of a Corsican blue tit population optimal? Oecologia 117:80-89.[CrossRef]

Both C, Tinbergen JM, Visser ME, 2000. Adaptive density dependence in clutch size. Ecology 81:3391-3403.[CrossRef]

Boyce MS, Perrins CM, 1987. Optimizing great tit clutch-size in a fluctuating environment. Ecology 68:142-153.[CrossRef]

Charnov EL, Krebs JR, 1974. On clutch size and fitness. Ibis 116:217-219.

Crawley MJ, 1993. GLIM for ecologists. Oxford: Blackwell Scientific Publications.

Daan S, Dijkstra C, Tinbergen JM, 1990. Family planning in the kestrel (Falco tinnunculus): the ultimate control of variation in laying date and clutch size. Behaviour 114:83-116.

Dhondt AA, Adriaensen F, Matthysen E, Kempenaers B, 1990. Nonadaptive clutch sizes in tits. Nature 348:723-725.[CrossRef]

Dias PC, Blondel J, 1996. Local specialization and maladaptation in the Mediterranean blue tit (Parus caeruleus). Oecologia 107:79-86.

Endler JA, 1986. Natural selection in the wild. Princeton: Princeton University Press.

Fisher RA, 1930. The genetical theory of natural selection. Oxford: Clarendon Press.

Gustafsson L, Sutherland WJ, 1988. Cost of reproduction in the collared flycatcher Ficedula albicollis. Nature 335:813-817.[CrossRef]

Kawecki TJ, 1995. Demography of source-sink populations and the evolution of ecological niches. Evolutionary Ecology 9:38-44.

Kruuk LEB, Slate J, Pemberton JM, Brotherstone S, Guinnes F, Clutton-Brock TH, 2002. Antler size in red deer: heritability and selection but no evolution. Evolution 56:1683-1695.[CrossRef][ISI][Medline]

Lessells CM, 1993. The cost of reproduction: do experimental manipulations measure the edge of the options set? Etologia 3:95-111.

Lindén M, 1990. Reproductive investment and its fitness consequences in the great tit Parus major. PhD thesis. University of Uppsala, Uppsala, Sweden.

McNamara JM, 1997. Optimal life histories for structured populations in fluctuating environments. Theoretical Population Biology 51:94-108.[CrossRef]

Merilä J, Sheldon BC, Kruuk LEB, 2001. Explaining stasis: microevolutionary studies in natural populations. Genetica 112–113:199-222.[CrossRef]

Monaghan P, Nager RG, 1997. Why don't birds lay more eggs. Trends Ecol Evol 12:270-274.[CrossRef]

Nur N, 1984a. The consequences of brood-size for breeding blue tits. 1. Adult survival, weight change and the cost of reproduction. J Anim Ecol 53:479-496.[CrossRef]

Nur N, 1984b. The consequences of brood-size for breeding blue tits. 2. Nestling weight, offspring survival and optimal brood size. J Anim Ecol 53:497-517.[CrossRef]

Nur N, 1986. Is clutch size variation in the blue tit (Parus caeruleus) adaptive? An experimental study. J Anim Ecol 55:983-999.[CrossRef]

Perrins CM, Moss D, 1975. Reproductive rates in the great tit. J Anim Ecol 44:695-706.

Pettifor RA, 1993a. Brood-manipulation experiments. I. The number of offspring surviving per nest in blue tits (Parus caeruleus). J Anim Ecol 62:131-144.[CrossRef]

Pettifor RA, 1993b. Brood-manipulation experiments. II. A cost of reproduction in blue tits (Parus caeruleus)? J Anim Ecol 62:145-159.[CrossRef]

Pettifor RA, Perrins CM, McCleery RH, 1988. Individual optimization of clutch size in great tits. Nature 336:160-162.[CrossRef]

Pettifor RA, Perrins CM, McCleery R, 2001. Individual optimization of fitness: variation in reproductive output, including clutch size, mean nestling mass and offspring recruitment, in manipulated broods of great tits Parus major. Journal of Animal Ecology 70:62-79.[CrossRef]

Price T, Kirkpatrick M, Arnold SJ, 1988. Directional selection and the evolution of breeding date in birds. Science 240:798-799.[Abstract/Free Full Text]

Sanz JJ, Tinbergen JM, 1999. Energy expenditure, nestling age, and brood size: an experimental study of parental behavior in the great tit Parus major. Behav Ecol 10:598-606.[Abstract/Free Full Text]

Tinbergen JM, Both C, 1999. Is clutch size individually optimized? Behav Ecol 10:504-509.[Abstract/Free Full Text]

Tinbergen JM, Daan S, 1990. Family planning in the great tit (Parus major): optimal clutch size as integration of parent and offspring fitness. Behaviour 114:161-190.

van Balen JH, 1973. A comparative study of the breeding ecology of the great tit Parus major in different habitats. Ardea 61:1-93.

Verhulst S, 1995. Reproductive decisions in great tits: an optimality approach. PhD thesis. University of Groningen, Groningen, The Netherlands.

Verhulst S, Perrins CM, Riddington R, 1997. Natal dispersal of great tits in a patchy environment. Ecology 78:864-872.[CrossRef]

Visser ME, Lessells CM, 2001. The cost of egg production and incubation in great tit. Proc R Soc Biol Sci B 268:1271-1277.[CrossRef]

White GC, Burnham KP, 1999. Program MARK. Bird Study 46:120-138.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Behav EcolHome page
A. P. Moller
Interval between clutches, fitness, and climate change
Behav. Ecol., January 1, 2007; 18(1): 62 - 70.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Lay Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (7)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Tinbergen, J. M.
Right arrow Articles by Sanz, J. J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Tinbergen, J. M.
Right arrow Articles by Sanz, J. J.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?