Karl Inne Ugland,
University of Oslo,
Marine Zoology,
Pb 1066, 0316
Oslo, Norway. Introduction
Classical and molecular genetics are concerned with the nature and transmission of genetic information, and how this information is translated into phenotypes. Population genetics theory is most successful when dealing with simply inherited traits  traits whose transmission follows simple Mendelian rules. Yet many of the most interesting and important traits are not so simply inherited: they depend on several genes, which often interact in complex ways with one another and with the environment.
Population genetics includes a large body of mathematical theory; one of the most richest and most successful bodies of mathematics in biology. The useful application of this theory has been greatly enhanced in recent years by new molecular techniques.
Usually, the first step is to study the frequency of different genotypes or phenotypes in a sample of the population. The concept locus is used to designate a chromosomal location, and the concept allele designates an alternative form of the gene occupying the considered position. Usually we are more interested in the frequencies of the different alleles than in the frequencies of the different genotypes. This is because allele frequency is a more economical way to characterise the population. The number of possible genotypes is enormous. Consider for example 100 loci, each with 4 segregating alleles. With 4 alleles A, B, C, D the total number of possible genotypes at each locus is 10: (1) four homozygotes: AA, BB, CC and DD plus (2) six heterozygotes: AB, AC, AD, BC, BD and CD. For all 100 loci there are thus 10^{100} possible genotypes, i.e. more than the number of atoms in our universe! It is, therefore, much better to stick to allelefrequencies.
The characterisation of a population in terms of allele frequencies rather than genotype frequencies has another advantage. In a Mendelian population, the genotypes are scrambled every generation by segregation and recombination. New combinations are put together only to be taken apart in later generations. If we are to study a population over a time period of more than a few generations, then the stable entity is the gene.
In nature it is the organism that survives and reproduces, so natural selection acts on the organisms and it is the fitness of the organism that determines the likelihood of survival and reproduction. But what an individual actually transmits to future generations is a random sample of its genes, and those genes that increase fitness will be more represented in future generations. Thus, there is selection of genes determining the properties like vigour and fertility.
Any description of nature, whether verbal or mathematical, will only be a caricature and therefore necessarily incomplete. In the previous century, scientists realised that we cannot ask whether a mathematical model is true, we can only ask if it gives a good or bad description of our data, and so may be used for prediction of future observations and experiments. Some models are only rough caricatures, but the advantage with this class of models is that they are easy to understand. An example is the successful theory of thermodynamics where the gas molecules are regarded as elastic spheres.
The gene pool
One of the most useful conventions in population genetics is the model of a gene pool. We assume that each parent contributes equally to a large pool of gametes. Each offspring is regarded as a random sample of one egg and one sperm from this pool. Each gene is chosen randomly, as if we were drawing different coloured beans from a bag. This oversimplified model building is called "beanbag genetics" and has a surprisingly strong predictive power. In 1964, the great British geneticist, J.B.S. Haldane, wrote an amusing and spirited article "A Defence of Beanbag Genetics" (Persp. Biol. & Med., 7: 343359).
Clearly the beanbag model is a crude representation of nature. A real population has individuals of all ages, some dying, some choosing mates, some giving birth, etc, etc. Keeping track of all this information is not only impractical, but often is neither interesting nor important. Most questions of genetic or evolutionary interest do not require minute details. The gene pool model will give us the same kind of insight that simple models in the physical sciences do, and often with predictions that are sufficiently accurate for most uses. Further it is important to be aware of that departure from the simple models may usually be treated with appropriate modifications.
Suppose now that a population consist of individuals with an extremely long life span. In the tradition of "beanbag genetics" the population is homozygous at the locus determining life span; i.e. all genotypes are AA and code for longlife. However, in any population there will be some random death rate basically determined from outer sources like biological (virus, sickness, predation or accidents due, for example, to nature catastrophes). Thus, rather living forever, the longlived strategy has to tolerate the operation of a random death rate; that is, the survival rate, S, is close to 1. This strategy is called evolutionary stable if the population rejects all kinds of invading alleles that reduce the life span. In order to investigate this stability, we assume that a small fraction ε of the gene pool consists of an allele 'a' such that the genotypes Aa induces a shorter life span via changes in the birth rates and survival rates. We must find the requirement that the new allele will increase in the gene pool.
A HardyWeinberg type model for analysing the evolutionary stability of infinite longlife
Since we are only interested in the general principles that govern the stability of the longlife strategy we simplify the population dynamics of the mutant allele as follows: On average they live t years and reproduces b_{Aa} offspring of which l_{Aa} survives the next t years. The probability of surviving t years for the longlived individuals is L = S^{t} = (S multiplied by itself t times). Finally, we ignore the homozygous mutant allele (of genotype aa), since their abundance is of the order ε∗ ε in the gene pool.
Thus the initial genotype frequencies are (1ε) of AA and ε of Aa where the abundant first group of individuals (AA) have a very long life span and the second potential invader group (Aa) only live t years. In the gene pool the frequency of A is (1ε)*1 + ε*(1/2) = 1ε/2 while the frequency of a is .
We shall now try to answer the question: Under what circumstances is the fraction of the allele a expected to increase in the gene pool?
In order to solve this problem we need another assumption: The population size is so large that the mating may be considered as a random mixture of the genes. We must consider all kind of matings between the genotypes:
(1) Consider first a longlived female (AA) that mates with a longlived male (AA). During t years she produces b_{AA} offspring (of genoype AA) that have a chance of l_{AA} to survive the next t years. In the population the probability of such matings is (1ε)*(1ε).
(2) Consider next a longlived female (AA) that mates with a mutant shortlived male (Aa). During t years she produces b_{AA} offspring of which 50% are longlived (AA) with a chance of l_{AA} to survive the next t years, and 50% are shortlived (Aa) with a chance of l_{Aa} to survive the next t years . In the population the probability of such matings is (1ε)*ε.
(3) Then consider a mutant shortlived female (Aa) that mates with a longlived male (AA). During t years she produces b_{Aa} offspring of which 50% are longlived (AA) with a chance of l_{AA} to survive the next t years, and 50% are shortlived (Aa) with a chance of l_{Aa} to survive the next t years . In the population, the probability of such matings is ε*(1ε).
(4) Finally consider a mutant shortlived female (Aa) that mates with a mutant shortlived male (Aa). During t years she produces b_{Aa} offspring of which 25% are longlived (AA) with a chance of l_{AA} to survive the next t years, and 50% are shortlived (Aa) with a chance of l_{Aa} to survive the next t years. In the population the probability of such matings is ??ε*ε.
After t years the number of longlived and shortlived recruits are respectively
W_{AA} = (1ε)^{2} b_{AA} l_{AA} + (1ε)*ε (1/2) b_{AA} l_{AA} + ε*(1ε)(1/2) b_{Aa} l_{AA} + ??ε*ε(1/4) b_{Aa} l_{AA}
W_{Aa} = (1ε)*ε (1/2) b_{AA} l_{Aa} + ε*(1ε)(1/2) b_{Aa} l_{Aa} + ??ε*ε(1/2) b_{Aa} l_{Aa}
Now, the number of adult longlived individuals surviving a period of t years will be (1ε)N^{ }*S where N is the total population size. In order to keep the population at a stable size of N individuals, we have to introduce a common death rate m for the recruits of all genotypes:
(1m)*N*W_{AA} + (1m)*N*W_{Aa} + (1ε)N^{ }*L = N
After t years the new fraction of the shortlived allele a in the gene pool is
Thus we see that the fraction of the allele a will increase in the gene pool if
,
i.e.
which is equivalent to
Reordering the terms, finally gives the requirement that the mutant allele coding for shortlived may invade the gene pool:
Thus, if the mutant allele a codes for a birth rate b_{Aa} and survival rate l_{Aa} such that this inequality is satisfied it will invade the gene pool. But what does this inequality essentially tell us? In order to see the essential requirement for invasion we may use standard techniques in numerical analysis. First, since ε is a small quantity in the order of 1/100, we may use the standard approximations in the final expressions: _{}
so the inequality is practically the same as
Then we take a closer look at the extra nonspecific mortality rate m that was introduced in order to keep the population at its equilibrium value N. In order to understand the magnitude of this quantity we consider the population without the mutant allele a (that is all genotypes are AA), so the equilibrium equation reduces to N*W_{AA} + N^{ }*L = N which is seen to imply
Equilibrium condition:
W_{AA} + L = 1
The equilibrium equation simply states that in order to stabilise the population, the recruitment to the stock must be balanced with the death rates. Now according to classical Darwinian thinking, successful new mutants only accomplish small changes. In our case, this means that the new allele a will only slightly modify the recruitment, so the final adjustment of the new shortlived recruits will be small. Hence
_{
}so we finally obtain the following essential criteria for successful invasion:
Criteria for successful invasion of the gene pool:
In words this criteria simply says that the new allele a will successfully invade the gene pool if the number of surviving mutant offspring to adulthood exceeds 1. Note that offspring of genotype Aa is created in two different ways: Either the mother is of genotype AA and mates with a male of genotype Aa in which case 50% of the offspring will be Aa, or the mother is of genotype Aa in which case also 50% of the offspring will be Aa. This explains why it is the average litter size of the two genotypes that enters the invasion criteria. Now the expected number of mutant offspring that reaches adulthood is found by multiplying with their survival rate l_{Aa} . Since the recruitment of the longlived genotypes exactly matches the adult mortality 1  L, the total recruitment to the adult group is exactly 1 (i.e. the equilibrium condition). Note that this is precisely where the longlife enters the equations: the total recruitment is the recruitment of offspring that manage to live up to adulthood PLUS the number of adults surviving the period (as mentioned above, only nongenetic accidents cause a small death rate among the longlived individuals). So the criteria for successful invasion of the gene pool may be restated in words as follows:
Criteria for successful invasion of the gene pool:
The recruitment of the shortlived individuals must exceed 1
Note that this is good classical Darwinism: The genotype with the highest recruitment will at the end fill up the whole gene pool.
It is now easy to see why Jim Baen´s hypothesis is supported by this genetic model. The whole issue about why we must die may be formulated the question of whether a longlived strategy is evolutionary stable. A simple, first approximation, genetic model says that the longlived strategy will only be evolutionary stable if it is difficult for mutant alternatives to establish a recruitment larger than 1. But as Jim Baen has pointed on classical biological knowledge: there are all sorts of difficulties with inbreeding. Therefore almost any alternative strategy that reduces energy devoted to longlife and put it into clutch size or vitality of offspring will have a better total recruitment than the longlived. Therefore, the strategy of eternal life is evolutionary unstable. For human beings this has the consequence that 125 years is the maximum age.
Cohort analysis of the evolutionary stability of infinite longlife
In this section we derive the same result in an alternative model that takes all the yearclasses (cohorts) into account.
Suppose there are M adult longlived individuals consisting of 1000 year classes, and let the clutch size per individual be b_{A} , the survival up to maturity is l_{m} and the survival rate of adults is s_{A} . At equilibrium Mb_{A} offspring are produced each year by 1000 year classes:
which is seen to reduce to
Equilibrium condition
i.e. the population will be in equilibrium if the recruitment to the adult population (i.e. birth rate multiplied by immature survival) balances the adult death rate.
Consider now a potential invading allele a that codes for living only 100 years as adults, i.e. 10% of the prevailing strategy of 1000 years. Initially there are individuals of these new genotypes having a birth rate of b_{a} , a survivorship of the immature period equal l_{a} and an adult survivorship of s_{a} . With these parameters the iteration equation for the adult mutants are
Substituting a general solution of the form gives the characteristic equation for the growth of the mutant
In order to see where the root lies, define the two functions
and notice that
,
but since tm > tm1 , f(r) will eventually become greater than g(r) . Thus somewhere they intercept, and the crucial question is whether they intercept to the left of r = 1 , in which case the number of invaders will decline, or to the right of r = 1 , in which case the number of short lived individuals will increase. Since
it follows that the condition for interception to the right of r = 1 , i.e. f(1) < g(1) ,is
which finally gives
Criteria for successful invasion of the gene pool:
In words: The recruitment to the mutant shortlived individuals (=average birth rate of long and short lived genotypes multiplied with the survival probability up to maturity) must be larger than the adult death rate of the shortlived individuals. Note that this condition is equivalent to that found in our previous HardyWeinberg type model, so the same discussion applies.
Let us finally see that if the following approximation of the population dynamics of the mutant
is substituted into the formula of the HardyWeinberg model, we arrive at the same criteria for εM mutant to invade a homozygous population of (1ε)M longlived adults:
where the second term in the denominator is the number of longlived adult that survived the total life span of the short lived (i.e. 100 years) and the third term is the number of adult longlived that recruits to the adult population of long lived (since these are prevailing we ignore the contribution of longlived from the shortlived). Using the equilibrium criterion,
,
we obtain
which is equivalent to
Thus, we arrive at precisely the same inequality, demonstrating an inner consistency in our models.
A competition model with tradeoff between reproduction and survival.
While the two previous models started with the assumption that longlife was an evolutionary stable strategy, we now develop a model where the different life history strategies are in direct competition with each other. However, if a new mutant code for better survival we shall now explicitly regard this as a strategy to allocate more energy to survival at the cost of energy devoted to reproduction. The simplest equation for such a tradeoff model is
A first approximation to the population growth will then be
which has the solution
From this we see that natural selection will favour a strategy which maximises the function
Putting the derivative equal to zero,
,
gives the optimum survival rate as
_{
}implying an evolutionary stable expected life span of
So also in this model will the longlived be selected against because they do not provide the optimal combination of birth rate and survival rate. Another consequence of tremendous importance is that given a strategy near the optimal, there will be no selection against deleterious genes at ages over T^{*}. This allows death genes, racingcar analogy construction of organs, etc. Now human beings have evolved from species where these mechanisms already have been operating, so this is an evolutionary blind route with respect to longlife unless it is possible to enter our genes directly and reconstruct them. We are an animal that has evolved over many millions of years to be short lived.
