Skip to main content

Scaling, growth and cyclicity in biology: a new computational approach

Abstract

Background

The Phenomenological Universalities approach has been developed by P.P. Delsanto and collaborators during the past 2–3 years. It represents a new tool for the analysis of experimental datasets and cross-fertilization among different fields, from physics/engineering to medicine and social sciences. In fact, it allows similarities to be detected among datasets in totally different fields and acts upon them as a magnifying glass, enabling all the available information to be extracted in a simple way. In nonlinear problems it allows the nonscaling invariance to be retrieved by means of suitable redefined fractal-dimensioned variables.

Results

The main goal of the present contribution is to extend the applicability of the new approach to the study of problems of growth with cyclicity, which are of particular relevance in the fields of biology and medicine.

Conclusion

As an example of its implementation, the method is applied to the analysis of human growth curves. The excellent quality of the results (R2 = 0.988) demonstrates the usefulness and reliability of the approach.

Background

Scaling, growth and cyclicity are basic "properties" of all living organisms and of many other biological systems, such as tumors. The search for scaling laws and universal growth patterns has led G.B. West and collaborators to the discovery of remarkably elegant results, applicable to all living organisms [14], and extensible to, e.g., tumors [57]. Cyclicity seems to be an almost unavoidable consequence of the feedback of every active biosystem from its environment. In the context of the present contribution we wish to extend the applicability of the Phenomenological Universalities (PUN) approach [8, 9], which allows the scaling invariance lost in nonlinear problems to be recovered, to growth phenomena that also involve cyclicities. The latter are particularly relevant in biology and medicine. We wish to make clear from the beginning that only mathematical universalities, as provided e.g. by partial differential equations, represent "true" universalities. As such, in a "top-down" approach, they have been used for centuries. However, we are often challenged, as in the present context, by observational or experimental datasets, from which we wish to "infer" some (more or less) general "laws" using a "bottom-up" approach. PUNs represent a paradigm for performing perform such a task on themost general level.

In muchthe same way that integers are defined as the 'Inbegriff' of a group of objects, when their nature is completely disregarded, PUN's may be defined as the 'Inbegriff' of a given body of phenomenology when the field of application and the nature of the variables involved are completely disregarded. They have been developed [8, 9] as a new epistemological tool for discovering, directly from the experimental data, formal similarities in totally different contexts and fields ranging from physics to biology and social sciences. This PUN "classification" can be made conspicuous by means of a simple test based on plots in the plane (a,b), where a and b are variables defined in Eq.(1) and (2), respectively.

Model and methods

In the present context PUN's are formulated as a method f for solving the following problem: given a string of data y i (t i ) and assuming that they refer to a phenomenology, which can be reduced to a first order ODE

y ˙ ( t ) = a ( y , t ) y ( t ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaiaacqGGOaakcqWG0baDcqGGPaqkcqGH9aqpcqWGHbqycqGGOaakcqWG5bqEcqGHSaalcqWG0baDcqGGPaqkcqWG5bqEcqGGOaakcqWG0baDcqGGPaqkcqGHSaalaaa@3E1F@
(1)

we search for a solution y(t), based not on simple numerical fitting, but on a universal (i.e. absolutely general) framework. The problem, of course, can be generalized to higher order ODE's, PDE's and/or to vectorial, rather than scalar relations, but we prefer here to keep the formalism at its simplest level.

If Eq. (1) refers to growth of any given organism or biological object, y is the mass (or length, height, etc.) of the body and t is time. To solve such a problem, let us start by assuming that a is a function solely of z = ln y and that its derivative with respect to z may be expanded as a set of powers of a. It follows that

b = a ˙ = d a d z z ˙ = n = 1 . α n a n ( z ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOyaiMaeyypa0JafmyyaeMbaiaacqGH9aqpjuaGdaWcaaqaaiabdsgaKjabdggaHbqaaiabdsgaKjabdQha6baakiqbdQha6zaacaGaeyypa0ZaaabCaeqaleaacqWGUbGBcqGH9aqpcqaIXaqmaeaacqGHEisPa0GaeyyeIuoacaaMc8UccqaHXoqydaWgaaWcbaGaemOBa4gabeaakiabdggaHnaaCaaabeWcbaGaemOBa4gaaOGaeiikaGIaemOEaONaeiykaKIaeyOla4caaa@4C4A@
(2)

If a satisfactory fit of the experimental data is obtained by truncating the set at the N-th term (or power of a), then we state that the underlying phenomenology belongs to the Universality Class UN.

It can be easily shown [8] that the Universality Class U 1 (i.e. with N = 1) represents the well known 'Gompertz' law [10], which has been used for more than a century to study diverse growth phenomena. The class U 2 includes, besides Gompertz as a special case, most of the commonly used growth models proposed to date in several fields of research, i.e., besides the already mentioned model of West and collaborators [1, 4], the exponential, logistic, theta-logistic, potential, von Bertalanffy, etc. models (for a review see Ref. [11]).

Restricting our attention to the class U 2, by solving the differential equations z ˙ = a MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOEaONbaiaacqGH9aqpcqWGHbqyaaa@2FAC@ and a ˙ = b MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyyaeMbaiaacqGH9aqpcqWGIbGyaaa@2F7D@ , and writing b, for brevity as

and with the (normalized) initial conditions y(0) = 1 and a(0) = 1,

a = β γ [ 1 ( 1 + β γ ) e β t ] 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaeyypa0JaeyOeI0scfa4aaSaaaeaacqaHYoGyaeaacqaHZoWzaaGcdaWadaqaaiabigdaXiabgkHiTmaabmaabaGaeGymaeJaey4kaSscfa4aaSaaaeaacqaHYoGyaeaacqaHZoWzaaaakiaawIcacaGLPaaacqWGLbqzdaahaaqabSqaaiabgkHiTiabek7aIjabdsha0baaaiaawUfacaGLDbaakmaaCaaaleqabaGaeyOeI0IaeGymaedaaaaa@45FE@
(4)

and

y = e z = { 1 + γ β [ 1 e x p ( β t ) ] } ( 1 / γ ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyEaKNaeyypa0Jaemyzau2aaWbaaeqaleaacqWG6bGEaaGccqGH9aqpdaqfGaqabeqaleaacqGGOaakcqGHsislcqaIXaqmcqGHVaWlcqaHZoWzcqGGPaqkaOqaamaacmaabaGaeGymaeJaey4kaSscfa4aaSaaaeaacqaHZoWzaeaacqaHYoGyaaGcdaWadaqaaiabigdaXiabgkHiTiabcwgaLjabcIha4jabcchaWjabcIcaOiabek7aIjabdsha0jabcMcaPaGaay5waiaaw2faaaGaay5Eaiaaw2haaaaaaaa@4DDC@
(5)

It is interesting to observe that Eq.(5) can be written as

which shows that the scaling invariance, which was lost due to the nonlinearity of a(z), may be recovered if the fractal-dimensioned variable u = y-γand the new variable τ = exp(βt) are considered. In fact, γ is, in general, non integer. In Eq. (6), c1 and c2 are constants: c 2 = γ β MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4yam2aaSbaaSqaaiabikdaYaqabaGccqGH9aqpcqGHsisljuaGdaWcbaqaaiabeo7aNbqaaiabek7aIbaaaaa@3427@ , c1 = 1 - c2.

It may also be useful to note that y is the solution of the Ordinary Differential Equation (ODE)

y ˙ = γ 1 y p γ 2 y , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaiaacqGH9aqpcqaHZoWzdaWgaaWcbaGaeGymaedabeaakiabdMha5naaCaaabeWcbaGaemiCaahaaOGaeyOeI0Iaeq4SdC2aaSbaaSqaaiabikdaYaqabaGccqWG5bqEcqGHSaalaaa@3AB0@
(7)

where p = 1 + γ ; γ1 and γ2 are constants: γ2 = β/γ and γ1 = 1 - γ2. Their sum is equal to 1, because of the normalization chosen (y(0) = 1). Equation (7) coincides with West's universal growth equation, except that here p may be totally general, while West and collaborators adopt Kleiber's prescription (p = 3/4) [12], which seems to be well supported by animal growth data. For other systems, different choices of p may be preferable: in particular C. Guiot et al. suggest a dynamical evolution of p in the transition from an avascular phase to an angiogenetic stage in tumors [13].

Equation (7) may have (as in Refs. [1] to [6]) a very simple energy balance interpretation, with γ1yprepresenting the input energy (through a fractal branched network), γ2y the metabolism and y ˙ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaiaaaaa@2D5A@ the asymptotically vanishing growth. In fact all UN (at least up to N = 3) fulfil energy conservation (or, equivalently, follow the first Principle of Thermodynamics). However, in U 1 there is no fractal dimensionality and both input energy and metabolism are proportional to y. In U 2, as we have seen, the energy input term has a fractal dimensionality. In U 3, one more term with a fractal exponent, again equal to p, contributes to growth. In fact, it can be shown (proof omitted here for brevity) that the U 3 ODE can be written as:

y ˙ + δ 1 d d t ( y p ) = δ 2 y p δ 3 y , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaiaacqGHRaWkcqaH0oazdaWgaaWcbaGaeGymaedabeaajuaGdaWcaaqaaiabdsgaKbqaaiabdsgaKjabdsha0baakiabcIcaOiabdMha5naaCaaabeWcbaGaemiCaahaaOGaeiykaKIaeyypa0JaeqiTdq2aaSbaaSqaaiabikdaYaqabaGccqWG5bqEdaahaaqabSqaaiabdchaWbaakiabgkHiTiabes7aKnaaBaaaleaacqaIZaWmaeqaaOGaemyEaKNaeyilaWcaaa@47DB@
(8)

where δ1, δ2 and δ3 are constants related to the coefficients β, γ, δ, of the truncated U3 series expansion of b

Let us now consider the case

a ( z , t ) = a ¯ ( z ) + a ˜ ( t ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaeiikaGIaemOEaONaeyilaWIaemiDaqNaeiykaKIaeyypa0JafmyyaeMbaebacqGGOaakcqWG6bGEcqGGPaqkcqGHRaWkcuWGHbqygaacaiabcIcaOiabdsha0jabcMcaPaaa@3DE9@
(10)

in which a is assumed to be the sum of two contributions to the growth rate, one ( a ¯ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyyaeMbaebaaaa@2D39@ ), that depends only on z (or y), while the other ( a ˜ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyyaeMbaGaaaaa@2D30@ ) is solely time-dependent. Then, by writing,

y = y ¯ ( t ) y ˜ ( t ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyEaKNaeyypa0JafmyEaKNbaebacqGGOaakcqWG0baDcqGGPaqkcuWG5bqEgaacaiabcIcaOiabdsha0jabcMcaPiabgYcaSaaa@38EB@
(11)

it follows

y ˙ = ( a ¯ + a ˜ ) y ¯ y ˜ , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaiaacqGH9aqpdaqadaqaaiqbdggaHzaaraGaey4kaSIafmyyaeMbaGaaaiaawIcacaGLPaaacuWG5bqEgaqeaiqbdMha5zaaiaGaeyilaWcaaa@37D6@
(12)

Eq. (1) can consequently be split into a system of two uncoupled equations

y ¯ ˙ = a ¯ ( z ) y ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaeHbaiaacqGH9aqpcuWGHbqygaqeaiabcIcaOiabdQha6jabcMcaPiqbdMha5zaaraaaaa@34E9@
(13)

and

y ˜ ˙ = a ˜ ( t ) y ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaGGbaiaacqGH9aqpcuWGHbqygaacaiabcIcaOiabdsha0jabcMcaPiqbdMha5zaaiaaaaa@34C2@
(14)

Eq. (13) can be solved as before (for the case a(z)) giving rise to the classes UN. A general solution of Eq. (14) can be written as

a ˜ = z ˜ ˙ = n = 1 . A n E n MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyyaeMbaGaacqGH9aqpcuWG6bGEgaacgaGaaiabg2da9maaqahabeWcbaGaemOBa4Maeyypa0JaeGymaedabaGaeyOhIukaniabggHiLdGaaGPaVRGaemyqae0aaSbaaSqaaiabd6gaUbqabaGccqWGfbqrdaWgaaWcbaGaemOBa4gabeaaaaa@3F0B@
(15)

where E n = exp[i(nωt + Ψ n )]. Then, if the sum in Eq.(15) can be truncated to the M-th term, we will state that the corresponding phenomenology belongs to the class UN/TM. It may be interesting to remark that the class U 0/TM and its phenomenology, involving the appearance of hysteretic loops and other effects, has been analyzed in detail, in a completely different context (Slow and Fast Dynamics [9]), under the name of Nonclassical Nonlinearity.

From Eq.(15) it is easy to obtain z ˜ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOEaONbaGaaaaa@2D62@ , y ˜ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaGaaaaa@2D60@ and b ˜ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOyaiMbaGaaaaa@2D32@ by simple integration or differentiation, but it is no longer possible, from a simple fitting of the experimental curve b(a), to obtain the relevant UN or TM parameters analytically, keeping, of course, only the real part of Eq.(15). However, assuming that N = 0, i.e. a = a ˜ ( t ) MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaeyypa0JafmyyaeMbaGaacqGGOaakcqWG0baDcqGGPaqkaaa@32A4@ , (case U0/T1) we can easily see that the curve b vs. a becomes

a 2 + ( b ω ) 2 = A 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyae2aaWbaaeqaleaacqaIYaGmaaGccqGHRaWkdaqfGaqabeqaleaacqaIYaGmaOqaamaabmaajuaGbaWaaSaaaeaacqWGIbGyaeaacqaHjpWDaaaakiaawIcacaGLPaaaaaGaeyypa0Jaemyqae0aaWbaaeqaleaacqaIYaGmaaGccqGHSaalaaa@3A27@
(16)

i.e. it represents an ellipse (see Fig. 1a), with ω being the ratio between the two semi-axes. In the case U 0/T 2, the "interference" between the ellipses generated by the first and second harmonics gives rise to plots, which include two complete ellipses (Fig.1b) or, according to whether A2 <<A1 or A2 > A1, one complete and one collapsed in a cusp (Fig. 1c) or in a knot (Fig. 1d), respectively. The plot b vs. a also depends, of course, also on the phase shift between the two harmonics and more complex curves may result (Fig. 1e) if its value is not close to 0 or to π. For N > 2, the plots obviously become more complex, nevertheless they may often be relatively easy to decipher, as in the U 0/T 3 case shown in Fig. 1f.

Figure 1
figure 1

Examples of b ( a ) curves belonging to the class U 0/ TM. Examples of b(a) curves belonging to the class U 0/TM, as described by Eq. (15). We have assumed for all the plots ω = 2, A1 = 1 and Ψ1 = 0. As predicted by Eq. (16), in the case M = 1 we obtain an ellipse with a ratio ω between the two semi-axes. Examples of the M = 2 case are shown in the plots (b), (c), (d) and (e), for different choices of the parameters A2 and Ψ2. As expected, three ellipses appear in the case U 0/T 3 (plot f).

When N ≠ 0 the additional problem of interference between a ¯ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyyaeMbaebaaaa@2D39@ and a ˜ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyyaeMbaGaaaaa@2D30@ and between b ¯ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOyaiMbaebaaaa@2D3B@ and b ˜ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOyaiMbaGaaaaa@2D32@ arises. However, in the case UN/T 1, if for brevity we write Φ1 = ω t + Ψ1, from

a = a ¯ + A 1 c o s Φ 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaeyypa0JafmyyaeMbaebacqGHRaWkcqWGbbqqdaWgaaWcbaGaeGymaedabeaakiabcogaJjabc+gaVjabcohaZjabfA6agnaaBaaaleaacqaIXaqmaeqaaaaa@39A2@
(17)

and

b = b ¯ A 1 ω s i n Φ 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOyaiMaeyypa0JafmOyaiMbaebacqGHsislcqWGbbqqdaWgaaWcbaGaeGymaedabeaakiabeM8a3jabcohaZjabcMgaPjabc6gaUjabfA6agnaaBaaaleaacqaIXaqmaeqaaaaa@3B88@
(18)

it follows that

( a a ¯ ) 2 + ( b b ¯ ) 2 ω 2 = A 1 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeiikaGIaemyyaeMaeyOeI0IafmyyaeMbaebacqGGPaqkdaahaaqabSqaaiabikdaYaaakiabgUcaRKqbaoaalaaabaGaeiikaGIaemOyaiMaeyOeI0IafmOyaiMbaebacqGGPaqkdaahaaqabeaacqaIYaGmaaaabaGaeqyYdC3aaWbaaeqabaGaeGOmaidaaaaakiabg2da9iabdgeabnaaDaaaleaacqaIXaqmaeaacqaIYaGmaaGccqGHSaalaaa@4277@
(19)

i.e. we still have an ellipse, whose centre, however, moves alongside the b ¯ ( a ¯ ) MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOyaiMbaebacqGGOaakcuWGHbqygaqeaiabcMcaPaaa@3050@ curve. As a consequence, a number (not necessarily integer) n e = ΔT/T of deformed ellipses is generated (ΔT is the time interval considered and T = 2π/ω). In the case UN/T 0, of course, only one ellipse is visible, since, in the plot b(a), the ellipse is retraced upon itself any number of times.

To illustrate the interference between an U2 curve (Fig. 2a) and the cyclicity contribution, we show in Fig. 2b and Fig. 2c the b(a) plots in the cases n e = 1 and 5, respectively. In spite of the ellipses' deformation, due to the curvature of the b ¯ ( a ¯ ) MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOyaiMbaebacqGGOaakcuWGHbqygaqeaiabcMcaPaaa@3050@ line, the approximate values of n e , ω and A1 can be retrieved and used as initial values for a fitting of y = y ˜ y ¯ MathType@MTEF@5@5@+=feaafiart1ev1aqatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyEaKNaeyypa0JafmyEaKNbaGaacuWG5bqEgaqeaaaa@3174@ , where

Figure 2
figure 2

Interference between U2 and the cyclical term. Interference between U2 and the cyclical term: (a) M = 0, i.e. no cyclical term; (b) M = 1 and n e = ΔT/T = 1, number of periods; (c) M = 1 and n e = 5

y ˜ = exp ( A 1 ω sin ω t ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaGaacqGH9aqpcyGGLbqzcqGG4baEcqGGWbaCdaqadaqaaKqbaoaalaaabaGaemyqae0aaSbaaeaacqaIXaqmaeqaaaqaaiabeM8a3baakiGbcohaZjabcMgaPjabc6gaUjabeM8a3jabdsha0bGaayjkaiaawMcaaiabcYcaSaaa@414D@
(20)

as it can be immediately obtained from Eq.(14).

In Eq. (20) it has been assumed Ψ1 = 0. Such an assumption is justified by the fact that cyclicity is usually due to an interaction of the system being considered with its environement (as a feedback from it), and t = 0 is chosen as the time at which the interaction starts.

Results

As an example of application of the proposed methodology, we consider in the following the curve of human weight development from birth to maturity. We refer to the classical work of Davenport [14] (nowadays an auxological standard), which suggests that the human growth rate exhibits three maxima: one intrauterine, a second one around the 6-th year and a third one other around the 16-th year. The last growth acceleration (adolescent spurt) seems to be activated by the secretions of the pituitary gland and/or the anterior lobe of the hypophysis, while no clear explanations have been proposed for the prenatal and the mid-childhood spurts.

Even if Davenport's finding are still actual, there has been a considerable debate over their interpretation. In fact for man, as for other social mammalians (e.g. elephants, lions, primates), growth development is greatly affected by cyclicity. It has been stated that it cannot be described by a simple curve, but that it requires at least two Gompertz or logistic-like curves (or three for humans [15]), describing the early growth and the juvenile phases separately. The period of extended juvenile growth is most marked in humans, for whom the total period of growth to mature size is very long in comparison with all other mammals.

In addition to the above main accelerations, many authors have observed short-term oscillations in longitudinal data. In the paper of Butler and McKie [16], 135 children were monitored at six monthly intervals from 2 to 18 years of age. Longitudinal studies reveal a cyclic, rhythmic pattern, as a sequence of spurts and lags occurring up to adolescence. In addition, in the paper of Wales [17], very short time cyclicities are reported, such as postural changes in height throughout the day, due to spinal disc compression. Variations in height velocity have also been described with the season of the year, possibly modulated through the higher central nervous system and secretion of melatonin and other hormones with circadian rhythmicity.

Our method allows us to fit Davenport's curve without the need to consider coupled logistic curves. In Figure 3 we show the original data, relative to human weight development from birth to maturity, and the fitting obtained with the curve U2/T1. The value of R2 = 0.998 confirms the correctness of the PUN classification and the accuracy and reliability of the approach. The presence of cyclicity is betrayed by the plot b(a) in Fig. 4, which clearly exhibits a loop (a very distorted ellipse). Since the curve of Fig. 3 was obtained from 'transversal', instead of 'longitudinal' data, it has been possible to detect only the overall "macroscopic" periodicity. In addition, by separately plotting the curves U2 and T1 vs. time, it is confirmed that the minima and maxima of the T1 curve fall at about 6 years and 17 years, i.e. where the complete U2/T1 curve has its inflection points (see Fig.5).

Figure 3
figure 3

Curve of human weight development from birth to maturity. Curve of human weight development from birth to maturity, based on Davenport data (dots) [14]. The solid line represents the fitting obtained with a U2/T1 curve (R2 = 0.998).

Figure 4
figure 4

Evidence of a cyclicity period in the human growth curve. Evidence of a cyclicity period in the human growth curve (Fig. 3) from the plot b(a).

Figure 5
figure 5

Separate plots of the curves U2 and T1. Separate plots of the curves U2 and T1 (see Eqs. (5) and (20)) for the case presented in Figs. 3 and 4. The minimum and maximum of the T1 curve coincide with the times for which the rate of growth is expected to have a minimum or a maximum, in correspondence with the inflection points in the y(t) curve.

Discussion and conclusion

After a short review of the Phenomenological Universalities (PUN) approach, we have proceeded to extend its range of applicability to problems of growth with cyclicity. We have analyzed in detail the case a ( z , t ) = a ¯ ( z ) + a ˜ ( t ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaeiikaGIaemOEaONaeyilaWIaemiDaqNaeiykaKIaeyypa0JafmyyaeMbaebacqGGOaakcqWG6bGEcqGGPaqkcqGHRaWkcuWGHbqygaacaiabcIcaOiabdsha0jabcMcaPaaa@3D9B@ , in which the growth rate is assumed to be separable in two terms, depending on z = ln y and t (time), respectively. y is the normalized mass (or height, length, etc.) of the body, the development of which is under analysis. As a result, we find that the UN classes, which have been defined and studied for problems without cyclicity, can be generalized as UN/TM classes, where TM represents the solution of the case with only the time dependent term a ˜ ( t ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyyaeMbaGaacqGGOaakcqWG0baDcqGGPaqkaaa@3052@ .

In the plots b vs. a (where a and b represent the first and second derivatives of z = ln y, respectively), the presence of cyclicity is betrayed by the appearance of "loops", which look like distorted ellipses, in a number that is equal to ΔT/T, where ΔT is the time range being considered and T is the cyclicity period. To be more specific, if we consider, e.g., the class U 2/T 1, we have two parameters (β and γ), which characterize the class U 2, and two more (A and ω), related to the cyclcity. From the appearance of the loops it is possible to obtain "initial" or "guess" values of A and ω, which allow us to fit the experimental or observational data using the U 2/T 1 general solution (Eqs. (11), (5) and (20)).

In order to demonstrate the reliability and accuracy of the method, the very important and yet not well understood problem of human growth has been considered. The classical transversal curve of Davenport [15] has been analyzed. The results (Figs. 3 and 5), with a value of R2 = 0.998 and the prediction of the acceleration spurts, demonstrate the validity of the approach. More information about human growth mechanisms may be obtained by analyzing longitudinal growth curves for individual or specific groups with the proposed methodology, thus leading to suggestions or evaluations of models incorporating suitable growth mechanisms. Many applications can, of course, be envisaged, such as the diagnosis of undernourishment or diseases, which affect the growth of an individual, or the comparative study of diverse growth patterns in different populations, or the correlation between mass and height development, etc.

An extension of the method presented in this paper to the case of coupled equations, or, more generally, vectorial relations (see e.g. [18]), is in progress.

References

  1. West GB, Brown JH, Enquist BJ: The fourth dimension of life: Fractal geometry and allometric scaling of organisms. Science. 1999, 284: 1677-1679. 10.1126/science.284.5420.1677.

    Article  CAS  PubMed  Google Scholar 

  2. West GB, Brown JH, Enquist BJ: A general model for otogenetic growth. Nature. 2001, 413: 628-631. 10.1038/35098076.

    Article  CAS  PubMed  Google Scholar 

  3. Gillooly JF, Charnov EL, West GB, Savage VM, Brown JH: Effects of Size and Temperature on Developmental Time. Nature. 2002, 417: 70-73. 10.1038/417070a.

    Article  CAS  PubMed  Google Scholar 

  4. West GB, Brown JH: Life's universal scaling laws. Physics Today. 2004, 57 (9): 36-43. 10.1063/1.1809090.

    Article  Google Scholar 

  5. Guiot C, DeGiorgis PG, Delsanto PP, Gabriele P, Deisboeck TS: Does tumor growth follow a 'universal law'?. J Theor Biol. 2003, 225: 147-283. 10.1016/S0022-5193(03)00221-2.

    Article  PubMed  Google Scholar 

  6. Delsanto PP, Guiot C, Degiorgis PG, Condat AC, Mansury Y, Desiboeck TS: Growth model for multicellular tumor spheroids. Appl Phys Lett. 2004, 85: 4225-4227. 10.1063/1.1812842.

    Article  CAS  Google Scholar 

  7. Delsanto PP, Griffa M, Condat CA, Delsanto S, Morra L: Bridging the gap between mesoscopic and macroscopic models: the case of multicellular tumor spheroids. Phys Rev Lett. 2005, 94: 148105-10.1103/PhysRevLett.94.148105.

    Article  CAS  PubMed  Google Scholar 

  8. Castorina P, Delsanto PP, Guiot C: Classification Scheme for Phenomenological Universalities in Growth Problems in Physics and Other Sciences. Phys Rev Lett. 2006, 96: 188701-10.1103/PhysRevLett.96.188701.

    Article  CAS  PubMed  Google Scholar 

  9. Delsanto PP, ed: Universality of Nonclassical Nonlinearity with applications to NDE and Ultrasonics Springer. 2007

    Google Scholar 

  10. Gompertz B: On the nature of the function expressive of the law of human mortality and on a new mode of determining life contingencies. Phil Trans Roy Soc Lond. 1825, 123: 513-10.1098/rstl.1825.0026.

    Article  Google Scholar 

  11. de Vladar HP: Density-dependence as a size-independent regulatory mechanism. J Theor biol. 2006, 238: 245-256.

    Article  PubMed  Google Scholar 

  12. Krieger RE, ed: The fire of life: an Introduction to Animal Energetics Huntington. 1975

    Google Scholar 

  13. Guiot C, Delsanto PP, Carpinteri A, Pugno N, Mansury Y, Deisboeck TS: The dynamic evolution of the power exponent in a universal growth model of tumors. J Theor Biol. 2006, 240: 459-10.1016/j.jtbi.2005.10.006.

    Article  PubMed  Google Scholar 

  14. Davenport CB: Human Growth Curve. J Gen Physiol. 1926, 10: 205-216. 10.1085/jgp.10.2.205.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Lozy M: A critical analysis of the double and triple logistic growth curves. Ann Hum Biol. 1978, 5: 389-94. 10.1080/03014467800003021.

    Article  PubMed  Google Scholar 

  16. Butler GE, McKie M, Ratcliffe SG: The cyclical nature of prepubertal growth. Ann Hum Biol. 1990, 17: 177-98. 10.1080/03014469000000952.

    Article  CAS  PubMed  Google Scholar 

  17. Wales JK: A brief history of the study of human growth dynamics. Ann Hum Biol. 1998, 25: 175-84. 10.1080/03014469800005552.

    Article  CAS  PubMed  Google Scholar 

  18. Banchio AJ, Condat CA: Seasonality and Harvesting, Revised. Universality of Nonclassical Nonlinearity with applications to NDE and Ultrasonics Springer. Edited by: Delsanto PP. 2007

    Google Scholar 

Download references

Acknowledgements

We wish to acknowledge the support of a Lagrange fellowship from the C.R.T. Foundation (for A.S.G.).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Antonio S Gliozzi.

Additional information

Authors' contributions

The first author PPD has developed the general formalism, the second ASG has developed the numerical tools and carried out the numerical analysis and the third CG has suggested and analyzed the applicative context of the paper.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Delsanto, P.P., Gliozzi, A.S. & Guiot, C. Scaling, growth and cyclicity in biology: a new computational approach. Theor Biol Med Model 5, 5 (2008). https://doi.org/10.1186/1742-4682-5-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-4682-5-5

Keywords