BACKGROUND
organisms that need to perform multiple tasks face a fundamental tradeoff: no phenotype can be optimal at all tasks  <cit> .
this tradeoff situation is reminiscent of tradeoffs in economics and engineering.
these fields analyze tradeoffs using pareto optimality theory  <cit> .
pareto optimality was recently used in biology to study tradeoffs in evolution  <cit> .
in contrast to the classic fitness-landscape approaches in which organisms maximize a single fitness function  <cit> , the pareto approach deals with several performance functions, one for each task, that all contribute to fitness .figure  <dig>
an overview of pareto theory for evolutionary tradeoffs.
the classical viewpoint of a fitness landscape: phenotypes are arranged along the slopes near the peak of a fitness hill maximum.
in contrast, the pareto viewpoint suggests a tradeoff between tasks.
for each task there is a performance function, which is maximal at a point known as the archetype for that task.
the fitness function in each niche is a combination of the different performance functions .
optimality in a niche in which task  <dig> is most important, is achieved near archetype  <dig> .
optimality in a niche in which all tasks are equally important, is achieved close to the middle of the pareto front .
the entire pareto front- the set of maxima of all possible fitness functions that combine these performances- is contained within the convex hull of the archetypes.
different numbers of tasks give various polygons or polyhedra, generally known as polytopes.
two tasks lead to a suite of variation along a line segment.
three tasks lead to a suite of variation on the triangle whose vertices are the three archetypes.
four archetypes form a tetrahedron.
this is true while there are enough traits measured: in lower dimensional trait spaces one finds projections of these polytopes.
pareto theory makes strong predictions on the range of phenotypes that evolve in such a multiple-objective situation: the evolved phenotypes lie in a restricted part of trait-space, called the pareto front.
the pareto front is defined as phenotypes that are the best possible compromises between the tasks; phenotypes on the pareto front can’t be improved at all tasks at once.
any improvement in one task comes at the expense of other tasks.
shoval et al.
<cit>  calculated the shape of the pareto front in trait space under a set of general assumptions.
evolved phenotypes were predicted to lie in a polygon or polyhedron in trait space, whose vertices are extreme morphologies, called archetypes, which are each optimal at one of the tasks .
thus, two tasks lead to phenotypes on a line that connects the two archetypes, three tasks to a triangle, four tasks to a tetrahedron and so on .
these polyhedra can have slightly curved edges in some situations  <cit> .
one does not need to know the tasks in advance: tasks can be inferred from the data, by considering the organisms closest to each archetype.
this theory can be rejected in principle by datasets which lie in a cloud without sharp vertices, and hence do not fall into well-defined polygons.
the shoval et al.
theory has been applied so far to datasets from animal morphology  <cit> , bacterial gene expression  <cit> , cancer  <cit> , biological circuits  <cit>  and animal behavior  <cit> .
in all of these cases, multi-dimensional trait data was found to be well described by low-dimensional polygons or polyhedra .
tasks were inferred based on the properties of the organisms  closest to the vertices.
an algorithm for detecting polyhedra in biological data and inferring tasks was recently presented  <cit> .
however, some of the fundamental predictions of the theory have not been tested yet.
the theory predicts that as long as the tasks stay more-or-less constant, for example dictated by biomechanical constraints, the vertices of the polygon also do not change.
moreover, the polygons in the theory are not necessarily due to phylogenetic history, but rather to convergent evolution to pareto-optimal solutions.
thus, for example, after a mass extinction which removes most of the species from a class  <cit> , survivor species are predicted to evolve to re-fill the same polygon as their ancestors  <cit> .
to test these predictions requires a class of organisms evolving over geological timescales, with mass extinctions, and whose geometry is well-defined and can be linked to function.
an excellent model system for this purpose is ammonoid fossil shells.
ammonoids were a successful and diverse group of species, which lived in the seas from  <dig> to  <dig> million years ago .
ammonoid shells can be described by a morphospace defined by three geometrical parameters, defined in the pioneering work of david raup   <cit> .
in this morphospace, the outer shell is a logarithmic spiral, whose radius grows with each whorl by a factor w, the whorl expansion rate.
there is a constant ratio between the inner and outer shell radii, denoted d. finally, the shell cross section can range from circular to elliptical, as described by s, the third parameter.
raup’s w-d-s parameterization can be robustly measured from fossils  <cit>  although the coiling axis changes throughout ontogeny and thus, the coiling axis is sometimes difficult to exactly locate in actual specimens  <cit> .
it has been the setting for extensive research on ammonoid morphology and evolution  <cit> , as well as the morphology of other shelled organisms  <cit> .figure  <dig>
raup morphospace coordinates.
ammonoid shell morphology can be described by three dimensionless geometrical parameters: w, the whorl expansion rate, is defined by a/b in the figure.
d, the internal to external shell ratio, is x/a.
s, the opening shape parameter, is y/z.
the shell diameter can also be related to the parameters in this figure as shown.
plotting each genus of ammonoids as a point in this morphospace, ignoring coiling axis changes, raup discovered that most of the theoretical morphospace is empty: many possible shell forms are not found.
the existing forms lie in a roughly triangular region in the w-d plane .
one reason for this distribution is geometric constraints.
researchers have suggested that ammonoids tend to lie to the left of the hyperbola w = 1/d  <cit> , because beyond this curve shells are gyroconic  .
such gyroconic shells are mechanically weaker and less hydrodynamically favorable  <cit> .
it is noteworthy, however, that shells to the right of the curve do exist in nature, for example in the bactritida or orthocerida lineages, which are probably ancestral to the ammonoids   <cit> , as well as in heteromorph ammonoids that occasionally occur in the mesozoic and more commonly in the cretaceous.
thus the w = 1/d curve is unlikely to be an absolute geometric constraint .figure  <dig>
ammonoids repeatedly filled the same triangle in d-w plane after mass extinctions.
all of the ammonoid data used in the present study.
red points are genera before the ff  mass extinction, genera after ff are denoted by blue points.
the green curve is w = 1/d.
ammonoids before the ff extinction, together with a schematic arrow for the direction of evolution from ancestral taxa.
genera between ff and dm mass extinctions fill out a triangle , surviving genera from the ff mass extinction are denoted by red bold points.
ammonoids between dm and pt mass extinctions fill a triangle, surviving genera from the dm mass extinction are denoted by red bold points.
ammonoids after the pt mass extinction fill a triangle, surviving genera from the pt mass extinction are denoted by red bold points.
ammonoids from different periods, together, genera between ff and dm are denoted blue, dm to pt in red and post pt in green.
the shell morphologies of the three archetypes at the vertices of the triangle are shown.
studies in recent years have considered a larger dataset of ammonoids than raup  <cit> .
work, saunders and nikolaeva  <cit>  show that after each mass extinction, ammonoid genera refill the same roughly triangular morphospace  <cit> .
the repeated convergence to the same suite of variation raises the question of the relation between ammonoid morphology and function.
most studies hypothesize a fitness function, which has an optimum in the middle of the triangular region  <cit>  .
the fitness function is often taken to be dominated by hydrodynamic drag; this assumption is compelling since the contours of hydrodynamic efficiency, experimentally measured by chamberlain  <cit> , show peaks at positions close to the most densely occupied regions of morphospace  <cit> .
the ammonoid genera are assumed to also occupy the slopes that descend from the fitness peaks, until bounded by the geometric constrains  <cit> .
interestingly, raup did not espouse the idea of a single task  dominating fitness, but rather noted that multiple tasks might be at play  <cit> .
in every niche, different tasks become important, leading to niche-dependent fitness functions with different maxima .
the idea of multiple tasks was elegantly employed by westermann  <cit> , who described ammonoid morphospace by mapping it to a triangle.
at the vertices are three ‘end member’ morphologies which correspond to three lifestyles.
each morphology is mapped to a point in the triangle, which is interpreted as portraying the relative distance from the end members and hence the relative weights of the three lifestyles.
the westermann morphospace was useful in comparing different datasets and in interpreting ammonoid lifestyles  <cit> .
the main drawback of the westermann morphospace is that, because it involves nonlinear dimensionality reduction, different morphologies can be mapped to the same point, and in some cases slight differences in shape can lead to large differences in the westermann projection.
thus, it is of interest to seek a relation between shape and tasks without such drawbacks.figure  <dig>
the performance contours of the three putative tasks for ammonoid shells.
contours for shell economy, defined as the ratio of internal volume to shell volume, with red denoting high values, and blue low values.
for gyroconic shells , this performance function becomes constant, and equal to the lowest contour shown .
the triangle encapsulating the entire ammonoid dataset is shown in black.
contours for the drag coefficient measured by chamberlain  <cit> , red lines denote lower drag or better hydrodynamics.
contours for the growth function defined in the main text, red lines denote quicker growth.
the contours of the three tasks give rise to a suite of variation denoted by blue points.
to address this, we explore evolutionary tradeoffs between tasks in the framework of pareto optimality theory, to quantitatively explain the suite of variation in direct morphospace , and to infer the putative tasks at play.
we find that ammonoid morphology in the w-d-s morphospace falls within a square pyramid, suggesting five tasks.
the triangular region observed by raup is the projection of this pyramid on the w-d plane, and the westermann morphospace is a dimensionality reduction of the three-dimensional pyramid to a two-dimensional triangle.
we propose putative tasks whose performance contours jointly lead to the observed suite of variations, including hydrodynamic efficiency, shell economy, compactness and rapid shell growth.
the position of each species in this pyramid, namely its distance from each vertex, indicates the relative importance of each task in the niche in which that species evolved.
after the ff and dm mass extinctions , surviving ammonoids refill essentially the same pyramid.
after the pt extinction , part of the pyramid is refilled.
these findings lend support to the pareto theory of evolutionary tradeoffs in the context of evolution on geological timescales.
RESULTS
ammonoid distributions in the w-d plane converge to a similar triangle after major extinctions
we begin by considering ammonoid morphology in the w-d plane, and later consider the three dimensional w-d-s space .
we combine the data of saunders, work and nikolaeva  <cit>  for paleozoic ammonoids , with the data of mcgowan  <cit>  for mesozoic ammonoids .
the data is classified into three sets between mass extinctions: from ff to dm , from dm to pt , and after pt .
we tested whether the ammonoid distribution in each set falls in a triangle more closely than randomized data, based on the statistical test of  <cit> .
we use an archetype analysis algorithm   <cit>  to find triangles, which enclose as much of the data as possible.
we find that a triangle describes each dataset much better than randomized datasets in which the w and d coordinates are randomly permuted .
randomized data rarely fill out a triangle as well as the real data .
we next tested how similar the triangles are for the three datasets.
we computed the ratio between the intersection area of the triangles to the union area as a measure for triangle similarity.
the three triangles show large ratios of intersection to union area ,  and  pairs respectively, p <10- <dig> compared to randomly generated triangles, see methods), indicating that the triangles are very similar.
we conclude that after each extinction, ammonoids re-populate essentially the same triangular region.
the vertices of the triangle describing the joint dataset of ammonoids after ff  are  <dig> d1*,w1*= <dig> , <dig> ~ <dig> , <dig>    <dig> d2*,w2*= <dig> , <dig> ~ <dig>    <dig> d3*,w3*= <dig> , <dig> ~ <dig> .
<dig>
we next ask which tasks might relate to each of the vertices.
economy of shell material may determine the first archetype
raup  <cit>  suggested that a possible need of the ammonoids is to maximize their internal volume relative to shell volume.
this is important if shell production is costly, and also in terms of buoyancy considerations.
ammonoids are thought to achieve neutral buoyancy by balancing shell weight with buoyancy from their air-filled chambers; high internal volume relative to shell material extends the range over which neutral buoyancy can be reached  <cit> .
to calculate shell material relative to internal volume at each point in morphospace, we follow raup and assume that shell thickness is a fixed fraction of radius, namely thickness/radius =  <dig> , as measured by  <cit> .
interestingly, this ratio is close to the optimal ratio of thickness/radius = <dig>  from calculations of mechanical strength in tube-like bones  <cit> .
we improve slightly on chamberlin and raup’s original calculation  <cit>  by numerically evaluating the necessary integrals rather than using the analytical approximations of  <cit>  , yielding corrections of about 10%.
the maximum of internal volume relative to shell thickness occurs at  = .
this point is close to archetype one d1*,w1*= <dig> , <dig> .
the calculated contours of internal volume relative to shell thickness-namely the performance contours of the task of economy- have a curving ridge that points towards the third archetype .
performance drops sharply on either side of this ridge.figure  <dig>
the three dimensional pareto front of the ammonoid dataset.
the rms error for pcha optimal polygons and polyhedra is computed for different numbers of possible vertices: line, triangle, tetrahedron, 5-vertex polyhedron, etc.
error decreases with increasing the number of archetypes up to  <dig>  increasing the number beyond  <dig> doesn't improve the fit by much .
the best fit 5-archetype polygon resembles a square pyramid.
blue points denote ff to dm ammonoids, red are dm to pt and green are post pt ammonoids.
archetypes are numbered, their morphology is shown, and the suggested tasks are listed in panel a.
the second archetype may optimize hydrodynamics
we conjecture that the second archetype maximizes the hydrodynamic efficiency of the ammonoids.
low drag is important for ammonoids in order to swim rapidly.
hydrodynamic efficiency is measured by the drag coefficient, which is a dimensionless number specific to each geometrical shape.
the drag coefficient is proportional to the force which should be applied in order to keep an object of a given surface area moving at a given velocity in water.
drag coefficients were measured by chamberlain  <cit>  using plexiglass models of shells  <cit> .
the contours of hydrodynamic efficiency are shown in figure 4b.
drag monotonically increases with d and w, hence we can conclude that the ammonoid morphology with minimal drag has the lowest possible values of d and w, namely  =  this is close to the second vertex of the triangle, archetype two at d2*,w2*= <dig> , <dig> .
the third archetype may optimize rapid shell growth
the remaining vertex of the triangle, archetype  <dig>  has a large value of w. thus the shell radius at this vertex increases rapidly with each revolution of the spiral .
there are different possible tasks that might relate to large w, including rapid growth, shell-orientation and swimming capabilities.
in westermann morphospace, large w compared to d and s is interpreted as nektonic .
here, we wish to demonstrate an essential approach, and thus focus on one of these potential tasks: rapid growth, and leave other possibilities to future study.
the fossil dataset we use does not contain information on growth.
however, if we assume that the ability to generate shell material  is proportional to body mass  we can predict the growth, or at least a function proportional to growth, using only the dimensionless parameters we have.
an evolute shell allows volume to grow rapidly with each whorl.
rapid growth may be important because predation tends to decrease with organism size.
this would select for increased w. however, the whorl expansion rate w cannot grow without bound in order to avoid cyrtoconic shells- the shell must close over itself at least once to provide space for the ammonoid body .
a coiled shell is also important in order to benefit from increased shell thickness, because until the ammonoid is closed, the thinner shell is exposed to the outside.
the small value of d is also reasonable for such a function, because when w is large, a small d is a must in order to benefit from the advantages of w < 1/d .
a similar function was suggested in snails where shell growth rate was found to be larger in snails in the presence of predators  <cit> .
to be concrete, we consider a putative performance function that penalizes the ammonoid for the smallness of its diameter, namely p3=∫0∞1diamtdt .
contours of this performance function are shown in figure 4c.
the function peaks at  =  close to the third archetype d3*,w3*= <dig> .
<dig>  at this archetype, ammonoids reach large diameters most rapidly.
one may ask if the advantage of growth comes from the increased diameter which might make the ammonoid too large for specific predators, or from the increased shell thickness which make it stronger.
it is difficult to distinguish between this two conjectures since from  <cit>  we know that this quantities are proportional to one another.
it is likely that both diameter and shell thickness contribute to fitness.
the three putative performance functions, shell economy, hydrodynamic efficiency, and shell growth together give rise to a triangular shaped pareto front.
the pareto front boundaries are given by the points of tangency of the contours of the different performance functions.
figure 4d shows the computed pareto front, which resembles a slightly curved triangle, and is similar to the observed suite of variation.
ammonoid data is enclosed by a pyramid in w-d-s morphospace
up to now, we considered ammonoid morphology projected on the w-d plane.
we now turn to the analysis of the data in the three-dimensional morphospace, given by w,d and s—whorl expansion, radii ratio and the shape of the shell opening.
low values of the parameter s correspond to oblate elliptical openings, giving rise to compressed shells .
an s value of  <dig> corresponds to a circular shell opening; high values corresponding to depressed shell morphologies .figure  <dig>
size is enriched at some of the archetypes.
ammonoid shell diameter as a function of distance from each archetype shows that small diameters are prevalent near archetypes  <dig> and  <dig>  data includes diameter for  <dig> genera   <cit> , divided into  <dig> bins with equal number of genera according to the distance from each archetype.
average diameter for each bin is plotted .
for convenience, a fit of the averages to a line is shown.
no diameter enrichment near archetype  <dig> .
positive diameter enrichment near archetype  <dig>   positive diameter enrichment near archetype  <dig> .
negative diameter enrichment near archetype  <dig>   negative diameter enrichment near archetype  <dig> .
we attempted to enclose the 3d dataset by polygons with  <dig> to  <dig> vertices.
we evaluated the extent to which each polygon explains the data, by calculating the rms distance of points outside the polyhedron.
we find that beyond  <dig> vertices, the rms error does not decrease significantly : shapes with  <dig> or more vertices do not improve the closeness of fit appreciably.
hence a 5-vertex polygon is a parsimonious description of the data .
this 5-vertex shape has four vertices that lie approximately on a plane.
we thus consider this shape as a pyramid.
a square pyramid encloses the data better than randomly permuted dataset with p <10- <dig> .
the five vertices of the pyramid suggest five archetypes, whose coordinates are given in table  <dig>  the square base of the pyramid has two vertices at low s , and two others, which match them for w and d values, but have higher s values .
the apex of the pyramid has a thin opening with s =  <dig> .table  <dig>
summary of suggested ammonoid archetypes
suggested task
w
s
d
archetype number
coordinates of the archetypes found by the archetype analysis algorithms with 5-archetypes, along with their putative tasks.
projecting the pyramid on the w-d plane, we find that the apex of the pyramid matches the ‘growth’ archetype described above; the ‘economy’ and ‘hydrodynamic’ archetypes each corresponds to the projection of two 3d archetypes: the economy archetype corresponds to archetypes  <dig> and  <dig>  and the hydrodynamic archetype to archetypes  <dig> and  <dig> .
economy, hydrodynamic and growth performance functions are maximized near three of the pyramid vertices
we repeated the calculation of economy performance  in three dimensions.
the 2d contours shown previously  were evaluated at s =  <dig>  by varying s, we find that the maximal economy is found at  = .
this is reasonably close to vertex  <dig> of the observed pyramid d1*,s1*,w1*= <dig> , <dig> , <dig> .the internal volume to shell volume ratio in this vertex is 96% of the optimum value.
for comparison, this ratio drops to nearly zero near vertices  <dig> and  <dig> of the pyramid.
the hydrodynamic efficiency measured in  <cit>  includes data at values of s other than s =  <dig>  this indicates that optimal hydrodynamic efficiency is at low s values, i.e.
s→ <dig>  the resulting optimum is thus close to vertex  <dig> of the pyramid, which is d2*,s2*,w2*= <dig> , <dig> , <dig> : note the low values of d,s and w.
archetype  <dig> has an s value close to  <dig> .
the dependence of the growth performance function on s comes only implicitly through the volume-to-surface ratio.
it is unclear from the present simplified model for the growth performance function why  <dig>   is selected as the optimal s value for archetype  <dig>  this s value might be due, for example, to diminishing returns of shell production per body mass.
in other words, the assumption that shell material production is constantly proportional to body mass might be imprecise.
if shell production grows slower than linearly with body mass , this will favor smaller-volume ammonoids with smaller value of s that will increase diameter faster.
the last two pyramid archetypes may be related to size
two pyramid vertices remain to be explained, vertices  <dig> and  <dig>  these vertices have large values of s, and correspond to depressed shells .
we find that these shapes have the smallest ratio of surface area to volume .
they are therefore the most globular in the suite of variation, in the sense that their height is most similar to their width and depth.
one feature of globular ammonoids is small size for a given internal volume, because spherical shapes have the minimal diameter of all shapes with the same volume.
up to now, we did not consider the absolute size of the ammonoids, only on dimensionless shape traits w, d and s. to address this, we correlated data by mcgowan  <cit>  on ammonoid size  with distance from the five vertices of the pyramid.
we find an enrichment of small ammonoids most strongly near archetypes  <dig> and 5: the genera nearest to these vertices have the smallest diameters .
archetypes  <dig> and  <dig> are enriched with large ammonoids and archetype  <dig> has weak enrichment since its s value  is relatively larger than archetypes  <dig> and  <dig> .
archetypes  <dig> and  <dig> may thus correspond to economy and hydrodynamic tasks respectively, combined with a need for smallness.
this relation between diameter and globularity is in line also with  <cit> , which used a different dataset.figure  <dig>
archetypes in trait space and in westermann space.
ammonoids in westermann space.
the red ellipses are the pyramid archetypes projected on the westermann space, with 5% error in the archetype positions.
three of the pyramid archetypes lies near the vertices of the triangle, another archetype is near the edge because both d and s are dominant.
the region near archetype  <dig> is severely warped  because d,s and w are all relatively small.
the pyramid in the w-d-s trait space, the ellipsoids are 5% errorbars around each vertex.
these small ellipsoids translates to the red ellipses of subfigure a when switching to westermann morphospace.
we further compared the way ammonoids from different periods fill out the pyramid.
the main difference between periods is between paleozoic and mesozoic genera.
mesozoic ammonoids tend to have lower s values than paleozoic ones, as found by mcgowan  <cit> .
in the pyramid, they are more densely arrayed near the face defined by vertices  <dig>   <dig> and  <dig>  and away from  <dig> and especially from  <dig>  this may be interpreted in the present framework as a shift in the niches occupied by later ammonoids, in which tasks corresponding to archetypes  <dig> and  <dig> contribute less to fitness than they did in the paleozoic niches.
finally, we mapped the five archetypes of the pyramid to the westremann morphospace.
we find that that three archetypes,  <dig>   <dig> and  <dig>  map near the three vertices of the westermann triangle .
the two other vertices of the pyramid map closer to the edges of the triangle.
some of the archetypes map slightly outside of the triangle since they are exptrapolated points which lie outside of the ammonite dataset.
we also asked about the sensitivity of this transformation, by testing a small region around each archetype .
we find that one of the archetypes, archetype  <dig>  lies in a region of morphospace which is severely warped by the westermann transformation, and maps to a wide region in the triangle.
the other archetypes are less sensitive and map to relatively small regions of the triangle .
discussion
this study explored how tradeoffs between multiple tasks may have contributed to the evolution of ammonoid shell morphology.
ammonoid shell data on  <dig> genera were studied in raup’s three parameter morphospace.
the data is well described by a square pyramid.
this finding is interpreted in light of pareto theory on tradeoffs between tasks.
the five vertices of the pyramid may be interpreted as archetype morphologies optimal for a single task, and morphologies in the middle of the pyramid are generalists which compromise between the tasks.
we propose candidate tasks for the archetypes.
hydrodynamic efficiency is a good candidate for one of these tasks, and is maximized near vertex  <dig> of the pyramid .
other putative tasks can be inferred from the position of the other vertices of the pyramidal shell distribution.
we propose that economy of shell material  is a second task, quantitated by the ratio of internal volume and total shell material.
the maximum of this function matches one of the vertices of the pyramid.
a third task may be rapid growth.
a performance function relating to rapid growth of ammonoid diameter is maximal near the apex of the pyramid, at shells with high w. two other tasks may relate to small spherical-like shells combined with low drag and high economy.
it is interesting to relate this study with previous work by westermann and ritterbush based on the idea that ammonoids face tradeoffs between different tasks, which determine their morphologies  <cit> .
westermann proposed a morphospace which, instead of working in d-w-s space, works in a 2-dimensional projection which consists of ratios of related measurements.
westermann morphospace has many advantages.
as a method to reduce 3-dimensional data into 2-dimensional one, it helps visualize data in order to achieve better understanding of the geometry.
it is also useful in understanding the different niches that ammonoids occupy and infer the various tasks they face  <cit> .
westermann's 2-dimensional representation also has drawbacks.
as a dimensionality reduction method, it loses information about the data.
ammonoid shells with very different geometries can be mapped to the same point in westermann morphospace.
moreover, because the westemann map is nonlinear, there are regions in morphospace that map to the triangle with relatively large errors.
for example, a small region around the point of minimal values of d, s and w  can mapped to the entire westermann triangle  depending on slight variation in the values.
the present approach does not show these drawbacks because it works directly in w-d-s morphospace.
it thus distinguishes between morphologies which are mapped to the same westermann point.
the square pyramid identified here suggests two new tasks  in addition to westermanns three.
these are archetypes  <dig> and  <dig> which correspond to hydrodynamic drag and compact shell economy.
we also propose a different interpretation of the other three tasks.
for example the westermann oxycone endmember, which is linked to nektonic lifestyle, corresponds to our archetype  <dig> .
furthermore, the same oxyconic endmember relates to certain morphologies near our archetype  <dig> .
the present study also bears on the question of geometric constraints in evolution  <cit> .
the w = 1/d line in ammonoids is thought to be an outstanding example of a geometric constraint  <cit> , because of the disadvantages of the open shell morphology beyond this curve.
this assumption is challenged by the existence of organisms with w > 1/d, including lineages ancestral to ammonoids as well as several heteromorphs  <cit> .
the present approach can make the concept of geometric constraint more precise by relating it to biological tasks.
we consider the performance functions of tasks, some of which indeed show a decline beyond the w = 1/d line.
in particular, economy and hydrodynamics contours both begin to sharply decline when w > 1/d .
this provides a more principled explanation, replacing strict geometric constraint with the more subtle dependence of specific performance functions on geometry.
other taxa may perform a different set of tasks, including a task with an archetype in the ‘forbidden for ammonoid’ region, w > 1/d.
such tasks might explain the morphology of the taxa which show gyroconic shells.
an alternative view is that some characters states do not require a functional explanation, but rather were neutral enough for a clade to succeed for some time.
this study adds to previous studies that used the pareto approach to analyze other biological systems  <cit> .
these systems showed lines, triangles or tetrahedra in morphospace.
ammonoids are the first system in which a pyramidal pareto front is observed.
for this purpose, we find that the archetype analysis algorithm pcha  <cit>  is an efficient way to detect high order polyhedra in data  <cit> .
the present approach can be readily extended to other shelled organisms such as gastropods and bivalves.
one application of the present approach is a quantitative inference of which task is important for fitness in the particular niche of each genus.
the closer the shell morphology is to a given vertex of the pyramid, the more important the corresponding task.
since ammonoid shells are carried by currents and found in rocks far from the habitat of the living organism, it is challenging to connect morphology with behavior.
the present approach can offer quantitative inference about the relative contribution of tasks to fitness, to provide insight into the ecological niche of these extinct organisms.
more generally, this study supports basic predictions of the pareto theory for evolutionary tradeoffs  <cit> , which we hope will be useful also for other biological contexts.
CONCLUSIONS
this study supports fundamental predictions of the pareto theory of tradeoffs by shoval et al.
<cit>  that have not been previously tested on the scale of hundreds of millions of years of evolution.
ammonoid shell data on  <dig> genera is well-described by a square pyramid in morphospace.
the five vertices of the pyramid may be interpreted as archetype morphologies optimal for a single task.
inferred tasks include shell economy, rapid growth, compactness and hydrodynamic efficiency.
the vertices of the polygon do not change over the timescale of interest, as predicted in the case where the tasks stay more-or-less constant because they are dictated by biomechanical considerations.
moreover, the polygons and polyhedra in the theory are not necessarily due to phylogenetic history, but rather to convergent evolution to pareto-optimal solutions.
this agrees with the finding that after a mass extinction which removes almost all of the species, survivor species evolve to re-fill the same polygon as their ancestors.
this approach may be used to infer biological tasks from data in other biological contexts.

 organisms that need to perform multiple tasks face a fundamental tradeoff: no phenotype can be optimal at all tasks. this tradeoff situation is reminiscent of tradeoffs in economics and engineering. these fields analyze tradeoffs using pareto optimality theory. pareto optimality was recently used in biology to study tradeoffs in evolution. in contrast to the classic fitness-landscape approaches in which organisms maximize a single fitness function, the pareto approach deals with several performance functions, one for each task, that all contribute to fitness.figure an overview of pareto theory for evolutionary tradeoffs. the classical viewpoint of a fitness landscape: phenotypes are arranged along the slopes near the peak of a fitness hill maximum. in contrast, the pareto viewpoint suggests a tradeoff between tasks. for each task there is a performance function, which is maximal at a point known as the archetype for that task. the fitness function in each niche is a combination of the different performance functions. optimality in a niche in which task is most important, is achieved near archetype. optimality in a niche in which all tasks are equally important, is achieved close to the middle of the pareto front. the entire pareto front- the set of maxima of all possible fitness functions that combine these performances- is contained within the convex hull of the archetypes. different numbers of tasks give various polygons or polyhedra, generally known as polytopes. two tasks lead to a suite of variation along a line segment. three tasks lead to a suite of variation on the triangle whose vertices are the three archetypes. four archetypes form a tetrahedron. this is true while there are enough traits measured: in lower dimensional trait spaces one finds projections of these polytopes. pareto theory makes strong predictions on the range of phenotypes that evolve in such a multiple-objective situation: the evolved phenotypes lie in a restricted part of trait-space, called the pareto front. the pareto front is defined as phenotypes that are the best possible compromises between the tasks; phenotypes on the pareto front can’t be improved at all tasks at once. any improvement in one task comes at the expense of other tasks. shoval et al. calculated the shape of the pareto front in trait space under a set of general assumptions. evolved phenotypes were predicted to lie in a polygon or polyhedron in trait space, whose vertices are extreme morphologies, called archetypes, which are each optimal at one of the tasks. thus, two tasks lead to phenotypes on a line that connects the two archetypes, three tasks to a triangle, four tasks to a tetrahedron and so on. these polyhedra can have slightly curved edges in some situations. one does not need to know the tasks in advance: tasks can be inferred from the data, by considering the organisms closest to each archetype. this theory can be rejected in principle by datasets which lie in a cloud without sharp vertices, and hence do not fall into well-defined polygons. the shoval et al. theory has been applied so far to datasets from animal morphology, bacterial gene expression, cancer, biological circuits and animal behavior. in all of these cases, multi-dimensional trait data was found to be well described by low-dimensional polygons or polyhedra. tasks were inferred based on the properties of the organisms closest to the vertices. an algorithm for detecting polyhedra in biological data and inferring tasks was recently presented. however, some of the fundamental predictions of the theory have not been tested yet. the theory predicts that as long as the tasks stay more-or-less constant, for example dictated by biomechanical constraints, the vertices of the polygon also do not change. moreover, the polygons in the theory are not necessarily due to phylogenetic history, but rather to convergent evolution to pareto-optimal solutions. thus, for example, after a mass extinction which removes most of the species from a class, survivor species are predicted to evolve to re-fill the same polygon as their ancestors. to test these predictions requires a class of organisms evolving over geological timescales, with mass extinctions, and whose geometry is well-defined and can be linked to function. an excellent model system for this purpose is ammonoid fossil shells. ammonoids were a successful and diverse group of species, which lived in the seas from to million years ago. ammonoid shells can be described by a morphospace defined by three geometrical parameters, defined in the pioneering work of david raup. in this morphospace, the outer shell is a logarithmic spiral, whose radius grows with each whorl by a factor w, the whorl expansion rate. there is a constant ratio between the inner and outer shell radii, denoted d. finally, the shell cross section can range from circular to elliptical, as described by s, the third parameter. raup’s w-d-s parameterization can be robustly measured from fossils although the coiling axis changes throughout ontogeny and thus, the coiling axis is sometimes difficult to exactly locate in actual specimens. it has been the setting for extensive research on ammonoid morphology and evolution, as well as the morphology of other shelled organisms.figure raup morphospace coordinates. ammonoid shell morphology can be described by three dimensionless geometrical parameters: w, the whorl expansion rate, is defined by a/b in the figure. d, the internal to external shell ratio, is x/a. s, the opening shape parameter, is y/z. the shell diameter can also be related to the parameters in this figure as shown. plotting each genus of ammonoids as a point in this morphospace, ignoring coiling axis changes, raup discovered that most of the theoretical morphospace is empty: many possible shell forms are not found. the existing forms lie in a roughly triangular region in the w-d plane. one reason for this distribution is geometric constraints. researchers have suggested that ammonoids tend to lie to the left of the hyperbola w = 1/d, because beyond this curve shells are gyroconic. such gyroconic shells are mechanically weaker and less hydrodynamically favorable. it is noteworthy, however, that shells to the right of the curve do exist in nature, for example in the bactritida or orthocerida lineages, which are probably ancestral to the ammonoids, as well as in heteromorph ammonoids that occasionally occur in the mesozoic and more commonly in the cretaceous. thus the w = 1/d curve is unlikely to be an absolute geometric constraint.figure ammonoids repeatedly filled the same triangle in d-w plane after mass extinctions. all of the ammonoid data used in the present study. red points are genera before the ff mass extinction, genera after ff are denoted by blue points. the green curve is w = 1/d. ammonoids before the ff extinction, together with a schematic arrow for the direction of evolution from ancestral taxa. genera between ff and dm mass extinctions fill out a triangle, surviving genera from the ff mass extinction are denoted by red bold points. ammonoids between dm and pt mass extinctions fill a triangle, surviving genera from the dm mass extinction are denoted by red bold points. ammonoids after the pt mass extinction fill a triangle, surviving genera from the pt mass extinction are denoted by red bold points. ammonoids from different periods, together, genera between ff and dm are denoted blue, dm to pt in red and post pt in green. the shell morphologies of the three archetypes at the vertices of the triangle are shown. studies in recent years have considered a larger dataset of ammonoids than raup. work, saunders and nikolaeva show that after each mass extinction, ammonoid genera refill the same roughly triangular morphospace. the repeated convergence to the same suite of variation raises the question of the relation between ammonoid morphology and function. most studies hypothesize a fitness function, which has an optimum in the middle of the triangular region. the fitness function is often taken to be dominated by hydrodynamic drag; this assumption is compelling since the contours of hydrodynamic efficiency, experimentally measured by chamberlain, show peaks at positions close to the most densely occupied regions of morphospace. the ammonoid genera are assumed to also occupy the slopes that descend from the fitness peaks, until bounded by the geometric constrains. interestingly, raup did not espouse the idea of a single task dominating fitness, but rather noted that multiple tasks might be at play. in every niche, different tasks become important, leading to niche-dependent fitness functions with different maxima. the idea of multiple tasks was elegantly employed by westermann, who described ammonoid morphospace by mapping it to a triangle. at the vertices are three ‘end member’ morphologies which correspond to three lifestyles. each morphology is mapped to a point in the triangle, which is interpreted as portraying the relative distance from the end members and hence the relative weights of the three lifestyles. the westermann morphospace was useful in comparing different datasets and in interpreting ammonoid lifestyles. the main drawback of the westermann morphospace is that, because it involves nonlinear dimensionality reduction, different morphologies can be mapped to the same point, and in some cases slight differences in shape can lead to large differences in the westermann projection. thus, it is of interest to seek a relation between shape and tasks without such drawbacks.figure the performance contours of the three putative tasks for ammonoid shells. contours for shell economy, defined as the ratio of internal volume to shell volume, with red denoting high values, and blue low values. for gyroconic shells, this performance function becomes constant, and equal to the lowest contour shown. the triangle encapsulating the entire ammonoid dataset is shown in black. contours for the drag coefficient measured by chamberlain, red lines denote lower drag or better hydrodynamics. contours for the growth function defined in the main text, red lines denote quicker growth. the contours of the three tasks give rise to a suite of variation denoted by blue points. to address this, we explore evolutionary tradeoffs between tasks in the framework of pareto optimality theory, to quantitatively explain the suite of variation in direct morphospace, and to infer the putative tasks at play. we find that ammonoid morphology in the w-d-s morphospace falls within a square pyramid, suggesting five tasks. the triangular region observed by raup is the projection of this pyramid on the w-d plane, and the westermann morphospace is a dimensionality reduction of the three-dimensional pyramid to a two-dimensional triangle. we propose putative tasks whose performance contours jointly lead to the observed suite of variations, including hydrodynamic efficiency, shell economy, compactness and rapid shell growth. the position of each species in this pyramid, namely its distance from each vertex, indicates the relative importance of each task in the niche in which that species evolved. after the ff and dm mass extinctions, surviving ammonoids refill essentially the same pyramid. after the pt extinction, part of the pyramid is refilled. these findings lend support to the pareto theory of evolutionary tradeoffs in the context of evolution on geological timescales.  ammonoid distributions in the w-d plane converge to a similar triangle after major extinctions we begin by considering ammonoid morphology in the w-d plane, and later consider the three dimensional w-d-s space. we combine the data of saunders, work and nikolaeva for paleozoic ammonoids, with the data of mcgowan for mesozoic ammonoids. the data is classified into three sets between mass extinctions: from ff to dm, from dm to pt, and after pt. we tested whether the ammonoid distribution in each set falls in a triangle more closely than randomized data, based on the statistical test of. we use an archetype analysis algorithm to find triangles, which enclose as much of the data as possible. we find that a triangle describes each dataset much better than randomized datasets in which the w and d coordinates are randomly permuted. randomized data rarely fill out a triangle as well as the real data. we next tested how similar the triangles are for the three datasets. we computed the ratio between the intersection area of the triangles to the union area as a measure for triangle similarity. the three triangles show large ratios of intersection to union area, and pairs respectively, p <10- compared to randomly generated triangles, see methods), indicating that the triangles are very similar. we conclude that after each extinction, ammonoids re-populate essentially the same triangular region. the vertices of the triangle describing the joint dataset of ammonoids after ff are d1*,w1*=, ~, d2*,w2*=, ~ d3*,w3*=, ~. we next ask which tasks might relate to each of the vertices. economy of shell material may determine the first archetype raup suggested that a possible need of the ammonoids is to maximize their internal volume relative to shell volume. this is important if shell production is costly, and also in terms of buoyancy considerations. ammonoids are thought to achieve neutral buoyancy by balancing shell weight with buoyancy from their air-filled chambers; high internal volume relative to shell material extends the range over which neutral buoyancy can be reached. to calculate shell material relative to internal volume at each point in morphospace, we follow raup and assume that shell thickness is a fixed fraction of radius, namely thickness/radius = , as measured by. interestingly, this ratio is close to the optimal ratio of thickness/radius = from calculations of mechanical strength in tube-like bones. we improve slightly on chamberlin and raup’s original calculation by numerically evaluating the necessary integrals rather than using the analytical approximations of, yielding corrections of about 10%. the maximum of internal volume relative to shell thickness occurs at  =. this point is close to archetype one d1*,w1*=,. the calculated contours of internal volume relative to shell thickness-namely the performance contours of the task of economy- have a curving ridge that points towards the third archetype. performance drops sharply on either side of this ridge.figure the three dimensional pareto front of the ammonoid dataset. the rms error for pcha optimal polygons and polyhedra is computed for different numbers of possible vertices: line, triangle, tetrahedron, 5-vertex polyhedron, etc. error decreases with increasing the number of archetypes up to increasing the number beyond doesn't improve the fit by much. the best fit 5-archetype polygon resembles a square pyramid. blue points denote ff to dm ammonoids, red are dm to pt and green are post pt ammonoids. archetypes are numbered, their morphology is shown, and the suggested tasks are listed in panel a. the second archetype may optimize hydrodynamics we conjecture that the second archetype maximizes the hydrodynamic efficiency of the ammonoids. low drag is important for ammonoids in order to swim rapidly. hydrodynamic efficiency is measured by the drag coefficient, which is a dimensionless number specific to each geometrical shape. the drag coefficient is proportional to the force which should be applied in order to keep an object of a given surface area moving at a given velocity in water. drag coefficients were measured by chamberlain using plexiglass models of shells. the contours of hydrodynamic efficiency are shown in figure 4b. drag monotonically increases with d and w, hence we can conclude that the ammonoid morphology with minimal drag has the lowest possible values of d and w, namely = this is close to the second vertex of the triangle, archetype two at d2*,w2*=,. the third archetype may optimize rapid shell growth the remaining vertex of the triangle, archetype has a large value of w. thus the shell radius at this vertex increases rapidly with each revolution of the spiral. there are different possible tasks that might relate to large w, including rapid growth, shell-orientation and swimming capabilities. in westermann morphospace, large w compared to d and s is interpreted as nektonic. here, we wish to demonstrate an essential approach, and thus focus on one of these potential tasks: rapid growth, and leave other possibilities to future study. the fossil dataset we use does not contain information on growth. however, if we assume that the ability to generate shell material is proportional to body mass we can predict the growth, or at least a function proportional to growth, using only the dimensionless parameters we have. an evolute shell allows volume to grow rapidly with each whorl. rapid growth may be important because predation tends to decrease with organism size. this would select for increased w. however, the whorl expansion rate w cannot grow without bound in order to avoid cyrtoconic shells- the shell must close over itself at least once to provide space for the ammonoid body. a coiled shell is also important in order to benefit from increased shell thickness, because until the ammonoid is closed, the thinner shell is exposed to the outside. the small value of d is also reasonable for such a function, because when w is large, a small d is a must in order to benefit from the advantages of w < 1/d. a similar function was suggested in snails where shell growth rate was found to be larger in snails in the presence of predators. to be concrete, we consider a putative performance function that penalizes the ammonoid for the smallness of its diameter, namely p3=∫0∞1diamtdt. contours of this performance function are shown in figure 4c. the function peaks at = close to the third archetype d3*,w3*=. at this archetype, ammonoids reach large diameters most rapidly. one may ask if the advantage of growth comes from the increased diameter which might make the ammonoid too large for specific predators, or from the increased shell thickness which make it stronger. it is difficult to distinguish between this two conjectures since from we know that this quantities are proportional to one another. it is likely that both diameter and shell thickness contribute to fitness. the three putative performance functions, shell economy, hydrodynamic efficiency, and shell growth together give rise to a triangular shaped pareto front. the pareto front boundaries are given by the points of tangency of the contours of the different performance functions. figure 4d shows the computed pareto front, which resembles a slightly curved triangle, and is similar to the observed suite of variation. ammonoid data is enclosed by a pyramid in w-d-s morphospace up to now, we considered ammonoid morphology projected on the w-d plane. we now turn to the analysis of the data in the three-dimensional morphospace, given by w,d and s—whorl expansion, radii ratio and the shape of the shell opening. low values of the parameter s correspond to oblate elliptical openings, giving rise to compressed shells. an s value of corresponds to a circular shell opening; high values corresponding to depressed shell morphologies.figure size is enriched at some of the archetypes. ammonoid shell diameter as a function of distance from each archetype shows that small diameters are prevalent near archetypes and data includes diameter for genera, divided into bins with equal number of genera according to the distance from each archetype. average diameter for each bin is plotted. for convenience, a fit of the averages to a line is shown. no diameter enrichment near archetype. positive diameter enrichment near archetype positive diameter enrichment near archetype. negative diameter enrichment near archetype negative diameter enrichment near archetype. we attempted to enclose the 3d dataset by polygons with to vertices. we evaluated the extent to which each polygon explains the data, by calculating the rms distance of points outside the polyhedron. we find that beyond vertices, the rms error does not decrease significantly: shapes with or more vertices do not improve the closeness of fit appreciably. hence a 5-vertex polygon is a parsimonious description of the data. this 5-vertex shape has four vertices that lie approximately on a plane. we thus consider this shape as a pyramid. a square pyramid encloses the data better than randomly permuted dataset with p <10-. the five vertices of the pyramid suggest five archetypes, whose coordinates are given in table  the square base of the pyramid has two vertices at low s, and two others, which match them for w and d values, but have higher s values. the apex of the pyramid has a thin opening with s = .table summary of suggested ammonoid archetypes suggested task w s d archetype number coordinates of the archetypes found by the archetype analysis algorithms with 5-archetypes, along with their putative tasks. projecting the pyramid on the w-d plane, we find that the apex of the pyramid matches the ‘growth’ archetype described above; the ‘economy’ and ‘hydrodynamic’ archetypes each corresponds to the projection of two 3d archetypes: the economy archetype corresponds to archetypes and and the hydrodynamic archetype to archetypes and. economy, hydrodynamic and growth performance functions are maximized near three of the pyramid vertices we repeated the calculation of economy performance in three dimensions. the 2d contours shown previously were evaluated at s =  by varying s, we find that the maximal economy is found at  =. this is reasonably close to vertex of the observed pyramid d1*,s1*,w1*=,,.the internal volume to shell volume ratio in this vertex is 96% of the optimum value. for comparison, this ratio drops to nearly zero near vertices and of the pyramid. the hydrodynamic efficiency measured in includes data at values of s other than s =  this indicates that optimal hydrodynamic efficiency is at low s values, i.e. s→ the resulting optimum is thus close to vertex of the pyramid, which is d2*,s2*,w2*=,,: note the low values of d,s and w. archetype has an s value close to. the dependence of the growth performance function on s comes only implicitly through the volume-to-surface ratio. it is unclear from the present simplified model for the growth performance function why is selected as the optimal s value for archetype this s value might be due, for example, to diminishing returns of shell production per body mass. in other words, the assumption that shell material production is constantly proportional to body mass might be imprecise. if shell production grows slower than linearly with body mass, this will favor smaller-volume ammonoids with smaller value of s that will increase diameter faster. the last two pyramid archetypes may be related to size two pyramid vertices remain to be explained, vertices and these vertices have large values of s, and correspond to depressed shells. we find that these shapes have the smallest ratio of surface area to volume. they are therefore the most globular in the suite of variation, in the sense that their height is most similar to their width and depth. one feature of globular ammonoids is small size for a given internal volume, because spherical shapes have the minimal diameter of all shapes with the same volume. up to now, we did not consider the absolute size of the ammonoids, only on dimensionless shape traits w, d and s. to address this, we correlated data by mcgowan on ammonoid size with distance from the five vertices of the pyramid. we find an enrichment of small ammonoids most strongly near archetypes and 5: the genera nearest to these vertices have the smallest diameters. archetypes and are enriched with large ammonoids and archetype has weak enrichment since its s value is relatively larger than archetypes and. archetypes and may thus correspond to economy and hydrodynamic tasks respectively, combined with a need for smallness. this relation between diameter and globularity is in line also with, which used a different dataset.figure archetypes in trait space and in westermann space. ammonoids in westermann space. the red ellipses are the pyramid archetypes projected on the westermann space, with 5% error in the archetype positions. three of the pyramid archetypes lies near the vertices of the triangle, another archetype is near the edge because both d and s are dominant. the region near archetype is severely warped because d,s and w are all relatively small. the pyramid in the w-d-s trait space, the ellipsoids are 5% errorbars around each vertex. these small ellipsoids translates to the red ellipses of subfigure a when switching to westermann morphospace. we further compared the way ammonoids from different periods fill out the pyramid. the main difference between periods is between paleozoic and mesozoic genera. mesozoic ammonoids tend to have lower s values than paleozoic ones, as found by mcgowan. in the pyramid, they are more densely arrayed near the face defined by vertices and and away from and especially from this may be interpreted in the present framework as a shift in the niches occupied by later ammonoids, in which tasks corresponding to archetypes and contribute less to fitness than they did in the paleozoic niches. finally, we mapped the five archetypes of the pyramid to the westremann morphospace. we find that that three archetypes, and map near the three vertices of the westermann triangle. the two other vertices of the pyramid map closer to the edges of the triangle. some of the archetypes map slightly outside of the triangle since they are exptrapolated points which lie outside of the ammonite dataset. we also asked about the sensitivity of this transformation, by testing a small region around each archetype. we find that one of the archetypes, archetype lies in a region of morphospace which is severely warped by the westermann transformation, and maps to a wide region in the triangle. the other archetypes are less sensitive and map to relatively small regions of the triangle. discussion this study explored how tradeoffs between multiple tasks may have contributed to the evolution of ammonoid shell morphology. ammonoid shell data on genera were studied in raup’s three parameter morphospace. the data is well described by a square pyramid. this finding is interpreted in light of pareto theory on tradeoffs between tasks. the five vertices of the pyramid may be interpreted as archetype morphologies optimal for a single task, and morphologies in the middle of the pyramid are generalists which compromise between the tasks. we propose candidate tasks for the archetypes. hydrodynamic efficiency is a good candidate for one of these tasks, and is maximized near vertex of the pyramid. other putative tasks can be inferred from the position of the other vertices of the pyramidal shell distribution. we propose that economy of shell material is a second task, quantitated by the ratio of internal volume and total shell material. the maximum of this function matches one of the vertices of the pyramid. a third task may be rapid growth. a performance function relating to rapid growth of ammonoid diameter is maximal near the apex of the pyramid, at shells with high w. two other tasks may relate to small spherical-like shells combined with low drag and high economy. it is interesting to relate this study with previous work by westermann and ritterbush based on the idea that ammonoids face tradeoffs between different tasks, which determine their morphologies. westermann proposed a morphospace which, instead of working in d-w-s space, works in a 2-dimensional projection which consists of ratios of related measurements. westermann morphospace has many advantages. as a method to reduce 3-dimensional data into 2-dimensional one, it helps visualize data in order to achieve better understanding of the geometry. it is also useful in understanding the different niches that ammonoids occupy and infer the various tasks they face. westermann's 2-dimensional representation also has drawbacks. as a dimensionality reduction method, it loses information about the data. ammonoid shells with very different geometries can be mapped to the same point in westermann morphospace. moreover, because the westemann map is nonlinear, there are regions in morphospace that map to the triangle with relatively large errors. for example, a small region around the point of minimal values of d, s and w can mapped to the entire westermann triangle depending on slight variation in the values. the present approach does not show these drawbacks because it works directly in w-d-s morphospace. it thus distinguishes between morphologies which are mapped to the same westermann point. the square pyramid identified here suggests two new tasks in addition to westermanns three. these are archetypes and which correspond to hydrodynamic drag and compact shell economy. we also propose a different interpretation of the other three tasks. for example the westermann oxycone endmember, which is linked to nektonic lifestyle, corresponds to our archetype. furthermore, the same oxyconic endmember relates to certain morphologies near our archetype. the present study also bears on the question of geometric constraints in evolution. the w = 1/d line in ammonoids is thought to be an outstanding example of a geometric constraint, because of the disadvantages of the open shell morphology beyond this curve. this assumption is challenged by the existence of organisms with w > 1/d, including lineages ancestral to ammonoids as well as several heteromorphs. the present approach can make the concept of geometric constraint more precise by relating it to biological tasks. we consider the performance functions of tasks, some of which indeed show a decline beyond the w = 1/d line. in particular, economy and hydrodynamics contours both begin to sharply decline when w > 1/d. this provides a more principled explanation, replacing strict geometric constraint with the more subtle dependence of specific performance functions on geometry. other taxa may perform a different set of tasks, including a task with an archetype in the ‘forbidden for ammonoid’ region, w > 1/d. such tasks might explain the morphology of the taxa which show gyroconic shells. an alternative view is that some characters states do not require a functional explanation, but rather were neutral enough for a clade to succeed for some time. this study adds to previous studies that used the pareto approach to analyze other biological systems. these systems showed lines, triangles or tetrahedra in morphospace. ammonoids are the first system in which a pyramidal pareto front is observed. for this purpose, we find that the archetype analysis algorithm pcha is an efficient way to detect high order polyhedra in data. the present approach can be readily extended to other shelled organisms such as gastropods and bivalves. one application of the present approach is a quantitative inference of which task is important for fitness in the particular niche of each genus. the closer the shell morphology is to a given vertex of the pyramid, the more important the corresponding task. since ammonoid shells are carried by currents and found in rocks far from the habitat of the living organism, it is challenging to connect morphology with behavior. the present approach can offer quantitative inference about the relative contribution of tasks to fitness, to provide insight into the ecological niche of these extinct organisms. more generally, this study supports basic predictions of the pareto theory for evolutionary tradeoffs, which we hope will be useful also for other biological contexts.  this study supports fundamental predictions of the pareto theory of tradeoffs by shoval et al. that have not been previously tested on the scale of hundreds of millions of years of evolution. ammonoid shell data on genera is well-described by a square pyramid in morphospace. the five vertices of the pyramid may be interpreted as archetype morphologies optimal for a single task. inferred tasks include shell economy, rapid growth, compactness and hydrodynamic efficiency. the vertices of the polygon do not change over the timescale of interest, as predicted in the case where the tasks stay more-or-less constant because they are dictated by biomechanical considerations. moreover, the polygons and polyhedra in the theory are not necessarily due to phylogenetic history, but rather to convergent evolution to pareto-optimal solutions. this agrees with the finding that after a mass extinction which removes almost all of the species, survivor species evolve to re-fill the same polygon as their ancestors. this approach may be used to infer biological tasks from data in other biological contexts. 
their evolutionary history includes repeated mass extinctions. after mass extinctions, surviving species evolve to refill essentially the same pyramid, suggesting that the tasks are unchanging. organisms that need to perform multiple tasks face a fundamental tradeoff: no design can be optimal at all tasks at once. ammonoid fossil shells provide an excellent model system for this purpose. we infer putative tasks for each archetype, related to economy of shell material, rapid shell growth, hydrodynamics and compactness. these results support pareto optimality theory as an approach to study evolutionary tradeoffs, and demonstrate how this approach can be used to infer the putative tasks that may shape the natural selection of phenotypes. recent theory based on pareto optimality showed that such tradeoffs lead to a highly defined range of phenotypes, which lie in low-dimensional polyhedra in the space of traits.
 organisms that need to perform multiple tasks face a fundamental tradeoff: no design can be optimal at all tasks at once. recent theory based on pareto optimality showed that such tradeoffs lead to a highly defined range of phenotypes, which lie in low-dimensional polyhedra in the space of traits. the vertices of these polyhedra are called archetypes- the phenotypes that are optimal at a single task. to rigorously test this theory requires measurements of thousands of species over hundreds of millions of years of evolution. ammonoid fossil shells provide an excellent model system for this purpose. ammonoids have a well-defined geometry that can be parameterized using three dimensionless features of their logarithmic-spiral-shaped shells. their evolutionary history includes repeated mass extinctions.  we find that ammonoids fill out a pyramid in morphospace, suggesting five specific tasks - one for each vertex of the pyramid. after mass extinctions, surviving species evolve to refill essentially the same pyramid, suggesting that the tasks are unchanging. we infer putative tasks for each archetype, related to economy of shell material, rapid shell growth, hydrodynamics and compactness.  these results support pareto optimality theory as an approach to study evolutionary tradeoffs, and demonstrate how this approach can be used to infer the putative tasks that may shape the natural selection of phenotypes. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords multi-objective optimalityrepeated evolutionpareto frontdiversityperformancegoalissue-copyright-statement© the author 2015 