The developing bird pelvis passes through ancestral dinosaurian conditions.

Living birds (Aves) have bodies substantially modified from the ancestral reptilian condition. The avian pelvis in particular experienced major changes during the transition from early archosaurs to living birds1,2. This stepwise transformation is well documented by an excellent fossil record2-4; however, the ontogenetic alterations that underly it are less well understood. We used embryological imaging techniques to examine the morphogenesis of avian pelvic tissues in three dimensions, allowing direct comparison with the fossil record. Many ancestral dinosaurian features2 (for example, a forward-facing pubis, short ilium and pubic 'boot') are transiently present in the early morphogenesis of birds and arrive at their typical 'avian' form after transitioning through a prenatal developmental sequence that mirrors the phylogenetic sequence of character acquisition. We demonstrate quantitatively that avian pelvic ontogeny parallels the non-avian dinosaur-to-bird transition and provide evidence for phenotypic covariance within the pelvis that is conserved across Archosauria. The presence of ancestral states in avian embryos may stem from this conserved covariant relationship. In sum, our data provide evidence that the avian pelvis, whose early development has been little studied5-7, evolved through terminal addition-a mechanism8-10 whereby new apomorphic states are added to the end of a developmental sequence, resulting in expression8,11 of ancestral character states earlier in that sequence. The phenotypic integration we detected suggests a previously unrecognized mechanism for terminal addition and hints that retention of ancestral states in development is common during evolutionary transitions.

A new role for joint mobility in reconstructing vertebrate locomotor evolution.

Reconstructions of movement in extinct animals are critical to our understanding of major transformations in vertebrate locomotor evolution. Estimates of joint range of motion (ROM) have long been used to exclude anatomically impossible joint poses from hypothesized gait cycles. Here we demonstrate how comparative ROM data can be harnessed in a different way to better constrain locomotor reconstructions. As a case study, we measured nearly 600,000 poses from the hindlimb joints of the Helmeted Guineafowl and American alligator, which represent an extant phylogenetic bracket for the archosaurian ancestor and its pseudosuchian (crocodilian line) and ornithodiran (bird line) descendants. We then used joint mobility mapping to search for a consistent relationship between full potential joint mobility and the subset of joint poses used during locomotion. We found that walking and running poses are predictably located within full mobility, revealing additional constraints for reconstructions of extinct archosaurs. The inferential framework that we develop here can be expanded to identify ROM-based constraints for other animals and, in turn, will help to unravel the history of vertebrate locomotor evolution.

Inner workings of the alligator ankle reveal the mechanistic origins of archosaur locomotor diversity.

Major transformations in the locomotor system of archosaurs (a major clade of reptiles including birds, crocodiles, dinosaurs, and pterosaurs) were accompanied by significant modifications to ankle anatomy. How the evolution of such a complex multi-joint structure is related to shifts in ankle function and locomotor diversity across this clade remains unclear and weakly grounded in extant experimental data. Here, we used X-ray Reconstruction of Moving Morphology to reconstruct skeletal motion and quantify the sources of three-dimensional ankle mobility in the American alligator, a species that retains the ancestral archosaur ankle structure. We then applied the observed relationships between joint excursion and locomotor behaviors to predict ankle function in extinct archosaurs. High-resolution reconstructions of Alligator skeletal movement revealed previously unseen regionalized coordination among joints responsible for overall ankle rotation. Differences in joint contributions between maneuvers and steady walking parallel transitions in mobility inferred from the ankle structure of fossil taxa in lineages with more erect hind limb postures. Key ankle structures related to ankle mobility were identified in the alligator, which permitted the characterization of ancestral archosaur ankle function. Modifications of these structures provide morphological evidence for functional convergence among sublineages of bird-line and crocodylian-line archosaurs. Using the dynamic insight into the internal sources of Alligator ankle mobility and trends among locomotor modes, we trace anatomical shifts and propose a mechanistic hypothesis for the evolution of ankle structure and function across Archosauria.

A proposed standard for quantifying 3-D hindlimb joint poses in living and extinct archosaurs.

The last common ancestor of birds and crocodylians plus all of its descendants (clade Archosauria) dominated terrestrial Mesozoic ecosystems, giving rise to disparate body plans, sizes, and modes of locomotion. As in the fields of vertebrate morphology and paleontology more generally, studies of archosaur skeletal structure have come to depend on tools for acquiring, measuring, and exploring three-dimensional (3-D) digital models. Such models, in turn, form the basis for many analyses of musculoskeletal function. A set of shared conventions for describing 3-D pose (joint or limb configuration) and 3-D kinematics (change in pose through time) is essential for fostering comparison of posture/movement among such varied species, as well as for maximizing communication among scientists. Following researchers in human biomechanics, we propose a standard methodological approach for measuring the relative position and orientation of the major segments of the archosaur pelvis and hindlimb in 3-D. We describe the construction of anatomical and joint coordinate systems using the extant guineafowl and alligator as examples. Our new standards are then applied to three extinct taxa sampled from the wider range of morphological, postural, and kinematic variation that has arisen across >250 million years of archosaur evolution. These proposed conventions, and the founding principles upon which they are based, can also serve as starting points for measuring poses between elements within a hindlimb segment, for establishing coordinate systems in the forelimb and axial skeleton, or for applying our archosaurian system more broadly to different vertebrate clades.

Paleobiological reconstructions of articular function require all six degrees of freedom.

Paleobiologists typically exclude impossible joint poses from reconstructions of extinct animals by estimating the rotational range of motion (ROM) of fossil joints. However, this ubiquitous practice carries the assumption that osteological estimates of ROM consistently overestimate true joint mobility. Because studies founded on ROM-based exclusion have contributed substantially to our understanding of functional and locomotor evolution, it is critical that this assumption be tested. Here, we evaluate whether ROM-based exclusion is, as currently implemented, a reliable strategy. We measured the true mobilities of five intact cadaveric joints using marker-based X-ray Reconstruction of Moving Morphology and compared them to virtual osteological estimates of ROM made allowing (a) only all three rotational, (b) all three rotational and one translational, and (c) all three rotational and all three translational degrees of freedom. We found that allowing combinations of motions in all six degrees of freedom is necessary to ensure that true mobility is always successfully captured. In other words, failing to include joint translations in ROM analyses results in the erroneous exclusion of many joint poses that are possible in life. We therefore suggest that the functional and evolutionary conclusions of existing paleobiological reconstructions may be weakened or even overturned when all six degrees of freedom are considered. We offer an expanded methodological framework for virtual ROM estimation including joint translations and outline recommendations for future ROM-based exclusion studies.

Alligators employ intermetatarsal reconfiguration to modulate plantigrade ground contact.

Feet must mediate substrate interactions across an animal's entire range of limb poses used in life. Metatarsals, the 'bones of the sole', are the dominant pedal skeletal elements for most tetrapods. In plantigrade species that walk on the entirety of their sole, such as living crocodylians, intermetatarsal mobility offers the potential for a continuum of reconfiguration within the foot itself. Alligator hindlimbs are capable of postural extremes from a belly sprawl to a high walk to sharp turns - how does the foot morphology dynamically accommodate these diverse demands? We implemented a hybrid combination of marker-based and markerless X-ray reconstruction of moving morphology (XROMM) to measure 3D metatarsal kinematics in three juvenile American alligators (Alligator mississippiensis) across their locomotor and maneuvering repertoire on a motorized treadmill and flat-surfaced arena. We found that alligators adaptively conformed their metatarsals to the ground, maintaining plantigrade contact throughout a spectrum of limb placements with non-planar feet. Deformation of the metatarsus as a whole occurred through variable abduction (twofold range of spread) and differential metatarsal pitching (45 deg arc of skew). Internally, metatarsals also underwent up to 65 deg of long-axis rotation. Such reorientation, which correlated with skew, was constrained by the overlapping arrangement of the obliquely expanded metatarsal bases. Such a proximally overlapping metatarsal morphology is shared by fossil archosaurs and archosaur relatives. In these extinct taxa, we suggest that intermetatarsal mobility likely played a significant role in maintaining ground contact across plantigrade postural extremes.

Contrast-enhanced XROMM reveals in vivo soft tissue interactions in the hip of Alligator mississippiensis.

Extant archosaurs exhibit highly divergent articular soft tissue anatomies between avian and crocodilian lineages. However, the general lack of understanding of the dynamic interactions among archosaur joint soft tissues has hampered further inferences about the function and evolution of these joints. Here we use contrast-enhanced computed tomography to generate 3D surface models of the pelvis, femora, and hip joint soft tissues in an extant archosaur, the American alligator. The hip joints were then animated using marker-based X-Ray Reconstruction of Moving Morphology (XROMM) to visualize soft tissue articulation during forward terrestrial locomotion. We found that the anatomical femoral head of the alligator travels beyond the cranial extent of the bony acetabulum and does not act as a central pivot, as has been suggested for some extinct archosaurs. Additionally, the fibrocartilaginous surfaces of the alligator's antitrochanter and femoral neck remain engaged during hip flexion and extension, similar to the articulation between homologous structures in birds. Moreover, the femoral insertion of the ligamentum capitis moves dorsoventrally against the membrane-bound portion of the medial acetabular wall, suggesting that the inner acetabular foramen constrains the excursion of this ligament as it undergoes cyclical stretching during the step cycle. Finally, the articular surface of the femoral cartilage model interpenetrates with those of the acetabular labrum and antitrochanter menisci; we interpret such interpenetration as evidence of compressive deformation of the labrum and of sliding movement of the menisci. Our data illustrate the utility of XROMM for studying in vivo articular soft tissue interactions. These results also allow us to propose functional hypotheses for crocodilian hip joint soft tissues, expanding our knowledge of vertebrate connective tissue biology and the role of joint soft tissues in locomotor behavior.

A coordinate-system-independent method for comparing joint rotational mobilities.

Three-dimensional studies of range of motion currently plot joint poses in a 'Euler space' whose axes are angles measured in the joint's three rotational degrees of freedom. Researchers then compute the volume of a pose cloud to measure rotational mobility. However, pairs of poses that are equally different from one another in orientation are not always plotted equally far apart in Euler space. This distortion causes a single joint's mobility to change when measured based on different joint coordinate systems and precludes fair comparison among joints. Here, we present two alternative spaces inspired by a 16th century map projection - cosine-corrected and sine-corrected Euler spaces - that allow coordinate-system-independent comparison of joint rotational mobility. When tested with data from a bird hip joint, cosine-corrected Euler space demonstrated a 10-fold reduction in variation among mobilities measured from three joint coordinate systems. This new quantitative framework enables previously intractable, comparative studies of articular function.

It's in the loop: shared sub-surface foot kinematics in birds and other dinosaurs shed light on a new dimension of fossil track diversity.

The feet of ground-dwelling birds retain many features of their dinosaurian ancestry. Experiments with living species offer insights into the complex interplay among anatomy, kinematics and substrate during the formation of Mesozoic footprints. However, a key aspect of the track-making process, sub-surface foot movement, is hindered by substrate opacity. Here, we use biplanar X-rays to image guineafowl walking through radiolucent substrates of different consistency (solid, dry granular, firm to semi-liquid muds). Despite substantial kinematic variation, the foot consistently moves in a looping pattern below ground. As the foot sinks and then withdraws, the claws of the three main toes create entry and exit paths in different locations. Sampling these paths at incremental horizons captures two-dimensional features just as fossil tracks do, allowing depth-based zones to be characterized by the presence and relative position of digit impressions. Examination of deep, penetrative tracks from the Early Jurassic confirms that bipeds had an equivalent looping response to soft substrates approximately 200 Ma. Our integration of extant and extinct evidence demonstrates the influence of substrate properties on sinking depth and sub-surface foot motion, both of which are significant sources of track variation in the fossil record of dinosaurs.

A biplanar X-ray approach for studying the 3D dynamics of human track formation.

Recent discoveries have made hominin tracks an increasingly prevalent component of the human fossil record, and these data have the capacity to inform long-standing debates regarding the biomechanics of hominin locomotion. However, there is currently no consensus on how to decipher biomechanical variables from hominin tracks. These debates can be linked to our generally limited understanding of the complex interactions between anatomy, motion, and substrate that give rise to track morphology. These interactions are difficult to study because direct visualization of the track formation process is impeded by foot and substrate opacity. To address these obstacles, we developed biplanar X-ray and computer animation methods, derived from X-ray Reconstruction of Moving Morphology (XROMM), to analyze the 3D dynamics of three human subjects' feet as they walked across four substrates (three deformable muds and rigid composite panel). By imaging and reconstructing 3D positions of external markers, we quantified the 3D dynamics at the foot-substrate interface. Foot shape, specifically heel and medial longitudinal arch deformation, was significantly affected by substrate rigidity. In deformable muds, we found that depths measured across tracks did not directly reflect the motions of the corresponding regions of the foot, and that track outlines were not perfectly representative of foot size. These results highlight the complex, dynamic nature of track formation, and the experimental methods presented here offer a promising avenue for developing and refining methods for accurately inferring foot anatomy and gait biomechanics from fossil hominin tracks.

Mandibular and dental characteristics of Late Triassic mammaliaform Haramiyavia and their ramifications for basal mammal evolution.

As one of the earliest-known mammaliaforms, Haramiyavia clemmenseni from the Rhaetic (Late Triassic) of East Greenland has held an important place in understanding the timing of the earliest radiation of the group. Reanalysis of the type specimen using high-resolution computed tomography (CT) has revealed new details, such as the presence of the dentary condyle of the mammalian jaw hinge and the postdentary trough for mandibular attachment of the middle ear-a transitional condition of the predecessors to crown Mammalia. Our tests of competing phylogenetic hypotheses with these new data show that Late Triassic haramiyids are a separate clade from multituberculate mammals and are excluded from the Mammalia. Consequently, hypotheses of a Late Triassic diversification of the Mammalia that depend on multituberculate affinities of haramiyidans are rejected. Scanning electron microscopy study of tooth-wear facets and kinematic functional simulation of occlusion with virtual 3D models from CT scans confirm that Haramiyavia had a major orthal occlusion with the tallest lingual cusp of the lower molars occluding into the lingual embrasure of the upper molars, followed by a short palinal movement along the cusp rows alternating between upper and lower molars. This movement differs from the minimal orthal but extensive palinal occlusal movement of multituberculate mammals, which previously were regarded as relatives of haramiyidans. The disparity of tooth morphology and the diversity of dental functions of haramiyids and their contemporary mammaliaforms suggest that dietary diversification is a major factor in the earliest mammaliaform evolution.

3-D range of motion envelopes reveal interacting degrees of freedom in avian hind limb joints.

Measuring range of motion (ROM) is a valuable technique that can link bone morphology to joint function in both extant and extinct taxa. ROM results are commonly presented as tables or graphs of maxima and minima for each rotational degree of freedom. We investigate the interactions among three degrees of freedom using X-ray reconstruction of moving morphology (XROMM) to measure ROM of the main hind limb joints of Helmeted Guineafowl (Numida meleagris). By plotting each rotation on an axis, we generate three-dimensional ROM volumes or envelopes composed of hundreds of extreme joint positions for the hip, knee, and intertarsal joints. We find that the shapes of ROM volumes can be quite complex, and that the contribution of long-axis rotation is often substantial. Plotting in vivo poses from individual birds executing different behaviors shows varying use of potential rotational combinations within their ROM envelopes. XROMM can provide unprecedented high-resolution data on the spatial relationship of skeletal elements and thereby illuminate/elucidate the complex ways in which soft and hard tissues interact to produce functional joints. In joints with three rotational degrees of freedom, two-dimensional representations of ROM (flexion/extension and abduction/adduction) are incomplete.

The birth of a dinosaur footprint: subsurface 3D motion reconstruction and discrete element simulation reveal track ontogeny.

Locomotion over deformable substrates is a common occurrence in nature. Footprints represent sedimentary distortions that provide anatomical, functional, and behavioral insights into trackmaker biology. The interpretation of such evidence can be challenging, however, particularly for fossil tracks recovered at bedding planes below the originally exposed surface. Even in living animals, the complex dynamics that give rise to footprint morphology are obscured by both foot and sediment opacity, which conceals animal-substrate and substrate-substrate interactions. We used X-ray reconstruction of moving morphology (XROMM) to image and animate the hind limb skeleton of a chicken-like bird traversing a dry, granular material. Foot movement differed significantly from walking on solid ground; the longest toe penetrated to a depth of ∼5 cm, reaching an angle of 30° below horizontal before slipping backward on withdrawal. The 3D kinematic data were integrated into a validated substrate simulation using the discrete element method (DEM) to create a quantitative model of limb-induced substrate deformation. Simulation revealed that despite sediment collapse yielding poor quality tracks at the air-substrate interface, subsurface displacements maintain a high level of organization owing to grain-grain support. Splitting the substrate volume along "virtual bedding planes" exposed prints that more closely resembled the foot and could easily be mistaken for shallow tracks. DEM data elucidate how highly localized deformations associated with foot entry and exit generate specific features in the final tracks, a temporal sequence that we term "track ontogeny." This combination of methodologies fosters a synthesis between the surface/layer-based perspective prevalent in paleontology and the particle/volume-based perspective essential for a mechanistic understanding of sediment redistribution during track formation.

Validation of XMALab software for marker-based XROMM.

Marker-based XROMM requires software tools for: (1) correcting fluoroscope distortion; (2) calibrating X-ray cameras; (3) tracking radio-opaque markers; and (4) calculating rigid body motion. In this paper we describe and validate XMALab, a new open-source software package for marker-based XROMM (C++ source and compiled versions on Bitbucket). Most marker-based XROMM studies to date have used XrayProject in MATLAB. XrayProject can produce results with excellent accuracy and precision, but it is somewhat cumbersome to use and requires a MATLAB license. We have designed XMALab to accelerate the XROMM process and to make it more accessible to new users. Features include the four XROMM steps (listed above) in one cohesive user interface, real-time plot windows for detecting errors, and integration with an online data management system, XMAPortal. Accuracy and precision of XMALab when tracking markers in a machined object are ±0.010 and ±0.043 mm, respectively. Mean precision for nine users tracking markers in a tutorial dataset of minipig feeding was ±0.062 mm in XMALab and ±0.14 mm in XrayProject. Reproducibility of 3D point locations across nine users was 10-fold greater in XMALab than in XrayProject, and six degree-of-freedom bone motions calculated with a joint coordinate system were 3- to 6-fold more reproducible in XMALab. XMALab is also suitable for tracking white or black markers in standard light videos with optional checkerboard calibration. We expect XMALab to increase both the quality and quantity of animal motion data available for comparative biomechanics research.

Guineafowl with a twist: asymmetric limb control in steady bipedal locomotion.

In avian bipeds performing steady locomotion, right and left limbs are typically assumed to act out of phase, but with little kinematic disparity. However, outwardly appearing steadiness may harbor previously unrecognized asymmetries. Here, we present marker-based XROMM data showing that guineafowl on a treadmill routinely yaw away from their direction of travel using asymmetrical limb kinematics. Variation is most strongly reflected at the hip joints, where patterns of femoral long-axis rotation closely correlate to degree of yaw divergence. As yaw deviations increase, hip long-axis rotation angles undergo larger excursions and shift from biphasic to monophasic patterns. At large yaw angles, the alternately striding limbs exhibit synchronous external and internal femoral rotations of substantial magnitude. Hip coordination patterns resembling those used during sidestep maneuvers allow birds to asymmetrically modulate their mediolateral limb trajectories and thereby advance using a range of body orientations.

Long-axis rotation: a missing degree of freedom in avian bipedal locomotion.

Ground-dwelling birds are typically characterized as erect bipeds having hind limbs that operate parasagittally. Consequently, most previous research has emphasized flexion/extension angles and moments as calculated from a lateral perspective. Three-dimensional (3D) motion analyses have documented non-planar limb movements, but the skeletal kinematics underlying changes in foot orientation and transverse position remain unclear. In particular, long-axis rotation of the proximal limb segments is extremely difficult to measure with topical markers. Here, we present six degree of freedom skeletal kinematic data from maneuvering guineafowl acquired by marker-based XROMM (X-ray Reconstruction of Moving Morphology). Translations and rotations of the hips, knees, ankles and pelvis were derived from animated bone models using explicit joint coordinate systems. We distinguished sidesteps, sidestep yaws, crossover yaws, sidestep turns and crossover turns, but birds often performed a sequence of blended partial maneuvers. Long-axis rotation of the femur (up to 38 deg) modulated the foot's transverse position. Long-axis rotation of the tibiotarsus (up to 65 deg) also affected medio-lateral positioning, but primarily served to either re-orient a swing phase foot or yaw the body about a stance phase foot. Tarsometatarsal long-axis rotation was minimal, as was hip, knee and ankle abduction/adduction. Despite having superficially hinge-like joints, birds coordinate substantial long-axis rotations of the hips and knees to execute complex 3D maneuvers while striking a diversity of non-planar poses.

Articular surface interactions distinguish dinosaurian locomotor joint poses.

Our knowledge of vertebrate functional evolution depends on inferences about joint function in extinct taxa. Without rigorous criteria for evaluating joint articulation, however, such analyses risk misleading reconstructions of vertebrate animal motion. Here we propose an approach for synthesizing raycast-based measurements of 3-D articular overlap, symmetry, and congruence into a quantitative "articulation score" for any non-interpenetrating six-degree-of-freedom joint configuration. We apply our methodology to bicondylar hindlimb joints of two extant dinosaurs (guineafowl, emu) and, through comparison with in vivo kinematics, find that locomotor joint poses consistently have high articulation scores. We then exploit this relationship to constrain reconstruction of a pedal walking stride cycle for the extinct dinosaur Deinonychus antirrhopus, demonstrating the utility of our approach. As joint articulation is investigated in more living animals, the framework we establish here can be expanded to accommodate additional joints and clades, facilitating improved understanding of vertebrate animal motion and its evolution.

Technical note: A volumetric method for measuring the longitudinal arch of human tracks and feet.

Fossil footprints (i.e., tracks) were believed to document arch anatomical evolution, although our recent work has shown that track arches record foot kinematics instead. Analyses of track arches can thereby inform the evolution of human locomotion, although quantifying this 3-D aspect of track morphology is difficult. Here, we present a volumetric method for measuring the arches of 3-D models of human tracks and feet, using both Autodesk Maya and Blender software. The method involves generation of a 3-D object that represents the space beneath the longitudinal arch, and measurement of that arch object's geometry and spatial orientation. We provide relevant tools and guidance for users to apply this technique to their own data. We present three case studies to demonstrate potential applications. These include, (1) measuring the arches of static and dynamic human feet, (2) comparing the arches of human tracks with the arches of the feet that made them, and (3) direct comparisons of human track and foot arch morphology throughout simulated track formation. The volumetric measurement tool proved robust for measuring 3-D models of human tracks and feet, in static and dynamic contexts. This tool enables researchers to quantitatively compare arches of fossil hominin tracks, in order to derive biomechanical interpretations from them, and/or offers a different approach for quantifying foot morphology in living humans.

Three-dimensional skeletal kinematics of the shoulder girdle and forelimb in walking Alligator.

Crocodylians occupy a key phylogenetic position for investigations of archosaur locomotor evolution. Compared to the well-studied hindlimb, relatively little is known about the skeletal movements and mechanics of the forelimb. In this study, we employed manual markerless XROMM (X-ray Reconstruction Of Moving Morphology) to measure detailed 3-D kinematics of the shoulder girdle and forelimb bones of American alligators (Alligator mississippiensis) walking on a treadmill. Digital models of the interclavicle, scapulocoracoid, humerus, radius and ulna were created using a 3-D laser scanner. Models were articulated and aligned to simultaneously recorded frames of fluoroscopic and standard light video to reconstruct and measure joint motion. Joint coordinate systems were established for the coracosternal, glenohumeral and elbow joints. Our analysis revealed that the limb joints only account for about half of fore/aft limb excursion; the remaining excursion results from shoulder girdle movements and lateral bending of the vertebral column. Considerable motion of each scapulocoracoid relative to the vertebral column is consistent with coracosternal mobility. The hemisellar design of the glenohumeral joint permits some additional translation, or sliding in the fore-aft plane, but this movement does not have much of an effect on the distal excursion of the bone.