Publications

The musculoskeletal configuration of the mammalian pectoral limb has been heralded as a key anatomical feature leading to the adaptive radiation of mammals, but limb function in the non‐mammaliaform cynodont outgroup remains unresolved. Conflicting reconstructions of abducted and adducted posture are based on mutually incompatible interpretations of ambiguous osteology. We reconstruct the pectoral limb of the Triassic non‐mammaliaform cynodont Massetognathus pascuali in three dimensions, by combining skeletal morphology from micro‐computed tomography with muscle anatomy from an extended extant phylogenetic bracket. Conservative tests of maximum range of motion suggest a degree of girdle mobility, as well as substantial freedom at the shoulder and the elbow joints. The glenoid fossa supports a neutral pose in which the distal end of the humerus points 45° posterolaterally from the body wall, intermediate between classically ‘sprawling’ and ‘parasagittal’ limb postures. Massetognathus pascuali is reconstructed as having a near‐mammalian complement of shoulder muscles, including an incipient rotator cuff (m. subscapularis, m. infraspinatus, m. supraspinatus, and m. teres minor). Based on close inspection of the morphology of the glenoid fossa, we hypothesize a posture‐driven scenario for the evolution of the therian ball‐and‐socket shoulder joint. The musculoskeletal reconstruction presented here provides the anatomical scaffolding for more detailed examination of locomotor evolution in the precursors to mammals.

Traditional Hill-type muscle models, parameterized using high-quality experimental data, are often “too weak” to reproduce the joint torques generated by healthy adults during rapid, high force tasks. This study investigated whether the failure of these models to account for different types of motor units contributes to this apparent weakness; if so, muscle-driven simulations may rely on excessively high muscle excitations to generate a given force. We ran a series of forward simulations that reproduced measured ankle mechanics during cycling at five cadences ranging from 60 to 140 RPM. We generated both “nominal” simulations, in which an abstract ankle model was actuated by a 1-element Hill-type plantar flexor with a single contractile element (CE), and “test” simulations, in which the same model was actuated by a 2-element plantar flexor with two CEs that accounted for the force-generating properties of slower and faster motor units. We varied the total excitation applied to the 2-element plantar flexor between 60 and 105% of the excitation from each nominal simulation, and we varied the amount distributed to each CE between 0 and 100% of the total. Within this test space, we identified the excitation level and distribution, at each cadence, that best reproduced the plantar flexor forces generated in the nominal simulations. Our comparisons revealed that the 2-element model required substantially less total excitation than the 1-element model to generate comparable forces, especially at higher cadences. For instance, at 140 RPM, the required excitation was reduced by 23%. These results suggest that a 2-element model, in which contractile properties are “tuned” to represent slower and faster motor units, can increase the apparent strength and perhaps improve the fidelity of simulations of tasks with varying mechanical demands.

Existing "off-the-shelf" musculoskeletal models are problematic when simulating movements that involve substantial hip and knee flexion, such as the upstroke of pedalling, because they tend to generate excessive passive fibre force. The goal of this study was to develop a refined musculoskeletal model capable of simulating pedalling and fast running, in addition to walking, which predicts the activation patterns of muscles better than existing models. Specifically, we tested whether the anomalous co-activation of antagonist muscles, commonly observed in simulations, could be resolved if the passive forces generated by the underlying model were diminished. We refined the OpenSim™ model published by Rajagopal et al. (IEEE Trans Biomed Eng 63:1-1, 2016) by increasing the model's range of knee flexion, updating the paths of the knee muscles, and modifying the force-generating properties of eleven muscles. Simulations of pedalling, running and walking based on this model reproduced measured EMG activity better than simulations based on the existing model-even when both models tracked the same subject-specific kinematics. Improvements in the predicted activations were associated with decreases in the net passive moments; for example, the net passive knee moment during the upstroke of pedalling decreased from 36.9 N m (existing model) to 6.3 N m (refined model), resulting in a dramatic decrease in the co-activation of knee flexors. The refined model is available from SimTK.org and is suitable for analysing movements with up to 120° of hip flexion and 140° of knee flexion.

Birds fly effectively and maneuver nimbly by dynamically changing the shape of their wings during each wingbeat. These shape changes have yet to be quantified automatically at high temporal and spatial resolution. Therefore, we developed a custom 3D surface reconstruction method, which uses a high-speed camera to identify spatially encoded binary striped patterns that are projected on a flying bird. This non-invasive structured-light method allows automated 3D reconstruction of each stand-alone frame and can be extended to multiple views. We demonstrate this new technique by automatically reconstructing the dorsal surface of a parrotlet wing at 3200 fps during flapping flight. From this shape we analyze key parameters such as wing twist and angle of attack distribution. While our binary ‘single-shot’ algorithm is demonstrated by quantifying dynamic shape changes of a flying bird, it is generally applicable to moving animals, plants and deforming objects.

Digit reduction is a major trend that characterizes horse evolution, but its causes and consequences have rarely been quantitatively tested. Using beam analysis on fossilized centre metapodials, we tested how locomotor bone stresses changed with digit reduction and increasing body size across the horse lineage. Internal bone geometry was captured from 13 fossil horse genera that covered the breadth of the equid phylogeny and the spectrum of digit reduction and body sizes, from Hyracotherium to Equus. To account for the load-bearing role of side digits, a novel, continuous measure of digit reduction was also established—toe reduction index (TRI). Our results show that without accounting for side digits, three-toed horses as late as Parahippus would have experienced physiologically untenable bone stresses. Conversely, when side digits are modelled as load-bearing, species at the base of the horse radiation through Equus probably maintained a similar safety factor to fracture stress. We conclude that the centre metapodial compensated for evolutionary digit reduction and body mass increases by becoming more resistant to bending through substantial positive allometry in internal geometry. These results lend support to two historical hypotheses: that increasing body mass selected for a single, robust metapodial rather than several smaller ones; and that, as horse limbs became elongated, the cost of inertia from the side toes outweighed their utility for stabilization or load-bearing.

Hill-type models are ubiquitous in the field of biomechanics, providing estimates of a muscle’s force as a function of its activation state and its assumed force-length and force-velocity properties. However, despite their routine use, the accuracy with which Hill-type models predict the forces generated by muscles during submaximal, dynamic tasks remains largely unknown. This study compared human gastrocnemii forces predicted by Hill-type models to the forces estimated from ultrasound-based measures of tendon length changes and stiffness during cycling, over a range of loads and cadences. We tested both a traditional model, with one contractile element, and a differential model, with two contractile elements that accounted for independent contributions of slow and fast muscle fibres. Both models were driven by subject-specific, ultrasound-based measures of fascicle lengths, velocities, and pennation angles and by activation patterns of slow and fast muscle fibres derived from surface electromyographic recordings. The models predicted on average, 54 % the time-varying gastrocnemii forces estimated from the ultrasound-based methods. However, differences between predicted and estimated forces were smaller under low speed-high activation conditions, with models able to predict nearly 80 % of the gastrocnemii force over a complete pedal cycle. Additionally, the predictions from the Hill-type muscle models tested here showed that a similar pattern of force production could be achieved for most conditions with and without accounting for the independent contributions of different muscle fibre types.

Flying animals must successfully contend with obstacles in their natural
environments. Inspired by the robust manoeuvring abilities of flying animals,
unmanned aerial systems are being developed and tested to improve flight
control through cluttered environments. We previously examined steering
strategies that pigeons adopt to fly through an array of vertical obstacles
(VOs). Modelling VO flight guidance revealed that pigeons steer towards
larger visual gaps when making fast steering decisions. In the present experiments,
we recorded three-dimensional flight kinematics of pigeons as they
flew through randomized arrays of horizontal obstacles (HOs). We found
that pigeons still decelerated upon approach but flew faster through a
denser array of HOs compared with the VO array previously tested. Pigeons
exhibited limited steering and chose gaps between obstacles most aligned to
their immediate flight direction, in contrast to VO navigation that favoured
widest gap steering. In addition, pigeons navigated past the HOs with more
variable and decreased wing stroke span and adjusted their wing stroke
plane to reduce contact with the obstacles. Variability in wing extension,
stroke plane and wing stroke path was greater during HO flight. Pigeons
also exhibited pronounced head movements when negotiating HOs, which
potentially serve a visual function. These head-bobbing-like movements
were most pronounced in the horizontal (flight direction) and vertical directions,
consistent with engaging motion vision mechanisms for obstacle
detection. These results show that pigeons exhibit a keen kinesthetic sense of
their body and wings in relation to obstacles. Together with aerodynamic flapping
flight mechanics that favours vertical manoeuvring, pigeons are able to
navigate HOs using simple rules, with remarkable success.

Mechanistically linking movement behaviors and ecology is key to understanding the adaptive evolution of locomotion. Predator evasion, a behavior that enhances fitness, may depend upon short bursts or complex patterns of locomotion. However, such movements are poorly characterized by existing biomechanical metrics. We present methods based on the entropy measure of randomness from Information Theory to quantitatively characterize the unpredictability of non-steady-state locomotion. We then apply the method by examining sympatric rodent species whose escape trajectories differ in dimensionality. Unlike the speed-regulated gait use of cursorial animals to enhance locomotor economy, bipedal jerboa (family Dipodidae) gait transitions likely enhance maneuverability. In field-based observations, jerboa trajectories are significantly less predictable than those of quadrupedal rodents, likely increasing predator evasion ability. Consistent with this hypothesis, jerboas exhibit lower anxiety in open fields than quadrupedal rodents, a behavior that varies inversely with predator evasion ability. Our unpredictability metric expands the scope of quantitative biomechanical studies to include non-steady-state locomotion in a variety of evolutionary and ecologically significant contexts.