# Tools for quantitative form description; an evaluation of different software packages for semi-landmark analysis

Léo Botton-Divet, Alexandra Houssaye, Anthony Herrel, Anne-Claire Fabre, Raphael Cornette

► To cite this version:

Léo Botton-Divet, Alexandra Houssaye, Anthony Herrel, Anne-Claire Fabre, Raphael Cornette. Tools for quantitative form description; an evaluation of different software packages for semi-landmark analysis. PeerJ, 2015, 3, pp.e1417. 10.7717/peerj.1417. hal-01277381

HAL Id: hal-01277381

https://hal.sorbonne-universite.fr/hal-01277381

Submitted on 22 Feb 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Tools for quantitative form description; an evaluation of different software packages for semi-landmark analysis

Léo Botton-Divet<fnref n="1" />,<fnref n="2" />, Alexandra Houssaye<fnref n="1" />, Anthony Herrel<fnref n="1" />,<fnref n="3" />, Anne-Claire Fabre<fnref n="4" /> and Raphael Cornette<fnref n="5" />

## ABSTRACT

The challenging complexity of biological structures has led to the development of several methods for quantitative analyses of form. Bones are shaped by the interaction of historical (phylogenetic), structural, and functional constrains. Consequently, bone shape has been investigated intensively in an evolutionary context. Geometric morphometric approaches allow the description of the shape of an object in all of its biological complexity. However, when biological objects present only few anatomical landmarks, sliding semi-landmarks may provide good descriptors of shape. The sliding procedure, mandatory for sliding semi-landmarks, requires several steps that may be time-consuming. We here compare the time required by two different software packages ('Edgewarp' and 'Morpho') for the same sliding task, and investigate potential differences in the results and biological interpretation. 'Morpho' is much faster than 'Edgewarp,' notably as a result of the greater computational power of the 'Morpho' software routines and the complexity of the 'Edgewarp' workflow. Morphospaces obtained using both software packages are similar and provide a consistent description of the biological variability. The principal differences between the two software packages are observed in areas characterized by abrupt changes in the bone topography. In summary, both software packages perform equally well in terms of the description of biological structures, yet differ in the simplicity of the workflow and time needed to perform the analyses.

**Subjects** Zoology, Anatomy and Physiology, Computational Science

**Keywords** Geometric morphometrics, Sliding semi-landmark, Software comparison

## INTRODUCTION

Because the interaction of form and function impacts performance, biological shape is under direct selection. The shape of a bone is influenced by several parameters. The evolutionary history of an organism plays an important role and often a strong phylogenetic signal is detected in bone shape (Morgan, 2009; Fabre et al., 2013a;

<footnote n="1">¹ UMR 7179, Mécadev, Muséum National d'Histoire Naturelle, Centre National de la Recherche Scientifique, Paris, France</footnote>

<footnote n="2">² Université Denis Diderot (Paris VII), Paris, France</footnote>

<footnote n="3">³ Evolutionary Morphology of Vertebrates, Ghent University, Ghent, Belgium</footnote>

<footnote n="4">⁴ Animal Locomotion Laboratory, Department of Evolutionary Anthropology, Duke University, Durham, NC, USA</footnote>

<footnote n="5">⁵ UMR 7205, Institut de Systématique, Évolution, Biodiversité, Centre National de la Recherche Scientifique, Muséum National d'Histoire Naturelle, Université Pierre et Marie Curie, Ecole Publique des Hautes Études, Paris, France</footnote>
Fabre et al., 2014a; Álvarez, Ercoli & Prevosti, 2013). Moreover, structural factors can shape bones as physical laws constrain potential shapes and their morphogenesis (e.g., Cubo, 2004; Cubo et al., 2008). For example, functional constraints imposed through differences in locomotor behavior (e.g., Fabre et al., 2013b; Álvarez, Ercoli & Prevosti, 2013) and prey capture strategy or prey size (Andersson, 2004; Meloro et al., 2008; Cornette et al., 2013) can affect bone size and shape. Additionally, mechanical loads during development (e.g., Beaupre, Orr & Carter, 1990; Bass et al., 2002) and during the lifetime of an individual (Lanyon et al., 1982; Currey, 2003) can shape bones because of bone remodeling.

Morphometric approaches have been used extensively for quantifying shape variation in biological objects including bones (Adams, Rohlf & Slice, 2004). The first paradigm, often referred to as ‘traditional morphometrics’ was based on a statistical treatment of distances or ratios (Marcus, 1990). Yet, this approach is limited as the information captured by linear measurements does not describe the whole geometry. As a consequence, the visualization of the geometry subsequent to statistical analysis is not possible and objects with different shapes can theoretically give rise to similar measurements. Additionally, many structures are difficult to measure in practice, leading to the use of measures such as maximal length which are functionally but not anatomically homologous.

These limitations have driven the development of geometric morphometrics in the 1990's (Bookstein, 1991; Dryden & Mardia, 1993; Adams, Rohlf & Slice, 2004; Zelditch, Swiderski & Sheets, 2012). Geometric morphometrics use outline descriptors (Fourier analysis) or point coordinates (called landmarks) rather than distances to describe the geometry of biological objects. In this approach, size is explicitly defined (centroid size, sum of the squared distances of each landmark to the centroid; Bookstein, 1991) and the relative position of landmarks is conserved (Adams, Rohlf & Slice, 2004; Zelditch, Swiderski & Sheets, 2012).

Since a big part of biological variability cannot be assessed by using anatomical landmarks (biologically homologous landmarks) only, sliding semi-landmarks were developed to quantify complex shapes devoid of landmarks. Sliding semi-landmarks can be placed on curves (Bookstein, 1997; Gunz, Mitteroecker & Bookstein, 2005) and surfaces (Gunz, Mitteroecker & Bookstein, 2005). Due to the impossibility to define anatomically homologous points on curves and surfaces this approach generates landmarks that are spatially homologous after sliding (Parr et al., 2012). Sliding semi-landmarks are allowed to move on curves and surfaces in order to optimize a pre-defined criterion (Gunz & Mitteroecker, 2013). However, the choice of this criterion has been subject to debate. The most commonly used criterion is the bending energy (Skinner & Gunz, 2010; Pizzo et al., 2011; Cornette et al., 2013; Fabre et al., 2014b). The alternative involves the minimization of the Procrustes distance (Perez, Bernal & Gonzalez, 2006; Gunz & Mitteroecker, 2013). Sliding semi-landmarks are particularly well suited for the study of bones, providing descriptors of crests or outlines (curve sliding semi-landmarks; Morgan, 2009; De Groote, Lockwood & Aiello, 2010; Monteiro & Nogueira, 2010; Álvarez, Ercoli & Prevosti, 2013) and surfaces such as articular surfaces and the diaphysis of long bones (surface sliding semi-landmarks; Fabre et al., 2013a; Fabre et al., 2013b; Fabre et al., 2014b; Schlager, 2013a; Cornette et al., 2013; Morita et al., 2014; Cornette, Tresset & Herrel, 2014).
Several software packages are currently available that allow one to perform the sliding procedure in three dimensions. One of the first packages that was developed and that is still used in many studies (e.g., *Kulemeyer et al., 2009*; *Fabre et al., 2013a*; *Fabre et al., 2013b*; *Cornette et al., 2013*; *Cornette, Tresset & Herrel, 2014*) is 'Edgewarp' (*Bookstein & Green, 1994*). The 'EVAN toolbox' (http://evan-society.org) also performs sliding of semi-landmarks amongst other operations. Some authors have developed their own routines such as the Mathematica (Wolfram Inc., Modesto, California, USA) code used by *Gunz & Mitteroecker (2013)*. Recently two R packages have been published: geomorph (*Adams & Otárola-Castillo, 2013*), featuring landmark placement, treatment and analysis, and 'Morpho' (*Schlager, 2013b*) featuring landmark importation from several other software packages and performing geometric morphometric treatment and analysis.

Because specimen digitization, the first step of any three dimensional surface analysis, is a time-consuming step that cannot be easily shortened, it could be of great interest to reduce the duration of the second step: the sliding procedure. The aim of the present study is to compare the workflow and the results obtained with two different software packages for the same three dimensional sliding task. We selected 'Edgewarp' (*Bookstein & Green, 1994*), an established reference for these types of analyses and 'Morpho' (*Schlager, 2013b*) a recently published set of R routines. We focused on a practical and biological test comparing a long bone in different mustelids. We followed the sliding procedure detailed in *Gunz, Mitteroecker & Bookstein (2005)* using both software packages. Differences between the results were examined and discussed in the light of the biological structures involved. The time required was measured for several steps and the global workflow (e.g., file handling, external software requirements) was compared.

# MATERIALS AND METHODS

## Material

We used the humeri from 10 specimens belonging to five mustelid species (2 specimens per species; Table 1): *Meles meles* (Linnaeus, 1758), *Mustela putorius* (Linnaeus, 1758), *Gulo gulo* (Linnaeus, 1758), *Martes martes* (Linnaeus, 1758), and *Enhydra lutris* (Linnaeus, 1758). These species were chosen because they are widely distributed across the mustelid phylogeny and illustrate a relatively large diversity in both size, shape and ecology (*Nowak, 2005*; *Schutz & Guralnick, 2007*; *Wilson et al., 2009*; *Hunter & Barrett, 2011*). By choosing these specimens we attempted to encounter the largest number of potential pitfalls and difficulties during the sliding procedure. All specimens used are housed in the collections of 'Mammiferes et Oiseaux' from the Muséum National d'Histoire Naturelle in Paris, France.

## 3D modeling

Bones were digitized in three dimensions using a white light fringe Breuckmann 3D surface scanner (white light fringe StereoSCAN³D model with a camera resolution of five megapixels). Raw scans were treated using Geomagic (Geomagic Studio; Raindrop Geomagic, Research Triangle Park, North Carolina, USA) to fill-in remaining holes and to remove highly creased edges and spikes. Next, models were decimated to contain 100,000 triangles
<table><caption>**Table 1** Specimens used in this study.</caption><thead><tr><th>Species</th><th>Institution specimen</th></tr></thead><tbody><tr><td rowspan="2">Martes martes</td><td>CG 1994-806</td></tr><tr><td>CG 2005-232</td></tr><tr><td rowspan="2">Gulo gulo</td><td>CG 1983-946</td></tr><tr><td>CG 1995-1208</td></tr><tr><td rowspan="2">Mustela putorius</td><td>CG 1991-605</td></tr><tr><td>CG 2004-639</td></tr><tr><td rowspan="2">Meles meles</td><td>CG 2005-707</td></tr><tr><td>CG 1987-28</td></tr><tr><td rowspan="2">Enhydra lutris</td><td>CG 1935-124</td></tr><tr><td>A12503</td></tr></tbody></table>

**Notes.**

Institutional abbreviations are as follows: CG, Muséum National d'Histoire Naturelle Catalogue Général, Paris, France; A, Muséum National d'Histoire Naturelle Anatomie Comparée, Paris, France.

resulting in homogeneous 3D surface models. The number of triangles was arbitrarily fixed to generate models that are not too cumbersome, but without altering the geometry of the object. The models were converted to ‘.sur’ using the ‘obj2sur’ tool provided at ‘ftp://brainmap.stat.washington.edu/edgewarp/utils/’ in order to import them into ‘Edgewarp.’

## 3D anatomical landmark and curve digitization

Landmarks and curves were digitized on the surfaces of the scans using the Landmark software package (Wiley et al., 2005). Twenty-seven 3D homologous anatomical landmarks were chosen and are visible on all specimens (Fig. 1 and Table 2). Eighteen 3D curves were defined at the margins of articular surfaces and along crests. All curves are bordered by anatomical landmarks as recommended by Gunz, Mitteroecker & Bookstein (2005). The curves were digitized with a high density of points (40–160 points per curve depending on the curve length) and then sub-sampled to the number listed in Table 3. This dataset was then used to perform an identical sliding procedure using both ‘Edgewarp’ and ‘Morpho.’

A semi-automatic point placement was used to place surface semi-landmarks on the scans. A template was modeled using the Blender software (Blender Online Community, 2014) following the procedure described in Souter et al. (2010). We created a mesh of 598 points representing a simplified form of the humerus. Landmarks and curves were digitized on the template surface as it was done on the actual scans of the specimens. Then in both ‘Edgewarp’ and ‘Morpho,’ landmarks and curves were used to compute a thin plate spline (TPS) deformation of the template. Next, surface sliding landmarks were projected from the deformed template onto the bone surface (Gunz & Mitteroecker, 2013). In summary, the template used in this study contains a total of 817 points including 27 3D homologous anatomical landmarks, 192 sliding semi-landmarks on curves and 598 sliding semi-landmarks on surfaces.

The sliding procedure was performed following the algorithm detailed in Gunz, Mitteroecker & Bookstein (2005). Four TPS relaxations were performed, the first TPS
<figure><img src="image_2.png" /><figcaption>**Figure 1** 3D view of a *Martes martes* MNHN 2005-232 humerus showing the location of the 27 anatomical landmarks and 18 curves used to quantify the humeral shape. (A) caudal; (B) cranial; (C) lateral views. Red: anatomical landmarks. Blue: curve sliding semi-landmarks. Yellow: surface sliding semi-landmarks. See Table 2 for landmark definitions.</figcaption></figure>

relaxation was performed against the template, the three others against Procrustes consensus calculated using the data from the previous iteration.

## Time estimation

The time for surface pre-processing and initial projection was measured by the value displayed in 'Edgewarp' log. The time required by 'Edgewarp' for the three iterations of sliding against Procrustes consensus was measured using a timer. The time required to
<table><caption>**Table 2 Definition of the anatomical landmarks.**</caption><thead><tr><td>LM</td><td>Definition</td></tr></thead><tbody><tr><td>1</td><td>Most disto-medial point of the trochlea</td></tr><tr><td>2</td><td>Most medio-proximal point of the caudal side of the trochlea</td></tr><tr><td>3</td><td>Point of maximum of curvature of the olecranon fossa</td></tr><tr><td>4</td><td>Most latero-proximal point of the caudal side of the trochlea</td></tr><tr><td>5</td><td>Most distal point of contact between the trochlea and the capitulum</td></tr><tr><td>6</td><td>Most latero-proximal point of the cranial side of the capitulum</td></tr><tr><td>7</td><td>Point of maximum of curvature of the radial fossa</td></tr><tr><td>8</td><td>Maximum concavity of the cranial margin of the trochlea</td></tr><tr><td>9</td><td>Point of maximum of curvature of the coronoid fossa</td></tr><tr><td>10</td><td>Most medio-proximal point of the cranial side of the trochlea</td></tr><tr><td>11</td><td>Most distal point of the cranial side of the supracondylar foramen</td></tr><tr><td>12</td><td>Most proximal side of the cranial side of the supracondylar foramen</td></tr><tr><td>13</td><td>Most distal tip of the medial epicondyle</td></tr><tr><td>14</td><td>Most proximal tip of the medial epicondyle</td></tr><tr><td>15</td><td>Most distal point of the caudal side of the supracondylar foramen</td></tr><tr><td>16</td><td>Most proximal point of the caudal side of the supracondylar foramen</td></tr><tr><td>17</td><td>Most disto-cranial point of the lateral epicondylar crest</td></tr><tr><td>18</td><td>Most proximal point of the lateral epicondylar crest</td></tr><tr><td>19</td><td>Most distal point of the deltopectoral crest</td></tr><tr><td>20</td><td>Upper tip of the lesser tuberosity</td></tr><tr><td>21</td><td>Most disto-medial point of the lesser tuberosity</td></tr><tr><td>22</td><td>Most medio-caudal point of contact between the lesser tuberosity and humeral head</td></tr><tr><td>23</td><td>Disto-caudal tip of the humeral head</td></tr><tr><td>24</td><td>Latero-caudal point of contact between the greater tuberosity and the humeral head</td></tr><tr><td>25</td><td>Most antero-proximal point of the greater tuberosity crest</td></tr><tr><td>26</td><td>Tip of the tuberosita teres minor</td></tr><tr><td>27</td><td>Contact point between tricipital line and greater trochanter crest</td></tr></tbody></table>

**Notes.**

LM, landmark index.

run the ‘Morpho’ package in R was assessed by running the timestamp R function at the beginning and end of tasks and calculating the difference. These durations do not include time required for writing the script (scripts for data formatting for ‘Edgewarp’, and R scripting for ‘Morpho’). All analyses were run on a laptop computer (Asus k55vj) with a Intel® i7-3630QM CPU, 4 Gb of memory, running on Linux Ubuntu V14.04.

## Geometric morphometrics and visualization

Data analysis and visualization were performed using the R software (R Core Team, 2014). In order to superimpose geometries and isolate size and shape, a generalized Procrustes analysis (GPA) was performed (Gower, 1975; Rohlf & Slice, 1990; Dryden & Mardia, 1998) for the ‘Edgewarp’ and the ‘Morpho’ dataset separately using the ‘Rmorph’ package (Baylac, 2013). We reduced dimensionality of the datasets by keeping the nine first non-null axes of a Principal Component Analysis (PCA) performed on the coordinates
<table><caption>**Table 3** Curve designation and number of sliding semi-landmarks per curve; for anatomical landmark designation please refer to Table 2 and Fig. 1.</caption><thead><tr><td>Init LM</td><td>Term LM</td><td>NbSl</td></tr></thead><tbody><tr><td>1</td><td>2</td><td>5</td></tr><tr><td>2</td><td>4</td><td>10</td></tr><tr><td>4</td><td>5</td><td>15</td></tr><tr><td>5</td><td>6</td><td>15</td></tr><tr><td>6</td><td>8</td><td>10</td></tr><tr><td>8</td><td>10</td><td>3</td></tr><tr><td>10</td><td>1</td><td>15</td></tr><tr><td>11</td><td>12</td><td>5</td></tr><tr><td>12</td><td>11</td><td>5</td></tr><tr><td>15</td><td>16</td><td>5</td></tr><tr><td>16</td><td>15</td><td>5</td></tr><tr><td>17</td><td>18</td><td>23</td></tr><tr><td>20</td><td>21</td><td>5</td></tr><tr><td>22</td><td>23</td><td>10</td></tr><tr><td>23</td><td>24</td><td>9</td></tr><tr><td>25</td><td>27</td><td>26</td></tr><tr><td>26</td><td>27</td><td>21</td></tr><tr><td>27</td><td>19</td><td>5</td></tr></tbody></table>

**Notes.**

Abbreviations are as follows: init LM, initial anatomical landmark index; term LM, terminal anatomical landmark index; NbSl, number of curve sliding semi-landmarks.

in the tangent space. In order to assess differences between the dataset produced by each software, a PROTEST analysis was performed using the ‘vegan’ package (Oksanen et al., 2015) run with 100,000 iterations (Peres-Neto & Jackson, 2001). The ‘vegan’ package was then used to draw the Procrustean superimposition plot. The three dimensional visualization of bone surfaces, landmarks, and vectors between homologous points were performed using the ‘rgl’ (Adler & Murdoch, 2012) and ‘Morpho’ (Schlager, 2013b) packages in R.

# RESULTS

## Global shape variability

The PROTEST analysis highly supports the congruence of the two datasets ($m_{12} = 0.00325$; $P < 10^{-5}$). The Procrustean superimposition plot (Fig. 2) shows the differences between the two datasets. It becomes immediately clear that the between-methods variability is smaller than the inter-specimens variability, even within species. The Procrustes residuals of specimens analyzed with both methods vary from one specimen to the other (Table 4). The specimens showing the maximal Procrustes residuals are *Gulo gulo* 1983-946, *Mustela putorius* 1991-605, and *Enhydra lutris* 1935-124 (respectively c, e, and i; Fig. 2). The first axis separates *Enhydra lutris* (i and j; Fig. 2) and *Martes martes* (a and b; Fig. 2). On the two first axes *Mustela putorius* and *Gulo gulo* show intra-specific distances greater than the inter-specific ones with *Mustela putorius* 2004-639 (f; Fig. 2) being closer to *Gulo gulo*
<figure><img src="image_3.png" /></figure>

**Figure 2** Procrustes superimposition of the Morpho (circles) and Edgewarp (arrow heads) results. Arrows present Procrustes residuals. (a) *Martes martes* 2005-232; (b) *Martes martes* 1994-806; (c) *Gulo gulo* 1983-946; (d) *Gulo gulo* 1995-1208; (e) *Mustela putorius* 1991-605; (f) *Mustela putorius* 2004-639; (g) *Meles meles* 2005-707; (h) *Meles meles* 1987-28; (i) *Enhydra lutris* 1935-124; (j) *Enhydra lutris* A12503.

1995-1208 (d; Fig. 2) than to *Mustela putorius* 1991-605 (e; Fig. 2) and *Gulo gulo* 1983-946 (c; Fig. 2) being closer to *Mustela putorius* 1991-605 (e; Fig. 2) than to *Gulo gulo* 1995-1208 (d; Fig. 2).

## Differences between mean shapes

The middle part of the greater trochanteric crest (2 in Fig. 3) presents the strongest differences. The neighboring points on the surface of the greater trochanteric crest are also affected. The second most variable area is the curve sliding around the cranial side of the medial supracondylar foramen (1 in Fig. 3) at the latero-distal part. Additionally, a point on the diaphysis located on the disto-medial side of the greater trochanter crest varies considerably.

## The most variable specimens

As highlighted in Fig. 2 and Table 4, differences linked to the software are more important in some specimens than in others. *Gulo gulo* 1983-946, *Mustela putorius* 1991-605, and
<table><caption>**Table 4** Procrustes residuals per specimen from the superimposition of the 9 first principal components of the Edgewarp and Morpho slid datasets.</caption><thead><tr><td>Specimen</td><td>Procrustes residuals</td></tr></thead><tbody><tr><td>Enhydra lutris 1935-124</td><td>0.020489693</td></tr><tr><td>Enhydra lutris A12503</td><td>0.017785627</td></tr><tr><td>Gulo gulo 1983-946</td><td>0.037204213</td></tr><tr><td>Gulo gulo 1995-1208</td><td>0.010347460</td></tr><tr><td>Martes martes 1994-806</td><td>0.009211963</td></tr><tr><td>Martes martes 2005-232</td><td>0.015334387</td></tr><tr><td>Meles meles 1987-28</td><td>0.008102059</td></tr><tr><td>Meles meles 2005-707</td><td>0.011307574</td></tr><tr><td>Mustela putorius 1991-605</td><td>0.021816036</td></tr><tr><td>Mustela putorius 2004-639</td><td>0.005777793</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 3** Warped mesh calculated on the Procrustes mean shape showing the most varying areas. (A) caudal; (B) lateral; (C) cranial; (D) medial views. Sphere colors (from green to red) and size are proportional to the Euclidean distance between points in the Procrustes mean shape for each software output after sliding. Points showing the maximal distance between the results obtained by the two software packages are large red points (as opposed to minimal distances being represented by small green points). 1: cranial side of the medial supracondylar foramen; 2: greater trochanteric crest.</figcaption></figure>
<figure><img src="image_5.png" /></figure>

**Figure 4** Close up view of the most variable areas of some specimens after sliding. Yellow: 'Edgewarp'; blue: 'Morpho'; black vectors link the homologous points in the output of each software package; the red vectors highlight the most variable points. (A) *Gulo gulo* 1983-946 medio-cranial view of the greater trochanteric crest; (B) *Mustela putorius* 1991-605 medio-caudal view of the lesser trochanter (distal) crest along the diaphysis. (C) *Enhydra lutris* 1935-124 medial view of the distal part of the greater trochanteric crest.

*Enhydra lutris* 1935-124 display the largest differences between methods (Table 4). Beyond the differences noted in the mean shapes, we investigated the ones specific to each of these specimens. *Gulo gulo* 1983-946 (Figs. 4A) presents a displacement of points located on the deltoid crest, more specifically around the insertion of an ossified tendon (teres major/latissimus dorsi). The most variable point switches from one side of this structure to the other depending on the software used, suggesting that one of the methods did not succeed in sliding the point across this structure. *Mustela putorius* 1991-605 (Fig. 4B) presents a point shifting from one side to the other side of the crest on the diaphysis distal to the lesser trochanter. The 'Edgewarp' points (yellow) are located on the caudal side of the crest whereas the 'Morpho' points (blue) are located on the cranial side. *Enhydra lutris* 1935-124 (Fig. 4C) presents a variation associated with a small depression located medially to the deltoid crest. 'Edgewarp' points (yellow) slide on the distal margin of the foramen whereas the 'Morpho' points (blue) are close to the cranial margin of the depression. Surrounding points are also affected. The most variable points around the supracondylar foramen, at the latero-distal part of its cranial side are driven by the *Meles meles* 1987-28. For this specimen 'Edgewarp' points are concentrated on the most proximal part of this curve whereas 'Morpho' points are located more distally.

## Time

The times measured for surface pre-processing, the initial point projection, and the relaxations against the Procrustes mean shape, iterated three times, are given in Table 5. 'Morpho' requires a much shorter time compared to 'Edgewarp'. The time required for the sliding against Procrustes mean shape by 'Edgewarp' is moreover an active time, meaning that
<table><caption>**Table 5** Comparison of the duration of several tasks of the sliding procedure on the whole dataset with both ‘Edgewarp’ and ‘Morpho.’</caption><thead><tr><td></td><td>Edgewarp</td><td>Morpho</td><td>%</td></tr></thead><tbody><tr><td>Surface pre-processing (s)</td><td>1,320</td><td>4</td><td>0.3</td></tr><tr><td>Initial point projection (s)</td><td>600</td><td>20</td><td>3.3</td></tr><tr><td>Sliding against mean shape 3 iterations (min)</td><td>105</td><td>2</td><td>1.9</td></tr></tbody></table>

**Notes.**

%, percentage of time required by ‘Morpho’ compared to ‘Edgewarp.’

the user needs to run each step of the iterations manually. This part also requires the computation of the GPA which is not implemented internally in ‘Edgewarp’. Based on the results of the GPA, new input files must be generated prior to performing the next iteration.

## DISCUSSION

### General variability

The distribution of specimens in morphological space is similar, irrespective of the software used. The biological information is similarly described with both packages and thus both can be used to describe complex biological shapes without an inherent bias. Nevertheless, as some local divergence between methods can be observed, results could differ for intra-specific studies often describing more subtle shape variations. The morphospace obtained after sliding is structured by the most extreme ecological niches in our sample suggesting that functional constraints drive much of the shape variability. The first dimension of the Procrustean superimposition plot separates *Enhydra lutris* and *Meles meles* (respectively i, j and g, h Fig. 2) from the other species. The maximal distance on this axis is obtained between *Enhydra lutris* (i and j, Fig. 2) which is the most aquatic of the Mustelidae sampled here (Estes, 1980; Bodkin, 2001) and *Martes martes* (a and b, Fig. 2) which is the most arboreal (Nowak, 2005) species in our sample. The second dimension of the Procrustean superimposition plot separates *Meles meles* (g and h, Fig. 2), the only semi-fossorial species in our sample (Nowak, 2005) from *Enhydra lutris* (i and j). Phylogeny, size and/or body mass (allometry) have been shown to be important factors shaping the humerus in Carnivora (*Heinrich & Biknevicius, 1998; Fabre et al., 2013a; Martín-Serra, Figueirido & Palmqvist, 2014*). Our data suggest that life-style is indeed an important driver of forelimb shape in mustelids. Nevertheless, this is a small sample only and these results should be interpreted with caution.

### Differences

As they are computed from the same dataset, morphospaces obtained by the two software where expected to be the same. The two morphospaces (relative distances between specimens) are indeed generally similar, which is validated by the significant PROTEST. Nevertheless some specimens show slightly different positions in the morphospace depending on the software (*Gulo gulo* 1983-946, *Mustela putorius* 1991-605, and *Enhydra lutris* 1935-124, respectively c, e, and i; Fig. 2). Specimens showing the greatest variation in their relative position between the two software are bones with areas where a structure
interrupts the main curvature of the surface resulting in strong local deformations (Fig. 4). The fact that these local deformations are of a size similar to the distance between surface sliding landmarks appears to play a role. Indeed, highly curved structures can “puncture” the meshing constituted by the surface sliding semi-landmarks. Even if concerning only few points, these structures can lead to displacement of the specimen into the morphospace and could consequently have an impact in the biological interpretation of the morphospace (i.e., modifying phenetic affinities; Fig. 2). A greater number of surface sliding semi-landmarks in these specific areas could allow points to slide at the surface of a structure. Additionally, if the aforementioned structures present no biological meaning for the study, an a priori elimination could avoid giving too much weight to these structures in the analysis. Small foramina and ossified tendon insertions could be removed prior to the sliding procedure using dedicated tools in meshing software.

The differences observed between the two datasets thus likely result from differences during the sliding process at particular points of some specimens rather than reflect a global difference between methods. Such differences could be a consequence of differences between the two software packages in the computation of the relaxation as well as in the computation of the projection of the points on the surface of the specimen. Consequently mixing in a single analysis results obtained by the two software should be avoided.

## Time

'Edgewarp' appears more time-consuming than 'Morpho' for the same sliding task. 'Edgewarp' requires 30-300 times the time required by 'Morpho' for the same tasks. Two reasons can be put forward to explain these differences: computational efficiency and workflow complexity. 'Edgewarp' is a non-parallelized software whereas many 'Morpho' functions can be run on several cores, thus improving processing time. The effect of computation power remains low on a small dataset like the one presented here, but can lead to a considerable time difference in the case of the analysis of larger datasets. In addition to the number of steps constituting the workflow, the number of actions performed by the user should also be taken into consideration. The use of 'Edgewarp' requires one to open each surface, curve, and landmark file manually. The subsequent iterative relaxation requires several manual operations per iteration as does the data saving after sliding. The use of a programming language such as python (Van Rossum & De Boer, 1991), bash (Free Software Foundation, 2013), or even R (R Core Team, 2014) is useful to automate some parts of the workflow such as the `.sav` generation. Once written, scripts can be run several times or adapted, reducing the time required to perform these repetitive tasks. Furthermore 'Edgewarp' does not perform the Procrustes superimposition needed for the computation of the mean shape, entailing the use of another software package (here R with the 'Rmorph' library). Therefore the time needed cannot be spent doing other things. Conversely using 'Morpho,' the principal active component is the writing of the script. Computation can then be run as a background task or on a remote computer. This is an advantage as much for an initial analysis of the data as for modifications (for example adding specimens or correcting a landmark placement).
## Problems encountered

Deviations from the initial curves were observed in ‘Morpho.’ This deviation from the initial curve is a logical consequence of the relaxation procedure along tangents. Surface sliding semi-landmarks tend to move inside or outside the volume of the specimen depending on the local curvature as they are sliding on the tangent (Gunz, Mitteroecker & Bookstein, 2005). Curve sliding landmarks move outside of the curve on the opposite side of the center of curvature. This is why points are projected back onto the surface after relaxation. But in some cases such as for the greater trochanteric crest (2 in Fig. 3) of *Gulo gulo* and *Martes martes* this leads to deviations of the curve from its initial position. This effect is counter balanced in ‘Edgewarp’ by re-projecting points onto the ‘rail.’ This ‘rail’ (named ‘cur’ in ‘Edgewarp’) is constituted of a collection of points digitized on the bone surface, linked by segments. Consequently the definition of the ‘rail’ needs a large number of points to achieve a good description. However, when two curves are close, points can ‘jump’ from one curve to the other if the relaxed position is close to the second curve. Two strategies can counterbalance this effect in ‘Morpho’; first using a higher density of curve sliding landmarks (Gunz, Mitteroecker & Bookstein, 2005) and second, by reducing the relaxation step size. This last solution presents the advantage of being applicable on the initial dataset and does not require new data collection.

We also noticed points being projected on the wrong side of the supracondylar lateral crest and points projecting through bone from the olecranon fossa to the coronoid fossa while using ‘Edgewarp.’ The projection of surface sliding semi-landmarks from the template to the specimen surface is a crucial step preceding the actual sliding process. The main principle used in both software packages is a TPS deformation of the template constrained by anatomical and curve sliding landmarks followed by a projection onto the specimen surface. For some areas such as the lateral supracondylar crest, surface points can be projected onto the wrong side of the crest, if projected to the closest surface after TPS deformation. The ‘Morpho’ function ‘placePatch’ implements additional steps to avoid such effects. Inflation/deflation along the normals (vector giving the orientation of the mesh vertex) is performed; then the normals of projected points are compared to those of the template to avoid inside-out effects. In ‘Edgewarp,’ the graphical interface allows users to manually displace and correct wrongly projected points. The ‘Morpho’ strategy presents the advantage to be automated and integrated as a part of the ‘placePatch’ function. In contrast, ‘Edgewarp’ requires the manual displacement of the points, implying the correct identification of wrongly projected points prior to moving them to the correct side. This requires the investigation of all points on every specimen which is time consuming. The use of surface sliding semi-landmarks at the very vicinity of the curves appears to be a factor leading to projection on the wrong side of crests. Therefore, this problem could be circumvented by taking this fact into account during the design of the template.

Aside of the own properties of these software, the help file determines a part of the global ergonomy of any software package. Edgewarp's user manual (User's manual, EWSH3.19; available at ftp://brainmap.stat.washington.edu/edgewarp/) was last updated on March 2002. Additionally, it remains apparently unfinished as indicated by the author's comments
appearing along the text. On the other hand, Morpho's help file is frequently updated and
examples are provided for each function.

CONCLUSIONS

This study shows that the results obtained by two software packages for 3D shape
quantification, 'Edgewarp' (Bookstein & Green, 1994) and 'Morpho' (Schlager, 2013b),
are similar and that biological variability is similarly described by both software packages.
We highlight the fact that small structures jutting out of the surface as well as depressions
may have an impact on the results and underline the importance of identifying these
structures prior to the sliding operation. An increase in the density of the surface sliding
landmarks in these areas or the removal of structures that do not contain biologically
relevant information for the study are possible solutions for avoiding strong impact
of these structures on the morphospace. The present study focuses on a small number
of species, showing a large range of size and shape variation. This experimental design
highlights the kind of structures that drive divergence between the two software packages.
However, this design does not allow us to anticipate possible biases for different datasets
such as those used in intra-specific studies where the variability between specimens is often
quite subtle. Beyond the similarity in results, 'Morpho' is faster than 'Edgewarp' for the
same sliding task, notably as a result of the computational power of 'Morpho' and the
complexity of the workflow in 'Edgewarp.' 'Edgewarp,' on the other hand provides a visual
approach to the sliding procedure. Thanks to the exploitation of modern computational
power, 'Morpho' library, provides a powerful tool for the treatment of larger datasets.

ACKNOWLEDGEMENTS

The authors would like to thank S. Ivan Perez and an anonymous reviewer for their
comments on the original manuscript and their suggestions, S Schlager (University of
Freiburg, Germany) for advice and his reactivity, J Cuisin, G. Véron, J. Villemain, and
C Bens for the access to specimens from the collections Mammifères et Oiseaux, MNHN,
Paris and the UMS 2700 OMSI MNHN/CNRS for the access to the scanner. LB-D thanks
the doctoral school 'Frontières du vivant' and the Bettencourt Schueller foundation. LB-D
thanks M Depraetere (MNHN, Paris) for her useful comments on the manuscript.

ADDITIONAL INFORMATION AND DECLARATIONS

Funding

LB-D, AIH and AnH received financial support from the ANR-13-PDOC-001. A-CF
is supported by the Fondation Fyssen. The funders had no role in study design, data
collection and analysis, decision to publish, or preparation of the manuscript.

Grant Disclosures

The following grant information was disclosed by the authors:
Agence Nationale de la Recherche.
Fondation Fyssen.
**Competing Interests**

The authors declare there are no competing interests.

**Author Contributions**

* Léo Botton-Divet conceived and designed the experiments, performed the experiments, analyzed the data, contributed reagents/materials/analysis tools, wrote the paper, prepared figures and/or tables, reviewed drafts of the paper.

* Alexandra Houssaye and Anthony Herrel conceived and designed the experiments, reviewed drafts of the paper.

* Anne-Claire Fabre contributed reagents/materials/analysis tools, reviewed drafts of the paper.

* Raphael Cornette conceived and designed the experiments, contributed reagents/materials/analysis tools, reviewed drafts of the paper.

**Supplemental Information**

Supplemental information for this article can be found online at http://dx.doi.org/10.7717/peerj.1417#supplemental-information.

**REFERENCES**

**Adams DC, Otárola-Castillo E. 2013.** geomorph: an R package for the collection and analysis of geometric morphometric shape data. *Methods in Ecology and Evolution* **4**:393-399 DOI 10.1111/2041-210X.12035.

**Adams DC, Rohlf FJ, Slice DE. 2004.** Geometric morphometrics: ten years of progress following the “revolution”. *Italian Journal of Zoology* **71**:5-16 DOI 10.1080/11250000409356545.

**Adler D, Murdoch D. 2012.** rgl: 3D visualization device system (OpenGL). R package version 0.92. Available at https://cran.r-project.org/web/packages/rgl/index.html.

**Álvarez A, Ercoli MD, Prevosti FJ. 2013.** Locomotion in some small to medium-sized mammals: a geometric morphometric analysis of the penultimate lumbar vertebra, pelvis and hindlimbs. *Zoology* **116**:356-371 DOI 10.1016/j.zool.2013.08.007.

**Andersson KI. 2004.** Elbow-joint morphology as a guide to forearm function and foraging behaviour in mammalian carnivores. *Zoological Journal of the Linnean Society* **142**:91-104 DOI 10.1111/j.1096-3642.2004.00129.x.

**Bass SL, Saxon L, Daly RM, Turner CH, Robling AG, Seeman E, Stuckey S. 2002.** The effect of mechanical loading on the size and shape of bone in pre-, peri-, and postpubertal girls: a study in tennis players. *Journal of Bone and Mineral Research* **17**:2274-2280 DOI 10.1359/jbmr.2002.17.12.2274.

**Baylac M. 2013.** Rmorph, a morphometrics library for R. Available from the author. baylac@mnhn.fr.

**Beaupre GS, Orr TE, Carter DR. 1990.** An approach for time-dependent bone modeling and remodeling—application: a preliminary remodeling simulation. *Journal of Orthopaedic Research* **8**:662-670 DOI 10.1002/jor.1100080507.

**Blender Online Community. 2014.** Blender—a 3D modelling and rendering package. Blender Institute, Amsterdam: Blender Foundation.

**Bodkin JL. 2001.** Sea otters. In: *Encyclopedia of ocean sciences*. Oxford: Academic Press, 2614-2621.
**Bookstein FL. 1991.** *Morphometric tools for landmark data: geometry and biology*. Cambridge: Cambridge University Press.

**Bookstein FL. 1997.** Landmark methods for forms without landmarks: morphometrics of group differences in outline shape. *Medical Image Analysis* 1:225–243
DOI 10.1016/S1361-8415(97)85012-8.

**Bookstein FL, Green WD. 1994.** Edgewarp: a flexible program package for biometric image warping in two dimensions. In: *Visualization in biomedical computing*, vol. 2359. International society for optics and photonics, 135–147.

**Cornette R, Baylac M, Souter T, Herrel A. 2013.** Does shape co-variation between the skull and the mandible have functional consequences? A 3D approach for a 3D problem. *Journal of Anatomy* 223:329–336 DOI 10.1111/joa.12086.

**Cornette R, Tresset A, Herrel A. 2014.** The shrew tamed by Wolff's law: do functional constraints shape the skull through muscle and bone covariation?: the Shrew Tamed by Wolff'S Law. *Journal of Morphology* 276:301–309 DOI 10.1002/jmor.20339.

**Cubo J. 2004.** Pattern and process in constructional morphology. *Evolution & Development* 6:131–133 DOI 10.1111/j.1525-142X.2004.04018.x.

**Cubo J, Legendre P, De Ricqlès A, Montes L, De Margerie E, Castanet J, Desdevises Y. 2008.** Phylogenetic, functional, and structural components of variation in bone growth rate of amniotes. *Evolution & Development* 10:217–227 DOI 10.1111/j.1525-142X.2008.00229.x.

**Currey J. 2003.** The many adaptations of bone. *Journal of Biomechanics* 36:1487–1495
DOI 10.1016/S0021-9290(03)00124-6.

**De Groote I, Lockwood CA, Aiello LC. 2010.** Technical note: a new method for measuring long bone curvature using 3D landmarks and semi-landmarks. *American Journal of Physical Anthropology* 141:658–664 DOI 10.1002/ajpa.21225.

**Dryden IL, Mardia KV. 1993.** Multivariate shape analysis. *Sankhyā: The Indian Journal of Statistics, Series A* 55:460–480.

**Dryden IL, Mardia KV. 1998.** *Statistical analysis of shape*. Chichester: Wiley.

**Estes JA. 1980.** Enhydra lutris. *Mammalian Species* 133:1–8 DOI 10.2307/3503844.

**Fabre A-C, Cornette R, Peigné S, Goswami A. 2013a.** Influence of body mass on the shape of forelimb in musteloid carnivorans. *Biological Journal of the Linnean Society* 110:91–103
DOI 10.1111/bij.12103.

**Fabre A-C, Cornette R, Perrard A, Boyer DM, Prasad GV, Hooker JJ, Goswami A. 2014a.** A three-dimensional morphometric analysis of the locomotory ecology of Deccanolestes, a eutherian mammal from the Late Cretaceous of India. *Journal of Vertebrate Paleontology* 34:146–156 DOI 10.1080/02724634.2013.789437.

**Fabre A-C, Cornette R, Slater G, Argot C, Peigné S, Goswami A, Pouydebat E. 2013b.** Getting a grip on the evolution of grasping in musteloid carnivorans: a three-dimensional analysis of forelimb shape. *Journal of Evolutionary Biology* 26:1521–1535 DOI 10.1111/jeb.12161.

**Fabre A-C, Goswami A, Peigné S, Cornette R. 2014b.** Morphological integration in the forelimb of musteloid carnivorans. *Journal of Anatomy* 225:19–30 DOI 10.1111/joa.12194.

**Free Software Foundation. 2013.** Bash (4.3) Unix shell program. Available at http://ftp.gnu.org/gnu/bash/.

**Gower JC. 1975.** Generalized procrustes analysis. *Psychometrika* 40:33–51
DOI 10.1007/BF02291478.

**Gunz P, Mitteroecker P. 2013.** Semilandmarks: a method for quantifying curves and surfaces. *Hystrix*, the *Italian Journal of Mammalogy* 24:103–109.
**Gunz P, Mitteroecker P, Bookstein FL. 2005.** Semilandmarks in three dimensions. In: *Modern morphometrics in physical anthropology*. New York: Kluwer Academic/Plenum Publishers, 73–98.

**Heinrich RE, Biknevicius AR. 1998.** Skeletal allometry and interlimb scaling patterns in mustelid carnivorans. *Journal of Morphology* **235**:121–134 DOI 10.1002/(SICI)1097-4687(199802)235:2<121::AID-JMOR3>3.0.CO;2-C.

**Hunter L, Barrett P. 2011.** *Carnivores of the world*. Princeton: Princeton University Press.

**Kulemeyer C, Asbahr K, Gunz P, Frahnert S, Bairlein F. 2009.** Functional morphology and integration of corvid skulls—a 3D geometric morphometric approach. *Frontiers in Zoology* **6**:2 DOI 10.1186/1742-9994-6-2.

**Lanyon LE, Goodship AE, Pye C, MacFie J. 1982.** Mechanically adaptive bone remodelling. *Journal of biomechanics* **15**:141–154 DOI 10.1016/0021-9290(82)90246-9.

**Marcus LF. 1990.** *Traditional morphometrics*. Ann Arbor: University of Michigan Museum of Zoology, 77–122.

**Martín-Serra A, Figueirido B, Palmqvist P. 2014.** A three-dimensional analysis of the morphological evolution and locomotor behaviour of the carnivoran hind limb. *BMC Evolutionary Biology* **14**:129 DOI 10.1186/1471-2148-14-129.

**Meloro C, Raia P, Piras P, Barbera C, O'Higgins P. 2008.** The shape of the mandibular corpus in large fissiped carnivores: allometry, function and phylogeny. *Zoological Journal of the Linnean Society* **154**:832–845 DOI 10.1111/j.1096-3642.2008.00429.x.

**Monteiro LR, Nogueira MR. 2010.** Adaptive radiations, ecological specialization, and the evolutionary integration of complex morphological structures. *Evolution* **64**:724–744 DOI 10.1111/j.1558-5646.2009.00857.x.

**Morgan CC. 2009.** Geometric morphometrics of the scapula of South American caviomorph rodents (Rodentia: Hystricognathi): form, function and phylogeny. *Mammalian Biology*—Zeitschrift für Säugetierkunde **74**:497–506 DOI 10.1016/j.mambio.2008.09.006.

**Morita Y, Amano H, Nakatsukasa M, Kondo O, Ogihara N. 2014.** A geometric morphometric study of neurocranial shape variations in the crania of modern Japanese. In: *Dynamics of learning in neanderthals and modern humans*, vol. 2. Tokyo: Springer, 153–160.

**Nowak RM. 2005.** *Walker's carnivores of the world*. Baltimore: JHU Press.

**Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O'Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. 2015.** *vegan: community ecology package*. R package version 2.3-0.

**Parr WCH, Wroe S, Chamoli U, Richards HS, McCurry MR, Clausen PD, McHenry C. 2012.** Toward integration of geometric morphometrics and computational biomechanics: new methods for 3D virtual reconstruction and quantitative analysis of Finite Element Models. *Journal of Theoretical Biology* **301**:1–14 DOI 10.1016/j.jtbi.2012.01.030.

**Peres-Neto P, Jackson D. 2001.** How well do multivariate data sets match? The advantages of a Procrustean superimposition approach over the Mantel test. *Oecologia* **129**:169–178 DOI 10.1007/s004420100720.

**Perez SI, Bernal V, Gonzalez PN. 2006.** Differences between sliding semi-landmark methods in geometric morphometrics, with an application to human craniofacial and dental variation. *Journal of Anatomy* **208**:769–784 DOI 10.1111/j.1469-7580.2006.00576.x.

**Pizzo A, Mazzone F, Rolando A, Palestrini C. 2011.** Combination of geometric morphometric and genetic approaches applied to a debated taxonomical issue: the status of Onthophagus
massai (Coleoptera, Scarabaeidae) as an endemic species vicarious to Onthophagus fracticornis in Sicily. *Zoology* **114**:199–212 DOI **10.1016/j.zool.2011.03.003**.

**R Core Team. 2014.** *R: a language and environment for statistical computing*. Vienna: R Foundation for Statistical Computing.

**Rohlf FJ, Slice D. 1990.** Extensions of the Procrustes method for the optimal superimposition of landmarks. *Systematic Biology* **39**:40–59.

**Schlager S. 2013a.** Soft-tissue reconstruction of the human nose: population differences and sexual dimorphism. PhD Thesis, Universitätsbibliothek Freiburg.

**Schlager S. 2013b.** *Morpho: calculations and visualizations related to geometric morphometrics*. R package version 2.1-1.141011.

**Schutz H, Guralnick RP. 2007.** Postcranial element shape and function: assessing locomotor mode in extant and extinct mustelid carnivorans. *Zoological Journal of the Linnean Society* **150**:895–914 DOI **10.1111/j.1096-3642.2007.00303.x**.

**Skinner MM, Gunz P. 2010.** The presence of accessory cusps in chimpanzee lower molars is consistent with a patterning cascade model of development: molar cusp patterning in chimpanzees. *Journal of Anatomy* **217**:245–253 DOI **10.1111/j.1469-7580.2010.01265.x**.

**Souter T, Cornette R, Pedraza J, Hutchinson J, Baylac M. 2010.** Two applications of 3D semi-landmark morphometrics implying different template designs: the theropod pelvis and the shrew skull. *Comptes Rendus Palevol* **9**:411–422 DOI **10.1016/j.crpv.2010.09.002**.

**Van Rossum G, De Boer J. 1991.** Interactively testing remote servers using the Python programming language. *CWi Quarterly* **4**:283–303.

**Wiley DF, Amenta N, Alcantara DA, Ghosh D, Kil YJ, Delson E, Harcourt-Smith W, Rohlf FJ, St John K, Hamann B. 2005.** Evolutionary morphing. In: *Proceedings of IEEE visualization 2005*. Piscataway: IEEE, 431–438.

**Wilson DE, Mittermeier RA, Cavallini P, Llobet T. 2009.** *Handbook of the mammals of the world*. Barcelona: Lynx.

**Zelditch ML, Swiderski DL, Sheets HD. 2012.** *Geometric morphometrics for biologists: a primer*. San Diego: Elsevier.