- Short report
- Open Access
Demographic analysis reveals gradual senescence in the flatworm Macrostomum lignano
Frontiers in Zoology volume 6, Article number: 15 (2009)
Free-living flatworms ("Turbellaria") are appropriate model organisms to gain better insight into the role of stem cells in ageing and rejuvenation. Ageing research in flatworms is, however, still scarce. This is partly due to culture difficulties and the lack of a complete set of demographic data, including parameters such as median lifespan and age-specific mortality rate. In this paper, we report on the first flatworm survival analysis. We used the species Macrostomum lignano, which is an emerging model for studying the reciprocal influence between stem cells, ageing and rejuvenation. This species has a median lifespan of 205 ± 13 days (average ± standard deviation [SD]) and a 90th percentile lifespan of 373 ± 32 days. The maximum lifespan, however, is more than 745 days, and the average survival curve is characterised by a long tail because a small number of individuals lives twice as long as 90% of the population. Similar to earlier observations in a wide range of animals, in M. lignano the age-specific mortality rate increases exponentially, but levels off at the oldest ages. To compare the senescence of M. lignano with that of other ageing models, we determined the mortality rate doubling time, which is 0.20 ± 0.02 years. As a result, we can conclude that M. lignano shows gradual senescence at a rate similar to the vertebrate ageing models Rattus norvegicus and Mus musculus. We argue that M. lignano is a suitable model for ageing and rejuvenation research, and especially for the role of stem cells in these processes, due to its accessible stem cell system and regeneration capacity, and the possibility of combining stem cell studies with demographic analyses.
Flatworms have been an object of ageing studies since Child's initial investigations [1, 2]. Researchers tended to focus on the role of stem cells and cell renewal during ageing, and the causal effect of regeneration and starvation on rejuvenation [2–4]. Despite these fascinating themes, the extent of flatworm ageing research remained limited in comparison to that of other model organisms such as Caenorhabditis elegans, Drosophila melanogaster and rodents. The lack of detailed demographic data partly accounts for this, as the only available data are the maximum lifespans of several species. These data, however, include many discrepancies due to non-specified or non-standardised culture conditions or culture problems such as the presence of fungal and bacterial contaminations [3, 4]. Without a basic set of demographic data, the most fundamental question – at which age can an individual be considered old? – remains unanswered. As a result, it is hard to draw any conclusions about, for example, old-age regeneration capacity or the rate of cell renewal as a function of age. Previously published data about these issues are often contradictory or ambiguous and there is still little known about senescence, rejuvenation and the causes of death in flatworms [3, 4]. This demonstrates that establishing a survival curve, median lifespan and 90th percentile lifespan is a prerequisite for the experimental design of ageing studies and should be the first step in initiating ageing research with a new model organism. Emerging ageing models are often first described demographically, after which detailed studies follow, stemming from these initial descriptions [5, 6]. Because lifespan parameters indicate when individuals can be considered young or old, they allow for choosing age groups to study biomarkers as a function of age and for experiments in which young and old worms are studied comparatively. Furthermore, the survival curve indicates what proportion of the initial cohort is alive at a certain age. Therefore, it can be used to calculate how large an initial culture set-up is needed to retain individuals at a desired age to give the experiment enough statistical power. Besides lifespan parameters, data about the age-related changes in mortality rate provide a basic measure for the rate of senescence , and can be used to study rejuvenation by experimental manipulation, such as regeneration and caloric restriction.
In this manuscript, the first flatworm survival curve and demographic dataset are presented. We used Macrostomum lignano (Rhabditophora, Platyhelminthes), which is a new model for stem cell biology, development, regeneration and the study of sexual selection [8–13], as well as an emerging model for ageing and rejuvenation research, and especially for the role of stem cells in these processes . Egger et al. suggested that, in M. lignano, repeated regeneration induces a lifespan extension and possible rejuvenation [8, 9, 11], because individuals were cut up to 59 times over a period of 26.5 months and were still able to regenerate . The demographic data in this manuscript can be used to design experiments in which mortality rate parameters are used to conclude whether repeated regeneration slows down the ageing process or induces an actual rejuvenation in comparison to uncut individuals of M. lignano.
The data presented here were obtained in three replicate cohorts, maintained under controlled culture conditions. In M. lignano, the average (± SD) median lifespan is 205 ± 13 days, the 90th percentile lifespan 373 ± 32 days, and the maximum lifespan more than 745 days (2.04 years). While all individuals in replicate 1 and 2 were already dead at this age, 4% of the individuals in replicate 3 are still alive (figure 1). The lifespan data of the separate replicates are listed in table 1, and individuals of four, 27 and 76 weeks old are presented in figure 2 to illustrate the morphological changes as a function of age. The maximum lifespan here observed is considerably longer than 42 weeks (0.8 year), which was previously reported by Egger et al. . These authors noted, however, that culture limitations like bacteria or parasites may have shortened the lifespan of the specimens. It is remarkable that the maximum lifespan of 745 days is twice as long as the 90th percentile lifespan, which is visualised by the long tail in the survival curve (figure 1), and more than three times the median lifespan, indicating a levelling off of the mortality rate acceleration at advanced ages (see below). This phenomenon was previously described in large medfly cohorts by Carey et al. . Furthermore, these authors concluded that medflies appear not to have a characteristic lifespan, because they observed different maximum lifespans in the several cohorts. Because the maximum lifespan is based on only one individual in the cohort, it results in a great variability which can also be observed in M. lignano (table 1). Therefore it is better to use more reliable statistics such as the median lifespan and 90th percentile when demographic data are compared among laboratories, or in experiments studying lifespan extension and rejuvenation.
Lifespan itself, however, cannot indicate the characteristics of senescence, and therefore age-related mortality rate acceleration should be studied in addition. In a wide range of animal populations, including human, an exponential increase in the age-specific mortality rate is observed after maturation and is considered a hallmark of senescence. However, the increase in mortality rate decelerates at advanced ages in several species [5, 15]. A similar mortality pattern can be observed in M. lignano (figure 3). There are several possible explanations for this, the first one being the genetic heterogeneity of the population. As the population ages, frailer individuals with higher death rates will die out in greater numbers than those with lower death rates, thereby transforming the population into one consisting mostly of individuals with low death rates [5, 15, 16]. A second explanation is that the chance of death for individuals actually levels off at older ages, resulting in increasing life expectancies at the oldest ages [5, 15, 16].
To compare the rate of senescence of M. lignano with that of other ageing models, the Gompertz function was plotted (table 2). The Gompertz function describes the exponential increase in the age-specific mortality rate, and is the most common model used for interspecies comparison [17–19]. It has proven to be especially useful in simply but accurately describing the kinetics of ageing of populations . Furthermore, the Gompertz coefficient G allows for calculating the mortality rate doubling time (MRDT), which is held to be a fundamental measure of senescence [7, 20]. Based on the MRDT, Finch  characterised senescence using a continuum with three general subdivisions according to the rate of degenerative changes: rapid (0.005–0.10 years), gradual (0.10–more than 10 years), and negligible (no indication of senescence). Humans, for example, show gradual senescence, with a MRDT of 8 years  and maximum lifespan of 122 years . M. lignano has a MRDT of 0.20 ± 0.02 years and is therefore also a gradually ageing species. This is in contrast to the commonly used invertebrate ageing models, Caenorhabditis elegans and Drosphila melanogaster. Both show rapid senescence, with a MRDT of 0.02–0.04 years and a short maximum lifespan of 0.16 and 0.30 years respectively . The vertebrate models Mus musculus and Rattus norvegicus show gradual senescence, with MRDTs of 0.27 and 0.30 years respectively [7, 20]. The MRDT of M. lignano is thus much more similar to the rodent models than the invertebrate models. Although M. musculus and R. norvegicus have a similar rate of senescence, their lifespan is longer: 4.5 and 5.5 years respectively . The median lifespan ranges from about 20 to 30 months in M. musculus  and is around 30 months for several strains of R. norvegicus , while it is 205 ± 13 days (between 6 and 7 months) in M. lignano. This is due to a larger initial mortality rate (A0) in M. lignano (average ± SD: 0.42 ± 0.06/year) than in the lab mouse and rat (0.03 and 0.02/year, respectively) . However, we cannot exclude that this initial mortality rate may slightly decrease when culture conditions could be further optimised in this young model.
Studying the influence of treatment on the rate of senescence, e.g. induction of rejuvenation through regeneration, is possible by comparing mortality parameters among treatment replicates . Mortality parameters such the s parameter, describing the degree of deceleration in mortality at older ages, and the age-independent Makeham constant are however not included in the Gompertz function . Therefore, one should first determine the mortality model which best describes the data. Analysis of the data obtained here demonstrates that the average overall survival curve and replicate 1 and 3 are best described by the logistic function, which contains the s parameter, while replicate 2 is best described by the Gompertz function (table 3 and 4).
The short lifespan in comparison to other gradually ageing models and the ease of culturing makes M. lignano a suitable ageing model, but the major advantage is the presence of a very accessible stem cell system. The stem cells can be visualised and quantified in vivo, and also manipulated in several ways. Stem cells can be arrested in mitosis or in S-phase by adding colchicine or hydroxyurea respectively to the medium, or eliminated through irradiation. In contrast, the proliferation rate can also be increased by injuring the individual, which induces regeneration [11, 24–28]. Furthermore, this manuscript demonstrates that it is possible to perform demographic analyses, which are necessary to draw reliable conclusions when stem cells are studied at different ages. Because the accessible stem cell population of M. lignano can be combined with demographic studies, this species has the potential to play a key role in obtaining a better understanding of stem cell biology in tissue homeostasis, ageing, and even rejuvenation.
M. lignano is a free-living, hermaphrodite flatworm with a five-day embryonic development and a generation time of about 2–3 weeks . M. lignano was first found in 1995 in Lignano, Italy, resulting in the first lab cultures. Afterwards, it was resampled several times, and cultures have been established in several labs. [10, 11]. The species and its natural habitat is described in detail by Ladurner et al. [10, 11].
In the lab, M. lignano is easily cultured in f/2, a nutrient-enriched artificial seawater medium at a salinity of 32‰ , and incubated at 20°C with a 60% relative humidity and a 13 h:11 h light: dark cycle . Individuals are fed ad libitum with the diatom Nitzschia curvilineata, which is grown under identical conditions as the worms and can be obtained from the culture collection of algae (SAG) at the University of Göttingen (strain 48.91, http://sagdb.uni-goettingen.de) .
Survival was followed in three replicate cultures consisting of 100 individuals each. To initiate the replicates, one-day-old juveniles were put in separate wells of a 12-well plate and, by maintaining the worms individually, reproduction and hence a mixture of individuals of different ages could be avoided. About every 30 days, individuals were put into new 12-well plates with new culture medium and diatoms to maintain ad libitum food resources. The number of survivors was counted about every 10 days. As in demographic studies of other model organisms [31, 32], animals that died due to age-independent injury (for example rupture in C. elegans) were censored to ensure that the analysis reflects the natural lifespan. In M. lignano, age-independent death can be caused by infection with Thraustochytrium caudivorum . This can be recognised by the characteristic dissolution of the tail plate . In previous ageing cultures, it was observed at different ages (ranging from three weeks old to the oldest individuals in culture), which allowed us to conclude that it is age-independent. In the meantime, we were able to establish parasite-free cultures by optimising the working procedures and using the Triton treatment presented by Schärer et al. . Kaplan-Meier survival curves were constructed (for completeness, non-censored survival curves are also given in additional file 1), and the median lifespan (50% mortality), 90th percentile lifespan (90% mortality), and maximum lifespan were determined.
Characterising senescence was done by calculating the age-specific mortality rate and determining the mortality parameters of both the Gompertz and the Logistic model. The equations of these models are y(t) = A0eGt and respectively . Identifying which demographic model best fits the data and determining the mortality parameters was done by using WinModest software. This software uses the maximum likelihood method, which is based on the age distribution of deaths. The method provides better parameter estimates that are more consistent and less influenced by technical aspects of the experimental design such as sample size than those of other methods [17, 18]. The WinModest software is very straightforward in use and is made freely available by Dr. Pletcher . To determine the parameters of the Logistic model, we used the complete survival dataset of the three replicates. Because the Gompertz model describes the exponential increase in mortality rates, we used the dataset until Day 367. We chose this subset because at Day 367 the average 90th percentile was reached, and after this day there is an obvious deceleration of the increasing mortality rate. The MRDT was calculated using the formula MRDT = ln(2)/G [7, 20].
Child CM: A study of senescence and rejuvenenscence based on experiments with Planaria dorotocephala. Archiv fiir Entwicklungsmechanik der Organismen. 1911, 31: 537-616. 10.1007/BF02162177.
Child CM: Senescence and Rejuvenescence. 1915, Chicago: The University of Chicago Press
Haranghy L, Balázs A: Ageing and rejuvenation in planarians. Experimental Gerontology. 1964, 1: 77-91. 10.1016/0531-5565(64)90009-9.
Lange CS: A possible explanation in cellular terms of the physiological ageing of the planarian. Experimental Gerontology. 1968, 3: 219-230. 10.1016/0531-5565(68)90005-3.
Carey JR, Liedo P, Orozco D, Vaupel JW: Slowing of mortality-rates at older ages in large medfly cohorts. Science. 1992, 258: 457-461. 10.1126/science.1411540.
Valdesalici S, Cellerino A: Extremely short lifespan in the annual fish Nothobranchius furzeri. Proceedings of the Royal Society of London Series B-Biological Sciences. 2003, 270: S189-S191. 10.1098/rsbl.2003.0048.
Finch CE: Longevity, Senescence and the Genome. 1990, Chicago: The University of Chicago Press
Egger B, Ladurner P, Nimeth K, Gschwentner R, Rieger R: The regeneration capacity of the flatworm Macrostomum lignano – on repeated regeneration, rejuvenation, and the minimal size needed for regeneration. Development Genes and Evolution. 2006, 216: 565-577. 10.1007/s00427-006-0069-4.
Egger B: Regeneration: Rewarding, but potentially risky. Birth Defects Research (Part C). 2008, 84: 257-264. 10.1002/bdrc.20135.
Ladurner P, Schärer L, Salvenmoser W, Rieger RM: A new model organism among the lower Bilateria and the use of digital microscopy in taxonomy of meiobenthic Platyhelminthes: Macrostomum lignano, n. sp. (Rhabditophora, Macrostomorpha). Journal of Zoological Systematics and Evolutionary Research. 2005, 43: 114-126. 10.1111/j.1439-0469.2005.00299.x.
Ladurner P, Egger B, De Mulder K, Pfister D, Kuales G, Salvenmoser W, Schärer L: The stem cell system of the basal flatworm Macrostomum lignano. Stem cells: from Hydra to man. Edited by: Bosch ThCG. 2008, Berlin – Heidelberg – New York: Springer, 75-94.
Morris J, Nallur R, Ladurner P, Egger B, Rieger R, Hartenstein V: The embryonic development of the flatworm Macrostomum sp. Development Genes and Evolution. 2004, 214: 220-239. 10.1007/s00427-004-0406-4.
Scharer L, Sandner P, Michiels NK: Trade-off between male and female allocation in the simultaneously hermaphroditic flatworm Macrostomum sp. Journal of Evolutionary Biology. 2005, 18: 396-404. 10.1111/j.1420-9101.2004.00827.x.
Mouton S, Willems M, Braeckman BP, Egger B, Ladurner P, Schärer L, Borgonie G: The free-living flatworm Macrostomum lignano: a new model organism for ageing research. Experimental Gerontology. 2009, 44: 243-249. 10.1016/j.exger.2008.11.007.
Vaupel JW, Carey JR, Christensen K, Johnson TE, Yashin AI, Holm NV, Iachine IA, Kannisto V, Khazaeli AA, Liedo P, Longo VD, Zeng Y, Manton KG, Curtsinger JW: Biodemographic trajectories of longevity. Science. 1998, 280: 855-860. 10.1126/science.280.5365.855.
Curtsinger JW, Fukui HH, Townsend DR, Vaupel JW: Demography of Genotypes – Failure of the Limited Life-Span Paradigm in Drosophila melanogaster. Science. 1992, 258: 461-463. 10.1126/science.1411541.
Pletcher SD: Model fitting and hypothesis testing for age-specific mortality data. Journal of Evolutionary Biology. 1999, 12: 430-439. 10.1046/j.1420-9101.1999.00058.x.
Promislow DEL, Tatar M, Pletcher S, Carey JR: Below threshold mortality: implications for studies in evolution, ecology and demography. Journal of Evolutionary Biology. 1999, 12: 314-328. 10.1046/j.1420-9101.1999.00037.x.
Yen K, Steinsaltz D, Mobbs CV: Validated analysis of mortality rates demonstrates distinct genetic mechanisms that influence lifespan. Experimental Gerontology. 2008, 43: 1044-1051. 10.1016/j.exger.2008.09.006.
Arking R: Biology of Aging. 1998, Sunderland: Sinauer Associates Inc, 2
Robine JM, Allard M: The oldest human. Science. 1998, 279: 1834-1835. 10.1126/science.279.5358.1831h.
Austad SN, Kristan DM: Are mice calorically restricted in nature?. Aging Cell. 2003, 2: 201-207. 10.1046/j.1474-9728.2003.00053.x.
Altun M, Bergman E, Edstrom E, Johnson H, Ulfhake B: Behavioral impairments of the aging rat. Physiology & Behavior. 2007, 92: 911-923. 10.1016/j.physbeh.2007.06.017.
Bode A, Salvenmoser W, Nimeth K, Mahlknecht M, Adamski Z, Rieger RM, Peter R, Ladurner P: Immunogold-labeled S-phase neoblasts, total neoblast number, their distribution, and evidence for arrested neoblasts in Macrostomum lignano (Platyhelminthes, Rhabditophora). Cell and Tissue Research. 2006, 325: 577-587. 10.1007/s00441-006-0196-2.
Ladurner P, Rieger R, Baguna J: Spatial distribution and differentiation potential of stem cells in hatchlings and adults in the marine platyhelminth Macrostomum sp.: A bromodeoxyuridine analysis. Developmental Biology. 2000, 226: 231-241. 10.1006/dbio.2000.9867.
Nimeth KT, Mahlknecht M, Mezzanato A, Peter R, Rieger R, Ladurner P: Stem cell dynamics during growth, feeding, and starvation in the basal flatworm Macrostomum sp (Platyhelminthes). Developmental Dynamics. 2004, 230: 91-99. 10.1002/dvdy.20035.
Pfister D, De Mulder K, Philipp I, Kuales G, Hrouda M, Eichberger P, Borgonie G, Hartenstein V, Ladurner P: The exceptional stem cell system of Macrostomum lignano: Screening for gene expression and studying cell proliferation by hydroxyurea treatment and irradiation. Frontiers in Zoology. 2007, 4: 9-10.1186/1742-9994-4-9.
Schärer L, Ladurner P, Rieger RM: Bigger testes do work more: experimental evidence that testis size reflects testicular cell proliferation activity in the marine invertebrate, the free-living flatworm Macrostomum sp. Behavioral Ecology and Sociobiology. 2004, 56: 420-425. 10.1007/s00265-004-0802-9.
Guillard R, Ryther J: Studies of marine planktonic diatoms. I. Cyclotella nana Hustedt, and Detonula confervacea (Cleve) Gran. Canadian Journal of Microbiology. 1962, 8: 229-239.
Egger B, Ishida S: Chromosome fission or duplication in Macrostomum lignano (Macrostomorpha, Plathelminthes) – remarks on chromosome numbers in 'archoophoran turbellarians'. Journal of Zoological Systematics and Evolutionary Research. 2005, 43: 127-132. 10.1111/j.1439-0469.2005.00300.x.
Andersen JB, Li XL, Judge CS, Zhou A, Jha BK, Shelby S, Zhou L, Silverman RH, Hassel BA: Role of 2–5A-dependent RNase-L in senescence and longevity. Oncogene. 2007, 26: 3081-3088. 10.1038/sj.onc.1210111.
Crawford D, Libina N, Kenyon C: Caenorhabditis elegans integrates food and reproductive signals in lifespan determination. Aging Cell. 2007, 6: 715-721. 10.1111/j.1474-9726.2007.00327.x.
Schärer L, Knoflach D, Vizoso DB, Rieger G, Peintner U: Thraustochytrids as novel parasitic protists of marine free-living flatworms: Thraustochytrium caudivorum sp nov parasitizes Macrostomum lignano. Marine Biology. 2007, 152: 1095-1104. 10.1007/s00227-007-0755-4.
We want to thank S.D. Pletcher for his useful suggestions concerning the use of WinModest. We also want to thank W. Govaerts and A.L. Kheibarshekan for helping us to get more insight into the mathematics of model fitting. S. Mouton, M. Willems, and P. Back received funding from the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen).
The authors declare that they have no competing interests.
SM designed the study, carried out the experiments, and wrote the manuscript. MW co-designed the study and participated in drafting the manuscript. PB helped to analyse and interpret the results. BPB helped to analyse the results and contributed to the manuscript drafting. GB critically revised the manuscript. All authors read, edited and approved the final manuscript.