Behavioral Ecology Advance Access originally published online on January 8, 2008
Behavioral Ecology 2008 19(2):344-352; doi:10.1093/beheco/arm136
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Evolution of parental favoritism among different-aged offspring
Department of Psychology, University of Texas at Austin, Austin, TX 78712, USA
Address correspondence to J. Jeon, who is now at the Division of EcoScience, Ewha Womans University, Seoul 120-750, Republic of Korea. E-mail: evopsy{at}gmail.com.
Received 16 May 2007; revised 16 November 2007; accepted 16 November 2007.
| ABSTRACT |
|---|
|
|
|---|
The theories of intrafamilial conflict and parental investment have yet to examine how parents' decisions about resource allocation are influenced by the fact that their offspring may be different in age. Two counteracting effects of offspring growth on parental allocation of resources have deterred the development of a formal model: A parent may favor older offspring due to their greater reproductive value or favor younger offspring due to their higher marginal returns from extra resources. Using an evolutionary invasion analysis in class-structured populations, I present a formal model that explores how a parent should allocate its resources among different-aged offspring from the viewpoint of the parent. The parent's evolutionarily stable strategy is to allocate its resources such that the marginal benefit to each offspring's survival, weighted by the survival probability to the reproductive age, is equal to the marginal cost to the parent's residual survival. Two general situations are considered in which younger offspring obtain higher marginal returns than older offspring. In nearly all circumstances, a parent is expected to bias its resources toward older offspring. The results may account for the widespread yet puzzling phenomenon of parental bias toward older offspring in view of previous theories of intrafamilial conflict.
Key words: evolutionarily stable strategy, hatching asynchrony, invasion analysis, offspring age, parental investment, reproductive value.
| INTRODUCTION |
|---|
|
|
|---|
Parents are expected to adaptively allocate their limited resources among offspring so as to maximize parental fitness (Smith and Fretwell 1974
Within-brood resource allocation, however, should be shaped not only by the offspring's cryptic internal condition (commonly called "need") but also by the offspring's external condition: information that parents may directly access, such as offspring age, sex, and competitive ability (Glassey and Forbes 2002
). In particular, age differences among offspring are common in asynchronously hatching birds and many mammals. Hatching asynchrony in altricial birds greatly influences food allocation among nestlings with earlier hatched (older) nestlings getting more food, leading to higher mortality of later-hatched (younger) nestlings (Magrath 1990
). These findings have been interpreted as the consequence of scramble competition among nestlings, where parents exert little or no control over resource distribution (Macnair and Parker 1979
; Parker and Macnair 1979
). Yet, parents may have full or partial control over food distribution and preferentially invest in older offspring. Even when parents seem to passively feed more competitive offspring, it may be in parental interests to do so. Therefore, to understand how both cryptic (internal) and noncryptic (external) conditions of offspring affect family dynamics, it is necessary to elucidate the optimal distribution of resources among different-aged offspring from the viewpoint of parents.
Unfortunately, no formal investigations have been made about how parents should allocate resources among offspring of different ages (Clutton-Brock 1991
). Making realistic predictions has been considered very difficult due to 2 potentially counteracting effects of offspring growth on the ES level of parental provisioning. Parents may favor older offspring because older offspring have higher reproductive value (i.e., should make a higher genetic contribution to future populations) (Fisher 1930
). This is both because older offspring are closer to reproductive maturity and because the instantaneous rate of juvenile mortality tends to decrease with increasing age (Clutton-Brock 1991
). On the other hand, parents may favor younger offspring because younger offspring have greater need for food than older offspring, so the effect of a unit of additional investment on offspring survival will be higher for younger offspring (Rubenstein 1982
; Sargent and Gross 1986
; Clutton-Brock 1991
).
The results of recent empirical studies of food allocation in asynchronously hatched broods pose a serious problem for current models of parent–offspring interactions (Lotem 1998
; Cotton et al. 1999
; Smiseth and Amundsen 2002
). According to the honest signaling model's definition of need, younger (smaller) nestlings are in greater need of food; hence, they should beg more vigorously and receive more food than older (larger) nestlings. However, many experimental studies have found that older nestlings actually receive more food, even though they beg less intensively than their younger siblings (Kilner 1995
; Price et al. 1996
; Lotem 1998
; Cotton et al. 1999
; Smiseth and Amundsen 2002
). This striking discrepancy between theoretical predictions and empirical findings may indicate that fundamental aspects of parent–offspring interactions have yet to be understood (Cotton et al. 1999
). Meanwhile, the scramble competition model assumes that offspring fully control food allocation within broods and predicts that older offspring with higher competitive ability should get more food than younger offspring, which is consistent with most observed allocation patterns in asynchronously hatched broods (Mock and Parker 1997
; Parker et al. 2002
). However, there remains the possibility that parents may have partial or full control over food distribution and actively prefer feeding older offspring. Indeed, recent empirical evidence suggests that parents may play an active role in food allocation within asynchronous broods (Kilner 1995
; Krebs et al. 1999
; Smiseth et al. 2003
).
To accurately assess the consequences of parental allocation of resources among offspring of different ages, it is necessary to construct evolutionary models of class-structured populations in which individuals fall into different age classes. In such models, Fisher's (1930)
concept of reproductive value acts as a weighting factor for comparing the effects of changes in age-specific survival and fecundity on a common scale (Taylor 1990
; Charlesworth 1994
; Taylor and Frank 1996
; Frank 1998
; Pen and Weissing 2002
). After envisioning a genetically homogeneous population with a fixed allocation strategy among offspring of different age classes, one can determine whether the resident population can be invaded by a rare mutant allocation strategy. This approach can reveal the outcome of long-term evolution and has been recently referred to as an evolutionary invasion analysis (for a clear introduction, see Otto and Day 2007
). Assuming parental control over resource allocation, I present a model that investigates how a parent should invest in offspring of different age classes. Within the model's life-history framework based on reproductive value, the effect of additional food on the survival of younger offspring will be allowed to be higher than the effect on survival of older offspring. The model presented here thus integrates the 2 counteracting effects of offspring growth on parental allocation strategies within broods.
| THE MODEL |
|---|
|
|
|---|
Dynamics of resident population
Consider a diploid resident population that is structured into parents (class P) and their offspring, which are classified into n age classes (class Oi : i = 1, ..., n). For the sake of simplicity, the model employs hermaphroditic individuals that reproduce as both males and females. During a time step, an individual parent mates with another adult, producing f newborn offspring, allocates its limited resources among its different-aged offspring, and survives to the next time step with probability
. The population is further structured into distinct broods, for which each parent provisions only its own offspring within the brood. The length of the time step corresponds to the average age gap between offspring age classes; it will be measured in days for asynchronously hatching birds, whereas it will typically be measured at a larger time scale for mammals (the model thus focuses on within-year dynamics for the case of asynchronously hatching birds). Given that offspring provisioning usually occurs over a period of time, it would be desirable to consider dynamic aspects of parent–offspring interactions as well. Unfortunately, such an analysis is extremely difficult to carry out. Following other authors (e.g., Parker et al. 1989
When a focal parent invests in its current brood composed of different-aged offspring, it suffers a cost in terms of reduced future survival. The average number of age-i offspring per adult in the population is denoted as ui. The per capita amount of resources that an age-i offspring receives is xi (it is assumed that parents can discriminate between offspring of differing ages with negligible cost). Thus, the total amount of resources that an average parent provides for its current brood per time step is proportional to Ex = 
uixi. An average parent's allocation strategy is denoted by a vector x = [x1, x2, ..., xn], where each trait value xi can be estimated by its genotypic value plus the unexplained residual. It is assumed that the n trait values are genetically uncorrelated with each other. An age-i offspring survives and grows into an offspring aged i + 1 with survival probability
i =
i(xi), an increasing function of xi. A focal parent's survival probability into the next time step is
(Ex), a decreasing function of Ex. Figure 1 depicts the life cycle for the resident population. The population transition matrix A corresponding to Figure 1 is
![]() | (1) |
|
The resident population characterized by the matrix A will eventually reach a constant rate of geometric growth: A*u* =
u*, where
is the dominant eigenvalue of the resident matrix A* at equilibrium and the column vector of stable age distribution u* = [u
,u
, ... ,u
,u
] is the dominant right eigenvector of A*. Because u* is determined up to a constant, I let the number of parental class u
equals 1.0 so that u
could represent the average number of age-i offspring per adult in the population at equilibrium. The row vector of individual reproductive values v* = [v
, v
, ..., v
, v
] is the dominant left eigenvector of A* : v*A* =
v*. The 2 vectors u* and v* are obtained in Appendix A.
It is unrealistic, however, to expect that the resident population will grow indefinitely: Density-dependent factors will ensure that the population size is no longer changing. The constraint
= 1 requires that the following equation be satisfied
|
| (2) |
= 
u
x
. Equation 2 means that, by the operation of density dependence, the rate at which new reproducing adults are introduced into the population (the left-hand side) is balanced by the death rate of existing adults (the right-hand side). The assumption of density dependence implies that it should be decided which vital rate in the life cycle will be influenced by the population density (Mylius and Diekmann 1995
1(x1)) can be expressed in terms of other vital rates such that Equation 2 is satisfied (in this system, it does not make any difference if I choose other vital rates to be density dependent).
One may suspect that the resident matrix A in Equation 1 does not accurately describe the actual population that it is meant to. In order for a newborn age-1 offspring to become an age-(n + 1) adult, it should be continually cared for by its living parent throughout n time steps (as implied in
1(x1), ...,
n(xn)), yet the matrix A also assumes that there will be a nonzero parental mortality 1 –
(Ex) for each time step. In other words, although the model assumes that each and every age-i offspring will receive xi from the parents and survive to the next time step with probability
i(xi), it seems that some offspring will never receive xi because they would lose their parents across time. Appendix B clarifies this issue and shows that each and every age-i offspring can be regarded as receiving xi from its living parent at any time step.
Invasion analysis
To find the ES allocation strategy among different-aged offspring, I seek to determine under what conditions a rare mutant allocation strategy x°( = [x
, ..., x
], where x
= x
+
x) can invade the resident population fixed for the resident strategy x*. The mutant strategy x° is assumed to be rare enough that the average number of age-i offspring per mutant adult is given by u
, that is, the offspring's age distribution in the resident population with x* (with the rarity assumption, the eigenvectors of the resident population, i.e., u* and v*, can be used to characterize the mutant dynamics. See Taylor and Frank [1996
, p. 36] for the proof). The transition matrix of the mutant subpopulation is
![]() | (3) |
The invasion fitness of a mutant strategy x° is sought as the initial growth rate of the mutant subpopulation within the resident population, that is, the dominant eigenvalue
(x°, x*) of A(x°, x*) (Metz et al. 1992
). Therefore, for a resident strategy x° to be ES,
(x°, x*) should have a local maximum at x° = x* such that
|
| (4) |
(x°, x*) is negative semidefinite (i.e., all its eigenvalues are less than or equal to zero) at the equilibrium point (x
, ..., x
) (Day and Taylor 2003|
| (5) |
Note that although the correct fitness measure of a rare mutant strategy x° is
(x°, x*), one can instead use W(x°, x*) = v*A(x°, x*)u* as the fitness function because 
/
x° and
W/
x° have the same sign (Taylor and Frank 1996
). From Equation 5 and Appendix A, a parent's ES allocation strategy x*( = [x
, ..., x
])is obtained by n equilibrium conditions as follows:
![]() | (6) |
= 
u
x
. The n equilibrium conditions emphasize that the n trait values comprising the vector of allocation strategy are coevolving. Equation 6 indicates that, at the ES allocation from a parent's viewpoint, the marginal benefit from additional resources for each offspring's survival
'i(x
), weighted by the survival probability of each offspring to reproductive age (i.e., 



), is equal to the marginal cost of current total expenditure to the parent's residual survival
' (E
). The second-order condition for evolutionary stability, as well as the condition for convergence stability, is verified in Appendix C. | RESULTS |
|---|
|
|
|---|
To derive concrete predictions about the parental ES allocation of resources among different-aged offspring, I allow the effect of additional resources on the survival of younger offspring to be greater than the effect on the survival of older offspring, in accordance with previous research (West-Eberhard 1975
The ES pattern of parental resource allocation can be obtained as follows. The survival probability of an offspring is assumed to increase with the resources provided, with diminishing returns (i.e.,
'(x) > 0,
''(x) < 0) (Smith and Fretwell 1974
). Consider 2 offspring spaced by 1 time step. At evolutionary equilibrium, each offspring will receive its ES level of resources (i.e., x
and x
, where i = 1,..., n – 1) from the parent. Because
'(x) decreases with x, whether x
is larger or smaller than x
can be determined by comparing the younger offspring's marginal return at x
with its marginal return at x
(i.e., comparing
'i(x
) with
'i(x
), respectively). Let us denote
'i(x
) –
'i(x
) by
(x
). Note that
(x
) can be regarded as a function of x
only; the equilibrium conditions in Equation 6 imply that I can rewrite
'i(x
) –
'i(x
) as
|
| (7) |
The sign of
(x
) informs whether x
is larger or smaller than x
. For instance, a plus sign indicates that x
is less than x
(i.e., the older offspring receives more resources than the younger offspring at ESS).
A younger offspring yields higher marginal returns for every level of resources provided
If a younger offspring could obtain a higher marginal return than an older offspring for every level of resources provided by the parent during a given feeding episode (i.e.,
'i(x) >
'i + 1(x) for all x), then it can be shown that
(x
) =
'i(x
) –
'i(x
) > 0 for all levels of x (see Appendix D for the proof). Because
'(x) decreases with x, it is revealed that x
< x
. Therefore, rather surprisingly, if a younger offspring yields higher marginal returns than an older offspring for every level of resources provided during a feeding bout, then the parent's ES allocation strategy x* (=[x
, ..., x
]) is to distribute its resources such that older offspring get more resources than younger offspring (i.e., x
< x
<
< x
).
I illustrate the above result by assuming explicit functions for an age-i offspring survival:
|
| (8) |
1) and c, a so-called "shape constant" (Parker et al. 1989
i(xi) rises to Ki. Offspring of distinct ages differ only in their asymptotes: Ki is assumed to be larger than Ki+1, which ensures
'i(x
) >
'i(x
) for every value of x (Figure 2a). In other words, the asymptote Ki of an age-i offspring survival curve will be inversely related to the offspring's age: K1 > K2 > ... > Kn. Throughout this paper, I assume that the parent's residual survival is an exponentially decreasing function of current parental expenditure (i.e., each increment of parental investment is more costly than the previous increment):
|
| (9) |
|
A younger offspring yields higher marginal returns only when the resources already provided are none or few
Suppose that a younger offspring yields a higher marginal return than an older offspring with none or few resources provided, but it yields a lower marginal return with high resources. In this situation, it can be proved that
(x
) =
'i(x
) –
'i(x
) is positive and hence x
< x
(see Appendix E for the proof). Therefore, as in the previous situation, the parent's ES allocation is to bias its resources in favor of older offspring (i.e., x
< x
<
<x
).
To illustrate the above result, let an age-i offspring's survival be given by
|
| (10) |
'i(x
) >
'i(x
) for small values of x and yet
'i(x
) >
'i(x
) for large values of x, ci is assumed to be larger than ci+1. In other words, the shape constant ci of an age-i offspring survival curve will be inversely related to the offspring's age: c1 > c2 > ... > cn (Figure 3a) (Sargent and Gross 1986
|
Figure 3b gives the ES allocation of resources between 2 different-aged offspring with the younger offspring's shape constant c1 being fixed at a certain value. As expected, despite the younger offspring's higher marginal return when no or few resources have been provided, the parent will preferentially invest in the older offspring throughout the range of the older offspring's shape constant c2 (i.e., whenever there is any age difference between the 2 offspring). The graph also reveals that as c2 decreases, whereas c1 is fixed (hence, as the older offspring's age increases relative to the younger offspring's fixed age), the older offspring's resource share increases and the younger offspring's resource share, to a lesser extent, decreases. In other words, the larger the age spacing between 2 siblings, the more favored the older offspring.
In terms of the explicit functions examined above, the survival curves of younger and older offspring may differ both in their asymptotic levels of survival (K) and in their shape constants (c). This case is, however, a variation of the second situation in which a younger offspring yields higher marginal returns only when the resources provided are none or few. Hence, the main conclusion of biased allocation toward older offspring remains the same (Appendix E). On the other hand, it can be argued that 2 different-aged offspring will differ in the minimum threshold levels of resources at which their survival probabilities start increasing, hereafter referred to as the "threshold level of resources." For instance, older offspring will tend to have larger body sizes, thus requiring a higher threshold level of resources than the younger offspring. I investigated the consequence of different threshold levels of resource requirements between siblings. The analysis revealed that parents will still bias their resource allocation in favor of the older offspring, as long as the older offspring's threshold level of resources is greater than or equal to the younger offspring's threshold. An "exceptional" allocation favoring the younger offspring may be favored if the younger offspring's threshold is sufficiently greater than the older offspring's one; yet it is difficult to imagine how a younger offspring's threshold level of resources could be higher than the older offspring's one in natural populations.
| DISCUSSION |
|---|
|
|
|---|
Using an evolutionary invasion analysis of class-structured populations, the current model examined how a parent should distribute its limited resources among offspring of different ages. At the parent's ESS, the marginal benefit to each offspring's survival, weighted by the survival probability to the reproductive age, is equal to the marginal cost to the parent's residual survival. After previous discussions, I considered 2 general situations in which younger offspring obtain a higher marginal return from additional resources than older offspring. The model predicted that a parent should bias its resources in favor of older offspring in nearly all circumstances and that the degree of parental favoritism will increase as the age spacing between 2 offspring increases.
Comparison with traditional allocation theory
Traditional theory for the allocation of limited resources among alternative options suggests that, when offspring vary in quality (e.g., offspring age, sex, or condition), a parent's optimal allocation is to adjust its investment such that the marginal returns in each offspring are all equal (Lloyd 1985
, 1988
; Temme 1986
; Haig 1990
; Lessells 2002
). Accordingly, when offspring are of different ages, a parent's optimal allocation of resources would be to allocate resources in a way that the marginal gains of each offspring's survival are all identical at evolutionary equilibrium:
'1(x
) =
' 2(x
) =
=
'n(x
) = –
'(x
+
+ x
), where
'i (x
) is the marginal gain in terms of offspring's survival from investment in an age-i offspring and
'(x
+
+ x
) is the marginal cost of current expenditure to the parent's residual survival.
By contrast, the class-structured model presented here shows that the parental ESS is to equalize the marginal gains in each offspring's survival weighted by the probability of survival to reproductive age (see Equation 6). In this scheme, the consequences of the parental actor's investment in recipient offspring structured into different age classes are measured on a common scale, that is, the extent to which current investment in an age-i offspring eventually contributes to that offspring's growth into a reproducing adult. Interestingly, due to the weighting by the survival probability to reproductive age, the marginal gain from investment in a younger offspring must be higher than that from investment in an older offspring at the parental ESS:
'1(x
) >
'2(x
) >
>
'n(x
) = –
'(E
) (see Equation 7, which shows that
'i(x
) =
'i + 1(x
) /
i + 1(x
) >
'i + 1(x
)). Contrary to traditional allocation theory, which insists on the equilibration of marginal gains across all offspring at ESS, the model presented here sheds new light on the evolution of parental investment allocation: the marginal gains from each offspring are not necessarily equal at ESS. Rather, when offspring are structured into different age classes, the marginal gains from investment in younger offspring will be greater than those from investment in older offspring at the parental ESS. What are equalized in this case are the marginal gains from each offspring weighted by the offspring's probability of survival to reproductive age. Although the model presented here does not take into account the role of offspring begging, the fact that younger offspring obtain a higher marginal return than older offspring (i.e., they are needier) at evolutionary equilibrium may partly explain why younger offspring tend to beg more vigorously than older offspring in real situations (Kilner 1995
; Price et al. 1996
; Lotem 1998
; Cotton et al. 1999
; Smiseth and Amundsen 2002
).
One may wonder to what extent the main result of the current model, that is, biased allocation of parental resources toward older offspring, depends on the assumption of higher marginal benefits for younger offspring. Remarkably, the result is fairly robust even under the opposite assumption. If investment in an older offspring yields higher marginal returns for any level of resources already provided, then the model predicts that the parent will favor the older offspring. If investment in the older offspring yields higher marginal returns only when no or few resources have already been provided, then a parent will also favor the older offspring on the condition that the age gap between 2 offspring is relatively small (Jeon J, unpublished data). These insights suggest that the results presented here may be widely applicable across a variety of family systems.
Implications for models of intrafamilial conflict
The 2 counteracting effects of offspring growth on parental resource allocation were thought to be 1) older offspring's higher reproductive value (hence, parents should prefer older offspring) and 2) younger offspring's higher fitness increment from extra resources (hence, parents should prefer younger offspring) (Clutton-Brock 1991
). After constructing a structured population model based on reproductive value, I allowed younger offspring to gain a higher marginal return from additional resources than older offspring under 2 different sets of conditions. Consequently, the model presented here disentangles and integrates the 2 counteracting effects of offspring growth on parental allocation of resources.
Under nearly all circumstances, the model generates a clear-cut prediction that a parent will preferentially invest in older offspring over younger offspring. An exceptional allocation favoring younger offspring is possible only when the younger offspring's threshold level of resources is sufficiently greater than that of the older offspring, which is likely to be rare in real populations. Nevertheless, there are a handful of studies reporting equal allocation of resources or biased allocation toward younger offspring in asynchronously hatched broods (e.g., Leonard and Horn 1996
; Krebs et al. 1999
; Stamps et al. 1985
). It would be interesting to investigate if the precondition of exceptional allocation is fulfilled in such species.
Provided that parents win the parent–offspring conflict and, hence, can impose the parental optimum, the main prediction of biased allocation toward older offspring may well account for a widespread yet enigmatic phenomenon observed in most asynchronously hatching birds: Younger nestlings indeed receive less food than older nestlings, even though the former beg more intensively (Kilner 1995
; Price and Ydenberg 1995
; Price et al. 1996
; Lotem 1998
; Cotton et al. 1999
; Smiseth and Amundsen 2002
). Existing models of parent–offspring interactions, lacking in formal analysis of class-structured populations, have failed to explain this phenomenon. In particular, the honest signaling model predicts that parents should bias their resources toward the needier (i.e., younger) offspring, which runs counter to abundant empirical data. The model presented here provides a straightforward solution for this hitherto unresolved problem.
In addition to the lack of formal analysis of class-structured populations, the honest signaling model suffers from a conceptual problem in defining offspring's cryptic condition (or need) in terms of the shape constant c of offspring's survival functions. In such models, needier offspring (hungrier offspring in the short term; younger [smaller] offspring in the long term) correspond to those yielding a higher marginal gain from additional resources (i.e., 
/
x) (Price et al. 1996
). Unfortunately, in presenting an explicit example for their analysis, the honest signaling theorists assumed that needier offspring correspond to those having a lower value of c in their survival functions compared with their siblings (Godfray 1991
, 1995
). In reality, as the Results section shows, the c values of survival curves do not reliably reflect marginal gains from extra resources (
/
x): A survival curve with a higher c value may yield either a higher or lower marginal gain than a curve with a lower c value, depending on the amount of resources, x, already provided (see Figure 3a) (Parker et al. 2002
; Royle et al. 2002
). Hence, in the honest signaling model's terminology, younger offspring should not be regarded as being always needier because they will obtain a higher marginal return only when the resources already provided are none or few. Future research should focus on clarifying ambiguous terms like offspring need and "condition" found in current literature.
Limitations and conclusions
Several limitations of the present study should be noted. First, the current model assumes that parents always win the parent–offspring conflict and, hence, that the amount of resources transferred in nature is at the parental optimum. Yet the optimal resource allocations from the perspective of each different-aged offspring will differ, and it is possible that offspring may (at least partly) win the parent–offspring conflict. Future research should address how the evolutionary conflict among various family members is resolved within an age-dependent life-history context. Second, the present study should be distinguished from Lack's (1947
, 1954
) interpretation of hatching asynchrony as a low-cost means of adjusting brood size by imposing within-brood competitive asymmetries. The present study simply asked how parents should allocate their resources within asynchronously hatched broods but not how adaptive brood reduction is made possible due to hatching asynchrony. Third, like other models of parent–offspring interactions (e.g., Godfray 1995
; Parker et al. 2002
), the present model implicitly assumes that environmental conditions are favorable throughout the investment period and, hence, that enough food is available to be divided among all offspring. Yet empirical evidence suggests that parental provisioning decisions may be contingent on food availability (e.g., Boland et al. 1997
; Kilner 2002
; Bize et al. 2006
). To address this issue, the current model should be expanded to include the effect of food availability. Fourth, the current model assumed hermaphroditic populations to avoid the complexity resulting from dioecy. Although the results obtained here should apply well to dioecious populations, it must be recognized that most birds have biparental care and the 2 parents can follow different allocation rules (Lessells 2002
).
Theories of intrafamilial conflict, by their very nature, deal with social interactions between individuals of different classes; yet researchers have been slow to explicitly model the evolution of class-structured populations. The class-structured model presented here asked how a parental actor's behavior is influenced by the fact that its recipient offspring are not identical but classified into distinct age classes. The model's central prediction appears to well account for a hitherto puzzling phenomenon observed in nature, which suggests that the study of intrafamilial interactions should benefit from the application of evolutionary invasion analysis in class-structured populations.
| Appendix |
|---|
|
|
|---|
A. Reproductive values and stable age distribution
Assuming that
remains at 1 due to density-dependent regulation, the vector of individual reproductive values v* = [v
, ..., v
, v
] in the resident population can be determined up to a scalar constant from the recurrence equation v*A* = v*. Because the eigenvector v* can be scaled at will, I choose the reproductive value of parents v
to equal 1. Then, the reproductive value of each offspring class is
|
| (A1) |
Likewise, the vector of stable age distribution u* = [u
, ..., u
, u
] can be determined up to a constant from the recurrence equation u*A* = u*. By setting the relative number of parental class u
at 1, I obtain the relative number of each offspring class:
|
| (A2) |
B. The internal consistency of the resident matrix A
For simplicity, I consider a population in which offspring are classified into 2 age classes. Suppose that the age-structured population has just reached a stable age structure at time t. Let the total number of adults in the population be N (i.e., the population is divided into N broods). The numbers of individuals in age-1 and age-2 offspring classes are Nf and Nf
1, respectively (see Equation A2). Clearly, each age-i offspring in this population will receive xi resources.
At time t + 1, N
adults will survive and continue to provision their offspring. By contrast, N(1 –
) adults will die, and hence, N(1 –
) territories will become vacant. The model reasonably assumes that the death of a parent at any time step will immediately result in the death of all immature offspring in the brood. The newly recruited Nf
1
2 adults, who were age-2 offspring in the previous time step, will occupy the N(1 –
) vacant territories (recall that Nf
1
2 = N(1 –
) due to density-dependent factors, see Equation 2). Each of the newly recruited adults will produce f newborn offspring and provision them. Thus, 2 subpopulations can be identified in the population at time t + 1. Subpopulation A is composed of N
broods, each of which is occupied by an adult who survived from the previous time unit with probability
. The numbers of age-1 and age-2 offspring in subpopulation A are N
f and Nf
1, respectively. Subpopulation B is composed of Nf
1
2 broods, each of which is occupied by an adult who was an age-2 offspring in the previous time step. The number of age-1 offspring in subpopulation B is Nf2
1
2, and there are no age-2 offspring in subpopulation B.
Because the population has already reached the evolutionary equilibrium at time t, the total number of age-2 offspring in the population at time t + 1 should remain at Nf
1, which is actually the number of age-2 offspring in subpopulation A. Likewise, the number of age-1 offspring in the entire population at time t + 1 should remain at Nf, which is the number of age-1 offspring in subpopulation A (=N
f) plus the number of age-1 offspring in subpopulation B (=Nf2
1
2 = Nf(1 –
)). Consequently, each and every age-i (i = 1, 2) offspring in the entire population can be regarded as receiving x
resources from its living parent at any given time step, with the average parental mortality holding constant at 1 –
. The same can be said of a population with n offspring age classes.
C. Evolutionary stability and convergence stability
The second-order condition for the equilibrium point (x
, ..., xn
) to be ES is that the Hessian matrix of
(x°, x*) is negative semidefinite at that point (Day and Taylor 2003
). In order to obtain the derivative of
(x°, x*) with respect to a mutant strategy x°, I use the following equation shown by Taylor and Frank (1996)
:
|
| (A3) |
(x°, x*) with respect to a mutant strategy x° and evaluate at x° = x*:
![]() | (A4) |
Let M be the n x n Hessian, real symmetric matrix of
(x°, x*), with elements
|
| (A5) |
In order for M to be negative semidefinite at the equilibrium point (x1*, ..., xn*), all its eigenvalues should be zero or negative. The Routh–Hurwitz condition for n-dimensional matrices can be used to determine the local stability of the n x n Hessian matrix M. The actual calculation, however, is very cumbersome, and hence, the local maximum at the equilibrium point (x1*, ..., xn*) was numerically verified.
Another stability criterion called convergence stability ensures evolutionary attainability through a series of small steps (Eshel and Motro 1981). Denote by
(=
x° + (1 –
)x*) the mean of parental allocation strategies in a population with a rare mutant strategy x° at frequency
and the ES strategy x* at frequency (1 –
) with
being small. Multidimensional convergence stability requires that the so-called Jacobian matrix of selection gradient, with elements
|
| (A6) |
D. The pattern of ES allocation when
'i(x)>
'i+1(x) for all positive x
Equation 8 reveals that
can be regarded as a function of x
only. Differentiating
(x
) with respect to x
, I obtain the following:
![]() | (A7) |
(x
) is a decreasing function of x
.
Because the survival probability must lie in between 0 and 1, I assume that 
(x) approaches an asymptote of Ki and 
(x) approaches Ki+1 as x approaches infinity (0 < Ki+1 < Ki
1 because
'
(x) >
'
(x) for all positive x). In other words, both
'
(x) and
'
(x) approach zero as x approaches infinity, which implies that
(x
) approaches zero as x
approaches infinity. Because
(x
) is a decreasing function of x
,
(x
) has a positive value for all x
. Thus, x
is lower than x
for all x
.
E. The pattern of ES allocation when
'i(x)>
'i+1(x) only for small values of x
Denote by
the specific value of x at which
'i + 1(x) becomes equal to
'i(x). That is,
'i(x) >
'i + 1(x) for x <
but
'i(x) <
'i + 1(x) for x>
. First, consider x
, the value of which is greater than
. Given that
i(x
) >
i + 1(x
) and
'i(x
) <
'i + 1(x
), I can derive the following inequality:
|
| (A8) |
Hence, x
is lower than x
for x
>
.
Second, consider x
, the value of which is less than or equal to
. The decreasing function
(x
) has its minimum value at x
=
. The minimum value is positive as follows:
|
| (A9) |
Hence,
(x
) is positive, which means that x
is lower than x
for x
. In sum, x
is lower than x
for all x
.
| ACKNOWLEDGEMENTS |
|---|
I thank David Buss, Jae Choe, Troy Day, Jeremy Field, Mark Kirkpatrick, Sang-Im Lee, Olof Leimar, Douglas Mock, Geoff Parker, Ido Pen, Peter Taylor, and 2 anonymous reviewers for helpful comments on earlier versions of the manuscript. I am especially grateful to Susan Lappan for improving my English.
| REFERENCES |
|---|
|
|
|---|
Bize P, Piault R, Moureau B, Heeb P. A UV signal of offspring condition mediates context-dependent parental favouritism. Proc R Soc Lond Ser B (2006) 273:2063–2068.[Medline]
Boland CRJ, Heinsohn R, Cockburn A. Experimental manipulation of brood reduction and parental care in cooperatively breeding white-winged choughs. J Anim Ecol (1997) 66:683–691.[CrossRef]
Charlesworth B. Evolution in age-structured populations (1994) 2nd ed. Cambridge (UK): Cambridge University Press.
Charnov EL. Life history invariants: some explorations of symmetry in evolutionary ecology. (1993) Oxford: Oxford University Press.
Clutton-Brock TH. The evolution of parental care (1991) Princeton (NJ): Princeton University Press.
Cotton PA, Wright J, Kacelnik A. Chick begging strategies in relation to brood hierarchies and hatching asynchrony. Am Nat (1999) 153:412–420.[CrossRef][Web of Science]
Day T, Taylor PD. Evolutionary dynamics and stability in discrete and continuous games. Evol Ecol Res (2003) 5:605–613.
Eshel I, Motro A. Kin selection and strong evolutionary stability of mutual help. Theor Pop Biol (1981) 19:420–433.[CrossRef][Web of Science][Medline]
Fisher RA. The genetical theory of natural selection (1930) Oxford: Oxford University Press.
Frank SA. Foundations of social evolution (1998) Princeton (NJ): Princeton University Press.
Glassey B, Forbes S. Begging and asymmetric nestling competition. In: The evolution of begging: competition, cooperation and communication—Wright J, Leonard ML, eds. (2002) Dordrecht (The Netherlands): Kluwer Academic Publishers. 269–282.
Godfray HCJ. Signalling of need by offspring to their parents. Nature (1991) 352:328–330.[CrossRef][Web of Science]
Godfray HCJ. Signalling of need between parents and young: parent-offspring conflict and sibling rivalry. Am Nat (1995) 146:1–24.[CrossRef][Web of Science]
Haig D. Brood reduction and optimal parental investment when offspring differ in quality. Am Nat (1990) 136:550–566.[CrossRef][Web of Science]
Johnstone RA. Begging and sibling competition: how should offspring respond to their rivals? Am Nat (2004) 163:388–406.[CrossRef][Web of Science][Medline]
Kilner R. When do canary parents respond to nestling signals of need? Proc R Soc Lond Ser B (1995) 260:343–348.[CrossRef]
Kilner RM. Sex differences in canary (Serinus canaria) provisioning rules. Behav Ecol Sociobiol (2002) 52:400–407.[CrossRef][Web of Science]
Kilner RM, Johnstone RA. Begging the question: are offspring solicitation behaviors signals of need? Trends Ecol Evol (1997) 12:11–15.[CrossRef]
Krebs EA, Cunningham RB, Donnelly CF. Complex patterns of food allocation in asynchronously hatching broods of crimson rosellas. Anim Behav (1999) 57:753–763.[CrossRef][Web of Science][Medline]
Lack D. The significance of clutch size. Ibis (1947) 89:302–352.[CrossRef]
Lack D. The natural regulation of animal numbers. (1954) Oxford: Oxford University Press.
Leonard M, Horn A. Provisioning rules in tree swallows. Behav Ecol Sociobiol (1996) 38:341–347.[CrossRef][Web of Science]
Lessells CM. Parentally biased favouritism: why should parents specialize in caring for different offspring? Philos Trans R Soc Lond Ser B Biol Sci (2002) 357:381–403.
Lloyd DG. Parallels between sexual strategies and other allocation strategies. Experientia (1985) 41:1277–1285.[CrossRef][Web of Science][Medline]
Lloyd DG. A general principle for the allocation of limited resources. Evol Ecol (1988) 2:175–187.[CrossRef]
Lotem A. Higher levels of begging behavior by small nestlings: a case of a negatively correlated handicap. Isr J Zool (1998) 44:29–45.[Web of Science]
Macnair MR, Parker GA. Models of parent-offspring conflict III. Intra-brood conflict. Anim Behav (1979) 26:111–122.[CrossRef][Web of Science]
Magrath RD. Hatching asynchrony in altricial birds. Biol Rev (1990) 65:587–622.
Metz JAJ, Nisbet RM, Geritz SAH. How should we define fitness for general ecological scenarios? Trends Ecol Evol (1992) 7:198–202.[CrossRef]
Mock DW, Parker GA. The evolution of sibling rivalry. (1997) Oxford: Oxford University Press.
Mylius SD, Diekmann O. On evolutionarily stable life histories, optimization and the need to be specific about density dependence. Oikos (1995) 74:218–224.[CrossRef][Web of Science]
Otto S, Day T. A biologist's guide to mathematical modeling in ecology and evolution. (2007) Princeton (NJ): Princeton University Press.
Parker GA, Macnair MR. Models of parent-offspring conflict IV. Suppression: evolutionary retaliation by the parent. Anim Behav (1979) 27:1210–1235.[CrossRef][Web of Science]
Parker GA, Mock DW, Lamey TC. How selfish should stronger sibs be? Am Nat (1989) 133:846–868.[CrossRef][Web of Science]
Parker GA, Royle NJ, Hartley IR. Begging scrambles with unequal chicks: interactions between need and competitive ability. Ecol Lett (2002) 5:206–215.[CrossRef][Web of Science]
Pen I, Weissing FJ. Towards a unified theory of cooperative breeding: the role of ecology and life history re-examined. Proc R Soc Lond Ser B (2000) 267:2411–2418.[Medline]
Pen I, Weissing FJ. Optimal sex allocation: steps towards a mechanistic theory. In: Sex ratios: concepts and research methods—Hardy I, ed. (2002) New York: Cambridge University Press. 26–45.
Price K, Harvey H, Ydenberg R. Begging tactics of nestling yellow-headed blackbirds, Xanthocephalus xanthocephalus, in relation to need. Anim Behav (1996) 51:421–435.[CrossRef][Web of Science]
Price K, Ydenberg R. Begging and provisioning in broods of asynchronously-hatched yellow-headed blackbird nestlings. Behav Ecol Sociobiol (1995) 37:201–208.[Medline]
Ricklefs RE. Density dependence, evolutionary optimization, and the diversification of avian life histories. Condor (2000) 102:9–22.
Rodríguez-Gironés MA, Cotton PA, Kacelnik A. The evolution of begging: signaling and sibling competition. Proc Natl Acad Sci USA (1996) 93:14627–14641.
Royle NJ, Hartley IR, Parker GA. Begging for control: when are offspring solicitation behaviors honest? Trends Ecol Evol (2002) 17:434–440.[CrossRef]
Rubenstein DI. Reproductive value and behavioral strategies: coming of age in monkeys and horses. Perspect Ethol (1982) 5:469–487.
Sargent RC, Gross MR. William's principle: an explanation of parental care in teleost fishes. In: The behavior of teleost fishes—Pitcher T, ed. (1986) London: Croom Helm Publishers. 275–293.
Smiseth PT, Amundsen T. Senior and junior nestlings in asynchronous bluethroat broods differ in their effectiveness of begging. Evol Ecol Res (2002) 4:1177–1189.
Smiseth PT, Bu RJ, Erikenaes AK, Amundsen T. Food limitation in asynchronous bluethroat broods: effects on food distribution, nestling begging, and parental provisioning rules. Behav Ecol (2003) 14:793–801.
Smith CC, Fretwell SD. The optimal balance between size and number of offspring. Am Nat (1974) 108:499–506.[CrossRef][Web of Science]
Stamps J, Clark A, Arrowood P, Kus B. Parent-offspring conflict in budgerigars. Behaviour (1985) 94:1–40.[CrossRef]
Taylor PD. Allele frequency change in a class-structured population. Am Nat (1990) 135:95–106.[CrossRef][Web of Science]
Taylor PD, Frank SA. How to make a kin selection model. J Theor Biol (1996) 180:27–37.[CrossRef][Web of Science][Medline]
Temme DH. Seed size variability: a consequence of variable genetic quality among offspring? Evolution (1986) 40:414–417.[CrossRef][Web of Science]
Trivers R. Parent-offspring conflict. Am Zool (1974) 14:249–264.[Web of Science]
West-Eberhard MJ. The evolution of social behavior by kin selection. Q Rev Biol (1975) 50:1–33.[CrossRef]
Winkler DW. A general model for parental care. Am Nat (1987) 130:526–543.[CrossRef][Web of Science]
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
F. Mas, K. F. Haynes, and M. Kolliker A chemical signal of offspring quality affects maternal care in a social insect Proc R Soc B, August 7, 2009; 276(1668): 2847 - 2853. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||








