Kornelius Kupczik

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.

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.

Mandibular molar root morphology in Neanderthals and Late Pleistocene and recent Homo sapiens.

Neanderthals have a distinctive suite of dental features, including large anterior crown and root dimensions and molars with enlarged pulp cavities. Yet, there is little known about variation in molar root morphology in Neanderthals and other recent and fossil members of Homo. Here, we provide the first comprehensive metric analysis of permanent mandibular molar root morphology in Middle and Late Pleistocene Homo neanderthalensis, and Late Pleistocene (Aterian) and recent Homo sapiens. We specifically address the question of whether root form can be used to distinguish between these groups and assess whether any variation in root form can be related to differences in tooth function. We apply a microtomographic imaging approach to visualise and quantify the external and internal dental morphologies of both isolated molars and molars embedded in the mandible (n=127). Univariate and multivariate analyses reveal both similarities (root length and pulp volume) and differences (occurrence of pyramidal roots and dental tissue volume proportion) in molar root morphology among penecontemporaneous Neanderthals and Aterian H. sapiens. In contrast, the molars of recent H. sapiens are markedly smaller than both Pleistocene H. sapiens and Neanderthals, but share with the former the dentine volume reduction and a smaller root-to-crown volume compared with Neanderthals. Furthermore, we found the first molar to have the largest average root surface area in recent H. sapiens and Neanderthals, although in the latter the difference between M(1) and M(2) is small. In contrast, Aterian H. sapiens root surface areas peak at M(2). Since root surface area is linked to masticatory function, this suggests a distinct occlusal loading regime in Neanderthals compared with both recent and Pleistocene H. sapiens.

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.

The mechanical function of the periodontal ligament in the macaque mandible: a validation and sensitivity study using finite element analysis

Whilst the periodontal ligament (PDL) acts as an attachment tissue between bone and tooth, hypotheses regarding the role of the PDL as a hydrodynamic damping mechanism during intraoral food processing have highlighted its potential importance in finite element (FE) analysis. Although experimental and constitutive models have correlated the mechanical function of the PDL tissue with its anisotropic, heterogeneous, viscoelastic and non‐linear elastic nature, in many FE simulations the PDL is either present or absent, and when present is variably modelled. In addition, the small space the PDL occupies and the inability to visualize the PDL tissue using μCT scans poses issues during FE model construction and so protocols for the PDL thickness also vary. In this paper we initially test and validate the sensitivity of an FE model of a macaque mandible to variations in the Young’s modulus and the thickness of the PDL tissue. We then tested the validity of the FE models by carrying out experimental strain measurements on the same mandible in the laboratory using laser speckle interferometry. These strain measurements matched the FE predictions very closely, providing confidence that material properties and PDL thickness were suitably defined. The FE strain results across the mandible are generally insensitive to the absence and variably modelled PDL tissue. Differences are only found in the alveolar region adjacent to the socket of the loaded tooth. The results indicate that the effect of the PDL on strain distribution and/or absorption is restricted locally to the alveolar bone surrounding the teeth and does not affect other regions of the mandible.

Brief communication: Contributions of enamel-dentine junction shape and enamel deposition to primate molar crown complexity

Molar crown morphology varies among primates from relatively simple in some taxa to more complex in others, with such variability having both functional and taxonomic significance. In addition to the primary cusps, crown surface complexity derives from the presence of crests, cuspules, and crenulations. Developmentally, this complexity results from the deposition of an enamel cap over a basement membrane (the morphology of which is preserved as the enamel-dentine junction, or EDJ, in fully formed teeth). However, the relative contribution of the enamel cap and the EDJ to molar crown complexity is poorly characterized. In this study we examine the complexity of the EDJ and enamel surface of a broad sample of primate (including fossil hominin) lower molars through the application of micro-computed tomography and dental topographic analysis. Surface complexity of the EDJ and outer enamel surface (OES) is quantified by first mapping, and then summing, the total number of discrete surface orientation patches. We investigate the relative contribution of the EDJ and enamel cap to crown complexity by assessing the correlation in patch counts between the EDJ and OES within taxa and within individual teeth. We identify three patterns of EDJ/OES complexity which demonstrate that both crown patterning early in development and the subsequent deposition of the enamel cap contribute to overall crown complexity in primates.

Enamel Thickness in Bornean and Sumatran Orangutan Dentitions

Dental enamel thickness has received considerable attention in ecological models of the adaptive significance of primate morphology. Several authors have theorized that the degree of enamel thickness may reflect selective pressures related to the consumption of fallback foods (dietary items that may require complex processing and/or have low nutritional value) during times of preferred food scarcity. Others have speculated that enamel thickness reflects selection during mastication of foods with particular material properties (i.e., toughness and hardness). Orangutans prefer ripe fruit when available, but show interspecific and sex differences in the consumption of fallback foods (bark, leaves, and figs) and other preferred foods (certain seeds). Bornean orangutans (Pongo pygmaeus) have also been reported to masticate more mechanically demanding foods than Sumatran orangutans (Pongo abelii). To test these ecological models, we assessed two-dimensional enamel thickness in orangutan full dentitions using established histological and virtual quantification methods. No significant differences in average enamel thickness (AET) were found between species. We found significant differences in the components of enamel thickness indices between sexes, with males showing greater enamel-dentine junction lengths and dentine core areas, and thus relatively thinner enamel than females. Comparisons of individuals of known sex and species revealed a dentition-wide trend for Bornean females to show greater AET than Sumatran females. Differences between small samples of males were less evident. These data provide only limited support for ecological explanations of enamel thickness patterns within great ape genera. Future studies of dietary ecology and enamel thickness should consider sex differences more systematically.

Long anterior mandibular tooth roots in Neanderthals are not the result of their large jaws

Tooth root length has been shown to taxonomically distinguish Neanderthals from modern humans. However, this may result from differences in jaw size between both taxa, although most previous studies have revealed a very low or non-existent correlation between tooth size and jaw size in recent modern humans. We therefore investigated, within a broader taxonomical frame, to what extent measurements on the anterior tooth roots and the symphyseal region covary. Our samples comprise permanent mandibular incisors and canines from Mauer, Neanderthals, and extant and fossil modern humans sensu lato. Using micro-computed tomography, we took linear and cross-sectional surface area measurements of the roots and the symphyseal region and calculated the root volume. We also measured 3D landmarks to quantify the overall size of the mandible using centroid size. Furthermore, we analyzed the relationship between root size and symphyseal shape, based on Procrustes shape variables of semi-landmarks along the symphyseal outline. Our results show that Neanderthals have significantly larger anterior tooth roots than recent modern humans in terms of root length, mid-sagittal surface area and volume, even after correction for mandibular size. In contrast, symphyseal height and width do not differ significantly between both taxa, whereas, without scaling, the mid-sagittal symphyseal surface area and the centroid size of the mandible do differ. Importantly, no significant correlation was found between any of the root and symphyseal measurements after correction for overall mandibular size. The shape analyses revealed that Neanderthals have a vertical symphyseal profile with an evenly-thick symphysis, whereas recent modern humans display an unevenly-thick symphysis, comprising a pronounced incurvatio mandibularis and a bony chin. These results suggest a negative evolutionary allometry for the recent modern human anterior root size. Therefore, root length and other root dimensions can be considered taxonomically relevant for distinguishing Neanderthals from modern humans.

Reconstruction of muscle fascicle architecture from iodine-enhanced microCT images: A combined texture mapping and streamline approach.

Skeletal muscle models are used to investigate motion and force generation in both biological and bioengineering research. Yet, they often lack a realistic representation of the muscles internal architecture which is primarily composed of muscle fibre bundles, known as fascicles. Recently, it has been shown that fascicles can be resolved with micro-computed tomography (µCT) following staining of the muscle tissue with iodine potassium iodide (I2KI). Here, we present the reconstruction of the fascicular spatial arrangement and geometry of the superficial masseter muscle of a dog based on a combination of pattern recognition and streamline computation. A cadaveric head of a dog was incubated in I2KI and µCT-scanned. Following segmentation of the masseter muscle a statistical pattern recognition algorithm was applied to create a vector field of fascicle directions. Streamlines were then used to transform the vector field into a realistic muscle fascicle representation. The lengths of the reconstructed fascicles and the pennation angles in two planes (frontal and sagittal) were extracted and compared against a tracked fascicle field obtained through cadaver dissection. Both fascicle lengths and angles were found to vary substantially within the muscle confirming the complex and heterogeneous nature of skeletal muscle described by previous studies. While there were significant differences in the pennation angle between the experimentally derived and µCT-reconstructed data, there was congruence in the fascicle lengths. We conclude that the presented approach allows for embedding realistic fascicle information into finite element models of skeletal muscles to better understand the functioning of the musculoskeletal system.

Structure–function relations of primate lower incisors: a study of the deformation of Macaca mulatta dentition using electronic speckle pattern interferometry (ESPI)

Teeth adopt a variety of different morphologies, each of which is presumably optimized for performing specific functions during feeding. It is generally agreed that the enamel cap is a crucial element in controlling the mechanical behavior of mammalian teeth under load. Incisors are particularly interesting in terms of structure–function relations, as their role in feeding is that of the ‘first bite’. However, little is known how incisor cap morphology is related to tooth deformation. In the present paper we examine the mechanical behavior of mandibular central incisors in the cercopithecine primate Macaca mulatta under loads similar to those encountered during ingestion. We map three‐dimensional displacements on the labial surface of the crown as it is compressed, using electronic speckle pattern interferometry (ESPI), an optical metrology method. In addition, micro‐computed tomography is used to obtain data regarding the morphology of the enamel cap, which in the M. mulatta lower incisors exhibits missing or very little enamel on the lingual face. The results showed that although compressed along a longitudinal axis, deformation in the incisors mostly occurred in the lingual direction and orthogonal to the direction of the applied load. Both isolated, embedded teeth and teeth in the mandible showed considerable lingual deformation. Incisor deformation in the mandible was generally greater, reflecting the additional freedom of movement enabled by the supporting structures. We show that the association with adjacent teeth in the arch is significant for the behavior of the tooth under load. Finally, loading two teeth simultaneously in the mandible showed that they work as one functional unit. We suggest that these results demonstrate the importance of enamel cap morphology in directing deformation behavior; an ability stemming from the stiffness of the enamel cap overlying the more pliable dentin.