| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|
| ||||||||||||||||||||||||||||||||
1 Institute of Control Sciences, Moscow, Russia.
2 Max Planck Institute for Demographic Research, Rostock, Germany.
3 Wayne State University, Detroit, Michigan.
4 Department of Entomology, University of California at Davis.
5 Duke University Center for Demographic Studies, Durham, North Carolina.
Address correspondence to A. I. Yashin, Duke University Center for Demographic Studies, 2117 Campus Drive, Box 90408, Durham, NC 27708. E-mail: yashin{at}cds.duke.edu
| Abstract |
|---|
|
|
|---|
Despite this general agreement on the descriptive aspects of this relationship, the prior work did not shed much light on the mechanisms underlying these widespread alterations in fecundity. In large part, this was because a robust analysis needed individual female, as opposed to population-based fecundity data. Although such individual female experiments have been done for some time (1,4,1215), their intensive study began only recently (14,1619).
Pretzlaff and Arking (20) were the first to note an exponential decrease in a mean-population fecundity in Drosophila. A theoretical prediction of optimal fecundity pattern in Drosophila was presented by Stearns and colleagues (21), but unfortunately that study did not report how the prediction was yielded. Muller and colleagues (22) successfully simulated the declining part of the age-related fecundity in individual Mediterranean fruit fly (medfly) females by an exponential function randomly beginning at an age of egg-laying peak, but unfortunately that study did not analyze the initial and mid-life part of the mean-population or individual fecundity pattern. Kindlmann and colleagues (23, p. 837) developed a model in which "inclusion of the senescence function ... resulted in a triangular shaped fecundity function." Nonetheless, these authors do not use any quantitative criteria to estimate the "triangularity" of the pattern. Dixon and Agarwala (24) compared such a form with a mean-population fecundity pattern in ladybird beetles, but they also did not go beyond the triangular appearance of the pattern. Finally, Rauser and colleagues (25) presented evolutionary reasons why mean fecundity values plateau at very advanced ages. None of these studies used a rigorous quantitative analysis.
Two other studies did attempt to quantitatively describe the overall fecundity pattern, but each yielded only rough approximations. Shanley and Kirkwood (26) used a cosine presentation of an overall fecundity pattern in a mouse population, and Cichon (27) calculated schematic curves of reproductive rates in Drosophila, presumably formed by optimal resource partitioning. There was no comparison with experimental data in either case.
These prior studies may generally be summarized as yielding a qualitative but not a quantitative understanding of the relationship between fecundity and longevity in this organism. As a consequence, no rigorous mechanistic understanding was possible.
As described in our preceding articles, we have developed and verified a mathematical model of Drosophila melanogaster female fecundity which includes the three phases of maturation, maturity, and reproductive senescence. Initially, Novoseltsev and colleagues (28) analyzed the theoretical intrinsic mechanisms that relate egg production to longevity in individual fly females. This model assumed that the observed fecundity for any individual female represented the interactions and trade-offs between an age-related reproductive program and an age-related decline in the homeostatic capacity of the reproductive system (presumably due to increasing oxidative vulnerability). We showed that the individual rate of egg-laying followed a genetically prescribed constant maximal value. This means that, until the onset of senescence, the individual fecundity pattern is flat and has no maximum. After that time, it declines exponentially. To describe such a pattern in adult females, we assumed the existence of three age-related stages: maturation, maturity, and senescence. These phases are described by four parameters, namely Xonset, RC, T, and
sen (Figure 1).
|
sen is the rate of senescence. Maturity presumably ends at the onset of senescence Xsen = Xonset + T, after which the senescence lasts S days, S = (LS T Xonset). Here LS means life span. The model was quantitated using a preexisting set of individual female longevity and fecundity data. This was subsequently verified as being applicable to other independently published Drosophila data sets (29). The model was then tested on independent data sets of female Medflies and was found to generally apply, although questions about the nature of the reproductive senescent phase in the medfly were left unanswered (e.g., the expected trade-off between fecundity and longevity was not observed in the medfly populations).
Thus our model generally mirrors the genetic and physiological features of the Diptera species, although the analysis to date suggests that the reproductive senescent phase of the medfly might involve alterations in different parameters in the medfly as compared to Drosophila (30).
Thus, we now extend our prior analysis by focusing on the manner in which aging and senescence develop in Drosophila and medflies. We use a different analytical method to detect the parameters associated with longevity for each species. Having identified the species-specific parameters, we then proceed to analyze the manner by which reproductive senescence develops in each species. We conclude that reproductive senescence in these two species involves different parameters, and this finding may explain the different patterns of reproductive senescence observed in these two laboratory model organisms.
| METHODS |
|---|
|
|
|---|
The Drosophila data used involved three control (i.e., normal-lived) lines in an artificial selection procedure for shortened and prolonged life span (1,4,31). Thus the united population data set consists of Ra, Rb, and Rc flies (149, 172, and 172 flies, relatively). Individual numbers from 1 to 493 were given to these flies.
A newly eclosed male and female fly were placed in vials which were changed daily, and fecundity and survival were tallied daily. Age at death was documented and longevity calculated. In the event of the death of a male in any of the assay vials, the dead fly was removed and replaced by another male of the same age from a cohort of reserve flies (1,4,20). However, the disturbed patterns were taken into account in our study.
Previously, egg-laying scores in medfly were approximated by regular three-stage patterns via the least-mean-square procedure (30). To test if it is possible to use this approach in Drosophila, we applied this fitting procedure to the individual fecundity scores of Drosophila females. This procedure allows us to divide the fecundity score into three parts, the first referred to as maturation, the second maturity, and the third senescence.
The methods used for the detection and analysis of critical points in the onset of senescence are described in conjunction with the results, as this approach is more informative to the reader.
| RESULTS |
|---|
|
|
|---|
After this point, the fecundity rate starts to decrease, presumably at the exponentially decelerated rate.
sen is the rate of senescence, which has a dimension of (1/day) and means that, at senescence, daily egg laying diminishes by exp[
sen] each following day as related to the current day. Mathematically,
sen is the exponent of the exponential function that describes an age-related decrease in fecundity at older ages. Previously, we and others (28,31) used another parameter,
tail = (
sen)1, instead of
sen.
tail is the time constant of the reproductive "tail" in the individual fecundity pattern. The senescence stage lasts for S days and ends when death occurs, at age LS.
Nonetheless, the random character of egg scores calls the existence of discrete times of onset of senescence into question. Indeed, in some patterns the location of this point is obvious, but in others, it is not. To test if there really exists a critical point that manifests the onset of senescence, we apply the following procedure.
The least-mean-square procedure allows dividing maturity stage from senescence stage in a fly female by fitting the individual scores with a three-stage regular pattern. To the time axis we allocate all individual egg scores so as to put together individual change points predicted by least-mean-squares to be the onsets of senescence. An example of how this is done is shown in Figure 2. The egg scores for the jth fly are presented with the presumable onset of senescence placed in a central point. This point is positioned exactly between two time intervals. To the left of the point, the maximal length of the plateau (Tmax) is placed, and to the right, the maximal length of senescence tail in the studied population, Smax.
|
The graphs below schematically show the changes in the number of flies that are to be averaged at each i, where the centered age is i = 0, 1, ... , (Tmax + Smax). The number of flies is referred to as Ni. The graphs Ni(i) start from Ni = 1 at point i = 0 (i.e., from the fly with the maximal length of maturity stage) and at the point i = Tmax arrives to the maximum value, which is equal to the population quantity, N. Then it moves downward to Ni = 1 at point i = Tmax + Smax (the fly with the maximal S length).
Let r(i,j) denote a number of eggs laid by fly j at day i. The difference Di between the scores and the respective plateau levels, averaged in series at all ages, yields
|
|
The pattern Di for the Drosophila population and the Ni pattern are presented in Figure 3. Having been averaged, the deviations of the individual plateau levels Di calculated from 493 individual scores yields a two-part pattern. Left of the central point the differences are averaged to give 0, whereas to the right they produce the curve, which is close to the exponential curve.
|
Fecundity Rate at Maturity in Drosophilaand Medfly Is Constant
Egg scores in Figure 1 represent the "noisy" components, over which the regular individual patterns are superimposed. But how flat is the noisy reproductive plateau at maturity in individuals, on average? Figure 3 is too rough to see the details of the plateau structure. To answer that question, we estimate individual errors between the random scores in each fly and its plateau level, and then average them over the population. This procedure allows us to evaluate how close these scores are, on average, to their regular patterns.
For each day of maturity we average only the errors of the mature individuals. Thus, the average gives
|
|
Errors EMi for 493 Drosophila and 1170 Medflies are presented in Figure 4. It can be seen that there are no maxima in the patterns on these graphs. This means that maximal values in the scores are evenly scattered over the maturity period. Thus individual "noisy" fecundity allows for a flat approximation at maturity. The patterns show the number of mature flies, and this number smoothly diminishes from the initial value because the flies leave the "steady-state" mature group after having achieved the onset of senescence. Thus, we present experimental confirmation of our theoretical finding that the flat maturity plateau exists in the life course of individual fruit flies of the species Drosophila and medfly.
|
To answer the question, one needs to calculate the individual errors, compare them at each age with the respective value of the decreasing exponential curve, and then average the errors. The results will be conditional, based on the premise that the fly is at the senescence stage of life on the day of averaging. Thus one has
|
|
sen is the rate of senescence in the jth fly. The other notations are the same as in Equation 2. The errors ESi are presented in Figure 5, with the x-axis showing the day of senescence. All onsets of senescence are placed at the origin. No specific features are seen in the patterns for both species. This finding definitely confirms our previous assumption about the exponential deceleration of individual fecundity at senescence.
|
Individual Longevity and Fecundity
The relationship between individual fecundity and senescence is intriguing. Indeed, it would be interesting to prospectively know which flies will live long. What are the individual parameters that differ in flies and thus make it possible to distinguish the long-living ones from the others in a population?
Unfortunately, individual data on Drosophila and medfly fecundity are stochastic and heterogeneous, thus preventing a direct exploration of the relationships between fecundity and longevity. The correlation coefficients between the parameters in both populations (Drosophila and 1000-fly Carey's medfly) are presented in Table 1.
|
sen, in both cases, is rather high: r = 0.632 in Drosophila and r = 0.423 in medfly [the correlation coefficients were calculated for the reversed parameters,
tail = (
sen)1]. Additionally, r = 0.474 between the life span and length of maturity period in medfly. To overcome difficulties related to the extreme stochasticity in the data, we divide both total populations into 10-decile clusters each. Each cluster contains one-tenth of the total scores (50 in Drosophila and 94 in medfly, with the exception of the first cluster in each) sorted in ascending order of LS. Then we calculated the correlation coefficients between the cluster parameters.
Drosophila and medfly data are presented in Table 2. One can see that a correlation arises, resulting in higher values of all coefficients. The most striking are 0.96 and 0.95 values of the correlation coefficients between
sen and life span in Drosophila and medfly; another is the 0.98 value of the medfly correlation between T and LS. The other large value (correlation coefficient 0.91 between T and
sen in medfly) is not essential for our analysis.
|
The overall dependencies between the mean decile values of the individual parameters and the mean decile life spans are shown in Figure 6. We can see that, indeed, neither RC nor T is related to prolongation of life span in Drosophila. Only the exponent
sen is clearly related to prolonged senescence. The smaller this exponent is, the slower is the decrease in egg-laying after the onset of senescence, and the longer is life span. The maturity period and the reproductive capacity tend to be constant in all clusters. This means that longevity in Drosophila is associated only with the rate of senescence.
|
sen, and the length of the maturity stage, T.
|
Representative examples taken from different clusters of Drosophila and medflies confirm this finding. Four such patterns from the Drosophila population are presented in Figure 8. Flies from clusters 2, 5, 8, and 10 have a similar plateau period but drastically differ in the lengths of the senescence stage. Analogous examples from medflies (data not shown) demonstrate that the length of both the maturity stage and the senescence stage increases from one cluster to the following one.
|
| DISCUSSION |
|---|
|
|
|---|
To describe the three-stage pattern, we used a set of parameters. These are the onset of reproduction, Xonset, the height of a fecundity plateau at maturity, RC, the length of the maturity period, T, and the rate of senescence,
sen. The last parameter characterizes the rate of the exponential decrease of fecundity in a senescing individual. The egg-laying rate diminishes each day as related to the current one by exp[
sen]. In this study, the rate of the senescence parameter was used instead of
tail, namely
sen = (
tail)1 used previously.
Due to the high degree of heterogeneity, a special technique was applied to test if the onset of senescence really exists, namely, whether the rate of egg laying at the maturity stage is constant, and if fecundity at senescence exponentially decelerates.
That analysis confirmed that individual fecundity plateaus, on average, are flat despite the existence of the high "noise" level observed in individual scores of fecundity. Because egg-laying shows a rhythmic activity (34), at least part of this noise may be caused by the stroboscopic effect of daily periodic cycles. An obvious maximum was not observed at maturity, on average. This means that sporadic maxima in the individual patterns are evenly scattered over the maturity ages.
The changes in the mathematical parameters described herein have predictive value and are interesting, but one would like to know the nature of the genetic and physiological events underlying these mathematically described transitions. At the onset of senescence, Xsen, the slow diminution in fecundity undergoes a drastic decrease, and this can likely be related to an alteration of the signals regulating the development of the ripening eggs in the ovaries. It is known that fecundity and life span in Drosophila are under the joint control of the insulin-like signaling pathway (35,36) and the nutritional status of the organism (37). These two inputs modulate the animal's relative hormone levels (e.g., juvenile hormone [JH] and 20-hydroxyecdysone [20E]). Well fed flies with a good nutritional status have higher levels of JH and lower levels of 20E, whereas starved flies have the opposite situation (38). These hormones regulate the expression of various genes, including the various isoforms of the Broad-Complex (BR-C) gene (38). The high JH/low 20E levels characteristic of well fed flies activate the Z1 isoform of the BR-C gene during stages 89 of egg development; this isoform represses apoptosis and allows the eggs to develop normally. The reverse situation, characteristic of flies with a nutritional shortage, activates BR-C isoforms Z1, Z2, and Z3 during that same checkpoint stage; this expression pattern subsequently induces apoptosis in that egg chamber. Thus it seems that the hormones known to affect fly fecundity (39) operate by changing the number of immature eggs allowed to develop in each egg chamber (ovariole) of the ovary. Individual egg chambers autonomously sense the nutritional status of the animal, some entering apoptosis and others not. This sensing then permits a graded fecundity response from each individual female. Integrating our observations with the existing genetic and physiological data allows the generalization that the individual fecundity plateau occurs when nutrition is good, the insulin-like signaling pathway is activated, JH level is high, ecdysone level is lower (but present), BR-C Z1 isoform is expressed, and apoptosis is repressed. Under these conditions, characteristic of the mature female adult, almost all eggs that enter the ovarioles are allowed to develop. However, when functional decline of nutritional mechanisms and degeneration of the soma as a whole arise due to individual senescence, then the nutritional status of the organism could change to yield internal conditions somewhat similar to those observed in the experimentally starved flies in the Terashima and Bownes study (38). Then JH level drops, ecdysone level increases, BR-C isoform Z1 is activated, and apoptosis is activated. Because a constantly increasing proportion of developing eggs are signaled to undergo apoptosis, significantly fewer eggs are laid, and the reproductive senescence stage begins.
The exponentially decelerated decline of fecundity is characterized by the exponent, or the rate of senescence,
sen. Mathematically, it means that fecundity in senescing flies is described by the exponential function RC · exp[
sen · (x Xsen)], where x is age (x > Xsen), RC is steady-state fecundity at maturity, and Xsen is the onset of senescence. The rate of senescence,
sen, is the exponent in the exponential function.
No specific dynamic properties were observed in the senescence-related parts of the fecundity patterns. This statement is even stronger than the thesis about the constancy of the egg-laying rate at the plateau stage. It calls into question the claim by Rauser and colleagues (25) that a late-life fecundity plateau exists. In any case, this question must be studied more thoroughly, especially when keeping in mind that Rauser and colleagues have studied population and not individual fecundity, and have presented serious evolutionary reasons for why mean fecundity values must plateau at very advanced ages.
Individual Senescence in Drosophila and Medfly
The specificity of the individual fecundity patterns in the two species poses a number of questions. What are the genetic properties related to senescence in Drosophila and medfly? Are aging and senescence similar in the two species?
In search of an answer to those questions, let us draw schematically the fecundity patterns in each species (Figure 9). Figure 9A shows six cluster-related patterns in Drosophila, and Figure 9B shows six such patterns for medflies. The schemes differ in that in Drosophila only the senescence phase changes with a cluster number whereas the maturity stage is constant. In medflies, both the length of the maturity phase and the exponent,
sen, change. The height of the reproductive plateau at maturity is shown as constant.
|
Thus it is clear that longevity is associated with different parameters; namely, in Drosophila with
sen only, and in medflies with T and
sen. These parameters may present a solid basis for prediction of life spans in individual flies, especially in Drosophila.
Senescence in Drosophila and medfly follows different scenarios. At the moment, the difference between the two species cannot be empirically explained. However, if we assume that the same genetic and physiological mechanisms known to be operative in Drosophila are also operative in medfly, albeit with species-specific variation in particular regulatory steps, then we may consider a possible molecular basis for this difference in their kinetics of reproductive senescence. Both mechanisms probably begin to work at the onset of senescence. The Drosophila process probably involves the mechanism determining the ability of each egg in the ovary to choose between development and apoptosis (38). The medfly process might well be a diminution of the number of eggs entering the ovary from germarium. A different role of these mechanisms can produce the observed effects in the two species. This topic will be the subject of a future investigation.
Thus it is possible that, in Drosophila, genetic mechanisms might regulate the length of the maturity stage, sustaining it at some constant period. In humans, the length of the period between menarche and menopause is under strict genetic control (40). Should our supposition about Drosophila be confirmed, then such a fact would add Drosophila to the two species known for such a trait, namely humans and pilot whales (41). If so, then perhaps we might in the future use Drosophila not only as a model for human diseases (42) but also for reproductive processes.
Mair and colleagues (43) have demonstrated that the life-span extension observed in Drosophila females subjected to dietary restriction does not seem to be caused by a reduction in vitellogenesis or ovarian activity. Thus the processes regulating reproductive senescence in Drosophila seem to be separable from those influencing longevity and somatic senescence. Given that finding, it is reasonable to assume that reproductive senescence is regulated in part by those processes centered on the expression of the BR-C and related genes, whereas the longevity and somatic senescence processes may be influenced by the organism's stress-resistance processes centered on the expression of the dFOXO transcription factor and its downstream targets. If so, then the reproductive and somatic senescent processes may be indirectly related within the individual organism by the power available for reproduction, Smax [see Figure 3 of (28)], which is itself related to overall oxygen consumption rate and oxidative vulnerability (29).
Individual and Population Aging and Senescence
Aging, senescence, and reproduction are inevitable features of living beings (43,44). Following the definition by Rose (44), aging is a persistent decline in age-specific components of an organism due to internal physiological deterioration, which is accompanied by an increase in the risk of death. Decline of reproduction (4548) and physical activity (19,4951) and age-related changes in accumulation of reserves (52,53) usually characterizes aging in individuals. Generally, such functional decline is used to present senescence.
In this study we definitely discriminate mortality rate and individual senescing. The Gompertz equation,
|
|
It is common belief that the mortality-related Gompertz parameter b may be considered to be a measure of the demographic "rate of aging." The slope of the logarithm of the mortality curve is often associated with changes in individual aging. Nonetheless, such interpretation may be incorrect: The changes in the slope of this curve do not necessarily correspond to the individual changes in aging. For example, the improvement in survival in developed countries in the first half of the past century is characterized by a rectangularization of the survival curve. This trend corresponds to an increase in the slope of the respective logarithm of the mortality rate. Naturally, this does not mean that the rate of individual aging increases as well (54).
To describe individual aging, we use the homeostatic approach, in which aging is treated as an age-related decrease of the homeostatic capacity (55,56). Such a decrease continues until the death of the individual. The homeostatic capacity S(x), a basic notion of the approach, describes the ability of the systemic mechanisms in the organism to convert into Adenosine triphosphate (ATP) substances delivered from external sources. Aging is governed by a quasiexponential function
|
|
|
|
|
|
This means that the rate of individual aging of the reproductive system becomes approximately constant, equal to
sen. The aging process continues all the time, but only the onset of reproductive senescence made it explicit. A diminution of reproduction after the onset of senescence can be observed in individual fecundity patterns.
Thus, in fruit flies we relate senescence with reproductive changes. Senescence arises when aging results in an inability of the mechanisms supporting the reproductive system to fulfill the endowed reproductive program, sustaining the constant rate of egg-laying. Being closely related to mechanistic causes of senescence (28), the
sen parameter can be treated as a general characteristic of senescence in fruit flies.
A large difference between the Gompertz-based measure of aging, b, and the parameter
sen is that the first one describes the demographic rate of population aging, whereas the second is directly related to senescing in an individual. The first one characterizes the processes at all ages, whereas
sen describes only the last period of life history, its senescence stage.
In general, the parallel study of the two species suggests that the three-stage approach is applicable individually to a wide range of animals. The parameters describing individual fecundity change during laboratory selection. Existing genetic and physiological data identify the molecular mechanisms most likely to be responsible for the modeled changes in fecundity and longevity of these two Dipteran species. It is possible to use the technique to analyze individual fecundity in other Drosophila studies as well as in differing species such as ladybird beetles (24).
| Acknowledgments |
|---|
We thank James W. Vaupel for the opportunity to complete this work at the Max-Planck Institute for Demographic Research. We also thank A. I. Mikhalski, A. A. Butov, and A. A. Romaniukha for discussion on the paper and for valuable comments.
| Footnotes |
|---|
Received October 22, 2004
Accepted February 8, 2005
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
| ||||||||||||||||||||||||||||||||
| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|