A novel approach to the analysis of human growth
© Gliozzi et al.; licensee BioMed Central Ltd. 2012
Received: 2 November 2011
Accepted: 17 May 2012
Published: 17 May 2012
Skip to main content
© Gliozzi et al.; licensee BioMed Central Ltd. 2012
Received: 2 November 2011
Accepted: 17 May 2012
Published: 17 May 2012
Several formulations have been proposed in order to model human growth from birth to maturity. They are usually based on “ad hoc” heuristic assumptions. In the present contribution we adopt, as an alternative, a completely general (interdisciplinary) approach, based on the formalism of the Phenomenological Universalities (PUN).
The main PUN class investigated to date, i.e. UN, can only account for the overall growth pattern. For a realistic description it is necessary to add to it one or more “spurts”, as expected on biological grounds, due to the stimulation of growth and sex hormones.
A new PUN class (UN + FM) is generated and shown to be able to provide excellent agreement with standard auxological datasets. The accuracy of the fitting and reliability of the model suggest applications both at the diagnostic and therapeutic level.
The developed formalism can be suitably related to the biological description of bone plate growth under selective hormonal stimulation on the bone epiphysis; i.e., the additional increase of stature is the “macroscopic” response to a well defined biological signal.
Growth is an extremely complex and non-linear biological process, driven by hormonal mechanisms and characterized by an intrinsic variability reflecting environmental and genetic influences and individual adaptive responses.
Individual human growth has been widely investigated [1–3]. In his brief history of human growth dynamics, Wales  reports that periodic accelerations and decelerations in growth had already been observed in 1777 by Montbeillard in his child, while seasonal and diurnal changes were detected by various authors during the 1960s and in later years . More recently, growth variability has been analyzed in detail leading to the concept of “saltatory growth” . The occurrence of “spurts” or discontinuities in longitudinal growth is now known to occur over small time intervals, and may therefore be detected only by using appropriate time scales.
When observations are made over infrequently collected data in “longitudinal” studies, or when growth charts are produced as “transversal” datasets, combining large numbers of data from a given homogeneous population (e.g. in terms of geographical area and of temporal range), growth variability is largely missed. Nevertheless the median curve from a growth chart usually exhibits spurts at mid-childhood (around 6–7 yr in both sexes) and at puberty (around 12 yr for females and 14 yr for males). Growth charts are very helpful as reference tools for pediatricians, in order to monitor individual growth and provide therapeutic interventions when children’s height is well below median values.
Similarly, animal growth has been shown to follow irregular “saltatory” patterns , which again are barely observable in studies based on infrequently collected data over large populations. In fact, at a low resolution time-scale, most animal species exhibit a very regular mass (m, or weight) growth pattern, characterized by an initial exponential growth and a progressive decrease of the growth rate a(t), until a final saturation value, the carrying capacity M, is reached. Such “universal” behaviour has been analyzed by West et al. [8, 9], who showed that by renormalizing the time by the effective growth rate and considering as dependent variable the ratio (m/M)0.25, most species share the same growth curve. West and collaborators have also investigated the dependence of the exponent value on the fractal nature of tissue perfusion, and the relation between the carrying capacity and some basic biological parameters, such as the cellular mass and metabolic rate of the species. For related literature, see [10–13].
For what concern the biological mechanisms underlying human longitudinal growth, a detailed explanation is still missing, and will probably entail multifactorial contributions to the main effect of selective hormonal stimulation on bone epiphyses. Basically, bone epiphyses are the regions where chondrocyte proliferation in the growth plate occurs; these cells can be seen as “stem-like cells”, whose growth is stimulated, at some specific point, by specific hormones . The most important of these hormones, i.e. the growth hormone GH, is secreted by the adenohypophysis and is controlled by the hypothalamus by means of various GH-stimulating mediators, e.g. GHRH (Growth Hormone Releasing Hormone) and GHR-IH (Growth Hormone Release-Inhibiting Hormone). Sleep, physical training, traumata and nutritional conditions are known to stimulate secretion of GH.
The GH effect on bone growth is known to be mediated by specific receptors, such as the IGF-I (Insulin-like growth factor) receptor . Sex hormones (SH, mainly estrogens, which in both boys and girls stimulate GH and IGF-I secretion) are also known to affect the endochondrial bones. Initially they induce a rapid growth in cartilage, producing the growth peak, then cartilages are calcified and height growth is blocked . Besides constitutional primitive and/or familial causes, any impairment in the production/effectiveness of such hormones may be responsible for pathological short stature [17, 18] (and can be counteracted with hormonal therapies, see Ranke et al. ). It would therefore be valuable to obtain direct indications relating the growth “spurts” and the biological effects of hormones in order to derive therapeutic guidelines. Such hormonal influence on the basic mechanisms of growth is fully in agreement with West’s assumption that the saturation value (or carrying capacity) is related to the cellular metabolic rate.
With regard to human growth description, many models have been proposed to date, mainly based on heuristic assumptions. A number of mathematical functions have been proposed, starting from the linear-logarithmic model [20, 21] and the first exponent function model . The logistic model was subsequently applied to human growth by Ozaki  and by Nelder . The Gompertz function  was used to describe the height of 24 male and 24 female subjects by Deming . In order to obtain a more realistic description, Thissen et al.  derived the double logistic function and, for a higher fitting accuracy, Bock and Thissen  proposed the triple logistic function. Similarly, the Jolicoeur model uses three terms , while Preece and Baines only use two , with a total of 5 free parameters. For a review, see [31, 32].
The goal of the present paper is to propose an alternative method for the fitting and modelling of human growth charts and individual (longitudinal) growth curves, based on the formalism of the Phenomenological Universalities (PUN). The PUN approach was inspired by casual observation of striking similarities between growth curves (and other phenomenologies) in different and independent fields and contexts, such as in Elastodynamics, Physics, Biology and Economics. It is explained in some detail in the next Section. In the succeeding section the problem of the exogenous vs. endogenous nature of the spurts is discussed and a specific formalism for their treatment is proposed. In the final section, several datasets of great auxological relevance are analysed in order to assess the applicability of the proposed approach and the validity and significance of the results.
Given a set of experimental data, the first thing that one may want to do, in order to extract from it as much information as possible, is to find a suitable fitting function. As a next step it may be desirable to construct a model out of it. For this purpose one should restrict one’s attention to the raw data and analyse them independently of the field of application. Such an unbiased procedure may be provided by the Phenomenological Universalities (PUN) approach, recently proposed by P.P. Delsanto and collaborators [33, 34] and applied to a wide range of topics (auxology , tumor growth [36, 37], nonlinear elasticity , and others [39–41]).
E.g. for N = 2 we have, in addition to the linear term βa, also the quadratic term γa2.
where a0 = a(t = 0) is the initial growth rate. Here and in the following we normalize the variable y(t), so that y(0) = 1.
Equation (5) represents the Gompertz law , which has been extensively used in all kinds of growth problems for almost two centuries. The parameter β represents a retardation factor, which may be defined as , where is the asymptotic value of z(t).
The formalism of the previous section is perfectly suitable to describe the evolution of an endogenous system, i.e. one in which all the “seeds” of growth are already present at the initial time (t = 0). There are, however, situations in which the “rules of the game” change, suddenly or slowly in a way that was not implicitly foreseeable at the time t = 0, possibly as a consequence of “forces” external to the system: we shall call these processes “exogenous”.
It is easy to find examples of exogenous processes in all scientific fields: e.g., a severe sickness or accident (such as the loss of a limb) during childhood, or the inception of angiogenesis in oncology, or a change of state in a thermodynamic process, or a “September 11” for the NYSE or OTC, etc. If the disruption of the system takes place at a well defined time t1, then the PUN formalism may be applied to both periods (before and after t1) and the change in the PUN parameters may be used to characterize quantitatively the variation which has occurred. If the change is sudden but its precise timing is not known, then t1 may be used as an adjustable parameter, in order to pinpoint the precise time of the perturbing event.
However, the distinction between exogenous and endogenous processes may be rather subtle. E.g., in the case of human development we can very well theorize that all the “seeds” are present in the DNA of the embryonic cells, so that the whole process is, at least implicitly, endogenous and a unique PUN curve should describe the entirety of perinatal growth (i.e. both before and after birth). In fact this turns out to be correct, at least as a first approximation . Likewise for the growth spurts mentioned in the Background, which occur both at mid-childhood and at a pre-puberal age due to hormonal stimulation.
A more accurate description of the human growth process can, however, be achieved by an alternative treatment, in which we superimpose on an overall U1 or U2 model one or more “mathematical spurts”, modelled as small additional corrections that are also to be described in the framework of a UN class. The resulting formalism is briefly described in the following and will be adopted in the remainder of this paper.
Let us then assume that the overall growth curve can be fitted, up to a certain level of approximation, by Eq. (5). We have restricted our treatment to the class U1 in order to limit the number of “free” parameters of the model, since the corresponding results are already quite satisfactory (see next section); the improvements achievable with U2 or U3 would be only marginal. It should be noted that, although in Eq. (5) there are three parameters, only one (β) should be considered as a “free” parameter, since y0 and a0 may be directly evaluated from the first few data.
The corresponding PUN classes will be called (U 1 + GM) and (U 1 + HM), respectively, with M = 1, 2, … being the number of spurts assumed. Obviously, since we are dealing with a steadily evolving process, results will be better with the former class (U 1 + GM), although the latter may be preferable when it is important to decrease the number of free parameters. A good compromise may be to adopt a class, which we shall call (U 1 + FM), which is identical to (U 1 + GM) except that the a0m parameters for the spurts are assumed to be equal to the U1 parameter a0. Such an assumption may easily be justified on biological grounds, since bone cells have their own duplication time, which is expected to remain the same both under ‘endogenous’ and ‘exogenous’ (i.e. hormonally induced) growth conditions.
Fitting parameter and R2 values corresponding to the longitudinal curve shown in Figure 1a (code IVS3)
the “historical” growth charts published by  and based on the London County Council survey updated to 1965 [45, 46]. In this survey a random sample of London schoolchildren was taken, including approximately 1,000 boys at each year of age;
the growth charts of Romano et al. , referring to Noonan boys and girls, based on the results of the National Cooperative Growth Study (which includes a large number of Noonan children). Noonan syndrome is a genetically-based pathology whose incidence in the general population is estimated between 1:1,000 and 1:2,500 . It is a clinically heterogeneous disorder predominantly characterized by dysmorphic facial features, congenital heart disease, proportionate post-natally short stature and other deformities .
Such data are normally used to define a “median” growth curve, which is the one considered in the following, and the percentiles, which are useful to pediatricians in order to establish threshold levels for underdevelopment.
U1 + H1
U1 + F1
t m (yr)
t m (yr)
t m (yr)
t m (yr)
Finally it can be observed that β1 is larger (in absolute value) for boys (−0.59) than for girls (−0.38) in the general population, while β1 for Noonan girls (−0.93) is larger than for Noonan boys (−0.43). Something similar happens with σ.
All the datasets analysed up to here are transversal, as stated. However, in order to evaluate the effect of Growth Hormone (GH) administration therapy to selected cases of critically short stature in infancy, “longitudinal” curves should be investigated. As an example, the dataset referring to a single young male, denoted by the code IVS3 in ref. , is analysed in Figure 1: it shows the two spurts, the first corresponding to the therapeutic administration of GH and the second to prepuberal growth (Table 1 reports the corresponding parameters and R2 values). In this case, as well as in similar ones, idiopathic short stature is corrected by the therapeutic administration of GH, which induces an additional ‘artificial’ spurt. A (U1 + F2) model is therefore required to account for the two spurts.
Growth charts represent a very useful reference tool for pediatricians to monitor the growth of individual children (i.e. size, rate of growth and effects of an eventual treatment). They are transversal in the sense that they yield avarage values (divided into percentiles) over large, well defined populations, e.g. boys (or girls) in a given country and time frame. Several, mostly heuristic, formulations have been proposed in order to model the growth, i.e. to fit e.g. the median growth curve (out of the charts) from early childhood to maturity. In this contribution we have presented an alternative approach, based on the formalism of the Phenomenological Universalities (PUN) [33, 34], which represent a new, completely general and interdisciplinary methodology.
As a result, we have found that the main PUN class studied to date, i.e. U1 and U2 [36–38, 41], can only predict the overall human growth pattern. For a more realistic description it is necessary to add to it one or more “spurts”, as also suggested by other authors  and well justified on biological grounds. Consequently a new PUN class (UN + FM) has been developed and shown to provide excellent agreement with some standard auxological datasets. Its formalism can be suitably related to the biological description of bone plate growth under selective hormonal stimulation of the bone epiphysis. In other words, the additional increase of stature is the “macroscopic” response to a well defined biological signal.
The accuracy of the fitting (see Figures 1a, 2c, 3c, and 4) is very important not only “per se”, but also because it adds a real significance to the model parameters. E.g., by comparing Figures 2c and 32c, we not only find that boys grow much later but faster than girls (as is well known), but we are also able to quantify the difference: growth spurt inception time (13 yr vs. 9 yr) and duration (2 yr vs. 3.2 yr) and total accretion due to the spurt (almost the same: 17.7 cm vs. 18 cm).
Likewise we may quantify the differences between the growth patterns of healthy vs. growth impaired children (e.g. Noonan [50–52]). E.g., as we can see in Table 2, the time of pubertal growth onset is quite different in Noonan vs. control girls (about two years later for the latter), while it is approximately the same for males. Likewise for the pubertal growth duration. As a consequence the difference in “added” stature is relatively small for boys (about 2 cm), while it is dramatic for girls (about 8 cm!). This suggests (even if other effects should also be considered, such as the time of reference for measurements) that pubertal growth is much more severely affected by Noonan syndrome in females than in males.
To conclude, the formalism presented in this contribution may also be applied to the analysis of other datasets of auxological interest, e.g. referring to variables such as the Body Mass Index (BMI): see e.g. . It may also be used for the systemic analysis of two or more variables assembled as a complex or vectorial quantity . In this case the goal is not only to investigate their time evolution, but also the degree of correlation and mutual dependence. As an example it would be instructive to study how mass growth follows stature in terms of relative increment of the various body components occurring differently with age, sex and life style.
We wish to acknowledge Gianpaolo Perego for his help on the numerical work.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.