Mandibular corpus bone strain in goats and alpacas: implications for understanding the biomechanics of mandibular form in selenodont artiodactyls

Does the model reflect the system? When two-dimensional biomechanics is not 'good enough'

Models are mathematical representations of systems, processes or phenomena. In biomechanics, finite-element modelling (FEM) can be a powerful tool, allowing biologists to test form-function relationships in silico, replacing or extending results of in vivo experimentation. Although modelling simplifications and assumptions are necessary, as a minimum modelling requirement the results of the simplified model must reflect the biomechanics of the modelled system. In cases where the three-dimensional mechanics of a structure are important determinants of its performance, simplified two-dimensional modelling approaches are likely to produce inaccurate results. The vertebrate mandible is one among many three-dimensional anatomical structures routinely modelled using two-dimensional FE analysis. We thus compare the stress regimes of our published three-dimensional model of the chimpanzee mandible with a published two-dimensional model of the chimpanzee mandible and identify several fundamental differences. We then present a series of two-dimensional and three-dimensional FE modelling experiments that demonstrate how three key modelling parameters, (i) dimensionality, (ii) symmetric geometry, and (iii) constraints, affect deformation and strain regimes of the models. Our results confirm that, in the case of the primate mandible (at least), two-dimensional FEM fails to meet this minimum modelling requirement and should not be used to draw functional, ecological or evolutionary conclusions.

Mechanical properties of food and masticatory behavior in llamas, Llama glama

Mammals typically process food items more extensively in their oral cavities than do other vertebrates. Dental morphology, jaw-muscle activity patterns, mandibular movements, and tongue manipulation work to facilitate oral fragmentation of dietary items. While processing mechanically challenging foods, mammals modulate mandibular movements and bite forces via recruitment of greater jaw-adductor muscle forces and protracted biting or chewing. Because jaw-loading patterns are influenced by magnitude; frequency; and duration of muscular, bite, and reaction forces during routine feeding behaviors, relatively larger jaws are thought to be more characteristic of mammals that experience higher masticatory loads due to the processing of mechanically challenging foods. The ease of food fracture during post-canine biting and chewing is mainly determined by food stiffness and toughness. Such foods have been associated with increased loading magnitude and/or greater amounts of cyclical loading (i.e., chewing duration). Dietary properties are thought to modulate cyclical loading through changes in chewing frequency and chewing investment. On the other hand, chewing frequency has been found to be independent of dietary properties in rabbits and primates; however, little evidence exists regarding the influence of dietary properties on these parameters in a broader range of mammals. Here, we assessed chewing behavior in seven adult llamas (Llama glama) processing foods with a wide range of mechanical properties (grain, hay, carrots, and dried corn). Each subject was filmed at 60 frames/s, with video slowed for frame-by-frame computer analysis to obtain length of feeding bout and number of chewing cycles for each food type. These parameters were used to calculate chewing frequency (chews/s), chewing investment (chews/g), and chewing duration (s/g). Chewing frequency was not significantly related to mechanical properties of food, but chewing investment and chewing duration were significantly related to dietary stiffness and toughness. Therefore, cyclical loading is positively influenced by stiff and tough foods. This suggests that variation in jaw morphology in extinct and extant mammals is positively related to dietary stiffness and toughness, which requires greater chewing investment and increased chewing duration.

Mandibular symphyseal fusion in fossil primates: Insights from correlated patterns of jaw shape and masticatory function in living primates

OBJECTIVES

Variation in primate masticatory form and function has been extensively researched through both morphological and experimental studies. As a result, symphyseal fusion in different primate clades has been linked to either the recruitment of vertically directed balancing-side muscle force, the timing and recruitment of transversely directed forces, or both. This study investigates the relationship between jaw muscle activity patterns and morphology in extant primates to make inferences about masticatory function in extinct primates, with implications for understanding the evolution of symphyseal fusion.

MATERIALS AND METHODS

Three-dimensional mandibular landmark data were collected for 31 extant primates and nine fossil anthropoids and subfossil lemur species. Published electromyography (EMG) data were available for nine of the extant primate species. Partial least squares analysis and phylogenetic partial least squares analysis were used to identify relationships between EMG and jaw shape data and evaluate variation in jaw morphology.

RESULTS

Primates with partial and complete symphyseal fusion exhibit shape-function patterns associated with the wishboning motor pattern and loading regime, in contrast to shape-function patterns of primates with unfused jaws. All fossil primates examined (except Apidium) exhibit jaw morphologies suggestive of the wishboning motor pattern demonstrated in living anthropoids and indriids.

DISCUSSION

Partial fusion in Catopithecus, similar to indriids and some subfossil lemurs, may be sufficient to resist, or transfer, some amounts of transversely directed balancing-side muscle force at the symphysis, representing a transition to greater reliance on transverse jaw movement during mastication. Furthermore, possible functional convergences in physiological patterns during chewing (i.e., Archaeolemur) are identified.

Mandibular biomechanics after marginal resection: Correspondences of simulated volumetric strain and skeletal resorption

Serious mandibular diseases such as tumor or osteonecrosis often require segmental or marginal mandibulectomy, the latter with improved outcome thanks to preserved mandibular continuity. Nevertheless, gradual osteolytic and/or osteosclerotic skeletal changes frequently indicate repetitive resections. Based on the fundamental adaptivity of bone to mechanical loads, the question arose whether resection-related anatomical alterations trigger relevant pathological skeletal adaptations. For a clinical case after mandibular box resection due to progressive osteoradionecrosis (ORN), routine biomechanical loading was simulated by finite element method, respecting pathology-related anatomy, tissue properties, and biting capacity. By 3D-visualization of the mandible's pathological development from follow-up-CT's over four years, remarkable correspondences of skeletal resorptions and increased unphysiological strain were revealed. Higher unphysiological load was correlated with more serious and earlier skeletal alterations. Three months post-operatively, serious buccal destruction at the distal resection corner occurred in correspondence with dominant tensile strain. At the resection, elevated strain caused by reduced alveolar height corresponded to skeletal compromise, observed 8-9 months post-operatively. ORN-related lesions, diagnosed before resection, entailed unphysiological strain coinciding with local skeletal alterations. Simulations with "healthy" instead of pathological tissue coefficients induced quantitative improvements of 25-33%, but without fundamental change. These results suggest a decisive contribution of resection-related biomechanical skeletal adaptations to this patient's mandibular decline with hemimandibulectomy about 2.5 years after the first resection. However, mechanical stress concentrations in sharp angles as the distal resection corner and reduced stability due to decreased alveolar height generally bear the danger of pathological biomechanics and severe skeletal adaptations for patients after mandibular box resection.

An ontogenetic perspective on symphyseal fusion, occlusion and mandibular loading in alpacas (Vicugna pacos)

A primary hypothesis for the evolution of mandibular symphyseal fusion in some mammals is that it functions to resist loads incurred during routine mastication. Anecdotal support for this hypothesis is based on the fact that when the symphysis fuses, it typically does so early during postnatal ontogeny prior to or around the time of weaning. However, little is known about the process of fusion, particularly relative to feeding behaviors and the dynamics of mastication, including occlusion and masticatory loading. In the present study, we investigate the timing and process of symphyseal fusion in alpacas (Vicugna pacos) in the context of maturation of the oral apparatus and oral behavior. We also report on in vivo strains from the symphysis and corpus in young alpacas prior to and following full fusion and M₁ occlusion. Results show that fusion begins rostrally by 1 month and is complete by 6-7 months whereas all deciduous premolars and M₁ come into occlusion by 6 months. Although symphyseal loading patterns are maintained throughout ontogeny, in young alpacas symphyseal strain magnitudes are low compared with adults but corpus strain magnitudes are comparable to those found in adults. Reduced symphyseal loading in young individuals is contrary to what might be predicted given that the symphysis is still fusing. When considered in light of the development of occlusion and rumination, strain magnitudes may be necessarily low and reflect an overall delay in the maturation of masticatory dynamics.

In vitro bone strain distributions in a sample of primate pelves

The pelvis is a critical link in the hindlimb locomotor system and has a central role in resisting loads associated with locomotion, but our understanding of its structural biomechanics is quite limited. Empirical data on how the pelvis responds to the loads it encounters are important for understanding pelvic adaptation to locomotion, and for testing hypotheses regarding how the pelvis is adapted to its mechanical demands. This paper presents in vitro strain gauge data on a sample of monkey and ape cadaveric specimens (Macaca, Papio, Ateles, Hylobates), and assesses strain magnitudes and distributions through the bones of the pelvis: the ilium, ischium and pubis. Pelves were individually mounted in a materials testing system, loads were applied across three hindlimb angular positions, and strains were recorded from 18 locations on the pelvic girdle. Peak principal strains range from 2000 to 3000 με, similar to peak strains recorded from other mammals in vivo. Although previous work has suggested that the bones of the pelvis may act as bent beams, this study suggests that there are likely additional loading regimes superimposed on bending. Specifically, these data suggest that the ilium is loaded in axial compression and torsion, the ischium in torsion, the pubic rami in mediolateral bending, and the pubic symphysis is loaded in a combination of compression and torsion. Compressive strains dominate the pelves of all species representatives. Shear strains change with limb position; hip flexion at 45° induces smaller shear strains than mid-stance (90°) or hip extension (105°). The pelvic girdle is a complex structure that does not lend itself easily to modeling, but finite element analyses may prove useful to generate and refine hypotheses of pelvic biomechanics.

In vivo bone strain and finite element modeling of the mandible of Alligator mississippiensis

Forces experienced during feeding are thought to strongly influence the morphology of the vertebrate mandible; in vivo strain data are the most direct evidence for deformation of the mandible induced by these loading regimes. Although many studies have documented bone strains in the mammalian mandible, no information is available on strain magnitudes, orientations or patterns in the sauropsid lower jaw during feeding. Furthermore, strain gage experiments record the mechanical response of bone at a few locations, not across the entire mandible. In this paper, we present bone strain data recorded at various sites on the lower jaw of Alligator mississippiensis during in vivo feeding experiments. These data are used to understand how changes in loading regime associated with changes in bite location are related to changes in strain regime on the working and balancing sides of the mandible. Our results suggest that the working side mandible is bent dorsoventrally and twisted about its long-axis during biting, and the balancing side experiences primarily dorsoventral bending. Strain orientations are more variable on the working side than on the balancing side with changes in bite point and between experiments; the balancing side exhibits higher strain magnitudes. In the second part of this paper, we use principal strain orientations and magnitudes recorded in vivo to evaluate a finite element model of the alligator mandible. Our comparison demonstrates that strain orientations and mandibular deformation predicted by the model closely match in vivo results; however, absolute strain magnitudes are lower in the finite element model.

The role of the sutures in biomechanical dynamic simulation of a macaque cranial finite element model: implications for the evolution of craniofacial form

The global biomechanical impact of cranial sutures on the face and cranium during dynamic conditions is not well understood. It is hypothesized that sutures act as energy absorbers protecting skulls subjected to dynamic loads. This hypothesis predicts that sutures have a significant impact on global patterns of strain and cranial structural stiffness when analyzed using dynamic simulations; and that this global impact is influenced by suture material properties. In a finite element model developed from a juvenile Rhesus macaque cranium, five different sets of suture material properties for the zygomaticotemporal sutures were tested. The static and dynamic analyses produced similar results in terms of strain patterns and reaction forces, indicating that the zygomaticotemporal sutures have limited impact on global skull mechanics regardless of loading design. Contrary to the functional hypothesis tested in this study, the zygomaticotemporal sutures did not absorb significant amounts of energy during dynamic simulations regardless of loading speed. It is alternatively hypothesized that sutures are mechanically significant only insofar as they are weak points on the cranium that must be shielded from unduly high stresses so as not to disrupt vitally important growth processes. Thus, sutural and overall cranial form in some vertebrates may be optimized to minimize or otherwise modulate sutural stress and strain.

A preliminary analysis of correlations between chewing motor patterns and mandibular morphology across mammals

The establishment of a publicly-accessible repository of physiological data on feeding in mammals, the Feeding Experiments End-user Database (FEED), along with improvements in reconstruction of mammalian phylogeny, significantly improves our ability to address long-standing questions about the evolution of mammalian feeding. In this study, we use comparative phylogenetic methods to examine correlations between jaw robusticity and both the relative recruitment and the relative time of peak activity for the superficial masseter, deep masseter, and temporalis muscles across 19 mammalian species from six orders. We find little evidence for a relationship between jaw robusticity and electromyographic (EMG) activity for either the superficial masseter or temporalis muscles across mammals. We hypothesize that future analyses may identify significant associations between these physiological and morphological variables within subgroups of mammals that share similar diets, feeding behaviors, and/or phylogenetic histories. Alternatively, the relative peak recruitment and timing of the balancing-side (i.e., non-chewing-side) deep masseter muscle (BDM) is significantly negatively correlated with the relative area of the mandibular symphysis across our mammalian sample. This relationship exists despite BDM activity being associated with different loading regimes in the symphyses of primates compared to ungulates, suggesting a basic association between magnitude of symphyseal loads and symphyseal area among these mammals. Because our sample primarily represents mammals that use significant transverse movements during chewing, future research should address whether the correlations between BDM activity and symphyseal morphology characterize all mammals or should be restricted to this "transverse chewing" group. Finally, the significant correlations observed in this study suggest that physiological parameters are an integrated and evolving component of feeding across mammals.

Intraspecific scaling of chewing cycle duration in three species of domestic ungulates

In mammals, chewing cycle duration (CCD) increases with various measures of size, scaling with body mass0.13–0.28 and jaw length0.55. Proposed explanations for these scaling relationships include the allometry of body size, basal metabolic rate and tooth size, on the one hand, and pendular mechanics treating the jaw as a gravity-driven pendulum, on the other. Little is known, however, about the relationship between CCD and size within species. Recent research in dogs demonstrates altogether different scaling exponents and weaker correlations. This research suggests that breed-specific growth rates influence the maturation of the neural networks generating chewing rhythm, which may be altered because of changes in jaw mass during early postnatal growth. Here, we explored the intraspecific scaling of CCD within a sample of adult horses ranging from miniatures to draft breeds and an ontogenetic sample of goats and alpacas from infants to adults. In horses, CCD scales with body mass0.19 and jaw length0.57, although in neither case is the correlation significant. In the ontogenetic samples of goats and alpacas, CCD is significantly correlated with body mass, scaling as CCD∝body mass0.37 in both species. In goats, but not alpacas, CCD is also significantly correlated with jaw length, scaling as jaw length1.032. As in dogs, the scaling of CCD in horses may reflect the influence of selective breeding on growth trajectories of different breeds, resulting in reduced body and jaw size differences among infants, when CCD is established, compared with adults. However, the allometric scaling of tooth size in horses of different breeds may be a potential influence on the scaling of CCD. The scaling of CCD with body and jaw size in goats, and to a lesser extent in alpacas, also suggests that the development of peripheral masticatory structures such as the teeth and occlusal relations may play a role in changes in CCD during the earliest stages of postnatal ontogeny.

Allometry of masticatory loading parameters in mammals

Considerable research on the scaling of loading patterns in mammalian locomotor systems has not been accompanied by a similarly comprehensive analysis of the interspecific scaling of loading regimes in the mammalian masticatory complex. To address this deficiency, we analyzed mandibular corpus bone strain in 11 mammalian taxa varying in body size by over 2.5 orders of magnitude, including goats, horses, alpacas, pigs, and seven primate taxa. During alert chewing and biting of hard/tough foods, bone-strain data were collected with rosette gauges placed along the lateral aspect of the mandibular corpus below the molars or premolars. Bone-strain data were used to characterize relevant masticatory loading parameters: peak loading magnitudes, chewing cycle duration, chewing frequency, occlusal duty factor, loading rate, and loading time. Interspecific analyses indicate that much as observed in limb elements, corpus peak-strain magnitudes are similar across mammals of disparate body sizes. Chewing frequency is inversely correlated with body size, much as with locomotor stride frequency. Some of this allometric variation in chewing frequency appears to be due to a negative correlation with loading time, which increases with body size. Similar to the locomotor apparatus, occlusal duty factor, or the duration of the chewing cycle during which the corpus is loaded, does not vary with body size. Peak principal-strain magnitudes are most strongly positively correlated with loading rate and only secondarily with loading, with this complex relationship best described by a multiple regression equation with an interaction term between loading rate and loading time. In addition to informing interpretations of craniomandibular growth, form, function, and allometry, these comparisons provide a skeleton-wide perspective on the patterning of osteogenic stimuli across body sizes.

The daily habitual in vivo strain history of a non-weight-bearing bone

Daily mechanical loading strongly influences the architecture and composition of bone tissue. Throughout the day, the amplitudes, rates, frequencies, and the dispersion over time of these loads vary. Nevertheless, most experimental and descriptive studies on the aforementioned relationship consider only cyclic loading and, in addition, focus on weight-bearing bones. A more complete assessment of the daily loading of bone might lead to a better understanding of the natural everyday stimulus for bone maintenance or adaptive responses. In the present study, we measured the daily habitual strain history of the non-weight-bearing mandible bone in the rabbit. Long-term continuous strain recordings were made using an implantable telemetry device able to read out bone-bonded strain gauges. The lateral surface of the rabbit mandibular corpus was chosen as the bone surface of interest. During the recordings, which lasted up to 33 h, the rabbits (N=7) were able to move unrestrictedly in their cages, performing their habitual behaviours. Analysis of the recordings revealed that the measured bone surface was subjected to 2.9 (+/-1.4)x10(3) strain events per hour of which 1.8 (+/-1.0)x10(3) had amplitudes < or =10 microstrains (muvarepsilon). Larger strain amplitudes occurred less often and principal strains fell within the range of -517 (+/-118) muvarepsilon to 298 (+/-81) muvarepsilon. Strain rates never exceeded 10,000 muvarepsilon/s and only 8.9% (+/-7.2%) of the habitual strain rates were higher than 1000 muvarepsilon/s. Strain frequency spectra displayed clear peaks at 4-5 and 9 Hz. The wirelessly recorded daily strain history of the rabbit mandible featured peak strain amplitudes resembling those of other mammalian mandibles, but much smaller than those found in many long-bone strain measurements.