|
|
||||||||
a Max Planck Institute for Demographic Research, Rostock, Germany
b Research Center for Medical Genetics, Russian Academy of Medical Sciences, Moscow
c Department of Cell Biology, University of Calabria, Rende, Italy
d N. N. Petrov Institute of Oncology, St. Petersburg, Russia
e Ulyanovsk State University, Russia
f Department of Molecular Pathology, University of Bologna, Italy
g Center for Demographic Studies, Duke University, Durham, North Carolina
Anatoli I. Yashin, Max Planck Institute for Demographic Research, Doberaner Strasse 114, 18057 Rostock, Germany E-mail: yashin{at}demogr.mpg.de.
Decision Editor: John Faulkner, PhD
| Abstract |
|---|
|
|
|---|
RECENT demographic studies of aging and survival reveal three important features characterizing mortality rates and survival dynamics in populations of developed parts of the world. The first feature is the deceleration of the age-specific mortality rate at old ages (1). The second feature deals with the age pattern of mortality and survival improvement (2) (3). The third characterizes changes in such a pattern during the 20th century. Fig. 1 shows two main patterns of change in survival of Swedish females in the 20th century.
|
An explanation of the mortality-deceleration effect was given by Vaupel and colleagues (1) in terms of a frailty model. Long debates about the compression of morbidity (2) (5) (6) resulted in an admission of the fact that both the effect of rectangularization of the survival curve and the lengthening of the tail of survival distribution do take place in modern survival patterns. In this paper we show that the rectangularization trend in survival improvement that dominated in the first part of the 20th century was replaced by a near-parallel shift of the survival curve to the right in the second part. We also show that the deceleration of old-age mortality, the rectangularization trend, and observed changes in survival patterns with years can be explained by dynamic properties of mortality and survival in a heterogeneous population. The goal of this paper is not to fit model to the data, but to show that all three features can be explained by using the model of population heterogeneity in which the logarithms of mortality rates in subpopulations cross over. We also provide evidence that such populations exist in nature.
For simplicity, we consider the mixture of two such subpopulations, which we call "labile" and "stable." We also assume that survival in these two groups has different responses to environmental changes. An improvement in living standards increases the proportion of initially frail (labile) individuals in such a mixture in the old part of the population. This is an interesting property, because it shows that today's centenarians might not be among the longest-lived individuals if they were born two centuries ago. We discuss the plausibility of the latter model, in the light of recent findings about the aging process in humans and laboratory animals and of the wide data on genetics and physiology of extremely old individuals (centenarians). We discuss empirical evidence of the presence of subpopulations of individuals with different slopes of the mortality curve, and we describe possible mechanisms by which an initially frail organism may become robust (or vise versa) later in life. These mechanisms include the following: (i) the antagonistic gene action; (ii) slower aging, which may be associated with not optimal health parameters in some ages; and (iii) the antagonistic change in stress resistance during an individual life.
Finally, we develop a mathematical model of mortality and aging, which illustrates one possible mechanism of antagonistic change in relative risk in the course of an individual's life. This model explains the possibility of a centenarian's origin from an initially frail part of a cohort. It describes survival of two kinds of individuals, which originally have different sensitivities to changes in environmental conditions and different abilities to adapt to these changes. The model establishes interrelations between these characteristics and survival. We show that the trends in survival patterns experienced by a mixture of such individuals resemble those observed in humans during the previous century.
| The Model |
|---|
|
|
|---|
Antagonistically Changing Frailty in the Process of an Individual's Life
Let us assume that the human population is a mixture of two subpopulations of individuals, both with Gompertz mortality rates, but with different slopes of a logarithm of mortality. At the beginning of the century the mortality rate for one subpopulation with lower slope is higher than that of the other. The intersection of such curves happens at an age that is beyond the observed range of the human life span (say, at 130 years, or more). The improvements in the standards of living produce a significant reduction of the mortality rate for the first (labile) subpopulation, and less significant changes in the survival chances in the second (stable) subpopulation (see
Fig. 2).
|
85 years. After this time the two curves evolve together (shift to the right) with almost the same speed (
Fig. 3).
|
Fig. 4 shows how survival function corresponding to Swedish females in 1861 can be approximated as a mixture of two survival functions. The lower survival function corresponds to labile individuals; the higher survival function represents stable individuals. Fig. 5 shows how the survival function for Swedish females in 1995 can be approximated as a mixture of survival functions for labile and stable individuals with the same properties. In this figure the survival function for labile individuals intersects that of stable individuals around age 70 years. This intersection happens as a result of faster progress in the mortality reduction in the labile subpopulation.
|
|
Model of Difference in Stress Response That Produces an Intersection of Mortality Rates
We suggest a model of physiological aging and survival in which adaptation to the stresses of life has a metabolic cost. The model shows that today's centenarians may originate from an initially frail part of a cohort. The model describes survival of two kinds of individuals, the labile and stable. These individuals have different sensitivities to changes in environmental conditions, and so different amplitudes of the response to stress. The model establishes interrelations between amplitude and survival. The low amplitude of the response to stress corresponds to stable individuals; the higher amplitude corresponds to labile ones. The mathematical development is reported in the Appendix; here a qualitative description is given.
For simplicity, we assume that at any age an organism may be in one of two possible states. One is the state of normal functioning. The other is a state of stress, disease, or other tension (arousal), induced by some external or internal conditions, in which the systems of an organism are functioning with an "overload"
(7). A long stay in this state must have a cost for an organism: when it accumulates, the survival chances of an organism decline. The response of an organism to stress is regulated by the parameter
, called the rate of recovery, or the adaptation rate (see the Appendix). The lability property is associated with small values of
; stability is associated with high ones. Because those who cannot adapt properly die first, the survival function of individuals with a small adaptation rate (i.e., labile individuals) is initially lower than that of individuals with high values of this rate (i.e., stable individuals). This means that stable individuals have a survival advantage earlier in life. The model establishes interrelations between the rate of adaptation and the accumulation of damage: a high rate yields a metabolic cost included in the definition of allostatic load; see 3 of the Appendix. Then, at older ages, individuals with smaller values of the adaptation rate may have higher survival chances than those with higher values of
. That is why the survival curves for these two groups of individuals may intersect. Another important parameter of the model
characterizes permanently acting environmental disturbances. The higher the
, the higher the level of disturbances. The life-span distribution in this model is obtained by simulation of individual life spans for a large sample of individuals (see the Appendix).
Fig. 6 shows survival patterns, generated by the model, in hypothetical populations of labile (small
) and stable (large
) individuals, for two values of
: the two values of
represent the transition from low standards of living (high value of
) to improved ones (lower value of
).
|
, respective survival functions intersect. The age at intersection moves to the left when the level of disturbances declines. Perhaps the most striking result shown by the simulation is that the variation of
(an increase in living standards) has a small impact on the stable population, whereas it substantially increases the survival of the labile one, which exceeds that of the stable population at older ages. Fig. 7 shows the survival patterns of a mixture of labile and stable individuals (with survival functions shown in Fig. 6) before and after an improvement in living standards. One can see from this figure that the resulting survival function shows both features. It became more rectangular, and it has a longer tail than the initial one.
|
| The Genes: Experimental Data |
|---|
|
|
|---|
values) with respect to those of younger individuals (low
), the cohort of centenarians should be formed by both labile and stable individuals (
Fig. 6). The model discussed above can be checked by comparing the levels of genetic homogeneity between centenarians and younger individuals (2060 years old) extracted from the same population. Let us consider a stress responder gene, with A and a alleles conferring lability (small values of
) and stability (high values of
), respectively. Let us assume that the frequencies of A and a at the birth of the cohort are pA = pa, and that the population is in HardyWeinberg equilibrium. When survival selection operates under high
(thin lines in
Fig. 4), the allelic pool will tend toward an increase of a, with a consequent decrease of heterozygosity. However, as
decreases (thick lines in
Fig. 6), A will tend to remain in the gene pool, and the level of heterozygosity will increase in the cohort of centenarians. Therefore centenarians are expected to show increasing heterozygosity with respect to youths for stress-responder genes affecting labilitystability. We checked our hypothesis by estimating heterozygosity at 12 autosomal loci in centenarians and younger individuals, after verification of HardyWeinberg equilibrium. The results are shown in Table 1 .
|
| Discussion |
|---|
|
|
|---|
Paradoxes in Centenarians
The progress in health and living standards experienced by the human population today tends to increase the proportion of originally frail individuals in successive generations
(11)
(12). The proportion of centenarians in human populations of developed countries also increases faster than before
(4). It is clear that both effects may mirror an influence of the increasing standards of living on the mortality curves of homogeneous subcohorts, comprising heterogeneous groups of individuals. This, however, does not explain the observed trends in survival patterns. The explanation would be easier to obtain if we assume that a substantial part of the group of centenarians originates from an initially vulnerable, frail part of a generation, whose deaths at adult age were prevented by the improvements in health and living standards gained by industrial progress. Because industrial progress is likely to increase survival chances for all individuals in the population, one should find a good reason why some originally frail individuals get survival advantages at old age. Several examples discussed below provide compelling evidence about the possibility of the original frailty of today's centenarians.
Recent genetic studies in centenarians pointed out some unexpected findings, as the presence of genetic risk factors in their gene pool (13). For example, several alleles, known to be associated with increased cardiovascular risk in middle age, are present in centenarians at the same frequency as in younger individuals (14) (15) (16) (17). In some cases, such risk factors occur even more frequently in centenarians than in younger individuals (18) (19) (20) (21) (22). Such counterintuitive accumulation of originally harmful alleles may be easier to understand if the biological and physiological role of respective alleles in survival is taken into account.
For example, let us consider the guanine insertion/deletion polymorphism 4G/5G in the promoter of the plasminogen activator inhibitor 1 (PAI-1) gene, a predictor of the risk of atherothrombotic disease. The 4G4G genotype is associated with a high plasma level of plasminogen activator inhibitor and therefore with an increased risk of atherothrombosis and myocardial infarction in adulthood (19). Unexpectedly, the frequency of the 4G4G genotype is higher in centenarians than in younger individuals (20). Our model can explain this apparent paradox on the basis of the physiological role of the PAI-1 hemostatic protein. The improvement in medical care and living conditions in early and in middle ages (e.g., by saving lives and preventing deaths of those who would otherwise die) promotes an increase in the proportion of individuals carrying the harmful genotype. However, this genotype, which is a risk genotype at middle age, may turn out to be advantageous at older ages, when the rate of metabolic processes decelerates, and many processes including a recovery from injury go slower. In such ages, a higher rate of blood coagulation may be beneficial for faster recovery (e.g., for stopping bleeding). Thus individuals with potentially harmful alleles, who yet survived under the pressure of selection and reached old age, may get an advantage from the same alleles later in life.
Centenarians may also be those individuals whose development (and aging) processes go slower than those in other individuals of the same generation. Such individuals may have a higher risk of death earlier in life because their health parameters at young and adult ages may substantially deviate from those "optimized" by the evolution. This is confirmed by the fact that the risk of death at a given age regarded as a function of physiological parameters can often be approximated by a U-shaped curve (23). For example, being both underweight and overweight is associated with an increased risk of mortality. A person who grows old slower may have a body mass index as well as other physiological parameters (such as metabolic rate) that are not optimal for survival at early and adult life. But if such a person survives this period, he or she may have a survival advantage at older ages because of a slower aging phenotype (see Fig. 8). Progress in medicine and health care together with improvements in living standards increases survival of such individuals earlier in life. This yields an increase in the proportion of such individuals in old ages, where they get survival advantages, and it contributes to the increase in the proportion of centenarians in a population.
|
The same considerations can be applied to model organisms. Stress experiments with nematode worms Caenorhabditis elegans reveal an intersection of survival curve in the population exposed to 4 hours of heat shock with that in the control group (30). The survival curves corresponding to a population of some long-lived mutants of C. elegans intersect those of wild-type worms under normal conditions (31).
Experiments with rodents give one more confirmation of the existence of the labile phenotype. In aging studies on rodents, rats and mice exposed to high doses (0.5 mg/rat) of pineal gland preparation epithalamin (known as geroprotector) show lower survival earlier in life and better survival later in life than those in the control group ( Fig. 9). These rodents also manifested both increased maximal life span and postponed aging (32) (33) (34). That is, treated rodents lived longer and aged slower (by a number of physiological parameters). At the same time, they are frailer at their young and adult ages, and they are more robust at old ages than rodents in the control group. Furthermore, in populations with lower survival earlier in life and higher survival later in life, tumorigenesis is going slower (32). This observation is expected in the population of labile individuals, in which the accumulation of damage associated with the excess of metabolic activity is slower than in the stable individuals. A similar antagonistic change of survival in adult and in old ages was observed in some experiments on rodents exposed to caloric restriction, a treatment known to increase longevity and postpone aging (35).
|
Evolutionary Adaptation May Produce an Intersection of Mortality Curves
The intersection of survival curves may have evolutionary reasons. The evolutionary theories suggest that genes that would enhance fitness earlier in life may be selected for, even if they produce detrimental effects later in life
(37). Because the notion of fitness involves fertility as well as mortality, one may assume that mortality rates for genotypes may have a spectrum of possible shapes (the average fertility rates for the populations of these genotypes must also be different). These theories link the life span with the quality of maintenance and repair systems in the cells of an organism. However, they do not provide an explicit description of possible mechanisms, which relate characteristics of physiological and biological aging with demographic survival curves. The idea of disposable soma suggested by Kirkwood
(38) looks promising; however, it requires additional specification of resource-allocation strategies that may depend on environmental conditions. In addition, none of the theories explains why some individuals in a population are more sensitive to industrial progress than the others.
Conclusions
Mortality decline in the human population forms a survival pattern, which shows signs of both the rectangularization of the survival curve and the lengthening of the tails of survival distributions
(5)
(39)
(40). In this paper we propose a mixed stochasticdeterministic mathematical model of survival, which incorporates parameters directly related to qualitative features such as adaptive capacity and sensitivity to environmental changes of an individual. The qualitative analysis and the simulation experiments show that the observed pattern of mortality as well as trends in human survival can be explained in terms of a mixture of two subpopulations. An important feature of these subpopulations is that their mortality rates intersect, thus giving a solid interpretation to the intersection of mortality curves often observed in the studies of aging and survival, including human and nonhuman subjects.
| Acknowledgments |
|---|
Received January 29, 2001
Accepted May 11, 2001
| Appendix ENDIX |
|---|
|
|
|---|
, µ, and A. The first, X = (Xt)t
0, characterizes the rate of metabolic processes. It may be regarded as a difference between metabolic rates in the state of normal functioning and in the arousal state. The dynamics of this process reflects changes in metabolic characteristics during the adaptation to the new state of an organism induced by external or internal conditions. For simplicity, we assume that there are only two possible states"0" and "1." State 0 is characterized by the absence of stress, disease, or other conditions that require substantial adaptational efforts. State 1 is an arousal state. It is characterized by elevated metabolic activity (as in the case of stress, disease, etc.). The process N = (Nt)t
0 describes random changes between these states. The process
t characterizes the intensity of transitions (incidence of disease, stress hits) resulting in an arousal state
![]() | (1) |
is a nonnegative parameter. This process depends on At, the fifth process contributing to aging, which we associate with the price for staying in the arousal state. The process µt describes the rate of recovery from the arousal state. We assume that
![]() | (2) |
Here m and b are nonnegative parameters. The process A = (At)t
0 describes the accumulation of damage, that is, a cost paid by an organism associated with its stay in the arousal state (e.g., an allostatic load;
(7),
(41),
(42)). It seems reasonable to assume that A = (At)t
0 depends on the history of the process X up to time t, for example,
![]() | (3) |
are nonnegative parameters. We also assume that the process X is related to process N by the equation
![]() | (4) |
Here W = (Wt)t
0 is the standard Wiener process, and parameters
and
> 0 are fixed for all genetically identical individuals. This equation contains deterministic and stochastic components. The deterministic part represents the homeostatic feedback mechanism by which process X adapts to the level of N. The stochastic component characterizes permanently acting disturbances, for example, oxidative stress or immunological stress, which disturbs homeostasis. It is important to note that parameter
influences both the X process and the A process. In the first case it characterizes the rate of adaptation of an organism to changes in the state. The higher the value of
, the faster the organism adapts. The variance of fluctuations of the process Xt is
2/2
; that is, it is inversely proportional to parameter
. Thus the weaker the adaptive capacity of an organism, the higher the fluctuations of X. The high rate of adaptation is, however, not for free. It has a metabolic cost included in the definition of allostatic load
(7)
(41). Note that according to Fries
(3), even for the same trajectories of X the allostatic load will be accumulated faster for an individual with larger values of
. Thus any jump of Nt from 0 to 1 induces an increase of an average level of X. An increase in Xt causes faster accumulation of allostatic load At and finally results in a reduction of a stress resistance µt and in an increase of intensity of diseases
t.
Three possible causes of death were considered in this model. The first is associated with a decrease of vital characteristics below the admissible level (in the example given below, Xt < -1). The second arises when the recovery rate falls below its admissible level (in the example given below, µt < 0.01). The third arises when X exceeds the admissible level (in the example given below, Xt > 2).
Simulation of the Modern Pattern of Survival
We simulated life-span data by using a computer experiment with our model and calculated survival functions for different sets of parameter values. The control set (F1) of the parameters is
= 370,
=550,
1 = 6.4,
2 = 6.0 a = 0.00066,
= 0.2, b = 0.5, and m = 10. The statistics is based on 2000 independent random realizations consisting of two cohorts: 25% (500 individuals) with
= 370 and 75% with
= 550 (1500 individuals); T* = 40 years. The sample size of simulated data is 2000 individuals.
The parameter
in this model is supposed to be genetically determined and fixed for all individuals in a given birth cohort. As it defines the rate of adaptation and is inversely proportional to the variance of Xt, an increase in
leads to an increase in survival probability for the younger ages, where the allostatic load is small and most of the deaths occur because of fluctuations of Xt. Because the allostatic load increases faster with age for such individuals, the survival curve for them goes steeper to zero. Thus we can simulate the phenomenon of the intersection of the survival curves by assigning different values of parameter
to respective populations of individuals. Note that the old-age survival probability is significantly higher for the organisms with a lower level of feedback rigidity (with smaller
). The main cause of the mortality among individuals with strong feedback (large values of
) is the fast decrease in the recovery rate caused by accumulation of damage (allostatic load). The individuals with low feedback (weak homeostatic mechanisms, small values of
) die mostly because of crossing the admissible boundary by the process Xt.
Thus individuals with the large values of
usually do not die from the fluctuations of Xt, because they are small. However, the accumulation of allostatic load for them is going faster. Individuals with small
have higher levels of fluctuations in Xt (are more labile). Thus they die more often at the beginning of life. However, the allostatic load in these individuals increases at a slower rate, which gives them better survival chances later in life compared with more stable individuals.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. M. Kulminski, S. V. Ukraintseva, I. V. Culminskaya, K. G. Arbeev, K. C. Land, L. Akushevich, and A. I. Yashin Cumulative Deficits and Physiological Indices as Predictors of Mortality and Long Life J. Gerontol. A Biol. Sci. Med. Sci., October 1, 2008; 63(10): 1053 - 1059. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Rockwood and A. Mitnitski Frailty in Relation to the Accumulation of Deficits J. Gerontol. A Biol. Sci. Med. Sci., July 1, 2007; 62(7): 722 - 727. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. De Benedictis and C. Franceschi The Unusual Genetics of Human Longevity Sci. Aging Knowl. Environ., June 28, 2006; 2006(10): pe20 - pe20. [Abstract] [Full Text] |
||||
![]() |
J. E. Morley Editorial: Hot Topics in Geriatrics J. Gerontol. A Biol. Sci. Med. Sci., January 1, 2003; 58(1): M30 - 36. [Full Text] [PDF] |
||||
![]() |
J. E. Morley Editorial: Citations, Impact Factor, and the Journal J. Gerontol. A Biol. Sci. Med. Sci., December 1, 2002; 57(12): M765 - 769. [Full Text] [PDF] |
||||
![]() |
J. E. Morley, H. M. Perry III, and D. K. Miller Editorial: Something About Frailty J. Gerontol. A Biol. Sci. Med. Sci., November 1, 2002; 57(11): M698 - 704. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|
| All GSA journals | The Gerontologist |
| Journals of Gerontology Series B: Psychological Sciences and Social Sciences | |