Musculoskeletal networks reveal topological disparity in mammalian neck evolution
BMC Evolutionary Biology volume 17, Article number: 251 (2017)
The increase in locomotor and metabolic performance during mammalian evolution was accompanied by the limitation of the number of cervical vertebrae to only seven. In turn, nuchal muscles underwent a reorganization while forelimb muscles expanded into the neck region. As variation in the cervical spine is low, the variation in the arrangement of the neck muscles and their attachment sites (i.e., the variability of the neck’s musculoskeletal organization) is thus proposed to be an important source of neck disparity across mammals. Anatomical network analysis provides a novel framework to study the organization of the anatomical arrangement, or connectivity pattern, of the bones and muscles that constitute the mammalian neck in an evolutionary context.
Neck organization in mammals is characterized by a combination of conserved and highly variable network properties. We uncovered a conserved regionalization of the musculoskeletal organization of the neck into upper, mid and lower cervical modules. In contrast, there is a varying degree of complexity or specialization and of the integration of the pectoral elements. The musculoskeletal organization of the monotreme neck is distinctively different from that of therian mammals.
Our findings reveal that the limited number of vertebrae in the mammalian neck does not result in a low musculoskeletal disparity when examined in an evolutionary context. However, this disparity evolved late in mammalian history in parallel with the radiation of certain lineages (e.g., cetartiodactyls, xenarthrans). Disparity is further facilitated by the enhanced incorporation of forelimb muscles into the neck and their variability in attachment sites.
The increase in locomotor and metabolic performance was one of the most important innovations in the evolution of mammals [1,2,3,4,5,6]. This innovation, however, was accompanied by an exceptionally low variability in the number of presacral vertebrae compared to other tetrapods (e.g., [7,8,9,10,11,12,13,14]). In fact, the number of cervical vertebrae in mammals is limited to seven, except in extant manatees and sloths . As mammals evolved a new locomotor mode based on an increase in sagittal axial motions, their back and nuchal muscles underwent an anatomical reorganization [15, 16]. The epaxonic muscles (particularly the iliocostalis system) were reduced along with the decrease of lateral axial motion [15, 16]. With the predominance of girdle-limb system as the main propeller in mammals, pectoral muscles also expanded into the dorsal region  and were integrated into the head/neck functional unit. Studies on the evolution of the mammalian neck usually have focused on the role of those muscles emigrating from the cervical region during early development [17,18,19,20]. In contrast, the muscles that expanded into the neck have been solely investigated for their impact on shoulder and forelimb mechanics (e.g., [21,22,23,24]).
Differences in ecology and size resulted in interspecific differences in the posture and mobility of the head in mammals during standing, locomotion, foraging, oral grooming, and other daily activities (e.g., [25,26,27,28,29,30]). The morphological basis of these differences, however, is poorly understood. Variation of the cervical column length as a whole has recently been shown to be an important factor in generating morphological disparity of the neck in mammals . As a consequence of the limited variability in the number of vertebrae [9, 12] and in vertebral shape [17, 32,33,34,35], the disparity of the cervical skeleton alone is still low. Hence, we suggest that interspecific variation in the arrangement of the neck muscles plus their attachment sites on the cervical vertebrae, the skull and other bones (i.e., the variability of the musculoskeletal organization) should be an important source of morphological disparity of the neck across mammals. Although there are numerous descriptions of the myology of the neck region for almost every mammalian family, only a few studies compared the neck muscle arrangement interspecifically in an evolutionary context (e.g., [36,37,38,39]). Moreover, these studies compared neck muscles only qualitatively, which prevents the quantification of the differentiation of the neck muscles arrangement. As a consequence, it is currently unknown whether the interspecific variation in muscle attachments actually affected the changes of the musculoskeletal organization of the neck across mammals.
Anatomical Network Analysis (AnNA) provides a novel framework to study the organization of the anatomical arrangement of bones and muscles of anatomical structures (i.e., the connectivity pattern) [40, 41]. Within this framework, bones and muscles are formalized as the nodes of the network, and the physical contacts among them are formalized as the links that connect the network’s nodes. Anatomical network models thus offer a mathematical description of the organization of the body . Through such mathematical formalism we can identify and quantify structural patterns, such as anatomical modules, without a priori assumptions about the developmental or functional factors causing them  (for a recent review on morphological modularity, see ). This allows direct phylogenetic comparisons. In this context, AnNA has been used, for example, to infer evolutionary trends in the skull of tetrapods  and the phylogenetic relation between morphological complexity and modularity in the skull of primates . AnNA formalization also allows combining information on skeletal and muscular tissues as it was done to study congenital musculoskeletal malformations , secondary injuries  in humans, and hindlimb functional integration in frogs .
We can also use network parameters as proxies to infer the morphological organization of the body. Table 1 summarizes the network parameters used in this study and their most common morphological interpretation. Further details on the interpretation of network concepts in a morphological context and their historical roots have been given elsewhere (see, e.g., [41, 43, 48,49,50]). In short, every interpretation derives from the biological role of connections in the network model. Broadly speaking, the connections we have modeled among the bones and muscles of the neck embody functional interactions (e.g., determining the range of neck motion), as well as developmental factors (e.g., those inducing muscles to attach to a specific vertebrae and not to other). Thus, for example, the number of such connections for a given element, or for the entire network (i.e., K), represents the amount of functional and developmental dependences of this element or of the whole network. Functional and developmental dependences are often associated to Rupert Rield’s concept of burden, or more generally, to the concept of constrain of body parts [40, 51,52,53]. Morphological interpretation of more elaborated network parameters, such as the density of connections (D), the mean clustering coefficient (C), and the mean shortest path length (L), follow a similar logic. Because connections represent biological interactions among anatomical parts, their relative amount (D) serves as a proxy of the complexity available to the system (e.g., to perform complex functions). In addition to quantifying the amount of connections, the way connections are set (e.g., creating intertwined patterns such as 3-node loops (C)) and their effects on topology (e.g., increasing the effective proximity of elements to interact together (L)) also have consequences in the overall integration of the anatomy. Thus, the greater the intertwining, the greater the integration; the closer the elements, the greater the integration. Moreover, differences in the amount of connections among the elements of the network (some have many, most have a few) introduces heterogeneity in the organization of the network. Such heterogeneity can be related before to structural disparity (or anisomerism sensu Gregory [43, 54]). Finally, the overall patterns of integration and heterogeneity among the parts of a network often results in the emergence of new properties, for example, modularity  (see [55,56,57,58] for general reviews on the origin and macroevolution of modularity at a morphological level). The more specific details of the modular organization of a network need a closer observation, but the overall degree of parcellation of the network into large, uniform modules (P) captures how much modular the neck is.
Here we applied AnNA to a phylogenetically broad dataset of mammalian necks, including all bones (cervical vertebrae, cranium, sternum, hyoid and pectoral girdle) and muscles involved in the motion of the head and neck (Fig. 1a-f). First, we modeled the neck of each species as a network in which nodes represented the aforementioned bones and muscles and links represented their physical arrangement or contacts. Then, we quantified seven network parameters that serve as proxies for the morphological organization of the neck anatomy (Table 1). Finally, we explored the disparity of neck musculoskeletal organization through time to answer, specifically: 1) whether the musculoskeletal organization (as captured by network parameters) of the neck really differs among mammals; 2) whether closely related species share similar network organization; 3) whether there is a consistent pattern of modularity across mammalian necks; 4) how extreme elongation and deviating vertebral numbers alter neck organization; and 5) how disparity in neck anatomical organization changed during mammalian evolutionary history.
Network parameters and phylogenetic signal
The values of the network parameters used as proxies for the musculoskeletal organization of the neck for individual taxa are listed in Table 2. The phylogenetic signal is statistically significant for the multivariate dataset of all network parameters (Kmult = 0.89, p < 0.001) as well as for five of the seven parameters (N, K, D, H, P; Table 3, Additional file 1: AF1). This suggests similar variation in neck organization in closely related species (Fig. 2). A Brownian motion model of evolution best explains trait evolution of the network parameters (see Additional file 1: AF1). Relative variability in connectivity K, complexity D, and integration by co-dependency C is high among the species examined (significant higher coefficients of variation; Msrl = 139.58, p < 0.001) (Table 3, Additional file 1: AF1). In contrast, relative variability of integration by effective proximity L, anisomerism H, and the degree of modularity P are low across all mammals (Table 3, Additional file: AF1). There is no significant relationship between the network parameters and either body mass (F = 0.170, p = 0.68), absolute neck length (F = 0.005, p = 0.94), relative neck length (F = 0.661, p = 0.42), or predatory behavior (i.e., predatory vs. non-predatory; F = 0.191, p = 0.68).
The number of anatomical elements N is low in monotremes compared to the general pattern of therians (Fig. 2a). However, most xenarthrans, the chiropterans, and the Pygmy sperm whale (Kogia breviceps) also have a decreased number of elements in their neck network. The number of anatomical connections K is uniformly high in marsupials in contrast to most other mammals (Fig. 2b). K is decreased in monotremes, xenarthrans, chiropterans, the Pygmy sperm whale and the Bactrian camel (Camelus bactrianus). Morphological complexity D is high in monotremes, intermediate in marsupials and xenarthrans and tends to decrease in most of the other placental lineages (Fig. 2c). For H, the largest contrast can be found between monotremes (very low H) and therians in general, whereas the pattern within therians is not uniform (Fig. 2d). The degree of modularity P is relatively invariable and only slight decreases can be found in diprotodonts, Pen-tailed treeshrew (Ptilocercus lowii), chiropterans and eulipotyphlans (Fig. 2e). The phylomorphospace (Fig. 3; Additional file 1: AF1) highlights the differences in neck musculoskeletal organization between monotreme and therian mammals (F = 1.695, p < 0.001). Marsupials cluster closely together. They slightly overlap with the distribution of placental mammal, which occupy a huge part of the morphospace. The Pygmy sperm whale is far from the other placental mammals.
Community structure and phenotypic modules in neck networks
The number and constitution of the connectivity modules in the neck varies among the mammalian species examined (a summary is given in Table 4; for detailed information see Additional file 1: AF1). In many cases however, module number varies based on a left-right split of modules that are united in other species (e.g., the pectoral elements are separated in left and right modules). Overall, five principal connectivity modules were detected: 1) cranio-pectoral, 2) ventral, 3) mid-cervical, 4) lower cervical, and 5) thoracic. This pattern is exemplified here for the lesser grison (Galictis cuja) (Fig. 1g). The cranio-pectoral module (present in 32 out of 48 taxa) groups the cranium, the C1, and the bones of the pectoral girdle (scapulae, claviculae, and, if included, humeri), as well as the suboccipital, cleidocephalic (or cephalohumeralis), atlantoscapularis, capital longus, and capital rhomboid muscles. The ventral module groups the sternum, hyoid, thyroid, mandible (when included), and the sternocephalic and infrahyoid muscles. In 32 out of 48 cases, ribs and the related scalenii muscles are also included in this module. The ventral module is combined with (parts of) the pectoral bones and muscles in 15 species, none of which are aclaviculate (Fig. 4b). The mid-cervical module groups C2 to C4 as well as the longus cervicis, spinalis, and their related interspinal, intertransversarii, and multifidii muscles. The lower cervical module groups C5 to C7 with the cervical longissimus, spinalis, and their related interspinal, intertransversarii, and multifidi muscles. These mid-cervical and lower cervical modules are present in all networks (Fig. 4a-d). The border between them is, however, shifted in some species (e.g., C4/C5 to C3/C4). If the attachments of scalenii muscles are limited to few specific cervical vertebrae, these muscle and the ribs are also included in the mid or lower cervical module. The thoracic module groups the thoracic spine, nuchal ligament, semispinalis (complexus + biventer cervicis), capital longissimus, cervical rhomboid and trapezius muscle.
In several species, the pectoral bones (plus the related muscles) are not grouped together with the cranium and C1 but separated; otherwise they are included in the ventral or thoracic module, respectively. For instance, in the long-necked camel the pectoral bones and muscles constitute a distinct module together with the nuchal ligament (Fig. 4c). In the giraffe (Giraffa camelopardalis), attachment sites of the pectoral muscles are shifted proximally to the midcervical module. In contrast, the pectoral bones and muscles are combined with the ventral elements into a ventro-pectoral unit in the parti-colored bat (Vespertilio murinus) (Fig. 4b). The sloths differ in the organization of their neck due to their aberrant number of cervical vertebrae. In the two-toed sloth (Choloepus didactylus; six cervical vertebrae), the pectoral elements form a separate module. The C5, C6, thoracic spine, and related muscles are grouped within one module. In the three-toed sloth (Bradypus tridactylus; nine cervical vertebrae; Fig. 4d), there is an upper (C2-C5) and lower (C6-C7) midcervical module. The evolutionary ‘new’ vertebrae C8 and C9, however, are grouped with the thoracic spine, scapulae, and related muscles. The claviculae are included in the ventral module. The ribs are grouped with the cranium and atlas in the Pygmy sperm whale.
Disparity in neck organization through time
The mean subclade disparity values for the observed and simulated data were plotted against node age (Fig. 2f). Subclade disparity through time is low, which is particularly obvious in the first two-thirds of mammalian evolution (i.e., during the Mesozoic). However, it is not significantly different from the expectation under a Brownian motion model of neck organizational evolution (morphological disparity index = 0.016, p = 0.65). Nevertheless, major shifts in disparity rate occurred in the middle to late Paleocene and in the middle to late Eocene. These shifts resulted in disparity peaks exceeding the 95% confidence interval of the simulated data. The disparity in neck musculoskeletal organization decreased after the Eocene-Oligocene border. A larger sample size would be required to infer significant results for the post-Eocene ages (i.e., more divergence events are needed).
Variation in neck organization across phylogeny
The more conserved network parameters (L, H, P) represent measurements of the neck’s integration by effective proximity, anisomerism, and degree of modularity [41, 43]. These measurements capture how distant parts of the neck (e.g., head and trunk, lower and upper cervical column) are integrated, that anatomical connections are not evenly distributed across bone (e.g. vertebrae) and the way the neck is modularized. Thus, our results indicate a basic constructional set-up of the neck across mammals determined by morphological regionalization (see further below). These conserved features likely arise due to shared developmental [8, 9, 17, 19, 59,60,61] and/or biomechanical/constructional constraints [31, 62, 63]. However, the musculoskeletal organization of the neck is not uniform for other morphological features captured by network proxies. Specifically, mammalian necks considerably vary in terms of their morphological burden, complexity, and integration by co-dependency (quantified by K, D, and C, respectively). These features describe the grade of specialization in the neck due to reduction of elements or enhancement of passive structures (e.g., the nuchal ligament) and the way the neck is structurally constrained by the setup of its muscular connections (i.e., its evolvability) [43, 49]. The variation arises from the major trends of epaxonic muscle modification during mammalian evolution, leading to differences in nuchal muscle organization among monotremes, marsupials and therians. Our findings confirm the plesiomorphic pattern of epaxonic neck muscle arrangement in monotremes [64,65,66,67]. It results in a musculoskeletal organization that is distinctively different from that of therian mammals (Fig. 3). Their low differentiation within the three longitudinal systems (longissimus, iliocostalis, transversospinalis; in particular the deep intervertebral muscles; low number of muscles) and muscle attachments that are evenly distributed among the vertebrae (high complexity but low irregularity) suggest low specialization to specific neck motion patterns. In marsupials, epaxonic muscles are more differentiated in deep and superficial layers. Moreover, most of the superficial muscles are attached to every cervical vertebra [66, 68,69,70,71,72]. This high connectivity results in structural constraints in the neck (i.e., morphological burdens) and low musculoskeletal disparity among marsupials in comparison to placentals (Fig. 3). Attachments of epaxonic neck muscles are very variable among placental mammals  and thus network parameters are as well. However, two major trends in neck evolution have been shown: First, there is a reduction of attachment sites of the neck muscles to only a few vertebrae/the skull; and second, there is an increased bracing of the head-trunk distance by ligamentous structures to accommodate for increasing head weight and neck length [15, 73,74,75,76,77,78]. This results in placental mammals generally having necks that are less complex compared to monotreme and marsupial mammals. For instance, few but specialized muscles have the small attachment sites and are able to induce a similar motion or the superficial epaxonic muscles attach only secondarily to the head via the nuchal ligament. In addition, neck variation in placental mammals is also highly influenced by variation in the organization of pectoral bones and muscles (see below).
The phylogenetic signal of most of the network parameters reveals that phylogenetic relationship accounts for much of the variation in neck organization. At the same time, network parameters also discriminate between monotreme, marsupial, and placental mammals. Within placental mammals, however, variation of network parameters is mostly limited to xenarthrans, chiropterans, and some cetartiodactyls. Although relatively species-poor, xenarthrans, show highly specialized neck morphologies related to their diverse fossorial or arboreal ecologies [79,80,81] and unique development [33, 82, 83]. In chiropterans, back muscles contribute only marginally to the stabilization of the head because of the lack of most cranial and cervical attachments [84,85,86]. Sagittal stability is instead achieved by the modified morphology of the cervical vertebrae in accordance with roosting behaviors . Cetartiodactyls have recently been show to exhibit the highest disparity in neck morphology across mammals . It ranges from the very short necks of cetaceans up to the extreme long ones of camelids and giraffids. As a consequence, neck musculoskeletal organization is similarly diverse. Several muscles with cervical attachment are reduced in the Pygmy sperm whale (and other cetaceans) or their attachment is shifted to the skull (e.g. scaleni muscles) and thus head stabilization is increased [88, 89]. On the other hand, cranial attachments of the dorsal neck muscle are mainly reduced in long necked species, such as the camel and giraffe. Muscle force is instead transferred by the modified nuchal ligament [74, 75, 77]. Surprisingly, neck network parameters in the dugong (Dugong dugon), although also being fully aquatic, do not show a similar alteration as in the Pygmy sperm whale. Instead, it closely resembles the Asian elephant (Elephas maximus) and other afrotherians (Table 1, Fig. 2).
In accordance with our findings, [90, 91] also showed that the effect of size and prey capture behavior is low in the neck compared to the thoracolumbar region. However, functional interpretations of the results of the analysis of the topological arrangement of parts needs to be inferred on a one to one basis and taking into account the specific ecological context of each taxa.
Regionalization and modularity in the mammalian neck
Despite the relative low and invariant number of neck vertebrae in mammals, several studies have uncovered a tripartite regionalization of the cervical spine based on developmental, morphological, allometric, and functional evidence [17, 25, 31, 32, 91]. Our results have now uncovered a corresponding regionalization of the musculoskeletal organization of the mammalian neck into an upper (cranium, C1), mid (C2-C4), and lower cervical module (C5-C7, in some species also the thoracic spine). This modularity pattern is conserved across mammals despite variations in size, feeding mechanisms, and locomotor modes (indicated by a uniform grade of modularity). This conserved pattern probably arose from the high number of connections between the vertebrae of the same module (or the cranium and C1) resulting in an increase of structural constraints and integration of these elements [40, 49, 52]. However, the boundaries between adjacent modules/regions are not consistent across different studies analyzing the morphology of the neck using different criteria. For example, vertebrae C1 and C2 are not part of the same connectivity module despite their close developmental, functional, and evolutionary relationship [25, 92, 93]. A similar dissociation of C1 and C2 into different regions has been shown for their scaling properties  and highlights the role of C2 as a functional mediator between the head joint and the postaxial column (see also ). In addition to the three ‘inner’ axial modules, there are two additional ‘outer’ modules bridging the distance between the trunk and the head (or the hyoid or upper vertebrae), with a muscular cuff on the dorso-lateral (pectoral) and ventro-lateral side. Many of these muscles were crucial for the evolutionary origin of the vertebrate neck [95, 96].
Neck organization in sloths
In general, a similar regionalization of the neck is observed in both genera of sloths, despite their variation in the number of cervical vertebrae. The evolutionary new C8 and C9 in Bradypus and their associated muscles are grouped together with the thoracic spine. This agrees with their thoracic origin and ossification sequence . Conversely, the evolutionary new Th1 provides the basis for the close association of the thoracic vertebral region to C5 and C6 in Choloepus neck organization. Divergence of the sloths’ necks becomes obvious when including their pectoral bones and muscles in the comparison. Their neck-shoulder arrangement represents two different solutions of locomotor possibilities under common functional constraints [80, 97]. Neck organization and modularity of Bradypus resembles those of other long necked species. The resemblance stems from its nearly complete lack of cervical and cranial attachments of the neck/shoulder muscles and its unusual clavicular and shoulder morphology [80, 98,99,100] (see the unusual lower cervical-thoracic-scapular module in Table 4). The muscles of Choloepus, in contrast, are so placed as to offer the greatest possible support dorsally and ventrally to the head as well as to the scapula [80, 98, 101, 102]. Thus, neck organization and modularity is closer to the general pattern as seen in the lesser grison (i.e., pronounced head support, functional connection of head and forelimb) .
Evolutionary integration of the neck and the forelimb
The enduring evolutionary and developmental relationship between the neck and the forelimb in mammals (and other amniotes) is well documented (e.g., [17, 19, 20, 95, 96, 103]). This relationship is most obvious in the brachial plexus innervating shoulder and forelimb muscles [104, 105]. However, there is also a strong functional integration between the neck and the forelimb, with several muscles connecting the pectoral girdle and the head/neck (often with repeated slips). Based on this functional connection, the posture and movements of the neck have a crucial influence on the mechanics of the forelimb in terms of gait efficiency, balance, stabilization, ground reaction forces, and kinematics [30, 106,107,108,109,110,111,112]. Our findings now highlight the consequences of this integration on the musculoskeletal organization of the neck. Although there is a conserved tripartition of the cervical spine and its associated muscles, the varying contribution of pectoral bones and muscles to different connectivity modules accounts for much of the observable neck disparity across mammals (e.g., see results on sloths). Major shifts in forelimb morphology and function (e.g., mobilization of the pectoral girdle, reduction of the claviculae) [113,114,115] are associated with increasing decoupling of the pectoral elements from the ventral module and their connection to the cranium and upper cervical region. This coincides with the increased role of head/neck movements on forelimb mechanics during fast and endurance running (i.e., cursoriality) [30, 108]. In mammals with extreme long necks (camel, giraffe), although also being capable of enduring walking, the pectoral bones and muscles are separated from the cranial module. However, it has recently been shown that neck pendulum mechanics and function in long-necked mammals is different compared to actual cursors .
Implications for the evolution and disparity of the mammalian neck
The differences in neck organization between monotremes and therian mammals is one of the striking findings of this study. They result in high disparity between them whereas their within-subclade disparity is low during the first two-thirds of mammalian evolution. Accordingly, the disparity of the neck of mammals was low during the Mesozoic (Fig. 2f) . Figure 5 illustrates the major grades of musculoskeletal organization during mammalian neck evolution. Monotreme, marsupial, and placental mammals differ in their degree of epaxonic muscle differentiation and the varying integration of the forelimb muscles. In general, the morphological complexity of the neck decreases from monotremes to placentals (Fig. 5) but disparity increases (Fig. 3).
The therian radiation that followed the K/T mass extinction and the appearance of most of the therian and marsupial (supra)orders during the Paleocene was accompanied by an abrupt increase in neck musculoskeletal disparity after a long period of low neck disparity. This was even associated with the appearance of locomotor and foraging specializations  and the diversification in body size (e.g., ). The second disparity peak during the Eocene coincides with the radiation and increased diversity of modern placental orders, like cetartiodactyls, perissodactyls, carnivorans, and xenarthrans (see  and references therein). Thus, disparity in neck organization emerged relatively late in the long mammalian history and is associated with the origin and radiation of specific lineages.
One of our crucial findings is that the musculoskeletal organization of the neck differs between monotreme, marsupial, and placental mammals (Fig. 5). Moreover, particularly the necks of placental mammals are characterized by a reduced complexity despite their increased disparity in musculoskeletal organization and length. Our network analyses revealed a mosaic complexity and disparity in the musculoskeletal organization of the mammalian neck despite the more obvious meristic (and other) constraints on the cervical spine. Musculoskeletal irregularity, effective proximity, degree of modularity, and the occurrence of three inner/axial regions are conserved features among mammalian necks. Thus, a shared biomechanical construction and common developmental interrelationships not only constrain variation in the cervical spine, but are similarly likely to limit musculoskeletal variability in the neck. The conservation of these traits contrasts, however, with the high variability in morphological burden, integration by co-relation, morphological complexity, and the configuration of the ventral and (cranio-)pectoral module in the neck. The expansion of limb muscles in the cervical region not only facilitated enhanced forelimb mechanics but also increased structural disparity (and thus derived motor patterns and mechanics) in the neck. Thus, we highlight the close integration of the neck and the forelimb during mammalian evolution. The disparity in neck musculoskeletal organization evolved late in mammalian history and in parallel with the radiation of some lineages (e.g., cetartiodactyls, xenarthrans). Finally, our findings show that the limited number of vertebrae in the cervical spine does not necessarily result in low musculoskeletal disparity during mammalian evolutionary diversification.
Data collection and anatomical network modeling
We collected the topographic data of the musculoskeletal system of the neck in 48 mammalian species through an extensive literature review (Table 2, Additional file 2: AF2). The sample represents all major monotreme, marsupial, and placental clades, as well a diversity of locomotor and feeding strategies. Rodent diversity is represented by members of all suborders (sciuromorphs, myomorphs, hystricomorphs, castorimorphs, and anumaluromorphs, respectively). Representatives of both genera of extant sloth (Bradypus and Choloepus) were included to examine the influence of their deviating number of cervical vertebrae on neck organization. We documented the number and specific connections/attachments of all skeletal structures and muscles constituting the neck motion system in these taxa.
In contrast to the domestic mammals, for which veterinary textbooks as well as anatomical publications are available, information for exotic species is scarce and literature descriptions are in many cases old. To overcome this problem, in all but two cases we always consulted two or more references to compare the data. We made sure that at least one detailed anatomical monograph was included. As the network models used here focus on the presence/absence of attachments among bones and muscles rather than more specific information on the nature and area of such attachments, the literature reviewed was of enough quality for our modeling approach.
We included all muscles originating from the cervical vertebrae, skull (cranium or mandible), nuchal ligament, or hyoid/thyroid, and inserting on the (cervical or thoracic) vertebrae, sternum, pectoral girdle (scapulae, claviculae, humeri), or ribs (see details in Additional file 3: AF3). Accordingly, we excluded the masticatory, facial, laryngeal, pharyngeal, and suprahyoid muscles from the analysis. A calibrated phylogenetic tree was constructed using the data from the Timetree of Life database [119, 120] (Fig. 2).
We built anatomical network models of the necks’ musculoskeletal systems, which comprised all anatomical units as well as the different types of physical interaction among them. Network nodes represented all bones (cranium, cervical vertebrae, thoracic spine, left and right ribs, hyoid, left and right claviculae, left and right scapulae, sternum, and, if involved, left and right humeri, and mandible), other passive elements (nuchal ligament, if present, thyroid), and all cervical muscles, as described above. Network connections represented all physical articulations between bones and other passive elements described, as well as the fleshy and tendinous attachments of the muscles onto the bones (see adjacency matrices in Additional file 4: AF4). Network models were analyzed using the igraph package  in R .
Network parameter analyses
The mathematical definitions and calculations of the network parameters examined here (N, K, D, C, L, H, P) are provided in Table 1. The degree of modularity (parcellation P) was measured from the connectivity modules identified using a spin-glass model and simulated annealing algorithm implemented in the R package netcarto [123, 124]. A connectivity module is defined as a group of nodes highly connected among them and poorly connected to nodes outside the group. P is 0 when all nodes are in a same module, and tends to 1 when nodes are evenly distributed within many modules. We tested the phylogenetic signal in the multivariate dataset of all network parameters by calculating Kmult  using R package phylocurve . We additionally tested the phylogenetic signal of the individual parameters using the Abouheif’s test with 1000 permutations  and Blomberg’s test  in the R package phytools . Mode of trait evolution was explored by comparing multivariate fits of Brownian motion, Ornstein-Uhlenbeck (single adaptive optimum), and Early burst models using Akaike information criterion (AIC) weights in R package mvMORPH . Distribution of the network parameters was visualized with a phylomorphospace of the first two principal components using R package phytools .
Relative variability of the network parameters was analyzed by statistical comparison of their coefficients of variation (CVs). 95% confidence intervals of the CVs were calculated by 10,000 bootstrap resampling. Significant differences among the parameters’ CVs were tested using the modified signed-likelihood ratio test (MSLRT) for equality of CVs (all parameters) and the asymptotic test for the equality of CVs (pair-wise comparisons, Bonferroni corrected) in the R package cvequality . In order to test for allometric effects on network parameters they were regressed against logtransformed body mass, absolute neck length, and relative neck length. Body mass and neck length data were taken from . Relative neck lengths were calculated by dividing absolute neck length by body mass1/3. Allometric analyses were done using phylogenetic generalized least square regressions in the R package caper . The effect of predatory behavior on logtransformed network parameters was examined by testing for significant differences between predatory and non-predatory mammals using a multivariate distance-based phylogenetic generalized least square regression (D-PGLS) with 1000 permutations in the R package geomorph . Species were classified as predatory when food intake involves head-neck movements to hold the food counteracting its resisting movements (carnivorous, insectivorous, piscivorous species). In contrast, species were classified as non-predatory when food intake does not involve such head-neck movements (food is just picked or harvested, e.g., browsers, grazers, but also myrmecophagous species). Differences in network parameters between monotreme, marsupial, and placental mammals were tested using a multivariate D-PGLS regression with 1000 permutations.
We calculated the quality of the partitions identified by the community detection algorithm using the optimization function Q . According to Newman and Girvan , if the number of connections within modules is not different from that expected at random, then Q will be close to 0. The higher the Q the stronger the modular pattern of the network (Qmax = 1). In practice, strongly modular networks show Q values ranging from 0.3 to 0.7 . Thus, we considered that an anatomical network has a strongly modular structure if Q − Qerror > 0.3. The expected error of Q was calculated using a jackknife procedure, where every link was taken as an independent observation  (more details are provided in Additional file 3: AF3). Finally, we estimated the statistical significance of each module using a two-sample Wilcoxon rank-sum test on the internal vs. external links of the module’s nodes. The null hypothesis was that the number of connections is the same inside as outside the module (i.e., as expected if the module were created at random); the alternative hypothesis was that the number of connections is higher inside than outside the module (i.e., the definition of connectivity module). An extensive account of these methods has been given elsewhere [41, 48, 135].
Disparity through time analysis
We carried out a disparity through time (DTT) analysis using the R package geiger  to trace the variation in neck organization through the evolution of mammals. First, we performed a principal component analysis (PCA) of the network parameters used as proxies of the morphological organization of the neck (i.e., N, K, C, D, L, H, and P) to account for their co-variation structure. Mean subclade disparity through time for the PC scores were calculated [116, 137]. Observed disparity in neck organization across our phylogeny was compared with that expected under a Brownian motion process performing 10,000 iterations. High disparity values indicate high variance within subclades; low disparity values indicate conservation within subclades and high variance among subclades. Finally, we calculated the morphological disparity index to quantify the overall difference in relative disparity of a clade compared to that expected under the null Brownian motion model [116, 137].
Abbreviations of network elements
Left and right side are indicated by l and r, respectively, added to the abbreviations.
asd atlantoscapularis dorsalis; asv atlantoscapularis ventralis; bc biventer cervicis; C1-C9 cervical vertebrae; cal longus capitis; cc cleidocervicalis; cl clavicle; cm cleidomastoideus; co cleidooccipitalis; col longus colli; cx complexus; cr cranium; hu humerus; hy hyoid; icc iliocostalis cervicis; id1-id8 intertransversarii cervicis dorsales; im1-im4 intertransversarii cervicis mediales; is1-is8 interspinalis; iv/iv1-iv8 intertransversarii cervicis ventrales (fused/separate); lat longus atlantis; lca longissimus capitis; lce longissimus cervicis; ln nuchal ligament; m1-m6 multifidi; md mandible; oca obliquus capitis caudalis; ocr obliquus capitis cranialis; oh omohyoideus; rca rhomboideus capitis; rce rhomboideus cervicis; rci rectus capitis dorsalis intermedius; rh rhomboideus (undifferentiated); ri ribs; rl rectus capitis lateralis; rma rectus capitis dorsalis major; rmi rectus capitis dorsalis minor; rv rectus capitis ventralis; sce spinalis cervicis; sc scapula; scm sternocleidomastoideus (sternal and clavicle part not separate); sd scalenus dorsalis; sh sternohyoideus; sm scalenus medius; so sternooccipitalis; spca splenius capitis; spce splenius cervicis; sp splenius (undifferentiated); ssca semispinalis capitis; ssce semispinalis cervicis; st sternum; sth sternothyroideus; stm sternomastoideus; stx sternomaxillaris; svc serratus ventralis cervicis; sv scalenus ventralis; tr trapezius; ts thoracic spine; ty thyroid
Alexander R, Dimery NJ, Ker R. Elastic structures in the back and their role in galloping in some mammals. J Zool. 1985;207(4):467–82.
Bramble DM, Carrier DR. Running and breathing in mammals. Science. 1983;219(4582):251–6.
Fischer MS. Crouched posture and high fulcrum, a principle in the locomotion of small mammals: the example of the rock hyrax (Procavia capensis)(Mammalia: Hyracoidea). J Hum Evol. 1994;26(5-6):501–24.
Gambaryan PP, Hardin H. How mammals run: anatomical adaptations. New York: Wiley; 1974.
Hildebrand M. Motions of the running cheetah and horse. J Mammal. 1959;40(4):481–95.
Magne de la Croix P. The evolution of locomotion in mammals. J Mammal. 1936;17(1):51–4.
Asher R, Lin K, Kardjilov N, Hautier L. Variability and constraint in the mammalian vertebral column. J Evol Biol. 2011;24(5):1080–90.
Buchholtz EA: Flexibility and constraint: patterning the axial skeleton in mammals. In: From clone to bone: the synergy of morphological and molecular tools in Palaeobiology. Edited by Asher R, Müller J. Cambridge: Cambridge University Press; 2012: 230-256.
Galis F. Why do almost all mammals have seven cervical vertebrae? Developmental constraints, Hox genes, and cancer. J Exp Zool. 1999;285(1):19–26.
Galis F, Carrier DR, Van Alphen J, Van Der Mije SD, Van Dooren TJ, Metz JA, ten Broek CM. Fast running restricts evolutionary change of the vertebral column in mammals. Proc Natl Acad Sci. 2014;111(31):11401–6.
Müller J, Scheyer TM, Head JJ, Barrett PM, Werneburg I, Ericson PG, Pol D, Sánchez-Villagra MR. Homeotic effects, somitogenesis and the evolution of vertebral numbers in recent and fossil amniotes. Proc Natl Acad Sci. 2010;107(5):2118–23.
Narita Y, Kuratani S. Evolution of the vertebral formulae in mammals: a perspective on developmental constraints. J Exp Zool B Mol Dev Evol. 2005;304(2):91–106.
Sánchez-Villagra MR, Narita Y, Kuratani S. Thoracolumbar vertebral number: the first skeletal synapomorphy for afrotherian mammals. Syst Biodivers. 2007;5(1):1–7.
Todd TW. Numerical significance in the thoracicolumbar vertebrae of the mammalia. Anat Rec. 1922;24(5):260–86.
Jouffroy F. Evolution of the dorsal muscles of the spine in light of their adaptation to gravity effects. In: Berthoz A, Graf W, Vidal PP, editors. The head-neck sensory motor system. Oxford: Oxford University Press; 1992. p. 22–35.
Slijper EJ. Comparative biological-anatomical investigations of the vertebral column and spinal musculature. K Ned Akad Wet, Verh (Tweede Sectie). 1946;42:1–128.
Buchholtz EA, Bailin HG, Laves SA, Yang JT, Chan MY, Drozd LE. Fixed cervical count and the origin of the mammalian diaphragm. Evol Dev. 2012;14(5):399–411.
Burke AC, Nowicki J. A new view of patterning domains in the vertebrate mesoderm. Dev Cell. 2003;4(2):159–65.
Hirasawa T, Fujimoto S, Kuratani S. Expansion of the neck reconstituted the shoulder–diaphragm in amniote evolution. Develop Growth Differ. 2016;58(1):143–53.
Hirasawa T, Kuratani S. A new scenario of the evolutionary derivation of the mammalian diaphragm from shoulder muscles. J Anat. 2013;222(5):504–17.
Carrier DR, Deban SM, Fischbein T. Locomotor function of forelimb protractor and retractor muscles of dogs: evidence of strut-like behavior at the shoulder. J Exp Biol. 2008;211(1):150–62.
Deban SM, Schilling N, Carrier DR. Activity of extrinsic limb muscles in dogs at walk, trot and gallop. J Exp Biol. 2012;215(2):287–300.
English AWM. An electromyographic analysis of forelimb muscles during overground stepping in the cat. J Exp Biol. 1978;76(1):105–22.
Payne R, Veenman P, Wilson A. The role of the extrinsic thoracic limb muscles in equine locomotion. J Anat. 2004;205(6):479–90.
Graf W, De Waele C, Vidal P. Functional anatomy of the head-neck movement system of quadrupedal and bipedal mammals. J Anat. 1995;186(Pt 1):55.
Graf W, de Waele C, Vidal P, Wang D, Evinger C. The orientation of the cervical vertebral column in unrestrained awake animals. II. Movement strategies. Brain Behav Evol. 1995;45(4):209–31.
Hart BL, Hart LA, Mooring MS, Olubayo R. Biological basis of grooming behaviour in antelope: the body-size, vigilance and habitat principles. Anim Behav. 1992;44(4):615–31.
Mooring MS, Blumstein DT, Stoner CJ. The evolution of parasite-defence grooming in ungulates. Biol J Linn Soc. 2004;81(1):17–37.
Vidal PP, Graf W, Berthoz A. The orientation of the cervical vertebral column in unrestrained awake animals. Exp Brain Res. 1986;61(3):549–59.
Zsoldos RR, Licka TF. The equine neck and its function during movement and locomotion. Zoology. 2015;118(5):364–76.
Arnold P, Amson E, Fischer MS. Differential scaling patterns of vertebrae and the evolution of neck length in mammals. Evolution. 2017;71(6):1587–1599.
Arnold P, Forterre F, Lang J, Fischer MS. Morphological disparity, conservatism, and integration in the canine lower cervical spine: insights into mammalian neck function and regionalization. Mamm Biology-Zeitschrift für Säugetierkunde. 2016;81(2):153–62.
Buchholtz EA, Stepien CC. Anatomical transformation in mammals: developmental origin of aberrant cervical anatomy in tree sloths. Evol Dev. 2009;11(1):69–79.
Buchholtz EA, Wayrynen KL, Lin IW. Breaking constraint: axial patterning in Trichechus (Mammalia: Sirenia). Evol Dev. 2014;16(6):382–93.
Johnson D, McAndrew T, Oguz Ö. Shape differences in the cervical and upper thoracic vertebrae in rats (Rattus norvegicus) and bats (Pteropus poiocephalus): can we see shape patterns derived from position in column and species membership? J Anat. 1999;194(2):249–53.
Bekele A. The comparative functional morphology of some head muscles of the rodents Tachyoryctes splendens and Rattus rattus. II. Cervical muscles. Mammalia. 1983;47(4):549–72.
Diogo R. The head and neck muscles of the Philippine colugo (Dermoptera: Cynocephalus volans), with a comparison to tree-shrews, primates, and other mammals. J Morphol. 2009;270(1):14–51.
Diogo R, Wood B. Soft-tissue anatomy of the primates: phylogenetic analyses based on the muscles of the head, neck, pectoral region and upper limb, with notes on the evolution of these muscles. J Anat. 2011;219(3):273–359.
Filan SL. Myology of the head and neck of the bandicoot (Marsupialia, Peramelemorphia). Aust J Zool. 1990;38(6):617–34.
Esteve-Altava B, Marugán-Lobón J, Botella H, Bastir M, Rasskin-Gutman D. Grist for Riedl's mill: a network model perspective on the integration and modularity of the human skull. J Exp Zool B Mol Dev Evol. 2013;320(8):489–500.
Rasskin-Gutman D, Esteve-Altava B. Connecting the dots: anatomical network analysis in morphological EvoDevo. Biol Theory. 2014;9(2):178–93.
Esteve-Altava B. In search of morphological modules: a systematic review. Biol Rev. 2017;92(3):1332–47.
Esteve-Altava B, Marugán-Lobón J, Botella H, Rasskin-Gutman D. Structural constraints in the evolution of the tetrapod skull complexity: Williston’s law revisited using network models. Evol Biol. 2013;40(2):209–19.
Esteve-Altava B, Boughner JC, Diogo R, Villmoare BA, Rasskin-Gutman D. Anatomical network analysis shows decoupling of modular lability and complexity in the evolution of the primate skull. PLoS One. 2015;10(5):e0127653.
Diogo R, Esteve-Altava B, Smith C, Boughner JC, Rasskin-Gutman D. Anatomical network comparison of human upper and lower, newborn and adult, and normal and abnormal limbs, with notes on development, pathology and limb serial homology vs. homoplasy. PLoS One. 2015;10:e0140030.
Murphy AC, Muldoon SF, Baker D, Lastowka A, Bennett B, Yang M, Bassett DS. Structure, function, and control of the musculoskeletal network. In: arXiv preprint arXiv:161206336; 2016.
Dos Santos DA, Fratani J, Ponssa ML, Abdala V. Network architecture associated with the highly specialized hindlimb of frogs. PLoS One. 2017;12(5):e0177819.
Esteve-Altava B. Challenges in identifying and interpreting organizational modules in morphology. J Morphol. 2017;278(7):960–74.
Rasskin-Gutman D. Boundary constraints for the emergence of form. In: Müller GB, Newman SA, editors. Origination of Organismal form: beyond the gene in developmental and evolutionary biology. Camebridge: MIT Press; 2003. p. 305.
Rasskin-Gutman D, Buscalioni AD. Theoretical morphology of the Archosaur (Reptilia: Diapsida) pelvic girdle. Paleobiology. 2001;27(1):59–78.
Riedl R, Jefferies RPS. Order in living organisms: a systems analysis of evolution. New York: Wiley; 1978.
Schoch RR. Riedl's burden and the body plan: selection, constraint, and deep time. J Exp Zool B Mol Dev Evol. 2010;314(1):1–10.
Wagner GP, Laubichler MD. Rupert Riedl and the re-synthesis of evolutionary and developmental biology: body plans and evolvability. J Exp Zool B Mol Dev Evol. 2004;302(1):92–102.
Gregory WK. Polyisomerism and anisomerism in cranial and dental evolution among vertebrates. Proc Natl Acad Sci. 1934;20(1):1–9.
Eble GJ. Morphological modularity and macroevolution: conceptual and empirical aspects. In: Callebaut W, Rasskin-Gutman D, editors. Modularity: understanding the development and evolution of natural complex systems. Camebridge: The MIT Press; 2005. p. 221–38.
Goswami A, Smaers J, Soligo C, Polly P. The macroevolutionary consequences of phenotypic integration: from development to deep time. Philos Trans R Soc B: Biol Sci. 2014;369(1649):20130254.
Klingenberg CP. Morphological integration and developmental modularity. Annu Rev Ecol Evol Syst. 2008;39:115–32.
Wagner GP, Pavlicev M, Cheverud JM. The road to modularity. Nat Rev Genet. 2007;8(12):921–31.
Buchholtz EA. Crossing the frontier: a hypothesis for the origins of meristic constraint in mammalian axial patterning. Zoology. 2014;117(1):64–9.
Galis F, Metz JA. Anti-cancer selection as a source of developmental and evolutionary constraints. BioEssays. 2003;25(11):1035–9.
Galis F, Van Dooren TJ, Feuth JD, Metz JA, Witkam A, Ruinard S, Steigenga MJ, Wunaendts LC. Extreme selection in humans against homeotic transformations of cervical vertebrae. Evolution. 2006;60(12):2643–54.
Kummer B. Bauprinzipien des Säugerskeletes Stuttgart. Germany: Georg Thieme Verlag; 1959.
Kummer B: Biomechanik des Säugetierskeletts. In: Handbuch der Zoologie, Band 8 (6). Edited by Helmcke J-G, Legerken Hv, Starck D, vol. 24. Berlin: Walter de Gruyter & Co. Verlag; 1959.
Coues E. On the myology of the Ornithorhynchus. Proc Essex Inst. 1871;6:127–73.
Jouffroy F, Lessertisseur J, Saban R. Particularités musculaires des Monotrémes – musculature post-craniénne. In: Grassé P-P, editor. Traité de Zoologie Mamifères, Tome XVI, Fac III. Paris: Masson et Cie editeurs; 1971. p. 679–836.
Nishi S. Zur vergleichenden Anatomie der eigentlichen (genuinen) Rückenmuskeln. Gegenbaurs Morphol Jahrb. 1916;50:167–318.
Virchow H. Die tiefen Rückenmuskeln des Ornithorhynchus. Gegenbaurs Morphol Jahrb. 1929;60:481–559.
Barbour RA. The musculature and limb plexuses of Trichosurus vulpecula. Aust J Zool. 1963;11(4):488–610.
Coues E, Wyman J. On the osteology and myology of Didelphys virginiana. Mem Boston Soc Nat Hist. 1872;2:41–154.
Jouffroy FK. Musculature épisomatique. In: Grassé PP, editor. Traité de Zoologie Mamifères, Tome XVI, Fac II. Paris: Masson et Cie editeurs; 1968. p. 479–548.
Jüschke S. Untersuchungen zur funktionellen Anpassung der Rückenmuskulatur und der Wirbelsäule quadrupeder Affen und Känguruhs. Anat Embryol. 1972;137(1):47–85.
Osgood WH, Herrick CJ, Obenchain JB. A monographic study of the American marsupial Caenolestes, with a description of the brain of Caenolestes by C. Judson Herrick. Field Mus Nat Hist (Zool Ser). 1921;14:1–162.
Bianchi M. The thickness, shape and arrangement of elastic fibres within the nuchal ligament from various animal species. Anat Anz. 1988;169(1):53–66.
Dimery NJ, Alexander R, Deyst KA. Mechanics of the ligamentum nuchae of some artiodactyls. J Zool. 1985;206(3):341–51.
Endo H, Yamagiwa D, Fujisawa M, Kimura J, Kurohmaru M, Hayashi Y. Modified neck muscular system of the giraffe (Giraffa camelopardalis). Ann Anat-Anat Anz. 1997;179(5):481–5.
Gellman K, Bertram J. The equine nuchal ligament 1: structural and material properties. Vet Comp Orthop Traumatol. 2002;15(1):1–6.
Mobarak A, Fouad S. A study on Lig. Nuchae of the one-humped camel (Camelus dromedarius). Anat Histol Embryol. 1977;6(2):188–90.
Preuschoft H, Klein N. Torsion and bending in the neck and tail of sauropod dinosaurs and the function of cervical ribs: insights from functional morphology and biomechanics. PLoS One. 2013;8:e78574.
Galliari FC, Carlini AA, Sánchez-Villagra MR. Evolution of the axial skeleton in armadillos (Mammalia, Dasypodidae). Mamm Biology-Zeitschrift für Säugetierkunde. 2010;75(4):326–33.
Miller RA. Functional adaptations in the forelimb of the sloths. J Mammal. 1935;16(1):38–51.
VanBuren CS, Evans DC. Evolution and function of anterior cervical vertebral fusion in tetrapods. Biol Rev. 2017;92:608–26.
Hautier L, Weisbecker V, Sánchez-Villagra MR, Goswami A, Asher RJ. Skeletal development in sloths and the evolution of mammalian vertebral patterning. Proc Natl Acad Sci. 2010;107(44):18903–8.
Varela-Lasheras I, Bakker AJ, van der Mije SD, Metz JA, van Alphen J, Galis F. Breaking evolutionary and pleiotropic constraints in mammals: on sloths, manatees and homeotic mutations. EvoDevo. 2011;2:11.
Macalister A. The myology of the Cheiroptera. Philos Trans R Soc Lond. 1872;162:125–71.
Maisonneuve P. Traité de l'ostéologie et de la myologie du Vespertilio murinus, précédé d'un exposé de la classification des chéiroptèi et de considérations sur les moeurs de ces animaux. Paris: O. Doin; 1878.
Mori M. Muskulatur des Pteropus edulis. Okajimas Folia Anat Jpn. 1960;36(3-4):253–307.
Fenton MB, Crerar LM. Cervical vertebrae in relation to roosting posture in bats. J Mammal. 1984;65(3):395–403.
Howell A. Aquatic mammals. Springfield: Charles C Thomas; 1930.
Schulte HW, Smith MDF. The external characters, skeletal muscles, and peripheral nerves of Kogia breviceps (Blainville). Bull Am Mus Nat Hist. 1918;38:7–72.
Randau M, Cuff AR, Hutchinson JR, Pierce SE, Goswami A. Regional differentiation of felid vertebral column evolution: a study of 3D shape trajectories. Organ Divers Evol. 2017;17(1):305–19.
Randau M, Goswami A, Hutchinson JR, Cuff AR, Pierce SE. Cryptic complexity in felid vertebral evolution: shape differentiation and allometry of the axial skeleton. Zool J Linnean Soc. 2016;178:183–202.
Evans FG. The morphology and functional evolution of the atlas-axis complex from fish to mammals. Ann N Y Acad Sci. 1939;39(1):29–104.
Jenkins FA. The evolution and development of the dens of the mammalian axis. Anat Rec. 1969;164(2):173–84.
Bogduk N, Mercer S. Biomechanics of the cervical spine. I: normal kinematics. Clin Biomech. 2000;15(9):633–48.
Matsuoka T, Ahlberg PE, Kessaris N, Iannarelli P, Dennehy U, Richardson WD, McMahon AP, Koentges G. Neural crest origins of the neck and shoulder. Nature. 2005;436(7049):347–55.
Ericsson R, Knight R, Johanson Z. Evolution and development of the vertebrate neck. J Anat. 2013;222(1):67–78.
Nyakatura JA. The convergent evolution of suspensory posture and locomotion in tree sloths. J Mamm Evol. 2012;19(3):225–34.
Humphry G. The myology of the limbs of the Unau, the Ai, the two-toed anteater, and the pangolin. J Anat Physiol. 1869;3:2–78.
Macalister A. On the myology of Bradypus Tridactylus; with remarks on the general mucular anatomy of the Edentata. Ann Mag Nat Hist. 1869;4(19):51–67.
Mackintosh H. On the myology of the genus Bradypus. Proc R Ir Acad. 1870;1:517–29.
Lucae JCG. Statik und Mechanik der Quadrupeden an dem Skelett und den Muskeln des Lemur und eines Choloepus. Abhandlungen der Senckenbergischen Naturforschenden Gesellschaft Frankfurt. 1884;13:1–92.
Mackintosh H. On the muscular anatomy of Choloepus didactylus. Proc R Ir Acad. 1875;1:66–67^69.
Gross MK, Moran-Rivard L, Velasquez T, Nakatsu MN, Jagla K, Goulding M. Lbx1 is required for muscle precursor migration along a lateral pathway into the limb. Development. 2000;127(2):413–24.
Howell AB. Morphogenesis of the shoulder architecture. Part VI. Therian Mammalia. Q Rev Biol. 1937;12(4):440–63.
Miller RA. Comparative studies upon the morphology and distribution of the brachial plexus. Am J Anat. 1934;54(1):143–75.
Dunbar DC, Macpherson JM, Simmons RW, Zarcades A. Stabilization and mobility of the head, neck and trunk in horses during overground locomotion: comparisons with humans and other primates. J Exp Biol. 2008;211(24):3889–907.
Hirasaki E, Kumakura H. Head movements during locomotion in a gibbon and Japanese macaques. Neuroreport. 2004;15(4):643–7.
Loscher DM, Meyer F, Kracht K, Nyakatura JA. Timing of head movements is consistent with energy minimization in walking ungulates. Proc R Soc B. 2016;283:20161908.
Runciman RJ, Richmond FJ. Shoulder and forelimb orientations and loading in sitting cats: implications for head and shoulder movement. J Biomech. 1997;30(9):911–9.
Weishaupt M, Wiestner T, Peinen K, Waldern N, Roepstorff L, Weeren R, Meyer H, Johnston C. Effect of head and neck position on vertical ground reaction forces and interlimb coordination in the dressage horse ridden at walk and trot on a treadmill. Equine Vet J. 2006;38(S36):387–92.
Zsoldos R, Kotschwar A, Kotschwar A, Groesel M, Licka T, Peham C. Electromyography activity of the equine splenius muscle and neck kinematics during walk and trot on the treadmill. Equine Vet J. 2010;42(s38):455–61.
Zubair HN, Beloozerova IN, Sun H, Marlinski V. Head movement during walking in the cat. Neuroscience. 2016;332:101–20.
Eaton TH. Modifications of the shoulder girdle related to reach and stride in mammals. J Morphol. 1944;75(1):167–71.
Jenkins FA. The movement of the shoulder in claviculate and aclaviculate mammals. J Morphol. 1974;144(1):71–83.
Schmidt M, Voges D, Fischer MS. Shoulder movements during quadrupedal locomotion in arboreal primates. Z Morphol Anthropol. 2002;83:235–42.
Harmon LJ, Schulte JA, Larson A, Losos JB. Tempo and mode of evolutionary radiation in iguanian lizards. Science. 2003;301(5635):961–4.
Janis C. Tertiary mammal evolution in the context of changing climates, vegetation, and tectonic events. Annu Rev Ecol Syst. 1993;24(1):467–500.
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.
Timetree of Life Database [http://timetreebeta.igem.temple.edu/]. Accessed 22.11.2016.
Hedges SB, Marin J, Suleski M, Paymer M, Kumar S. Tree of life reveals clock-like speciation and diversification. Mol Biol Evol. 2015;32:835–45.
Csardi G, Nepusz T. The igraph software package for complex network research. Int J Complex Syst. 2006;1695(5):1–9.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2016.
Network cartography - Netcarto [http://seeslab.info/downloads/network-cartography-netcarto/]. Accessed 02.10.2016.
Guimera R, Amaral LAN. Functional cartography of complex metabolic networks. Nature. 2005;433(7028):895–900.
Adams DC. A generalized K statistic for estimating phylogenetic signal from shape and other high-dimensional multivariate data. Syst Biol. 2014;63(5):685–97.
Goolsby EW. Phylogenetic comparative methods for evaluating the evolutionary history of function-valued traits. Syst Biol. 2015;64(4):568–78.
Abouheif E. A method for testing the assumption of phylogenetic independence in comparative data. Evol Ecol Res. 1999;1(8):895–909.
Blomberg SP, Garland T Jr, Ives AR, Crespi B. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 2003;57(4):717–45.
Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3(2):217–23.
Clavel J, Escarguel G. Merceron G: mvMORPH: an R package for fitting evolutionary models to morphometric data. Methods Ecol Evol. 2015;6(11):1311–9.
Marwick B, Krishnamoorthy K: cvequality: Tests for the Equality of Coefficients of Variation from Multiple Groups. 2016. R package version 0.1.1. https://CRAN.R-project.org/package=cvequality.
Orme D, Freckleton R, Thomas G, Petzoldt T, Fritz S, Isaac N, Pearse W: caper: Comparative Analyses of Phylogenetics and Evolution in R. 2013. R package version 0.5.2. https://CRAN.R-project.org/package=caper.
Adams DC, Collyer ML, Kaliontzopoulou A, Sherratt E: Geomorph: Software for geometric morphometric analyses. 2013. R package version 3.0.5. https://CRAN.R-project.org/package=geomorph.
Newman ME, Girvan M. Finding and evaluating community structure in networks. Phys Rev E. 2004;69(2):026113.
Esteve-Altava B, Marugán-Lobón J, Botella H, Rasskin-Gutman D. Network models in anatomical systems. J Anthropol Sci. 2011;89:1–10.
Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24(1):129–31.
Slater GJ, Price SA, Santini F, Alfaro ME. Diversity versus disparity and the radiation of modern cetaceans. Proc R Soc Lond B Biol Sci. 2010;277(1697):3097–104.
Ercoli MD, Álvarez A, Busker F, Morales MM, Julik E, Smith HF, Adrian B, Barton M, Bhagavatula K, Poole M. Myology of the head, neck, and thoracic region of the lesser Grison (Galictis cuja) in comparison with the red panda (Ailurus fulgens) and other carnivorans: phylogenetic and functional implications. J Mamm Evol. 2017;24(3):289–322.
Ercoli MD, Álvarez A, Stefanini MI, Busker F, Morales MM. Muscular anatomy of the forelimbs of the lesser grison (Galictis cuja), and a functional and phylogenetic overview of Mustelidae and other Caniformia. J Mamm Evol. 2015;22(1):57–91.
We thank Heiko Stark for helpful discussions, Jan Wölfer for assistance with the phylogentic comparative methods, and Mikaela Lui for language improvement. We are grateful to associated editor Leandro Monteiro and two anonymous reviewers who improved the manuscript with their helpful suggestions.
PA is funded by the Max Planck Society. BE-A is funded by the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No 654155. The funding supported the experimental design and execution of the work described in this study. The funding institution had no influence on the analysis and interpretation of the data.
Availability of data and materials
All data generated and analyzed during this study are included in this published article and its additional information files. Raw data (i.e., the adjacency matrices) are provided in the additional information files.
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.
AF1 Additional information on the results of phylogenetic, network, and modularity analyses. (PDF 1108 kb)
AF2 Systematics and references of investigated species. (PDF 414 kb)
AF3 Additional information on methods. (PDF 210 kb)
AF4 Adjacency matrices coding for the topological information of neck’s musculoskeletal organization of investigated species. (XLS 677 kb)
About this article
Cite this article
Arnold, P., Esteve-Altava, B. & Fischer, M.S. Musculoskeletal networks reveal topological disparity in mammalian neck evolution. BMC Evol Biol 17, 251 (2017). https://doi.org/10.1186/s12862-017-1101-1