Skip to main content

Impact of the terrestrial-aquatic transition on disparity and rates of evolution in the carnivoran skull



Which factors influence the distribution patterns of morphological diversity among clades? The adaptive radiation model predicts that a clade entering new ecological niche will experience high rates of evolution early in its history, followed by a gradual slowing. Here we measure disparity and rates of evolution in Carnivora, specifically focusing on the terrestrial-aquatic transition in Pinnipedia. We analyze fissiped (mostly terrestrial, arboreal, and semi-arboreal, but also including the semi-aquatic otter) and pinniped (secondarily aquatic) carnivorans as a case study of an extreme ecological transition. We used 3D geometric morphometrics to quantify cranial shape in 151 carnivoran specimens (64 fissiped, 87 pinniped) and five exceptionally-preserved fossil pinnipeds, including the stem-pinniped Enaliarctos emlongi. Range-based and variance-based disparity measures were compared between pinnipeds and fissipeds. To distinguish between evolutionary modes, a Brownian motion model was compared to selective regime shifts associated with the terrestrial-aquatic transition and at the base of Pinnipedia. Further, evolutionary patterns were estimated on individual branches using both Ornstein-Uhlenbeck and Independent Evolution models, to examine the origin of pinniped diversity.


Pinnipeds exhibit greater cranial disparity than fissipeds, even though they are less taxonomically diverse and, as a clade nested within fissipeds, phylogenetically younger. Despite this, there is no increase in the rate of morphological evolution at the base of Pinnipedia, as would be predicted by an adaptive radiation model, and a Brownian motion model of evolution is supported. Instead basal pinnipeds populated new areas of morphospace via low to moderate rates of evolution in new directions, followed by later bursts within the crown-group, potentially associated with ecological diversification within the marine realm.


The transition to an aquatic habitat in carnivorans resulted in a shift in cranial morphology without an increase in rate in the stem lineage, contra to the adaptive radiation model. Instead these data suggest a release from evolutionary constraint model, followed by aquatic diversifications within crown families.


Understanding factors which influence tempo and mode in evolution is an important theme in evolutionary biology [1-3]. One factor which may influence the evolutionary patterns is ecology [4,5]. In particular, the ‘adaptive radiation’ model of Simpson (1944) suggests that when organisms enter a new adaptive zone, that is a niche with relatively few competitors, there will be an initial burst in evolutionary rate [6,7]. This ‘early burst’ may be followed by a slowing of morphological diversification as the niche begins to become filled [2,8,9]. This model provides a potential link between ecological transitions and evolutionary rate.

Support for adaptive radiations has been found using empirical data [10-17], the best-known examples including Darwin’s finches [18-20], Hawaiian silverswords [21] and African lake cichlids [22,23]. Further, in terms of marine mammals, cetaceans underwent a rapid increase in body size disparity early in the clade’s history despite the lack of a rapid initial taxonomic diversification [24].

The mammalian order Carnivora is an ecologically and taxonomically diverse group of mammals which have been the source for many studies of morphological variation, though they display relatively low cranial disparity relative to other mammal orders [25-33]. Arguably the largest ecological transition in carnivoran evolution was the shift from terrestrial to aquatic lifestyle in the evolution of the Pinnipedia (seals, sea lions and walruses) [34]. This study investigates how this extreme ecological shift has influenced disparity (morphological diversity) and rates of evolution in pinniped skulls, in comparison to their fissiped relatives. The Carnivora provides an ideal case study for the influence of the terrestrial-aquatic transition on disparity and rates of evolution because both aquatic pinniped carnivorans (seals, sea lions, and walruses) and terrestrial (fissiped) carnivorans (dogs, bears, weasels, cats, hyaenas, mongooses and allies) have a large extant taxonomic diversity [33,35-37]. Moreover, pinnipeds are less divergent from their closest living fissiped relatives in cranial morphology than other marine mammal groups (e.g., cetaceans or sirenians) [38] and thus may be directly compared with fissipeds. Though increased rates of body size evolution in pinnipeds relative to fissipeds were not supported [39], rates of evolution in the cranium have not been compared.

In order to test hypotheses of cranial evolution in the Carnivora, we employed analyses that work within a phylogenetic framework. The phylogenetic relationships of fissiped and pinniped carnivorans used here are shown in Figure 1 [40,41]. There are 21 genera and 34–36 species divided into the three pinniped families: Phocidae (seals), Otariidae (fur seals and sea lions) and Odobenidae (walruses), which diverged ~29ma [40,42-44]. Molecular data generally indicate that odobenids are most closely related to otariids, forming the Otarioidea (Figure 1) whereas morphological data links the odobenids with phocids in the clade Phocomorpha (Figure 1) [37,44,45], so both hypotheses were employed here. The most basal pinnipeds were the enaliarctines, a stem radiation first known from California around 28ma [45], of which one representative fossil with excellent cranial preservation is included in this study (see Methods).

Figure 1
figure 1

Composite phylogeny used in this study. Extant relationships and branch lengths from [40], placement of fossils according to [41]. This shows the Otarioidea topology, with Odobenidae as sister taxon to the Otariidae. Analyses were also run on the same tree but with a Phocidae-Odobenidae sister grouping, following the Phocomorpha hypothesis. Branch colors: Feliformia, orange; non-pinniped Caniformia, red; stem pinnipeds and allodesmines, dark blue; Phocidae, mid-blue; Odobenidae, teal; Otariidae, pale blue.

The paraphyletic fissipeds are a group of mainly terrestrial (including arboreal and fossorial) carnivorans which consist of ten families, 105 genera and over 241 species (Figure 1); [46]. Pinnipeds are caniform carnivorans, and molecular and morphological evidence support the placement of pinnipeds within the arctoids (bears, raccoons, weasels, and allies). Within Arctoidea, there has been long-standing disagreement over whether musteloids or ursids are the sister group of pinnipeds (Figure 1) [35,36,45]. Most [35,40], but not all [37,47], recent evidence supports a closer relationship between musteloids and pinnipeds. Pinnipeds were also previously thought to be diphyletic, with otariids linked to ursids and phocids linked to mustelids [48,49]. However, a monophyletic origin of pinnipeds is now well supported by both morphological and molecular evidence [36,50].

To investigate the influence of the terrestrial-aquatic transition on morphological diversity, we tested evolutionary models on the carnivoran phylogeny. By combining a phylogeny with species’ information, statistical models can be used to infer the evolutionary past [51]. Such statistical models use specific parameters to determine how traits change in phylogenetic space. The most commonly used statistical model is Brownian Motion (‘BM’), which assumes that traits evolve in each instant of unit of time with a mean change of zero and unknown and constant variance. Within BM, the evolution of a continuous trait X, along a branch over time increment t, is quantified as dX(t) = σdB(t), where σ constitutes the magnitude of undirected, stochastic evolution (σ 2 is generally presented as the BM rate parameter) and dB(t) is Gaussian white noise.

Although generally agreed to be an unrealistic assumption for most analyses, the advantage of BM is that it is mathematically tractable. Recent phylogenetic comparative methods have continued using BM as a baseline model, but incorporate additional parameters to reflect more nuanced assumptions about the evolutionary process. Recent advances include the development of methods based on Ornstein-Uhlenbeck (‘OU’) assumptions. The OU model incorporates stabilizing selection as a constraint and quantifies the evolution of a continuous trait X as dX(t) = α[θX(t)]dt + σdB(t), where σ captures the stochastic evolution of BM, and α determines the rate of adaptive evolution towards an optimum trait value θ [52,53]. This standard OU model can be modified into multiple-optima OU models that allow optima to vary across the phylogeny [54]. In these implementations, the optima are defined a priori to allow testing of alternative parameterizations and therefore alternative biological hypotheses [54]. An added advantage of these OU model fitting approaches is that they allow proper multivariate models rather than fitting variables one at a time [55]. Other recent methods do not fix the number of shifts or their locations on the phylogeny, but instead implement algorithms that estimate them (e.g. reversible-jump), while jointly sampling OU parameters [56]. Such an approach was recently implanted in the R package bayou [57], and allows an inference of the location, magnitude, and number of adaptive shifts for univariate models. We hypothesize that the evolution of an aquatic lifestyle in pinnipeds will cause a selective regime shift at the basal node of pinnipeds which could be detected by these OU models.

Although OU-based methods are a powerful tool for testing alternative biological interpretations and ‘painting’ regime shifts onto branches of a phylogeny [54,57], they do not provide estimations of ancestral states nor of variable rates for individual branches. To overcome this empirical hurdle while avoiding increased model parameterization, [58] and [59] developed an approach (‘Independent Evolution’) that relies on similar assumptions as a multiple regime OU model, but requires fewer parameters. This approach assumes that population phenotypes are affected by the wandering adaptive peaks of adaptive surfaces (aligning with an Adaptive Peak model of evolution), and is therefore similar to an OU model with shifting locations in assuming that different regimes may occur at different locations in phylogenetic space. However, the formalization of IE assumptions differs considerably to the OU approaches, as the IE method utilizes a geometric approach that consists of two main steps: (1) quantifying an expected value for each internal node based only on phylogenetic and trait data, assuming a pure gradual mode of evolution; (2) quantifying the deviation of internal nodes to the expected value based on gradual evolution using a triangulation between the expected value and the two descendant values of the internal node in question. The triangulation between the gradual mode expectation of the ancestor and the observed descendants results in a rescaling of branches such that the barycentre among them provides the best fit for the data [58]; a procedure akin to Farris optimization [60,61]. The rate parameter is the distance between the barycentre of the triangulation procedure and the descendant value of the branch. This method therefore provides variable rate estimates for individual branches, as well as ancestral values for individual nodes. Under an adaptive radiation model, we would predict an increase in evolutionary rates on the branches at the base of Pinnipedia, relative to those higher up the pinniped phylogeny.

Measuring rates of evolution captures an important aspect of shape or trait evolution, but high rates of evolution may not translate simply into high diversity, or vice versa. For example, if taxa are constrained developmentally or ecologically to a particular range of shapes, they may show high rates of evolution and high amounts of convergence in form, but low overall morphological disparity [17,62,63]. In this scenario, analyses could accurately recover high rates of evolution but this would not show that the taxa of interest are repeatedly exploring the same range of morphospace rather than expanding into new morphologies. Alternatively, a clade could achieve high disparity through slow evolution if each shift moved into new regions of morphospace. Combining analyses of rates and morphological disparity thus allows for a more complete analysis of tempo and mode in the evolution of diversity. Under a model of diffusive evolution, we would predict that fissipeds would have greater disparity than fissipeds because fissipeds are the more taxonomically diverse and ancient clade.


Data set

This study included 151 specimens of adult carnivorans, including 64 fissiped and 87 pinniped specimens. These specimens constituted 34 fissiped species and 28 extant pinniped species (Table 1). Fissipeds are more taxonomically diverse than pinnipeds (241 and 36 species respectively), so fissipeds were relatively less densely sampled [46]. Specifically, we measured 14% of fissiped species, and 30% of fissiped genera. In contrast, 77% of pinniped species and all genera were included [46]. We were testing if the greater taxonomic diversity of fissipeds was reflected in a larger morphological disparity, so we sampled as broadly as possible within the group. Fissiped species were selected to encompass their full range of phylogenetic diversity, including representatives from every extant family and major ecological group. Although fossil fissipeds were unavailable for this study, previous studies suggest that most fossil fissipeds (except sabre-toothed cats) fall within the range of morphospace of extant clades [33]. The difference in sampling between fissipeds and pinnipeds may influence the results of variance-based disparity analyses as follows. Variance-based disparity measures may be overestimated in fissipeds relative to pinnipeds, because more dissimilar taxa were sampled for fissipeds, in order to cover their extant phylogenetic and ecological breadth without fully sampling their species diversity, than for pinnipeds, for which most extant and some fossil species were sampled. However, as our null hypothesis is that fissiped disparity should exceed pinniped disparity (due to their greater alpha diversity), this will not lead to a false rejection of the null hypothesis (type 1 error), but will make it more difficult to reject (type 2 error). Further, rarefaction will be used to statistically account for the differences in sampling (see below for details).

Table 1 Species list, sampling and taxonomic assignment

Several fossil pinnipeds were also available for study. Stem pinniped fossils were particularly important for the comparative analyses (disparity analyses were extant-only) because they provide information about the ancestral morphology of the group. These taxa may improve estimations of rates on the pinniped stem branches which spanned the terrestrial-aquatic transition. Cranial material of stem musteloids (sister taxa to pinnipeds) would have also been useful in this capacity, but were unavailable.

Enaliarctos emlongi (USNM 250345) is a stem pinniped from the Miocene (~20 mya) of California and represents an early radiation of enaliarctine pinnipeds. Allodesmus sp. (USNM 335445) is also from the Miocene of California but this specimen has not yet been carefully assessed in terms of its phylogenetic and taxonomic position. It is currently referred to the genus Allodesmus, which is a member of Desmatophocidae, an extinct pinniped family that has been variably related to either the Otarioidea or the Phocoidea [45,64-66]. Pontolis magnus (USNM 335567) is an extremely large fossil odobenid from the lower Pliocene Empire Formation of California and is thought to be closely related to the extinct Dusignathinae [67,68]. The two fossil phocids, Acrophoca longirostris (USNM 421632) and Piscophoca pacifica (USNM 360406) are both from the Pisco Formation of the Pliocene of Peru and are thought to be relatives of the monachine seals [69,70].

Specimens were obtained from the University of Cambridge Zoology Museum, the Natural History Museum (London), the United States National Museum of Natural History (Washington, D.C.), the Field Museum of Natural History (Chicago), and the American Museum of Natural History (New York). Specific details on landmark collection, unification, and mirroring (to fill in missing data) are described in [71,72] and [31]. From the datasets detailed in those studies, a subset of 11 overlapping landmarks, observable and with clear homology in both fissipeds and pinnipeds, and identifiable in all the fossils, was selected. These cranial landmarks are shown in Figure 2 and described in Table 2, and the complete dataset is available at This number is reduced from the original datasets because of the need to capture equivalent landmarks on disparate skull shapes and incomplete material. However, the common landmarks still include information from most regions of the skull and provide much more information about skull proportions than would equivalent linear measures.

Figure 2
figure 2

11 landmarks used in data analysis, shown on a skull of Arctocephalus gazella . Landmarks 4–7 were taken bilaterally. Landmark descriptions can be found in Table 2. Wireframe used to present shape variation from PCA shown in red.

Table 2 Description of landmarks taken

Two specimens of Ailuropoda melanoleuca were digitized by both observers to test inter-observer error. The average Procrustes distance between repeats (0.055) was smaller than the than the distance between different specimens of A. melanoleuca (0.081), and much less than the average Procrustes distance between fissiped specimens and the mean (0.114).


Three-dimensional landmark data were analyzed using the software MorphoJ [73]. Specimens were brought into the same shape space (removing all non-shape elements) via translation, rotation and scaling by Procrustes superimposition. Further, the Procrustes co-ordinates were symmetrized and the symmetric component of variation used in subsequent analyses. Previous work indicates that shape variation of the carnivoran cranium is limited enough to apply the tangent-space approximation [74]. Species averages were then calculated from specimen data so that each tip on the phylogenetic tree is represented by a single shape point [3]. Though most species were represented by multiple individuals sampling ranged from 1 to 7 individuals. Quantifying intraspecific variation is not critical here because we are measuring large-scale cross-taxonomic morphological variation. Any fossil specimens with high asymmetry were removed prior to analysis, as unusually high asymmetry may reflect distortion that occurred during or after fossilization. Principal components analysis was used to reduce the dataset to orthogonal principal components (PCs) for subsequent analysis. These PC data were used in disparity and rate analyses, as opposed to Procrustes co-ordinate scores, to provide more convenient and interpretable graphing of the evolutionary morphospace through time, and to overcome computing limits on variable numbers in some analyses. However, both data sets provide the same shape data when all the PCs are included.

Comparisons of disparity

The variation in shape within a group can be expressed in many different ways, and the disparity measure chosen may impact interpretation of the data. Hence, multiple measures of disparity should be used to fully understand the patterns of variation observed [3,75]. Range-based methods identify the maximal morphospace occupation of each group, and therefore are more sensitive to outliers and sample size bias. However they do inform about the maximum differences in shape within a group. Here, the sum of ranges across all PCs was used as a measure of range. Variance-based disparity metrics also include information about the spread of species within the morphospace, but are much more robust to sampling bias [3]. For example, a tightly clustered group with a few divergent species might have the same range as a widely dispersed group, but would have a lower variance. The mean multivariate distance to the group centroid was used as the variance-based disparity metric. Estimates of the total variance of each group should be less sensitive to uneven sampling than the convex hull area [75,76].

These metrics were calculated using the MDA package in Matlab [76,77] for fissipeds and extant pinnipeds. Fossil pinnipeds were excluded from this analysis to provide a fairer comparison to fissipeds, for which fossil specimens were unavailable. Disparity analyses were carried out using all principal component scores in order to take into account all of the variation in the dataset. To compute standard deviations and confidence intervals on the disparity measures, a bootstrap procedure was used. Specimens were re-sampled randomly with replacement 1000 times, and the mean (disparity value) and deviation calculated. In addition, to take into account the uneven sampling between fissipeds and pinnipeds, rarefaction was used. The pinniped sample was rarefied to a sample-size of five using bootstrapping. This represents an equivalent species-sampling to that of fissipeds.

Evolutionary rates

The phylogenetic relationships and branch lengths used for estimating evolutionary rate were taken from [40]. Fossil species were then added based on positions indicated by [41] and [70], with branch lengths equivalent to their earliest appearance in the fossil record (Figure 1). Due to controversy over the position of the walrus, and the potentially large impact of this key node to our interpretations, rate analyses were run using both the Otarioidea and Phocomorpha hypotheses.

A summary of comparative analyses run in this study can be found in Table 3. To investigate the relative likelihood of alternative evolutionary scenarios of shifts in cranial morphology in the carnivoran phylogeny, we used the ouch R-package [78] to fit Ornstein-Uhlenbeck models to multivariate data. To investigate the likelihood of an evolutionary scenario in which pinnipeds indicate a higher rate of change at their basal branch, we tested four alternative models. The first model is a neutral Brownian motion model, which assumes that evolutionary change and selection follow a random walk, and therefore does not align with adaptive radiation assumptions. The second model considers a single optimum along the entire phylogeny, consistent with a shared selective regime across carnivorans. The third model considers two separate optima, for terrestrial and aquatic carnivorans. This hypothesizes regime shifts both at the base of Pinnipedia and on the branch leading to the otter. Finally, the fourth model considers an adaptive shift at the base of pinnipedia, in concert with the reconstructed position of the terrestrial-aquatic transition. The analysis was run on the top 9 PCs which together account for 95.2% of variation (multivariate hypothesis test). The analyses included the top 9 PCs rather than all PCs due to the extensive time required to run the models.

Table 3 Rate and comparative analyses summary

To describe morphological changes along individual branches of the phylogeny we employed the R packages evomap [79] and bayou [57]. The IE method (available in evomap) provides both ancestral states and variable rates, allowing for a detailed description of how morphospace changed through evolutionary time along individual branches of the phylogeny. These changes can be visualized into an evolutionary morphospace that captures changes through time by taking snapshots of morphospace along intervals of time. This approach is different from the more widely used ‘phylo-morphospace’ approach [17] in that the evo-morphospace approach consists of visualizing the evolutionary changes in the morphospace through time, rather than the projection of a phylogeny into a morphospace. The evolutionary morphospace approach, available in evomap [79], thus fully captures evolutionary trends by displaying how morphospace is inferred to have changed over time through phylogenetic space. Rates and reconstructed node values were calculated for all PCs individually using evomap, of which the top 4 are presented (univariate IE analysis). From the node values, the Euclidean distance between ancestor–descendant pairs across all PCs was calculated, providing an estimate of total evolutionary change along each branch of the tree (multivariate IE analysis).

The method implemented by bayou allows the inference of the location, magnitude and number of adaptive shifts within a multiple-optima OU framework, hereby providing a detailed inference of selective shifts in the phylogeny (Univariate OU analysis). We applied this analysis to the top 4 PC’s as an independent validity check of results obtained by the IE method.


Cranial shape of fissiped and pinniped carnivorans

Figure 3 depicts the first and second principal components which account for 60.12% of the total variance. Illustration of shape variation is provided by the wireframes in Figure 4 and the comparison plate of skull photographs in Figure 5. PC1 distinguished phocid and odobenid pinnipeds from otariids and fissipeds. Phocids, with high scores on PC1, displayed large nasal openings, dorso-ventrally tall and mediolaterally wide crania. This is in comparison to the relatively dorsoventrally flat and mediolaterally narrow terrestrial skull that has more anteriorly positioned nasal bones. The aquatic mustelid Enhydra lutris (sea otter, superior-most orange star) fell close to otariid space and in the most positive position on PC1 among fissipeds, although some ursids (orange trefoil) also displayed relatively high PC1 scores. The stem pinniped E. emlongi and the desmatophocid Allodesmus sp. were similar to the otariids on this axis, whereas the fossil phocids P. pacifica and A. longirostris fell within or close to extant phocids.

Figure 3
figure 3

Scatterplots showing variation on PC1 - PC4. These axes represent 39.76%, 20.36%, 14.92% and 6.14% of variance. Based on species means. Fossil pinnipeds are as follows: En, Enaliarctos emlongi; Al, Allodesmus sp.; Po, Pontolis magnus; Pi, Piscophoca pacifica; Ac, Acrophoca longirostris. Polygons connect fissipeds (red) and extant pinnipeds (blue) and reflect groupings used in the disparity analyses. Extremal shapes are shown in Figure 4.

Figure 4
figure 4

Wireframes showing shape variation on PC1-PC4 in lateral and dorsal views. Anterior is to the right of the image. Landmarks the wireframe was based on are shown in Figure 2.

Figure 5
figure 5

Morphological variation in carnivoran skulls. A-C: representatives of the three families of pinnipeds. D-F: Examples of shape variation within fissipeds. Note the enlarged nasal opening typical of positive PC1 scores, found in phocids, odobenids, and the fissiped otter. Dog and cat represent dolichocephalic and brachiocephalic extremes respectively, reflected by PC2 score.

PC2 reflected snout elongation across Carnivora. Short-faced (brachycephalic) species, such as Enhydra, displayed positive scores on PC2, whereas long-snouted (dolichocephalic) species, like most canids, had negative scores. The fossil taxa A. longirostris (the long-necked seal), and Pontolis magnus (a stem odobenid), as well as the extant leopard seal, Hydrurga leptonyx, had crania that were very dolichocephalic.

Figure 3 also illustrates PC3 and 4, which combined represent 21.06% of the variance in the analyses. PC3 distinguished otariids from other pinniped and fissiped groups. Otariids differed in the more anterior placed parietal-occipital suture relative to the bullae, longer ventral portion of the jugal, and more anteriorly place nasofrontal relative to the jugal-maxillary suture. Both E. emlongi and Allodesmus sp. were near to, but just outside of, otariid morphospace on this axis. PC4 represents around 6% of variation and did not distinguish at all fissipeds from pinnipeds, which overlap fully on this axis.

Disparity analyses

Results of disparity analyses of fissiped and pinniped carnivorans can be found in Table 4. For pinnipeds, results are shown both for the full number of specimens, and for a sample which has been rarefied down to the equivalent fissiped sampling-level (proportional to their taxonomic diversity) using bootstrapping (five specimens). The sum of ranges, a range-based disparity metric, is more sensitive to sample size. When all pinnipeds are included, pinniped disparity was larger, but the pinniped confidence interval overlapped with the fissiped value, indicating that there was not a significant difference in disparity between the two groups. However, when pinnipeds are bootstrapped the range-based disparity falls drastically, such that it is below the lower confidence limit of fissipeds. The mean distance to the centroid, a variance-based disparity metric, was much more robust to sampling changes. In this case, both original- and rarefied-pinniped disparity was above the 95% confidence interval for fissipeds, indicating that pinnipeds had higher variance, which was robust to sampling.

Table 4 Results of disparity analyses

Evolutionary rates

Evolutionary rates on the carnivoran tree were assessed using both Independent Evolution (evomap) and Ornstein-Uhlenbeck (bayou and ouch) models. Results for the IE model can be found in Figures 6, 7, and Additional files 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10, bayou results in Additional files 11, 12, 13, 14, 15, 16, 17 and 18, and ouch results in Table 5.

Figure 6
figure 6

Evolutionary morphospace showing the reconstructed evolution of Carnivora on PC1 and PC3 through time. Based on IE analysis. Phylogeny and colors as shown in Figure 1. Points at the zero time point represent both nodes and tip values.

Figure 7
figure 7

Multivariate IE analysis results showing shape change on each branch. Multivariate evolutionary distances are calculated from node estimates based on IE analysis. Thicker branches represent greater morphological change on that branch.

Table 5 Results of multivariate hypothesis test in ouch

The AIC values for the four hypotheses tested using the multivariate ouch model-fitting approach are shown in Table 5. The BM model is supported for both Phocomorpha and Otarioidea. This suggests that there is no single selective regime influencing all carnivorans, nor is there a drastic shift in selective regime at the base of pinnipeds, but instead support a randomly varying adaptive peak.

The IE and bayou univariate rate analyses based on the Otarioidea hypothesis paint a similar picture of trait evolution for the first 4 PCs (Additional files 1, 11, 12, 13 and 14). On PC1, the highest evolutionary rates within pinnipeds are found in crown phocids and on the O. rosmarus branch, and similarly bayou finds regime shifts at the base of Phocidae and on the O. rosmarus branch. On PC2, which represents snout length, the canids have high evolutionary rates and high likelihood of a regime shift, associated with their dolichocephalic morphology, as does the long-snouted fossil pinniped A. longirostris. PC3 is associated with high evolutionary rates on, and a regime shift at, the base of the Otariidae. These patterns support the ouch results as they show that selection is acting in different morphological directions for the pinniped families (phocids and odobenids versus otariids), and cannot be united by a single selection regime at the base of the group. The results based on Phocomorpha are extremely similar to those based on Otarioidea, with the exception that the bayou model places an additional morphological shift on the stem Phocomorpha branch. Results for the Phocomorpha analyses can be found in Additional files 2, 15, 16, 17 and 18. In short, neither the evomap method nor the bayou method provides support for a rate or regime shift at the basal branch of pinnipeds.

Univariate IE analyses of individual PC scores can be mapped to show the morphological changes through time along individual branches of the phylogeny, producing an evolutionary morphospace graph. Figure 6 shows the evolutionary morphospace of PC1 against PC3, the two axes linked with fissiped-pinniped variation, for the Otarioidea hypotheses. Similar plots for PC1 and PC2, can be found in Additional file 3. Videos showing the evolutionary morphospace through time are presented in Additional files 4 and 5. Figure 7 shows that the fissipeds (orange and red) continue to evolve in the lower left region of the morphospace (negative PC1 and PC3) throughout the 60 million years. However, the appearance of crown pinnipeds marks the invasion of new areas of the morphospace. Crown phocids and otariids evolve quickly into the positive regions of PC1 and PC3 respectively, with many points showing similar trajectories. Results for the Phocomorpha analysis are essentially the same, and can be found in Additional files 6, 7, 8 and 9.

The total evolutionary change along each branch of the phylogeny was estimated with a multivariate IE analysis using all individual PCs. Results of this multivariate analysis are shown in Figure 7 for the Otarioidea hypothesis. The amount of shape change on each branch of the tree is indicated by its line thickness. Results of the same analysis using the Phocomorpha hypothesis are similar, and so are included in Additional file 10. These results do not indicate a marked morphological shift at the basal branch of the pinnipeds, but match the results of the previous analyses in suggesting evolutionary shifts higher in the pinniped tree. Within pinnipeds, the highest evolutionary shifts are found on branches leading to the walrus (Odobenus rosmarus), the hooded seal (Cystophora cristata), and within crown monachine seals. Branches connecting crown group pinnipeds, which reflect the original radiation of pinnipeds, generally indicate lower amount of shape change (Figure 7, Additional file 10). In the fissipeds, the highest amount of shape change tends to be at terminal branches, including those leading to Enhydra lutra, Mephitis mephitis, Proteles cristata, Nasua nasua and Potos flavus, but are also quite high on the stem branch of the Canidae.


During the ~43ma of their evolution, crown-group carnivorans have diversified into many ecological niches and geographical areas [33,80-82]. However, invasion of the aquatic environment is arguably the largest ecological transition in carnivoran evolution [44]. This study has demonstrated that although pinniped carnivorans have distinctive cranial morphology to fissipeds, and relatively large disparity for their taxonomic age; this high disparity was not achieved by rapid evolution at the base of Pinnipedia, as would be predicted under an adaptive radiation model. Instead, pinnipeds evolved at low-moderate rates at the base of the group, followed by later spikes of increased rate on terminal branches.

Fossil pinnipeds

Principal component analysis on cranial shape demonstrated that pinniped and fissiped carnivorans occupy distinctive regions of morphospace. In particular, phocids and odobenids are primarily distinguished from fissipeds on PC1, representing a dorsoventrally tall, mediolaterally wide skull, with high and posteriorly-placed nasal bones increasing the size of the nasal opening associated with dorsally-placed nares. Similarly, otariids are distinguished from fissipeds on PC3 by their anteriorly placed nasofrontal suture and long jugals. The inclusion of fossil pinnipeds in this study provides promising preliminary results concerning the extinct diversity and evolutionary history in this group. The stem pinniped, E. emlongi, shared similar cranial morphology to otariids, confirming previous qualitative work [48,83,84]. This also demonstrates that, despite the fact that E. emlongi is phylogenetically intermediate between fissipeds and pinnipeds, in terms of its skull morphology it can be fully distinguished from extant fissipeds on PC3, supporting inferences from post-crania that enaliarctines were likely fully adapted for an aquatic lifestyle [83,85].

Allodesmus sp. is intermediate between phocids and otariids, reflecting its phylogenetic position as a member of the Desmatophocidae, an enigmatic clade that is likely the sister group to one of those two families [45,64]. In contrast, the stem odobenid P.magnus lies far from the extant walrus in morphospace, with a relatively dolichocephalic morphology compared to the extant walrus. This suggests that the vast taxonomic diversity of extinct, stem odobenids may have encompassed a much larger amount of morphological disparity than other pinniped or fissiped clades and did not progress linearly from an Enaliarctos-type ancestral morphology to that of modern walruses [41]. Inclusion of more fossil walruses into future analyses could clarify evolutionary patterns in this diverse group.

One of the most unusual fossil pinnipeds, A. longirostris, the long-necked seal, has a score on PC1 similar to that of other phocids; however, it is more dolichocephalic than any living phocid. Dolichocephaly is very common in fissipeds and is thought to be related to the variable expression of the Runx2 gene [74,86]. This illustrates that some aspects of phenotypic variation are important in both fissipeds and pinnipeds, and is reflected by strong overlap between the two groups on PC2, representing dolichocephalic-brachiocephalic morphologies. In fact, A. longirostris is more dolichocephalic than any extant fissiped, including canids, raising questions about the functional implications of this divergent skull morphology. This species is known from the Pliocene of Peru and although it’s feeding ecology is not well understood, it has features consistent with the generalized pierce feeding typical of the group, but with interdigitate tooth cusps that have been linked with filter feeding in extant groups [87].

Fissiped and pinniped cranial disparity

Disparity analyses revealed that fissipeds do not have greater morphological disparity than extant pinnipeds, despite their greater taxonomic diversity and earlier divergence age. In this study, pinnipeds were sampled much more densely than fissipeds, which can lead to discrepancies in the measurement of disparity. In particular, fissiped specimens were selected broadly from across the group, to encompass their full range of diversity. However, this sampling scheme is likely to results in an inflated variance-based disparity, because taxa selected were more taxonomically diverse. Therefore, it is even more surprising that pinnipeds had greater variance than fissipeds. To further account for these sampling differences, the pinniped sample was rarefied down to a similar sampling level. While between-group patterns remained the same for variance-based disparity, range-based disparity was much lower in the rarefied sample. This reflects that range-based disparity is very sensitive to sample size, whereas variance-based metrics are not.

These data suggest that despite their younger taxonomic age and constituting relatively fewer species, pinnipeds have more variable skull morphology than the sampled fissipeds. This supports some previous data which showed fissipeds have relatively low cranial disparity compared to other mammalian groups [25], but is surprising given the large diversity of dietary specialization in the group [81]. Pinnipeds, on the other hand, are generally opportunistic feeders, but include a number of cranial specializations for prey capture and sexual display or combat [31]. These range from tusks and suction feeding specializations in the walrus to the bizarre inflatable nasal balloon in the hooded seal, and demonstrate the remarkable plasticity of cranial form in this group.

Heterogeneous disparity patterns such as these have been found in other carnivoran groups [74]. Domesticated dogs (a single species) have cranial disparity equivalent to the rest of the order, likely related to strong selection pressures during artificial selection. Those data showed that changing selection regimes (for features favored by dog breeders) have resulted in evolution of novel skull shapes, outside those found in wild Carnivora [74]. Despite this, modularity and integration patterns, which reflect the mechanisms generating natural variation, are not altered. Pinnipeds have also attained increased cranial diversity and evolved novel morphology outside the range of the other carnivorans. Hence both artificial and natural selection environments can influence cranial disparity in the Carnivora.

Adaptive radiation models in Pinnipedia

We did not find support for an increase in the rate of cranial evolution in stem pinniped lineages, nor a systematic shift in selective regime associated with the terrestrial-aquatic transition. Additionally, there was only a small amount of multivariate shape change on pinniped basal branches. All these results point to more gradual evolution across the terrestrial-aquatic transition. In contrast, the highest rates of evolution and strongest selective regime shifts were found within the crown group. The adaptive radiation model predicts that evolutionary rates should be high at the base of the radiation, then gradually slow as the niche fills. Therefore, these results do not provide support for an adaptive radiation model in the evolution of pinnipeds based on the current sample. This result was recovered using both Phocomorpha and Otarioidea arrangements, demonstrating that it is robust to the phylogenetic hypothesis used.

This result is in agreement with some other large-scale comparative studies which have shown that adaptive radiations may be quite rare [7,88]. In particular, very few groups were found which matched an early burst, followed by stasis, pattern. Those authors found that even though patterns tended to match constrained evolution, groups evolved too slowly for stabilizing selection [7]. Instead they suggest either oscillating selection or developmental constraints in shaping long-term evolutionary patterns. Similarly, [39] found similar rates of body size evolution in fissipeds and pinnipeds, suggesting that a large body size range in pinnipeds was achieved without an increase in evolutionary rate. Interestingly, angular measurements of the carnassials in carnivore evolution do display adaptive radiation-type patterns of evolution, establishing strong between-clade differences early in Carnivora [11,89,90]. However, the strong functional links between carnassial morphology and prey acquisition in fissiped carnivores may make this example a stronger candidate for an adaptive radiation model.

Despite the fact that evolutionary rates were low at the base of the pinniped radiation, increases in evolutionary rate were noted at other places within Pinnipedia. In particular, there were increased evolutionary rates in phocid terminal branches and the walrus branch. This matched the regime shifts detected by the OU model in phocids and Odobenus on PC1, indicating repeated selection and convergence toward a positive PC1 optimum. This suggests multiple instances of rapid evolution toward a tall cranium, with a highly expanded nasal opening and recessed nasal bones. This supports previous work based on similar data, which hypothesized convergent evolution of recessed nasals associated with sexual displays and sediment feeding in phocids [31]. For example, M. leonina and C. cristata have distinctive cranial morphology related to sexual dimorphism specializations, and both have relatively high rates of evolution [31]. Similarly, O. rosmarus is highly specialized for suction feeding, and has high estimated evolutionary rates, though fossil walruses were poorly sampled here. This suggests specialization for marine ecologies may be driving increased morphological rates in phocids, and that they may show relatively high plasticity with regard this type of morphological variation.

The bayou analysis, based on the Ornstein-Uhlenbeck model, and the evomap analysis, based on the independent evolution model, both suggest that the regime shifts for pinniped families were associated with different selective optima (positive PC1 vs PC3). This indicates selective divergence of the crown groups, and does not support a single unifying selective regime for the group. Phocids, odobenids and otariids differ in their ecology in several ways. For example, phocids inhabit higher latitudes, spend more time in water, use different swimming mechanisms and have more diverse dietary and reproductive strategies than otariids [44,85,91]. Therefore, there may have been bursts of evolution within crown Pinnipedia associated with exploiting new ecological opportunities within the aquatic realm.

As well as extrinsic factors, such as environment, it is also important to consider intrinsic factors [88]. One explanation for the high pinniped disparity is that evolution of an aquatic lifestyle released certain constraints acting on the terrestrial carnivoran skull. For example, mastication produces high loadings in the maxillary dentition which are transmitted through the rest of the skull [92-96]. Pinnipeds, and some other aquatic groups (e.g., cetaceans), do not masticate their food [44,97]. In concert with the loss of mastication, pinnipeds also evolved a shortened tooth row, homodont dentition and lost the typical mammalian tooth structure [34,87,98-101]. Instead of masticating with posterior teeth, pinniped behaviors (such as pierce feeding or male-male combat) involve biting with the anterior dentition. Loss of posterior tooth biting may therefore have changed the types of loadings experienced on the cranium as a whole and may have allowed pinnipeds to evolve new and disparate skull forms not accessible to fissipeds. The relaxation of functional constraints is one potential hypothesis for pinniped skull diversity which warrants future exploration. This hypothesis could be further investigated using functional analysis (e.g., FEA) of biting in the pinniped and fissiped cranium.

This theory is in line with the idea of an ‘escape from the terrestrial adaptive zone’. The concept of a terrestrial ‘adaptive zone’, or constrained area of morphospace associated with function, was originally suggested by [32] for the pinniped ankle. In order to examine the idea of a terrestrial adaptive zone in the cranium, we must first establish convergence of fissipeds within a confined region of morphospace and second demonstrate that morphological variation in the cranium is functionally related to their environment [32].

On the first point, there is strong overlap of feliform and fissiped caniform morphology on the first four PCs, suggesting convergence between species in both clades. Within the Caniformia, there is some overlap between mustelids and ursids, though canids have distinctive dolichocephalic cranial shapes. On the second point, there is ample evidence that skull morphology is strongly influenced by functional pressures. The morphology of the fissiped cranium is influenced by diet [28,33,81,102-106], sensory evolution [107] and thermoregulatory demands [108]. The pinniped skull is adapted for marine functions, such as aquatic prey acquisition, and polygynous reproductive strategies in island colonies [31,87,109]. Relative to their terrestrial ancestors, pinnipeds have specialized feeding systems [87,100,110-112], thermoregulation [113,114], sub-aquatic vision and hearing [115-119], diving [120,121] and intra-specific displays [31,122,123]. Therefore, we suggest that, similarly to the ankle, pinniped skull morphology may have diverged out of a constrained region of morphospace occupied by fissipeds by diffusive evolution, followed by stronger selection in different directions within crown pinniped families.


The evolution of aquatic habits in pinniped carnivorans presents an ideal case study of the effect of ecological transitions on evolutionary processes because both terrestrial and aquatic groups are well represented in modern faunas. This study suggests that pinnipeds are more diverse in cranial morphology than are fissipeds, despite their more recent evolutionary origin. However, this increased disparity was not achieved through higher rates of evolution, or selective regime shifts, associated with the terrestrial-aquatic transition. Instead pinniped evolution is characterized initially by low to moderate evolutionary rates. Later in pinniped evolution, there were several bursts of more rapid evolution, associated with shifts in selection regime toward divergent optima, which may reflect ecological diversification within the aquatic realm. Thus we demonstrate that very high evolutionary rates in stem lineages are not necessary to produce morphological diversity. Instead ecological diversification within the aquatic realm, and perhaps the release of terrestrial constraints, seem to be important in driving pinniped cranial disparity.



Principal components analysis


Principal components


United states national museum


Independent evolution model


Ornstein-Uhlenbeck model


  1. 1.

    Gould SJ. Tempo and mode in the macroevolutionary reconstruction of darwinism. Proc Natl Acad Sci U S A. 1994;91(15):6764–71.

    PubMed Central  CAS  PubMed  Google Scholar 

  2. 2.

    Simpson GG. Tempo and Mode in Evolution. New York: Columbia University Press; 1944.

    Google Scholar 

  3. 3.

    Foote M. The evolution of morphological diversity. Annu Rev Ecol Syst. 1997;28:129–52.

    Google Scholar 

  4. 4.

    Schluter D. The Ecology of Adaptive Radiation. Oxford: Oxfrod University Press; 2000.

    Google Scholar 

  5. 5.

    Losos JB. Adaptive radiation, ecological opportunity, and evolutionary determinism. Am Nat. 2010;175(6):623–39.

    PubMed  Google Scholar 

  6. 6.

    Losos JB, Miles DB. Testing the hypothesis that a clade has adaptively radiated: Iguanid lizard clades as a case study. Am Nat. 2002;160(2):147–57.

    PubMed  Google Scholar 

  7. 7.

    Harmon LJ, Losos JB, Davies TJ, Gillespie RG, Gittleman JL, Jennings WB, et al. Early bursts of body size and shape evolution are rare in comparative data. Evolution. 2010;64(8):2385–96.

    PubMed  Google Scholar 

  8. 8.

    Schluter D. Introduction to the symposium: Species interactions and adaptive radiation. Am Nat. 2000;156:S1–3.

    Google Scholar 

  9. 9.

    Foote M. Morphological disparity in Ordovician-Devonian crinoids and the early saturation of morphological space. Paleobiology. 1994;20(3):320–44.

    Google Scholar 

  10. 10.

    Raia P, Carotenuto F, Passaro F, Piras P, Fulgione D, Werdelin L, et al. Rapid action in the Palaeogene, the relationship between phenotypic and taxonomic diversification in Coenozoic mammals. Proc R Soc Lond Ser B Biol Sci. 2013;280(1750):20122244.

    CAS  Google Scholar 

  11. 11.

    Meloro C, Raia P. Cats and dogs down the tree: the tempo and mode of evolution in the lower carnassial of fossil and living Carnivora. Evol Biol. 2010;37(4):177–86.

    Google Scholar 

  12. 12.

    Harmon LJ, Schulte JA, Larson A, Losos JB. Tempo and mode of evolutionary radiation in iguanian lizards. Science. 2003;301(5635):961–4.

    CAS  PubMed  Google Scholar 

  13. 13.

    Mahler DL, Revell LJ, Glor RE, Losos JB. Ecological opportunity and the rate of morphological evolution in the diversification of Greater Antillean anoles. Evolution. 2010;64(9):2731–45.

    PubMed  Google Scholar 

  14. 14.

    Gillespie R. Community assembly through adaptive radiation in Hawaiian spiders. Science. 2004;303(5656):356–9.

    CAS  PubMed  Google Scholar 

  15. 15.

    Givnish TJ. Adaptive radiation, dispersal, and diversification of the Hawaiian lobeliads. In: Kato M, editor. The Biology of Biodiversity. Tokyo: Springer-Verlag; 1999. p. 67–90.

    Google Scholar 

  16. 16.

    Burbrink FT, Pyron RA. How does ecological opportunity influence rates of speciation, extinction, and morphological diversification in New World ratsnakes (tribe Lampropeltini)? Evolution. 2010;64(4):934–43.

    PubMed  Google Scholar 

  17. 17.

    Sidlauskas B. Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach. Evolution. 2008;62(12):3135–56.

    PubMed  Google Scholar 

  18. 18.

    Grant BR, Grant PR. How and Why Species Multiply: The Radiation of Darwin's Finches. Princeton: Princeton University Press; 2008.

    Google Scholar 

  19. 19.

    Grant PR, Grant BR. Evolution of character displacement in Darwin's finches. Science. 2006;313(5784):224–6.

    CAS  PubMed  Google Scholar 

  20. 20.

    Schluter D. Character displacement and the adaptive divergence of finches on islands and continents. Am Nat. 1988;131(6):799–824.

    Google Scholar 

  21. 21.

    Baldwin BG, Sanderson MJ. Age and rate of diversification of the Hawaiian silversword alliance (Compositae). Proc Natl Acad Sci U S A. 1998;95(16):9402–6.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. 22.

    Seehausen O. African cichlid fish: a model system in adaptive radiation research. Proc R Soc Lond Ser B Biol Sci. 2006;273(1597):1987–98.

    Google Scholar 

  23. 23.

    Kocher TD. Adaptive evolution and explosive speciation: The cichlid fish model. Nat Rev Genet. 2004;5(4):288–98.

    CAS  PubMed  Google Scholar 

  24. 24.

    Slater GJ, Price SA, Santini F, Alfaro ME. Diversity versus disparity and the radiation of modern cetaceans. Proc R Soc Lond Ser B Biol Sci. 2010;277(1697):3097–104.

    Google Scholar 

  25. 25.

    Marcus LF, Hingst-Zaher E, Zaher H. Application of landmark morphometrics to skulls representing the orders of living mammals. Hystrix. 2000;11(1):27–47.

    Google Scholar 

  26. 26.

    Wesley-Hunt GD. The morphological diversification of carnivores in North America. Paleobiology. 2005;31(1):35–55.

    Google Scholar 

  27. 27.

    Valkenburgh BV, Ruff CB. Canine tooth strength and killing behaviour in large carnivores. J Zool. 1987;212(3):379–97.

    Google Scholar 

  28. 28.

    Radinsky LB. Evolution of skull shape in carnivores 1: Representative modern carnivores. Biol J Linn Soc. 1981;15(4):369–88.

    Google Scholar 

  29. 29.

    Meloro C, Raia P, Piras P, Barbera C, O'Higgins P. The shape of the mandibular corpus in large fissiped carnivores: allometry, function and phylogeny. Zool J Linn Soc. 2008;154:832.

    Google Scholar 

  30. 30.

    Goswami A, Polly PD. The influence of modularity on cranial morphological disparity in Carnivora and Primates (Mammalia). PLoS ONE. 2010;5(3):e9517.

    PubMed Central  PubMed  Google Scholar 

  31. 31.

    Jones KE, Goswami A. Morphometric analysis of cranial shape in pinnipeds (Mammalia, Carnivora): convergence, ecology, ontogeny, and dimorphism. In: Goswami A, Friscia A, editors. Carnivoran Evolution: new Views on Phylogeny, Form, and Function. Cambridge: Cambridge University Press; 2010. p. 342–73.

    Google Scholar 

  32. 32.

    Polly PD. Adaptive zones and the pinniped ankle: a three-dimensional quantitative analysis of carnivoran tarsal evolution. In: Sargis EJ, Dagosto M, editors. Mammalian Evolutionary Morphology. Netherlands: Springer; 2008. p. 167–96.

    Google Scholar 

  33. 33.

    Goswami A, Milne N, Wroe S. Biting through constraints: cranial morphology, disparity and convergence across living and fossil carnivorous mammals. Proc R Soc Lond Ser B Biol Sci. 2011;278(1713):1831–9.

    Google Scholar 

  34. 34.

    Repenning CA. Adaptive evolution of sea lions and walruses. Syst Biol. 1976;25:375–90.

    Google Scholar 

  35. 35.

    Eizirik E, Murphy WJ, Koepfli KP, Johnson WE, Dragoo JW, Wayne RK, et al. Pattern and timing of diversification of the mammalian order Carnivora inferred from multiple nuclear gene sequences. Mol Phylogen Evol. 2010;56(1):49–63.

    CAS  Google Scholar 

  36. 36.

    Flynn JJ, Finarelli JA, Zehr S, Hsu J, Nedbal MA. Molecular phylogeny of the Carnivora (Mammalia): Assessing the impact of increased sampling on resolving enigmatic relationships. Syst Biol. 2005;54(2):317–37.

    PubMed  Google Scholar 

  37. 37.

    Agnarsson I, Kuntner M, May-Collado LJ. Dogs, cats, and kin: A molecular species-level phylogeny of Carnivora. Mol Phylogen Evol. 2010;54(3):726–45.

    CAS  Google Scholar 

  38. 38.

    Monteiro ELD, Monteiro LR, dos Reis SF. Skull shape and size divergence in dolphins of the genus Sotalia: A tridimensional morphometric analysis. J Mammal. 2002;83(1):125–34.

    Google Scholar 

  39. 39.

    Slater GJ, Harmon LJ, Wegmann D, Joyce P, Revell LJ, Alfaro ME. Fitting models of continuous trait evolution to incompletely sampled comparative data using approximate bayesian computation. Evolution. 2012;66(3):752–62.

    PubMed  Google Scholar 

  40. 40.

    Nyakatura K, Bininda-Emonds ORP. Updating the evolutionary history of Carnivora (Mammalia): a new species-level supertree complete with divergence time estimates. BMC Biol. 2012;10:12.

    PubMed Central  PubMed  Google Scholar 

  41. 41.

    Deméré TA, Berta A, Adam PJ. Pinnipedimorph evolutionary biogeography. Bull Am Mus Nat Hist. 2003;279:32.

    Google Scholar 

  42. 42.

    Koretsky IA, Sanders AE. Paleontology of the Late Oligocene Ashley and Chandler Bridge Formations of South Carolina, 1:Paleogene pinniped remains; the oldest known seal (Carnivora: Phocidae). In: Cenozoic Mammals of the Land and Sea: Tributes to the Career of Clayton E Ray. Smithsonian: Cont Paleo; 2002. p. 179–83.

    Google Scholar 

  43. 43.

    Arnason U, Gullberg A, Janke A, Kullberg M, Lehman N, Petrov EA, et al. Pinniped phylogeny and a new hypothesis for their origin and dispersal. Mol Phylogen Evol. 2006;41:345–54.

    Google Scholar 

  44. 44.

    Berta A, Sumich J, Kovaks K. Marine Mammals Evolutionary Biology. San Diego: Academic Press; 2006.

    Google Scholar 

  45. 45.

    Berta A, Wyss AR. Pinniped phylogeny. In: Contributions in marine mammal paleontology honoring Frank C Whitmore, Jr. Proceedings of the San Diego Society of Natural History. 1994;29:33–56.

    Google Scholar 

  46. 46.

    Wilson DE, Reeder DM. Mammal Species of the World: A Taxonomic and Geographic Reference. 3rd ed. Baltimore: Johns Hopkins University Press; 2005.

    Google Scholar 

  47. 47.

    Luan P, Ryder OA, Davis H, Zhang Y, Yu L. Incorporating indels as phylogenetic characters: Impact for interfamilial relationships within Arctoidea (Mammalia: Carnivora). Mol Phylogen Evol. 2012;66(3):748–56.

    Google Scholar 

  48. 48.

    Tedford RH. Relationship of pinnipeds to other carnivores (Mammalia). Syst Zool. 1976;25(4):363–74.

    Google Scholar 

  49. 49.

    Mclaren IA. Are the Pinnipedia biphyletic? Syst Zool. 1960;9(1):18–28.

    Google Scholar 

  50. 50.

    Wyss AR. Evidence from flipper structure for a single origin of pinnipeds. Nature. 1988;334:427.

    Google Scholar 

  51. 51.

    Pagel M. Inferring the historical patterns of biological evolution. Nature. 1999;401(6756):877–84.

    CAS  PubMed  Google Scholar 

  52. 52.

    Hansen TF. Stabilizing selection and the comparative analysis of adaptation. Evolution. 1997;51(5):1341–51.

    Google Scholar 

  53. 53.

    Hansen TF. Adaptive landscapes and macroevolutionary dynamics. In: Svensson E, Calsbeek R, editors. The Adaptive Landscape in Evolutionary Biology. Oxford, UK: Oxford University Press; 2012. p. 205–26.

    Google Scholar 

  54. 54.

    Butler MA, King AA. Phylogenetic comparative analysis: A modeling approach for adaptive evolution. Am Nat. 2004;164(6):683–95.

    Google Scholar 

  55. 55.

    Butler MA, King AA. Multivariate comparative analysis using OUCH. Integr Comp Biol. 2009;49:E24–4.

    Google Scholar 

  56. 56.

    Uyeda JF, Harmon LJ. A novel Bayesian method for inferring and interpreting the dynamics of adaptive landscapes from phylogenetic comparative data. Syst Biol. 2014;63(6):902–18.

    PubMed  Google Scholar 

  57. 57.

    Uyeda JF, Eastman J, Harmon LJ. bayou: Bayesian fitting of Ornstein-Uhlenbeck models to phylogenies. R package version 1.0.1. 2014.

    Google Scholar 

  58. 58.

    Smaers JB, Vinicius L. Inferring macro-evolutionary patterns using an adaptive peak model of evolution. Evol Ecol Res. 2009;11(7):991–1015.

    Google Scholar 

  59. 59.

    Smaers JB, Dechmann DKN, Goswami A, Soligo C, Safi K. Comparative analyses of evolutionary rates reveal different pathways to encephalization in bats, carnivorans, and primates. Proc Natl Acad Sci U S A. 2012;109(44):18006–11.

    PubMed Central  CAS  PubMed  Google Scholar 

  60. 60.

    Farris JS. Estimating phylogenetic trees from distance matrices. Am Nat. 1972;106(951):645.

    Google Scholar 

  61. 61.

    Farris JS, Farris JS. Methods for computing Wagner trees. Syst Zool. 1970;19(1):83.

    Google Scholar 

  62. 62.

    Goswami A, Smaers JB, Soligo C, Polly PD. The macroevolutionary consequences of phenotypic integration: from development to deep time. Philos Trans R Soc Lond, Ser B: Biol Sci. 2014;369(1649):20130254.

    CAS  Google Scholar 

  63. 63.

    Slater GJ. Phylogenetic evidence for a shift in the mode of mammalian body size evolution at the Cretaceous-Palaeogene boundary. Methods Ecol Evol. 2013;4(8):734–44.

    Google Scholar 

  64. 64.

    Demere TA, Berta A. The Miocene pinniped Desmatophoca oregonensis condon, 1906 (Mammalia: Carnivora), from the Astoria Formation, Oregon. Smithson Contrib Paleobiol. 2002;93:113–47.

    Google Scholar 

  65. 65.

    Repenning CA, Tedford RH. Otarioid seals of the Neogene. US Geol Survey Prof Pap. 1977;992:1–93.

    Google Scholar 

  66. 66.

    Barnes LG, Kiyoharu H. Miocene pinnipeds of the otariid subfamily Allodesminae in the north pacific ocean: systematics and relationships. Island Arc. 1994;3:329–60.

    Google Scholar 

  67. 67.

    Kohno N. A new Miocene odobenid (Mammalia : Carnivora) from Hokkaido, Japan, and its implications for odobenid phylogeny. J Vert Paleontol. 2006;26(2):411–21.

    Google Scholar 

  68. 68.

    Demere TA. The family Odobenidae:a phylogenetic analysis of fossil and living taxa. Proc San Diego Soc Nat Hist. 1994;29:99–123.

    Google Scholar 

  69. 69.

    Walsh S, Naish D. Fossil seals from Late Neogene deposits in South America: A new pinniped (Carnivora, Mammalia) assemblage from Chile. Palaeontology. 2002;45:821–42.

    Google Scholar 

  70. 70.

    Berta A, Kienle S, Bianucci G, Sorbi S: A re-evaluation of Pliophoca etrusca (Pinnipedia, Phocidae) from the Pliocene of Italy: phylogenetic and biogeographic implications. J Vert Paleontol 2014. in press.

  71. 71.

    Goswami A. Cranial modularity shifts during mammalian evolution. Am Nat. 2006;168(2):270–80.

    PubMed  Google Scholar 

  72. 72.

    Goswami A. Morphological integration in the carnivoran skull. Evolution. 2006;60(1):169–83.

    PubMed  Google Scholar 

  73. 73.

    Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Res. 2011;11(2):353–7.

    Google Scholar 

  74. 74.

    Drake AG, Klingenberg CP. Large-scale diversification of skull shape in domestic dogs: Disparity and modularity. Am Nat. 2010;175(3):289–301.

    PubMed  Google Scholar 

  75. 75.

    Ciampaglio CN, Kemp M, McShea DW. Detecting changes in morphospace occupation patterns in the fossil record: characterization and analysis of measures of disparity. Paleobiology. 2001;27(4):695–715.

    Google Scholar 

  76. 76.

    Navarro N. MDA: a MATLAB-based program for morphospace-disparity analysis. Comput Geosci. 2003;29(5):655–64.

    Google Scholar 

  77. 77.

    Navarro N. MDA: Morphospace-disparity analysis for Matlab. 12th ed. Dijon: Biogeosciences; 2001.

    Google Scholar 

  78. 78.

    Butler MA, King AA. ouch: Ornstein-Uhlenbeck models for phylogenetic comparative hypotheses. R package. 2009.

    Google Scholar 

  79. 79.

    Smaers JB: evomap: R package for the evolutionary mapping of continuous traits. In.:; 2014.

  80. 80.

    Van Valkenburgh B. Locomotor diversity within past and present guilds of large predatory mammals. Paleobiology. 1985;11(4):406–28.

    Google Scholar 

  81. 81.

    Van Valkenburgh B. Deja vu: the evolution of feeding morphologies in the Carnivora. Integr Comp Biol. 2007;47(1):147–63.

    PubMed  Google Scholar 

  82. 82.

    Christiansen P, Wroe S. Bite forces and evolutionary adaptations to feeding ecology in carnivores. Ecology. 2007;88(2):347–58.

    PubMed  Google Scholar 

  83. 83.

    Berta A. New Enaliarctos (Pinnipedimorpha) from the Oligocene and Miocene of Oregon and the Role of "Enaliarctids" in Pinniped Phylogeny. Washington D.C: Smithsonian Inst. Press; 1991.

    Google Scholar 

  84. 84.

    Hunt RM, Barnes LG. Basicranial evidence for ursid affinities of the oldest pinnipeds. In: Berta A, Demere T, editors. Marine Mammal Paleontology Honoring Frank C Whitmore, Jr. San Diego: Proceedings of the San Diego Museum of Natural History; 1994. p. 268.

    Google Scholar 

  85. 85.

    Bebej RM. Swimming mode inferred from skeletal proportions in the fossil pinnipeds Enaliarctos and Allodesmus (Mammalia, Carnivora). J Mamm Evol. 2009;16(2):77–97.

    Google Scholar 

  86. 86.

    Sears KE, Goswami A, Flynn JJ, Niswander LA. The correlated evolution of Runx2 tandem repeats, transcriptional activity, and facial length in Carnivora. Evol Dev. 2007;9(6):555–65.

    CAS  PubMed  Google Scholar 

  87. 87.

    Adam PJ, Berta A. Evolution of prey capture strategies and diet in the Pinnipedimorpha (Mammalia, Carnivora). Oryctos. 2002;4:83–107.

    Google Scholar 

  88. 88.

    Wagner CE, Harmon LJ, Seehausen O. Ecological opportunity and sexual selection together predict adaptive radiation. Nature. 2012;487(7407):366–U124.

    CAS  PubMed  Google Scholar 

  89. 89.

    Crusafont-Pairo M, Truyols-Santonja J. A biometric study of evolution of fissiped carnivores. Evolution. 1956;10:314–32.

    Google Scholar 

  90. 90.

    Werdelin L. Carnivoran ecomorphology: A phylogenetic perspective. In: Gittleman JL, editor. Carnivore Behavior, Ecology and Evolution, vol. 1. Ithica, NY: Cornell University Press; 1996. p. 582–624.

    Google Scholar 

  91. 91.

    Ferguson SH. How seals divide up the world: environment, life history, and conservation. Oecologia. 2006;150:318.

    PubMed  Google Scholar 

  92. 92.

    Hylander WL, Picq PG, Johnson KR. Masticatory-stress hypotheses and the supraorbital region of primates. Amer J Phys Anthrop. 1991;86(1):1–36.

    CAS  PubMed  Google Scholar 

  93. 93.

    Smith M, Savage RJ. The mechanics of mammalian jaws. School Sci Rev. 1959;141:289–301.

    Google Scholar 

  94. 94.

    Ross CF, Patel BA, Slice DE, Strait DS, Dechow PC, Richmond BG, et al. Modeling masticatory muscle force in finite element analysis: Sensitivity analysis using principal coordinates analysis. Anat Rec. 2005;283A(2):288–99.

    Google Scholar 

  95. 95.

    Vinyard CJ, Wall CE, Williams SH, Hylander WL. Comparative functional analysis of skull morphology of tree-gouging primates. Amer J Phys Anthrop. 2003;120(2):153–70.

    PubMed  Google Scholar 

  96. 96.

    Kupczik K, Dobson CA, Crompton RH, Phillips R, Oxnard CE, Fagan MJ, et al. Masticatory loading and bone adaptation in the supraorbital torus of developing macaques. Amer J Phys Anthrop. 2009;139(2):193–203.

    CAS  PubMed  Google Scholar 

  97. 97.

    Berta A, Adam PJ. Evolutionary Biology of Pinnipeds. In: de Buffrenil V, Mazin J-M, editors. Secondary adaptation of tetrapods to a life in water. Munchen Germany: Verlag Dr Frederich Pfeil; 2001. p. 235–60.

    Google Scholar 

  98. 98.

    Boessenecker RW. New records of the fur seal Callorhinus (Carnivora: Otariidae) from the Plio-Pleistocene Rio Dell Formation of Northern California and comments on otariid dental evolution. J Vert Paleontol. 2011;31(2):454–67.

    Google Scholar 

  99. 99.

    Jones KE, Ruff CB, Goswami A. Morphology and biomechanics of the pinniped jaw: mandibular evolution without mastication. Anat Rec. 2013;296(7):1049–63.

    Google Scholar 

  100. 100.

    Miller EH, Sung HC, Moulton VD, Miller GW, Finley JK, Stenson GB. Variation and integration of the simple mandibular postcanine dentition in two species of phocid seal. J Mammal. 2007;88(5):1325–34.

    Google Scholar 

  101. 101.

    Salazar-Ciudad I, Jernvall J. A computational model of teeth and the developmental origins of morphological variation. Nature. 2010;464(7288):583–U138.

    CAS  PubMed  Google Scholar 

  102. 102.

    Tseng ZJ, Wang XM. Cranial functional morphology of fossil dogs and adaptation for durophagy in borophagus and epicyon (carnivora, Mammalia). J Morphol. 2010;271(11):1386–98.

    PubMed  Google Scholar 

  103. 103.

    Sacco T, Van Valkenburgh B. Ecomorphological indicators of feeding behaviour in the bears (Carnivora: Ursidae). J Zool. 2004;263:41–54.

    Google Scholar 

  104. 104.

    Slater GJ, Dumont ER, Van Valkenburgh B. Implications of predatory specialization for cranial form and function in canids. J Zool. 2009;278(3):181–8.

    Google Scholar 

  105. 105.

    Wroe S. Bite club: comparative bite force in big biting mammals and the prediction of predatory behaviour in fossil taxa. Proc R Soc Lond Ser B Biol Sci. 2005;272:619.

    Google Scholar 

  106. 106.

    Tseng ZJ, Binder WJ. Mandibular biomechanics of Crocuta crocuta, Canis lupus, and the late Miocene Dinocrocuta gigantea (Carnivora, Mammalia). Zool J Linn Soc. 2009;158(2):638–96.

    Google Scholar 

  107. 107.

    Finarelli JA, Goswami A. The evolution of orbit orientation and encephalization in the Carnivora (Mammalia). J Anat. 2009;214(5):671–8.

    PubMed Central  PubMed  Google Scholar 

  108. 108.

    Van Valkenburgh B, Friscia A, Theodor J: Respiratory turbinals in felids and canids: a quantitative comparison. J Zool 2004, 264(3).

  109. 109.

    Jones KE, Goswami A. Quantitative analysis of the influences of phylogeny and ecology on phocid and otariid pinniped (Mammalia; Carnivora) cranial morphology. J Zool. 2010;280(3):297–308.

    Google Scholar 

  110. 110.

    Kastelein RA. Oral suction of a Pacific walrus (Odobenus rosmarus divergens) in air and under water. Mamm Biol. 1994;59:105.

    Google Scholar 

  111. 111.

    Marshall CD. Feeding kinematics, suction and hydraulic jetting capabilities in bearded seals (Erignathus barbatus). J Exp Biol. 2008;211:699.

    PubMed  Google Scholar 

  112. 112.

    King JE. The feeding mechanism and jaws of the crabeater seal (Lobodon carcinophagus). Annls Mus Hist nat Paris. 1961;25(4):462–6.

    Google Scholar 

  113. 113.

    Van Valkenburgh B, Curtis A, Samuels JX, Bird D, Fulkerson B, Meachen-Samuels J, et al. Aquatic adaptations in the nose of carnivorans: evidence from the turbinates. J Anat. 2011;218(3):298–310.

    PubMed Central  PubMed  Google Scholar 

  114. 114.

    Mellish JAE, Horning M, York AE. Seasonal and spatial blubber depth changes in captive harbor seals (Phoca vitulina) and Steller's sea lions (Eumetopias jubatus). J Mammal. 2007;88(2):408–14.

    Google Scholar 

  115. 115.

    Kastelein RA, Zweypfenning CVJ, Spekreijse H, Dubbeldam JL, Born EW. The Anatomy of the walrus head (Odobenus rosmarus). Part 3: The eyes and their function in Walrus ecology. Aquat Mamm. 1993;19(2):61–9.

    Google Scholar 

  116. 116.

    Hunt RM. Auditory bulla in Carnivora: Anatomical basis for reappraisal of carnivore evolution. J Morphol. 1974;143(1):21–75.

    PubMed  Google Scholar 

  117. 117.

    Nummela S, Thewissen JGM. Sensory Evolution on the Threshold: Adaptations in Secondarily Aquatic Vertebrates. Berkeley: University of California Press; 2008.

    Google Scholar 

  118. 118.

    Kastak D, Schusterman RJ. Low-frequency amphibious hearing in pinnipeds: Methods, measurements, noise, and ecology. J Acoust Soc Am. 1998;103(4):2216–28.

    CAS  PubMed  Google Scholar 

  119. 119.

    Kastak D, Schusterman RJ. In-air and underwater hearing sensitivity of a northern elephant seal (Mirounga angustirostris). Can J Zool. 1999;77(11):1751–8.

    Google Scholar 

  120. 120.

    Berg L, Pyenson N. Osteological correlates and phylogenetic analysis of deep diving in living and extinct pinnipeds: what good are big eyeballs? J Vert Paleontol. 2008;28:51A.

    Google Scholar 

  121. 121.

    Stenfors LE, Sade J, Hellstrom S, Anniko M. How can the hooded seal dive to a depth of 1000 m without rupturing its tympanic membrane? A morphological and functional study. Acta Oto-Laryngol. 2001;121(6):689–95.

    CAS  Google Scholar 

  122. 122.

    Sanvito S, Galimberti F, Miller EH. Having a big nose: Structure, ontogeny, and function of the elephant seal proboscis. Can J Zool. 2007;85(2):207–20.

    Google Scholar 

  123. 123.

    Bartholomew GA. A model for evolution of pinniped polygyny. Evolution. 1970;24(3):546.

    Google Scholar 

Download references


We would like to thank Dr. Leandro Monteiro and three anonymous reviewers for their helpful suggestions regarding both analyses and text. For access to collections, we would like to thank Nick Pyenson, Matt Lowe, Richard Sabin, Bill Stanley, Jean Spence, Paula Jenkins, and Linda Gordon. We are also grateful to Matthew Carrano for use of his digitizer whilst at the USNM, Dave Polly for helpful discussion, and Abbey Drake for advice on analyses. Funding during fissiped data collection was provided by the National Science Foundation DDIG #0308765, the University of Chicago Hinds Fund, and an AMNH collections study grant. We would like to thank Gavin Thomas and one anonymous reviewer for their comments on an earlier version of this manuscript.

Supporting data

3D landmark data used in these analyses are available from

Author information



Corresponding author

Correspondence to Katrina E Jones.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KEJ designed the study, collected pinniped dataset, performed shape and disparity analyses, and wrote the manuscript. JBS performed rate analyses and created figures. AG designed the study, and collected fissiped dataset. All authors contributed to manuscript preparation and approved the final manuscript.

Authors’ information

KEJ, Research Scholar, Museum of Comparative Zoology, Harvard University.

JS, Assistant Professor, Department of Anthropology, Stony Brook University.

AG, Reader in Paleobiology, Department of Genetics, Evolution & Environment and Department of Earth Sciences, University College London.

Additional files

Additional file 1:

UnivariateIE_Otarioidea.pdf, figure, PC1-4 from IE analysis. Thickness of branch is proportional to rate. Green, positive direction; red, negative direction.

Additional file 2:

UnivariateIE_Phocomorpha.pdf, figure, PC1-4 from IE analysis. Thickness of branch is proportional to rate. Green, positive direction; red, negative direction.

Additional file 3:

IE_Evomorphospace_PC1PC2_otarioidea.pdf, figure, Evolutionary morphospace for PC1 vs PC2 for otarioidea.

Additional file 4:, movie, Evolutionary morphospace movie for PC1 vs PC3 for otarioidea.

Additional file 5:, movie, Evolutionary morphospace movie for PC1 vs PC2 for otarioidea.

Additional file 6:

IE_Evomorphospace_PC1PC2_phocomorpha.pdf, figure, Evolutionary morphospace for PC1 vs PC2 for phocomorpha.

Additional file 7:

IE_Evomorphospace_PC1PC3_phocomorpha.pdf, figure, Evolutionary morphospace for PC1 vs PC3 for phocomorpha.

Additional file 8:, movie, Evolutionary morphospace movie for PC1 vs PC2 for phocomorpha.

Additional file 9:, movie, Evolutionary morphospace movie for PC1 vs PC3 for phocomorpha.

Additional file 10:

MultivariateIE_phocomorpha.pdf, figure, Multivariate distances calculated from nodes from IE model. Thickness of branches represents multivariate distance between ancestor–descendant.

Additional file 11:

Bayou_PC1_otarioidea.pdf, figure, Results of bayou analysis on PC1 for otarioidea. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Additional file 12:

Bayou_PC2_otarioidea.pdf, figure, Results of bayou analysis on PC2 for otarioidea. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Additional file 13:

Bayou_PC3_otarioidea.pdf, figure, Results of bayou analysis on PC3 for otarioidea. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Additional file 14:

Bayou_PC4_otarioidea.pdf, figure, Results of bayou analysis on PC4 for otarioidea. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Additional file 15:

Bayou_PC1_phocomorpha.pdf, figure, Results of bayou analysis on PC1 for phocomorpha. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Additional file 16:

Bayou_PC2_phocomorpha.pdf, figure, Results of bayou analysis on PC2 for phocomorpha. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Additional file 17:

Bayou_PC3_phocomorpha.pdf, figure, Results of bayou analysis on PC3 for phocomorpha. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Additional file 18:

Bayou_PC4_phocomorpha.pdf, figure, Results of bayou analysis on PC4 for phocomorpha. Circles at the node represent the likelihood value of a shift occurring at that node. Nodes with a likelihood value of 0.2 or greater were mapped with a regime shift, represented by a change in color, on the subsequent branches.

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Jones, K.E., Smaers, J.B. & Goswami, A. Impact of the terrestrial-aquatic transition on disparity and rates of evolution in the carnivoran skull. BMC Evol Biol 15, 8 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Disparity
  • Carnivora
  • Fissiped
  • Pinniped
  • Cranial morphology
  • Radiation
  • Ecological transition
  • Aquatic mammal