Michael J. Fagan

The earliest evidence for anatomically modern humans in northwestern Europe

The earliest anatomically modern humans in Europe are thought to have appeared around 43,000–42,000 calendar years before present (43–42 kyr cal bp), by association with Aurignacian sites and lithic assemblages assumed to have been made by modern humans rather than by Neanderthals. However, the actual physical evidence for modern humans is extremely rare, and direct dates reach no farther back than about 41–39 kyr cal bp, leaving a gap. Here we show, using stratigraphic, chronological and archaeological data, that a fragment of human maxilla from the Kent’s Cavern site, UK, dates to the earlier period. The maxilla (KC4), which was excavated in 1927, was initially diagnosed as Upper Palaeolithic modern human. In 1989, it was directly radiocarbon dated by accelerator mass spectrometry to 36.4–34.7 kyr cal bp. Using a Bayesian analysis of new ultrafiltered bone collagen dates in an ordered stratigraphic sequence at the site, we show that this date is a considerable underestimate. Instead, KC4 dates to 44.2–41.5 kyr cal bp. This makes it older than any other equivalently dated modern human specimen and directly contemporary with the latest European Neanderthals, thus making its taxonomic attribution crucial. We also show that in 13 dental traits KC4 possesses modern human rather than Neanderthal characteristics; three other traits show Neanderthal affinities and a further seven are ambiguous. KC4 therefore represents the oldest known anatomically modern human fossil in northwestern Europe, fills a key gap between the earliest dated Aurignacian remains and the earliest human skeletal remains, and demonstrates the wide and rapid dispersal of early modern humans across Europe more than 40 kyr ago.

Assessing mechanical function of the zygomatic region in macaques: validation and sensitivity testing of finite element models.

Crucial to the interpretation of the results of any finite element analysis of a skeletal system is a test of the validity of the results and an assessment of the sensitivity of the model parameters. We have therefore developed finite element models of two crania of Macaca fascicularis and investigated their sensitivity to variations in bone material properties, the zygomatico‐temporal suture and the loading regimen applied to the zygomatic arch. Maximum principal strains were validated against data derived from ex vivo strain gauge experiments using non‐physiological loads applied to the macaque zygomatic arch. Elastic properties of the zygomatic arch bone and the zygomatico‐temporal suture obtained by nanoindentation resulted in a high degree of congruence between experimental and simulated strains. The findings also indicated that the presence of a zygomatico‐temporal suture in the model produced strains more similar to experimental values than a completely separated or fused arch. Strains were distinctly higher when the load was applied through the modelled superficial masseter compared with loading an array of nodes on the arch. This study demonstrates the importance of the accurate selection of the material properties involved in predicting strains in a finite element model. Furthermore, our findings strongly highlight the influence of the presence of craniofacial sutures on strains experienced in the face. This has implications when investigating craniofacial growth and masticatory function but should generally be taken into account in functional analyses of the craniofacial system of both extant and extinct species.

Assessment of the role of sutures in a lizard skull: a computer modelling study

Sutures form an integral part of the functioning skull, but their role has long been debated among vertebrate morphologists and palaeontologists. Furthermore, the relationship between typical skull sutures, and those involved in cranial kinesis, is poorly understood. In a series of computational modelling studies, complex loading conditions obtained through multibody dynamics analysis were imposed on a finite element model of the skull of Uromastyx hardwickii, an akinetic herbivorous lizard. A finite element analysis (FEA) of a skull with no sutures revealed higher patterns of strain in regions where cranial sutures are located in the skull. From these findings, FEAs were performed on skulls with sutures (individual and groups of sutures) to investigate their role and function more thoroughly. Our results showed that individual sutures relieved strain locally, but only at the expense of elevated strain in other regions of the skull. These findings provide an insight into the behaviour of sutures and show how they are adapted to work together to distribute strain around the skull. Premature fusion of one suture could therefore lead to increased abnormal loading on other regions of the skull causing irregular bone growth and deformities. This detailed investigation also revealed that the frontal–parietal suture of the Uromastyx skull played a substantial role in relieving strain compared with the other sutures. This raises questions about the original role of mesokinesis in squamate evolution.

Combining geometric morphometrics and functional simulation: an emerging toolkit for virtual functional analyses

The development of virtual methods for anatomical reconstruction and functional simulation of skeletal structures offers great promise in evolutionary and ontogenetic investigations of form–function relationships. Key developments reviewed here include geometric morphometric methods for the analysis and visualization of variations in form (size and shape), finite element methods for the prediction of mechanical performance of skeletal structures under load and multibody dynamics methods for the simulation and prediction of musculoskeletal function. These techniques are all used in studies of form and function in biology, but only recently have they been combined in novel ways to facilitate biomechanical modelling that takes account of variations in form, can statistically compare performance, and relate performance to form and its covariates. Here we provide several examples that illustrate how these approaches can be combined and we highlight areas that require further investigation and development before we can claim a mature theory and toolkit for a statistical biomechanical framework that unites these methods.

Functional Evolution of the Feeding System in Rodents

The masticatory musculature of rodents has evolved to enable both gnawing at the incisors and chewing at the molars. In particular, the masseter muscle is highly specialised, having extended anteriorly to originate from the rostrum. All living rodents have achieved this masseteric expansion in one of three ways, known as the sciuromorph, hystricomorph and myomorph conditions. Here, we used finite element analysis (FEA) to investigate the biomechanical implications of these three morphologies, in a squirrel, guinea pig and rat. In particular, we wished to determine whether each of the three morphologies is better adapted for either gnawing or chewing. Results show that squirrels are more efficient at muscle-bite force transmission during incisor gnawing than guinea pigs, and that guinea pigs are more efficient at molar chewing than squirrels. This matches the known diet of nuts and seeds that squirrels gnaw, and of grasses that guinea pigs grind down with their molars. Surprisingly, results also indicate that rats are more efficient as well as more versatile feeders than both the squirrel and guinea pig. There seems to be no compromise in biting efficiency to accommodate the wider range of foodstuffs and the more general feeding behaviour adopted by rats. Our results show that the morphology of the skull and masticatory muscles have allowed squirrels to specialise as gnawers and guinea pigs as chewers, but that rats are high-performance generalists, which helps explain their overwhelming success as a group.

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.

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 Youngs modulus of the PDL material.

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.

Predicting skull loading: applying multibody dynamics analysis to a macaque skull.

Evaluating stress and strain fields in anatomical structures is a way to test hypotheses that relate specific features of facial and skeletal morphology to mechanical loading. Engineering techniques such as finite element analysis are now commonly used to calculate stress and strain fields, but if we are to fully accept these methods we must be confident that the applied loading regimens are reasonable. Multibody dynamics analysis (MDA) is a relatively new three dimensional computer modeling technique that can be used to apply varying muscle forces to predict joint and bite forces during static and dynamic motions. The method ensures that equilibrium of the structure is maintained at all times, even for complex statically indeterminate problems, eliminating nonphysiological constraint conditions often seen with other approaches. This study describes the novel use of MDA to investigate the influence of different muscle representations on a macaque skull model (Macaca fascicularis), where muscle groups were represented by either a single, multiple, or wrapped muscle fibers. The impact of varying muscle representation on stress fields was assessed through additional finite element simulations. The MDA models highlighted that muscle forces varied with gape and that forces within individual muscle groups also varied; for example, the anterior strands of the superficial masseter were loaded to a greater extent than the posterior strands. The direction of the muscle force was altered when temporalis muscle wrapping was modeled, and was coupled with compressive contact forces applied to the frontal, parietal and temporal bones of the cranium during biting. Anat Rec, 291:491–501, 2008.

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.