Biomechanical changes in the rabbit masticatory system during postnatal development

Lagomorpha as a Model Morphological System

Due to their global distribution, invasive history, and unique characteristics, European rabbits are recognizable almost anywhere on our planet. Although they are members of a much larger group of living and extinct mammals [Mammalia, Lagomorpha (rabbits, hares, and pikas)], the group is often characterized by several well-known genera (e.g., Oryctolagus, Sylvilagus, Lepus, and Ochotona). This representation does not capture the extraordinary diversity of behavior and form found throughout the order. Model organisms are commonly used as exemplars for biological research, but there are a limited number of model clades or lineages that have been used to study evolutionary morphology in a more explicitly comparative way. We present this review paper to show that lagomorphs are a strong system in which to study macro- and micro-scale patterns of morphological change within a clade that offers underappreciated levels of diversity. To this end, we offer a summary of the status of relevant aspects of lagomorph biology.

Dietary variation and mechanical properties of articular cartilage in the temporomandibular joint: implications for the role of plasticity in mechanobiology and pathobiology

Due to their nature as tissue composites, skeletal joints pose an additional challenge in terms of evaluating the functional significance of morphological variation in their bony and cartilaginous components in response to altered loading conditions. Arguably, this complexity requires more direct means of investigating joint plasticity and performance than typically employed to analyze macro- and micro-anatomical phenomena. To address a significant gap in our understanding of the plasticity of the mammalian temporomandibular joint (TMJ), we investigated the histology and mechanical properties of condylar articular cartilage in rabbits subjected to long-term variation in diet-induced masticatory stresses, specifically cyclical loading. Three cohorts of male weanlings were raised for six months on different diets until adulthood. Following euthanasia, the TMJ condyles on one side were dissected away, fixed, decalcified, dehydrated, embedded and sectioned. Safranin O staining was employed to identify variation in proteoglycan content, which in turn was used to predict patterns of articular cartilage stiffness in contralateral condylar specimens for each treatment group. Hematoxylin and eosin staining was used to quantify diet-induced changes in chondrocyte hypertrophy and cellularity. Mechanical tests document significant decreases in articular cartilage stiffness corresponding to patterns of extracellular matrix relative protein abundance in rabbits subjected to greater cyclical loading. This indicates that TMJs routinely subjected to higher masticatory stresses due to a challenging diet eventually develop postnatal decreases in the ability to counter compressive loads during postcanine biting and chewing. These findings provide novel information regarding TMJ performance, with broader implications about the costs and benefits of phenotypic plasticity as well as implications for how such biological processes affect connective tissue mechanobiology and pathobiology.

Shape Variation in the Craniomandibular System and Prevalence of Dental Problems in Domestic Rabbits: A Case Study in Evolutionary Veterinary Science

In contrast to wild lagomorphs, pet rabbits exhibit a noticeably high frequency of dental problems. Although dietary habits are considered as a major factor contributing to acquired malocclusions, the exact causes and interrelationships are still under debate. In this regard, an important aspect that has not been considered thoroughly to date is the effect of diet-induced phenotypic plasticity in skull morphology. Therefore, we conducted a geometric morphometric analysis on skull radiological images of wild and pet rabbits in order to quantify intraspecific variation in craniomandibular morphology. The statistical analyses reveal a significant morphological differentiation of the craniomandibular system between both groups. Furthermore, the analysis of covariance shows that the force-generating modules (cranium and mandible) vary independently from the force-receiving module (hypselodont teeth) in pet rabbits, which is in contrast to their wild relatives. Our findings suggest that the phenotypic changes in domestic rabbits impact mastication performance and, consequently, oral health. An adequate close-to-nature nutrition throughout the whole life and especially beginning early parallel to weaning (phase of increased phenotypic plasticity) is necessary to ensure a normal strain on the teeth by promoting physiological lateral gliding movements and avoiding direct axial loads.

Masticatory biomechanics in the rabbit: a multi-body dynamics analysis

Multi-body dynamics is a powerful engineering tool which is becoming increasingly popular for the simulation and analysis of skull biomechanics. This paper presents the first application of multi-body dynamics to analyse the biomechanics of the rabbit skull. A model has been constructed through the combination of manual dissection and three-dimensional imaging techniques (magnetic resonance imaging and micro-computed tomography). Individual muscles are represented with multiple layers, thus more accurately modelling muscle fibres with complex lines of action. Model validity was sought through comparing experimentally measured maximum incisor bite forces with those predicted by the model. Simulations of molar biting highlighted the ability of the masticatory system to alter recruitment of two muscle groups, in order to generate shearing or crushing movements. Molar shearing is capable of processing a food bolus in all three orthogonal directions, whereas molar crushing and incisor biting are predominately directed vertically. Simulations also show that the masticatory system is adapted to process foods through several cycles with low muscle activations, presumably in order to prevent rapidly fatiguing fast fibres during repeated chewing cycles. Our study demonstrates the usefulness of a validated multi-body dynamics model for investigating feeding biomechanics in the rabbit, and shows the potential for complementing and eventually reducing in vivo experiments.

Functional geometric morphometric analysis of masticatory system ontogeny in papionin primates

The three-dimensional configuration of the primate masticatory system is constrained by the need to maximize bite forces while avoiding distraction of the temporomandibular joint (TMJ). Within these bounds, shape variation has predictable effects on functional capacities such as mechanical advantage and gape. In this study, geometric morphometric analysis is used to investigate the ontogeny of masticatory function in papionin monkeys and test the hypothesis that biomechanical constraints determine the location of molar eruption. This "constrained eruption hypothesis" predicts that the distalmost molar (DMX) will occupy a consistent location anterior to the TMJ and that jaw adductor muscles will maintain consistent positions relative to both DMX and TMJ. Craniometric landmarks were digitized on cross-sectional ontogenetic series of nine papionin species. Form-space PCA of Procrustes residuals, visualization of Bookstein shape coordinates, and nonparametric ANOVA were used to identify ontogenetic shape trends and test for significant ontogenetic changes in relative landmark positions. In most taxa, DMX maintains a consistent position relative to the TMJ while the anterior dentition migrates anteriorly. Where significant intraspecific ontogenetic differences occur, they involve anterior migration of DMX in later dental stages, likely due to late adolescent growth of the posterior palate. Attachments of the anterior temporalis and deep masseter also maintain consistent positions relative to the TMJ; however, the superficial masseter migrates anteriorly throughout ontogeny. All muscle attachments migrate laterally relative to the TMJ, reflecting positive scaling of adductor PCSA. Overall, results support the constrained eruption hypothesis and suggest mechanisms by which functional capacity is maintained during ontogeny.

Masticatory loading, function, and plasticity: a microanatomical analysis of mammalian circumorbital soft-tissue structures

In contrast to experimental evidence regarding the postorbital bar, postorbital septum, and browridge, there is exceedingly little evidence regarding the load-bearing nature of soft-tissue structures of the mammalian circumorbital region. This hinders our understanding of pronounced transformations during primate origins, in which euprimates evolved a postorbital bar from an ancestor with the primitive mammalian condition where only soft tissues spanned the lateral orbital margin between frontal bone and zygomatic arch. To address this significant gap, we investigated the postorbital microanatomy of rabbits subjected to long-term variation in diet-induced masticatory stresses. Rabbits exhibit a masticatory complex and feeding behaviors similar to primates, yet retain a more primitive mammalian circumorbital region. Three cohorts were obtained as weanlings and raised on different diets until adult. Following euthanasia, postorbital soft tissues were dissected away, fixed, and decalcified. These soft tissues were divided into inferior, intermediate, and superior units and then dehydrated, embedded, and sectioned. H&E staining was used to characterize overall architecture. Collagen orientation and complexity were evaluated via picrosirius-red staining. Safranin-O identified proteoglycan content with additional immunostaining performed to assess Type-II collagen expression. Surprisingly, the ligament along the lateral orbital wall was composed of elastic fibrocartilage. A more degraded organization of collagen fibers in this postorbital fibrocartilage is correlated with increased masticatory forces due to a more fracture-resistant diet. Furthermore, the lack of marked changes in the extracellular composition of the lateral orbital wall related to tissue viscoelasticity suggests it is unlikely that long-term exposure to elevated masticatory stresses underlies the development of a bony postorbital bar.

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.

Phenotypic plasticity and function of the hard palate in growing rabbits

Morphological variation related to differential loading is well known for many craniomandibular elements. Yet, the function of the hard palate, and in particular the manner in which cortical and trabecular bone of the palate respond to masticatory loads, remains more ambiguous. Here, experimental data are presented that address the naturalistic influence of biomechanical loading on the postweaning development and structure of the hard palate. A rabbit model was used to test the hypothesis that variation in the morphology of the hard palate is linked to variation in masticatory stresses. Rabbit siblings were divided as weanlings into soft and hard/tough dietary treatment groups of 10 subjects each and were raised for 15 weeks until subadulthood. MicroCT analyses indicate that rabbits subjected to elevated masticatory loading developed hard palates with significantly greater bone area, greater cortical bone thickness along the oral lamina, and thicker anterior palates. Such diet-induced levels of palatal plasticity are comparable to those for other masticatory elements, which likely reflect osteogenic responses for maintaining the functional integrity of the palate vis-à-vis elevated stresses during unilateral mastication. These data support a role for mechanical loading in the determination of palatal morphology, especially its internal structure, in living and fossil mammals such as the hominin Paranthropus. Furthermore, these findings have potential implications for the evolution of the mammalian secondary hard palate as well as for clinical considerations of human oral pathologies.

Pushing the limit: masticatory stress and adaptive plasticity in mammalian craniomandibular joints

Excessive, repetitive and altered loading have been implicated in the initiation of a series of soft- and hard-tissue responses or ;functional adaptations' of masticatory and locomotor elements. Such adaptive plasticity in tissue types appears designed to maintain a sufficient safety factor, and thus the integrity of given element or system, for a predominant loading environment(s). Employing a mammalian species for which considerable in vivo data on masticatory behaviors are available, genetically similar domestic white rabbits were raised on diets of different mechanical properties so as to develop an experimental model of joint function in a normal range of physiological loads. These integrative experiments are used to unravel the dynamic inter-relationships among mechanical loading, tissue adaptive plasticity, norms of reaction and performance in two cranial joint systems: the mandibular symphysis and temporomandibular joint (TMJ). Here, we argue that a critical component of current and future research on adaptive plasticity in the skull, and especially cranial joints, should employ a multifaceted characterization of a functional system, one that incorporates data on myriad tissues so as to evaluate the role of altered load versus differential tissue response on the anatomical, cellular and molecular processes that contribute to the strength of such composite structures. Our study also suggests that the short-term duration of earlier analyses of cranial joint tissues may offer a limited notion of the complex process of developmental plasticity, especially as it relates to the effects of long-term variation in mechanical loads, when a joint is increasingly characterized by adaptive and degradative changes in tissue structure and composition. Indeed, it is likely that a component of the adaptive increases in rabbit TMJ and symphyseal proportions and biomineralization represent a compensatory mechanism to cartilage degradation that serves to maintain the overall functional integrity of each joint system. Therefore, while variation in cranial joint anatomy and performance among sister taxa is, in part, an epiphenomenon of interspecific differences in diet-induced masticatory stresses characterizing the individual ontogenies of the members of a species, this behavioral signal may be increasingly mitigated in over-loaded and perhaps older organisms by the interplay between adaptive and degradative tissue responses.

Heterogeneous postnatal transitions in myosin heavy chain isoforms within the rabbit temporalis muscle

Postnatal changes in the fiber type composition and fiber cross-sectional area were investigated in the superficial (TEM1) and deep (TEM23) temporalis of male rabbits. It was hypothesized that, due to the transition from suckling to chewing during early postnatal development, the proportion of fast fiber types would decrease, while the proportion of fibers positive for myosin heavy chain (MyHC) cardiac alpha would increase, and that, due to the influence of testosterone during late postnatal development, the proportion of these alpha fibers would decrease again. Classification of the fibers types was performed by immunohistochemistry according to their MyHC content. The proportion of alpha fiber types significantly increased in both muscle portions from 2% and 8% for TEM1 and TEM23 at week 1 to 29% and 54% at week 8, respectively,. While in TEM1 the proportion of this fiber type did not change thereafter, it decreased again to 27% in TEM23 at week 20. The change for the fast fiber types was opposite to that of the alpha fiber types. Significantly more MyHC IIX fibers were found in TEM1 than in TEM23 in adult rabbits. In the first 8 weeks, the cross-sectional areas of all fibers increased. After this period, only MyHC cardiac alpha + I fibers continued to increase significantly. It was concluded that there are developmental differences in the myosin heavy chain transitions of the two portions of the temporalis muscle.

Postnatal transitions in myosin heavy chain isoforms of the rabbit superficial masseter and digastric muscle

We investigated the early (< 8 weeks) and late (> 8 weeks) postnatal development of the fibre type composition and fibre cross-sectional area in the superficial masseter and digastric muscle of male rabbits. It was hypothesized, first, that due to the transition between suckling and chewing, during early postnatal development the increase in the proportion of slow fibre types and in fibre cross-sectional areas would be larger in the masseter than in the digastric; and second, that due to the supposed influence of testosterone during late postnatal development, the proportion of slow fibre types in both muscles would decrease. Fibre types were classified by immunostaining according to their myosin heavy chain (MyHC) content. The proportion of slow fibre types significantly increased in the masseter, from 7% at week 1 to 47% at week 8, and then decreased to 21% at week 20, while in the digastric it increased from 5% in week 1 to 19% at week 8 and remained the same thereafter. The changes in the proportion of fast fibre types were the opposite. The remarkable increase and decrease in the proportion of slow fibre types in the masseter was attributed predominantly to MyHC-cardiac alpha fibres. During early development, the cross-sectional area of all fibres in both muscles increased. However, only the fast fibre types in the masseter continued to grow further after week 8. Before weaning, the fast fibre types in the digastric were larger than those in the masseter, but after week 8, they became larger in the masseter than in the digastric. In adult animals, masseter and digastric had the same percentage of fast fibre types, but these fibres were almost twice as large in masseter as in digastric.

Daily activity of the rabbit jaw muscles during early postnatal development

Early postnatal development of the jaw muscles is characterized by the transition from suckling to chewing behavior. As chewing develops the jaw closing muscles become more powerful compared with the jaw openers. These changes are likely to affect the amount of daily muscle activity. Therefore, the purpose of this study was to characterize for a jaw opener (digastric) and jaw closer (masseter) the total duration of daily muscle activity (i.e. the duty time), and the daily burst numbers and lengths during early postnatal development. Using radiotelemetry the activity of these muscles was recorded in 10 young New Zealand White rabbits between three and eight weeks of age. Fiber-type composition was analyzed at eight weeks of age by determining the myosin heavy chain content of the fibers. During postnatal development both muscles showed no significant decrease or increase in their daily activity. However, the interindividual variation of the duty time and burst number significantly decreased. There were no significant differences between the digastric and masseter except for the most powerful activities at eight weeks of age, where the masseter showed a significantly higher duty time and burst number than the digastric. The masseter contained a higher number of slow-type fibers expressing myosin heavy chain-I and myosin heavy chain-cardiac alpha than the digastric. The present results suggest that the amount of jaw muscle activation is already established early during postnatal development, before the transition from suckling to chewing behavior. This amount of activation seems to be related to the number of slow-type fibers.

Mechanical strain on the human skull in a humanoid robotic model

Patterns of strain were analyzed in a dry human skull at 15 different regions on the lateral and medial surfaces of the mandible. The strains were induced with a human robotic system that represented each of 8 bilateral muscles by a DC servomotor connected to a wire and pulley. The tractions of the simulated muscles (masseter, medial pterygoid, anterior temporalis, and posterior temporalis) were increased from 1x to 4x with each representing different levels of traction or force (5, 3, 4, and 4 N, respectively). The study was done with the teeth in maximum intercuspal occlusion. Bite forces were also measured with a transducer and reached a maximum of 40 N on the posterior teeth with less force on the anterior dentition. The smallest traction level (1x) developed some small strains. At 2x, compressive strains developed more on the medial (lingual) side beneath the molars through the corpus and radiated into the anterior ramus. Strains at 3x to 4x significantly increased both the tensile and compressive strains throughout the mandible with more strains developing in the ramus. The increased bilateral traction and loading developed significant compressive forces on both sides of the mandible. Evaluation of disparities between compressive and tensile strains at one site, and comparison between the medial and lateral sides of strain, suggested some visible distortion of portions of the mandible under the higher loads.

Synchronous recording of load-deformation behaviour and polarized light-microscopic images of the rabbit incisor periodontal ligament during tensile loading

Tooth-periodontal ligament-bone segments were cut in the form of rectangular prisms (1.5 mm wide, 0.65 mm thick, and long enough to allow anchorage of the bone and tooth-end portions in a stretching jig) from the mandibular incisors of 10 rabbits. The experimental set-up enabled simultaneous recording, on video, of the changing image brightness under polarizing optics together with extension across the periodontal ligament. Specimens were stretched until failure at a velocity of 0.5 mm/min. The tensile load-deformation curve of the ligament exhibited an initial, non-linear region that was followed by a linear region, a subsequent yielding region preceding the maximum point, and a final descending region. Gradual increases in the intensity of birefringence in the linear and yielding regions indicated that stress concentrations occur in the supporting fibres attached to mineralized tissues. In the final descending region of the curve, progressive breakages of individual fibre bundles occurred, mainly in the middle zone of the ligament. Analysis of the polarized light-microscopic images showed that the increases in brightness and area of birefringent collagen fibre bundles occurred in parallel with the stress generated. These results suggest that the collagen fibre bundles became aligned with the direction of loading and the intensity of their birefringence increased according to the applied tensile force.

The role of passive muscle tensions in a three-dimensional dynamic model of the human jaw

The role of passive muscle tensions in human jaw function are largely unknown. It seems reasonable to assume that passive muscle-tension properties are optimized for the multiple physiological tasks the jaw performs in vivo. However, the inaccessibility of the jaw muscles is a major obstacle to measuring their passive tensions, and understanding their effects. Computer modelling offers an alternative method for doing this. Here, a three-dimensional, dynamic model was used to predict active and passive jaw-muscle tensions during simulated postural rest, jaw opening and chewing. The model included a rigid mandible, two temporomandibular joints, multiple dental bite points, and an artificial food bolus located between the right first molars. It was driven by 18 Hill-type actuators representing nine pairs of jaw muscles. All anatomical forms, positions and properties used in the model were based on previously published, average values. Two states were stimulated, one in which all optimal lengths for the length-tension curves in the closing muscles were defined as their fibre-component lengths when the incisor teeth were 2 mm apart (S2), and another in which the optimal lengths were set for a 12.0 mm interincisal separation (S12). At rest, the jaw attained 3.6 mm interincisal separation in S2, and 14.8 mm in S12. Activation of the inferior lateral pterygoid (ILP) and digastric (DG) muscles in various combinations always induced passive jaw-closer tensions, and compressive condylar loads. Maximum midline gape (from maximum bilateral co-activation of DG and ILP) was 16.2 mm in S2, and 32.0 mm in S12. When both model states were driven with muscle patterns typical for human mastication, recognizable unilateral and vertical "chopping" chewing cycles were produced. Both states revealed condylar loading in the opening and closing phases of mastication. During unilateral chewing, compressive force on the working-side condyle exceeded that on the balancing side. In contrast, during the "chopping" cycle, loading on the balancing side was greater than that on the working side. In S2, chewing was limited in both vertical and lateral directions. These results suggest that the assumptions used in S12 more closely approximated human behaviour than those in S2. Despite its limitations, modelling appears to provide a useful conceptual framework for developing hypotheses regarding the role of muscle tensions during human jaw function.

Effect of maturation on 31P magnetic resonance spectroscopy of the rabbit masseter muscle

This work studies the dynamic metabolic changes of the rabbit masseter muscle during post-natal development. The composition and proportion of oxidative and glycolytic muscle fibers alter during maturation. The masseter muscle, as most muscles of the craniofacial region, exhibits unusual development in composition of isoforms of myosin. The effect of this unusual composition on the dynamic metabolic properties of the masseter muscle have not been assessed. The metabolism of the rabbit masseter muscle was studied by means of 31P-nuclear magnetic resonance (NMR) spectroscopy. Contraction was elicited by electrical stimulation of the muscle in the anesthetized animal. Five animals were studied at 8 weeks and 24 weeks so that both the juvenile and adult stages could be evaluated. The dynamic biochemical changes in the masseter muscle were studied by the analysis of NMR spectra. A single-turn surface coil (copper) was used, and the original signal was treated with Fourier transforms to obtain 31P spectra. The low signal-to-noise ratio required averaging 16 acquisitions (acquisition time = 400 msec, repetition rate = 1.8 sec) in 30 sec and then obtaining continuous spectra for 27 min. Each averaged spectrum demonstrated five peaks: inorganic phosphate (Pi), creatine phosphate (PCr), and three peaks related to adenosine triphosphate (ATP). The protocol involved recording an initial three-minute rest period, stimulating the muscle at 5 Hz for 3 min twice, separated by three-minute rest periods, and stimulating the muscle at 50 Hz twice for 3 min separated by rest periods. The Pi/PCr ratio increased significantly in the adult masseter during both 5-Hz stimulations, evoking twitching, and the first 50-Hz stimulation, evoking tetany (repeated ANOVA, P < 0.05). The resting pH (6.96 +/- 0.13) was significantly lowered during both twitching (6.85 +/- 0.10; P < 0.0038) and tetany (6.55 +/- 0.13; P < 0.0001), but only in the adult masseter muscle. These finding suggest that the adult masseter muscle possesses more glycolytic fibers as it modifies its metabolism during postnatal development.

Plasticity of craniomandibular muscle function: 31P magnetic resonance spectroscopy of the rabbit masseter muscle

The masseter muscle was studied during postnatal development of the rabbit from the juvenile to adult stage in which the oral function was altered during maturation by modifying the diet to soft food. The muscle was assessed using phosphate magnetic resonance (31P NMR) spectroscopy with a single-turn copper surface coil to study potential changes in phosphate metabolism. The 31P NMR spectra consisted of five peaks related to unbound forms of inorganic phosphate (Pi), creatine phosphate (PCr), and three peaks related to the adenosine triphosphate (ATP). The masseter was assessed in one group of five rabbits at 8 weeks postnatally (juvenile) and after 4 months of this experimental group masticating on soft food. They were compared with a control group of five rabbits raised on a normal hard diet. The Pi/PCr ratio increased in the adult masseter much higher during twitching, tetany, and periodic contraction than in the juvenile regardless as to whether the adult animal had been raised from the juvenile period on soft or hard diet. There were relatively few differences between the experimental adult animals raised on a soft diet and the normal adult animals despite the soft diet animals demonstrating a significantly lower weight and smaller muscle mass. These findings suggest that chronic underuse of the masseter muscle by decreasing the masticatory loads has a minimal effect on the phosphate metabolism of the maturing masseter.

31P-magnetic resonance spectroscopy of the rabbit masseter muscle

Dynamic biochemical changes in the masseter muscle were studied in 14 New Zealand adult male rabbits by 31P-nuclear magnetic resonance (NMR) spectroscopy. NMR spectra were obtained during rest and electrical stimulation of the muscle in the anaesthetized animal at 33 recording sessions. Electrical stimulation was applied by a pair of copper wires placed separately with hypodermic needles into the muscle. NMR spectra were acquired with a 2 x 3 cm, double-turn, copper transmit/receive coil. Sixteen spectra were averaged over 30 s to obtain averaged spectra continuously during a 30-min recording. The spectra were processed automatically using a non-linear 'least-squares' fitting program on the spectrometer. A Lorentzian line shape was assumed for the peaks, and values of peak height, area and chemical shifts were generated. Each averaged spectrum consisted of five peaks: inorganic phosphate (Pi), creatine phosphate (PCr), and three peaks related to ATP. Data were analysed as to absolute changes in Pi and PCr, in the ratio of Pi/PCr, and the shift of Pi to PCr to estimate pH. Several protocols were used in which ranges of frequency, intensity and duration of electrical stimulation were tested. The protocol for detailed studies involved stimulating the muscle twice at 5 Hz for 3 min separated by a 3-min rest period, then stimulating twice at 50 Hz for 3 min separated by a rest period. During contraction of the muscle, there was a significant increase in the Pi/PCr ratio (p < 0.05) as compared to the resting level. The ratio reached a plateau over a 3-min contraction using 5-Hz stimulation, then increased significantly more with the 50-Hz stimulation but decayed during the 3 min. Sustained stimulation with 50 Hz for 15-45 min evoked an initial sharp change in Pi/PCr, which then reached a steady plateau that remained over the entire stimulation. These findings indicate that the rabbit masseter muscle is relatively fatigue resistant in maintaining a steady-state equilibrium in the relation of Pi to PCr.