Form-function relationships underlie rapid dietary changes in a lizard.

Macroevolutionary changes such as variation in habitat use or diet are often associated with convergent, adaptive changes in morphology. However, it is still unclear how small-scale morphological variation at the population level can drive shifts in ecology such as observed at a macroevolutionary scale. Here, we address this question by investigating how variation in cranial form and feeding mechanics relate to rapid changes in diet in an insular lizard (Podarcis siculus) after experimental introduction into a new environment. We first quantified differences in the skull shape and jaw muscle architecture between the source and introduced population using three-dimensional geometric morphometrics and dissections. Next, we tested the impact of the observed variation in morphology on the mechanical performance of the masticatory system using computer-based biomechanical simulation techniques. Our results show that small differences in shape, combined with variation in muscle architecture, can result in significant differences in performance allowing access to novel trophic resources. The confrontation of these data with the already described macroevolutionary relationships between cranial form and function in these insular lizards provides insights into how selection can, over relatively short time scales, drive major changes in ecology through its impact on mechanical performance.

Modeling and control of anterior-posterior and medial-lateral sways in standing posture.

To study essential anterior-posterior and medial-lateral sways of the stance caused by rotational movements about the ankle and hip joints, a mathematical model is developed for the 3D postural kinematics and dynamics. The model is in the form of nonlinear differential-algebraic equations corresponding to a biomechanical system with holonomic constraints. A nonlinear feedback control law is further derived for stabilizing the upright stance, whilst eliminating internal torques induced by the constraints on postural movements. Numerical simulations of the model parametrized with experimental data of human body segments illustrate the performance of postural balancing with the proposed control.

Assessments of bilateral asymmetry with application in human skull analysis.

As a common feature, bilateral symmetry of biological forms is ubiquitous, but in fact rarely exact. In a setting of analytic geometry, bilateral symmetry is defined with respect to a point, line or plane, and the well-known notions of fluctuating asymmetry, directional asymmetry and antisymmetry are recast. A meticulous scheme for asymmetry assessments is proposed and explicit solutions to them are derived. An investigation into observational errors of points representing the geometric structure of an object offers a baseline reference for asymmetry assessment of the object. The proposed assessments are applicable to individual, part or all point pairs at both individual and collective levels. The exact relationship between the developed treatments and the widely used Procrustes method in asymmetry assessment is examined. An application of the proposed assessments to a large collection of human skull data in the form of 3D landmark coordinates finds: (a) asymmetry of most skulls is not fluctuating, but directional if measured about a plane fitted to shared landmarks or side landmarks for balancing; (b) asymmetry becomes completely fluctuating if one side of a skull could be slightly rotated and translated with respect to the other side; (c) female skulls are more asymmetric than male skulls. The methodology developed in this study is rigorous and transparent, and lays an analytical base for investigation of structural symmetries and asymmetries in a wide range of biological and medical applications.

A biomechanical analysis of prognathous and orthognathous insect head capsules: evidence for a many-to-one mapping of form to function.

Insect head shapes are remarkably variable, but the influences of these changes on biomechanical performance are unclear. Among 'basal' winged insects, such as dragonflies, mayflies, earwigs and stoneflies, some of the most prominent anatomical changes are the general mouthpart orientation, eye size and the connection of the endoskeleton to the head. Here, we assess these variations as well as differing ridge and sclerite configurations using modern engineering methods including multibody dynamics modelling and finite element analysis in order to quantify and compare the influence of anatomical changes on strain in particular head regions and the whole head. We show that a range of peculiar structures such as the genal/subgenal, epistomal and circumocular areas are consistently highly loaded in all species, despite drastically differing morphologies in species with forward-projecting (prognathous) and downward-projecting (orthognathous) mouthparts. Sensitivity analyses show that the presence of eyes has a negligible influence on head capsule strain if a circumocular ridge is present. In contrast, the connection of the dorsal endoskeletal arms to the head capsule especially affects overall head loading in species with downward-projecting mouthparts. Analysis of the relative strains between species for each head region reveals that concerted changes in head substructures such as the subgenal area, the endoskeleton and the epistomal area lead to a consistent relative loading for the whole head capsule and vulnerable structures such as the eyes. It appears that biting-chewing loads are managed by a system of strengthening ridges on the head capsule irrespective of the general mouthpart and head orientation. Concerted changes in ridge and endoskeleton configuration might allow for more radical anatomical changes such as the general mouthpart orientation, which could be an explanation for the variability of this trait among insects. In an evolutionary context, many-to-one mapping of strain patterns onto a relatively similar overall head loading indeed could have fostered the dynamic diversification processes seen in insects.

New insights into the biomechanics of Legg-Calvé-Perthes' disease: The Role of Epiphyseal Skeletal Immaturity in Vascular Obstruction.

OBJECTIVES: Legg-Calvé-Perthes' disease (LCP) is an idiopathic osteonecrosis of the femoral head that is most common in children between four and eight years old. The factors that lead to the onset of LCP are still unclear; however, it is believed that interruption of the blood supply to the developing epiphysis is an important factor in the development of the condition. METHODS: Finite element analysis modelling of the blood supply to the juvenile epiphysis was investigated to understand under which circumstances the blood vessels supplying the femoral epiphysis could become obstructed. The identification of these conditions is likely to be important in understanding the biomechanics of LCP. RESULTS: The results support the hypothesis that vascular obstruction to the epiphysis may arise when there is delayed ossification and when articular cartilage has reduced stiffness under compression. CONCLUSION: The findings support the theory of vascular occlusion as being important in the pathophysiology of Perthes disease.Cite this article: M. Pinheiro, C. A. Dobson, D. Perry, M. J. Fagan. New insights into the biomechanics of Legg-Calvé-Perthes' disease: The Role of Epiphyseal Skeletal Immaturity in Vascular Obstruction. Bone Joint Res 2018;7:148-156. DOI: 10.1302/2046-3758.72.BJR-2017-0191.R1.

Cranial sutures work collectively to distribute strain throughout the reptile skull.

The skull is composed of many bones that come together at sutures. These sutures are important sites of growth, and as growth ceases some become fused while others remain patent. Their mechanical behaviour and how they interact with changing form and loadings to ensure balanced craniofacial development is still poorly understood. Early suture fusion often leads to disfiguring syndromes, thus is it imperative that we understand the function of sutures more clearly. By applying advanced engineering modelling techniques, we reveal for the first time that patent sutures generate a more widely distributed, high level of strain throughout the reptile skull. Without patent sutures, large regions of the skull are only subjected to infrequent low-level strains that could weaken the bone and result in abnormal development. Sutures are therefore not only sites of bone growth, but could also be essential for the modulation of strains necessary for normal growth and development in reptiles.

Masticatory loadings and cranial deformation in Macaca fascicularis: a finite element analysis sensitivity study.

Biomechanical analyses are commonly conducted to investigate how craniofacial form relates to function, particularly in relation to dietary adaptations. However, in the absence of corresponding muscle activation patterns, incomplete muscle data recorded experimentally for different individuals during different feeding tasks are frequently substituted. This study uses finite element analysis (FEA) to examine the sensitivity of the mechanical response of a Macaca fascicularis cranium to varying muscle activation patterns predicted via multibody dynamic analysis. Relative to the effects of varying bite location, the consequences of simulated variations in muscle activation patterns and of the inclusion/exclusion of whole muscle groups were investigated. The resulting cranial deformations were compared using two approaches; strain maps and geometric morphometric analyses. The results indicate that, with bite force magnitude controlled, the variations among the mechanical responses of the cranium to bite location far outweigh those observed as a consequence of varying muscle activations. However, zygomatic deformation was an exception, with the activation levels of superficial masseter being most influential in this regard. The anterior portion of temporalis deforms the cranial vault, but the remaining muscles have less profound effects. This study for the first time systematically quantifies the sensitivity of an FEA model of a primate skull to widely varying masticatory muscle activations and finds that, with the exception of the zygomatic arch, reasonable variants of muscle loading for a second molar bite have considerably less effect on cranial deformation and the resulting strain map than does varying molar bite point. The implication is that FEA models of biting crania will generally produce acceptable estimates of deformation under load as long as muscle activations and forces are reasonably approximated. In any one FEA study, the biological significance of the error in applied muscle forces is best judged against the magnitude of the effect that is being investigated.

Improving the validation of finite element models with quantitative full-field strain comparisons.

The techniques used to validate finite element (FE) models against experimental results have changed little during the last decades, even though the traditional approach of using single point measurements from strain gauges has major limitations: the strain distribution across the surface is not captured and the accurate determination of strain gauge positions on the model surface is difficult if the 3D surface topography of the bone surface is not measured. The full-field strain measurement technique of digital speckle pattern interferometry (DSPI) can overcome these problems, but the potential of this technique has not yet been fully exploited in validation studies. Here we explore new ways of quantifying and visualising the variation in strain magnitudes and orientations within and between repeated DSPI measurements as well as between the DSPI measurements and FEA results. We show that our approach provides a much more comprehensive and accurate validation than traditional methods. The measurement repeatability and the correspondence between measured and predicted strains vary to a great degree within and between measurement areas. The two models used in this study predict the measured strain directions and magnitudes surprisingly well considering that homogeneous and isotropic mechanical properties were assigned to the models. However, the full-field comparisons also reveal some discrepancies between measured and predicted strains that are most probably caused by inaccuracies in the models' geometries and the degree of simplification of the modelled material properties.

Intraluminal thrombus has a selective influence on matrix metalloproteinases and their inhibitors (tissue inhibitors of matrix metalloproteinases) in the wall of abdominal aortic aneurysms.

BACKGROUND: The influence of intraluminal thrombus (ILT) on the proteolytic environment within the wall of an abdominal aortic aneurysm (AAA) is unknown. This is the first study to examine the correlation between ILT thickness and the levels of matrix metalloproteinases (MMPs) and their natural inhibitors (tissue inhibitors of matrix metalloproteinases [TIMPs]) within the adjacent AAA wall. METHODS: Thirty-five patients undergoing elective repair of AAAs were studied. A single full-thickness infrarenal aortic sample was obtained uniformly from the arteriotomy site from each patient. All samples were snap frozen and analyzed for total and active MMP 2, 8, and 9 and TIMP 1 and 2. Thrombus thickness at the specimen site was measured on the preoperative contrast computed tomographic angiograms. RESULTS: There was a statistically significant correlation between ILT thickness, concentration of TIMP 1, and active concentration of MMP 9. MMP 2 (active and total) and TIMP 2 demonstrated a positive correlation with ILT thickness, although not statistically significant. CONCLUSION: In this novel study, we found a significant positive correlation of ILT thickness with active MMP 9 and TIMP 1 concentration in the adjacent AAA wall, and this may have implications for AAA expansion and eventual rupture.

Finite element modelling of squirrel, guinea pig and rat skulls: using geometric morphometrics to assess sensitivity.

Rodents are defined by a uniquely specialized dentition and a highly complex arrangement of jaw-closing muscles. Finite element analysis (FEA) is an ideal technique to investigate the biomechanical implications of these specializations, but it is essential to understand fully the degree of influence of the different input parameters of the FE model to have confidence in the model's predictions. This study evaluates the sensitivity of FE models of rodent crania to elastic properties of the materials, loading direction, and the location and orientation of the models' constraints. Three FE models were constructed of squirrel, guinea pig and rat skulls. Each was loaded to simulate biting on the incisors, and the first and the third molars, with the angle of the incisal bite varied over a range of 45°. The Young's moduli of the bone and teeth components were varied between limits defined by findings from our own and previously published tests of material properties. Geometric morphometrics (GMM) was used to analyse the resulting skull deformations. Bone stiffness was found to have the strongest influence on the results in all three rodents, followed by bite position, and then bite angle and muscle orientation. Tooth material properties were shown to have little effect on the deformation of the skull. The effect of bite position varied between species, with the mesiodistal position of the biting tooth being most important in squirrels and guinea pigs, whereas bilateral vs. unilateral biting had the greatest influence in rats. A GMM analysis of isolated incisor deformations showed that, for all rodents, bite angle is the most important parameter, followed by elastic properties of the tooth. The results here elucidate which input parameters are most important when defining the FE models, but also provide interesting glimpses of the biomechanical differences between the three skulls, which will be fully explored in future publications.

Validation of a morphometric reconstruction technique applied to a juvenile pelvis.

Three-dimensional reconstructions of bone geometry from microCT (computed tomography) data are frequently used in biomechanical and finite element analyses. Digitization of bone models is usually a simple process for specimens with a complete geometry, but in instances of damage or disarticulation it can be very challenging. Subsequent to digitization, further imaging techniques are often required to estimate the geometry of missing bone or connecting cartilage. This paper presents an innovative approach to the reconstruction of incomplete scan data, to reproduce proper anatomical arrangements of bones, including absent connecting cartilaginous elements. Utilizing geometric morphometric tools, the reconstruction technique is validated through comparison of a reconstructed 9 year old pelvis, to the original CT data. A principal component analysis and an overlay of the two pelves provide a measure of the accuracy of the reconstructed model. Future work aims to investigate the biomechanical effects of any minor positional error on the bone's predicted structural properties through the use of finite element analysis.

The effects of the periodontal ligament on mandibular stiffness: a study combining finite element analysis and geometric morphometrics.

It is generally accepted that the periodontal ligament (PDL) plays a crucial role in transferring occlusal forces from the teeth to the alveolar bone. Studies using finite element analysis (FEA) have helped to better understand this role and show that the stresses and strains in the alveolar bone are influenced by whether and how PDL is included in FE models. However, when the overall distribution of stresses and strains in crania and mandibles are of interest, PDL is often not included in FE models, although little is known about how this affects the results. Here we study the effect of representing PDL as a layer of solid material with isotropic homogeneous properties in an FE model of a human mandible using a novel application of geometric morphometrics. The results show that the modelling of the PDL affects the deformation and thus strain magnitudes not only of the alveolar bone around the biting tooth, but that the whole mandible deforms differently under load. As a result, the strain in the mandibular corpus is significantly increased when PDL is included, while the strain in the bone beneath the biting tooth is reduced. These results indicate the importance of the PDL in FE studies. Thus we recommend that the PDL should be included in FE models of the masticatory apparatus, with tests to assess the sensitivity of the results to changes in the Young's modulus of the PDL material.

Isolated word recognition of silent speech using magnetic implants and sensors.

There are a number of situations where individuals wish to communicate verbally but are unable to use conventional means-so called 'silent speech'. These include speakers in noisy and covert situations as well as patients who have lost their voice as a result of a laryngectomy or similar procedure. This paper focuses on those who are unable to speak following a laryngectomy and assesses the possibility of speech recognition based on a magnetic implant/sensors system. Permanent magnets are placed on the tongue and lips and the changes in magnetic field resulting from movement during speech are monitored using a set of magnetic sensors. The sensor signals are compared to sets of pre-recorded templates using the dynamic time warping (DTW) method, and the best match is identified. Experimental trials are reported for subjects with intact larynx, typically using 500-1000 utterances used for speaker dependant training and testing. It is shown that recognition rates of over 90% are achievable for vocabularies of at least 57 isolated words: sufficient to drive command-and-control applications.

Feedback control from the jaw joints during biting: an investigation of the reptile Sphenodon using multibody modelling.

Sphenodon, a lizard-like reptile, is the only living representative of a group that was once widespread at the time of the dinosaurs. Unique jaw mechanics incorporate crushing and shearing motions to breakdown food, but during this process excessive loading could cause damage to the jaw joints and teeth. In mammals like ourselves, feedback from mechanoreceptors within the periodontal ligament surrounding the teeth is thought to modulate muscle activity and thereby minimise such damage. However, Sphenodon and many other tetrapods lack the periodontal ligament and must rely on alternative control mechanisms during biting. Here we assess whether mechanoreceptors in the jaw joints could provide feedback to control muscle activity levels during biting. We investigate the relationship between joint, bite, and muscle forces using a multibody computer model of the skull and neck of Sphenodon. When feedback from the jaw joints is included in the model, predictions agree well with experimental studies, where the activity of the balancing side muscles reduces to maintain equal and minimal joint forces. When necessary, higher, but asymmetric, joint forces associated with higher bite forces were achievable, but these are likely to occur infrequently during normal food processing. Under maximum bite forces associated with symmetric maximal muscle activation, peak balancing side joint forces were more than double those of the working side. These findings are consistent with the hypothesis that feedback similar to that used in the simulation is present in Sphenodon.

Comparison between in vivo and theoretical bite performance: using multi-body modelling to predict muscle and bite forces in a reptile skull.

In biomechanical investigations, geometrically accurate computer models of anatomical structures can be created readily using computed-tomography scan images. However, representation of soft tissue structures is more challenging, relying on approximations to predict the muscle loading conditions that are essential in detailed functional analyses. Here, using a sophisticated multi-body computer model of a reptile skull (the rhynchocephalian Sphenodon), we assess the accuracy of muscle force predictions by comparing predicted bite forces against in vivo data. The model predicts a bite force almost three times lower than that measured experimentally. Peak muscle force estimates are highly sensitive to fibre length, muscle stress, and pennation where the angle is large, and variation in these parameters can generate substantial differences in predicted bite forces. A review of theoretical bite predictions amongst lizards reveals that bite forces are consistently underestimated, possibly because of high levels of muscle pennation in these animals. To generate realistic bites during theoretical analyses in Sphenodon, lizards, and related groups we suggest that standard muscle force calculations should be multiplied by a factor of up to three. We show that bite forces increase and joint forces decrease as the bite point shifts posteriorly within the jaw, with the most posterior bite location generating a bite force almost double that of the most anterior bite. Unilateral and bilateral bites produced similar total bite forces; however, the pressure exerted by the teeth is double during unilateral biting as the tooth contact area is reduced by half.

Validating a voxel-based finite element model of a human mandible using digital speckle pattern interferometry.

Finite element analysis is a powerful tool for predicting the mechanical behaviour of complex biological structures like bones, but to be confident in the results of an analysis, the model should be validated against experimental data. In such validation experiments, the strains in the loaded bones are usually measured with strain gauges glued to the bone surface, but the use of strain gauges on bone can be difficult and provides only very limited data regarding surface strain distributions. This study applies the full-field strain measurement technique of digital speckle pattern interferometry to measure strains in a loaded human mandible and compares the results with the predictions of voxel-based finite element models of the same specimen. It is found that this novel strain measurement technique yields consistent, reliable measurements. Further, strains predicted by the finite element analysis correspond well with the experimental data. These results not only confirm the usefulness of this technique for future validation studies in the field of bone mechanics, but also show that the modelling approach used in this study is able to predict the experimental results very accurately.

Masticatory loading and bone adaptation in the supraorbital torus of developing macaques.

Research on the evolution and adaptive significance of primate craniofacial morphologies has focused on adult, fully developed individuals. Here, we investigate the possible relationship between the local stress environment arising from masticatory loadings and the emergence of the supraorbital torus in the developing face of the crab-eating macaque Macaca fascicularis. By using finite element analysis (FEA), we are able to evaluate the hypothesis that strain energy density (SED) magnitudes are high in subadult individuals with resulting bone growth in the supraorbital torus. We developed three micro-CT-based FEA models of M. fascicularis skulls ranging in dental age from deciduous to permanent dentitions and validated them against published experimental data. Applied masticatory muscle forces were estimated from physiological cross-sectional areas of macaque cadaveric specimens. The models were sequentially constrained at each working side tooth to simulate the variation of the bite point applied during masticatory function. Custom FEA software was used to solve the voxel-based models and SED and principal strains were computed. A physiological superposition SED map throughout the face was created by allocating to each element the maximum SED value from each of the load cases. SED values were found to be low in the supraorbital torus region throughout ontogeny, while they were consistently high in the zygomatic arch and infraorbital region. Thus, if the supraorbital torus arises to resist masticatory loads, it is either already adapted in each of our subadult models so that we do not observe high SED or a lower site-specific bone deposition threshold must apply.

Combined finite element and multibody dynamics analysis of biting in a Uromastyx hardwickii lizard skull.

Lizard skulls vary greatly in shape and construction, and radical changes in skull form during evolution have made this an intriguing subject of research. The mechanics of feeding have surely been affected by this change in skull form, but whether this is the driving force behind the change is the underlying question that we are aiming to address in a programme of research. Here we have implemented a combined finite element analysis (FEA) and multibody dynamics analysis (MDA) to assess skull biomechanics during biting. A skull of Uromastyx hardwickii was assessed in the present study, where loading data (such as muscle force, bite force and joint reaction) for a biting cycle were obtained from an MDA and applied to load a finite element model. Fifty load steps corresponding to bilateral biting towards the front, middle and back of the dentition were implemented. Our results show the importance of performing MDA as a preliminary step to FEA, and provide an insight into the variation of stress during biting. Our findings show that higher stress occurs in regions where cranial sutures are located in functioning skulls, and as such support the hypothesis that sutures may play a pivotal role in relieving stress and producing a more uniform pattern of stress distribution across the skull. Additionally, we demonstrate how varying bite point affects stress distributions and relate stress distributions to the evolution of metakinesis in the amniote skull.

Rigid-body analysis of a lizard skull: modelling the skull of Uromastyx hardwickii.

Lizard skulls vary greatly in their detailed morphology. Theoretical models and practical studies have posited a definite relationship between skull morphology and bite performance, but this can be difficult to demonstrate in vivo. Computer modelling provides an alternative approach, as long as hard and soft tissue components can be integrated and the model can be validated. An anatomically accurate three-dimensional computer model of an Uromastyx hardwickii skull was developed for rigid-body dynamic analysis. The Uromastyx jaw was first opened under motion control, and then muscle forces were applied to produce biting simulations where bite forces and joint forces were calculated. Bite forces comparable to those reported in the literature were predicted, and detailed muscular force information was produced along with additional information on the stabilizing role of temporal ligaments in late jaw closing.