Arie Van Der Meijden

Julio César González-Gómez

María D Pulido-Osorio

Anthony Herrel

CIBIO/InBio Biopolis

JEB -Methods and Techniques

Keywords: Automated voluntary bite force measurement Bite force, carnivora, Environmental enrichment, feeding

Bite force is a key performance trait of the feeding system, but maximal in vivo bite force has been measured in few large mammals. The alternative, modelling of bite force from anatomy, cannot be validated without in vivo measurements. To overcome existing limitations of ethics, safety, and animal well-being, we here propose a semi-automated method to obtain voluntary maximum bite forces from large mammals using bite plates that automatically dispense a food reward if an incrementally increasing threshold force value is reached. We validated our method using two Malayan sun bears, two Andean spectacled bears and a lioness. We show that voluntary bite force measurement using positive reinforcement is a non-invasive and reliable method to record maximum voluntary bite force performance in large mammals. Our results further show that in vivo data are critical as modeling efforts from osteology have greatly underestimated bite forces in Andean bears.

Introduction

Feeding is among the most essential interactions of an animal with its environment [START_REF] Schwenk | Feeding: form, function and evolution in tetrapod vertebrates[END_REF]. Specializations of the feeding apparatus allow specific trophic niches to be exploited, often at the expense of taking advantage of others due to biomechanical trade-offs [START_REF] Aguirre | Ecomorphological analysis of trophic niche partitioning in a tropical savannah bat community[END_REF][START_REF] Herrel | Force-velocity trade-off in Darwin's finch jaw function: A biomechanical basis for ecological speciation?[END_REF][START_REF] Konuma | Trade-offs between force and fit: Extreme morphologies associated with feeding behavior in carabid beetles[END_REF][START_REF] Van Valkenburgh | Déjà vu: The evolution of feeding morphologies in the Carnivora[END_REF]. The biomechanics and comparative functional ecology of the feeding apparatus has therefore been studied extensively to understand the relative roles of adaptation and functional trade-offs in driving the diversity of the feeding system [START_REF] Aguirre | Ecomorphological analysis of trophic niche partitioning in a tropical savannah bat community[END_REF][START_REF] Christiansen | Bite forces and evolutionary adaptations to feeding ecology in carnivores[END_REF][START_REF] Grossnickle | Feeding ecology has a stronger evolutionary influence on functional morphology than on body mass in mammals[END_REF][START_REF] Herrel | Evolution of bite performance in turtles[END_REF][START_REF] Wroe | Convergence and remarkably consistent constraint in the evolution of carnivore skull shape[END_REF][START_REF] Wroe | Bite club: comparative bite force in big biting mammals and the prediction of predatory behaviour in fossil taxa[END_REF]. The performance of the feeding apparatus in vertebrates includes a variety of traits such as processing rate, gape, and bite force. The latter has been shown to be a limiting factor in the feeding performance of many vertebrates [START_REF] Aguirre | The implications of food hardness for diet in bats[END_REF][START_REF] Maestri | Diet, bite force and skull morphology in the generalist rodent morphotype[END_REF]. Beyond feeding and prey-capture, the jaw-cranial system is an essential functional complex in defensive and reproductive behavior. Quantifying the performance of this complex is therefore necessary to understand the ecological role and evolution of the whole organism [START_REF] Anderson | Bite force in vertebrates: opportunities and caveats for use of a nonpareil whole-animal performance measure[END_REF], and can be informative in terms of conservation biology [START_REF] Erickson | Comparison of bite-force performance between long-term captive and wild American alligators ( Alligator mississippiensis )[END_REF][START_REF] Fleming | Body Size and Bite Force of Stray and Feral Cats-Are Bigger or Older Cats Taking the Largest or More Difficult-to-Handle Prey[END_REF].

Maximum bite force is a key performance metric of the jaw-cranial system, determining dietary diversity and defensive ability in a wide range of taxa. Despite its importance, bite force has been measured in surprisingly few large vertebrates with the exception of crocodiles [START_REF] Dumont | The effects of gape angle and bite point on bite force in bats[END_REF][START_REF] Erickson | The ontogeny of bite-force performance in American alligator (Alligator mississippiensis)[END_REF]. Large, and sometimes dangerous animals are notoriously difficult to measure the bite force of, likely due to limited access and the inability to manipulate these animals easily. This is particularly true for large mammals. There are understandable limitations to the handling and coercion necessary to produce strong defensive bites in these animals, as such provoked defensive bites involve considerable risk to the experimenter and the animal. Furthermore, it may negatively affect the association the animal has with humans, sometimes permanently. Thus, considerations of animal well-being, ethics, conservation, and safety have severely limited the availability of larger mammals for the measurement of bite forces.

To fill this hiatus many researchers have used a modeling approach, yet these models cannot be validated if no actual bite force measurements exist [START_REF] Christiansen | Evolutionary implications of bite mechanics and feeding ecology in bears[END_REF][START_REF] Heethoff | A new use for synchrotron X-ray microtomography: threedimensional biomechanical modeling of chelicerate mouthparts and calculation of theoretical bite forces[END_REF][START_REF] Thomason | Cranial strength in relation to estimated biting forces in some mammals[END_REF][START_REF] Westneat | A biomechanical model for analysis of muscle force, power output and lower jaw motion in fishes[END_REF][START_REF] Wroe | Bite club: comparative bite force in big biting mammals and the prediction of predatory behaviour in fossil taxa[END_REF]. Any ecologically relevant comparative study of the feeding apparatus of larger mammals therefore requires measurement of bite force in vivo to get a direct estimate of maximum bite force to validate existing models. Once the principles are better understood and modelled in extant species, inferences can be then made on the performance and trophic niche of fossil species, allowing better insights into the evolution of the ecomorphological relationships in emblematic groups such as large carnivores. To measure in vivo bite force, several methods have thus far been employed. Bite force can be measured by either making an animal bite a device, and recording the force, or by implanting a sensor in the teeth. The latter method requires anesthesia, dental surgery and tethering the animals with wires, and is highly invasive. It has been used to measure the bite force of domestic dogs [START_REF] Brunski | In vivo forces on dental implants: Hard-wiring and telemetry methods[END_REF], domestic pigs [START_REF] Bousdras | A novel approach to bite force measurements in a porcine model in vivo[END_REF] and macaques [START_REF] Hylander | An in vivo strain-gauge analysis of the squamosal-dentary joint reaction force during mastication and incisal biting in Macaca mulatta and Macaca fascicularis[END_REF]. Less invasive is to induce the animal to bite on a device, and recording the force. The bite on the device can be enticed behaviorally or through electrostimulation of the jaw muscles. The latter method has been used on various anesthetized animals: domestic dogs [START_REF] Ellis | Calibration of estimated biting forces in domestic canids: comparison of post-mortem and in vivo measurements[END_REF], domestic pigs [START_REF] Ström | Bite-force development, metabolic and circulatory response to electrical stimulation in the canine and porcine masseter muscles[END_REF], possums [START_REF] Thomason | Forces of biting, body size, and masticatory muscle tension in the opossum Didelphis virginiana[END_REF] and rhesus monkeys [START_REF] Dechow | A method of bite force measurement in primates[END_REF].

Animals can be behaviorally enticed to bite a force sensing device through defensive behavior, feeding, play, or by training. Defensive bites are the easiest to elicit, particularly in wild animals.

This method has been used on mammals ranging from bats [START_REF] Aguirre | The implications of food hardness for diet in bats[END_REF][START_REF] Dumont | The effects of gape angle and bite point on bite force in bats[END_REF] and rodents [START_REF] Becerra | Another one bites the dust: Bite force and ecology in three caviomorph rodents (Rodentia, Hystricognathi)[END_REF] to ferrets [START_REF] Dessem | Jaw-muscle activity in ferrets,Mustela putorius furo[END_REF], wild pigs [START_REF] Sicuro | Quantifying bite force in coexisting tayassuids and feral suids: a comparison between morphometric functional proxies and in vivo measurements[END_REF] and canids including foxes and dogs [START_REF] Brassard | Bite force and its relation to jaw shape in domestic dogs[END_REF][START_REF] Brassard | Masticatory system integration in a commensal canid: interrelationships between bones, muscles and bite force in the red fox[END_REF]. Utilizing feeding and/or play behavior, [START_REF] Lindner | Measurement of Bite Force in Dogs: A Pilot Study[END_REF] measured bite forces of domestic dogs with a rawhide covered force sensor, but could only elicit 6-9 bites of even the most enthusiastic participating dogs. A combination of training and play was used to measure bite forces in spotted hyenas [START_REF] Binder | Development of bite strength and feeding behaviour in juvenile spotted hyenas ( Crocuta crocuta )[END_REF] and manual training was also used to measure the bite forces of dogs [START_REF] Brassard | Bite force and its relation to jaw shape in domestic dogs[END_REF]. Manual training with a food reward was used to train restrained macaques and galagos to bite a force sensor [START_REF] Hylander | Mandibular function in Galago crassicaudatus and Macaca fascicularis: An in vivo approach to Stress Analysis of the mandible[END_REF]. However, most of the above-mentioned studies were invasive or time-consuming, and bites were part of an interaction with humans in non-feeding contexts. To overcome limitations of ethics, safety, and animal well-being, we here propose an automated method to obtain voluntary maximum bite forces from large mammals, and validate it as part of environmental enrichment in a zoo setting.

We validated our method using two bear species and a lion, and show that in vivo data are critical as current modeling efforts may underestimate bite forces.

Materials and Methods

Description of the bite force device

To record bite forces in bears we designed sturdy bite plates that have less than 2.5 mm bending deflection under the maximum expected loads. In order to reduce bending by increasing the second moment of area of the bite plate beams, the plates were attached to stainless steel ubeams. These were further stiffened, where possible, with square steel tubing (see Figure 1). The force is transmitted from the bite plates to a piezo sensor with a range of 20 kN (Kistler model 9331B). In voluntary biting, the position of the teeth on the plates cannot be precisely controlled [START_REF] Lappin | Reliable quantification of bite-force performance requires use of appropriate biting substrate and standardization of bite out-lever[END_REF]. This may result in variation in the effective in-lever length of the plates.

To minimize this source of variation, we exposed only the first 45 mm of the plates to the animals, and lengthened the in-lever arms to 435 mm, resulting in a mechanical advantage of 1.37 in the middle of the plates ± 0.083 to the front and rear extremes. This design has the added advantage of lowering the fulcrum reaction forces. See supplementary file 1 for details of the parts.

The charge signal of the piezo sensor was transformed to a voltage using a Kistler 5018a charge amplifier, which was then sampled using a LabJack U3 I/O board at a frequency of approximately 2 kHz. This high sampling rate is necessary, as we did not know the loading and unloading rates and bite durations of bear bites. Under-sampling may unnecessarily increase the variance in measured bite forces, and can lead to an underestimation of the bite force. The samples were converted to Newtons, corrected for the mechanical advantage of the middle of the plates, and recorded using a laptop computer.

The exposed bite plates were sanded to reduce sharp corners. Then, to protect the teeth of the animals from the stainless-steel plates, we covered the plates with three layers of "Aramid Armor" Kevlar mesh tubing (bears) or two layers of fire hose material (lion). We chose these materials for their high tensile strength to withstand pulling, and their high wear resistance. We found the latter to be more resistant to attempts at destruction.

To provide a reward for each bite, we placed a thin flexible plastic hose between the bite plates through which honey, diluted with lukewarm water for fluidity, could be delivered directly into the mouth of the bears. The lioness was rewarded with cow blood. The delivery of the liquid reward was controlled with a solenoid valve, which could be controlled from the laptop through the output of the LabJack I/O board, and by hand using a remote-control button. When the solenoid valve was activated, the liquid reward was forced through the tube by gravity, as the reward container was placed higher above the ground than the bite plates.

To record which side of the jaw was used, how the jaws were applied to the plates, and which individual was doing the biting, a video of each bite with a force above the trigger threshold was recorded using a USB camera (ELP, Ailipu, Shenzhen, China, see figure 2).

To achieve maximal bite forces, we programmed the reward to be triggered above an incrementally increasing threshold, initially set at 20N. The peak bite force was then stored only if it was above the threshold. The reward threshold was then increased to be the mean of the previous bites, with a maximum of 20 previous bites. In this manner, the reward threshold increased a little with each subsequent bite above the previous threshold, as the mean of a sliding window of the last 20 measurements. This was intended to mimic the natural slow depletion of a hard object with a desired food, such as chewing on a piece of wood with a honeycomb, or on an animal bone.

Testing and validation of the method

To test, improve and validate our method, we introduced the setup to captive specimens of two bear species: the Malayan sun bear (Helarctos malayanus (Raffles, 1821)) and the Andean spectacled bear (Tremarctos ornatus (Cuvier, 1825)). To show the method also works for felids, we tested the device a lioness (Panthera leo (Linnaeus, 1758)).

We started with two female specimens of the Malayan sun bear named Linh (born 8-4-2011) and Cindy (born 16-7-2000) at the Ouwehands Dierenpark Rhenen in the Netherlands. Both bears were born in captivity, and are estimated to weigh about 50kg. Initially, the bears were separated in order to avoid any competition. The device was attached to the bars of the indoor part of the enclosure, with only the bite plates protruding beyond the bars (see supplementary file 1). The bears were introduced to the device by their daily caretakers. Interest in the device was rewarded with peanuts and walnuts as external rewards. Bear Linh rapidly learned to bite the plates, and was rewarded both with the diluted honey automatically dispensed from the device, as well as the external rewards. After this initial training, the bears were left with the device available in their enclosure without supervision. Linh regularly used the biteforce plates correctly to obtain the reward over a period of several days, but Cindy had to be re-trained multiple times. This allowed us to test different training approaches, including using dummy plates and external rewards. We found that training by rewarding desired behavior using manual control of the internal reward system only, was most effective. This training method was therefore used subsequently.

At the DierenPark Amersfoort in the Netherlands we measured bite forces of two young male specimens of the Andean spectacled bear, two brothers named Ariba (91kg) and Teju (80kg, both born in captivity 24-12-2017). Both bears were trained exclusively using manual control of the internal reward system, and each learned to bite the plates correctly, using only the molars, within a single session of less than 20 minutes. The bears were then allowed access to the device in their enclosure for several hours per day in the subsequent days, with the reward being triggered automatically. Re-training with manual control of the internal reward system was done only when necessary. Lioness Savannah at Maia Zoo in Portugal was trained over 3-4 sessions of 20-30 minutes. Savannah is estimated to weigh 100-120kg and was born 13-5-2014 in captivity.

Since the animals can choose how they reach forces high enough to trigger a reward, there was a potential for them to use their body mass to achieve a higher downward force on the upper plate.

To check whether the recorded bite forces were unintentionally inflated in this manner, we mounted the device upside down in the last session with the spectacled bears. In this manner, the lower plate is recording the force, and any body mass force exerted on the upper plate is not transferred to the transducer.

Analysis of bite force data

For each bite event, the force data, the time, and the accompanying video were recorded on the laptop computer using a custom Python script (available from the lead author upon request). After each day of trials, the videos were reviewed to identify the bear that was biting, and bites were classified as valid (unilateral bites using molars) or invalid (anything else). For valid bites, we noted which side of the jaw the bear used. This allowed us to test for fatigue and left-right preferences.

The lioness only used the canines or incisors to bite the device.

Before analysis, the bite force traces were filtered using a 4 th order low-pass Butterworth filter with a cutoff frequency of 20Hz. To quantify the duration of the bites, we used the find_peaks function in the Scripy python package to determine the peak widths with a relative height setting of 0.1 and a minimum peak height of 100N.

To analyze the bear data, we performed linear models and ANOVAs on the log 10 -transformed data in R version 21.09.1 (R Development Core Team 3.0.1., 2013) using the lmer function of the lme4 package [START_REF] Kuznetsova | lmerTest Package: Tests in Linear Mixed Effects Models[END_REF]. We performed an ANOVA to test for differences in force using To compare our in vivo measurements at the rear molar with the estimated values of bite force at the carnassials and the canines by [START_REF] Christiansen | Evolutionary implications of bite mechanics and feeding ecology in bears[END_REF], we measured the distances of the molar, carnassial and canine tip to the joint axis from 3D scans of a skull of T. ornatus (MNHN specimen CG1990-696_F) and H. malayanus (MNHN specimen GC1914-360). To calculate the force estimated at the rear molar, we multiplied the force at the carnassial by the ratio of the distances from the joint axis to the carnassial and the rear molar, and similarly with the ration of out-lever distances of the canine and rear molar. Since the lioness used only the canines and incisors to bite, these values could be compared to modeled values directly.

Results

We recorded 4,068 valid bites in four recording sessions/days for the Andean bears, 1,937 bites in five recording sessions/days for the Malayan bears, and 852 bites over seven recording days for the lioness (see Table 1).

Behavioral observations

We found that the animals approached and used the device repeatedly, despite being fed their normal diet and having access to their complete indoor and outdoor enclosures. Wear marks and video recordings indicate that most bites were in the middle of the plates. Bears quickly favored unilateral bites with the rear molars, while the lioness settled on mostly using the incisors. Bites were usually concentrated in bouts of many repeated bites until the force necessary to trigger a reward became very high. We noted large individual differences in the necessity for re-training between the individuals in both species. In the Malayan sun bears, Linh was able to keep using the device correctly over a period of several days. Bear Cindy however, had to be re-trained several times to use the device correctly. Similarly, in the Andean spectacled bears, Teju consistently used the device correctly, while Ariba had to be re-trained after a hiatus of a day or two. In both bear species, the bear that retained the training best was also the bear that most frequently used the device. Lioness Savannah was slower to learn to use the device, but did not require re-training.

Statistical analysis

There was no difference in force between left and right in the bears, even when taking species and specimen into account (species F= 0.0040, p=0.96, side F=0.026, p= 0.87).

We also found no evidence for fatigue. The relationship between bite force and time since the previous bite was significant (F=153, p<2e-16), but negative, suggesting that the higher forces were recorded as part of a burst of repeated bites (as in Figure 2).

We found no evidence for influence of the bite plate orientation, suggesting that bears do not use their body mass to increase downward force (F=0.87, p=0.35). In fact, the single highest bite force (2196.8N) was recorded with the plates inverted.

Across all bear specimens, we found that of the valid bites, 13.6% were over 80% of the maximal recorded bite force, reducing to 2.8% and 0.47% for bites over 90% and 95% of maximal bite force respectively (but note the large individual variation in Table 1). The lioness had a lower percentage of peak forces. For the bears, at an average rate of one valid bite every 3.3 seconds during a biting bout (here defined as a string of bites with a hiatus between bites of less than 100s), a bite > 95%

of the maximum force occurred on average every 12 minutes, and a bite of >90% of the maximum on average every 2 minutes.

Correspondence to models

We found that the recorded maximal bite forces at the rear molar for the Malayan sun bears was close to the estimated maximal force of 1722±423N, based on the model of [START_REF] Christiansen | Evolutionary implications of bite mechanics and feeding ecology in bears[END_REF].

The canine bite force of lions is estimated at 1104-1483N, based on the dry skull method and regression analysis respectively [START_REF] Thomason | Cranial strength in relation to estimated biting forces in some mammals[END_REF], and at 1768N by [START_REF] Wroe | Bite club: comparative bite force in big biting mammals and the prediction of predatory behaviour in fossil taxa[END_REF], which is in line with the maximal incisor bite force of 1594N we recorded. The recorded maximum in vivo force for the Andean spectacled bears of 2197N was much higher than the estimated force for this species at 1236±296N, even when corrected for the better mechanical advantage of the rear molar.

Discussion

Bite force is one of the key performance variables of the feeding system and may constrain dietary diversity, prey size and type, and the capacity of an animal to defend itself in intra-or interspecific contexts. Understanding bite force variation and how it relates to diet may be particularly important in rare or endangered species, especially species involved in head-start programs or rereleased into the wild. Our method seems to have overcome the ethical and safety limitations that have hitherto largely precluded the measurement of bite forces in large captive carnivores, many of which are vulnerable.

We found that the maximal bite forces of the Andean bears and the lioness we recorded are well above those predicted for these species based on the dry skull method [START_REF] Christiansen | Evolutionary implications of bite mechanics and feeding ecology in bears[END_REF]. This underestimation is a known shortcoming of the dry skull method, in part due to an overestimation of the physiological cross section area (PCSA) of the masseter and underestimation of the PCSA of the temporalis muscle [START_REF] Law | Dry versus wet and gross: Comparisons between the dry skull method and gross dissection in estimations of jaw muscle cross-sectional area and bite forces in sea otters[END_REF][START_REF] Thomason | Cranial strength in relation to estimated biting forces in some mammals[END_REF][START_REF] Wroe | Bite club: comparative bite force in big biting mammals and the prediction of predatory behaviour in fossil taxa[END_REF]. In the case of the Malayan bears, when corrected for the mechanical advantage of the rear molar, the in vivo values of 1907-2021N are within the upper range of the estimated maximal force of 1722±423N.

For the Andean bears, even when corrected for the mechanical advantage of the rear molar, the maximal predicted bite force at 1236±296N remains well below the in vivo forces. We thus find that the one-sided bite forces we recorded are higher than the estimates reported in the literature. Such modeled bite forces may lead us to underestimate the potential diet range and defensive capability. Our results thus highlight the need for further in vivo measurement to validate and improve biomechanical models.

To train the bears to use the bite plates correctly, rather than using their paws, pulling, or biting at ineffective orientations (see supplementary video), we found that using only the manual control of the internal reward system of the device was most effective. We found our initial training regime using both internal and external rewards focused the attention of the bears on the trainers rather than their interaction with the device, leading to a longer training time. By concealing the manual reward trigger button in a pocket and observing the interaction of the animals with the device from a distance, we were able to train the correct bite orientation rapidly without the animals associating the reward with humans. Eventually this may even be done from a separate chamber, as the device is fitted with a camera to allow observation of the behavior remotely. Once trained, two of the bears and the lioness were able to use the device correctly for several subsequent days without any human presence, while the other two required a brief daily re-training. We found the fraction of incorrect bites to be relatively low (5-18%), and these were generally of a much lower force than valid bites. Since the device provided novelty and stimulated exploration, zoos regarded its use as habitat enrichment.

To be able to use these data for biomechanical modeling, it is necessary to know which teeth were used during the biting . Our videos show that once trained, the bears quickly favored using the rear molars to achieve the high bite forces necessary to receive a reward. Thus, by rewarding increasing force, our reward protocol naturally led to the use of the rear molars, which have the best mechanical advantage to achieve high forces. A session of an hour is sufficient to record several bites near the maximum bite force in a well-trained and sufficiently motivated animal.

Overall, we found the duration of the force peaks to be short, with mean peak widths of 48-69ms per bite, and peak force plateaus of ~10ms. Such rapid bite peaks make the measurement of bite forces using slower sampling rates difficult. It also suggests that previously recorded bite forces in mammals using lower sampling rates may be underestimates of the maximal forces.

This method may not produce bite forces as may be produced during aggressive or defensive behavior, as the motivation is of a different nature. Individual differences in motivation by the reward or ability to learn may also lead to underestimates of the maximal bite force. We therefore recommend to measure multiple specimens of the same species to obtain better estimates of the maximal capacity.

Our reward protocol is an attempt to mimic a natural situation where a hard or tough food item requires increasing chewing force to extract food. The biting behavior is driven by the desire to obtain the food reward. The recorded bite forces are therefore relevant in a feeding context, rather than the defensive bites often measured using coercion methods. To measure voluntary bite forces, we left the device present in the animals' normal enclosures without human presence.

By simultaneously recording force and video, we were able to identify the individual and assess

  the left and the right jaw while controlling for species and specimen (model: logforce_N ~ side + species / (1|individual)). To test for fatigue, we ran a linear model of bite force against the time since the previous bite in seconds, controlling for individual and species as nested factors (model: logforce_N ~ logtimeDiff + species / (1|individual)). Finally, we tested for any effect of the orientation of the bite plates while controlling for the individual bear (model: logforce_N ~ plateUD + (1|individual)).

Figure 1

 1 Figure 1 Schematic drawing of the bite plates. Scale bar is 20cm. On the left is a cross section of the device at the dashed line. The reward tube (not shown) opens between the bite plates on the left side of the device. Scale bar: 20cm.

Figure 2

 2 Figure 2 Exemplary bite force events for each species. Left column shows the recorded forces, the right column shows video stills of the corresponding event. From top to bottom: Spectacled bear Teju, Malayan bear Cindy, and lioness Savannah. The red line shows the unfiltered data, the blue line represents the filtered data.

Table 1

 1 Overview of results per individual. The lioness only used canines or incisors to bite. "i" indicates invalid bites and other force generating events, such as paw presses.

		Tables
	Species Individual Bites	Maximal force (N)	Mean Force (N)	Peak width (ms)	# of bites over 80% of max	L	# of bites R	i	L	Max force (N) R	i
	Andean	Ariba	352	2196.8	692.1	69	35 (9.9%)	151	137	64	2128	2197 1857
	Andean	Teju	3962	2167.4	996.0	62	767 (19.4%) 2640 1140 182	2133	2167 1966
	Malayan Linh	1696	1907.3	539.8	48	34 (2%)	724	872 100	1907	1811 1706
	Malayan Cindy	602	2020.6	1057.7	48	62 (10.3%)	310	249	43	1944	2021 1689
	Lion	Savannah	852	1593.8	378.9	51	4 (0.47%)

the validity and orientation of each bite. Our camera angle did not allow us to measure the jaw opening angle, but we would recommend recording this if the data are to be used for detailed biomechanical modelling. Although we did not record the body mass of the individuals in this pilot study, placing a force plate near the bite device would allow the body mass to be measured quickly and unobtrusively. We expect this method to be applicable to other large carnivores and omnivores. Herbivores, with their lower gape angle and narrower mouth opening, may require a different design of the bite plates and reward system.

Acknowledgements

We are greatly indebted to José Kok and Dorine van Appeldoorn of Ouwehands Zoo Rhenen, the Netherlands, to Saskia Visser-van Soest, Rob Saris and Kiki Klomp of DierenPark Amersfoort, the Netherlands, and to Nuno Alvura, Pedro Mendes and Paula Telinhos of Maia Zoo, Portugal, for access to the animals, and their generous support, time, and assistance.

Thanks to Anne-Claire Fabre for 3D scanning the bear skulls at the MNHN collection.

Funding

JCGG is currently supported by a doctoral scholarship from MinCiencias -Colfuturo "Fondo de Ciencia Tecnología e Innovación del Sistema General de Regalías para la conformación de una lista de proyectos elegibles para ser movilizados, priorizados y aprobados por el OCAD en el programa de Becas de Excelencia Doctoral del Bicentenario". AvdM is financed through FCT (Fundação para a Ciência e a Tecnologia, Portugal), I.P. under contract number DL57/2016/CP1440/CT0009.

Data Availability

All data is enclosed in the manuscript and the supplementary materials. The custom Python script is available from the lead author upon request.

Appendix

Supplementary file 1 -Schematic and photograpsh of the device Supplementary file 2 -Temporal distribution of bites over a day Supplementary file 3 -Data in tab-delimited format Supplementary video: https://youtu.be/BKEvEc1eal4

Competing interests

No competing interests declared

Author Contributions

12 Arie van der Meijden, Julio César González-Gómez, María D. Pulido-Osorio and Anthony Herrel conceived the idea. Arie van der Meijden designed the methodology, constructed the device and collected and analyzed the data. Arie van der Meijden led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.

Ethical statement

This research does not require ethical approval under Portuguese law (Decreto-Lei nº 113/2013, de 7 de Agosto).