A universal power law for modelling the growth and form of teeth, claws, horns, thorns, beaks, and shells
BMC Biology volume 19, Article number: 58 (2021)
A major goal of evolutionary developmental biology is to discover general models and mechanisms that create the phenotypes of organisms. However, universal models of such fundamental growth and form are rare, presumably due to the limited number of physical laws and biological processes that influence growth. One such model is the logarithmic spiral, which has been purported to explain the growth of biological structures such as teeth, claws, horns, and beaks. However, the logarithmic spiral only describes the path of the structure through space, and cannot generate these shapes.
Here we show a new universal model based on a power law between the radius of the structure and its length, which generates a shape called a ‘power cone’. We describe the underlying ‘power cascade’ model that explains the extreme diversity of tooth shapes in vertebrates, including humans, mammoths, sabre-toothed cats, tyrannosaurs and giant megalodon sharks. This model can be used to predict the age of mammals with ever-growing teeth, including elephants and rodents. We view this as the third general model of tooth development, along with the patterning cascade model for cusp number and spacing, and the inhibitory cascade model that predicts relative tooth size. Beyond the dentition, this new model also describes the growth of claws, horns, antlers and beaks of vertebrates, as well as the fangs and shells of invertebrates, and thorns and prickles of plants.
The power cone is generated when the radial power growth rate is unequal to the length power growth rate. The power cascade model operates independently of the logarithmic spiral and is present throughout diverse biological systems. The power cascade provides a mechanistic basis for the generation of these pointed structures across the tree of life.
The discovery of general models and mechanisms that create the phenotypes of organisms is a major goal of evolutionary developmental biology [1,2,3,4,5]. Very few such fundamental growth patterns exist, including logarithmic spiral growth [6, 7]. These growth patterns are important because they significantly influence the diversity of life by making some phenotypes very common while constraining or even prohibiting others, essentially favouring specific evolutionary trajectories [1, 8,9,10,11,12,13].
The vertebrate dentition, with its panoply of morphological diversity, is a superb system in which to investigate models of growth. All vertebrate teeth grow from the tip downwards towards the base to form the main body and individual cusps of each tooth. Teeth are often described as being ‘conical’ [14, 15]—this term may refer to the strict mathematical shape of a cone (Fig. 1a) or perhaps a single-pointed structure that folds down on all sides. While we have made great strides in determining the genetic influences on cusp formation and variation [16, 17], currently, we do not know the main determinants of cusp shape.
Following the discovery of the logarithmic spiral , in 1659, Sir Christopher Wren  first proposed that shells grow as a cone expanding along a logarithmic spiral (Fig. 1a, b), so that the trajectory of the midline of the shell forms a logarithmic spiral. This approach has since been used to model shell growth [7, 20]. Thompson  concluded that teeth follow this conical pattern of growth; here, we will test the suggestion of Wren  and Thompson  that teeth grow in the shape of a cone. At the same time, we strive to ascertain whether there are high-level, simple models governing how teeth grow and to determine if these patterns extend beyond the dentition to other pointed structures in nature.
Power cascade model simulates tooth growth
We represent the shape of the tooth in the manner in which it grows, from tip to base, by measuring the rate of lateral expansion of the tooth as the length increases. To measure this rate for a tooth, we take a 3D digitised surface of the tooth and place 10 equally spaced cross-sections perpendicular to its midline (Fig. 2a). The average radius of each cross-section is Radius = √ (cross-sectional area/π). We then plot log10Distance from the tip against log10Radius. The example Tyrannosaurus rex tooth illustrated in Fig. 2a shows a very close relationship to a straight line (Fig. 2b), fitting a linear model with R2 = 0.997.
All other teeth measured show a strong linear or close to linear relationship between log(Distance) and log(Radius) (here termed Log Distance-Radius plots; Fig. 2b; R2 range 0.9654–0.9998; see Additional file 1: Tables S1-S5 for sample sizes). This appears to indicate a general model of growth for single-cusped teeth, including bony fish, sharks, amphibians, reptiles, non-avian dinosaurs, birds, and mammals (Additional file 1: Figure S1).
The linear relationship, explicitly written log10(Radius) = Slope × log10(Distance) + Intercept, allows the growth and shape of a tooth to be characterised by its Slope and Intercept. This relationship can be written as Radius = 10Intercept × DistanceSlope. Therefore, the relationship is a ‘power law’  (monomial) with varying exponent (Slope) and multiplier (10Intercept).
When Slope = 1, there is a linear relationship between the raw measurements of Distance and Radius. Revolving this straight line around the x-axis to form a surface of revolution generates a cone (Fig. 1a; Additional file 1: Figure S2a). This shape is the cone expected by Wren . When Slope = 0.5, the surface of revolution is a paraboloid (Fig. 1c; Additional file 1: Figure S2a). Values below 0.5 are increasingly blunt at the very tip. The shape of the unicuspid teeth measured above therefore matches the surface of revolution of a power function, which we call a ‘power cone’ and is the same as the ‘power series’ of airplane nose cone designs . We will use the term ‘cone’ to refer to only the straight-sided conventional cone, while ‘power cone’ is the more general shape with a profile of varying curvature. Because Slope is always less than 1 (range 0.25–0.95), teeth are therefore not cone-shaped and do not match the conical model of Wren . In order to describe the folding of the tooth shape cascading down from the tip of the tooth according to the power function, we term this model the ‘power cascade’. The power cascade defines a new family of shapes that vary in Slope and Intercept.
The Intercept represents a scaling factor for the width of the tooth, with higher values resulting in wider teeth for the same length (Additional file 1: Figure S2b). The right-hand end of the Log Distance-Radius curve represents the maximum length of the tooth (Additional file 1: Figure S2c). Since all teeth plotted in Fig. 2 have been measured at ten equally spaced points along the tooth, the first distance for each tooth (x minimum) is 1/10 of the maximum distance (x maximum), and so these lines are the same length along the log-scaled x-axis (the only exception being the mammoth tusk; Additional file 1: Figure S2c). Ten intervals appears to be minimally required to represent the shape of a tooth; however, more points, whether or not equally spaced, can be used (see Additional file 1: Figure S3 for alternative sampling intervals) and do not substantially affect the calculation of Slope or Intercept (see Supplementary Discussion).
Tooth crowns, rows, and cusps follow the power cascade
The junction between the enamel and dentine layers in a tooth initially forms as the interface between the epithelium and mesenchyme during development . In most vertebrates, the enamel or enameloid has an equal thickness over the tooth surface, and so the outer enamel surface is an adequate approximation of the initial shape of the developing tooth. In some groups such as hominids, the enamel is of uneven thickness over the tooth (modelled by ). However, both the enamel-dentine junction (EDJ) and the outer enamel surface of the human canine show a tight fit to the power cascade (R2 > 0.997), although with different slopes (Fig. 3c).
Many reptiles and marine mammals have rows of unicuspid teeth along each jaw, sometimes showing variation in tooth shape along the row . To investigate whether the power cascade sufficiently represents growth in all teeth within a row, we measured 11 teeth in the lower jaw of the theropod dinosaur Tyrannosaurus rex. All teeth are represented by the power function (0.992 < R2 < 0.999; Fig. 3a), with Slope between 0.47 and 0.54.
Most mammal teeth have multiple cusps on the surface. Each of these cusps is generated from the position of an enamel knot that forms a local maximum of the epithelium-mesenchyme interface during early tooth development [16, 23]. Leopard seal Hydrurga leptonyx postcanine teeth are trident-shaped, with a large central cusp and smaller anterior and posterior cusps (Fig. 3d). Measuring each of these cusps separately shows a linear Log Distance-Radius relationship for each cusp. The Slope and Intercept are similar for all five postcanine teeth, with the central cusp tending to have a higher Intercept but similar Slope (Fig. 3d). More rounded cusps also show power cascade tooth growth, as seen in the wave-like shapes of the crab-eater seal Lobodon carcinophaga postcanine teeth (Fig. 3b) and in the enamel-dentine junction of the individual cusps of human molars (Fig. 3c). The entire heterodont tooth row of the leopard seal shows strong linear patterns in the incisors, canines, and main postcanine cusps (R2 > 0.997; Fig. 3e).
Power cascade model can predict tooth length and age
The attributes of the power cascade model have important implications for the growth and characterisation of teeth. Teeth are often worn or broken at the tip through use in life or during preservation and fossilisation. If teeth generally show growth according to a power cascade model, we can use this to estimate how much of the tip of a tooth has been lost. By sequentially adding length to the tip and fitting a linear model, the distance corresponding to the original tooth length should be closest to linear and have the highest R2 value for the regression. To demonstrate this concept, we can artificially remove part of the tip of a leopard seal canine and then estimate how much has been lost (Additional file 1: Figure S4b,c), resulting in an overestimate of 5.8 to 8.0% of the original length of the tooth (or less than 5% overestimate if accounting for apparent minor changes in growth early in tooth formation; see Supplementary Information). This method has been used to reconstruct the missing tip of some teeth in Fig. 2 such as the pliosaur fossil tooth (Additional file 1: Figure S4d).
A new method for estimating the age of mammals with ever-growing teeth, such as elephants and rodents, can be derived from this ability to reconstruct tooth length. Through use, a substantial amount of an ever-growing tooth is worn away (e.g. tusk wear on trees and gnawing in rodents), and so the complete length of the tooth must first be estimated using the power cascade. Then, we can use the rate of tooth growth (microns per day, based on crown extension rate [26, 27]) to calculate the length of time taken to grow the tooth. The result is a minimum estimate of the age of the individual.
In African elephants Loxodonta africana , Radius increases with Distance along the tusk following the power cascade relationship. Intercept is higher in males than females (Fig. 4a), thereby giving the ability to determine sex from an isolated tusk according to the rate of increase in radius. As expected, Radius at the base of a tusk increases with Age in years according to the power cascade (Fig. 4b), therefore giving an estimated Age for a given tusk Radius. A similar power cascade pattern is found in the incisors of the Zaisan mole vole Ellobius tancrei  during juvenile growth (Fig. 4c) but at a much smaller scale. While incisor circumference or width have been used to estimate age previously based solely on regression [30, 31], the power cascade describes the underlying pattern, and allows for estimation of age from growth parameters alone. Additional work on the limitations and accuracy of predictions derived from this approach is currently underway.
Close inspection of the Log Distance-Radius curves in Fig. 2 shows some deviations from the expected linear pattern of the power cascade. This is most obvious in the acrodont snake teeth, where the tooth fuses to the bone rather than sitting in a bony socket. The initial shaft of the tooth follows the power cascade model, but close to its base, the radius increases faster than expected, deviating from linear. This is likely due to additional widening of the dental epithelium at the base of the tooth to aid the fusion of the tooth to the bone in these snakes. Deviation from the power cascade pattern may also occur as the developing tooth is curved to fit within the jaw before mineralisation has occurred.
Another apparent cause of deviation from the power cascade is the presence of grooves running down part of the length of the tooth, such as those found in snake fangs (Fig. 2) and felid canines (Additional file 1: Figure S1). The infolding of the grooves reduces the cross-sectional area and so appears to cause deviation from the expected linear pattern. Wear on the tip or lateral surface of the teeth will also cause deviation, and so slightly worn or broken teeth were only included where they could be confidently reconstructed.
Tooth shape is defined by power cascade morphospace of Slope and Aspect Ratio
Power cones can be considered self-similar curves in that all curves of the same power can be stretched or rescaled to be the same shape (using affine transformations of translation and scaling akin to self-affine fractals; Additional file 1: Figure S5). Therefore, the only distinguishing feature of teeth with the same Slope is the relative stretching of the curve, i.e. the Aspect Ratio, calculated as the maximum length divided by the maximum diameter (Additional file 1: Figure S2).
We can use these two key parameters of the power cascade to define dimensions of a morphospace that illustrates the range of combinations found in vertebrate teeth (Fig. 5). Squat shapes are at the bottom of the graph—these are mostly individual cusps, as well as a marine reptile unicuspid tooth. Elongated teeth (‘tusks’) occupy the top of the graph, where the region with the highest aspect ratio is occupied by the tusks of the narwhal Monodon monoceros, woolly mammoth Mammuthus primigenius, African elephant Loxodonta africana, and walrus Odobenus rosmarus. Most mammal canines have a Slope between 0.35 and 0.60 (with a mean standard error Slope for each tooth of 0.011). The highest Slope is found in the human canine EDJ. Snake fangs have a Slope between 0.4 and 0.6 and a high Aspect Ratio (above 2.5). All of the Tyrannosaurus rex teeth along the row fall within a small range of Aspect Ratio and Slope. Our examination of tooth shape not only includes extant organisms, but also the greater morphological diversity of extinct organisms, increasing our confidence in the generality of the model .
Power cascade growth is not accounted for by existing models of tooth growth
There are currently two general models that each strive to describe or explain various aspects of tooth development. Enamel knots produce inhibitory signals that prevent new enamel knots forming close to an existing knot . The ‘patterning cascade’ model describes how this inhibition, along with the folding of the epithelial-mesenchyme interface, creates limitations on the size and position of successive cusps during development . First described in seal postcanine teeth, the patterning cascade model has since been extended to primate molars [34, 35]. The second model, the ‘inhibitory cascade’, describes the relative size of sequentially produced teeth, such as molars, as a linear change in size along a tooth row [2, 36]. Neither of these models addresses the shape of cusps.
The power cascade model proposed here is a third general model of tooth development complementary to the two existing models, indicating how the shapes of unicuspid teeth and individual cusps are generated. After determination of cusp shape by the power cascade model, we postulate that cusp spacing is dictated by inhibition of enamel knots according to the patterning cascade , and number of cusps is controlled by the number of enamel knots that can fit in the total area of the tooth. The sizes of sequential teeth are then directed by the inhibitory cascade [2, 36]. Therefore, cusp shape, cusp number, and tooth size can be simulated according to this trio of models to generate the main features of an entire tooth row.
The power function has been used to represent or measure a limited set of teeth in previous studies, including the tips of shapes designed for mechanical penetration testing , using an average Slope of 0.5. Felid canine profiles measured using power functions  showed that they generally had a Slope of ~ 0.55. Both of these studies are consistent with the current findings in many mammal canines, but they did not generalise this pattern to all teeth or cusps.
Detailed developmental computer simulations of tooth morphogenesis have used a 3D reaction-diffusion-like model that calculates bending stresses to form cusps and teeth [39, 40]. This model produces cusp positions that can have morphological variation similar to biological teeth [39, 41]. Here we tested whether the cusp shapes produced by that model conform to the power cascade model. Varying five parameters of the model that simulates the development of ringed seal postcanine teeth  shows that most of the cusp shapes produced do not closely resemble the expected power cascade, with R2 between 0.59 and 0.97 (Additional file 1: Figure S6). Therefore, the power cascade model describes cusp shape (or cross-sectional profile) substantially better than complex in silico models, although this may be a result of the limited number of cells in the simulations.
Given the power of this new model to define the limits of tooth shape in animals, we expanded our focus to compare it with existing models of growth in other morphological systems. Wren’s  model of shells growing as a cone bending to form a logarithmic spiral has since been used to model shells and teeth [6, 7, 20]. Starting with a cone, a logarithmic spiral is generated when one side grows faster than the other, causing the cone to bend to one side (Fig. 1b; Additional file 1: Figure S7d). A mechanism to generate a logarithmic spiral is the unequal growth rates of the two sides A and B. Logarithmic spirals have a formula in polar coordinates S = a eb θ, where θ is the angle of rotation around the origin, S is the resulting radius of the logarithmic spiral, and a and b are parameters affecting the size and rate of expansion of the spiral, respectively (Additional file 1: Figure S7a). The radius of the shell opening expands linearly with the angle of rotation (Radius = c θ, where c is a parameter affecting the rate of growth of the shell opening), which creates a cone spiralling around the central axis (Additional file 1: Figure S7b). This model was used to generate shell shapes of many types by modifying relative rates of growth [7, 42, 43].
The Raup  shell equation describes shell growth using a cone, which is the shape where Slope = 1 in our Log Distance-Radius plots (Fig. 6; Additional file 1: Figure S2). If this model accurately describes shell growth, all shells should fall on the right-hand edge of the morphospace in Fig. 6. The shells of molluscs (scaphopod Dentalium sp. and gastropod Bembicium auratum) and cephalopods (nautilus Nautilus pompilius and ram’s head squid Spirula spirula) each apparently form logarithmic spirals, but follow the power cascade with Slope between 0.37 and 0.88 (Fig. 6). This shows that power cones can bend to form logarithmic spirals in an analogous manner to that first proposed by Wren  for cones (a specific power cone; Fig. 1d). It also establishes that not all shell shapes can be generated by the existing model of development . In order to accommodate such shapes, the Raup  model must have the Slope parameter added, such that Radius = c θSlope. In the first description of the shell growth model, Raup  assumes that ‘the rate of expansion of the generating curve is approximately constant’, i.e. Slope = 1, and so this parameter was not included in his model. In contrast, Thompson  suggested that the growth may not be constant in some shells but in fact vary ‘in accordance with some simple law’, and Ackerly  showed that for some shells there is an allometric component to the change in radius. Our power cascade model accounts for this important feature of growth.
The long axis of each tooth grows as a logarithmic spiral [6, 46], which can be seen in an extreme form in the curved upper tusks of the babirusa Babyrousa celebensis. However, we find that the Slope of these tusks (0.25) is considerably less than 1, and therefore, they are not conical (Fig. 5): their high Aspect Ratio can make them appear more conical. This means that teeth cannot be modelled by the Raup  shell equation. The radius of the circle must change logarithmically with the angle of rotation to form a power cone, rather than a straight-sided cone with Slope = 1.
A general model of growth for horns, claws, spines, beaks, and thorns
Thompson  expected that pointed and spiral structures such as horns and claws would follow the same growth pattern as shells, which has been used to model some horn-like structures . If horns grow according to the shell model and are spiralled cones, then their Slope parameter will be 1. From measurements of bony horn cores from vertebrates including mammals, non-avian dinosaurs (referred to here as dinosaurs) and reptiles, we have found that Log Distance-Radius plots are linear and the Slope is typically between 0.4 and 0.8 (Fig. 6; Additional file 1: Figure S8), demonstrating that they do follow the power cascade but are not growing according to the original conical shell model.
Other structures throughout vertebrates also show power cascade growth, including mammal, bird and dinosaur claw and hoof bones (unguals), the bony beaks of birds and dinosaurs, and spines of fish (Fig. 6). Outside vertebrates, the power cascade model is also followed in arthropod fangs and cephalopod beaks. Beyond animals, it is found in thorns and prickles in plants (Fig. 6).
The rose prickle (generally called a thorn) represents an interesting exception. While the concave shape of a mature prickle does not follow the power cascade prediction, a young prickle does (Additional file 1: Figure S9). It appears that the prickle is initially generated following the power cascade growth with Slope = 0.6, but then as the stem to which it is attached grows, the base of the prickle is stretched along the long axis of the branch. The result is the typical concave shape of a rose prickle, where only the top half follows the power cascade, not the basal half that has been stretched (Additional file 1: Figure S9). In general, it appears that deviations from the power cascade are more likely in pointed structures controlled by multiple growth processes.
The power cascade model can be added to the logarithmic spiral model to generate a ‘power spiral’ that can simulate realistic shapes of pointed, curved structures (Additional file 1: Figure S7c). Figure 7 shows some comparisons between real teeth and power spiral models, using both circular cross-sections that would be generated in surfaces of revolution and other cross-sectional shapes (elliptical, lenticular, truncated circle) implemented in a Mathematica notebook (v. 12.0, Wolfram Research Inc., Champaign, IL) available in the Supplementary Information (see also Additional file 1: Figure S10).
The majority of the structures that are closely emulated by the power cascade grow from tip to base, including teeth, horns, thorns, and prickles. These shapes are presumably formed as each addition of material increases the radius by a constant proportion for a proportional increase in length. For example, bovid horns grow from tip to base, increasing in radius down the horn, and they generally follow the power cascade model. In contrast, cervid antlers grow from base to tip, with the growing antler branching, and the antler points being the last structures to form. Despite this directional difference in growth—and the antler starting from a wider base and narrowing towards the tip—antler points also follow the power cascade (Fig. 6). This shows that the proportional growth pattern can act both when increasing the radius of the structure as it cascades downwards from the tip to the base, and also when decreasing the radius to cascade upwards from base to tip. It appears that only the direction of radial growth differs between these two scenarios.
Since many of the structures examined here (including teeth and claws) are used to penetrate food or other materials, it may be argued that selection to maximise penetration ability or structural strength is the cause of the underlying similarity in shape as described by the power cascade model. However, many structures that are not for penetration (such as shells, rounded teeth or backward-curving horns) still follow the power cascade pattern. Given that structures that conform to the power cone can vary from sharp and long to blunt and short, we argue that the most parsimonious explanation for the model fit is an underlying biophysical or developmental mechanism rather than strong selection for shapes that coincidentally fit a power cascade-like pattern. The power cascade generates a base set of allowed variations (Fig. 5), and selection chooses from among these shapes, as occurs with the selection of relative tooth size in hominins according to the inhibitory cascade .
Mechanism and generality of power cascade
The log-log linear pattern of the power cascade can be compared with allometric plots of the relative sizes of body components during growth , such as head size versus body size in humans. A linear allometric relationship is produced when two components grow exponentially at different rates. The power cascade relationship shows that there is an allometric relationship within the same structure due to differential growth rates of Radius and Distance.
We can demonstrate this growth process by examining power function growth in Distance and Radius over time (Fig. 8a): Distance ∝ TimerD and Radius ∝ TimerR, where rD and rR are the growth rates for Distance and Radius, respectively. Power function growth is very common in biology, including for human height  and elephant tusks (Fig. 4b). When both axes of the growth over time curves are logged, the plot log(Distance) vs log(Time) is linear with slope rD (similarly for Radius and rR; Fig. 8b). By solving the log(Distance) equation for log(Time) and substituting into the log(Radius) equation, the relationship between log(Distance) and log(Radius) through time becomes apparent (Fig. 8c). If rD and rR are equal, then Radius increases linearly with Distance (Fig. 8d) and produces a conical shape (with Log Distance-Radius power cascade Slope of 1). If instead the rates of growth of Distance and Radius differ (e.g. rD = 2rR), then the log-log growth over time trajectories will not be parallel (Fig. 8e-f), and the result will be a power cone such as a paraboloid (Fig. 8h). The Log Distance-Radius power cascade Slope of such a structure will be rR/rD = 0.5 (Fig. 8g).
Therefore, the power cascade is an expression of allometry as a shape: power cones show unequal power growth within the same structure, or ‘constant differential growth-ratios’ in the terminology of Huxley . The cone is produced through isometric growth between Distance and Radius, while a power cone results from allometric growth (rD ≠ rR). The same shapes can also be generated through exponential (as opposed to power) growth of body parts, although this is not commonly found in organisms. Constant differential growth of the two sides of a structure must generate a logarithmic spiral (; Fig. 1b). In the same manner, differential power growth of Distance and Radius must generate a power cone (Fig. 1c). Both mechanisms could operate at the same time, forming a power cone on a logarithmic spiral, or a power spiral (Fig. 1d).
The power cascade, and likewise the logarithmic spiral, can be seen as ‘dynamical patterning modules’  that generate patterns and structures in metazoans and plants. Despite over three centuries of research , the specific molecules driving logarithmic spiral growth are not known (although recent work has begun to reveal some components in gastropod shells ). Likewise, the identity of signalling molecules and genes that influence the differential growth of the power cascade very likely must vary widely across animals and plants. Here we show that common growth patterns in animals and plants generate power cones. These shapes may be considered the default family of shapes for pointed structures, meaning they are more likely to independently evolve multiple times and will be a likely source of homoplasy in evolution.
Since the time of Wren , the logarithmic spiral has been considered a fundamental pattern of biological growth, generated by differential growth rates . The power cascade directs the shape of an immense range of structures and likely is as widespread and elementary as the logarithmic spiral in nature. Due to the huge breadth of structures and taxa in which this pattern is found, it appears that the power cascade is a fundamental pattern of growth in myriad organisms [6, 51].
Materials and methods
Specimens and 3D scanning
Specimens scanned for this study were sourced from the following institutions or private collections: Museums Victoria, Melbourne, Australia (NMV); Monash University Zoology Research Collection, Melbourne, Australia (MZRC); Australian Museum, Sydney, Australia (AM); American Museum of Natural History, New York, NY, USA (AMNH); Evans EvoMorph Collection, Monash University, Melbourne, Australia (EEM); Natural History Museum of Los Angeles County, Los Angeles, CA, USA (LACM); Natural History Museum UK, London, UK (NHMUK); John Canning, Victoria, Australia (JC); Monash University Integrated Morphology and Palaeontology Laboratory, Melbourne, Australia (MU-IMP); South Australian Museum, Adelaide, Australia (SAMA); Judith Pollock, Melbourne, Australia (JFP); Tasmanian Museum and Art Gallery, Hobart, Australia (TMAG). 3D surface models of teeth were generated using microCT (Zeiss Xradia 520 Versa XRM, Monash University XMFIG X-ray Microscopy Facility for Imaging Geomaterials; Skyscan 1174, Monash University), medical CT (Siemens, Monash Biomedical Imaging), or 3D surface scanning (Laser Design Surveyor 2025; Artec Space Spider). Scans of specimens from the following museums were obtained from researchers or online databases such as MorphoSource and Aves3D: Alabama Museum of Natural History, Tuscaloosa, AL, USA (ALMNH); Ditsong National Museum of Natural History, Pretoria, South Africa (DNMNH); Kronosaurus Korner, Richmond, Queensland, Australia (KK); Museum für Naturkunde, Berlin, Germany (ZMB); Museum of Comparative Zoology, Harvard, MA, USA (MCZ); Natural History Museum, London, UK (NHMUK); Ohio University Vertebrate Collection, Athens, Ohio, USA (OUVC); Royal Ontario Museum, Ontario, Canada (ROM); Smithsonian Institution, Washington, D.C., USA (USNM); Stony Brook University, NY, USA (SBU); University of Leipzig Anatomical Collection, Leipzig, Germany (ULAC); University of Washington Burke Museum, Seattle, USA (UWBM); Yale Peabody Museum, Newhaven, CT, USA (YPM). A total of 200 teeth/cusps from 120 specimens were examined, and 46 specimens for non-tooth structures (see Additional file 1: Tables S1-S5 for sample sizes). CT data were thresholded and/or segmented in Avizo v. 9.6 (Thermo Scientific, Waltham, MA, USA). Surface files were processed in Geomagic Wrap v. 2015 (3D Systems, Rock Hill, SC, USA). The dataset includes all structures that were measured and considered sufficiently complete (essentially unworn and not broken), and no measured structures were excluded from the study for any other reason. We strove to include as wide a taxonomic range as possible, such as teeth from all orders of vertebrates, and from as many families of mammals as possible. Where permitted by institutions/individuals, 3D models are available at  on MorphoSource (https://www.morphosource.org/).
Power cascade shape analysis
For each 3D surface model of a tooth, we estimated a 3D midline through the centre of the tooth from tip to base. First, the midline was approximated by eye in Rhinoceros 3D v. 5.0 or 6.0 (Robert McNeel & Associates, Seattle, WA, USA). At 10 equally spaced locations along the estimated midline, test cross-sections through the tooth surface were generated perpendicular to the estimated midline. The centroid of each test cross-section was calculated, and the final midline was fit as a 3D spline running through these centroids. To measure the rate of radius increase, we placed 10 equally spaced cross-sections perpendicular to its midline (Fig. 2a) and calculated the average radius of each cross-section (Radius = √ [cross-sectional area/π]) using Rhinoceros 3D. Cross-sectional area includes all area interior to the external surface of the tooth, including pulp cavity. The use of tooth cross-sectional area perpendicular to the midline essentially removes the effect of the spiral (logarithmic or otherwise) so that our measurements concentrate on the rate of expansion of the radius of the tooth, not the trajectory of the midline in 3D space. All measurement data are included in Additional file 2: Data S1. Circumferential measurements of elephant and mammoth tusks were obtained from the literature [28, 53] and transformed to radius (Radius = [circumference/(2π)]). We then plotted log10 Distance from tip vs log10 Radius. An ordinary least squares (OLS) linear model was fit to each log Distance vs log Radius plot, with the R2 of the regression indicating goodness of fit to the predicted linear power cascade pattern using R Statistical Computing v. 3.5.0  and RStudio v. 1.1.447. Tooth development simulations were generated using the ToothMaker model [40, 41], and parameters and measurement data are in Additional file 3: Data S2. Non-tooth structures were scanned and measured in an analogous manner to teeth, and the measurement data are in Additional file 4: Data S3. The power cascade model was implemented with a logarithmic spiral growth pattern (forming a power spiral) in Mathematica v. 12.0 (Wolfram Research Inc., Champaign, IL), incorporating various cross-sectional shapes (circle, ellipse, lens, truncated circle)—see Fig. 7, Additional file 1: Supplementary Equations and Figure S10.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files. The 3D scan data that support the findings of this study are available on MorphoSource project Power Cascade (C1206) (https://www.morphosource.org/projects/0000C1206)  where permitted by institutions/individuals.
Maynard Smith J, Burian R, Kauffman S, Alberch P, Campbell J, Goodwin B, Lande R, Raup D, Wolpert L. Developmental constraints and evolution. Q Rev Biol. 1985;60(3):265–87.
Kavanagh KD, Evans AR, Jernvall J. Predicting evolutionary patterns of mammalian teeth from development. Nature. 2007;449(7161):427–32.
Newman SA. Physico-genetic determinants in the evolution of development. Science. 2012;338(6104):217–9.
Young NM, Winslow B, Takkellapati S, Kavanagh K. Shared rules of development predict patterns of evolution in vertebrate segmentation. Nat Commun. 2015;6:6690.
Abzhanov A. The old and new faces of morphology: the legacy of D’Arcy Thompson’s ‘theory of transformations’ and ‘laws of growth’. Development. 2017;144(23):4284–97.
Thompson DAW. On growth and form. 2nd ed. Cambridge: Cambridge University Press; 1942.
Raup DM. Geometric analysis of shell coiling: general problems. J Paleontol. 1966;40(5):1178–90.
Beldade P, Koops K, Brakefield PM. Developmental constraints versus flexibility in morphological evolution. Nature. 2002;416(6883):844–7.
Alberch P, Gale EA. A developmental analysis of an evolutionary trend - digital reduction in amphibians. Evolution. 1985;39(1):8–23.
Nijhout HF, Emlen DJ. Competition among body parts in the development and evolution of insect morphology. Proc Natl Acad Sci U S A. 1998;95(7):3685–9.
Wagner GP, Chiu CH, Laubichler M. Developmental evolution as a mechanistic science: the inference from developmental mechanisms to evolutionary processes. Am Zool. 2000;40(5):819–31.
Polly PD. Development and phenotypic correlations: the evolution of tooth shape in Sorex araneus. Evol Dev. 2005;7(1):29–41.
Jablonski D. Developmental bias, macroevolution, and the fossil record. Evol Dev. 2020;22(1–2):103–25.
Ungar PS. Mammal teeth: origin, evolution, and diversity. Baltimore: Johns Hopkins University Press; 2010.
Swindler DR. Primate dentition: an introduction to the teeth of non-human primates. Cambridge: Cambridge University Press; 2002.
Jernvall J, Kettunen P, Karavanova I, Martin LB, Thesleff I. Evidence for the role of the enamel knot as a control center in mammalian tooth cusp formation: non-dividing cells express growth stimulating Fgf-4 gene. Int J Dev Biol. 1994;38(3):463–9.
Jernvall J, Thesleff I. Reiterative signaling and patterning during mammalian tooth morphogenesis. Mech Dev. 2000;92(1):19–29.
Descartes R. Œuvres de Descartes: Correspondance II. Paris: Adam and Tanner, Léopold Cerf; 1638. 
Wallis J. Tractatus Duo: Prior, de Cycloide; Posterior, de Cissoide. Oxford: Lichfield; 1659.
Huxley JS. Problems of relative growth. London: Methuen & Co. Ltd.; 1932.
West G. Scale: the universal laws of life, growth, and death in organisms, cities and companies. New York: Penguin Publishing Group; 2017.
Fu DB, Hao HP. Numerical and experimental analysis of projectile nose geometry. J Spacecr Rocket. 2015;52(5):1510–4.
Tucker A, Sharpe P. The cutting-edge of mammalian development; how the embryo makes teeth. Nat Rev Genet. 2004;5(7):499–508.
Hakkinen TJ, Sova SS, Corfe IJ, Tjaderhane L, Hannukainen A, Jernvall J. Modeling enamel matrix secretion in mammalian teeth. PLoS Comput Biol. 2019;15(5):e1007058.
D’Amore DC, Harmon M, Drumheller SK, Testin JJ. Quantitative heterodonty in Crocodylia: assessing size and shape across modern and extinct taxa. PeerJ. 2019;7:e6485.
Dean C, Leakey MG, Reid D, Schrenk F, Schwartz GT, Stringer C, Walker A. Growth processes in teeth distinguish modern humans from Homo erectus and earlier hominins. Nature. 2001;414(6864):628–31.
Laws RM. A new method of age determination for mammals. Nature. 1952;169(4310):972–3.
Hanks J. Growth of the African elephant (Loxodonta africana). E Afr Wildl J. 1972;10:251–72.
Smith KM, Fisher DC. Sexual dimorphism and inter-generic variation in proboscidean tusks: multivariate assessment of American mastodons (Mammut americanum) and extant African elephants. J Mamm Evol. 2013;20(4):337–55.
Kuprina KV, Smorkatcheva AV. Noninvasive age estimation in rodents by measuring incisors width, with the Zaisan mole vole (Ellobius tancrei) as an example. Mammalia. 2019;83(1):64–9.
Pilgram T, Western D. Inferring the sex and age of African elephants from tusk measurements. Biol Conserv. 1986;36(1):39–52.
Urdy S, Wilson LAB, Haug JT, Sanchez-Villagra MR. On the unique perspective of paleontology in the study of developmental evolution and biases. Biol Theory. 2013;8:293–311.
Jernvall J. Linking development with generation of novelty in mammalian teeth. Proc Natl Acad Sci U S A. 2000;97(6):2641–5.
Skinner MM, Gunz P. The presence of accessory cusps in chimpanzee lower molars is consistent with a patterning cascade model of development. J Anat. 2010;217(3):245–53.
Paul KS, Astorino CM, Bailey SE. The Patterning Cascade Model and Carabelli’s trait expression in metameres of the mixed human dentition: exploring a morphogenetic model. Am J Phys Anthropol. 2017;162(1):3–18.
Evans AR, Daly ES, Catlett KK, Paul KS, King SJ, Skinner MM, Nesse HP, Hublin J-J, Townsend GC, Schwartz GT, et al. A simple rule governs the evolution and development of hominin tooth size. Nature. 2016;530:477–80.
Evans AR, Sanson GD. The effect of tooth shape on the breakdown of insects. J Zool. 1998;246(4):391–400.
Freeman PW, Lemen CA. The trade-off between tooth strength and tooth penetration: predicting optimal shape of canine teeth. J Zool. 2007;273(3):273–80.
Salazar-Ciudad I, Jernvall J. A computational model of teeth and the developmental origins of morphological variation. Nature. 2010;464(7288):583–6.
Harjunmaa E, Seidel K, Hakkinen T, Renvoise E, Corfe IJ, Kallonen A, Zhang Z-Q, Evans AR, Mikkola ML, Salazar-Ciudad I, et al. Replaying evolutionary transitions from the dental fossil record. Nature. 2014;512(7512):44–8.
Savriama Y, Valtonen M, Kammonen JI, Rastas P, Smolander O-P, Lyyski A, Hakkinen TJ, Corfe IJ, Gerber S, Salazar-Ciudad I, et al. Bracketing phenogenotypic limits of mammalian hybridization. R Soc Open Sci. 2018;5(11):180903.
Fowler DR, Meinhardt H, Prusinkiewicz P. Modeling seashells. Comput Graph. 1992;26(2):379–87.
Urdy S, Goudemand N, Bucher H, Chirat R. Allometries and the morphogenesis of the molluscan shell: a quantitative and theoretical model. J Exp Zool B Mol Dev Evol. 2010;314B(4):280–302.
Raup DM. Geometry of coiling in gastropods. Proc Natl Acad Sci U S A. 1961;47(4):602–9.
Ackerly SC. Kinematics of accretionary shell growth, with examples from brachiopods and mollusks. Paleobiology. 1989;15(2):147–64.
Moss ML. The logarithmic spiral and allometry as expressions of odontogenic thermodynamics. In: Kurten B, editor. Teeth: Form, Function, and Evolution. New York: Columbia University Press; 1982. p. 75–89.
Skalak R, Dasgupta G, Moss M, Otten E, Dullemeijer P, Vilmann H. Analytical description of growth. J Theor Biol. 1982;94(3):555–77.
World Health Organisation. World Health Organisation Child Growth Standards. France: World Health Organisation; 2006.
Newman SA, Bhat R. Dynamical patterning modules: a “pattern language” for development and evolution of multicellular form. Int J Dev Biol. 2009;53(5–6):693–705.
Johnson AB, Fogel NS, Lambert JD. Growth and morphogenesis of the gastropod shell. Proc Natl Acad Sci U S A. 2019;116(14):6878–83.
Held LI Jr. Models for embryonic periodicity. Monogr Dev Biol. 1992;24:1–119.
Evans AR, Pollock TI, Cleuren SGC, Parker WMG, Richards HL, Hocking DP, Adams JW, McCurry MR, Chieu L, Courtney A et al: Power Cascade: 3D models of teeth, claws, horns, thorns, beaks, and shells. MorphoSource database, Project C1206. https://www.morphosource.org/projects/0000C1206. 2021.
Grigoriev SE, Fisher DC, Obada T, Shirley EA, Rountrey AN, Savvinov GN, Garmaeva DK, Novgorodov GP, Cheprasov MY, Vasilev SE, et al. A woolly mammoth (Mammuthus primigenius) carcass from Maly Lyakhovsky Island (New Siberian Islands, Russian Federation). Quat Int. 2017;445:89–103.
R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.
We thank A. Courtney, E. Delson, M. DeVeer, A. Farke, D. Field, A. Haque, M. McCurry, A. McDonald, O. Panagiotopoulou, R. Rozzi, J. Rule, E. Sherratt, M. Skinner, and G. Wilson for 3D scanning or existing 3D scan data. We thank B. Weir, D. Annese (Monash University), K. Rowe, K. Roberts, R.-L. Erickson, T. Ziegler, R. Schmidt, M. Gomon, and D. Bray (Museums Victoria) for access to specimens under their care. We thank N. Campione and M. Camaiti for assistance in translating the Descartes and Wallis works, respectively. The authors acknowledge the technical assistance of the Australian National Imaging Facility fellows and equipment at the Monash Biomedical Imaging facility and Monash X-ray Microscopy for Imaging Geomaterials (XMFIG) (ARC LIEF LE130100115 to A.R.E.). For comments and discussion, we thank A. Courtney, G. Sanson, J. Jernvall, D. Smyth, E.S. Daly, F. Smith, D. Jones, K. Brink, C. Jackson, E. Sibert, A. Hsiang, K. Tsyganov, and L. Wilson.
A.R.E. was supported by an Australian Research Council Future Fellowship (FT130100968), Linkage Project (LP150100403) and Discovery Project (DP180101797), and Monash University Science and Interfaculty Seed grants. Support was given by Australian Government Research Training Program Scholarships (W.M.G.P., T.I.P., H.L.R. and K.L.S.G.), Monash Graduate Scholarship and Monash International Tuition Scholarship (S.G.C.C.), Monash Graduate Excellence Scholarship (T.I.P. and K.L.S.G.), Robert Blackwood Partnership PhD Top-Up Scholarships (W.M.G.P., H.L.R. and K.L.S.G.), and Holsworth Wildlife Research Endowments (S.G.C.C. and T.I.P.).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Vertebrate teeth show power cascade growth. Figure S2. Power cascade shapes are characterized as surfaces of revolution for power functions, with variables Slope, Intercept and MaxDistance. Figure S3. Alternative sampling intervals along an elephant Loxodonta africana NMV C30765 tusk. Figure S4. Two ways in which tooth can deviate from linear power cascade: tip offset and missing tip. Figure S5. Power cascade shapes are self-similar curves. Figure S6. In silico tooth development models do not produce cusps that closely approximate power cascade found in natural teeth. Figure S7. Logarithmic spiral, shell model, power cascade model and power spiral model. Figure S8. Pointed structures in vertebrates, invertebrates and plants show power cascade growth. Figure S9. Prickle growth in roses causes deviation from power cascade growth. Figure S10. Power cascade interface implemented in Mathematica for generating biological shapes using power cascade and logarithmic spiral. Figure S11. Graphical abstract – Power cascade combined with the logarithmic spiral can generate many biological shapes. Table S1. Number of species, specimens and structures in each class for all structures (teeth and non-teeth) measured in this study. Table S2. Number of species, specimens and teeth/cusps in each class for all teeth measured in this study. Table S3. Number of species, specimens and teeth/cusps in each mammalian order for all teeth measured in this study. Table S4. Number of species, specimens and non-tooth structures in each class for all non-tooth structures measured in this study. Table S5. Number of species, specimens and structures for each type of structure measured in this study. Supplementary Discussion. Resampling of power cascade variables. Effect of tip offset on power cascade linear pattern. Supplementary Equations. Derivation of power cascade growth mechanism. Mathematica implementation of power cascade model.
Distance and cross-sectional area data for vertebrate teeth, cusps and tooth rows. Taxonomy, specimen number, institution, tooth type/position and measurements of cross-sectional area for 10 distances from tip of tooth.
Distance and cross-sectional area data for developmental simulations of seal teeth from model by Savriama et al. . Model parameters and measurements of cross-sectional area for 10 distances from tip of tooth. Starting from the ringed seal model of Savriama et al. , parameters were increased and decreased to show the effect of each parameter on tooth shape. The central cusp on each model was measured for the power cascade. Parameters: Act, activator; Boy, buoyancy; Deg, degradation of activator; Egr, epithelial growth; Inh, inhibitor.
Distance and cross-sectional area data for non-tooth structures: antlers, beaks, chelicera, claws, hooves, horns, prickles, shells, spine and thorn. Taxonomy, specimen number, institution, structure type and measurements of cross-sectional area for 10 distances from tip of structure.
About this article
Cite this article
Evans, A.R., Pollock, T.I., Cleuren, S.G.C. et al. A universal power law for modelling the growth and form of teeth, claws, horns, thorns, beaks, and shells. BMC Biol 19, 58 (2021). https://doi.org/10.1186/s12915-021-00990-w