Wonky whales: the evolution of cranial asymmetry in cetaceans
BMC Biology volume 18, Article number: 86 (2020)
Unlike most mammals, toothed whale (Odontoceti) skulls lack symmetry in the nasal and facial (nasofacial) region. This asymmetry is hypothesised to relate to echolocation, which may have evolved in the earliest diverging odontocetes. Early cetaceans (whales, dolphins, and porpoises) such as archaeocetes, namely the protocetids and basilosaurids, have asymmetric rostra, but it is unclear when nasofacial asymmetry evolved during the transition from archaeocetes to modern whales. We used three-dimensional geometric morphometrics and phylogenetic comparative methods to reconstruct the evolution of asymmetry in the skulls of 162 living and extinct cetaceans over 50 million years.
In archaeocetes, we found asymmetry is prevalent in the rostrum and also in the squamosal, jugal, and orbit, possibly reflecting preservational deformation. Asymmetry in odontocetes is predominant in the nasofacial region. Mysticetes (baleen whales) show symmetry similar to terrestrial artiodactyls such as bovines. The first significant shift in asymmetry occurred in the stem odontocete family Xenorophidae during the Early Oligocene. Further increases in asymmetry occur in the physeteroids in the Late Oligocene, Squalodelphinidae and Platanistidae in the Late Oligocene/Early Miocene, and in the Monodontidae in the Late Miocene/Early Pliocene. Additional episodes of rapid change in odontocete skull asymmetry were found in the Mid-Late Oligocene, a period of rapid evolution and diversification. No high-probability increases or jumps in asymmetry were found in mysticetes or archaeocetes. Unexpectedly, no increases in asymmetry were recovered within the highly asymmetric ziphiids, which may result from the extreme, asymmetric shape of premaxillary crests in these taxa not being captured by landmarks alone.
Early ancestors of living whales had little cranial asymmetry and likely were not able to echolocate. Archaeocetes display high levels of asymmetry in the rostrum, potentially related to directional hearing, which is lost in early neocetes—the taxon including the most recent common ancestor of living cetaceans. Nasofacial asymmetry becomes a significant feature of Odontoceti skulls in the Early Oligocene, reaching its highest levels in extant taxa. Separate evolutionary regimes are reconstructed for odontocetes living in acoustically complex environments, suggesting that these niches impose strong selective pressure on echolocation ability and thus increased cranial asymmetry.
Cetaceans (whales, dolphins, and porpoises) are the most taxonomically diverse aquatic mammal clade  and inhabit most major ocean basins and some rivers . Whales appear in the fossil record approximately 52.5 million years ago (Mya), with the two extant cetacean suborders, mysticetes (baleen whales) and odontocetes (toothed whales), diverging around 39 million years ago (Mya) . Mysticetes evolved large body sizes and specialisations for bulk filter feeding whilst odontocetes evolved echolocation (biosonar) [4,5,6] and employ various raptorial and/or suction feeding strategies [7, 8]. Cetaceans have undergone extensive morphological changes to adapt to a fully aquatic lifestyle [9, 10] and show extremely divergent morphologies compared to their terrestrial artiodactyl relatives. Some of the most striking changes have occurred in the skull, including the posterior displacement of the nares, maxilla and premaxilla, and a shortening of the nasals [11,12,13].
Odontocetes are well-known to have asymmetrical crania , whereas mysticetes have bilaterally symmetrical skulls and no asymmetry in the nasofacial region . Fahlke et al.  hypothesised that basilosaurids and protocetids (early cetaceans belonging to the archaeocetes) also have cranial asymmetry thought to be linked to aquatic directional hearing with the most conspicuous asymmetry occurring in the rostrum [15, 16]. Asymmetry in odontocetes is always unidirectional, with a posterior and sinistral shift in the bones, linked to the hypertrophied melon, phonic lips, and nasal sacs, all of which are associated with high-frequency sound production and echolocation [16, 17]. Most of this asymmetry appears in the dorsal opening of the nares [14, 15, 18] and appears to be correlated with the degree of elevation in the cranial vertex . Species with high cranial vertices such as physeterids, kogiids, and ziphiids tend to have the most asymmetrical crania, likely because a functional component of asymmetry pertains to soft facial anatomy and consequently drives evolution of the underlying bony structures .
Odontocete asymmetry is thought to have evolved as a result of an evolutionary hyperallometric investment into sound-producing structures to facilitate the production of high frequency vocalisations [11, 19,20,21,22], but alternative explanations have been put forward. MacLeod et al.  proposed that skull asymmetry is a by-product of the selection pressure for an asymmetrically positioned larynx, an aquatic adaptation which enables the swallowing of large prey underwater without mastication. However, this has been argued against because reduction of tooth size and loss of shearing occlusion started after asymmetry was well developed, suggesting that swallowing prey whole may not be the driver of asymmetry . Alternatively, cranial asymmetry in basilosaurids and protocetids is thought to be linked to aquatic directional hearing . The limited or lack of asymmetry in mysticetes, which do not echolocate and instead specialise in low and infrasonic frequencies [23,24,25], suggests directional cranial asymmetry is more likely related to echolocation than hearing .
Previous studies have focused on either odontocete cranial shape and function , archaeocete asymmetry , or mysticete symmetry with modern odontocetes and archaeocetes for comparison . There is, however, little resolution on how cranial asymmetry evolved during the transition from archaeocetes to modern whales (Neoceti) , and little is known about archaeocete asymmetry and its relationship, if any, to that of odontocetes . To assess when and how often asymmetry may have arisen, where and if it is present in the archaeocete skull, and how it relates to the evolution of echolocation, it is necessary to adopt a comparative approach by broadly sampling across living and extinct cetaceans. Here we use geometric morphometric techniques to quantify asymmetry in the skull across modern and fossil species of Cetacea. We then use these data to reconstruct the evolution of asymmetry across cetaceans and test for shifts and jumps in the rate of evolution of cranial asymmetry across the cetacean phylogeny. Finally, we use these results to test potential factors associated with the evolution of asymmetry in specific cetacean clades, including presence or absence of echolocation, echolocation frequency, and inhabiting acoustically complex or high-pressure environments such as shallow rivers, cluttered icy waters, and deep ocean.
Cranial asymmetry across cetaceans
Comparing the sum radii (Σρspec) for each specimen in our data set, we found that odontocetes, especially the monodontids, physeterids, and kogiids, are the most asymmetrical of the cetaceans (Table 1).
Of 172 specimens (162 cetaceans + 10 terrestrial artiodactyls for reference), the top 43 with the highest sum radii per specimen (Σρspec ) are odontocetes. The highest ranking mysticetes are a balaenopterid (MNNZ MM001630) (Σρspec = 0.300), Aglaocetus moreni (Σρspec = 0.298), and Janjucetus hunderi (Σρspec = 0.295), ranked 47th, 50th, and 51st, respectively (Additional file 1: Table S1). The highest ranking archaeocetes are Basilosaurus isis (Σρspec = 0.308) and Zygorhiza kochii (Σρspec = 0.306) ranked 44th and 45th, respectively. The highest-ranking terrestrial artiodactyls do not appear until the 129th (Capricornis sumatrensis, Σρspec = 0.205) and 139th (Bos sp., Σρspec = 0.195) positions. The mysticetes and terrestrial artiodactyls dominate the lower end of the ranking with eight of the last ten positions occupied by extant balaenids and balaenopterids and one fossil pelocetid (Additional file 1: Table S1). For the whole cetacean data set, the most asymmetric landmarks are the nasals, the maxilla at the sutures with the nasals and premaxilla, and the posterior, dorsal premaxilla (Table 2). This distribution is heavily influenced by the odontocete sample (n = 120, 74% of cetacean specimens). For this reason, each cetacean suborder and the terrestrial artiodactyls were analysed separately. The mean total cranial radii for odontocetes is the highest of all groups at x̄ρ = 0.290 (Table 2). The most asymmetric landmarks for odontocetes are the dorsal maxilla (suture with nasal and premaxilla), nasals, and the posterior-dorsal maxilla (Table 2). Terrestrial artiodactyls have the lowest average total radii value (asymmetry) across the skull (x̄ρ = 0.171), followed closely by mysticetes (x̄ρ = 0.191). Archaeocetes showed a moderately high level of asymmetry in the skull (x̄ρ = 0.251). Cetacean subgroups differ greatly in identity of the most asymmetric landmarks and magnitude of landmark asymmetry (Table 2). For example, the most asymmetric landmark in odontocetes is the dorsal medial maxilla (suture with nasal and premaxilla) with the average sum of radii for that landmark (x̄ρland) = 0.013, whereas the most asymmetric landmark in mysticetes is the posterior ventral lateral most point of the maxilla, with x̄ρland = 0.005 (Table 2). This difference is evident when comparing average landmark asymmetry across the groups (Fig. 1).
The basilosaurid and protocetid archaeocetes show a high level of asymmetry, akin to the levels seen in fossil and extant odontocetes (Additional file 1: Table S1). The contribution of rostral landmarks to overall cranial asymmetry in these archaeocete families ranges from 13.8% in Aegyptocetus tarfa to 31.3% in Artiocetus clavis. The average amount of asymmetry concentrated in the rostrum is higher in archaeocetes (19.3%) (this includes the families Kekenodontidae, Pakicetidae, and Remingtonocetidae, which are not commonly associated with asymmetry) than in mysticetes (14.2%) and odontocetes (14.7%) (Additional file 1: Table S2–4).
As deformation during fossil preservation may create nonbiological asymmetry and previous studies suggest this may be concentrated in the rostrum of some fossil cetaceans  (see Martínez-Cáceres et al.,  and Martínez-Cáceres and de Muizon ), we also ran analyses without any fossils and without rostral landmarks. When fossils were removed, there was a decrease in asymmetry in mysticetes (x̄ρ = 0.142) (extant mysticetes: n = 12, 7% of the data set). In contrast, fossil odontocetes are more symmetrical than most extant odontocetes (Additional file 1: Table S1). For this reason, when fossil odontocetes were removed, the level of average asymmetry in the odontocete skull increased marginally (x̄ρ = 0.292) (extant odontocetes: n = 72, 44% of the data set). Excluding rostral landmarks had the most impact on archaeocetes and mysticetes, as some of the highest levels of asymmetry in those clades are found in the rostral region (Table 2; Additional file 1: Tables S2–4). However, overall, removal of the rostral landmarks had only a minor effect on results (Additional file 1: Fig. S1-S3 , Table S5b). Principal component analysis of landmark asymmetries showed that odontocetes exhibit a wide range of cranial asymmetry (Fig. 2 ) and cranial shape (Additional file 1: Fig. S4). Mysticetes and terrestrial artiodactyls overlap in asymmetry morphospace, whilst archaeocetes have a higher level of asymmetry, similar to more moderately asymmetric odontocetes (Fig. 2). See Additional file 1: Fig. S5 for identification of each specimen in the morphospace.
The basilosaurid and protocetid archaeocetes show a high amount of asymmetry (∑pspec) in the cranium (Fig. 3 (1)), similar to that observed in early odontocetes, such as xenorophids (Fig. 3 (4)). Asymmetry decreases towards the base of Neoceti, and mysticetes show the lowest level of cetacean asymmetry observed in the entire data set (Fig. 3 (2)), overlapping with terrestrial artiodactyls (Fig. 2). As expected, odontocetes show higher values of asymmetry but span nearly the full range of asymmetry morphospace (Fig. 2). The highest asymmetry is seen in the platanistids, monodontids, and physeteroids (Fig. 3 (6–8)). There are some exceptions within odontocetes, however, such as lower levels of asymmetry in the other extant river dolphins (Inia, Pontoporia, and Lipotes) [21, 31]. Lower asymmetry is also observed in the family Eurhinodelphinidae , the extant phocoenids [26, 33], and genus Sousa  (Fig. 3 (5)).
Evolutionary models of asymmetry
Reconstructing shifts in the rate of asymmetry evolution supported three shifts with a probability > 0.375, with several additional shifts at lower probabilities (Fig. 4). There is a probability (0.375) of a shift in asymmetry occurring in the family Xenorophidae during the Early Oligocene (~ 30 Mya); this represents the first, large shift in asymmetry in the cetacean phylogeny (Fig. 4). Another shift occurs in the physeteroids (~ 23 Mya; probability = 0.750), and a later shift (probability = 0.625) is seen in the Late Miocene/Early Pliocene in the Monodontidae (Fig. 4). There are two smaller shifts (probability = 0.250) in the Squalodelphinidae and Platanistidae in the Late Oligocene/Early Miocene and later in the ‘inioids’. There are no high probability shifts in asymmetry in the mysticete suborder, nor in the archaeocetes. Shifts with small (< 0.250) probabilities of occurrence are scattered throughout the phylogeny (Fig. 4), including one low probability shift at the root of Archaeoceti, but most shifts occur within the odontocetes. There is no measurable probability of a shift occurring in the archaeocete protocetids and basilosaurids. A slower or decreasing rate of asymmetry evolution is reconstructed within Mysticeti. Surprisingly, no shifts are reconstructed in the ziphiids, an odontocete family with bizarre asymmetrical premaxillary crests in most species (e.g. Ziphius cavirostris).
We found a similar pattern for ‘jumps’ (a temporary or rapid change in the trait) (Fig. 5) as we did for ‘shifts’, with the addition of several jumps occurring in the Mid-Late Oligocene. The largest jumps (probability = 0.750) occur in the physeteroids and the monodontids. Smaller jump probabilities (0.625) occur in the Delphinidae, specifically in the subfamily Globicephalinae (e.g. Globicephala spp., Pseudorca crassidens) and Platanistidae and Squalodelphinidae and also at around probability = 0.40 in the xenorophids and the kentriodontids (Fig. 5). The traces of the chains for the two models (shifts and jumps) show that a successful burn-in occurs before 25% of the model iterations are run, justifying the use of the default value (Additional file 1: Fig. S6, S7 and Model diagnostics [34,35,36]). All model diagnostics are provided in Additional file 1: Fig. S6–10; Table S6 and Model diagnostics section [34,35,36].
Evolutionary models of influence on asymmetry
Inclusion or exclusion of the rostrum made no difference to the ordering of the ‘goodness of fit’ of the models (Table 3; Additional file 1: Table S5b ). There was no difference in the ordering of the ‘goodness of fit’ of the top models when we ran all models with a phylogeny that includes only species that appear in a character matrix from Lloyd and Slater  (Additional file 1: Fig. S11-S13 Table S5c ). For this reason, the results focus on the analyses which include the rostral landmarks and the original ‘genus tree’ phylogeny.
The best fit model for both data sets is the ‘OUM-regime’ (AIC = − 448) (Table 3), which is the model with a selective regime suggesting that the monodontids, physeteroids, and platanistids are evolving under a single different regime (Table 3; Additional file 1: Table S5a ), under the assumption of Ornstein–Uhlenbeck (OU). The ‘OUM-regime-split’ model (the ‘regime’ model split into 4 separate ‘regimes’ i.e. one regime (evolutionary state) for the monodontids, one for physeteroids, and one platanistids) also received strong support (AIC = − 445). In both the ‘OUM-regime’ and ‘OUM-regime-split’ models, archaeocetes are placed into one regime, mysticetes in another, and the remaining odontocetes in a third. The third best fit model is the ‘OUM-echo-freq’ model (AIC = − 403) (Table 3), again under an OU assumption, with species categorised by their predominant echolocation/sound group.
Phylogenetic ANOVAs supported the ‘OUM-regime’, ‘OUM-regime-split’, and ‘OUM-echo-freq’ models as factors significantly associated with total cranial asymmetry (Σρspec) across cetaceans (F = 26.97, p < 0.001; F = 15.78, p < 0.001; F = 5.83, p < 0.001, respectively). Geological age, suborder, and presence/absence of echolocation were not significantly associated with cranial asymmetry (F = 1.10, p = 0.36; F = 1.57, p = 0.21; F = 1.44, p = 0.23, respectively). After correction for false discovery rate (using the Benjamini-Hochberg method [55, 56]), the regime, regime-split, and echolocation frequency models remained significant (p < 0.001, p < 0.01, p < 0.001, respectively) (Additional file 1: Table S7 [55, 56]). Hereafter, results with the Benjamini-Hochberg correction are discussed.
Our analyses of cranial asymmetry through the evolutionary history of whales suggests that the common ancestor of living whales (mysticetes and odontocetes) did not possess an asymmetric cranium, and thus, it is unlikely that echolocation was present at that stage of whale evolution or at any point in mysticete evolution. Cranial asymmetry is highest in crown odontocetes and first becomes a major feature of odontocete crania in the Early Oligocene soon after their divergence from mysticetes. This period has previously been identified as one of unusually high diversity and evolution in neocete skull morphology [13, 37, 57] alongside an explosive and rapid radiation of crown cetaceans [38, 57, 58].
Rostral asymmetry is observed in some archaeocetes and is potentially related to directional hearing, possibly increased by deformation in some cases. Fahlke et al.  suggest that Artiocetus clavis (GSP-UM 3458—the same specimen as used in this study) was found palate-up with no evident compression or deformation and further suggest that archaeocete asymmetry in the rostrum is consistent in direction. We found this same rostral asymmetry in this and other archaeocetes along with asymmetry in the jugal, orbit, and squamosal. This rostral asymmetry disappears in Neoceti and later arises in the nasofacial region in odontocetes. In archaeocetes, four of the ten most asymmetric landmarks (∑pland) were located in the rostrum (Table 2). This distribution could be inferred as torsion in the archaeocete rostrum as part of a complex of traits which led to directional hearing . This asymmetry then disappears during the transition from archaeocetes to early neocetes (Fig. 3; Additional file 1: Table S1). It is unclear whether this is due to an actual shift from a primitive form of aquatic directional hearing in specific archaeocetes (the basilosaurids and protocetids, as suggested by Fahlke et al. ) to a different regime (i.e. to high-frequency sound production in the odontocetes and low-frequency hearing in the mysticetes), or whether this is simply asymmetry unrelated to function. Asymmetry unrelated to function is reported for other mammals (e.g. dextral twist in the rostral region of some dogs ) or even brought on by developmental and environmental stressors [60, 61]. Further, it could be related to specific feeding strategies such as bottom-feeding or other lateralized behaviours. When looking at the primary landmarks displaying asymmetry in the basilosaurids and protocetids, there is no indication that these are dominated by rostral torsion more than in the other archaeocetes (Additional file 1: Table S2), and instead, asymmetry appears to be spread in no particular pattern across the jugal, squamosal (which are possibly more susceptible to deformation), rostrum, and orbit for these families. Rostral asymmetry in the archaeocetes is at least partly caused by fossil distortion in some specimens , but perhaps may also be biologically present in more archaecoete families than previously thought.
We found no high probability shifts (Fig. 4) in asymmetry occurring in the protocetids and basilosaurids, despite a rapid change from high asymmetry to a more symmetrical skull in the early mysticetes such as Coronodon havensteini (Fig. 3). We did, however, find evidence for small temporary and rapid change (jumps) in asymmetry in the later archaeocetes (Fig. 5). Echolocation, telescoping, and ecological specialisation rapidly evolved shortly after the divergence of Neoceti from Basilosauridae [4, 38], and there may have been a rapid regime change from directional hearing occurring at the same time, possibly with associated asymmetry.
Asymmetry is lowest in basal mysticetes such as Coronodon havensteini and the aetiocetids and remains low in mysticetes from the Oligocene to present. There are no high probability shifts in asymmetry in the mysticetes. Rather, Mysticeti largely display a slower or decreasing rate of the trait.
There are some increases in asymmetry observed in individual mysticetes, for example in Balaenopteridae indet NMNZ MM00163 and Aglaocetus moreni (FMNH P13407), but this likely represents taphonomic distortion in the rostrum rather than biological asymmetry. Balaenopteridae indet NMNZ MM00163 especially has some distortion in the supraoccipital, postorbital process, lateral posterior squamosal, and the parietal which likely account for its high Σρspec.
Quantifying cranial asymmetry in living and extinct mysticetes allows reconsideration of the evolution of echolocation in this clade. The consensus is that cranial asymmetry in whales evolved due to the production of high-frequency vocalisations [19,20,21]. The consistent level of symmetry in the mysticetes corroborates the hypothesis that mysticetes never evolved sophisticated echolocation [25, 62] and also contradicts the hypothesis that this suborder secondarily lost their echolocation capabilities . Our analysis further suggests that echolocation was likely not present in the common ancestor of mysticetes and odontocetes [25, 62] but evolved early in the common ancestor of odontocetes shortly after they diverged from mysticetes . As reported in Fahlke and Hampe , mysticete crania are similar in magnitude of asymmetry to terrestrial artiodactyls (Table 2; Fig. 2). In mysticetes, the highest level of cranial asymmetry was found across the rostrum (anterior and posterior maxilla and premaxilla), likely due to deformation. In some extant specimens, we observed that the tip of the rostrum has dried out and partly split apart. Even with drying-out and potential taphonomic deformation, the levels of asymmetry in mysticetes were lower than asymmetry seen in archaeocetes and much lower than that of odontocetes.
Cranial asymmetry first appears as a significant morphological trait in the Early Oligocene odontocetes Xenorophidae (Fig. 3, Fig. 4), suggesting that biosonar arose early in odontocete evolution. Odontocete asymmetry is overwhelmingly concentrated in the nasals including the posterior suture with the frontal, maxilla, and premaxilla. Most early odontocetes are less asymmetric (Fig. 3) compared to later extinct and modern forms , bar a few exceptions. The extant La Plata dolphin (Pontoporia blainvillei) is one of few living odontocetes with cranial symmetry but asymmetric nasal sacs , and it ranks here as the least asymmetric odontocete (∑pspec = 0.179, Additional file 1: Table S1). Other extant odontocetes with low cranial asymmetry include Sousa, Sotalia, and Steno (Fig. 3) which have been suggested to converge in skull morphology with kentriodontids , (Fig. 3). Phocoenids also exhibit a low level of cranial asymmetry (Fig. 3) . This low asymmetry is likely tied to their relatively low peak-power biosonar [22, 64]. Further, many descriptions of eurhinodelphinids have suggested that their crania are only slightly asymmetric [32, 65], as is supported here (Fig. 3). Thus, it should be considered that although some later fossil odontocetes had symmetrical skulls, they may have had asymmetrical nasal sacs as is observed in these extant species.
Macroevolutionary reconstruction of shifts and jumps in cranial asymmetry throughout cetacean evolution supported the first major positive shift (probability = 0.375) in asymmetry occurring in xenorophids during the Early Oligocene (~ 30 Mya) (Fig. 3; Fig. 4). This result adds further evidence to the idea that xenorophids and other odontocetes iteratively evolved specialisations for the production of high-frequency sounds necessary for echolocation [4,5,6, 39]. The distinct cranial morphology (and by inference, distinct soft tissue morphology) found in xenorophids (e.g. a deep rostral basin, a narrow premaxillary fossa, and a postnarial fossa) indicate a form of echolocation unique to the clade which interestingly, as it became more specialised, also became more asymmetrical, highlighting the importance of this trait for echolocation. The position of xenorophids as the earliest diverging clade within Odontoceti demonstrates that echolocation, telescoping, and ecological specialisation rapidly evolved shortly after the extinction of the Basilosauridae [5, 6, 38]. Since then, cranial asymmetry has increased and remained generally high throughout the odontocete lineage (Fig. 3), bar a few exceptions.
Later shifts are observed in the physeteroids in the Late Oligocene (~ 23 Mya) and in the Squalodelphinidae (this increase in asymmetry is also recently mentioned in Bianucci et al. ) and Platanistidae in the Late Oligocene/Early Miocene. The latter two families share marked asymmetry in the premaxillae with the right maxilla narrower than the left in dorsal view . Further, asymmetry is recorded in the frontal and maxillary crests of fossil platanistids such as Zarhachis flagellator  although the supraorbital crests are not as developed as the extreme maxillary crests in the extant Platanista gangetica which is one of the most asymmetric of all odontocete skulls . There is also marked skull asymmetry in the distantly related squalodelphinid, Notocetus vanbenedeni, which also sits within the superfamily, Platanistoidea [66, 67].
The best model fit, the regime model (p < 0.001) (model: 'regime'), assumes there is a distinct evolutionary regime for the most asymmetrical odontocete specimens (physeteroids, platanistids, and monodontids) indicating a single driver for their extreme asymmetry. We hypothesise that this regime may be linked to the pressures which arise from inhabiting acoustically complex environmental niches. The physeteroids were the first of the major odontocete crown lineages to rapidly diverge and are easily recognisable due to a highly asymmetric facial region and supracranial basin . Their large body size and hypertrophied nasal structures produce a low-frequency multi-pulsed sound , which facilitates long range detection of prey . This is highly advantageous when searching for patchy prey, especially as the physical properties of the water itself alter sound velocity and potentially constrain sensory morphology .
Platanista gangetica, the sole modern survivor of Platanistidae sits alone among river dolphins for having a highly asymmetric cranium and echolocating at broadband low frequency (BBLF). The unique, autapomorphic bony maxillary crests of Platanista may help achieve a higher directionality than expected for a cetacean that clicks nearly an octave lower than similar sized odontocetes , a feature that would be useful in the turbid, cluttered rivers they inhabit. Other species in this highly asymmetric model include both monodontids: belugas (Delphinapterus leucas) and the narwhal (Monodon monoceros). Monodon remains the most asymmetric skull in the sample, even when the rostrum is removed (Σρspec = 0.472) which rules out the possibility that the asymmetric tusk and residual teeth may be skewing the overall Σρspec (see below for details). Their unique sound repertoire (narrowband structured, NBS) is ideal for projecting and receiving signals in icy, shallow waters, where the animals can detect targets in high levels of ambient noise and backscatter  (Additional file 1: Table S8 [24, 40,41,42, 45,46,47,48,49,50,51,52,53,54, 64, 70,71,72]). Jumps detected in the delphinids all belong to the subfamily Globicephalinae (Fig. 5). In particular, the highly asymmetrical Globicephala (Table 1; Additional file 1: Table S1) has evolved a deep-dive pattern to target a deep-water niche occupied by large, calorific, and fast squid, and its acoustic behaviour is more akin to deep divers than to oceanic delphinids . The cochlea of Globicephala is also morphologically different to other delphinids , which could also represent adaptation to the extreme acoustic environment of the deep ocean. Further, studies suggest that Pseudorca (which also has a highly asymmetric cranium (Table 1)) echolocates with different vertical and horizontal plane patterns to other delphinids .
Surprisingly, no jumps or shifts are seen in the deep-diving ziphiids (beaked whales), an odontocete family with bizarre asymmetrical premaxillary crests and an asymmetric prenarial basin (Additional file 1: Fig. S14. The asymmetry of the beaked whale skull is marked [20, 75, 76], so much so that the right premaxilla, premaxillary crest, premaxillary sac fossa, and the nasal bone are around 30% larger than those on the left . Previous studies have suggested that the beaked whale genus Berardius (the most basal crown genus) shows the least bilateral asymmetry in the skull [78, 79], and we saw a similar result here. We attribute the underrepresentation of asymmetry in the ziphiid skull to the use of landmarks alone. Whilst detecting asymmetry in the shifting of the nasal, premaxilla, and maxilla to the left side of the skull, this method underrepresents the degree of asymmetry in the morphology of the bones themselves. The premaxilla is landmarked with points at the posterior dorsal premaxilla and the dorsal medial maxilla (suture with nasal and premaxilla) which accurately captures asymmetry in the positioning of the bone and its attachment but fails to capture the tapering of the highly asymmetric premaxillary crest itself (Additional file 1: Fig. S14). Future studies in this area should be done with curve sliding semi-landmarks and surface patches to more accurately capture the complex morphology  of the premaxillary crests and premaxillary sac fossae in ziphiids which are not represented using fixed landmarks alone.
‘Regime-split’ (model: ‘regime-split’) was the second-best model fit which had a significant effect (p < 0.01) on asymmetry in the cranium. This model suggests a different evolutionary regime for each of the most asymmetric groups. As above, it could be hypothesised that the highly asymmetric species live in unique, acoustically complex environments all of which have rather extreme specific environmental selection pressures. The reduction in the p value after phylogenetic correction for the regime and regime-split models suggests that the factors influencing asymmetry may be shared by closely related taxa.
Frequency of echolocation (model: ‘echo-freq’) also had a significant effect (p < 0.001) on the asymmetry in the cranium and was the third best model fit. Echolocation frequency has been widely suggested as a key driver of asymmetry in the cranium [16, 17] and soft tissues . Although not the best model fit, we suggest that this relationship be investigated in more detail, for example with a more detailed analysis of species-specific echolocation frequencies and associated categories across Cetacea . It is important to note that these methods assume a Brownian motion model, which oversimplifies the actual evolutionary model underlying the evolution of asymmetry (shown here to be better described by an OU model).
We found no support for several other potential drivers for observed patterns of cranial asymmetry, independent of phylogeny. There is no significant effect of geologic age of the specimen (e.g. Eocene, Miocene, extant) on sum radii in the skull (p = 0.36). This result is likely because, despite odontocete crania becoming more asymmetrical in most extant families, mysticetes do not. There is no significant effect of ‘suborder’ (p = 0.21) on the total sum radii across the cranium. This is not surprising as there is generally a clear phylogenetic relatedness in whether a cetacean is symmetric (mysticete) or asymmetric (odontocete). Presence or absence of echolocation (model: ‘echo’) has no significant effect (p = 0.23) on the sum of radii in the cranium. Again, this is not surprising as there is a clear phylogenetic relatedness in whether a cetacean can echolocate, i.e. the odontocete suborder, or not echolocate, i.e. the mysticete suborder.
There is a small chance that skulls used in this study may be more asymmetrical, i.e. deformed or distorted, than a standard skull of the species and therefore this is represented in the placed landmarks and the resulting ∑pspec. Where possible, we chose skulls based on their overall quality and representation of the species. This was not possible for fossils which are often represented by one specimen, but deformed skulls were removed from the study so as not to falsely imply there is biological asymmetry in the skull when there is none. Further, the sex of the specimen may slightly alter the degree of asymmetry in the skull. Female false killer whales, for example, have a slightly more asymmetrical skull than males , and this may partially explain why the individual in this study appears to have a higher level of asymmetry than the other delphinids. However, the sex of this specimen (USNM 11320) is listed as unknown. It is important to note that adult male narwhal exhibit an extreme form of asymmetry in the tusk and vestigial teeth . The specimen in this study (USNM 267959) is female and therefore lacks a highly asymmetric tusk, however, the paired tusks embedded in the maxillae may still exhibit asymmetry  and may affect the overlying bone structure. This has not skewed the results seen here as the top 6 landmarks of asymmetry in the Monodon skull are in the nasals and posterior premaxilla and maxilla (i.e. not the rostrum or anterior maxilla where imbedded tusks reside). Further, no landmarks were placed on tusks or teeth (see the “Methods” section: “Manually placed landmarks (Fn)”), which ensures that extreme asymmetry seen in some tusked species for example, Odobenocetops, is not captured in this study. The ‘skew’ in Globicephala and monodontid skulls has also been attributed to some asymmetry in the attachment of the neck muscles ; however, the asymmetry (Σρland) in eight landmarks associated with the condyle and posterior cranium do not differ between Globicephala and the monodontids compared to other closely related species (e.g. Feresa attenuata and Peponocephala electra).
Lastly, an argument against the hypothesis that echolocation drives asymmetry in the odontocete skull is that bats also echolocate and do not have cranial asymmetry as the natural condition . However, the extreme differences in the environments in which bats and cetaceans echolocate, as well as other ecological and morphological differences between the two clades, complicate any meaningful comparison . It should be noted that both odontocetes and bats share a remarkable convergence on narrow biosonar beams across species independent of body size [22, 86], with the ability to do this in odontocetes likely a result of cranial asymmetry.
With the most widely supported explanation of asymmetry being sound production, our results support the hypothesis that craniofacial asymmetry (along with concavity in the facial area, hypertrophied naso-facial muscles, air sacs, melon, and premaxillary sac fossa ) arose in odontocetes to support high-frequency echolocation. Further, echolocating in complex environments continues to be a primary factor driving the evolution of asymmetry in the odontocete skull, as supported by the independent evolutionary regimes for the most asymmetric odontocetes.
Our study represents the first comprehensive analysis of cranial asymmetry spanning the evolutionary history of cetaceans. We demonstrate that the common ancestor of living cetaceans had little cranial asymmetry and thus is unlikely to have possessed the ability to echolocate. Odontocetes display increasing cranial asymmetry from the Oligocene to present, reaching their highest levels in extant taxa. Separate evolutionary regimes are supported for three odontocete clades (monodontids, physeteroids, and platanistids) that inhabit acoustically complex environments, suggesting that echolocation and cranial asymmetry are continuing to evolve under strong selection in these niches. Surprisingly, no increases in asymmetry were recovered within the highly asymmetric ziphiids. We attribute this to the extreme, asymmetric shape of the premaxillary crests and sac fossae in these taxa not being captured by landmarks alone.
Mysticetes have maintained a low level of cranial asymmetry since their origin, and if asymmetry reflects ultrasonic sound production ability, it is unlikely that mysticetes were ever able to echolocate. Archaeocetes have a high level of asymmetry in the rostrum which could be linked to directional hearing, as reported by Fahlke et al. , but this rostral asymmetry disappears in early neocetes as the dichotomous hearing abilities of the two suborders became established.
Modelling the evolution of cranial asymmetry across living and extinct cetaceans recovered the highest probabilities of shifts in the trait at three main points: first, in the extinct odontocete xenorophids in the Early-Mid Oligocene, then in the physeteroids (Late Oligocene), and finally in the monodontids in the Late Miocene/Early Pliocene. Smaller shifts were found in the Squalodelphinidae and Platanistidae. This was also true for ‘jumps’ in the trait, with an additional jump in a branch of the delphinids (namely the Globicephalinae e.g. pilot whales and false killer whales). Additional episodes of rapid change were found in the Mid-Late Oligocene, a period of rapid evolution in cranial asymmetry in odontocetes. These results support studies suggesting that biosonar, the signature adaptation of odontocetes, and associated asymmetry were acquired at or soon after the origin of this clade [4,5,6, 39].
The data set comprises stem cetaceans (archaeocetes, n = 10) and both extant suborders: the baleen whales (mysticetes, n = 32) and toothed whales (odontocetes, n = 120). The final data set comprised 162 cetacean crania, of which 78 (48%) are extinct, ranging in age from 48.6 to 2.59 Mya. Additionally, 10 terrestrial artiodactyls (representing 7 of the 10 Arctiodactyla families) were included to provide a baseline for symmetry as cetaceans are nested within Artiodactyla. Specimen details (Additional file 1: Table S9) and museum abbreviations are provided in Additional file 1.
Specimens were selected to cover the widest possible phylogenetic spread, representing 38 families and 101 genera from the Eocene to the present.
The Early-Middle Eocene is represented by the land-dwelling family Pakicetidae through to semi-aquatic Ambulocetidae and Remingtonocetidae. The Pelagiceti are represented by the fully aquatic Basilosauridae of the Late Eocene through to the modern Neoceti. This includes representation of some early stem toothed mysticetes such as the Mammalodontidae and the Aetiocetidae. Three of the four extant mysticete families are represented. The odontocetes are represented by early stem families: the Xenorophidae and the Simocetidae of the Early-Mid Oligocene and the ‘Patriocetidae’ (phylogenetic position is still being clarified) of the Late Oligocene. The more crownward odontocetes of the Miocene are represented by the Eurhinodelphinidae, Kentriodontidae, Albireonidae, Squalodelphinidae, Squalodontidae, and Allodelphinidae among other extinct families. All ten extant odontocete families are represented. See Additional file 1: Table S9 for details.
Because many extant and all fossil specimens lack information on sex, sexual dimorphism could not be considered. All specimens are adult except for one subadult, Mesoplodon traversii. Specimens were selected based on completeness but some bones were broken (e.g. jugal) and were treated as missing data. Sixty-four (~ 39%) of the specimens, including some extant specimens, had missing data concentrated in the pterygoid, palatine, jugal, squamosal, and tip of the rostrum. For this reason and because fossils often have a higher proportion of missing data, we also ran analyses without any fossils and without rostral landmarks. Specimens with obvious taphonomic or other deformation were excluded from further analysis (Additional file 1: Table S10). Excluded specimens include the basilosaurid Cynthiacetus peruvianus which shows sinistral torsion in the rostrum. Although a potential natural feature in protocetids and basilosaurids [15, 16], it is suggested that rostral distortion in this particular specimen (MNHN.F.PRU10) is at least partly the original morphology of the skull and potentially a result of some taphonomic distortion . Inevitably, some fossil specimens have sections of reconstructed bone. Their inclusion in the study was based upon the extent and accuracy of the reconstruction and the unavailability of alternative specimens.
Skulls were scanned using a Creaform Go!SCAN 20, or Creaform Go!SCAN 50 depending on the size of the skull. Scans were cleaned, prepared, and merged in VX Elements v.6.0 and exported in ply format before being further cleaned and decimated in Geomagic Wrap software (3D Systems). Models were decimated down to 1,500,000 triangles, reducing computational demands without compromising on detail for further morphometric analyses. In many studies of morphology when the skull is incomplete, it is possible to digitally reconstruct bilateral elements by mirroring across the midline plane if preserved on one side [87,88,89]. However, due to the substantial asymmetry observed in many taxa in this study, mirroring a complete half of the skull was not possible (Fig. 6; Additional file 1: Fig. S15). For this reason, we limited mirroring to marginally damaged bones or easily mirrored missing bones only, where it was clear that mirroring would not mask any biological asymmetry, using the ‘mirror’ function in Geomagic Wrap (3D Systems).
Morphometric data collection
Manually placed landmarks (Fn)
We placed 123 anatomically defined landmarks over the surface of the skull in Stratovan Checkpoint (Stratovan, Davis, CA, USA) using the ‘single point’ option. We placed 57 landmarks on both the left-hand side (LHS) and right-hand side (RHS) of the skull, and 9 landmarks on the midline, totalling 123 landmarks covering both the dorsal and ventral sides of the skull (Fig. 7). Type I and II landmarks  were selected to comprehensively represent the full cranium (Fig. 7; Additional file 1: Table S11). ‘Landmark 15’ and the subsequent mirrored ‘landmark 79’ denote the back of the toothrow in most species. In some ziphiids, e.g. Mesoplodon carlhubbsi, the teeth (or tusks) erupt midway along the mandible  whilst other species present multiple pairs of tusks . In others (e.g. Hyperoodon ampullatus), teeth typically erupt as a single pair on the anterior mandible which often protrudes beyond the upper jaw . Without the mandible, it is challenging to pinpoint the positioning of the back of the toothrow, and even then, the presence and number of teeth is negligible in some species. Further, these tusks only erupt in adult males. For these reasons, and to avoid simply estimating where the true tooth row may be, ‘landmark 15’ and ‘landmark 79’ in specimens with mandibular prognathism, absent, maxillary-only, or vestigial dentition (including all ziphiids, narwhals (Monodon monoceros) and sperm whales (Physeter macrocephalus)) were consistently placed on the proximal lateral maxilla where the posterior end of a standard tooth row would be located (Additional file 1: Fig. S16).
As previously noted, some specimens have missing data. Geometric morphometric analyses and plotting functions implemented in geomorph v.3.1.0  require a full complement of landmarks . This complement can consist of actual landmarks and estimated positions for ‘missing’ landmarks. To estimate positions for missing landmarks, we placed ‘missing’ landmarks as close to the missing bone (areas that could not be digitally mirrored) as possible and then marked it as a ‘missing landmark’ in Checkpoint which automatically assigns a coordinate of − 9999. We then used the estimate.missing function in geomorph and the ‘TPS’ (thin plate spline) method to estimate the location of landmarks on incomplete specimens. A reference specimen which has a complete complement of landmarks is selected, and incomplete specimens are aligned against it using common landmarks . In a TPS-based estimation, missing landmarks are placed so that the overall bending energy between the reference and the incomplete specimen is smallest which creates a smooth deformation . TPS was chosen over regression-based methods (i.e. ‘Reg’ in geomorph) because it performs better in simulations with missing data .
Our study uses a phylogenetic framework to reconstruct macroevolutionary patterns of cranial asymmetry across Cetacea. To generate a tree that included all of our sampled taxa, we used the time-calibrated phylogeny from Lloyd and Slater . This ‘genus tree’ includes all species belonging to a genus that appear in a character matrix using taxonomic constraints to place taxa that lack data. We modified it as follows: First, we added several additional extant species (which were already represented to the genus level in the Lloyd and Slater phylogeny ) with position based on recently published studies. We placed Neophocaena asiaeorientalis in the same genus as Neophocaena phocaenoides , Sousa plumbea + Sousa teuszii + Sousa sahulensis in the same genus as Sousa chinensis , Orcaella heinsohni in the same genus as Orcaella brevirostris [96, 97], and Mesoplodon hotaula in the genus Mesoplodon next to Mesoplodon gingkodens . Finally, we placed Berardius minimus in the genus Berardius next to Berardius bairdii and Berardius arnuxii following its recent description by Yamada et al. . The following fossil species were directly swapped with their corresponding monophyletic congener as follows. We placed Balaneoptera sp. (SDNHM 83695) as a sister taxon to Balaenoptera siberi (although not present in our sample), close to extant Megaptera novaeangliae as in Martin , Balaenoptera floridana as a sister taxon to Balaenoptera davidsonnii [100, 101] (again, the latter species is not present in our sample), and Orycterocetus crocodilinus is placed in the physeterids according to Lambert et al. . We placed Globicephala sp. as a sister taxon to Globicephala etruriae [103,104,105] and Hemisyntrachelus cortesii in the same genus as Hemisyntrachelus oligodon according to Post and Bosselaers . We caution that Kentriodontidae is often considered a non-monophyletic ‘waste-basket’ for Late Oligocene and Miocene homodont odontocetes . Restrictions according to Peredo et al.  leave Tagicetus and Atocetus (previously referred to as Kentriodontidae) outside of the family (Additional file 1: Table S9). The positioning of Argyrocetus joaquinensis is also unclear . Two specimens (Xenorophus ChM PV7677 and Patriocetid or Waipatiid CCNHM 1078) were excluded from the analysis due to uncertainty in their position (Additional file 1: Table S10).
We generated mirrored landmarks for the right-hand side (RHS) of the skull and compared their positions to those of the original manually placed landmarks, measuring the amount of landmark displacement between the two. To do so, we used the 57 LHS landmarks and 9 midline landmarks (total = 66) (Additional file 1: Table S11) and mirrored the LHS landmarks onto the RHS using the mirrorfill function in the R package paleomorph v.0.1.4 . Before carrying out further analyses, we superimposed the specimens to remove all non-shape elements, i.e. size (scaling), translation, and rotation (positioning) from the data using Generalized Procrustes Analysis implemented in the gpagen function from the geomorph R package v.3.1.0 .
We used the R package landvR v.0.4  to calculate the Euclidean distances between a reference specimen (the computer -mirrored, landmarked specimen) (Rn) and a focal specimen (the manually landmarked specimen) (Fn). Both Rn and Fn are defined by three coordinates (x, y, z). The landmark displacements were measured for each landmark individually using the spherical coordinates system which measures between the nth landmark of the Fn and the Rn specimens respectively . This method provides 3 outputs (from Guillerme and Weisbecker ):
ρ, the Euclidean distance between Fn and Rn
ϕ, the azimuth angle formed by the projection of Rn on the equatorial plane (f(x) = 0)
θ (when using 3D data only), the polar angle formed by the projection Rn on the polar plane (f(y) = 0)
We estimated differences between Fn and Rn in the spherical coordinates system using the coordinates.difference function in landvR and extracted the ρ (radius) for each landmark, for each specimen. This provides a measure of the Euclidean distance between a manually placed landmark which accurately represented the specimen’s morphology (Fn) and a computer -mirrored landmark (Rn). If the specimen is asymmetric, the computer -mirrored landmark does not accurately reflect its morphology (Fig. 6).
The spherical coordinates system is preferable because it directly measures landmark displacement in any direction, and further, the values for each landmark displacement is discrete in space (i.e. independent from other landmarks) . This is important because it allows identification of asymmetry which may occur in discrete parts of the skull, e.g. the posterior nasal without blanket labelling all landmarks as ‘asymmetric’. We obtained ρ for each of the 123 landmarks for each specimen, including the terrestrial artiodactyls (21,156 radii values in total). The larger the radii (and consequently the larger the ϕ for a corresponding landmark) the more displacement between Fn and Rn. We interpret a higher ρ as an indication of more asymmetry in the skull (see Fig. 8 for a visualisation of this). Higher displacement means that there is a greater difference between the placement of Fn and Rn, indicating asymmetry in those landmarks . The closer the radius to 0, the more symmetrical the specimen as Fn has not displaced far from Rn. We took the averaged sum of radii for each landmark (x̄ρland) to find the most asymmetrical landmarks and identify their location on the skull for each group (archaeocetes, odontocetes, mysticetes, and terrestrial artiodactyls) (Table 2), as well as an average total cranial asymmetry (x̄ρ) for each group. We also took the sum radii for each of the individual specimens (Σρspec). We ran a principal component analysis (PCA) on the Procrustes aligned data using the ‘factoextra’ package v.1.0.7  in R to identify PC scores of maximum radii variation using the sum radii for each of the individual specimens (Σρspec) (Fig. 2). We then used overall asymmetry in each specimen (Σρspec) to reconstruct the evolution of cranial asymmetry.
For computational purposes, all polytomies in the tree were resolved by adding zero branch lengths using multi2di in ape v.5.0  prior to downstream phylogenetic analyses. All models were run with cetacean data only (i.e. no terrestrial artiodactyls), and model fit was assessed using AIC. We also conducted the analyses after removing the rostrum (NR) as it may be more easily deformed through both taphonomic deformation and drying out in extant specimens. We also ran all models with a phylogeny that includes only taxa that appear in a character matrix from Lloyd and Slater . We found no differences in the ordering of the best-fitting models with NR, nor with the phylogeny which only uses taxa that appear in the character matrix (see Additional file 1: Figs. 1-3, 11–13; Tables S5a-c); we thus focus the analyses on data including the rostral landmarks and the original ‘genus tree’.
Modelling the evolution of cranial asymmetry
We assessed how asymmetry has evolved across Cetacea using phylogenetic models of trait evolution. We investigated variation in rates of cranial asymmetry evolution using a relaxed Brownian motion (BM) process with the rjmcmc.bm function implemented in the R package geiger v.184.108.40.206 [114, 115]. This model uses a Bayesian framework with a reversible jump Markov Chain Monte Carlo (rjMCMC) algorithm to find the position and the amplitude of evolutionary rates changes and ‘shifts’ across the tree . ‘Jumps’ indicate a temporary or rapid change in the trait.
We ran the rjMCMC chain for 106 generations sampling each 10,000 generations with the combined ‘jump-rbm’ class model . We used the weakly informative Half-Cauchy distribution with scale parameter 25  as the prior density of the rate scalar and measurement error instead of the default exponential distribution and used default priors for the number of shifts (a Poisson distribution with mean equal to log(2) which places a 50% probability on a scenario with no shifts). For comparison and as a proposal mechanism for exploring the parameter space, we ran the same model with the ‘rbm’ class model and no jumps. We checked the effective sample size (ESS) and assessed convergence of the chains with Gelman and Rubin’s diagnostics [35, 36] using the ‘effectiveSize’ and ‘gelman.diag’ functions implemented in the R package coda  (Additional file 1: Figs. S6–10; Table S6 ).
Hypothesised evolutionary regimes for cranial asymmetry
Several state-dependent models were proposed as potential predictors for the level of asymmetry seen in the cetacean skull. For example, ‘echolocation’ (model: ‘echo’) (Table 4) is one model used to investigate whether the rate of evolution for skull asymmetry differs between species that can echolocate and those that cannot. We name two other models, the regime model (model: ‘regime’) and the ‘regime split’ model (model: ‘regime-split’). In these models we test whether evolutionary changes in asymmetry in the cetacean cranium (the studied trait) can be associated with the states of another discrete trait. By regime, we mean a particular condition or process that may be underlying the observed patterns of cranial asymmetry. We further fitted a ‘frequency echolocation’ model (model: ‘echo-freq’) (Table 4; Additional file 1: Table S8 [40,41,42, 45,46,47,48,49,50,51,52,53,54, 64, 70,71,72]). This study is not a specialist analysis of acoustics (nor behaviours affecting acoustics) in cetaceans, and we use these values to indicate potential drivers for the evolution of cranial asymmetry. We assigned species to several categories depending on how they predominantly produce sound (Table 4; Additional file 1: Table S8 [40,41,42, 45,46,47,48,49,50,51,52,53,54, 64, 70,71,72]). We fitted our several state-dependent models with different variants of the BM and OU model (Table 5). The multiple models and different rates are summarised below.
In addition, we evaluated the fit of 24 alternative models (all listed in Table 5) based on the states of a discrete character (Table 4) implemented in the mvMORPH package in R  using a maximum likelihood inference. We used the ‘fitDiscrete’ function in ‘geiger’ v.1.3-1 to fit various likelihood models for discrete character evolution. The model arguments tested were an ‘equal-rates’ model (ER) where all transitions occur at equal rates, a ‘symmetric transitions are equal’ model (SYM), and an ‘all rates different’ model (ARD) where each rate is a separate parameter [114, 117]. The ER model gave the best fit (Additional file 1: Table S12 [114, 117]) and was thus used in all of our alternative models using maximum-likelihood inference (Table 5).
The multiple models described above (Table 4) relax the assumption of a common dynamic for modelling the trait evolution by allowing the estimation of the model parameters that depend on the states of a discrete character. For these, we first had to ‘paint’ the evolutionary history (e.g. ancestral state) of the selective regime onto the tree. To do this, we used write.simmap in the phytools package v.0.6–99 . We ran these models under OUM and BMM assumptions (Table 5). We repeated these analyses for the data with the rostrum removed (Table 3) and with the phylogeny that includes only taxa that appear in a character matrix  (Additional file 1: Table S5b-c ). All analyses were done in R v.3.5.0 .
Although the relaxed BM process described above is very flexible and allows the investigation of changes in evolutionary rates across the tree without strong a priori, it is however limited for assessing and interpreting changes in evolutionary modes. Moreover, recovered changes in rates might result from long-term trends in the average asymmetry rather than actual changes in the pace of evolution. We considered multiple models including parameterizations of the BM and OU process (Table 5). Our models assume constant dynamics of trait evolution but a directional drift of the clade average value that might be interpreted as shifts in evolutionary rates in the relaxed BM model considered above. We also considered models with specific optimums (model “OUM” in ‘mvOU’), and ancestral states and/or rates (model “BMM” in ‘mvBM’) in different parts of the tree. The more parameterized and refined models allow for testing of evolutionary changes in the studied trait and can be associated with the states of another discrete trait. In this study, we consider different scenarios to assess whether the evolution of the skull asymmetry shows marked differences between the three major clades (archaeocetes, odontocetes, and mysticetes) and if it is related to the evolution of echolocation.
Analysis of variance
Lastly, we ran phylogenetically corrected ANOVAs on each of the different scenarios using the R package nlme (v.3.1-137)  and function ‘gls’ to test for correlations between the level of asymmetry seen in the skull and the potential scenarios (or regimes) hypothesised above (Table 4). ‘gls’ allows for a more flexible model with better power. Simulations were run using a ‘Pagel’s Lambda’ (λ) correlation structure (corPagel) in the ape package  by multiplying the off-diagonal of the phylogenetic covariance matrix by "lambda" constrained to be within [0, 1]. As we ran multiple models, we controlled for a false discovery rate using the Benjamini-Hochberg method [55, 56] (Additional file 1: Table S7 [55, 56]). The data sets generated and/or analysed during the current study are available in the Github repository: https://github.com/EllenJCoombs/Asymmetry-evolution-cetaceans .
Availability of data and materials
The data sets generated and/or analysed during the current study are available in the Github repository: https://github.com/EllenJCoombs/Asymmetry-evolution-cetaceans . All code and raw data are available to download, complete with MIT licence. Please cite this paper and the Zenodo doi: https://doi.org/10.5281/zenodo.3893943 when using the data or raw code. See Github for details.
Akaike information criterion
Analysis of variance
‘All rates different’ model
Frequency-modulated buzz clicks
- Fn :
Manually placed landmark
Million years ago
- p :
Reversible jump Markov Chain Monte Carlo
- Rn :
Computer -mirrored landmark
‘Symmetric transitions are equal’ model
- x̄ρ :
Average of the total radii (Σρ) values across the skull
- x̄ρland :
Averaged sum of radii for each landmark
- Σρspec :
Sum radii for each of the individual specimens
- Σρland :
Sum radii for specific landmarks
Fordyce RE, Barnes LG. The evolutionary history of whales and dolphins. Annu Rev Earth Planet Sci. 1994;22:419–55.
McGowen MR, Spaulding M, Gatesy J. Divergence date estimation and a comprehensive molecular tree of extant cetaceans. Mol Phylogenet Evol. 2009;53:891–906.
Marx FG, Fordyce RE. Baleen boom and bust: a synthesis of mysticete phylogeny, diversity and disparity. R Soc Open Sci. 2015;2:140434.
Geisler JH, Colbert MW, Carew JL. A new fossil species supports an early origin for toothed whale echolocation. Nature. 2014;508(7496):383–6.
Churchill M, Martínez-Cáceres M, de Muizon C, Mnieckowski J, Geisler J. The origin of high-frequency hearing in whales. Curr Biol Cell Press. 2016;26(16):2144–9.
Park T, Fitzgerald EMG, Evans AR. Ultrasonic hearing and echolocation in the earliest toothed whales. Biol Lett. 2016;12:20160060.
Hocking DP, Marx FG, Park T, Fitzgerald EM, Evans AR. A behavioural framework for the evolution of feeding in predatory aquatic mammals. Proc R Soc B Biol Sci. 2017;284(1850):20162750.
Hocking DP, Marx FG, Park T, Fitzgerald EM, Evans AR. Reply to comment by Kienle et al. 2017. Proc R Soc B Biol Sci. 2017;284(1863):20171836.
Fordyce RE, de Muizon C. Evolutionary history of cetaceans: a review. In: Secondary adaptation of tetrapods to life in water; 2001. p. 169–233.
Uhen MD. Evolution of marine mammals: Back to the sea after 300 million years. Anat Rec. 2007;290(6):514–22.
Heyning JE, Mead JG. Evolution of the nasal anatomy of cetaceans, in sensory abilities of cetaceans. US: Springer; 1990. p. 67–79.
Klima M. Development of the cetacean nasal skull. Adv Anat Embryol Cell Biol. 1999;149:1–143.
Churchill M, Geisler J, Beatty B, Goswami A. Evolution of cranial telescoping in echolocating whales (Cetacea: Odontoceti), evolution. Soc Stud Evol. 2018;72(5):1092–108.
Thompson P. The bottlenose dolphin. Trends Ecol Evol. 1990;5(11):390 Edited by Leatherwood. S and Reeves. R. Elsevier Science.
Fahlke JM, Hampe O. Cranial symmetry in baleen whales (Cetacea, Mysticeti) and the occurrence of cranial asymmetry throughout cetacean evolution. Naturwissenschaften. 2015;102(9-10):58.
Fahlke JM, Gingerich PD, Welsh RC, Wood AR. Cranial asymmetry in Eocene archaeocete whales and the evolution of directional hearing in water. Proc Natl Acad Sci U S A. 2011;108(35):14545–8.
Cranford TW, Amundin M, Norris KS. Functional morphology and homology in the odontocete nasal complex: implications for sound generation. J Morphol. 1996;228(3):223–85.
Macleod CD, Reidenberg JS, Weller M, Santos MB, et al. Breaking symmetry: the marine environment, prey size, and the evolution of asymmetry in cetacean skulls. Anat Rec. 2007;290(6):539–45.
Mead JG. Anatomy of the external nasal passages and facial complex in the Delphinidae (Mammalia: Cetacea). Smithsonian Contrib Zool. 1975;207:1–35.
Heyning JE. Comparative facial anatomy of beaked whales (Ziphiidae) and a systematic revision among the families of extant Odontoceti. Contrib Sci Na. 1989;405:1–64.
Huggenberger S, Leidenberger S, Oelschläger HHA. Asymmetry of the nasofacial skull in toothed whales (Odontoceti). J Zool. 2017;302(1):15–23.
Jensen FH, Johnson M, Ladegaard M, Wisniewska DM, Madsen PT. Narrow acoustic field of view drives frequency scaling in toothed whale biosonar. Curr Biol. 2018;28(23):3878–85.
Ketten DR. Structure and function in whale ears. Bioacoustics. 1997;8(1–2):103–35.
Reidenberg JS, Laitman JT. Discovery of a low frequency sound source in mysticeti (baleen whales): anatomical establishment of a vocal fold homolog. Anat Rec. 2007;290(6):745–59.
Park T, Evans AR, Gallagher SJ, Fitzgerald EMG. Low-frequency hearing preceded the evolution of giant body size and filter feeding in baleen whales. Proc R Soc B. 2017;284:20162528.
Marx FG, Lambert O, Uhen MD. Major Steps in the Evolution of Cetaceans, in: Cetacean Paleobiology. Wiley. 2016:157–97.
Martínez-Cáceres M, Lambert O, de Muizon C. The anatomy and phylogenetic affinities of Cynthiacetus peruvianus, a large Dorudon-like basilosaurid (Cetacea, Mammalia) from the late Eocene of Peru. Geodiversitas. 2017;39(1):7–163.
Martínez-Cáceres M, de Muizon C. A new basilosaurid (Cetacea, Pelagiceti) from the Late Eocene to Early Oligocene Otuma Formation of Peru. Comptes Rendus - Palevol. 2011;10(7):517–26 No longer published by Elsevier.
Lloyd GT and Slater GJ (in prep). A Total-Group Phylogenetic Metatree for Cetacea and the Importance of Fossil Data in Diversification Analyses. 2020. bioRxiv 2020.06.24.169078.
Creative Commons, 2020. Attribution-ShareAlike 3.0 Unported (CC BY-SA 3.0). Available at: https://creativecommons.org/licenses/by-sa/3.0/* [Accessed 06 Mar 2020].
Huggenberger S, Vogl TJ, Oelschlager HHA. Epicranial complex of the La Plata dolphin (Pontoporia blainvillei): topographical and functional implications. Mar Mamm Sci. 2010;26:471–81.
Lambert O. Phylogenetic affinities of the long-snouted dolphin Eurhinodelphis (Cetacea, Odontoceti) from the Miocene of Antwerp, Belgium. Palaeontol. 2005;48(3):653–79.
Racicot R. Dolphins, porpoises, and monodontids, Evolution. In: Wursig B, Thewissen J, Kovacs K, editors. Encyclopedia of Marine Mammals. San Diego: Academic Press/Elsevier; 2018. p. 271–4.
Plummer M, Best N, Cowles K, Vines K. CODA: convergence diagnosis and output analysis for MCMC. R News. 2006;6:7–1.
Gelman A. Prior distributions for variance parameters in hierarchical models. Bayesian Anal. 2006;1(3):515–33.
Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stati Sci Inst Math Stat. 1992;7(4):457–72.
Berta A, Sumich JL, Kovacs KM. Marine Mammals: Evolutionary Biology. 2015. Academic Press Inc; 3rd Revised edition edition (5 Jun. 2015).
Boessenecker RW, Fraser D, Churchill M, Geisler JH. A toothless dwarf dolphin (Odontoceti: Xenorophidae) points to explosive feeding diversification of modern whales (Neoceti). Proc R Soc B. 2017;284:20170531.
Racicot RA, Boessenecker RW, Darroch SAF, Geisler JH. Evidence for convergent evolution of ultrasonic hearing in toothed whales (Cetacea: Odontoceti). Biol Lett. 2019;15:20190083.
Kastelein RA, Bunskoek P, Hagedoorn M, Au WWL, et al. Audiogram of a harbor porpoise (Phocoena phocoena) measured with narrow-band frequency-modulated signals. J Acoust Soc Am. 2002;112(1):334–44.
Kyhn LA, Tougaard J, Jensen FH, Wahlberg M, et al. Feeding at a high pitch: source parameters of narrow band, high-frequency clicks from echolocating off-shore hourglass dolphins and coastal Hector’s dolphins. J Acoust Soc Am. 2009;125(3):1783–91.
Kyhn LA, Tougaard J, Beedholm K, Jensen FH, Ashe E, Williams R, Madsen PT. Clicking in a killer whale habitat: narrow-band, high-frequency biosonar clicks of harbour porpoise (Phocoena phocoena) and Dall’s porpoise (Phocoenoides dalli). PLoS One. 2013;8(5):e63763.
Jensen FH, Rocco A, Mansur RM, Smith BD, Janik VM, Madsen PT. Clicking in shallow rivers: short-range echolocation of Irrawaddy and Ganges River dolphins in a shallow, acoustically complex habitat. PLoS One. 2013;8(4):e59284.
Turl CW, Penner RH. Differences in echolocation click patterns of the beluga (Delphinapterus leucas) and the bottlenose dolphin (Tursiops truncatus). J Acoust Soc Am. 1989;86(2):497–502.
Madsen PT, Payne R, Kristiansen NU, Wahlberg M. Sperm whale sound production studied with ultrasound time/depth-recording tags. J Exp Biol. 2002;205:1899–906.
Johnson M, Hickmott LS, Aguilar Soto N, Madsen PT. Echolocation behaviour adapted to prey in foraging Blainville’s beaked whale (Mesoplodon densirostris). Proc R Soc B Biol Sci. 2008;275(1631):133–9.
Johnston DW, McDonald M, Polovina J, Domokos R, et al. Temporal patterns in the acoustic signals of beaked whales at Cross Seamount. Biol Lett. 2008;4(2):208–11.
Baumann-Pickering S. Baird’s beaked whale echolocation signals. J Acoust Soc Am. 2013;133(6):4321–31.
Baumann-Pickering S, Simonis AE, Wiggins SM. Aleutian Islands beaked whale echolocation signals. Mar Mamm Sci. 2013;29(1):221–7.
Moors-Murphy HB. Patterning in northern bottlenose whale (Hyperoodon ampullatus) click trains. Can Acoust. 2015;43(3):1–2.
Clark CW. Acoustic behavior of Mysticete whales. In: Sensory abilities of Cetaceans. US: Springer; 1990. p. 571–83.
Sjare BL, Smith TG. The vocal repertoire of white whales, Delphinapterus leucas, summering in the Cunningham Inlet, Northwest Territories. Can J Zool. 1986;64:407–15.
Turl CW. Echolocation abilities of the beluga, a review and comparison with the bottlenose dolphin Tursiops truncatus. In: Smith TG, St. Aubin DJ, Geraci JR, editors. Advances in research of the Beluga whale, Delphinapterus leucas, vol. 224; 1990. p. 119–28.
Racicot RA, Darroch SAF, Kohno N. Neuroanatomy and inner ear labyrinths of the narwhal, Monodon monoceros, and beluga, Delphinapterus leucas (Cetacea: Monodontidae). J Anat. 2018;233:421–39.
Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001;29(4):1165–88.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Fordyce RE. Whale evolution and Oligocene southern ocean environments. Palaeogeogr Palaeoclimatol Palaeoecol. 1980;31(C):319–36.
Slater GJ, Price SA, Santini F, Alfaro ME. Diversity versus disparity and the radiation of modern cetaceans. In: Proceedings of the Royal Society B: Biological Sciences; 2010. p. 3097–104.
Howell AB. Asymmetry in the skulls of mammals. Proc US Natl Museum. 1925;67(2599):1–18.
Willmore KE, Klingenberg CP, Hallgrímsson B. The relationship between fluctuating asymmetry and environmental variance in rhesus macaque skulls. Evolution. Soc Stud Evol. 2005;59(4):898–909.
Goswami A, Binder WJ, Meachen J, O’Keefe FR. The fossil record of phenotypic integration and modularity: a deep-time perspective on developmental and evolutionary dynamics. Proc Natl Acad Sci U S A. 2015;112(16):4891–6.
Fitzgerald EMG. A bizarre new toothed mysticete (Cetacea) from Australia and the early evolution of baleen whales. Proc R Soc B: Biol Sci R Soc. 2006;273(1604):2955–63.
Milinkovitch MC. Molecular phylogeny of cetaceans prompts revision of morphological transformations. Trends Ecol Evol. 1995;10:328–34.
Galatius A, Olsen MT, Steeman ME, Racicot RA, et al. Raising your voice: evolution of narrow-band high-frequency signals in toothed whales (Odontoceti). Biol J Linn Soc. 2019;126(2):213–24.
Benoit J, Adnet S, Welcomme J-L, Fabre P-H. New skull of Schizodelphis sulcatus Gervais, 1861 (Mammalia, Odontoceti, Eurhinodelphinidae) from the Lower Miocene of Pignan (Hérault, France) and its implications for systematics of Eurhinodelphinidae. Geobios. 2011;44(4):323–34.
Bianucci G , de Muizon C, Urbina M and Lambert O. Extensive diversity and disparity of the Early Miocene Platanistoids (Cetacea, Odontoceti) in the Southeastern Pacific (Chilcatay Formation, Peru). Life. 2020;10(3):27.
Lambert O, Bianucci G, Urbina M. Huaridelphis raimondii, a new early Miocene Squalodelphinidae (Cetacea, Odontoceti) from the Chilcatay Formation, Peru. J Vertebr Paleontol. 2014;34:987–1004.
Smith BD. In: Perrin WF, Wursig B, Thewissen JGM, editors. In Encyclopaedia of marine mammals. San Diego: Academic Press; 2002. p. 1208–13.
Park T, Mennecart B, Costeur L, Grohé C, Cooper N. Convergent evolution in toothed whale cochleae. BMC Evol Biol. 2019;19(1):1–11.
Sugimatsu H, Kojima J, Tamaki U, Bahl R, et al. Advanced technique for automatic detection and discrimination of a click train with short interclick intervals from the clicks of Ganges river dolphins (Platanista gangetica gangetica) recorded by a passive acoustic monitoring system using hydrophone arrays. Mar Technol Soc J. 2014;48(3):167–81.
Ladegaard M, Jensen FH, de Freitas M, et al. Amazon river dolphins (Inia geoffrensis) use a high-frequency short-range biosonar. J Exp Biol. 2015;218(19):3091–101.
Morisaka T, Connor RC. Predation by killer whales (Orcinus orca) and the evolution of whistle loss and narrow-band high frequency clicks in odontocetes. J Evol Biol. 2007;20(4):1439–58.
Aguilar Soto N, Johnson MP, Madsen PT, Díaz F, Domínguez I, et al. Cheetahs of the deep sea: deep foraging sprints in short-finned pilot whales off Tenerife (Canary Islands). J Anim Ecol. 2008;77:936–47.
Au WWL, Pawloski JL, Nachtigall PE. Echolocation signals and transmission beam pattern of a false killer whale (Pseudorca crassidens). J Acoust Soc Am. 1995;98(1):51–9.
Baker AN. Status, relationships, and distribution of Mesoplodon bowdoini Andrews, 1908 (Cetacea: Ziphiidae). Mar Mamm Sci. 2001;17(3):473–93.
Rommel S, Costidis AM, Fernandez A, Jepson PD, et al. Elements of beaked whale anatomy and diving physiology and some hypothetical causes of sonar-related stranding. J Cetacean Res Manag. 2005;7(3):189.
Rutherforde-Thorpe M. Notes on the osteology of a beaked whale. J Mammol JSTOR. 1938;19(3):354–62.
Yamada TK, Kitamura S, Abe S, et al. Description of a new species of beaked whale (Berardius) found in the North Pacific. Sci Rep. 2019;9:12723.
Kasuya T. Giant beaked whales. In: Perrin WF, Wursig B, Thewissen JGM, editors. Encyclopedia of Marine Mammals. London: Academic Press; 2002. p. 519–22.
Bardua C, Felice RN, Watanabe A, Fabre A-C, Goswami A. A practical guide to sliding and surface semi landmarks in morphometric analyses. Integr Organismal Biol. 2019;1(1):obz016.
Soldevilla MS, Henderson EE, Campbell GS, Wiggins SM, Hildebrand JA. Classification of Risso’s and Pacific white-sided dolphins using spectral properties of echolocation clicks. J Acoust Soc Am. 2008;124(1):609–24.
Kitchener DJ, Ross GJB, Caputi N. Variation in skull and external morphology in the false killer whale, Pseudorca crassidens, from Australia, Scotland and South Africa. Mammalia. 1990;54(1):119–36.
Nweeia MT, Nutarak C, Eichmiller FC, Eidelman N, et al. Considerations of anatomy, morphology, evolution, and function for narwhal dentition. In: Smithsonian at the poles: contributions to International Polar Year science; 2009. p. 223–40.
Ness AR. A measure of asymmetry of the skulls of odontocete whales. J Zool. 1967;153:209–21.
Cranford TW, Amundin M. Biosonar pulse production in odontocetes: the state of our knowledge. In: Thomas J, Moss CF, Vater M, editors. Echolocation in Bats and Dolphins. Chicago: University of Chigaco Press; 2004. p. 27–35.
Jakobsen L, Ratcliffe J, Surlykke A. Convergent acoustic field of view in echolocating bats. Nature. 2013;493:93–6.
Gunz P, Mitteroecker P, Neubauer S, Weber GW, et al. Principles for the virtual reconstruction of hominin crania. J Hum Evol. 2009;57(1):48–62.
Cardini A, Felizola Diniz Filho JA, Polly DP, Elton S. Biogeographic analysis using geometric morphometrics: clines in skull size and shape in a widespread African arboreal monkey. In: Lecture Notes in Earth Sciences; 2010. p. 191–217.
Gunz P, Mitteroecker P. Semilandmarks: a method for quantifying curves and surfaces. Hystrix. 2013;24:103–9.
Adams DC, Collyer ML, Kaliontzopoulou A. Geomorph: software for geometric morphometric analyses. In: R package version 3.1.0. See https://cran.r-project.org/package=geomorph; 2019.
Bookstein FL. Morphometric tools for landmark data: geometry and biology. Cambridge: Cambridge University Press; 1991.
Ellis R, Mead J. Beaked whales. A complete guide to their biology and conservation. Baltimore: John Hopkins University Press; 2017.
Reid JB, Evans PGH, Northridge SP. Atlas of cetacean distribution in northwest European waters. Peterborough, in press: Joint Nature Conservation Committee; 2003.
Jefferson TA, Wang JY. Revision of the taxonomy of finless porpoises (genus Neophocaena): the existence of two species. J Marine Anim Their Ecol. 2011;4(1):3–16.
Jefferson TA, Rosenbaum HC. Taxonomic revision of the humpback dolphins (Sousa spp.), and description of a new species from Australia. Marine Mamm Sci. 2014;30(4):1494–541.
Beasley I, Robertson KM, Arnold P. Description of a new dolphin, the Australian Snubfin dolphin Orcaella Heinsohni Sp. N. (Cetacea, Delphinidae). Marine Mamm Sci. 2005;21(3):365–400.
Vilstrup JT, Ho SY, Foote AD, et al. Mitogenomic phylogenetic analyses of the Delphinidae with an emphasis on the Globicephalinae. BMC Evol Biol. 2011;11:65.
Dalebout ML, Bake SC, Steel D, Thompson K, et al. Resurrection of Mesoplodon hotaula Deraniyagala 1963: a new species of beaked whale in the tropical Indo-Pacific. Marine Mamm Sci. 2014;30(3):1081–108.
Martin JA. From finbacks to humpbacks: investigation of the evolutionary history of the Balaenopteridae. Thesis: San Diego State University; 2014.
Demeré TA. The fossil whale, Balaenoptera Davidsonii (Cope 1872), with a review of other neogene species of Balaenoptera (Cetacea: Mysticeti). Marine Mamm Sci. 1986;2(4):277–98.
Ekdale EG. Morphological variation among the inner ears of extinct and extant baleen whales (Cetacea: Mysticeti). J Morphol. 2016;277(12):1599–615.
Lambert O, Bianucci G, Post K, et al. The giant bite of a new raptorial sperm whale from the Miocene epoch of Peru. Nature. 2010;466:105–8.
Bianucci G. The Odontoceti (Mammalia, cetacea) from Italian Pliocene. Systematics and phylogenesis of Delphinidae. Palaeontograpbia ltal. 1996;83:73–167.
Bianucci G, Sarti G, Catanzariti R, Santini U. Middle pliocene cetaceans from Monte Voltraio (Tuscany, Italy). Biostratigraphical, paleoecological and paleoclimatic observations’, Rivista Italiana di Paleontologia e Stratigrafia 1998; 104(1): 123–130.
Olson PA. Pilot whale Globicephala melas and G. muerorhynchus in Encyclopedia of Marine Mammals. Perrin W F, Wursig B, and Thewissen JGM. (eds.). Academic Press; 2nd edition. ISBN 0–12–551340-2. 2008; 847–52.
Post K, Bosselaers M. Late Pliocene occurrence of Hemisyntrachelus (Odontoceti, Delphinidae) in the southern North Sea – DEINSEA, vol. 11; 2005. p. 29–45. [ISSN 0923-9308].
Lambert O, Bianucci G, Urbina M, Geisler JH. A new inioid (Cetacea, Odontoceti, Delphinida) from the Miocene of Peru and the origin of modern dolphin and porpoise families. Zool J Linnean Soc. 2017;179:919–46.
Peredo CM, Uhen MH, Nelson MD. A new kentriodontid (Cetacea: Odontoceti) from the early Miocene Astoria Formation and a revision of the stem delphinidan family Kentriodontidae. J Vertebr Paleontol. 2018;38:e1411357.
Lambert O, de Muizon C, Bianucci G. A new archaic homodont toothed whale (Mammalia, Cetacea, Odontoceti) from the early Miocene of Peru. Geodiversitas. 2015;37:79–108.
Lucas T, Goswami A. paleomorph: geometric morphometric tools for paleobiology. In: R package version 0.1.4 (https://CRAN.R-project.org/package=paleomorph); 2017.
Guillerme T, Weisbecker V. landvR: tools for measuring landmark position variation: Zenodo; 2019. https://doi.org/10.5281/zenodo.2620785.
Kassambara A, Mundt F. factoextra: extract and visualize the results of multivariate data analyses. In: R package version 220.127.116.119; 2017. http://www.sthda.com/english/rpkgs/factoextra.
Paradis E, Schliep K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. 2018;35(3):526–8.
Harmon LJ, Weir JT, Brock CD, Glor RE, et al. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24(1):129–31.
Eastman JM, Alfaro ME, Joyce P, Hipp AL, Harmon LJ. A novel comparative method for identifying shifts in the rate of character evolution on trees. Evolution. 2011;65:3578–89.
Clavel J, Escarguel G, Merceron G. mvMORPH: an R package for fitting multivariate evolutionary models to morphometric data. Methods Ecol Evol. 2015;6(11):1311–9 Edited by Poisot T.
Yang Z. Computational molecular evolution. Oxford: Oxford University Press; 2006.
Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3(2):217–23.
Core TR. R: a language and environment for statistical computing. Vienna: R Foundation for statistical computing; 2017. URL. https://www.R-project.org/.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: linear and nonlinear mixed effects models. In: R package version 3.1–147, https://CRAN.R-project.org/package=nlme; 2020.
Coombs EJ, Clavel J, Park T, Churchill M and Goswami A. Wonky whales: the evolution of cranial asymmetry in cetaceans. Github. https://github.com/EllenJCoombs/Asymmetry-evolution-cetaceans. Accessed 11 June 2020. Release for publication. Zenodo: https://doi.org/10.5281/zenodo.3893943.
We thank Thomas Guillerme for invaluable help and guidance with his R package: landvR. We thank Ryan Felice and Carla Bardua for their help, insight, and skill with code and statistical analyses. Thank you to Graeme Lloyd and Graham Slater for giving us early access to their cetacean phylogenies. Thank you to Olivier Lambert for his knowledge and guidance on the manuscript. Thank you to the museum curators at the collections EJC and MC visited (see Additional file 1) to collect scans of skulls; we are very grateful for your wealth of knowledge and to be given access to such wonderful specimens. Finally, thank you to Eve Noirault for the scans of Pakicetus, Dorudon, and all of the terrestrial artiodactyls, and to João Leite for the scan of Notocetus. EJC was supported by the London Natural Environment Research Council Doctoral Training Partnership (London NERC DTP) training grant NE/L002485/1.
EJC was funded by the London Natural Environment Research Council Doctoral Training Partnership (London NERC DTP) training grant NE/L002485/1. JC was funded by a Marie Skłodowska-Curie Individual Fellowship (IF 797373-EVOTOOLS). TP was funded by a Maria Skłodowska-Curie individual Fellowship (748167/ECHO) and ERC Starting Grant (677774/TEMPO). MC was funded by the University of Wisconsin-Oshkosh Faculty Development Grant. AG was supported by a European Research Council Grant STG-2014–637171.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
All specimens ranked by sum radius (Σρspec). Table S2. Percentage of asymmetry in the rostrum – archaeocetes. Table S3. Percentage of asymmetry in the rostrum – mysticetes. Table S4. Percentage of asymmetry in the rostrum – odontocetes. Figure S1. Asymmetry in the cetacean skull with the rostrum removed. Figure S2. Reconstructed probability of shifts in cetacean cranial asymmetry (∑pspec) with the rostrum removed. Figure S3. Reconstructed jumps in the rate of cetacean cranial asymmetry (∑pspec) with the rostrum removed. Table S5. Akaike information criterion (AIC) rankings for each evolutionary model for asymmetry in the cetacean cranium. Figure S4. Additional morphospace occupation of cetacean crania used in this study. Figure S5. Principal Components plot with PC1 and PC2 for each specimen in the study. Model diagnostics. Figure S6. Trace of the chain for model 1. Figure S7. Trace of the chain for model 2. Figure S8. Further model diagnostics for chain 1. Figure S9. Further model diagnostics for chain 2. Table S6. Effective size (ES) for estimating the mean for each of the chains 1 and 2. Figure S10. Gelman diagnostics for the two chains. Figure S11. Asymmetry in the cetacean skull shown using a phylogeny that includes only taxa that appear in a character matrix. Figure S12. Reconstructed probability of shifts in cetacean cranial asymmetry (∑pspec) using a phylogeny that includes only taxa that appear in a character matrix. Figure S13. Reconstructed jumps in the rate of cetacean cranial asymmetry (∑pspec) using a phylogeny that includes only taxa that appear in a character matrix. Table S7. ANOVA results for each potential scenario for asymmetry in the cetacean cranium. Table S8. Frequency categories used to group all extant cetaceans for the ‘frequency echolocation’ model. Figure S14. Ziphiid skulls showing the marked asymmetry in the premaxillary crests. Table S9. List of specimens used in the study. Table S10. Skulls scanned but excluded from analysis. Figure S15. The landmark configuration with manually placed landmarks on half of the skull to be mirrored to the other half of the skull. Table S11. 123 landmarks added to the entire surface of the skull. Figure S16. The position of landmark 15 (to be mirrored as landmark 79). Table S12. Likelihood model results (AIC) for each potential scenario for asymmetry in the cetacean cranium.
About this article
Cite this article
Coombs, E.J., Clavel, J., Park, T. et al. Wonky whales: the evolution of cranial asymmetry in cetaceans. BMC Biol 18, 86 (2020). https://doi.org/10.1186/s12915-020-00805-4
- Trait evolution