#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Comparing growth patterns of three species: Similarities and differences


Authors: Norbert Brunner aff001;  Manfred Kühleitner aff001;  Werner Georg Nowak aff001;  Katharina Renner-Martin aff001;  Klaus Scheicher aff001
Authors place of work: University of Natural Resources and Life Sciences, Department of Integrative Biology and Biodiversity Research, Vienna, Austria aff001
Published in the journal: PLoS ONE 14(10)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0224168

Summary

Quantitative studies of the growth of dinosaurs have made comparisons with modern animals possible. Therefore, it is meaningful to ask, if extinct dinosaurs grew faster than modern animals, e.g. birds (modern dinosaurs) and reptiles. However, past studies relied on only a few growth models. If these models were false, what about the conclusions? This paper fits growth data to a more comprehensive class of models, defined by the von Bertalanffy-Pütter (BP) differential equation. Applied to data about Tenontosaurus tilletti, Alligator mississippiensis and the Athens Canadian Random Bred strain of Gallus gallus domesticus the best fitting growth curves did barely differ, if they were rescaled for size and lifespan. A difference could be discerned, if time was rescaled for the age at the inception point (maximal growth) or if the percentual growth was compared.

Keywords:

Birds – Curve fitting – Reptiles – Poultry – Chickens – Dinosaurs – Basal metabolic rate measurement – Chicken models

Introduction

Mathematical growth models aim at a simplified description of growth in terms of curves that fit well to size-at-age data [1]. As the growth of animals depends on multiple factors, the most-informative data came from controlled studies, e.g. for chicken [2] or for pigs [3]. By contrast, for wildlife and wild-caught fish, there remained considerable uncertainties about the proper choice of the growth model [4]. Dealing with extinct species the situation was even worse, as no weighing of body mass was possible for fossils. Nevertheless (e.g. Table 1), recent approaches led to mathematical growth models for dinosaurs [5] that have “revolutionized our understanding of dinosaur biology” [6]. For instance, it is now consensus that dinosaurs grew faster than modern reptiles.

Tab. 1. Age and mass data of Tenonotosaurus tilletti from [18].
Age and mass data of Tenonotosaurus tilletti from [<em class="ref">18</em>].

However, previous growth studies relied on few models only, whence model uncertainty may be an issue for the comparisons of growth curves of different species. (This paper compares several thousand models, as outlined in Fig 1). Another issue is scaling-up. [7] defined dimensionless mass and time ratios and concluded from a plot that the so rescaled growth data of 13 species were close to a “universal growth curve”. Thus, aside from the different scaling, all animals would grow in the same way. We therefore reconsider the conceptual question, how to compare the growth of species that differ in size and life span.

Fig. 1.
Named models (blue) and part of the search-region (yellow) for the exponent-pair of the best fitting growth model.

With respect to model uncertainty, this paper studies a differential Eq (1) of [8], that includes e.g. the models of Brody [9], von Bertalanffy [10], or Verhulst [11] as special cases. In view of its generality it allows a more accurate comparison of growth curves. Literature refers to this model as Bertalanffy-Pütter (BP) model; c.f. [12].


Eq (1) describes growth of mass m(t) at time t and it uses five free parameters that are optimized to obtain a growth curve with a best fit to given data. [13] related the non-negative exponent-pair a <b to the metabolism. The non-negative constants p and q are scaling constants. The fifth parameter is the initial condition of the differential equation that is needed to determine the function m(t); e.g. m(0) = m0, where m0 > 0 is an estimate for the hatching (or natal) mass.

In this paper we interpret Eq (1) as a definition for a class of models; the BP-class. Thereby, each exponent-pair defines a unique BP-model. For instance, the Verhulst-model is defined from the exponent-pair a = 1, b = 2. Each BP-models has three free parameters (p, q and m0). Fig 1 illustrates this approach and it indicates the exceptional character of these named models, when compared to the range of possible models, whereby for this paper we confined the search for the best-fit exponent-pairs to the yellow area. Richards’ model [14], and the generalized Bertalanffy model of Pauly [15] are represented as line segments (i.e. subclasses of the class of BP-models). Further, the model of Gompertz [16] in the following sense is a limit-case of the BP-class: Growth-curves of the Gompertz model are limits, for (a, b) → (1, 1), of BP-growth curves with exponents a, b [17].

Lee and Werning [18] compared the growth of Tenontosaurus tilletti with the growth of modern Alligator mississippiensis and they concluded that dinosaurs (more specifically: iguanodontids) were not scaled-up lizards, as they grew much faster. We revisit this issue and seek the best fitting BP-models for their data. These data are from [18] (Table 2) about Tenontosaurus tilletti (twelve data points, mass 23–1102 kg, age 1–26 years) and the data about Alligator mississippiensis (41 data points, mass 0.1–40.7 kg, age 1–42 years), which we retrieved from a plot in [18] (using DigitizeIt of Bormisoft®).

Tab. 2. Parameters of the best-fitting models.
Parameters of the best-fitting models.

We also verified the alligator-data from the original source [19], who over a time span of forty years captured and partly recaptured ca. 7000 alligators from Louisiana, USA. In order to explore the limits of dinosaur growth, we used data about modern avian dinosaurs, broiler chicken that were bred for fast growth and reared under optimal conditions. To this end we identified the best fitting BP-model for the data from [2] (Table 1) about the Athens Canadian Random Bred strain of Gallus gallus domesticus (28 data points, mass 0.04–2.23 kg, age 0–170 days).

Within the BP-class, model uncertainty was related to the variability of the exponents. To this end, the paper identified the region of near-optimal exponent-pairs. The exponent-pairs of this region could also be used to model growth without affecting the fit to the data significantly when the other parameters were optimized. We used them to explore the model uncertainty.

The best-fitting and the near optimal exponents were then used to compare the growth of different species despite their different scales in size and age. In addition, a dinosaur-year had more days, but these were shorter. As overall a year covered about the same time span as today, we used kg and years as units; e.g. weight gains in kg/year also for chicken.

Methods

The methods are explained in detail in our preprint [20] at BioRxiv. We therefore point out only the main issues. As was observed e.g. for chicken [2], the standard deviation of mass becomes higher for heavier animals, whence the method of least squares may not be suitable for data-fitting. Instead, as in [21] we minimized the sum of squared errors between the logarithm of the growth function and the logarithmically transformed data (SSLE). This defined the following function (2):


, assuming model (1) with exponents a, b. An exponent-pair was near-optimal, if its SSLEopt(a, b) exceeded the least SSLE by less than 5%. We did not minimize SSLE for each exponent-pair. Instead we considered exponent pairs of the search region (yellow area in Fig 1) on a grid (distance 0.01 in the x and y directions, respectively). Thereby, we searched 26,200 grid-points for the chicken, 88,730 for the alligators and 42,371 for Tenontosaurus.

Results

Table 2 summarizes the model parameters that minimized SSLE. The parameters for chicken are from [21]. In order to define dimensionless coordinates, asymptotic mass mmax was computed as the limit of the growth curve m(t), when time t approaches infinity. At “full age”, tfull, 90% of the asymptotic mass were reached; we used “full age” as a proxy for “adulthood”. (Further, we used 90% to avoid excessive extrapolations, if the asymptotic mass was larger than the observations.) The inflection point is defined by the maximal growth rate mmax′; it was attained at age tinfl with mass minfl. In comparisons between species the maximal growth rate (i.e. (tinfl)) is used as a proxy for the basal metabolic rate [22]. These parameters were all computed from the best fitting model.

In order to compare the growth curves, they were rescaled in dimensionless coordinates. Ideally, the dimensionless time coordinate corresponds to about the same stage of the biological development of the considered animals. We use a linear rescaling, assuming t = 0 has the same biological meaning for the considered animals and seeking a second point of time with the same meaning. Fig 2 of [20] used full age tfull; i.e. mass was reported as a fraction of the asymptotic mass (mmax) and time was expressed as a fraction of full age. In terms of these dimensionless coordinates, the best-fitting model curves were almost equal. The rescaled data, too, were close to these curves, except for a larger spread for tenontosaurs. Similar plots were obtained, if instead of tfull the half-weight age was used (i.e. m(t) = mmax/2) or any other fraction of the asymptotic weight (e.g. 15%). As this seemed to indicate that all animals would grow alike, aside from rescaling, we checked also other empirically meaningful ages.

In Fig 2 we used tinfl. In terms of this rescaling, a difference between the species could be discerned, although the rescaled growth curves remained close together: The graphical representation of the results uses red for chicken, green for alligators and blue for Tenontosaurus. Chicken grew faster than tenontosaurs, and these grew faster than alligators, whereby some tenontosaurs (blue points) grew even faster than chicken and slower than alligators.

Fig. 2.
Growth data and best fitting growth curves in dimensionless coordinates (mass as a fraction of mmax, time as a fraction of tinfl) for chicken (red), alligators (green), and tenontosaurs (blue).

As for another comparison of the growth, in Fig 3 we compared the relative growth rates (t)/m(t) using a dimensionless time scale. This graphical representation emphasized the differences in the growth rates best: well-fed broiler chicken grew more than ten times faster than alligators and Tenontosaurus and the latter grew somewhat faster than alligators. Fig 3 displays this for the rescaling using tinfl. In [20] this was also observed using tfull for rescaling.

Fig. 3.
Growth rates relative to body mass for chicken (red), alligators (green) and tenontosaurs (blue) with time as a fraction of tinfl.

With respect to model uncertainty, Fig 4 plots the optimal and near-optimal exponent-pairs. Despite the similarity of the data in dimensionless coordinates, the optimal exponent-pairs were different. However, due to the larger variance of the dinosaur-data the region of near-optimal exponents for dinosaurs was larger and it included the regions for alligators and for chicken. Thus, judging from the perspective of extinct dinosaurs, their growth data did not display a systematic difference to modern species, whence there was no fundamental change in the growth pattern. The regions of near-optimal exponents displayed fuzzy boundaries and points close to the diagonal were not connected to the regions. This was caused by the optimization strategy, aiming at a high accuracy for points next to the diagonal and at faster computations thereafter. However, despite these deficiencies the visualization of the near-optimal exponents verified the optimal character of the optimal exponent-pairs.

Fig. 4.
Optimal and near-optimal exponent-pairs for chicken (triangle and red area), alligators (upside triangle and green area) and tenontosaurs (circle and blue area). For better orientation, the exponent-pairs of three named models are plotted (blue).

Fig 5 used the near-optimal models to explore, how sensitive the maximal growth rate was to the choice of a model. The growth rate is a measure that cannot be observed directly from the data; it is derived from a growth model and depends on what model is selected. This was demonstrated for the maximal growth rate, which varied considerably even for growth curves that fitted well to the data. The clouds in Fig 5 display the values of m and at the inflection point of m(t), using near-optimal growth curves. Apparently, even well-fitting growth curves resulted in inaccurate estimates for the maximal . Nevertheless, regardless of the near-optimal model used, (tinfl) for chicken was much higher than (tinfl) for the larger alligators, and (tinfl) for dinosaurs was largest, whereas relative to body size, i.e. in terms of (tinfl)/m(tinfl), chicken grew fastest.

Fig. 5.
Maximal growth rate, (tinfl), and mass at the inflection point, minfl = m(tinfl), for near-optimal growth curves m(t) for chicken (red), alligators (green) and tenentosaurs (blue).

Discussion

For the data about three species of dinosaurs from [18] only Tenontosaurus provided feasible data. For the two other species, the plot of the near-optimal exponent-pairs (c.f. Fig 4) displayed large regions that almost covered the search grid. As a large region of near-optimal exponents indicates that data may not carry enough information to differentiate between growth models, the paper did not use them. However, in view of the inherent uncertainties of estimating the mass of dinosaurs [6], it was surprising that one in three datasets allowed to differentiate between the models.

Further, for all species the optimal exponent-pairs were quite remote from the exponent-pairs for the named models which are more common in growth studies. This indicates that BP-models provide a significantly better fit than the conventionally used models. In fish-biology it has long been accepted that exponent-pairs (a, b) with a < 1 and b = 1 might be better compatible with biological constraints for growth; e.g. the growth of gill surface area relative to mass growth [16]. Recently, also exponents b < 1 were considered as biologically meaningful [23]. In epidemiology, too, recent publications supported the use of BP-models to analyze outbreaks of diseases, e.g. [24].

The issue of rescaling proved to be tricky. Using the “full age” to define dimensionless coordinates did not allow to discern different growth patterns for different species. Using the age at the inflection point was more satisfactory and this age may have a biological meaning (phenomenologically, growth is fastest at this age). For the present data, this rescaling resulted in the expected outcome: Broiler chicken grew fastest and Tenontosaurus grew faster than modern reptiles. This pattern was confirmed under three different perspectives (Figs 2, 3 and 5). However, using a linear transformation for rescaling may be an oversimplification, as for different species the fraction t/tinfl may correspond to different stages of their biological development. Yet, using this linear transformation was a convenient tool to combine data and growth curves into one plot. Further, with respect to Fig 4 the faster growth of broiler chicken will also be observed for any nonlinear transformation of time that aims at a proper representation of biological development.

Conclusion

It is generally acknowledged that mass-at-age estimates for dinosaurs are highly uncertain. It was therefore surprising that data for Tenontosaurus allowed for the identification of a best fitting growth model within the comprehensive class of BP-models (1) with relatively low variability in the parameters (i.e. small region of the near-optimal exponent-pairs). However, data uncertainty did not allow to conclude that Tenontosaurus would need a different exponent-pair (model) than modern alligators or birds. On the contrary, depending on the rescaling, displaying the data in dimensionless coordinates did not always show notable differences. In order to display differences, we rescaled mass relative to the asymptotic limit (of the best fitting model) and time relative to the age at the inflection point (age of maximal growth). Using this rescaling, we obtained the expected results: Modern broiler chicken grew much faster than dinosaurs or alligators and dinosaurs grew faster than alligators.

Supporting information

S1 File [xlsx]
Age (years) and mass (kg) data of , retrieved from [] ().

S2 File [xlsx]
Computation of (, ), based on , for certain grid-points, namely exponents and , and for them the best fit-parameters (optimization results) initial mass , , , and .

S3 File [xlsx]
Age (days) and mass (g) data of male from [].

S4 File [xlsx]
Computation of (, ), based on , for certain grid-points, namely exponents and , and for them the best fit-parameters (optimization results) initial mass , , , and .

S5 File [xlsx]
Age (years) and mass (kg) data of from [].

S6 File [xlsx]
Computation of (, ), based on , for certain grid-points, namely exponents and , and for them the best fit-parameters (optimization results) initial mass , , , and .


Zdroje

1. Kahm M, Hasenbrink G, Lichtenberg-Frate H, Ludwig J, and Kschischo M. (2010) Fitting Biological Growth Curves with R. Journal of Statistical Software 33: 1–21.

2. Aggrey SE. (2002) Comparison of Three Nonlinear and Spline Regression Models for Describing Chicken Growth Curves. Poultry Sciences 81: 1782–1788.

3. Renner-Martin K, Kühleitner M, Brunner N, Hagmüller W. (2016) AIC-Based Selection of Growth Models: The Case of Piglets from Organic Farming. Open Journal of Modelling and Simulation 4: 17–23. doi: 10.4236/ojmsi.2016.42002

4. Renner-Martin K, Brunner N, Kühleitner M, Nowak WG, and Scheicher K. (2018) Optimal and near-optimal exponent-pairs for the Bertalanffy-Pütter growth model. PeerJ 6: e5973. doi: 10.7717/peerj.5973 30505634

5. Lee AH, Huttenlocker AK, Padian K, and Woodward HN. (2013) Analysis of Growth Rates. In Padian K. and Lamm E.T. (eds.) Bone histology of Fossil Tetrapods: Advancing Methods, Analysis and Interpretation. UCLA Press: Berkeley, USA, pp. 217–251.

6. Erickson GM. (2014) On Dinosaur Growth. Annual Review of Earth and Planetary Sciences 42: 675–697.

7. West GB, Brown JH, Enquist BJ. (2001) A general model for ontogenetic growth. Nature 413: 628–631. doi: 10.1038/35098076 11675785

8. Pütter A. (1920) Studien über physiologische Ähnlichkeit. VI. Wachstumsähnlichkeiten. Pflügers Archiv für die Gesamte Physiologie des Menschen und der Tiere 180: 298–340.

9. Brody S. (1945) Bioenergetics and growth. Hafner Publ. Comp.: New York, NY, USA.

10. von Bertalanffy L. (1949) Problems of organic growth. Nature 163: 156–158. doi: 10.1038/163156a0 18107581

11. Verhulst PF. (1838) Notice sur la loi que la population suit dans son accroissement, CorrespondenceMathematique et Physique (Ghent) 10: 113–121.

12. Ohnishi S, Yamakawa T, Akamine T. (2014) On the analytical solution for the Pütter-Bertalanffy growth equation. Journal of Theoretical Biology 343: 174–177. doi: 10.1016/j.jtbi.2013.10.017 24211257

13. von Bertalanffy L. (1957) Quantitative laws in metabolism and growth. Quarterly Reviews of Biology 32: 217–231.

14. Richards FJ. (1959) A Flexible Growth Function for Empirical Use. Journal of Experimental Botany 10: 290–300.

15. Pauly D. (1981) The relationship between gill surface area and growth performance in fish: a generalization of von Bertalanffy’s theory of growth. Reports on Marine Research (Berichte der deutschen wissenschaftlichen Kommission für Meeresforschung) 28:25–282.

16. Gompertz B. (1832) On the Nature of the Function Expressive of the Law of Human Mortality, and on a New Mode of Determining the Value of Life Contingencies. Philos. Trans. R. Soc. London 123: 513–585.

17. Marusic M, Bajzer Z. (1993) Generalized two-parameter equations of growth. Journal of Mathematical Analysis and Applications 179: 446–462.

18. Lee AH, and Werning S. (2008) Sexual maturity in growing dinosaurs does not fit reptilian growth models. Proceedings of the National Academy of Sciences USA 105: 582–587.

19. Rootes WL, Chabreck RH, Wright VL, Brown BW, and Hess TJ. (1991) Growth Rates of American Alligators in Estuarine and Palustrine Wetlands in Louisiana. Estuaries 14: 489–494.

20. Brunner N, Kühleitner M, Nowak WG, Renner-Martin K, Scheicher K. (2019) Growth Patterns of birds, dinosaurs and reptiles: Are differences real or apparent? Preprint (at BioRxiv). Published online: doi: 10.1101/597260v1

21. Renner-Martin K, Brunner N, Kühleitner M, Nowak WG, Scheicher K. (2019) Best-fitting growth curves of the von Bertalanffy-Pütter type. Poultry Science 98 (9): 3587–3592, doi: 10.3382/ps/pez122 30895317

22. Calder WA. III. (1985) Size, Function, and Life History. Harvard Univ. Press: Cambridge, USA.

23. Pauly D, and Cheung WWL. (2017) Sound physiological knowledge and principles in modeling shrinking of fishes under climate change. Global Change Biology. Published online: doi: 10.1111/gcb.13831 28833977

24. Viboud C, Simonsen L, Chowell G.(2016) A generalized-growth model to characterize the early ascending phase of infectious disease outbreaks. Epidemics 15: 27–37. doi: 10.1016/j.epidem.2016.01.002 27266847


Článek vyšel v časopise

PLOS One


2019 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#