* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Senescence as an Adaptation to Limit the Spread of
Survey
Document related concepts
Unilineal evolution wikipedia , lookup
Gene expression programming wikipedia , lookup
Sexual selection wikipedia , lookup
Hologenome theory of evolution wikipedia , lookup
The Selfish Gene wikipedia , lookup
Evolution of sexual reproduction wikipedia , lookup
Kin selection wikipedia , lookup
Viral phylodynamics wikipedia , lookup
Natural selection wikipedia , lookup
Theistic evolution wikipedia , lookup
Genetics and the Origin of Species wikipedia , lookup
Evolutionary landscape wikipedia , lookup
State switching wikipedia , lookup
Microbial cooperation wikipedia , lookup
Transcript
5 Senescence as an Adaptation to Limit the Spread of Disease 10 Josh Mitteldorf* and John Pepper† 15 *[email protected] Dept of Ecology & Evolutionary Biology University of Arizona, Tucson AZ 85720 †[email protected] Dept of Ecology & Evolutionary Biology University of Arizona, Tucson AZ 85720 20 submitted to Theoretical Population Biology 25 Running head: The Red Queen evolves senescence 30 Senescence as an Adaptation to Limit the Spread of Disease Josh Mitteldorf and John Pepper Abstract 1 5 10 15 The genetics and phenomenology of senescence suggest the hallmarks of an adaptation; but the effect of senescence on individual fitness is negative, so theorists have preferred to regard it as a non-adaptive byproduct of selection. We ask whether the possibility that senescence is an adaptation in its own right can be confidently excluded on theoretical grounds, and explore the viability of the Red Queen mechanism, transplanted from the evolution of sex. According to this hypothesis, pathogens provide a powerful motivation to maintain immunological diversity, (promoted at the population level by both sex and senescence). We study the behavior of a twodimensional spatial model, in which fatal epidemics spread through contiguous domains of genetically identical individuals. Heterogeneity creates barriers to the spread of disease, and senescence (in the form of programmed death) is observed to evolve on this basis. We conclude that the evolution of programmed death as an adaptation may be a more viable hypothesis than has been supposed, and should not be dismissed. Keywords: aging, senescence, programmed death, epidemics, group selection, red queen, population viscositySenescence as an Adaptation to Limit the Spread of Disease Josh Mitteldorf and John Pepper 20 Natural selection should ordinarily proceed towards lengthening life, not shortening it. Such selection, at the individual level, could conceivably be countered by selection at the population level, if selection somehow favored group-survival. - Williams (1957 p. 399) 25 Introduction In the half century since these words were written, naïve, ill-formed arguments from the good of the species gave way first to a healthy skepticism about group selection, and then to a solid science of multi-level selection, spearheaded by D.S. Wilson (1980; Sober and Wilson 1998), building on the pioneering work of Price (1970, 1972). But explaining organismal senescence as 30 a group-level adaptation remains a formidable challenge, for all the reasons that Williams (1966) appreciated: the costs of senescence are immediate, and directly affect reproductive potential; the benefits are diffuse, and act in the long term to favor broad communities that share a food source, 2 or ecological niche. Why would we not expect long-lived “cheaters” to take over a population of senescing individuals? It has long been believed that aging is a non-adaptive byproduct of natural selection, rather than an adaptation in its own right. Accumulated experimental evidence, however, indicates that 5 senescence appears to be a programmed process that proceeds under genetic controls that have been highly conserved over great spans of evolutionary time (Guarente and Kenyon 2000; Kenyon 2001). Programmed death has long been observed in semelparous organisms (Austad 2004, Steinbehrens and Sapolsky 1992), and the biochemical markers of apoptosis have been observed in populations of yeast under stress (Fabrizio et al. 2004). Protists suffer programmed 10 death in the form of replicative senescence (Clark 1998, 1999). For these and other reasons, it has been proposed that senescence has the hallmarks of an evolved adaptation (Mitteldorf 2004, Longo et al. 2005, Skulachev 1997, Bredesen 2004). Theorists have discounted experimental evidence for the adaptive nature of senescence, dismissing the idea as a logical impossibility. Some of the staunchest advocates for established 15 (pleiotropic) theories have been reluctant to engage in debate on the issue, declaring that it is a settled question. The model presented herein is offered as an existence proof. It is a counterexample to the widely-held view that no plausible evolutionary mechanism could select senescence as an adaptation in its own right. In fact, the present model is one realization of the Red Queen hypothesis, which is generally 20 regarded as one plausible explanation for the evolution of sex and the cost of males. According to the Red Queen hypothesis, fast-evolving parasites provide ongoing evolutionary pressure for a population to sustain genetic diversity. Sexual reproduction contributes to this end by mixing up 3 genomes; and senescence can also contribute to the same end by increasing the population turnover rate. The individual cost of maintaining males in a population can be as high as a factor of two in fitness; the individual cost of senescence is likely to be substantially lower. If we accept the Red Queen hypothesis as a viable explanation for the evolution of sex, why should we 5 not also regard it as a possible explanation for the evolution of senescence? Historical Review Strong between-group selection is required to overcome the self-evident within-group disadvantage. There is a hundred-year history of mechanisms proposed to meet this challenge, 10 but much of it is limited by the assumption (taking various forms, sometimes implicit) that individuals are inevitably weakened by age. Candidate mechanisms from the past have proposed that senescence exists for the purpose of eliminating damaged specimens, making room in the niche for younger, undamaged individuals. The first such model was proposed by Weismann (1889). 15 Suppose that an immortal individual could escape all fatal accidents, through infinite time – a supposition which is of course hardly conceivable. The individual would nevertheless be unable to avoid, from time to time, slight injuries to one or other part of its body. The injured parts could not regain their former integrity, and thus the longer the individual lived, the more defective and crippled it would become, and the less perfectly would it fulfill the purpose of its species. Individuals are injured by the operation of external forces and for this reason alone it is 20 necessary that new and perfect individuals should continually arise and take their place, and this necessity would remain even if the individuals possessed the power of living eternally. From this follows, on the one hand the necessity of reproduction and on the other the utility of death. Worn out individuals are not only valueless to the species, but they are even harmful, for they take the place of those which are sound. Hence by the operation of 4 natural selection, the life of our hypothetically immortal individual would be shortened by the amount which was useless to the species. This mechanism relies explicitly on population level selection, and Weismann did not here anticipate the difficulties associated with competition between different levels of selection. A 5 deeper logical flaw, which ultimately led Weismann to abandon this idea, is that there is no necessity for older individuals to be “defective and crippled.” If an organism is able to build itself entirely from seed, why should the capacity for regeneration and repair process be limited at all? Why should it be less costly to build a new individual from scratch than to repair injuries that may have accumulated in an older individual? 10 Starting in the middle of the 20th century, thinking about mechanisms of evolution was channeled toward exclusively individual selection, with individual reproductive success taken as the quantitative realization of Darwinian fitness. Within this paradigm, it is not possible for senescence to evolve as an adaptation; hence non-adaptive theories rose to the fore. Modern theories for evolution of senescence are rooted in the ideas of Medawar (1952) concerning the 15 declining force of natural selection with age, and Williams (1957), who proposed that organismal senescence was a maladaptive by-product of individual selection for high survival and fertility early in life. Empirical evidence discrediting the theories of aging as a maladaptive side-effect of evolution has motivated recent proposals for modeling aging as an adaptation. Libertini (1988) proposed 20 that the evolution of adaptive senescence may be possible if individuals dying of senescence are replaced by kin. He did not, however, specify a mechanism by which this might come about. Kirchner and Roy (1999) explore a model which is, like the present model, based on the utility of 5 senescence in avoiding epidemics. A key assumption of that model is that the disease being modeled is impose permanent sterility when it does not kill its victim, thus replacing a fertile member of the population with a dysfunctional individual that drains limited resources from the reproducing population. Were this to be proposed as a general mechanism for evolution of 5 senescence, it would suffer from the same essential limitation as the Weismann hypothesis: Why should it be less costly to create a new individual from seed than to restore those reproductive organs that were damaged by disease? Similarly, Travis (2004) models a mechanism for the evolution of programmed death as an adaptation, given that individual fecundity declines with age. Older individuals are eliminated in 10 favor of their neighbor kin because they are using communal resources, while not contributing reproductive output. Dytham and Travis (2006) further explore this model. The most limiting assumption has been relaxed somewhat in an updated model that substitutes a rapidly-changing environment, to which older individuals are progressively less well adapted (Travis et al 2006). The ubiquity of senescence suggests that its evolution should not require special or unusual 15 conditions. In searching for a general group-level force of sufficient intensity to compete with the individual cost of senescence, we are inspired by a growing recognition of the strength of population dynamics as a factor in natural selection (Gilpin 1975, Pepper & Smuts 2000; Pels et al 2002). In contrast to other group-level advantages, ecological interactions can be swift and deadly in their effects on a population. One of us (Mitteldorf 2006) has proposed a model for 20 evolution of senescence in a predator species, based on the ability of predator senescence to damp population fluctuations. Here, we suggest a different mechanism, based on host-pathogen population dynamics, and aimed at explaining senescence in the host (prey) species. 6 We are inspired by the parallel problem of the origin and maintenance of sexual reproduction, which also requires strong group selection by most accounts. The primary theoretical advantage of sex, due originally also to Weismann (1889), is that sexual recombination provides new genetic variation for selection to act upon. This idea was criticized by Williams (1966, 1975), as 5 being based upon group (or population) selection. However, many evolutionary theorists who remain skeptical of the potency of group selection in other contexts maintain an exception for evolution of sex. Some forms of group selection are more rapid and potent than others. In general, benefits to evolvability (Bell, 1982; Burt 2000) are slow, requiring an evolutionary time scale in order to 10 manifest. Over so many generations, the benefit is diffused, and the “cheater” problem becomes decisive. Population dynamics, in contrast, can include swift and sudden local extinctions. An adaptation that helps avoid such extinctions offers an advantage that can be competitive with individual reproductive value in its force (Mitteldorf 2006). This explains the strength of the Red Queen Hypothesis. Bell (1982) posits that in an 15 evolutionary arms race against parasites and pathogens, rapid evolutionary change is at a premium, and consequently the genetic variability promoted by sexual recombination affords a decisive advantage. (The phrase derives from an analogy by van Valen (1973) to the Red Queen in Charles Dodgson’s classic children’s story (Carroll 1871).) The appeal of the Red Queen hypothesis for the evolution of sexuality led us to wonder whether pathogen pressure could 20 similarly favor senescence, which also contributes to genetic diversity. Population dynamics are affected by spatial structure. Physical proximity and a common ecosystem provide the connections that link the fates of related altruists in predator/prey or 7 host/pathogen systems. Previous simulation studies have demonstrated that physical proximity can be a proxy for genetic relatedness, promoting the evolution of blind, local altruistic traits (Axelrod 1997, Mitteldorf and Wilson 2000, Dytham and Travis 1998, Pepper& Smuts 2001). When offspring are born and remain near their parents, spatial kinship correlations result. If 5 altruistic traits are directed toward spatial neighbors, then the beneficiaries have enhanced likelihood of carrying the relevant altruistic trait. The resulting inclusive fitness benefits of altruism can effectively oppose the within-group selective advantage of selfishness. But this is inclusive fitness in a broad sense, where relatedness cannot easily be calculated by descent (Hamilton 1964), or predicted analytically, but only measured empirically as a regression 10 coefficient (Hamilton 1970, 1972). Population structures: group selection and kin selection Hamilton’s (1964) theory of inclusive fitness has become a fundamental point of reference in explaining the evolution of altruism. Hamilton’s Rule states that a trait can evolve so long as its 15 benefit B and its cost C are related by rB > C. The r in this inequality stands for “relatedness”, and it is often modeled as a coefficient that derives from common descent. But, more generally, it is a correlation coefficient, related to the excess probability that the altruistic gene is shared by the beneficiary of the altruism. In cases that involve explicit kinship, Hamilton’s Rule is easily applied. Examples include risks 20 that parents incur to promote the survival of their offspring, and cooperative behavior within a colony of eusocial insects. But altruistic benefit may also be dispersed more widely in a geographic region, mediated by a common food stock or a common pool of parasites. Within a 8 compact geographic region, individuals are inbred and likely to share genes in common, but their kinship may be difficult to trace. In such circumstances, Hamilton’s Rule remains valid, but it has little predictive value when the value of r cannot be calculated a priori (Mitteldorf and Wilson 2000). Predator/prey and parasite/host dynamics are famously nonlinear, and here it may 5 be especially true that Hamilton’s Rule applies, but offers little predictive power. Whether such situations of spatially mediated genetic similarity are described as ‘group selection’ or ‘kin selection’ is a matter of taste (Pepper 2000). The formal equivalence of Hamilton’s (1964) inclusive fitness theory, and Price’s (1972) multilevel selection theory has been noted by various authors (Hamilton 1975; Wade 1980; Queller 1992), leading Bell (1997 p. 10 529) to remark that, “The distinction between group selection and kin selection, though conventional, seems somewhat artificial.” In the present study, we model the evolution of senescence as a promoter of diversity in a spatial array of organisms subject to infectious epidemics. From the literature on the evolution of sex, we borrow the idea that many organisms may be locked in a perpetual evolutionary arms race 15 with their parasites and pathogens. Because pathogens evolve rapidly, their biochemical attacks are tuned to the phenotype of the host, so that the host population is under constant selection pressure to change for change’s sake. Hosts typically have longer generations and slower evolution than their pathogens, so offspring will tend to be more resistant to disease than their parents, and therefore to have higher reproductive value. The idea that parents might be selected 20 to sacrifice themselves for the good of their offspring could be framed in terms of classical kin selection. However, in a spatially distributed population such as we model, ongoing population dynamics create a more diffuse distribution of relatedness, where the beneficiaries of altruistic death may include a range of near and more distant kin. The possibility of far-reaching contagion 9 greatly magnifies the communal advantage of altruistic senescence, but also makes its magnitude much more difficult to predict. In our model, when an individual dies, it is replaced by the offspring of a neighbor. Limited dispersal makes it likely that the neighbor is a close relative; but relatedness in this model is a 5 double-edged sword. Senescence as altruism can only evolve if the neighbor who benefits from the altruism is likely to share the altruistic (life span) gene; but the benefit of altruism accrues in the form of genetic diversity (at the locus that determines susceptibility to disease). The mechanism succeeds to the extent that, when an individual dies of old age, he is replaced by another individual who is likely to share his lifespan gene, but likely to have a variant 10 susceptibility gene. Hence, a principal assumption of our model is that we have separated the mutation probabilities at these two loci. Our model succeeds in regimes where mutations to the lifespan gene are rare compared to mutations in the susceptibility gene. There is good evidence that in nature, mutation probabilities vary across the genome, so diversity is promoted where beneficial, even as core systems are conserved (Wright 2004). Although our model is based on 15 clonal reproduction, the inclusion of sex in the model could conceivably strengthen it, since every birth provides a new opportunity for genetic diversity. A high rate of population turnover can help maintain the diversity that makes a population more resistant to epidemics. This is one part of the basis on which senescence is selected in our model; the other piece is that death itself keeps the grid porous, and epidemics propagate less widely at 20 the lower population density. 10 Description of the Model We have adapted a computational model developed by Socolar et al (2001). Sites are arrayed on a two-dimensional, n*n grid, with opposing edges identified to form a torus. Each site may be vacant or may be occupied by a single model organism. In each computational 5 cycle, a site is chosen at random, and if the site is occupied, one of three events may occur: (1) reproduction, (2) death, or (3) the origin of an epidemic. Reproduction occurs at a random neighbor site, creating a viscous population structure with relatedness that varies with physical proximity. Death may occur either because of background mortality, or because the organism’s genetically programmed life span has expired. Epidemics spread through a contiguous pathway 10 of occupied sites, where the occupants are genetically identical at the relevant locus. Epidemics are the source of the model’s highly non-linear behavior, generating selective effects that promote cooperation. The probabilities that govern the selection of one of these three events determine the dynamics of the model. The probability b of reproduction was controlled by a gene in some runs, but after we 15 determined that it tended to evolve to a maximum, we fixed b such that every selected organism which does not die or originate an epidemic is given an opportunity to reproduce. Similarly, we experimented with the background, random mortality rate m as a gene, and found that it evolves to zero when a lifespan gene was available as an alternative, better-regulated source of mortality. The probability of an epidemic is a parameter of the model, fixed as constant within any given 20 run. 11 Thus in most of our runs, the organisms’ genomes included just two loci: The first controls lifespan, implemented in our simple model as a threshold function. The lifespan gene is represented as a floating point number. An organism will vacate its site at the first opportunity after its age exceeds its genetically-programmed lifespan. The second is the susceptibility gene 5 described above, which implies susceptibility to a perfectly-matched pathogen, and is otherwise neutral. Epidemics kill the originator, and then explore all neighbor sites; if a genetically identical susceptibility gene is found, the epidemic kills the neighbor and continues to spread recursively through the grid. 10 (Evolution of pathogens is not modeled explicitly, but novel pathogens are assumed to appear stochastically and to be adapted to match the susceptibility allele of the organism in which they first appear.) The susceptibility gene is an 8-bit integer, and one bit will flip at random with a low probability, specified by the bit mutation rate. In variations of the model, less than 8 bits are used, so that genome space is effectively smaller. In some runs, mutations were disabled (“zero- 15 bit genome”), forcing a sparsely-populated grid, with shorter life spans. The lifespan gene is a floating-point number which mutates with a probability controlled by the lifespan mutation rate. Mutation multiplies the parent’s value by a random number distributed in a logarithmically uniform way from e – ½ range to e + ½ range, where the mutation range typically=0.03. We refer to n2 such cycles as one “time step”. Because each cycle begins with an independent, 20 random determination of a target site, some sites will not be active in any given time step, while other sites will have two or more opportunities to reproduce. On average, each site will be given one opportunity to reproduce or die or begin an epidemic in each time step. 12 Dynamics of the Model Beginning from any distribution of genotypes, the grid tends to self-organize into genetically 5 homogeneous domains, because the rule for reproduction keeps offspring close by. Epidemics take down a domain at a time. The larger a domain, the more likely it is to be wiped out by an epidemic. When either the background mortality or the mutation rate is high, domains can break apart, forming domain boundaries that limit the spread of epidemics. After a domain succumbs to an epidemic, it is re-filled from neighboring domains. If the 10 mutation rate is high, the grid will maintain sufficient diversity to avoid extinction; but the lower the mutation rate, the more turnover the population requires in order to sustain a level of diversity that will protect against epidemics. When mortality from background death rate plus programmed lifespan is too high, the population cannot sustain itself; but when mortality from these two sources is too low, the 15 population grows too dense and too homogeneous; each epidemic then takes down a larger domain. This, too, can lead to extinction. When either of the two sources of mortality is allowed to evolve, mortality is pushed lower by individual selection. But lower mortality in a region then triggers larger epidemics, creating countervailing selection pressure from a higher level. Depending on model parameters, the 20 result can be a stable balance at a modest mortality level, or wide fluctuations in both mortality and total population, or fluctuations that are wider still can drive the population to extinction. 13 Experiments with the Model Senescence vs Fertility limits vs Background Mortality The suppression of epidemics requires limits to population density that can only be achieved at a 5 cost in individual reproductive fitness. Our thesis is that this dynamic supplies a mechanism by which senescence may evolve. But there are other ways that population density may be limited, and epidemics thus suppressed: Individual fertility may be moderated; or the age-independent component of mortality (background mortality) may be enhanced. We conducted numerical experiments pitting senescence against each of these other two life history traits. 10 In the first set of experiments, we included a gene for fertility, defined as the probability that an organism will reproduce in the cycle when it is randomly selected. When the rates for senescent and background mortality were fixed at low levels, the fertility gene evolved a steady state value <1, successfully stabilizing the grid against epidemics. However, when fertility and life span were assigned separate genes, and permitted to evolve simultaneously, we found that fertility 15 evolved to its maximum value (=1) and that the life span gene was left to shoulder the entire burden of moderating population density and controlling epidemics. This was true even when the experiment was run with a single genotype, so that the enhanced diversity provided by a higher fertility rate was not a factor. Similarly, we included a gene for background mortality¸ defined as the probability that an 20 organism would die (regardless of its age) in the cycle when it is randomly selected. (This was the mode of death in the original model of Socolar et al (2000)) Again, when this was the only 14 option available for controlling epidemics, the mortality gene evolved to a functional finite value; however when both lifespan and mortality were permitted to evolve simultaneously, the mortality gene evolved to a low value, consistent with mutational load. Any of these three mechanisms was capable of generating the population turnover to support 5 diversity and limit the spread of epidemics. Given a choice, selection maintained a limited lifespan, and fixed the age-independent death rate at zero and fertility at its maximum value. This suggests that, at least within the assumptions of our model, a programmed life span provides special advantages as a mechanism for stabilizing population size and providing diversity. 10 In response to this result, we fixed fertility at 1 and fixed mortality at 0 in most subsequent experiments. (We did not always set extrinsic mortality to zero in deference to biological reality.) Does mutation dynamics affect evolution of life span? Two parameters of the model controlled the probability of mutation in lifespan, and the scale of 15 the distribution of mutated lifespans. A robust model should produce results independent of the particulars of inheritance, including mutation rate and mutation size. To check that our results did not depend unduly on the mutation dynamics, we allowed the scale of mutations to vary, and charted the lifespan that evolved in response. 15 Figure 1: Size of mutations at the life span locus has a modest effect on the life span that evolves. It also affects the “focus” of selection, determining how narrowly the evolved life span is defined. After experimenting with this parameter, we selected 0.03 as a mutation range for the remaining numerical experiments, being the value that best sharpens the effects of selection. The reported “range” of life span mutations scaled the (logarithmic) size of mutations at the lifespan locus. Mutation would produce an offspring between e – ½ range and e + ½ range times the parent’s lifespan. For computationally reasonable values of the range variable, the evolved 5 lifespan did not vary unduly with mutation dynamics. For small and especially for large values of the range, the (experimental) error bars on evolved lifespan were seen to widen. We chose range=0.03 as the value for which evolution seemed to be most efficient, in that the experiments produced a consistent result most reliably. 16 Parameters for this experiment were: 5 10 512*512 grid Lifespan mutation probability = 0.1 (with variable range) Epidemic Probability=3.16*10–5 Single genotype (no mutation at the susceptibility locus) Fertility = 1 Extrinsic mortality=0 Evolution of life span as a function of extrinsic mortality Traditional theories for the evolution of aging predict that high extrinsic mortality should be associated with the evolution of a shorter intrinsic lifespan. This prediction has been borne out in several studies (e.g. Austad 1993; Blanco and Sherman 2005; Dudycha and Tessier 1999; Ricklefs 1998; Stearns 2000). Intriguingly, there are circumstances – perhaps common in nature 15 – where the prediction appears to fail, and lower extrinsic mortality has the paradoxical result of promoting the evolution of a shorter intrinsic life span (Reznick et al 2004; Promislow 1991). Indeed, Ricklefs (1998) reports that the proportion of deaths attributable to senescence in the wild is negatively correlated with extrinsic mortality, in a review of bird and mammal field studies. The evidence may be strong enough to justify a radical hypothesis, based on population 20 dynamics rather than traditional population genetics: perhaps senescence is playing a complementary role, contributing to population homeostasis by maintaining the overall death rate in the face of fluctuating incidental mortality. The present experiment was designed to determine whether senescence in our model exhibits this behavior. We prescribed an extrinsic death rate globally, varying this parameter from run to run, and asked 25 what life span would evolve in response. The susceptibility gene and the gene for lifespan were the only two loci in this version of the model. 17 Parameters for this experiment were: 5 10 512*512 grid Lifespan mutation probability = 0.1 (with variable range) Epidemic Probability=3.16*10–5 8-bit susceptibility genotype, with mutation probability = 3.16*10–5 Fertility = 1 Extrinsic mortality varies Each run of the model comprised a “throwaway” period of 200 time steps, in which the imprint of initial conditions was permitted to fade; then the model continued to run until 108 total deaths were accumulated, and statistics were compiled from this sample. In Figure 2, two measures of evolved life span are plotted. Life span measured in time steps, is a model-dependent measure, the significance of which is relative to the particular parameters of 15 each run. The percentage of deaths attributable to senescence is a dimensionless, modelindependent measure of the fitness cost that senescence imposes. This latter measure has been computed for natural populations by several authors (Bonduriansky and Brassil 2002, Ricklefs 1998). Life span evolved to be longer when the extrinsic death rate was prescribed to be higher; hence the proportions of deaths from each of these causes had a complementary relationship. 18 Figure 2: When background mortality is low, there is more selection pressure to shorten life span (in order to limit epidemics). Senescence complements background mortality, contributing to population homeostasis. (The fact that senescent deaths appear to level off at a low value may be an artifact of the high mutation rate – the model equivalent of ‘mutational load’ as deleterious mutations shortening life 5 span accumulated faster than selection removed them.) We see that an extrinsic death rate that is sufficiently high obviates selection for limits on life span. For lower extrinsic mortality, the life span tends to adopt a complementary value, such that the total burden from senescence and background death rate is approximately constant. The proportion of deaths from the two causes combined decreased slowly as extrinsic death rate 10 increased. 19 Evolved limit to life span complements background mortality Evolved life span 1.4% 30 1.2% 1.0% 20 0.8% 0.6% 10 0.4% 0.2% 0 0.0% 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 Background mortality Life span % Deaths fr Senescence Figure 1: The X axis represents the prescribed rate of extrinsic senescence, as a probability per time step. Plotted against the Y axis is the proportion of deaths from aging and from the prescribed extrinsic mortality. (Epidemics accounted for the residual deaths.) A complementary relationship is seen, in agreement with the review by Ricklefs (1998). Evolved life span as a function of epidemic rate. Another model parameter that affects the evolved life span is the frequency with which novel pathogens are introduced (‘epidemic rate’). Frequent epidemics offer good feedback, so that life 5 span limits can evolve reliably; infrequent epidemics may fail to punish evolved longevity until it is too late, and extinction can result. More frequent epidemics evolved a shorter life span, but a smaller proportion of all deaths were attributable to senescence. Epidemics that were less frequent resulted in longer evolved life spans, but senescence claimed a larger percentage of all deaths. 10 This effect can be seen in the plot displayed as Figure 3. Once again, the statistics are compiled from an average of the first 108 individual deaths after an initial ‘throw-away’ period of 200 time steps. 20 Parameters for this experiment were: 5 512*512 grid Lifespan mutation probability = 0.1 with range=0.03 Epidemic probability varies 8-bit susceptibility genotype, with mutation probability = 3.16*10–5 Fertility = 1 Extrinsic mortality =0 21 22 Figure 3AB: As epidemic rate varies over two orders of magnitude, the limits to life span that evolve remain within a narrow range. The proportion of deaths attributable to senescence, however, varies widely. When epidemics are less frequent, senescent deaths become relatively more important. Senescence may be seen as a homeostatic mechanism regulating population density. Effect of the size of genotype space In the original model of Socolar et al (2000), there was a single genotype. Working with the 5 model in just this form, we were able to verify, given a lifespan gene, that senescence evolves reliably, and that selection prefers senescence to two other epidemic control mechanisms. But in the real world, epidemics are limited not only by limits on population density, but also by diversity of immune genotypes. (This is the essence of the Red Queen hypothesis.) To test this idea, we modified the model, allowing for 2n immune genotypes (“n-bit susceptibility gene”), 10 and specifying that epidemics will not spread between individuals with different susceptibility genotypes. Results for different values of n are plotted in Figure 4. Parameters for this experiment were: 15 512*512 grid Lifespan mutation probability = 0.1 with range=0.03 Epidemic probability = 3.16*10–5 Variable n-bit susceptibility genotype, with mutation probability = 3.16*10–5 Fertility = 1 Extrinsic mortality =0 23 Figure 4: If there is only a single immune genotype (0 genes) then epidemics must be controlled solely by interspersing voids between population patches. Life span evolves severe limits to accomplish this. As the number of genes increases, genetic diversity provides another mechanism whereby epidemics may be controlled, and larger genotype space permits longer life spans to evolve. With only a single genotype (n=0), the evolved limit to lifespan was strongest. More susceptibility genes implied more genotypes on the grid at any given time, so that epidemics could be avoided with weaker limits on lifespan. As n was increased and evolved lifespan grew, 5 the range of evolved lifespans expanded as well. The dynamics became more volatile (as measured by the time variance of population). Extinctions were a more frequent outcome, thus the range of parameter space in which we could gather data about evolved lifespan became narrower. In particular, higher rates of mutation in combination with large n led to instability and frequent extinctions. 24 Evolved life span as a function of mutation rate at susceptibility locus In our model, limits on life span evolve to enhance population turnover rate, supporting diversity. The product of mutation rate and population turnover rate is the source of that 5 diversity, so that the higher the mutation rate, the longer the life span that can be tolerated. This effect is seen in the unconnected square markers of Figure 5. In the same figure, the proportion of deaths attributable to senescence (solid line) is seen to decline as mutation rate increases. Both these measures vary within narrow limits, even as prescribed mutation rate is varied over 3 orders of magnitude. This is a tribute to the highly non-linear dynamics intrinsic to the model. 10 15 Parameters of these runs were as follows: 512*512 grid Lifespan mutation probability = 0.1 with range=0.03 Epidemic probability = 3.16*10–5 4-bit susceptibility genotype, with variable mutation probability Fertility = 1 Extrinsic mortality =0 25 Figure 5: Mutations at the susceptibility locus supply diversity which limits the spread of epidemics. The rate of such mutations was taken as a parameter, and lifespan was evolved in response to varying mutation rates. Lower mutation rates required a faster population turnover (shorter lifespan) to provide the same protection from epidemics (un-connected squares). The proportion of deaths attributable to senescence rose with shorter lifespan (downward-sloping line . Changes in both measures of mortality were modest, however, in comparison with the 1000-fold change in prescribed mutation rate. Effect of dispersal range Dytham and Travis (2006; Travis and Dytham 1999) have studied the interaction between 5 dispersal and evolution of cooperation, showing that altruism that is blind and local can evolve more easily when dispersal range is limited. When kin remain tightly clustered, the probability that the beneficiary of altruism will be another altruist is enhanced. In our model, dispersal distance corresponds to the size of the neighborhood in which offspring are distributed; hence 26 general considerations suggest that, in our model, smaller neighborhood size should lead to stricter limits on evolved lifespan. There is a second effect of dispersal range in our model: epidemics spread far more easily with the wider neighborhoods, and a corresponding definition of contiguity. In fact, the larger 5 neighborhoods sharply limit the population density that is robust to epidemics, facilitating the evolution of life span limits that is the subject of our investigation. All the data reported above were derived with nine-point neighborhoods, in the form of a 3*3 square. In order to study the effect of dispersal, we ran variants of the model with smaller crossshaped (von Neumann) neighborhoods of size 4 and larger 25-point neighborhoods (5*5 square). 10 When the neighborhood is expanded for the purpose of offspring dispersal, but epidemics are limited to spreading via von Neumann contiguity, the effect is to permit longer individual life spans to evolve. 27 Figure 6: Smaller neighborhood size affects dispersal rate, which should facilitate evolution of limits on life span. But neighborhood size also affects transmission of epidemics, and this turns out to be a stronger effect. Shorter life spans evolved with 9-point square neighborhoods than with 4-point von Neumann neighborhoods. For yet larger 25-point square neighborhoods, evolved life span was shorter yet, but frequent extinctions frustrated our data collection. The bottom line in the plot was extrapolated from incomplete computational results. Discussion The phenomenology of aging speaks of adaptation, and cries out for an evolutionary explanation. 5 Since aging offers the individual only costs and no benefits, the mechanism must depend on a higher level of selection. The ubiquity of aging in nature suggests a general evolutionary basis, 28 and we take this as justification for a model which spares detail. We have endeavored to capture the essence of a mechanism by which natural selection might have created organisms that are programmed to weaken and die. Sustainable population density has some finite value, and catastrophic, non-linear effects can 5 take down an entire local population when this limit is exceeded. This dynamic is a natural foil to the central mechanism of individual selection, i.e., the generational increase in gene frequency for any trait that makes an organism individually more competitive. The intensity of competition for individual reproduction is blunted, and this creates an opening for features such as sex and aging that carry substantial costs for individual reproductive fitness. 10 A demographic theory of aging But this alone cannot explain why aging has been so widely adopted. There are many ways that nature may hold population densities in check, and of these the simplest and most efficient is reproductive restraint. A promising result of our model is that programmed death at a specific 15 age is robustly selected in preference to either lower birth rate or higher (age-independent) death rate. In the original model of Socolar, et al (2000), it was shown that either lower birth rate or higher death rate was capable of stabilizing a population against epidemics. Our results demonstrate that programmed life span offers a superior combination of individual fitness and protection against epidemics than can be obtained with either of these other life history factors. 20 We speculate that this may reflect a property of the real biosphere as well, and may be a basis for natural selection in favor of aging. The advantage that both sex and aging offer over the simpler strategy of reproductive restraint is that they contribute as well to intraspecific diversity. It may 29 seem counter-intuitive that selection could favor genes that kill their bearers, simply to keep them form competing with kin that also carry the same allele and that will be better at replicating it; however, this phenomenon is well documented in other contexts (e.g., Hurst 1991, Travis 2004). 5 The advantage of aging may be understood qualitatively as “taming” stochastic mortality, and bringing it under genetic control. We imagine that local extinction has been an ever-present feature of species evolution, resulting in substantial selection pressure (at the group level) for demographic homeostasis (Mitteldorf, 2005). Epidemics are certainly an extreme example of disruptive patterns of death, but starvation and predation can also impose volatility on a 10 population. A programmed maximum life span offers a way for populations to take control of their demographics, albeit at substantial individual cost. The gradual weakening characteristic of senescence provides a further benefit: whatever environmental hardships are imposed on a population, weaker (older) individuals are likely to succumb first, leaving the younger and stronger ones to replenish the population. This leads to a 15 far more stable demographic than if all adults evinced equal (maximal) resistance to the hardship. This effect was not modeled herein, but suggests an extension of the present model in which epidemics attack a population with graded resistance. In real biological communities, epidemics frequently cull the old and the weak, while buttressing the immunity of the young and strong. 20 Diversity as a target of natural selection The major adaptive advantage of diversity accrues over the long haul, making a population more adaptable in the face of environmental change, and promoting the rate at which new adaptations 30 accrue in the genome. This advantage is far too slow to compete with the direct costs of reduced reproductive value. It is the brilliance of the Red Queen hypothesis that this objection is circumvented: fast-evolving pathogens provide an incentive for population diversity that is immediate and salient. We suggest a picture in which the Red Queen mechanism opens an 5 opportunity for aging and sex to evolve, but that they are maintained in the genome and shaped over the very long term by the contribution they make to the rate of evolution. Senescent death complements stochastic mortality One of our results (Fig. 2) suggests that selection may raise the level of senescent deaths when deaths from extrinsic causes are reduced, resulting in homeostasis with a constant total death 10 rate. Reznick et al (2004) have found that guppies evolved in Trinidad river pools with low predation rates have evolved more rapid senescence, compared to nearby sites where predators dominate the death rate. This phenomenon is difficult to explain in the context of standard theories for evolution of senescence (Mitteldorf and Pepper 2007), and may point to an adaptive role for aging, and a preferred death rate >0 that is optimized under genetic control. 15 In fact, Bryant and Reznick (2004) report anecdotal evidence that when guppies that have evolved high fertility and low rates of aging in high-predation sites are introduced into sites where predators are absent, epidemics of infectious disease have swept through the population. This may be a real-world example of the mechanism we have modeled herein. Conclusion 20 It may be considered a historic artifact that the Red Queen hypothesis enjoys wide acceptance as an explanation for evolution of sex, but has not previously been considered as an explanation for 31 aging. Both sex and aging contribute to the diversity that is essential to a population’s ability to evolve and adapt in the long run, and both exact a substantial cost at the individual level. The first theories of aging to gain wide acceptance invoked only individual selection, a properly conservative approach. Even as advances in genetic methodology and related experimental 5 techniques in the late 20th century have contradicted the principal predictions of these theories, the Red Queen, which is considered quite viable in another context, must surmount hurdles in our ingrained habits of thought in order to gain credibility as an explanation for aging. The prevalent paradigm for these many generations has been that natural selection is fundamentally a race to create more offspring, faster. This picture is incapable of 10 accommodating either the two-fold cost of males or the evolution of senescence as an adaptation. We suggest a new paradigm, in which demographic volatility presents a threat of local extinction that is sufficiently potent to compete head-on with individual selection for fertility. Demographic homeostasis becomes a primary target of natural selection, which necessarily limits individual reproductive potential. We see sexual reproduction, aging and other life history 15 traits as arising from an optimization process which opposes individual selection for reproductive success and selection at the community level, acting via local extinctions. References 20 AUSTAD, S. 1993. Retarded senescence in an insular population of Virginia opossums (Didelphis virginiana). J. Zool. Lond. 229: 695-708. AUSTAD, S. 2004. Is aging programmed? Aging cell 3:249-251 AXELROD, R. 1997. The complexity of cooperation. Princeton U Press. Princeton, NJ. 25 BELL, G. 1982. The Masterpiece of Nature: The evolution and genetics of sexuality. Croom Helm, London. 32 BELL, G. 1997. Selection: The Mechanism of Evolution. Chapman and Hall, New York. BLANCO, M.A. & SHERMAN, P.W. 2005. Maximum longevities of chemically protected and nonprotected fishes, reptiles, and amphibians support evolutionary hypotheses of aging. Mech Ageing Devel 126, 794-803. 5 BRONIKOWSKI, A.M. & PROMISLOW, D.E.L. 2005. Testing evolutionary theories of aging in wild populations. TREE, z, 271-273. BONDURIANSKY, R. AND BRASSIL, C. E. 2002. Senescence: rapid and costly ageing in wild male flies. Nature 420:377 10 BREDESEN, D.E. 2004. The non-existent aging program: how does it work? Aging Cell 3:255259 BRYANT, M. J. AND REZNICK, D. 2004. Comparative studies of senescence in natural populations of guppies. Am Nat 163:55-68 BURT, A. 2000. Perspective: sex, recombination, and the efficacy of selection – was Weismann right? Evolution 54(2):337-351 15 CARROLL, LEWIS [REV. C.L. DODGSON]; Through the Looking Glass and What Alice Found There, Macmillian & Co., 1871. CLARK, W. R. 1998. Reflections on an unsolved problem of biology: the evolution of senescence and death. Unpublished manuscript available at http://www.wrclarkbooks.com/evolution.html 20 ________. 1999. A means to an end: the biological basis of aging and death. Oxford University Press. New York. DUDYCHA, J.L. & TESSIER, A.J. 1999. Natural genetic variation of life span, reproduction, and juvenile growth in Daphnia. Evolution 53:1744-1756. DYTHAM, C & TRAVIS, J.M.J. 2006. Evolving dispersal and age at death. Oikos 113:530-538 25 FABRIZIO, P., BATTISTELLA, L., VARDAVAS, R., GATTAZZO, C., LIOU, L., DIASPRO, A. ET AL. 2004. Superoxide is a mediator of an altruistic aging program in Saccharomyces cerevisiae. J. Cell Biol., 166:1055–1067. GILPIN, M.L. 1975. Group Selection in Predator–Prey Communities. Princeton, NJ: Princeton University Press. 30 GORDEEVA, A. V. ET AL. 2004. Apoptosis in unicellular orgnims: Mechanism and evolution. Biochemistry (Moscow) 69(10): 1055-1066 GUARENTE, L AND KENYON, C. 2000. Genetic pathways that regulate ageing in model organisms. NATURE 408:255-262 HAMILTON, W.D. 1964. The evolution of social behavior. J Theor Biol 1:295-311. 35 ____ 1975. Innate social aptitudes of man: an approach from evolutionary genetics. In Biosocial Anthropology ed R. Fox. London, Malaby Press: 133-135. HAMILTON, W. D., AXELROD, R, AND TANESE, R. 1990. Sexual reproduction as an adaptation to resist parasites (A Review). Proc Nat Acad Sci 87(9): 3566-3573. HARDIN, G. 1998. The tragedy of the commons. Science 162, 1243-1248. 33 HURST, L. D. 1991, The incidences and evolution of cytoplasmic male killers. Proc. Royal Soc. Lond. B. 244(1310): 91-99 KENYON, C. 2001. A conserved regulatory system for aging. Cell 105:165-168 5 KIRCHNER, J.W. & ROY, B.A. 1999. The evolutionary advantages of dying young: epidemiological implications of longevity in metapopulations. Am Nat 154(2):140-159. LIBERTINI,G. 1988. An adaptive theory of the increasing mortality with increasing chronological age in populations in the wild." J. Theor. Biol. 132(2):145-162. MEDAWAR P. B. 1952. An unsolved problem of biology. London. H. K. Lewis. MITTELDORF, J. 2004. Ageing selected for its own sake. Evol. Ecol. Res., 6:937-953. 10 15 ____ 2006. Chaotic population dynamics and the evolution of ageing: proposing a demographic theory of senescence. Evol. Ecol. Res., 8:561-574 MITTELDORF , J. AND J.W. PEPPER. 2007. How can evolutionary theory accommodate recent empirical results on organismal senescence? Theory in Biosciences (accepted; forthcoming). MITTELDORF, J AND WILSON, D.S. 2000. Population viscosity and the evolution of altruism. J. Theor. Biol. 204:481-496 PELS, B., DEROOS, A.M. AND SABELIS, M.W. 2002. Evolutionary dynamics of prey exploitation in a metapopulation of predators. Am. Nat., 159: 172–189. PEPPER, J.W. 2000. Relatedness in group-structured models of social evolution. Journal of Theoretical 20 Biology 206(3):355-368. PEPPER, J.W. AND SMUTS, B.B. 2000. The evolution of cooperation in an ecological context: an agentbased model. In: Dynamics in Human and Primate Societies: Agent-Based Modeling of Social and Spatial Processes. T. A. Kohler and G. J. Gumerman, eds. Pp. 45-76. Oxford University Press: Oxford. 25 ____ 2001. Agent-based modeling of multilevel selection: the evolution of feeding restraint as a case study. Pp 57-68 in: W. C. Pitt, ed. Proceedings of the 4th Annual Swarm User Group Conference. Natural Resources and Environmental Issues, Volume XIII, S. J. and Jessie E. Quinney Natural Resources Research Library, Logan, UT. 30 ____ 2002. A mechanism for the evolution of altruism among non-kin: positive assortment through environmental feedback. Am. Nat., 160:205–213. PRICE, G.R. 1970. Selection and covariance. Nature 227:520-521. ____ 1972. Extension of covariance selection mathematics. Annals of human genetics 35(4): 485-490. 35 PROMISLOW, D.E.L. 1991. Senescence in natural populations of mammals: a comparative study Evolution 45:1869-1887 QUELLER, D. C. 1992. Quantitative genetics, inclusive fitness, and group selection. Am Nat 139(3): 540-558. REZNICK, D.N., BRYANT, M.J. ROFF, D., GHALAMBOR, C.K. & GHALAMBOR, D.E. 2004. Effect of extrinsic mortality on the evolution of senescence in guppies. Nature, 431:1095-1099. 40 RICKLEFS, R.E. 1998. Evolutionary theories of aging: Confirmation of a fundamental prediction, with implications for the genetic basis and evolution of life span. Am Nat. 152:24-44. 34 STEINBEHRENS, B. A. AND R. M. SAPOLSKY (1992). Stress, glucocorticoids, and aging. Aging Clin Expt Res 4(3):197-210. 5 SKULACHEV, V.P. 1997. Ageing is a specific biological function rather than the result of a disorder in complex living systems: biochemical evidence in support of Weismann’s hypothesis. Biochemistry (Moscow), 62:1191–1195. SOBER, E. & WILSON D.S. 1998. Unto Others. Cambridge, Harvard University Press. SOCOLAR, J, RICHARDS, S, AND WILSON, W. G. 2001. Evolution in a spatially structured population subject to rare epidemics. Phys Rev E 63:041908.1-8 10 STEARNS, S. 2000. Life history evolution: successes, limitations, and prospects. Naturwissenschaften 87:476–486 TRAVIS, J.M.J. 2004. The evolution of programmed death in a spatially structured population. J. of Gerontol. Biol. Sci. 59A:301-305 TRAVIS, J.M.J AND DYTHAM, C. 1999. The evolution of dispersal in a metapopulation: a spatially explicit, individual-based model. Proc. R. Soc. Lond. 265:17-23. 15 TRAVIS, J.M.J DYTHAM, C, AND MITTELDORF, J. 2006. Decreased longevity can evolve in response to lower rates of mortality. Submitted to Biology Letters. VAN VALEN, L. 1973. A new evolutionary law. Evolutionary Theory 1:1-30. WADE, M. J. 1980. Kin selection - its components. Science 210(4470): 665-667. 20 WEISMANN A. 1889. Essays upon heredity and kindred biological problems. Clarendon Press, Oxford. WILSON D. S. 1980. The natural selection of populations and communities. Benjamin Cummings. Menlo Park, CA. WILLIAMS, G. C. 1966. Adaptation and Natural Selection: A Critique of some Current Evolutionary Thought. Princeton University Press. 25 ____ 1975. Sex and Evolution. Princeton, NJ: Princeton University Press. WRIGHT, B. 2004. Stress-directed adaptive mutations and evolution. Molec. Microbiol. 52(3):643-650 35