Rat superficial masseter operates at long, unstable lengths during biting

Rat superficial masseter operates at long, unstable lengths during biting | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Rat superficial masseter operates at long, unstable lengths during biting Nicolai Konow, Brandon Reder, Daniel Bartlett, Devin Jenness, and 3 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-6739430/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 30 Oct, 2025 Read the published version in Scientific Reports → Version 1 posted 9 You are reading this latest preprint version Abstract The operating length ranges of mammalian jaw muscles have been estimated using twitch contractions or force measurements at the bite point, prompting a consensus that jaw muscles operate at short lengths on their force-length (FL) curve. However, since activation intensity truncates muscle optimal length (LO), we hypothesized that LO of rat superficial masseter (SM) would decrease with activation intensity, with high-force biting involving muscle shortening from long lengths on the FL curve. We measured muscle activation, strain, and force in vivo during biting on food with varying hardness and mapped the in vivo data from each muscle (N = 6) onto its FL relationship, measured in situ. Submaximal LO was approx. 12% shorter than twitch LO, and SM bite forces averaged 4.1 ± 3.9 N/cm2 (mean ± S.D.) and reached 10.6 N/cm2, corresponding to muscle activation and food hardness. Length operation ranged from 7% below LO (ascending FL plateau), to 27% beyond (descending limb). The finding that jaw muscles operate at long, potentially unstable lengths, particularly during hard food biting significantly expands our understanding of skeletal muscle function, with broad implications for craniofacial evolution, muscle mechanics and control, and healthy as well as pathological function of the jaw musculoskeletal system. Biological sciences/Evolution Biological sciences/Physiology Force-length operation length-tension relationship musculoskeletal function Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Introduction Skeletal movements in vertebrates are driven by muscle contractions, the force generation of which is limited by the overlap of actin and myosin filaments in the sarcomere, in accordance with the force-length (FL) relationship [1-4]. An emerging consensus has been that across scales of biological organization, sarcomeres, fibers, and the fascicles (bundles of fibers) they compose tend to operate at optimal lengths ( L O ), centered on the FL curve plateau, where skeletal muscle generates peak force ( F O ) [5]. However, under scenarios where stable muscle function is priority, fascicle operating length ranges may be truncated to shorter lengths of the ascending FL limb [6-8]. At such short lengths, additional shortening reduces muscle force as myofilament overlap increases [2, 9]. Conversely, muscle may also operate at long lengths, on the descending FL limb [8]. Here, muscle shortening up the descending FL limb increases force, but also involves risks of sarcomere disruption due to the reduced myofilament overlap [10]. Considering the risks of myofilament disruption, muscles have been thought to only operate at long lengths if they are unlikely to undergo abrupt lengthening. The prevailing consensus on muscle operating length ranges comes from studies of limb and cardiac muscle. Therefore, it remains unclear if the consensus about muscle operating lengths transfers directly to jaw muscle operation or not. Determining whether jaw muscles operate at short lengths, or at long and potentially unstable lengths would significantly expand our understanding of skeletal muscle function. This would also have broad implications for craniofacial evolution, muscle mechanics, and control, and for both healthy and pathological function of the craniofacial musculoskeletal system. Jaw muscles studies have suggested that force generation during biting at large gapes involve whole muscle and fascicle operating lengths that are short and centered on the ascending FL limb [11-16]. This mode of operation is consistent with the idea that jaw muscles have relatively large moment arms, which result in large muscle length-changes [17] that inherently could render sarcomeres mechanically unstable. However, jaw muscles rarely experience abrupt lengthening (i.e., eccentric contractions, except in cases of live prey struggling). This musculature would thus be expected to operate at FL plateau-lengths, because FL effects associated with muscle shortening as gape approaches occlusion would reduce force production dramatically [9]. If the SM operates at short ascending limb lengths, force loss during gape closing would be particularly challenging for incisor biting at large gape, as used by many mammals [18] including rodents. However, approaches in earlier jaw muscle studies, which include force measurements from the entire musculature at the bite point [13], estimations of muscle length from skeletal displacement, and elicitation of twitch contractions to measure FL from individual muscles [11, 14, 16], are unlikely to accurately reflect the patterns of individual muscle activation, contraction, and force generation used during feeding. Measurements of muscle activation, fascicle strains, and time-varying force production are critical for reliable interpretations of muscle function in vivo , for several reasons: (i) The length trajectory of muscle fascicles often differs considerably from the entire muscle-tendon complex, due to the elastic action of compliant tendinous elements [19, 20], and architectural gearing effects caused by fascicles rotating around their tendinous attachments [21-23]. (ii) In most skeletal systems, including jaws, force measured at the contact point can be generated by several muscles whose individual contributions cannot readily be disentangled from time-varying bite-point force [24]. (iii) Muscles are typically activated in a set sequence and at intensities corresponding to the mechanical task. Limb muscle studies have shown that muscle activation intensity impacts the operating length-range of its fascicles [25]. Tetanic activation results in a left-shift of the FL plateau by 10-15% relative to the L O of a FL curve generated with twitch contractions [26, 27]. However, a general lack of in vivo measurements of muscle activation, fascicle strains, and time-varying force generation from individual muscles of the jaw means that we do not know if they undergo the same activation-dependent operating length-shift. It also remains unknown what the FL effect on jaw muscle operation might be. We hypothesize that conclusions drawn in earlier jaw studies, suggesting that muscles of the jaw operate at relatively short lengths [11, 14, 16], were influenced by the twitch activation protocol used. We predict that operation at short twitch lengths actually correspond to operation on the FL plateau or longer ascending limb lengths under greater (submaximal tetanic) activation conditions, as necessary for generating sufficient bite force. To determine if earlier study protocols led to an imprecise consensus about the length-operation range of jaw muscles, we elicited incisor biting in rats feeding on an assortment of natural food items. We used food types that rats willingly accept with material properties, including systematic hardness differences that have been thoroughly characterized, to facilitate comparisons of our data with existing literature [28, 29]. We used direct instrumentation of the rat superficial masseter (SM) to collect direct measurements of time-varying force production from a major jaw closer muscle during feeding, along with electromyographical measurements of muscle activation and recruitment (EMG), and measurements of fascicle length-change dynamics (sonomicrometry). We then tested the same muscle preparations with an in situ ergometry experiment that measures twitch and submaximal tetanic FL. Finally, we mapped in vivo fascicle length-operating ranges and forces to the in situ FL curves [8, 20, 30]. For limb muscles, this approach has proven powerful by determining in vivo constraints on muscle force and fascicle length, which we have previously only indirectly achieved for jaw muscles [9]. This experimental design allowed us to test the hypothesis that large fascicle strains in SM, associated with incisor biting at large gape, lead to peak force generation at fascicle operating lengths that correspond to the descending limb on the submaximal tetanic FL curve, in contrast to the findings of earlier studies that used twitch contractions. Results Across the six rats tested in vivo in biting experiments, we were not always able to collect reliable bites for all food types, and a total of 247 bite recordings were used for analyses (see Supplemental Table 3). We found a strong relationship between recruitment of rat SM and its force production across incisor bites on foods with varying hardness. Food hardness correlated strongly with muscle recruitment intensity (GLM: F 2,236 = 10.8; p < 0.0001; η 2 = 0.32), measured as the integrated area under the rectified EMG signal (Fig. 2A). Harder food was associated with greater muscle force production (Fig. 2B), again with a statistically significant relationship (GLM: F 2,161 = 4.8; p < 0.01; η 2 = 0.07). Whereas EMG activity and force production generally increased with food hardness, this relationship plateaued at levels of hardness above 20 MPa. Rat chow, nominally the hardest food tested (50.44 MPa; [29]), elicited lower-than-expected EMG and force production compared to predictions based on muscle force production on softer foods. Stress developed in vivo by the rat SM increased considerably with food hardness and averaged 4.1 ± 3.9 N/cm 2 (mean ± S.D.). Greater stresses were consistently recorded for the two hardest food items (almond and chow) and peaked for almond at 10.6 N/cm 2 , with median stress for chow bites (7.2 ± 3.9 N/cm 2 ) being marginally lower than for almond bites (8.6 ± 5.3 N/cm 2 ) (F 2,125 = 4.966; p <0.05). Based on an average muscle PCSA of 0.51 ± 0.12 cm 2 (Supplemental Table 1), and specific tension data for vertebrate skeletal muscle [31, 32], the rat SM operated at submaximal recruitment during incisor biting. A Linear Least-Squares regression of EMG recruitment against force production indicated a strong EMG-force relationship (F 1,161 = 516.9; p < 0.0001), despite identifiable variation across subjects and food hardness (Fig. 2C). The operating length-ranges of fascicles in the rat SM were strongly correlated with food hardness. Peak relative lengths ( L MAX ) typically occurred at bite onset (Fig. 1B) and preceded as well as exceeded relative fascicle length at activation ( L @â ON ), likely due to earlier activation and action of different adductor muscles in the rat feeding system [33]. Increased food hardness elicited longer relative fascicle lengths at muscle activation onset ( L / L O @â ON ) (GLM: F 2,237 = 52.1; p < 0.0001; η 2 = 0.98) (Fig. 3A) and at occlusion ( L / L O @F MAX ) (F 2,237 = 74.4; p < 0.0001; η 2 = 0.68) (Fig. 3B). The strain range for contracting fascicles varied considerably (Fig. 3C) and averaged 10.5 ± 3.2%, with a statistically significant reduction in strain with increasing food hardness (F 2,237 = 24.8; p < 0.0001; η 2 = 0.69). Our in situ preparations of rat superficial masseters achieved peak tension development under twitch activation at optimal length ( L O ) reaching 12.3 ± 4.3% relative to tetanic activation ( F O ). Tetanic tension at L O averaged 10.1 ± 3.3 N/cm 2 , suggesting that our needle stimulation protocol elicited submaximal and not maximal tetanic contractions. Comparisons of L O for twitch and submaximal tetanic FL curves showed that the submaximal tetanic plateau was left-shifted by 11.3 ± 2.9% with respect to the twitch plateau (Fig. 4). Superimposition of in situ and in vivo data from rat SM (Fig. 5) revealed incisor biting to involve fascicle operating strain-ranges spanning from in excess of 25% over L O , on the descending FL limb (i.e., high ε) across the FL plateau to approximately 7% down the slow ascending FL limb (i.e., low ε). Both hard and soft food types elicited bites involving high fascicle peak strains (ε MAX ). However, fascicle operation at low strains were particularly evident for soft food types and fascicle operation at high strains were particularly common for hard food types, which also elicited the greatest tension generation as the jaw approached occlusion. Discussion Our in vivo study of the rat superficial masseter during incisorbiting generated the first available ensemble measurements of electromyographical activation, fascicle strain, and force from a jaw muscle during feeding. We mapped the in vivo stress and operating strain ranges of the same muscle onto its submaximal force-length (FL) relationship, as measured from the same muscle preparations in situ. Combined, these data revealed that as activation increases, the FL plateau of rat superficial masseter (SM) shifts leftward to a similar extent as previously reported for various limb muscles. Taken together, our data reveal that the rat SM operates in vivo across lengths that range from the upper FL ascending limb, across the plateau, and down the descending limb. This finding runs counter to the well-established paradigm that jaw muscles operate a short, stable lengths [11-16]. Our findings further develop our understanding of how skeletal muscle exploits the FL relationship in vivo to maximize performance for specific tasks, in this case biting. These results highlight the role of muscle operating length in maintaining musculoskeletal stability during routine motions. Our finding has broad implications, ranging from reconstruction of craniofacial evolution in vertebrates to the prevention and treatment of temporomandibular joint disorders (TMD) in humans. Muscle operating lengths depend on the intensity of muscle activation and recruitment. For the rat SM, we discovered a left-shift of the plateau of the muscle’s FL curve averaging 11.3 ± 2.9% L O as we increased stimulation voltage in situ, from twitch to submaximal tetanic stimulation. This shift in optimal length ( L O ) is consistent with earlier experiments that varied stimulation voltage to elicit twitch, submaximal tetanic, and maximal tetanic contractions in limb muscles, including toad plantaris (>15% L O ,[34]), Guinea Fowl lateral gastrocnemius (20% L O , [35]). rat medial gastrocnemius (12.6-11.2% L O , across muscle compartments, [27]) and mouse soleus and tibialis anterior (both 8.2% L O , [26]). For cat soleus, a similar shift wad demonstrated by varying stimulation frequency (25% L O , [36]). Our discovery of a plateau-shift for a jaw muscle that is comparable to measurements from limb muscles with a diversity of fascicle architectures, series elastic element prominences, and mechanical functions suggest that this trait is generalizable across striated skeletal muscle. The mechanistic explanation for the plateau-shift remains incomplete as both strain-dependent calcium sensitivity-changes affecting cross-bridge cycling [25, 36-38] and series elastic compliance [38-40] have been implicated. Regardless, our findings from rat SM suggest that studies using twitch activation to examine FL behavior in jaw muscles [14, 16] have generated an inaccurate consensus about the lengths across which jaw adductor muscles operate. Our conclusion is strengthened by the consideration that muscles rarely are twitch-activated nor supramaximally activated in vivo and typically operate at submaximal levels of activation and recruitment [27]. The central nervous system regulates muscle activation timing and both the intensity and completeness of muscle recruitment to changes in loading environment and force requirements [41-43], for instance, as elicited by changes in food hardness. In agreement with this principle, we found a strong positive relationship between the recruitment intensity (integrated EMG area) of the rat SM during in vivo incisor biting power-strokes with the hardness of the food substrate (LLS regression, p < 0001). A similar strong positive relationship was also established between muscle force output and food hardness, with few exceptions (discussed below). Overall, our findings indicate a generally strong agreement between EMG and force measurements, consistent with recordings from a diversity of limb musculoskeletal systems [44, 45]. This finding suggests that EMG holds great promise for understanding time-varying force production in other jaw muscles, where substrates for force-buckle attachments are unavailable, such as the pterygoid and temporalis. As food hardness was increased during rat in vivo biting, we found systematic changes in the operating length ranges of the muscle fascicles, as measured directly using sonomicrometry. As recruitment intensity and force production increased with food hardness, muscle fascicles operated across longer length-ranges but also contracted to relatively shorter lengths (Fig. 3B; GLM; both p < 0.0001). These findings underscore that muscles can exploit FL effects by tuning operating length to generate forces that match task-based requirements, contend with length-based limitations to force generation, and mitigate risks of contracting at long lengths [8, 9]. To better understand how the rat SM exploits its FL relationship during feeding, we superimposed our in vivo and in situ measurements. This is a validated and powerful method for rigorously interpreting operating requirements and constraints on skeletal muscles in limbs [20, 35, 46]. Mapping the in vivo operating length ranges of rat SM onto the muscle’s own FL relationships, as measured in situ allowed testing of the hypothesis that twitch FL studies inaccurately reflect the muscles’ contractile behavior. In line with earlier twitch studies of jaw muscles [14, 16], we determined that the in vivo operating length-ranges of the rat SM during incisor biting corresponded to the ascending limb (i.e. short operating lengths) of the twitch FL curve. However, when the in vivo operating length ranges were mapped to the significantly shorter plateau-lengths of the submaximal FL curve, with comparable stresses, we found that the realized operating lengths were typically distributed across the FL plateau. For bites on hard food (almond, chow) the muscle often operated at relatively high stresses on the descending limb (i.e., long lengths). This discovery is important for several reasons. Limb muscles will often, but not always, operate across the FL plateau where muscle force is maximized [8, 47]. Our data, therefore, suggest general similarities between limb and jaw muscles in this regard (see also [9]). We tested natural foods that rats would willingly accept and attempted to size-match food items to avoid food size-driven FL effects in our in vivo data. Evidence that our attempts to decouple muscle force from gape and fascicle length were not entirely successful is seen for our hardest food type, chow (E = 50.44 MPa [29]). Chow elicited significantly lower muscle recruitment intensity and force production (Figure 2) than our second-hardest food type, almond (E = 19.4 MPa). This result may be influenced by uncontrolled food properties such as moisture content, with which hardness likely co-varies [48] as well as the granular nature of chow that may afford the incisors easier purchase, and also fractures more easily between incisors than between materials testing plates. However, food-driven FL effects were likely not completely avoided as chow pieces invariably were larger than raisins. Consequently, length effects on muscle function, driven by variations of food size [9] likely also play a role because relative fascicle operating lengths were considerably longer for chow than for other food types (Fig. 3A, B). Long fascicle lengths may have prevented greater force generation on chow (Fig. 2A,B), although the associated depression in force likely was somewhat offset by force-velocity effects [46]. For example, muscle operating strains on chow were generally shorter (Fig. 3C) and slower (Fig. 1B) than for other food types. Conversely, our softest food type, raisin (E = 0.22 MPa) elicited greater-than-expected muscle forces (Fig. 1B) and more variable operating lengths (Fig. 2A,B) than other food types. This result may be caused by the greater compliance (i.e., stress to fracture) of raisins, requiring greater biting work than other food types tested (Supplemental Table 2). Together, the anomalous force-responses by SM on our nominally softest and hardest food type respectively may be the reason that the influence of food hardness on SM force production carried the lowest effect size (η 2 = 0.07) of our statistical comparisons. Future work building on ours could use manufactured food items that permit more careful control and parameterization of shape, size, and hardness [9, 49] than natural foods do. Analyzes could include single bite interactions on spherical food with controlled variations of size and/or hardness. Such analyses would also benefit from videofluoroscopic imaging [9] to precisely determine where the incisors make purchase on the food item and jaw gape. They should also incorporate skeletal and muscle moment arm analyses [17], especially considering that rodent jaw joint anatomy allows propalinal (caudal-to-rostral) jaw translation during gape closing, which likely affects muscle operating length independent from vertical gape movements. Foods used in our study are unlikely to have challenged the rat jaw apparatus, which for instance is capable of fracturing pistachio shell (E = 1.2 GPa; Pers. Obs.). Accordingly, forces measured in vivo were unlikely representative of the peak voluntary contractile capacity of the rat SM, which may be realized not only during feeding but also during defensive biting. Here, the maximal FL curve would be relevant, and our specific tension calculations suggest that the maximal FL plateau would be left-shifted by 10% from submaximal L o, placing peak bite force events on the deep descending limb. We were unable to elicit supramaximal contractions in situ due to the surgical infeasibility of accessing the nerve of masseter for mounting of a cuff electrode that would facilitate supramaximal stimulation of the rat SM. However, given the well-matched peak tension recorded in vivo during biting and elicited in our in situ tetanic contractions, the methodological limitations of needle electrode stimulation are unlikely to affect the conclusions of our study. Despite the good agreement between the specific tension values we recorded in vivo and the peak tension elicited in situ , future studies should explore alternative stimulation approaches [14] to elicit supramaximal stimulation and generate in situ FL data that match supramaximal behaviors. The discovery that the rat SM operated at long fascicle and sarcomere lengths to generate sufficient force to fracture hard food is an interesting finding from the perspectives of muscle as well as jaw joint stability and safety. Compared to short fascicle and sarcomere length operation, where muscle is stable [6], large, hard food items drove the muscle to develop high forces at long lengths where inherent sarcomere instability due to weak actomyosin overlap [2] might cause muscle damage [10]. The dangers of the descending FL limb have been highlighted for limb muscle [50] but may potentially be greater for jaw muscles due to their larger moment arms [17]. If food item slippage was to translate the jaw during high force production at large gape, the resulting stretch of the muscle might compromise its ability to counter dislocating torques at the temporomandibular joint (TMJ). Other potential deleterious results from such an event include damage to the dentition [16]. Jaw muscles in mammals tend to alternate between powering and balancing the jaw during feeding [24, 51, 52]. Accordingly, muscle instability might trigger jaw joint instability and our data suggest that biting on large and hard objects could be an overlooked trigger of dysphagia including temporomandibular joint disorders (TMD). We anticipate future FL operation analyses of mammalian jaw muscles to inform research on temporomandibular joint (TMJ) function and triggers of its dysfunction. Our data on the effects of food types with different material properties echo previous work by paleontologists and anthropologists seeking to reconstruct foraging and feeding in our mammalian and human forebearers [28-29]. However, such studies have previously lacked the analytical power that single muscle force measurements provide [20,51,55]. The force-length effects highlighted by our data are especially relevant for computational modeling of the craniofacial system, where the force-length relationships of individual muscles tend to be greatly simplified and certainly do not incorporate factors such as activation-dependent shifts of optimal length [53-55]. Computational models are a common tool given the challenging nature of experimental work on primates (including humans) and are the only tool available for studying craniofacial evolution in the fossil record. Several key transitions in vertebrate evolution involve changes to the mandible and jaw joint, where continued functioning and stability of the system would be paramount e.g., mandibular reduction in modern humans, and the transformation of the post-dentary bones into the mammalian middle-ear ossicles [56, 57]. Our experimental data demonstrate a clear need to fully understand jaw biomechanical systems to inform modeling approaches aiming to accurately reconstruct key functional transitions across major clades. Methods Sprague-Dawley rats ( N = 10, 50/50 M/F) were purchased from Charles River and same-sex pair-housed under a 12 H light : 12 H dark cycle at the Concord Field Station (Harvard University). All methods were carried out in accordance with relevant guidelines and regulations. Methods are reported in accordance with ARRIVE guidelines. Both husbandry and experimental procedures adhered to protocols approved by the Institutional Animal Care and Use Committees at Harvard University (FAS, 20-09_04) and UMass Lowell (23-03). Rats were provisioned with water and rodent chow ad libitum and at the time of experiments, body masses were 440 ± 113 g (mean ± S.D.) and ages were 13.5 ± 1.8 weeks (mean ± S.D.) (Supplemental Data Table 1). Power analysis based on muscle activation, strain, and force data from our pilot study suggested that N =5 would yield a Delta (mean difference over individual variation) for two-factored GLM designs of 1.8-2.5 and analysis power ranging between 87% and 96% at α = 0.01 (Systat v. 12). Instrumentation surgery Rats ( N = 6) were induced in a chamber on isoflurane (4-4.5%) via oxygen (2.0L/min) and transferred to a low-profile mask for anesthetic maintenance at 1.5-2.5% isoflurane via oxygen (0.8 L/min). Core temperature was monitored via a rectal thermocouple and adjusted using an infrared lamp and a controllable heat pad. Breathing rate and core temperature were continuously monitored and recorded at 10 min intervals throughout each procedure. After application of ophthalmic ointment and injection of 2ml sterile saline (SQ), incision areas were shaved and sterilized using three repeated alternating scrubs of chlorhexidine and betadine. The subject was draped with semitransparent plastic (Press’n Seal) to allow visual monitoring of breathing during the sterile implantation procedure that followed. Two skin incisions were made, one along the left lateral-to-ventral chin and one along the mid-dorsum of the neck to route wires subcutaneously through a sterile plastic straw, from a skin port (Microtech, Bronwyn PA) attached at the dorsal incision, to the chin incision, over the SM (Fig. 1A). A pair of sonomicrometry crystals (1 mm, Sonometrics Corporation, London ON) were implanted into the SM muscle along the length of a fascicle. One crystal was implanted at the myotendinous junction, and one at the mandibular origin. The distance between the crystals were 12 ± 3 mm. Crystals were anchored to the epimysium using 6-0 silk and oriented in the tissue to generate the strongest possible transducer pulse (monitored during surgery via an oscilloscope). To measure EMG, a bipolar, offset-hooked electrode [58] (California fine wire, Grover beach, CA) was implanted using a 25G hypodermic needle mid-belly into the muscle, along the same fascicle as the sonomicrometry. The EMG electrode was also sutured to the epimysium using 6-0 silk. A custom-made miniature strain buckle (FLK-1–11, Tokyo Sokki Kenkyujo) [59] was tied to the tendon of SM using 3-0 silk. The rat and its instrumentation was grounded using a large-gauge blunt monopolar steel electrode located under the skin port chassis. Both skin incisions were closed using 3-0 Vicryl. Flunixin Melgumine (30mg/kg, SQ) was administered as pain relief and analgesics intraoperatively, and at 12 hour intervals throughout recovery. Once awakened from surgery, subjects were fed soft (raisin) food within the hour to ensure recovery to normal behavior, and then fasted during recovery (overnight, 12-16 hours) to elicit optimal feeding performance during the in vivo feeding experiment. In vivo feeding experiment Our in vivo experiment was designed to test bite performance across a range of food types with known hardness (Range of Young’s Modulus; 0.22 – 50.44MPa; [28, 29, 60] (Supplemental data Table 2) and near-equivalent size (approx. 1.5 cm diameter). We avoided mechanical alteration of food size due to concerns that this would affect the opportunities for the subject’s incisors to make purchase on the food during biting. During experiments in vivo , the subject was placed in a transparent acrylic enclosure (25 x 25 x 25 cm). The data collection cable was connected to the skin port and the recording equipment and freely suspended above the subject’s neck using nylon string. The rat was presented with food items (Supplemental data Table 2) in a pseudorandomized sequence from a sturdily attached, spring-loaded metal (alligator) clip. Clipping of food was chosen to facilitate video recordings of the bite without the forelimbs obscuring the view of the incisors. Food clipping also prevented the rat from generating bite torque using the upper incisors and the forelimbs, which potentially would decouple jaw adductor muscle function from the biting action. A single bite series was allowed per food item to minimize FL effects in the muscle data caused by food size reduction due to repetitive biting. Figure 1B shows individual bites on different foods. To collect gape cycle data, we recorded high speed video of all bites in a quasi-ventral view (125fps, 1024x1024 pixels, Fastcam SA3, Photron, San Diego, CA). Sonomicrometry data were acquired at 520 Hz (transmit pulse: 406.25 ns; inhibit delay: 1.54 mm; range-gain: 25mm) using a digital sonomicrometer and a PC laptop running SonoLab (Sonometrics Inc., London, ON). Muscle EMG data were recorded at 100x amplification (Grass P511) with a 10 kHz high-cut and a 60Hz notch filter engaged. Strain gauge measurements were balanced and amplified (200x; Vishay Measurements Group 2120, Raleigh, NC) and all time-varying transducer data were digitized at 4kHz via a Powerlab 16/20 onto a PC laptop running LabChart 8.0 (ADInstruments). The animal ground was connected to the equipment ground via the sonomicrometer cable shield. A manually triggered square-wave TTL signal terminated video recordings and permitted resynchronization of time-varying transducer and gape cycle data. In situ FL ergometry Due to in situ preparation and sonomicrometry transduction difficulties for two of our in vivo rats, we used additional rats ( N = 4) to generate FL curves in situ . For these rats, sonomicrometry implantations were carried out during the ergometry preparation and EMG and force buckles were not mounted. Directly following the in vivo data collection, each subject was re-induced and maintained on isoflurane, including vitals controls and monitoring as described above. The skin incision on the chin was carefully reopened, and the tendon of the superficial masseter was gently dissected free from its insertion on the rostrum and connected to the lever-arm of a servomotor (Aurora Scientific 300C-LR) using a self-tightening knot of 1-0 silk and a dab pf cyanoacrylate gel (Duro) (Fig. 1C). Two pairs of stimulation electrodes were implanted and secured in place using 6-0 silk. Care was taken to place equivalent poles of each pair at the rostral (+) and caudal (-) ends of the muscle respectively. The insertion orientation of each electrode was normal to the fascicle plane, to avoid alterations in stimulation and recruitment as fascicles shortened against series elastic compliance. The stimulation electrode leads were flexibly suspended above the muscle preparation to minimize interference by motion during contractions. The left zygoma and ventral mandibular corpus were both firmly clamped to a customized stereotaxic frame onto which the servomotor was mounted. Muscle temperature was maintained as constant as possible using an infrared lamp and a dressing of saline soaked gauze prevented muscle desiccation during the procedure. We acquired data on muscle stimulation intensity (V), muscle fascicle length (measured via the sonomicrometry crystals, mm), and motor force (N) along with measurements from all implanted transducers at 4kHz via a NIDAQ (USB-6212) to a PC running IgorPro 8.0 (Wavemetrics, Lake Oswego, OR). Optimal stimulation parameters were determined for each preparation using a series of twitches with increasing voltage (6-12 volts depending on preparation; 125 pulses/s, 0.2 ms pulse width; Grass s48). Four to five 300 ms contractions from a passive tension that roughly corresponded to L O were used remove series compliance and collect time-varying data from the tendon strain-gauge (V) and the motor (N) for calibration purposes (see below) [59]. Video recordings of sample tetanic contractions in two preparations were used to verify that the fixed-end compliance of the preparation was minimal (measured to be 0.7 ± 0.2 mm). Using consecutive twitches, the preparation was shortened to minimal length ( L REF ), which was measured from the muscle origin on the mandible to the suture knot (i.e., the entire preparation) to the nearest millimeter using calipers. Starting at L REF , we then collected the twitch and submaximal tetanic FL curves for the preparation, and the fascicles, as measured via sonomicrometry (Fig. 1D). Three minutes of rest was allowed between each tetanic stimulation to reduce effects of muscle fatigue. A final tetanic contraction at L O was used to confirm that the preparation had lost no more than 10% of F O during testing, which was confirmed for all preparations. Following the in situ experiment, the subject was euthanized (Sodium Pentobarbital, IC). The muscle was excised to record fascicle length, including the segment of the fascicles that were not spanned by the crystals, so that a segment correction could be applied to the in vivo data. We also measured resting muscle pennation angle and muscle mass (without the tendon) to calculate muscle physiological cross sectional area (PCSA) [61], which was used to calculate muscle stress (N/cm 2 ). Data extraction & analysis Construction of force-length curves We constructed force-length (FL) curves to determine twitch and submaximal tetanic optimal length ( L O ) from our sonomicrometry measurements of fascicle length-change and peak force (P O ) from the motor and muscle force buckle for each SM preparation by fitting third order polynomials to the passive and active twitch and tetanic curves in IgorPRO v.8. Force (N) was expressed as stress (N/cm 2 ) by dividing with muscle PCSA (Supplemental table 1) and fascicle length (mm) was expressed as strain (( L - L O )/ L O ). To measure the percent left-shift of tetanic L O relative to twitch L O , we used tetanic L O to strain-normalize twitch FL curves [25, 26]. Passive FL curves were adjusted for shortening during stimulation [62]. Conditioning of in vivo measurements Gape cycles were digitized for usable bites in our custom confluence of DeepLabCut and XMAlab [63]. Time-varying traces of gape cycles were imported into IgorPro along with CSV data from Powerlab containing activation, length, and force measurements from the implanted transducers. Gape and instrument data were re-synchronized using the TTL square-wave signal. All subsequent signal conditioning and analyses were performed in IgorPro. Muscle EMG signals were filtered as needed using an FIR interpolation, base-line corrected and rectified. Sonomicrometry signals (V) were converted to length (mm) using a synchronized calibration recording (in SonoLab). These signals were then segment-corrected for the fascicle length not spanned by the crystal pair, and adjusted for a time-offset caused by the crystal coating epoxy [64]. The conditioned signal was smoothed using a quintic spline interpolation (s.d. = 0.01 – 0.05) and strain (ε = ( L - L O )/ L O ) was calculated using tetanic L O . To calibrate the tendon buckle signal (V) from in vivo recordings to force (N), we first regressed the time-varying gauge signal (V) against motor force (N) from the initial in situ tetanic contractions. Linear regressions of these two time-varying measurements for the period of tetanic force rise returned r 2 values ranging from 0.87 to 0.96, which were comparable to prior validations of this technology [45, 59, 65]. The average regression slope (V/N) was then used to calculate time-varying in vivo muscle-tendon force (N), which was quintic spline interpolated as needed and baseline-corrected to zero in a period of the in vivo trial where the animal was at rest and no feeding occurred. Extraction of in-vivo measurements Time-varying data from the conditioned and calibrated in vivo measurements of SM activation, strain, and force were graphed in IgorPro. For the period of gape closing, a semi-automated procedure was used to extract the following response variable: (1) Mean EMG (mV), (2) Fascicle starting length at the power stroke onset ( L MAX ), (3) Fascicle length at EMG onset ( L ACT ), (4) Fascicle length at peak force ( L MIN ), (5) Fascicle shortening strain (ε=[ L- L O ]/ L O ) during the power stroke (T ACT to T MIN ), (6) Force rise magnitude (N), and (7) Peak force (F MAX , N). Data were tabulated per individual ( N = 6), cycle number (to control for effects of food size reductions across repeated bites), and food hardness. Data extraction was performed by two researchers who were blinded to the treatment (food hardness) across trials. Hypothesis testing and statistics To test the hypothesis that activation and force production by the muscle fascicles would vary systematically with food hardness, we ran General Linear Models (GLM) in Systat v. 12. These analyses factored EMG (specifically the integrated area under the rectified electromyogram during jaw closing, in mV/s) or muscle force (N) as the response variables, individual as an independent effect, and food hardness (E) as the fixed factor. We interrogated the interaction term (IND*E) for statistical significant. Effect sizes of the interaction terms were calculated as the partial eta squared (η 2 ) in SPSS 30.0. We also probed the strength of the EMG-force relationship across individuals and food types using Linear Least-Squares regression in Systat v. 12. To determine if fascicle relative length ( L / L O ) at activation ( L ACT ) and at peak force generation ( L MIN ), as well as fascicle strain during the incisor biting power stroke (% L O ) increased with to food hardness (E), we ran GLM and effect-size analyses using the design described above. Stress-strain FL curves from the in situ experiments were used to determine if optimal length ( L O ) for rat SM varies with recruitment intensity, as previously reported for limb muscles. To determine this, we ran a GLM factoring L O as the response variable, with stimulation protocol (twitch or submax) as a fixed factor, individual as an independent effect with interrogation of the interaction-term and computation of effect-size as above. To determine the ranges of stresses and strains that the fascicles of the rat SM operate across during incisor biting, we plotted peak muscle stress against fascicle final strain (at peak force generation, Fig. 2B) across food items with varying hardness. Our in vivo fascicle operating strain data were recorded using the same sonomicrometer crystals used to measure fascicle strain in situ . This enabled us to superimpose in vivo stress-strain data onto the aggregate FL curves, expressed as stress versus strain, as measured in situ using submaximal stimulation. Declarations Data access statement: Data analyzed in this study are available from the corresponding author upon reasonable request. Supplemental information: Is available within this submission. Acknowledgements: B. Rimkus, T. Flint, S. Angelli-Nicholls, and K. Katugam-Dechenne, assisted with experiments. G. Crowell, and N. Raymondeau assisted with data extraction from the in vivo measurements. P. Ramirez assisted with husbandry, and A. A. Biewener provided access to facilities at the Concord Field Station. C. F. Ross and N. J. Gidmark both motivated and offered helpful discussions of the study. Funding: This work was funded by the National Science Foundation (award 2217246). Author contributions: Conceptualization: N.K., D. J., & T.P, Methodology: N.K. Investigation: N.K., B.R., D.B., D.J., T.P., J.M. & R.B. Data analysis: N.K, D.J., T. P. Writing—original draft: N.K. Writing—review and editing: N.K. B.R., D.B., D.J., T.P., J.M. & R.B. Competing interests: The authors declare that they have no competing influences. References Ramsey, R.W., and Street, S.F. (1940). The isometric length-tension diagram of isolated skeletal muscle fibers of the frog. Journal of Cellular and Comparative Physiology 15 , 11-34. Gordon, A.M., Huxley, A.F., and Julian, F.J. (1966). The variation in isometric tension with sarcomere length in vertebrate muscle fibres. J Physiol 184 , 170-192. Huxley, H., and Hanson, J. (1954). Changes in the Cross-Striations of Muscle during Contraction and Stretch and their Structural Interpretation. Nature 173 , 973-976. Huxley, A.F., and Niedergerke, R. (1954). Structural changes in muscle during contraction; interference microscopy of living muscle fibres. Nature 173 , 971-973. Rome, L.C., and Sosnicki, A.A. (1991). Myofilament overlap in swimming carp. II. Sarcomere length changes during swimming. Am J Physiol 260 , C289-296. Tu, M.S., and Daniel, T.L. (2004). Cardiac-like behavior of an insect flight muscle. J Exp Biol 207 , 2455-2464. Layland, J., Young, I.S., and Altringham, J.D. (1995). The length dependence of work production in rat papillary muscles in vitro. J Exp Biol 198 , 2491-2499. Azizi, E. (2014). Locomotor function shapes the passive mechanical properties and operating lengths of muscle. Proceedings of the Royal Society B: Biological Sciences 281 , 20132914. Gidmark, N.J., Konow, N., LoPresti, E., and Brainerd, E.L. (2013). Bite force is limited by the force–length relationship of skeletal muscle in black carp, Mylopharyngodon piceus . Biol Lett 9 , 20121181. Morgan, D.L., and Proske, U. (2006). Sarcomere popping requires stretch over a range where total tension decreases with length. J Physiol 574 , 627-628. Anapol, F., and Herring, S. (1989). Length-tension relationships of masseter and digastric muscles of miniature swine during ontogeny. J Exp Biol 143 , 1-16. Weijs, W.A., and van der Wielen-Drent, T.K. (1983). The relationship between sarcomere length and activation pattern in the rabbit masseter muscle. Arch Oral Biol 28 , 307-315. Nordstrom, S.H., and Yemm, R. (1974). The relationship between jaw position and isometric active tension produced by direct stimulation of the rat masseter muscle. Arch Oral Biol 19 , 353-359. Thexton, A.J., and Hiiemae, K.M. (1975). The twitch-contraction characteristics of opossum jaw musculature. Arch Oral Biol 20 , 743-748. Manns, A., Miralles, R., and Palazzi, C. (1979). EMG, bite force, and elongation of the masseter muscle under isometric voluntary contractions and variations of vertical dimension. The Journal of prosthetic dentistry 42 , 674-682. Mackenna, B.R., and Türker, K.S. (1978). Twitch tension in the jaw muscles of the cat at various degrees of mouth opening. Arch Oral Biol 23 , 917-920. Granatosky, M.C., and Ross, C.F. (2020). Differences in muscle mechanics underlie divergent optimality criteria between feeding and locomotor systems. Journal of Anatomy 237 , 1072-1086. Eng, C.M., Ward, S.R., Vinyard, C.J., and Taylor, A.B. (2009). The morphology of the masticatory apparatus facilitates muscle force production at wide jaw gapes in tree-gouging common marmosets (Callithrix jacchus). J Exp Biol 212 , 4040-4055. Rull, M., Solomon, J., and Konow, N. (2020). Elastic recoil action amplifies jaw closing speed in an aquatic feeding salamander. Proceeding of the Royal Society, series B Biology 287 , 32429804. Konow, N., and Roberts, T.J. (2015). The series elastic shock absorber: tendon elasticity modulates energy dissipation by muscle during burst deceleration. Proceedings of the Royal Society B: Biological Sciences 282 , 20142800. Laird, M.F., Granatosky, M.C., Taylor, A.B., and Ross, C.F. (2020). Muscle architecture dynamics modulate performance of the superficial anterior temporalis muscle during chewing in capuchins. Sci Rep-Uk 10 , 6410. Azizi, E., Brainerd, E.L., and Roberts, T.J. (2008). Variable gearing in pennate muscles. Proc Natl Acad Sci U S A 105 , 1745-1750. Laird, M.F., Iriarte-Diaz, J., Byron, C.D., Granatosky, M.C., Taylor, A.B., and Ross, C.F. (2023). Gape drives regional variation in temporalis architectural dynamics in tufted capuchins. 378 , 20220550. Hiiemae, K.M. (1978). Mammalian mastication : a review of the activity of the jaw muscles and the movements they produce in chewing. In Development, Function and Evolution of Teeth, P.M.B.a.K.A. Joysey, ed., pp. 359–398. Holt, N.C., and Azizi, E. (2014). What drives activation-dependent shifts in the force-length curve? Biology letters 10 . Sundar, S., Rimkus, B., Meemaduma, P.S., deLap, S., LaFave, N., Racca, A.W., Hettige, P., Moore, J., Gage, M., Shehaj, A., et al. (2022). Bridging the muscle genome to phenome across multiple biological scales. Journal of Experimental Biology. Tijs, C., Konow, N., and Biewener, A.A. (2021). Effect of muscle stimulation intensity on the heterogeneous function of regions within an architecturally complex muscle. Journal of Applied Physiology 130 , 941-951. Reed, D.A., and Ross, C.F. (2010). The influence of food material properties on jaw kinematics in the primate, Cebus . Archives of Oral Biology 55 , 946-962. Williams, S.H., Wright, B.W., Truong, V.d., Daubert, C.R., and Vinyard, C.J. (2005). Mechanical properties of foods used in experimental studies of primate masticatory function. Am J Primatol 67 , 329-346. Ahn, A.N., Konow, N., Tijs, C., and Biewener, A.A. (2018). Different segments within vertebrate muscles can operate on different regions of their force–length relationships. Integrative and Comparative Biology 58 , 219-231. Close, R.I. (1972). Dynamic properties of mammalian skeletal muscles. Physiol Rev 52 , 129-197. Buchanan, T.S. (1995). Evidence that maximum muscle stress is not a constant: differences in specific tension in elbow flexors and extensors. Med Eng Phys 17 , 529-536. Weijs, W.A., and Dantuma, R. (1975). Electromyography and mechanics of mastication in the albino rat. Journal of Morphology 146 , 1-33. Holt, N.C., and Azizi, E. (2016). The effect of activation level on muscle function during locomotion: are optimal lengths and velocities always used? Proceedings of the Royal Society of London B: Biological Sciences 283 , 20152832. Schwaner, M.J., Mayfield, D.L., Azizi, E., and Daley, M.A. (2024). Linking in vivo muscle dynamics to in situ force-length and force-velocity reveals that guinea fowl lateral gastrocnemius operates at shorter than optimal lengths. Journal of Experimental Biology. Rack, P.M., and Westbury, D.R. (1969). The effects of length and stimulus rate on tension in the isometric cat soleus muscle. J Physiol 204 , 443-460. Rassier, D.E., MacIntosh, B.R., and Herzog, W. (1999). Length dependence of active force production in skeletal muscle. Journal of Applied Physiology 86 , 1445-1457. Mayfield, D.L., and Holt, N.C. (2025). Does force depression resulting from shortening against series elasticity contribute to the activation dependence of optimum length? 328 , C528-C540. Holt, N.C., and Williams, C.D. (2018). Can Strain Dependent Inhibition of Cross-Bridge Binding Explain Shifts in Optimum Muscle Length? Integrative and Comparative Biology 58 , 174-185. Holt, N.C., and Mayfield, D.L. (2023). Muscle-tendon unit design and tuning for power enhancement, power attenuation, and reduction of metabolic cost. Journal of Biomechanics, 111585. HENNEMAN, E. (1957). Relation between Size of Neurons and Their Susceptibility to Discharge. Science 126 , 1345-1347. Holt, N.C., Wakeling, J.M., and Biewener, A.A. (2014). The effect of fast and slow motor unit activation on whole-muscle mechanical performance: the size principle may not pose a mechanical paradox. Proceedings of the Royal Society B: Biological Sciences 281 , 24695429. Llewellyn, M.E., Thompson, K.R., Deisseroth, K., and Delp, S.L. (2010). Orderly recruitment of motor units under optical control in vivo. Nat Med 16 , 1161-U1144. Hedrick, T.L., Tobalske, B.W., and Biewener, A.A. (2003). How cockatiels ( Nymphicus hollandicus ) modulate pectoralis power output across flight speeds. J Exp Biol 206 , 1363-1378. Konow, N., Collias, A., and Biewener, A.A. (2020). Skeletal muscle shape change in relation to varying foce requirements across locomotor conditions. Frontiers in Physiology 11 , 143. Gabaldon, A.M., Nelson, F.E., and Roberts, T.J. (2008). Relative shortening velocity in locomotor muscles: turkey ankle extensors operate at low V/Vmax. Am J Physiol Regul Integr Comp Physiol 294 , 200-210. Rome, L.C. (1998). Some advances in integrative muscle physiology. Comp Biochem Physiol B Biochem Mol Biol 120 , 51-72. Sinha, A., and Bhargav, A. (2020). Effect of state transition, drying kinetics and moisture content on Young's modulus variation during thermal drying of hygroscopic food materials. Journal of Food Engineering 279 , 109957. Püffel, F., Walthaus, O.K., Kang, V., and Labonte, D. (2023). Biomechanics of cutting: sharpness, wear sensitivity and the scaling of cutting forces in leaf-cutter ant mandibles. 378 , 20220547. Azizi, E., and Roberts, T.J. (2010). Muscle performance during frog jumping: influence of elasticity on muscle operating lengths. Proc Biol Sci 277 , 1523-1530. Langenbach, G.E.J., and van Eijden, T.M.G.J. (2001). Mammalian Feeding Motor Patterns. Amer. Zool. 41 , 1338-1351. Ram, Y., and Ross, C.F. (2020). Jaw Elevator Muscle Coordination during Rhythmic Mastication in Primates: Are Triplets Units of Motor Control? Brain, Behavior and Evolution 95 , 1-14. Peck, C.C., and Hannam, A.G. (2007). Human jaw and muscle modelling. Archives of Oral Biology 52 , 300-304. Shi, J., Curtis, N., Fitton, L.C., O'Higgins, P., and Fagan, M.J. (2012). Developing a musculoskeletal model of the primate skull: Predicting muscle activations, bite force, and joint reaction forces using multibody dynamics analysis and advanced optimisation methods. Journal of Theoretical Biology 310 , 21-30. Cox, P.G., and Watson, P.J. (2023). Masticatory biomechanics of red and grey squirrels (Sciurus vulgaris and Sciurus carolinensis) modelled with multibody dynamics analysis. 10 , 220587. Lautenschlager, S., Gill, P.G., Luo, Z.-X., Fagan, M.J., and Rayfield, E.J. (2018). The role of miniaturization in the evolution of the mammalian jaw and middle ear. Nature 561 , 533-537. Ledogar, J.A., Senck, S., Villmoare, B.A., Smith, A.L., Weber, G.W., Richmond, B.G., Dechow, P.C., Ross, C.F., Grosse, I.R., Wright, B.W., et al. (2022). Mechanical compensation in the evolution of the early hominin feeding apparatus. 289 , 20220711. Basmajian, J., Stecko, G. (1962). A new bipolar electrode for electromyography. Journal of Appplied Physiology 17 , 849. Richards, C.T., and Biewener, A.A. (2007). Modulation of in vivo muscle power output during swimming in the African clawed frog (Xenopus laevis). J Exp Biol 210 , 3147-3159. Afoakwa, E.O., Paterson, A., Fowler, M., and Vieira, J. (2009). Microstructure and mechanical properties related to particle size distribution and composition in dark chocolate. International Journal of Food Science & Technology 44 , 111-119. Vinyard, C.J., and Taylor, A.B. (2010). A Preliminary Analysis of the Relationship Between Jaw-Muscle Architecture and Jaw-Muscle Electromyography During Chewing Across Primates. The Anatomical Record: Advances in Integrative Anatomy and Evolutionary Biology 293 , 572-582. Otten, E. (1987). A Myocybernetic Model of the Jaw System of the Rat. Journal of Neuroscience Methods 21 , 287-302. DeLap, S., Tran, V., and Konow, N. (2024). DeepXROMM. A set of user-friendly tools for DeepLabCut-XMALab marker tracking. (https://pypi.org/project/deepxromm/). Marsh, R.L. (2016). Speed of sound in muscle for use in sonomicrometry. Journal of Biomechanics. Eng, C.M., Konow, N., Tijs, C., Holt, N.C., and Biewener, A.A. (2019). In vivo force–length and activation dynamics of two distal rat hindlimb muscles in relation to gait and grade. The Journal of Experimental Biology 222 , jeb205559. Tables Table 1. Summary data from in situ experiment on rat superficial masseter Rat Stim (V) F O (N) (submax) F O (N) (twitch) L O (mm) (submax) L O (mm) (twitch) Tension (submax) 1 10 6.5* - 12.8 - 10.5 2 9 6.7* - 8.7 - 14.1 3 10 5.7 0.4 6.1 6.7 12.8 4 11 6.6 # - - - 12.3 5 6 3.5 0.7 8.5 9.6 6.9 6 10 3.7 # - - - 7.3 7 6 6.2 1.9 5.2 6.2 8.1 8 12 6.2 1.5 3.1 3.4 11.1 9 8 7.3 2.2 5.8 6.1 9.0 *Twitch FL not measured. # FL curves from motor position, not fascicle length. Additional Declarations No competing interests reported. Supplementary Files SupplementaryMaterials.docx Cite Share Download PDF Status: Published Journal Publication published 30 Oct, 2025 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Revision requested 22 Sep, 2025 Editor assigned by journal 22 Sep, 2025 Reviews received at journal 18 Sep, 2025 Reviewers agreed at journal 18 Sep, 2025 Reviews received at journal 15 Sep, 2025 Reviewers agreed at journal 15 Sep, 2025 Reviewers invited by journal 12 Sep, 2025 Submission checks completed at journal 11 Sep, 2025 First submitted to journal 10 Sep, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-6739430","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":513587054,"identity":"6f365568-2469-47c7-8bd8-a9b95aee0ca0","order_by":0,"name":"Nicolai Konow","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA6klEQVRIie3QvQrCMBSG4U8Kcam6BryJlEBQ9E5cKg4udncQLQid/FkdxMtwVgK6iK6CS4urg5OTiElAx9TRIe/Q9sB5CCngcv1hTKpHCgq20WMfKKsXtZGaJqHaCWI9HgCSR5obQwCup0LyC9l7WRoOahDF6fZaWskWgbe9+BYiJOEs3FE0/WOHl9YySkA6DTvxCQ2JIrQnqpEhvqjaCDfkRSEMWWpSeVhJoEk7oeCGxOYUkkPUXdoz9ZMnB86fu26UeITXlxbCTjJL748h2H4SZItBI5oXx9n5ZiGfRvH30/th3eVyuVz23oqURO3uzNOWAAAAAElFTkSuQmCC","orcid":"","institution":"University of Massachusetts Lowell","correspondingAuthor":true,"prefix":"","firstName":"Nicolai","middleName":"","lastName":"Konow","suffix":""},{"id":513587055,"identity":"c5ddf128-df6e-4384-9b52-d7002246647c","order_by":1,"name":"Brandon Reder","email":"","orcid":"","institution":"University of Massachusetts Lowell","correspondingAuthor":false,"prefix":"","firstName":"Brandon","middleName":"","lastName":"Reder","suffix":""},{"id":513587056,"identity":"9a59f8cd-5cf2-43a9-9e36-124b9fa4c084","order_by":2,"name":"Daniel Bartlett","email":"","orcid":"","institution":"University of Massachusetts Lowell","correspondingAuthor":false,"prefix":"","firstName":"Daniel","middleName":"","lastName":"Bartlett","suffix":""},{"id":513587057,"identity":"0f1b1b64-05e3-49da-bdca-291e45d109e1","order_by":3,"name":"Devin Jenness","email":"","orcid":"","institution":"University of Massachusetts Lowell","correspondingAuthor":false,"prefix":"","firstName":"Devin","middleName":"","lastName":"Jenness","suffix":""},{"id":513587058,"identity":"9afc2c75-fb23-4a90-b0c9-a7b5eb0a94fc","order_by":4,"name":"Trushti Patel","email":"","orcid":"","institution":"University of Massachusetts Lowell","correspondingAuthor":false,"prefix":"","firstName":"Trushti","middleName":"","lastName":"Patel","suffix":""},{"id":513587059,"identity":"b0202197-6f7e-42f8-b220-9576c7e49d02","order_by":5,"name":"Jeffrey R. Moore","email":"","orcid":"","institution":"University of Massachusetts Lowell","correspondingAuthor":false,"prefix":"","firstName":"Jeffrey","middleName":"R.","lastName":"Moore","suffix":""},{"id":513587062,"identity":"e499a270-ceb4-45cf-9060-86b7d1801563","order_by":6,"name":"Robert J. Brocklehurst","email":"","orcid":"","institution":"University of Massachusetts Lowell","correspondingAuthor":false,"prefix":"","firstName":"Robert","middleName":"J.","lastName":"Brocklehurst","suffix":""}],"badges":[],"createdAt":"2025-05-24 13:53:15","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-6739430/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-6739430/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-025-21953-z","type":"published","date":"2025-10-30T15:58:13+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":91153112,"identity":"564b5cb5-c271-4204-8db5-d4d47063e081","added_by":"auto","created_at":"2025-09-12 07:29:31","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":158230,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eInstrumentation and experimental protocols. \u003c/strong\u003e(A) The unipennate rat superficial masseter (grey is tendon, red is fascicles), instrumented with EMG to measure muscle activation (â), sonomicrometry to measure length-change of muscle fascicles (ε), and a custom miniature strain-gauge buckle to measure muscle force ( \u003cem\u003ef \u003c/em\u003e). (B) Representative \u003cem\u003ein vivo \u003c/em\u003emeasurements of muscle activation, length, and force for incisor biting on food with varying hardness (color). Grey bars delineate biting power strokes. Note how starting length (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eMAX\u003c/sub\u003e) can exceed activation length (\u003cem\u003eL\u003c/em\u003e @\u003csub\u003eâON\u003c/sub\u003e) due to earlier activation of other adductor muscles. (C) The muscle preparations tested \u003cem\u003ein vivo \u003c/em\u003ewere also tested \u003cem\u003ein situ\u003c/em\u003e using an ergometry procedure to measure muscle FL with needle stimulation (\u003cstrong\u003e⚡\u003c/strong\u003e, 250ms for submax tetanic, and 2ms for twitch). (D) Representative contractions used to construct the muscle’s FL relationship. Grey bars indicate stimulation intervals.\u003c/p\u003e","description":"","filename":"Picture1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-6739430/v1/f84e39e030534283a4fb5512.jpg"},{"id":91153436,"identity":"ac3ba7f2-7f63-4d4b-bfb3-f5e00f494ece","added_by":"auto","created_at":"2025-09-12 07:37:31","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":114424,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eInfluence of food hardness on muscle recruitment and force production. \u003c/strong\u003eDuring rat incisor biting, greater food hardness increases superficial masseter (A) electromyographic recruitment, and (B) force production. (C) A strong relationship between EMG intensity and force production exists across rats (symbols) and food hardnesses (colors) (For Linear Least-Squares regression results, see methods. LCL, UCL, Lower and Upper Confidence Limits; LPL, UPL, Lower and Upper Probability Limits; ME, Model Estimate).\u003c/p\u003e","description":"","filename":"Picture2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-6739430/v1/9ee34542ae6f3184e1aa03e3.jpg"},{"id":91153111,"identity":"aefd6e38-5d50-4b55-b4f3-5891d90c45d8","added_by":"auto","created_at":"2025-09-12 07:29:31","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":97698,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eFascicle in rat SM operate at short relative lengths (A,B) and strains (C) during incisor-biting. \u003c/strong\u003eFascicle relative lengths (\u003cem\u003eL\u003c/em\u003e/\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) at activation onset (A) and coincident with peak force (B) increase with food hardness (color scheme as in Figure 2). Fascicle strains (Ԑ) averaging 5-7 % across food types are consistent with the pennate architecture of SM, and with operation beyond the FL plateau [5].\u003c/p\u003e","description":"","filename":"Picture3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-6739430/v1/e05730b25c7d93de0165a774.jpg"},{"id":91153113,"identity":"00671f29-7498-48a5-98bc-67ea01ca8238","added_by":"auto","created_at":"2025-09-12 07:29:31","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":40720,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eIn rat SM, the plateau of the relationship between relative force (\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eF\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e/\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eF\u003c/strong\u003e\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e\u003cstrong\u003e) and strain (%\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eL\u003c/strong\u003e\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e\u003cstrong\u003e) is left-shifted with increased muscle activation\u003c/strong\u003e. Submaximal tetanic \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e (dashed vertical line to black filled symbols) is left-shifted by 11.3 ± 2.9% from twitch \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO \u003c/sub\u003e(dashed vertical line to grey filled symbols), consistent with measurements from limb muscles [26, 27, 34, 36]. Open grey symbols show the passive FL relationship. The shaded area highlights how the operating length range of twitch-activated muscle on the ascending twitch FL limb is equivalent to the operating length range of submaximally-activated muscle on the descending submaximal FL limb.\u003c/p\u003e","description":"","filename":"Picture4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-6739430/v1/83bf2e14be7b63059f6db047.jpg"},{"id":91153108,"identity":"76699fe3-d499-4941-b0e6-4b2f3acca1ae","added_by":"auto","created_at":"2025-09-12 07:29:31","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":113633,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003ePlateau and descending FL limb operation of rat SM during incisor biting.\u003c/strong\u003e Superimposition of \u003cem\u003ein vivo\u003c/em\u003e stress-strain operation of rat SM (colored symbols are data from individual bites) onto the passive (grey) and active (black) stress-strain relationships obtained \u003cem\u003ein situ\u003c/em\u003e under submaximal tetanic stimulation. Bite data are peak stress (σ) plotted against fascicle strain (ε) at peak force production and occlusion, typically being the shortest fascicle length reached during biting (Fig. 3B). Note how bites on the hardest food type (chow, E = 50.4 MPa, orange) occur at long strains on the descending FL limb, whereas bites on the softest food type tested (raisin, E = 0.2 MPa, purple) involve low stress and occur at FL plateau lengths.\u003c/p\u003e","description":"","filename":"Picture5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-6739430/v1/1802791db9d4edc092e4e421.jpg"},{"id":95041285,"identity":"eaea22be-132f-48df-95a0-a0eb03c9bed2","added_by":"auto","created_at":"2025-11-03 16:11:06","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1394429,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-6739430/v1/c4181acf-efb1-441d-9ca5-71212e840955.pdf"},{"id":91154357,"identity":"5368db09-773c-4e5e-ab45-10e03d6862d1","added_by":"auto","created_at":"2025-09-12 07:45:31","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":16833,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryMaterials.docx","url":"https://assets-eu.researchsquare.com/files/rs-6739430/v1/e53b26ba7af07ff74fcb8ad4.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Rat superficial masseter operates at long, unstable lengths during biting","fulltext":[{"header":"Introduction","content":"\u003cp\u003eSkeletal movements in vertebrates are driven by muscle contractions, the force generation of which is limited by the overlap of actin and myosin filaments in the sarcomere, in accordance with the force-length (FL) relationship [1-4]. An emerging consensus has been that across scales of biological organization, sarcomeres, fibers, and the fascicles (bundles of fibers) they compose tend to operate at optimal lengths (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e), centered on the FL curve plateau, where skeletal muscle generates peak force (\u003cem\u003eF\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) [5]. However, under scenarios where stable muscle function is priority, fascicle operating length ranges may be truncated to shorter lengths of the ascending FL limb [6-8]. At such short lengths, additional shortening reduces muscle force as myofilament overlap increases [2, 9]. Conversely, muscle may also operate at long lengths, on the descending FL limb [8]. Here, muscle shortening up the descending FL limb increases force, but also involves risks of sarcomere disruption due to the reduced myofilament overlap [10]. Considering the risks of myofilament disruption, muscles have been thought to only operate at long lengths if they are unlikely to undergo abrupt lengthening. The prevailing consensus on muscle operating length ranges comes from studies of limb and cardiac muscle. Therefore, it remains unclear if the consensus about muscle operating lengths transfers directly to jaw muscle operation or not. Determining whether jaw muscles operate at short lengths, or at long and potentially unstable lengths would significantly expand our understanding of skeletal muscle function. This would also have broad implications for craniofacial evolution, muscle mechanics, and control, and for both healthy and pathological function of the craniofacial musculoskeletal system.\u003c/p\u003e\n\u003cp\u003eJaw muscles studies have suggested that force generation during biting at large gapes involve whole muscle and fascicle operating lengths that are short and centered on the ascending FL limb [11-16]. This mode of operation is consistent with the idea that jaw muscles have relatively large moment arms, which result in large muscle length-changes [17] that inherently could render sarcomeres mechanically unstable. However, jaw muscles rarely experience abrupt lengthening (i.e., eccentric contractions, except in cases of live prey struggling). This musculature would thus be expected to operate at FL plateau-lengths, because FL effects associated with muscle shortening as gape approaches occlusion would reduce force production dramatically [9]. If the SM operates at short ascending limb lengths, force loss during gape closing would be particularly challenging for incisor biting at large gape, as used by many mammals [18] including rodents. However, approaches in earlier jaw muscle studies, which include force measurements from the entire musculature at the bite point [13], estimations of muscle length from skeletal displacement, and elicitation of twitch contractions to measure FL from individual muscles [11, 14, 16], are unlikely to accurately reflect the patterns of individual muscle activation, contraction, and force generation used during feeding.\u003c/p\u003e\n\u003cp\u003eMeasurements of muscle activation, fascicle strains, and time-varying force production are critical for reliable interpretations of muscle function \u003cem\u003ein vivo\u003c/em\u003e, for several reasons: (i) The length trajectory of muscle fascicles often differs considerably from the entire muscle-tendon complex, due to the elastic action of compliant tendinous elements [19, 20], and architectural gearing effects caused by fascicles rotating around their tendinous attachments [21-23]. (ii) In most skeletal systems, including jaws, force measured at the contact point can be generated by several muscles whose individual contributions cannot readily be disentangled from time-varying bite-point force [24]. (iii) Muscles are typically activated in a set sequence and at intensities corresponding to the mechanical task. Limb muscle studies have shown that muscle activation intensity impacts the operating length-range of its fascicles [25]. Tetanic activation results in a left-shift of the FL plateau by 10-15% relative to the \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e of a FL curve generated with twitch contractions [26, 27]. However, a general lack of \u003cem\u003ein vivo\u003c/em\u003e measurements of muscle activation, fascicle strains, and time-varying force generation from individual muscles of the jaw means that we do not know if they undergo the same activation-dependent operating length-shift. It also remains unknown what the FL effect on jaw muscle operation might be. We hypothesize that conclusions drawn in earlier jaw studies, suggesting that muscles of the jaw operate at relatively short lengths [11, 14, 16], were influenced by the twitch activation protocol used. We predict that operation at short twitch lengths actually correspond to operation on the FL plateau or longer ascending limb lengths under greater (submaximal tetanic) activation conditions, as necessary for generating sufficient bite force.\u003c/p\u003e\n\u003cp\u003eTo determine if earlier study protocols led to an imprecise consensus about the length-operation range of jaw muscles, we elicited incisor biting in rats feeding on an assortment of natural food items. We used food types that rats willingly accept with material properties, including systematic hardness differences that have been thoroughly characterized, to facilitate comparisons of our data with existing literature [28, 29]. We used direct instrumentation of the rat superficial masseter (SM) to collect direct measurements of time-varying force production from a major jaw closer muscle during feeding, along with electromyographical measurements of muscle activation and recruitment (EMG), and measurements of fascicle length-change dynamics (sonomicrometry). We then tested the same muscle preparations with an \u003cem\u003ein situ\u003c/em\u003e ergometry experiment that measures twitch and submaximal tetanic FL. Finally, we mapped \u003cem\u003ein vivo\u003c/em\u003e fascicle length-operating ranges and forces to the \u003cem\u003ein situ\u003c/em\u003e FL curves [8, 20, 30]. For limb muscles, this approach has proven powerful by determining \u003cem\u003ein vivo\u003c/em\u003e constraints on muscle force and fascicle length, which we have previously only indirectly achieved for jaw muscles [9]. This experimental design allowed us to test the hypothesis that large fascicle strains in SM, associated with incisor biting at large gape, \u0026nbsp;lead to peak force generation at fascicle operating lengths that correspond to the descending limb on the submaximal tetanic FL curve, in contrast to the findings of earlier studies that used twitch contractions.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003eAcross the six rats tested \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003ein biting experiments, we were not always able to collect reliable bites for all food types, and a total of 247 bite recordings were used for analyses (see Supplemental Table 3).\u003c/p\u003e\n\u003cp\u003eWe found a strong relationship between recruitment of rat SM and its force production across incisor bites on foods with varying hardness. Food hardness correlated strongly with muscle recruitment intensity (GLM: F\u003csub\u003e2,236\u003c/sub\u003e = 10.8; p \u0026lt; 0.0001; η\u003csup\u003e2\u003c/sup\u003e = 0.32), measured as the integrated area under the rectified EMG signal (Fig. 2A). Harder food was associated with greater muscle force production (Fig. 2B), again with a statistically significant relationship (GLM: F\u003csub\u003e2,161\u003c/sub\u003e = 4.8; p \u0026lt; 0.01; η\u003csup\u003e2\u003c/sup\u003e = 0.07). Whereas EMG activity and force production generally increased with food hardness, this relationship plateaued at levels of hardness above 20 MPa. Rat chow, nominally the hardest food tested (50.44 MPa; [29]), elicited lower-than-expected EMG and force production compared to predictions based on muscle force production on softer foods.\u003c/p\u003e\n\u003cp\u003eStress developed \u003cem\u003ein vivo\u003c/em\u003e by the rat SM increased considerably with food hardness and averaged 4.1 ± 3.9 N/cm\u003csup\u003e2\u0026nbsp;\u003c/sup\u003e(mean ± S.D.). Greater stresses were consistently recorded for the two hardest food items (almond and chow) and peaked for almond\u0026nbsp;at 10.6 N/cm\u003csup\u003e2\u003c/sup\u003e, with median stress for chow bites (7.2\u0026nbsp;± 3.9 N/cm\u003csup\u003e2\u003c/sup\u003e) being marginally lower than for almond bites (8.6\u0026nbsp;±\u0026nbsp;5.3 N/cm\u003csup\u003e2\u003c/sup\u003e) (F\u003csub\u003e2,125\u003c/sub\u003e = 4.966; p \u0026lt;0.05).\u0026nbsp;Based on an average muscle PCSA of 0.51 ± 0.12 cm\u003csup\u003e2\u003c/sup\u003e (Supplemental Table 1), and specific tension data for vertebrate skeletal muscle [31, 32], the rat SM operated at submaximal recruitment during incisor biting. A Linear Least-Squares regression of EMG recruitment against force production indicated a strong EMG-force relationship (F\u003csub\u003e1,161\u003c/sub\u003e = 516.9; p \u0026lt; 0.0001), despite identifiable variation across subjects and food hardness (Fig. 2C).\u003c/p\u003e\n\u003cp\u003eThe operating length-ranges of fascicles in the rat SM were strongly correlated with food hardness. Peak relative lengths (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eMAX\u003c/sub\u003e) typically occurred at bite onset (Fig. 1B) and preceded as well as exceeded relative fascicle length at activation (\u003cem\u003eL\u0026nbsp;\u003c/em\u003e@â\u003csub\u003eON\u003c/sub\u003e), likely due to earlier activation and action of different adductor muscles in the rat feeding system [33]. Increased food hardness elicited longer relative fascicle lengths at muscle activation onset (\u003cem\u003eL\u003c/em\u003e/\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u0026nbsp;\u003c/sub\u003e@â\u003csub\u003eON\u003c/sub\u003e) (GLM: F\u003csub\u003e2,237\u003c/sub\u003e = 52.1; p \u0026lt; 0.0001; η\u003csup\u003e2\u003c/sup\u003e = 0.98) (Fig. 3A) and at occlusion (\u003cem\u003eL\u003c/em\u003e/\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e @F\u003csub\u003eMAX\u003c/sub\u003e) (F\u003csub\u003e2,237\u003c/sub\u003e = 74.4; p \u0026lt; 0.0001; η\u003csup\u003e2\u003c/sup\u003e = 0.68) (Fig. 3B). The strain range for contracting fascicles varied considerably (Fig. 3C) and averaged 10.5 ± 3.2%, with a statistically significant reduction in strain with increasing food hardness (F\u003csub\u003e2,237\u003c/sub\u003e = 24.8; p \u0026lt; 0.0001; η\u003csup\u003e2\u003c/sup\u003e = 0.69).\u003c/p\u003e\n\u003cp\u003eOur \u003cem\u003ein situ\u003c/em\u003e preparations of rat superficial masseters achieved peak tension development under twitch activation at optimal length (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) reaching 12.3 ± 4.3% relative to tetanic activation (\u003cem\u003eF\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e). Tetanic tension at \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e averaged 10.1 ± 3.3 N/cm\u003csup\u003e2\u003c/sup\u003e, suggesting that our needle stimulation protocol elicited submaximal and not maximal tetanic contractions. Comparisons of \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e for twitch and submaximal tetanic FL curves showed that the submaximal tetanic plateau was left-shifted by 11.3 ± 2.9% with respect to the twitch plateau (Fig. 4).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSuperimposition of \u003cem\u003ein situ\u003c/em\u003e and \u003cem\u003ein vivo\u003c/em\u003e data from rat SM (Fig. 5) revealed incisor biting to involve fascicle operating strain-ranges spanning from in excess of 25% over\u0026nbsp;\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e, on the descending FL limb (i.e., high\u0026nbsp;ε) across the FL plateau to approximately 7% down the slow ascending FL limb (i.e., low ε). Both hard and soft food types elicited bites involving high fascicle peak strains (ε\u003csub\u003e\u0026nbsp;MAX\u003c/sub\u003e). However, fascicle operation at low strains were particularly evident for soft food types and fascicle operation at high strains were particularly common for hard food types, which also elicited the greatest tension generation as the jaw approached occlusion.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eOur \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003estudy of the rat superficial masseter during incisorbiting generated the first available ensemble measurements of electromyographical activation, fascicle strain, and force from a jaw muscle during feeding. We mapped the \u003cem\u003ein vivo\u003c/em\u003e stress and operating strain ranges of the same muscle onto its submaximal force-length (FL) relationship, as measured from the same muscle preparations \u003cem\u003ein situ.\u003c/em\u003e Combined, these data revealed that as activation increases, the FL plateau of rat superficial masseter (SM) shifts leftward to a similar extent as previously reported for various limb muscles. Taken together, our data reveal that the rat SM operates \u003cem\u003ein vivo\u003c/em\u003e across lengths that range from the upper FL ascending limb, across the plateau, and down the descending limb. This finding runs counter to the well-established paradigm that jaw muscles operate a short, stable lengths [11-16]. Our findings further develop our understanding of how skeletal muscle exploits the FL relationship \u003cem\u003ein vivo\u003c/em\u003e to maximize performance for specific tasks, in this case biting. These results highlight the role of muscle operating length in maintaining musculoskeletal stability during routine motions. Our finding has broad implications, ranging from reconstruction of craniofacial evolution in vertebrates to the prevention and treatment of temporomandibular joint disorders (TMD) in humans.\u003c/p\u003e\n\u003cp\u003eMuscle operating lengths depend on the intensity of muscle activation and recruitment. For the rat SM, we discovered a left-shift of the plateau of the muscle’s FL curve averaging\u0026nbsp;11.3 ± 2.9% \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e as we increased stimulation voltage \u003cem\u003ein situ,\u003c/em\u003e from twitch to submaximal tetanic stimulation. This shift in optimal length (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) is consistent with earlier experiments that varied stimulation voltage to elicit twitch, submaximal tetanic, and maximal tetanic contractions in limb muscles, including toad plantaris (\u0026gt;15%\u003cem\u003e\u0026nbsp;L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e,[34]), Guinea Fowl lateral gastrocnemius (20%\u0026nbsp;\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e,\u0026nbsp;[35]). rat medial gastrocnemius (12.6-11.2%\u003cem\u003e\u0026nbsp;L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e, across muscle compartments,\u0026nbsp;[27]) and mouse soleus and tibialis anterior (both 8.2%\u003cem\u003e\u0026nbsp;L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e,\u0026nbsp;[26]). For cat soleus, a similar shift wad demonstrated by varying stimulation frequency (25%\u003cem\u003e\u0026nbsp;L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e,\u0026nbsp;[36]). Our discovery of a plateau-shift for a jaw muscle that is comparable to measurements from limb muscles with a diversity of fascicle architectures, series elastic element prominences, and mechanical functions suggest that this trait is generalizable across striated skeletal muscle. The mechanistic explanation for the plateau-shift remains incomplete as both strain-dependent calcium sensitivity-changes affecting cross-bridge cycling\u0026nbsp;[25, 36-38]\u0026nbsp;and series elastic compliance\u0026nbsp;[38-40]\u0026nbsp;have been implicated. Regardless, our findings from rat SM suggest that studies using twitch activation to examine FL behavior in jaw muscles\u0026nbsp;[14, 16]\u0026nbsp;have generated an inaccurate consensus about the lengths across which jaw adductor muscles operate. Our conclusion is strengthened by the consideration that muscles rarely are twitch-activated nor supramaximally activated \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003eand typically operate at submaximal levels of activation and recruitment\u0026nbsp;[27].\u003c/p\u003e\n\u003cp\u003eThe central nervous system regulates muscle activation timing and both the intensity and completeness of muscle recruitment to changes in loading environment and force requirements [41-43], for instance, as elicited by changes in food hardness. In agreement with this principle, we found a strong positive relationship between the recruitment intensity (integrated EMG area) of the rat SM during \u003cem\u003ein vivo\u003c/em\u003e incisor biting power-strokes with the hardness of the food substrate (LLS regression, \u003cem\u003ep\u003c/em\u003e \u0026lt; 0001). A similar strong positive relationship was also established between muscle force output and food hardness, with few exceptions (discussed below). Overall, our findings indicate a generally strong agreement between EMG and force measurements, consistent with recordings from a diversity of limb musculoskeletal systems [44, 45]. This finding suggests that EMG holds great promise for understanding time-varying force production in other jaw muscles, where substrates for force-buckle attachments are unavailable, such as the pterygoid and temporalis.\u003c/p\u003e\n\u003cp\u003eAs food hardness was increased during rat \u003cem\u003ein vivo\u003c/em\u003e biting, we found systematic changes in the operating length ranges of the muscle fascicles, as measured directly using sonomicrometry. As recruitment intensity and force production increased with food hardness, muscle fascicles operated across longer length-ranges but also contracted to relatively shorter lengths (Fig. 3B; GLM; both p \u0026lt; 0.0001). These findings underscore that muscles can exploit FL effects by tuning operating length to generate forces that match task-based requirements, contend with length-based limitations to force generation, and mitigate risks of contracting at long lengths [8, 9].\u003c/p\u003e\n\u003cp\u003eTo better understand how the rat SM exploits its FL relationship during feeding, we superimposed our \u003cem\u003ein vivo\u003c/em\u003e and \u003cem\u003ein situ\u003c/em\u003e measurements. This is a validated and powerful method for rigorously interpreting operating requirements and constraints on skeletal muscles in limbs [20, 35, 46]. Mapping the \u003cem\u003ein vivo\u003c/em\u003e operating length ranges of rat SM onto the muscle’s own FL relationships, as measured \u003cem\u003ein situ\u003c/em\u003e allowed testing of the hypothesis that twitch FL studies inaccurately reflect the muscles’ contractile behavior. In line with earlier twitch studies of jaw muscles [14, 16], we determined that the \u003cem\u003ein vivo\u003c/em\u003e operating length-ranges of the rat SM during incisor biting corresponded to the ascending limb (i.e. short operating lengths) of the twitch FL curve. However, when the \u003cem\u003ein vivo\u003c/em\u003e operating length ranges were mapped to the significantly shorter plateau-lengths of the submaximal FL curve, with comparable stresses, we found that the realized operating lengths were typically distributed across the FL plateau. For bites on hard food (almond, chow) the muscle often operated at relatively high stresses on the descending limb (i.e., long lengths). This discovery is important for several reasons. Limb muscles will often, but not always, operate across the FL plateau where muscle force is maximized [8, 47]. Our data, therefore, suggest general similarities between limb and jaw muscles in this regard (see also [9]).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWe tested natural foods that rats would willingly accept and attempted to size-match food items to avoid food size-driven FL effects in our \u003cem\u003ein vivo\u003c/em\u003e data. Evidence that our attempts to decouple muscle force from gape and fascicle length were not entirely successful is seen for our hardest food type, chow (E = 50.44 MPa [29]). Chow elicited significantly lower muscle recruitment intensity and force production (Figure 2) than our second-hardest food type, almond (E = 19.4 MPa). This result may be influenced by uncontrolled food properties such as moisture content, with which hardness likely co-varies [48]\u0026nbsp;as well as the granular nature of chow that may afford the incisors easier purchase, and also fractures more easily between incisors than between materials testing plates. However, food-driven FL effects were likely not completely avoided as chow pieces invariably were larger than raisins. Consequently, length effects on muscle function, driven by variations of food size\u0026nbsp;[9]\u0026nbsp;likely also play a role because relative fascicle operating lengths were considerably longer for chow than for other food types (Fig. 3A, B). Long fascicle lengths may have prevented greater force generation on chow (Fig. 2A,B), although the associated depression in force likely was somewhat offset by force-velocity effects\u0026nbsp;[46]. For example, muscle operating strains on chow were generally shorter (Fig. 3C) and slower (Fig. 1B) than for other food types. Conversely, our softest food type, raisin (E = 0.22 MPa) elicited greater-than-expected muscle forces (Fig. 1B) and more variable operating lengths (Fig. 2A,B) than other food types. This result may be caused by the greater compliance (i.e., stress to fracture) of raisins, requiring greater biting work than other food types tested (Supplemental Table 2). Together, the anomalous force-responses by SM on our nominally softest and hardest food type respectively may be the reason that the influence of food hardness on SM force production carried the lowest effect size (η\u003csup\u003e2\u003c/sup\u003e = 0.07)\u0026nbsp;of our statistical comparisons.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eFuture work building on ours could use manufactured food items that permit more careful control and parameterization of shape, size, and hardness [9, 49] than natural foods do. Analyzes could include single bite interactions on spherical food with controlled variations of size and/or hardness. Such analyses would also benefit from videofluoroscopic imaging [9] to precisely determine where the incisors make purchase on the food item and jaw gape. They should also incorporate skeletal and muscle moment arm analyses [17], especially considering that rodent jaw joint anatomy allows propalinal (caudal-to-rostral) jaw translation during gape closing, which likely affects muscle operating length independent from vertical gape movements.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eFoods used in our study are unlikely to have challenged the rat jaw apparatus, which for instance is capable of fracturing pistachio shell (E = 1.2 GPa; Pers. Obs.). Accordingly, forces measured \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003ewere unlikely representative of the peak voluntary contractile capacity of the rat SM, which may be realized not only during feeding but also during defensive biting. Here, the maximal FL curve would be relevant, and our specific tension calculations suggest that the maximal FL plateau would be left-shifted by 10% from submaximal \u003cem\u003eL\u003c/em\u003eo, placing peak bite force events on the deep descending limb. We were unable to elicit supramaximal contractions \u003cem\u003ein situ\u003c/em\u003e due to the surgical infeasibility of accessing the nerve of masseter for mounting of a cuff electrode that would facilitate supramaximal stimulation of the rat SM. However, given the well-matched peak tension recorded \u003cem\u003ein vivo\u003c/em\u003e during biting and elicited in our \u003cem\u003ein situ\u003c/em\u003e tetanic contractions, the methodological limitations of needle electrode stimulation are unlikely to affect the conclusions of our study. Despite the good agreement between the specific tension values we recorded \u003cem\u003ein vivo\u003c/em\u003e and the peak tension elicited \u003cem\u003ein situ\u003c/em\u003e, future studies should explore alternative stimulation approaches [14] to elicit supramaximal stimulation and generate \u003cem\u003ein situ\u003c/em\u003e FL data that match supramaximal behaviors.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe discovery that the rat SM operated at long fascicle and sarcomere lengths to generate sufficient force to fracture hard food is an interesting\u0026nbsp;finding from the perspectives of muscle as well as jaw joint stability and safety. Compared to short fascicle and sarcomere length operation, where muscle is stable [6], large, hard food items drove the muscle to develop high forces at long lengths where inherent sarcomere instability due to weak actomyosin overlap [2] might\u0026nbsp;cause muscle damage [10]. The dangers of the descending FL limb have been highlighted for limb muscle [50] but may potentially be greater for jaw muscles due to their larger moment arms [17]. If food item slippage was to translate the jaw during high force production at large gape, the resulting stretch of the muscle might compromise its ability to counter dislocating torques at the temporomandibular joint (TMJ). Other potential deleterious results from such an event include damage to the dentition [16]. Jaw muscles in mammals tend to alternate between powering and balancing the jaw during feeding [24, 51, 52]. Accordingly, muscle instability might trigger jaw joint instability and our data suggest that biting on large and hard objects could be an overlooked trigger of dysphagia including temporomandibular joint disorders (TMD). We anticipate future FL operation analyses of mammalian jaw muscles to inform research on temporomandibular joint (TMJ) function and triggers of its dysfunction.\u003c/p\u003e\n\u003cp\u003eOur data on the effects of food types with different material properties echo previous work by paleontologists and anthropologists seeking to reconstruct foraging and feeding in our mammalian and human forebearers [28-29]. However, such studies have previously lacked the analytical power that single muscle force measurements provide [20,51,55]. The force-length effects highlighted by our data are especially relevant for computational modeling of the craniofacial system, where the force-length relationships of individual muscles tend to be greatly simplified and certainly do not incorporate factors such as activation-dependent shifts of optimal length [53-55]. Computational models are a common tool given the challenging nature of experimental work on primates (including humans) and are the only tool available for studying craniofacial evolution in the fossil record. Several key transitions in vertebrate evolution involve changes to the mandible and jaw joint, where continued functioning and stability of the system would be paramount e.g., mandibular reduction in modern humans, and the transformation of the post-dentary bones into the mammalian middle-ear ossicles [56, 57]. Our experimental data demonstrate a clear need to fully understand jaw biomechanical systems to inform modeling approaches aiming to accurately reconstruct key functional transitions across major clades.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003eSprague-Dawley rats (\u003cem\u003eN\u003c/em\u003e = 10, 50/50 M/F) were purchased from Charles River and same-sex pair-housed under a 12 H light : 12 H dark cycle at the Concord Field Station (Harvard University). All methods were carried out in accordance with relevant guidelines and regulations. Methods are reported in accordance with ARRIVE guidelines. Both husbandry and experimental procedures adhered to protocols approved by the Institutional Animal Care and Use Committees at Harvard University (FAS, 20-09_04) and UMass Lowell (23-03). Rats were provisioned with water and rodent chow\u003cem\u003e\u0026nbsp;ad libitum\u003c/em\u003e and at the time of experiments, body masses were 440 ± 113 g (mean ± S.D.) and ages were 13.5 ± 1.8 weeks (mean ± S.D.) (Supplemental Data Table 1). Power analysis based on muscle activation, strain, and force data from our pilot study suggested that \u003cem\u003eN\u0026nbsp;\u003c/em\u003e=5 would yield a Delta (mean difference over individual variation) for two-factored GLM designs of 1.8-2.5 and analysis power ranging between 87% and 96% at α = 0.01 (Systat v. 12).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eInstrumentation surgery\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eRats (\u003cem\u003eN\u003c/em\u003e = 6) were induced in a chamber on isoflurane (4-4.5%) via oxygen (2.0L/min) and transferred to a low-profile mask for anesthetic maintenance at 1.5-2.5% isoflurane via oxygen (0.8 L/min). Core temperature was monitored via a rectal thermocouple and adjusted using an infrared lamp and a controllable heat pad. Breathing rate and core temperature were continuously monitored and recorded at 10 min intervals throughout each procedure. After application of ophthalmic ointment and injection of 2ml sterile saline (SQ), incision areas were shaved and sterilized using three repeated alternating scrubs of chlorhexidine and betadine. The subject was draped with semitransparent plastic (Press’n Seal) to allow visual monitoring of breathing during the sterile implantation procedure that followed.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTwo skin incisions were made, one along the left lateral-to-ventral chin and one along the mid-dorsum of the neck to route wires subcutaneously through a sterile plastic straw, from a skin port (Microtech, Bronwyn PA) attached at the dorsal incision, to the chin incision, over the SM (Fig. 1A). A pair of sonomicrometry crystals (1 mm, Sonometrics Corporation, London ON) were implanted into the SM muscle along the length of a fascicle. One crystal was implanted at the myotendinous junction, and one at the mandibular origin. The distance between the crystals were 12 ± 3 mm. Crystals were anchored to the epimysium using 6-0 silk and oriented in the tissue to generate the strongest possible transducer pulse (monitored during surgery via an oscilloscope). To measure EMG, a bipolar, offset-hooked electrode [58] (California fine wire, Grover beach, CA) was implanted using a 25G hypodermic needle mid-belly into the muscle, along the same fascicle as the sonomicrometry. The EMG electrode was also sutured to the epimysium using 6-0 silk. A custom-made miniature strain buckle (FLK-1–11, Tokyo Sokki Kenkyujo) [59] was tied to the tendon of SM using 3-0 silk. The rat and its instrumentation was grounded using a large-gauge blunt monopolar steel electrode located under the skin port chassis. Both skin incisions were closed using 3-0 Vicryl. Flunixin Melgumine (30mg/kg, SQ) was administered as pain relief and analgesics intraoperatively, and at 12 hour intervals throughout recovery. Once awakened from surgery, subjects were fed soft (raisin) food within the hour to ensure recovery to normal behavior, and then fasted during recovery (overnight, 12-16 hours) to elicit optimal feeding performance during the \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003efeeding experiment.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eIn vivo\u0026nbsp;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003efeeding experiment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eOur \u003cem\u003ein vivo\u003c/em\u003e experiment was designed to test bite performance across a range of food types with known hardness (Range of Young’s Modulus; 0.22 – 50.44MPa; [28, 29, 60] (Supplemental data Table 2) and near-equivalent size (approx. 1.5 cm diameter). We avoided mechanical alteration of food size due to concerns that this would affect the opportunities for the subject’s incisors to make purchase on the food during biting. During experiments\u003cem\u003e\u0026nbsp;in vivo\u003c/em\u003e, the subject was placed in a transparent acrylic enclosure (25 x 25 x 25 cm). The data collection cable was connected to the skin port and the recording equipment and freely suspended above the subject’s neck using nylon string. The rat was presented with food items (Supplemental data Table 2) in a pseudorandomized sequence from a sturdily attached, spring-loaded metal (alligator) clip. Clipping of food was chosen to facilitate video recordings of the bite without the forelimbs obscuring the view of the incisors. Food clipping also prevented the rat from generating bite torque using the upper incisors and the forelimbs, which potentially would decouple jaw adductor muscle function from the biting action. A single bite series was allowed per food item to minimize FL effects in the muscle data caused by food size reduction due to repetitive biting. Figure 1B shows individual bites on different foods.\u003c/p\u003e\n\u003cp\u003eTo collect gape cycle data, we recorded high speed video of all bites in a quasi-ventral view (125fps, 1024x1024 pixels, Fastcam SA3, Photron, San Diego, CA). Sonomicrometry data were acquired at 520 Hz (transmit pulse: 406.25 ns; inhibit delay: 1.54 mm; range-gain: 25mm) using a digital sonomicrometer and a PC laptop running SonoLab (Sonometrics Inc., London, ON). Muscle EMG data were recorded at 100x amplification (Grass P511) with a 10 kHz high-cut and a 60Hz notch filter engaged. Strain gauge measurements were balanced and amplified (200x; Vishay Measurements Group 2120, Raleigh, NC) and all time-varying transducer data were digitized at 4kHz via a Powerlab 16/20 onto a PC laptop running LabChart 8.0 (ADInstruments). The animal ground was connected to the equipment ground via the sonomicrometer cable shield. A manually triggered square-wave TTL signal terminated video recordings and permitted resynchronization of time-varying transducer and gape cycle data.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eIn situ\u0026nbsp;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003eFL ergometry\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eDue to \u003cem\u003ein situ\u0026nbsp;\u003c/em\u003epreparation and sonomicrometry transduction difficulties for two of our \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003erats, we used additional rats (\u003cem\u003eN\u003c/em\u003e = 4) to generate FL curves \u003cem\u003ein situ\u003c/em\u003e. For these rats, sonomicrometry implantations were carried out during the ergometry preparation and EMG and force buckles were not mounted.\u003c/p\u003e\n\u003cp\u003eDirectly following the \u003cem\u003ein vivo\u003c/em\u003e data collection, each subject was re-induced and maintained on isoflurane, including vitals controls and monitoring as described above. The skin incision on the chin was carefully reopened, and the tendon of the superficial masseter was gently dissected free from its insertion on the rostrum and connected to the lever-arm of a servomotor (Aurora Scientific 300C-LR) using a self-tightening knot of 1-0 silk and a dab pf cyanoacrylate gel (Duro) (Fig. 1C). Two pairs of stimulation electrodes were implanted and secured in place using 6-0 silk. Care was taken to place equivalent poles of each pair at the rostral (+) and caudal (-) ends of the muscle respectively. The insertion orientation of each electrode was normal to the fascicle plane, to avoid alterations in stimulation and recruitment as fascicles shortened against series elastic compliance. The stimulation electrode leads were flexibly suspended above the muscle preparation to minimize interference by motion during contractions. The left zygoma and ventral mandibular corpus were both firmly clamped to a customized stereotaxic frame onto which the servomotor was mounted. Muscle temperature was maintained as constant as possible using an infrared lamp and a dressing of saline soaked gauze prevented muscle desiccation during the procedure.\u003c/p\u003e\n\u003cp\u003eWe acquired data on muscle stimulation intensity (V), muscle fascicle length (measured via the sonomicrometry crystals, mm), and motor force (N) along with measurements from all implanted transducers at 4kHz via a NIDAQ (USB-6212) to a PC running IgorPro 8.0 (Wavemetrics, Lake Oswego, OR). Optimal stimulation parameters were determined for each preparation using a series of twitches with increasing voltage (6-12 volts depending on preparation; 125 pulses/s, 0.2 ms pulse width; Grass s48). Four to five 300 ms contractions from a passive tension that roughly corresponded to \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e were used remove series compliance and collect time-varying data from the tendon strain-gauge (V) and the motor (N) for calibration purposes (see below) [59]. Video recordings of sample tetanic contractions in two preparations were used to verify that the fixed-end compliance of the preparation was minimal (measured to be 0.7 ± 0.2 mm). Using consecutive twitches, the preparation was shortened to minimal length (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eREF\u003c/sub\u003e), which was measured from the muscle origin on the mandible to the suture knot (i.e., the entire preparation) to the nearest millimeter using calipers. Starting at \u003cem\u003eL\u003c/em\u003e\u003csub\u003eREF\u003c/sub\u003e, we then collected the twitch and submaximal tetanic FL curves for the preparation, and the fascicles, as measured via sonomicrometry (Fig. 1D). Three minutes of rest was allowed between each tetanic stimulation to reduce effects of muscle fatigue. A final tetanic contraction at \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e was used to confirm that the preparation had lost no more than 10% of \u003cem\u003eF\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e during testing, which was confirmed for all preparations. Following the \u003cem\u003ein situ\u0026nbsp;\u003c/em\u003eexperiment, the subject was euthanized (Sodium Pentobarbital, IC). The muscle was excised to record fascicle length, including the segment of the fascicles that were not spanned by the crystals, so that a segment correction could be applied to the \u003cem\u003ein vivo data.\u0026nbsp;\u003c/em\u003eWe also measured resting muscle pennation angle and muscle mass (without the tendon) to calculate muscle physiological cross sectional area (PCSA) [61], which was used to calculate muscle stress (N/cm\u003csup\u003e2\u003c/sup\u003e).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData extraction \u0026amp; analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eConstruction of force-length curves\u003c/em\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWe constructed force-length (FL) curves to determine twitch and submaximal tetanic optimal length (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) from our sonomicrometry measurements of fascicle length-change and peak force (P\u003csub\u003eO\u003c/sub\u003e) from the motor and muscle force buckle for each SM preparation by fitting third order polynomials to the passive and active twitch and tetanic curves in IgorPRO v.8. Force (N) was expressed as stress (N/cm\u003csup\u003e2\u003c/sup\u003e) by dividing with muscle PCSA (Supplemental table 1) and fascicle length (mm) was expressed as strain ((\u003cem\u003eL\u003c/em\u003e-\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e)/\u003cem\u003e\u0026nbsp;L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e). To measure the percent left-shift of tetanic \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u0026nbsp;\u003c/sub\u003erelative to twitch \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e, we used tetanic \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u0026nbsp;\u003c/sub\u003eto strain-normalize twitch FL curves [25, 26]. Passive FL curves were adjusted for shortening during stimulation [62].\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eConditioning of\u0026nbsp;\u003c/em\u003ein vivo\u003cem\u003e\u0026nbsp;measurements\u003c/em\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eGape cycles were digitized for usable bites in our custom confluence of DeepLabCut and XMAlab [63]. Time-varying traces of gape cycles were imported into IgorPro along with CSV data from Powerlab containing activation, length, and force measurements from the implanted transducers. Gape and instrument data were re-synchronized using the TTL square-wave signal. All subsequent signal conditioning and analyses were performed in IgorPro. Muscle EMG signals were filtered as needed using an FIR interpolation, base-line corrected and rectified. Sonomicrometry signals (V) were converted to length (mm) using a synchronized calibration recording (in SonoLab). These signals were then segment-corrected for the fascicle length not spanned by the crystal pair, and adjusted for a time-offset caused by the crystal coating epoxy [64]. The conditioned signal was smoothed using a quintic spline interpolation (s.d. = 0.01 – 0.05) and strain (ε = (\u003cem\u003eL\u003c/em\u003e-\u003cem\u003e\u0026nbsp;L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e)/\u003cem\u003e\u0026nbsp;L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) was calculated using tetanic \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e. To calibrate the tendon buckle signal (V) from \u003cem\u003ein vivo\u003c/em\u003e recordings to force (N), we first regressed the time-varying gauge signal (V) against motor force (N) from the initial \u003cem\u003ein situ\u0026nbsp;\u003c/em\u003etetanic contractions. Linear regressions of these two time-varying measurements for the period of tetanic force rise returned \u003cem\u003er\u003c/em\u003e\u003csup\u003e2\u003c/sup\u003e values ranging from 0.87 to 0.96, which were comparable to prior validations of this technology [45, 59, 65]. The average regression slope (V/N) was then used to calculate time-varying \u003cem\u003ein vivo\u003c/em\u003e muscle-tendon force (N), which was quintic spline interpolated as needed and baseline-corrected to zero in a period of the \u003cem\u003ein vivo\u003c/em\u003e trial where the animal was at rest and no feeding occurred.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eExtraction of in-vivo measurements\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eTime-varying data from the conditioned and calibrated \u003cem\u003ein vivo\u003c/em\u003e measurements of SM activation, strain, and force were graphed in IgorPro. For the period of gape closing, a semi-automated procedure was used to extract the following response variable: (1) Mean EMG (mV), (2) Fascicle starting length at the power stroke onset (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eMAX\u003c/sub\u003e), (3) Fascicle length at EMG onset (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eACT\u003c/sub\u003e), (4) Fascicle length at peak force (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eMIN\u003c/sub\u003e), (5) Fascicle shortening strain (ε=[\u003cem\u003eL- L\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e]/\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) during the power stroke (T\u003csub\u003eACT\u003c/sub\u003e to T\u003csub\u003eMIN\u003c/sub\u003e), (6) Force rise magnitude (N), and (7) Peak force (F\u003csub\u003eMAX\u003c/sub\u003e, N). Data were tabulated per individual (\u003cem\u003eN\u0026nbsp;\u003c/em\u003e = 6), cycle number (to control for effects of food size reductions across repeated bites), and food hardness. Data extraction was performed by two researchers who were blinded to the treatment (food hardness) across trials.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eHypothesis testing and statistics\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eTo test the hypothesis that activation and force production by the muscle fascicles would vary systematically with food hardness, we ran General Linear Models (GLM) in Systat v. 12. These analyses factored EMG (specifically the integrated area under the rectified electromyogram during jaw closing, in mV/s) or muscle force (N) as the response variables, individual as an independent effect, and food hardness (E) as the fixed factor. We interrogated the interaction term (IND*E) for statistical significant. Effect sizes of the interaction terms were calculated as\u0026nbsp;the partial eta squared (η\u003csup\u003e2\u003c/sup\u003e) in SPSS 30.0.\u0026nbsp;\u0026nbsp;We also probed the strength of the EMG-force relationship across individuals and food types using Linear Least-Squares regression in Systat v. 12.\u003c/p\u003e\n\u003cp\u003eTo determine if fascicle relative length (\u003cem\u003eL\u003c/em\u003e/\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) at activation (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eACT\u003c/sub\u003e) and at peak force generation (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eMIN\u003c/sub\u003e), as well as fascicle strain during the incisor biting power stroke (% \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) increased with to food hardness (E), we ran GLM and effect-size analyses using the design described above.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eStress-strain FL curves from the \u003cem\u003ein situ\u0026nbsp;\u003c/em\u003eexperiments were used to determine if optimal length (\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e) for rat SM varies with recruitment intensity, as previously reported for limb muscles. To determine this, we ran a GLM factoring \u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e as the response variable, with stimulation protocol (twitch or submax) as a fixed factor, individual as an independent effect with interrogation of the interaction-term and computation of effect-size as above.\u003c/p\u003e\n\u003cp\u003eTo determine the ranges of stresses and strains that the fascicles of the rat SM operate across during incisor biting, we plotted peak muscle stress against fascicle final strain (at peak force generation, Fig. 2B) across food items with varying hardness. Our \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003efascicle operating strain data were recorded using the same sonomicrometer crystals used to measure fascicle strain \u003cem\u003ein situ\u003c/em\u003e. This enabled us to superimpose \u003cem\u003ein vivo\u0026nbsp;\u003c/em\u003estress-strain data onto the aggregate FL curves, expressed as stress versus strain, as measured \u003cem\u003ein situ\u003c/em\u003e using submaximal stimulation.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eData access statement:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eData analyzed in this study are available from the corresponding author upon reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSupplemental information:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIs available within this submission.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements:\u0026nbsp;\u003c/strong\u003eB. Rimkus, T. Flint, S. Angelli-Nicholls, and K. Katugam-Dechenne, assisted with experiments. G. Crowell, and N. Raymondeau assisted with data extraction from the \u003cem\u003ein vivo\u003c/em\u003e measurements. P. Ramirez assisted with husbandry, and A. A. Biewener provided access to facilities at the Concord Field Station. C. F. Ross and N. J. Gidmark both motivated and offered helpful discussions of the study.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding:\u003c/strong\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThis work was funded by the National Science Foundation (award 2217246).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eConceptualization: N.K., D. J., \u0026amp; T.P, Methodology: N.K. Investigation: N.K., B.R., D.B., D.J., T.P., J.M. \u0026amp; R.B. Data analysis: N.K, D.J., T. P. Writing—original draft: N.K. Writing—review and editing: N.K. B.R., D.B., D.J., T.P., J.M. \u0026amp; R.B.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing influences.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n \u003cli\u003eRamsey, R.W., and Street, S.F. (1940). The isometric length-tension diagram of isolated skeletal muscle fibers of the frog. Journal of Cellular and Comparative Physiology \u003cem\u003e15\u003c/em\u003e, 11-34.\u003c/li\u003e\n \u003cli\u003eGordon, A.M., Huxley, A.F., and Julian, F.J. (1966). The variation in isometric tension with sarcomere length in vertebrate muscle fibres. J Physiol \u003cem\u003e184\u003c/em\u003e, 170-192.\u003c/li\u003e\n \u003cli\u003eHuxley, H., and Hanson, J. (1954). Changes in the Cross-Striations of Muscle during Contraction and Stretch and their Structural Interpretation. Nature \u003cem\u003e173\u003c/em\u003e, 973-976.\u003c/li\u003e\n \u003cli\u003eHuxley, A.F., and Niedergerke, R. (1954). Structural changes in muscle during contraction; interference microscopy of living muscle fibres. Nature \u003cem\u003e173\u003c/em\u003e, 971-973.\u003c/li\u003e\n \u003cli\u003eRome, L.C., and Sosnicki, A.A. (1991). Myofilament overlap in swimming carp. II. Sarcomere length changes during swimming. Am J Physiol \u003cem\u003e260\u003c/em\u003e, C289-296.\u003c/li\u003e\n \u003cli\u003eTu, M.S., and Daniel, T.L. (2004). Cardiac-like behavior of an insect flight muscle. J Exp Biol \u003cem\u003e207\u003c/em\u003e, 2455-2464.\u003c/li\u003e\n \u003cli\u003eLayland, J., Young, I.S., and Altringham, J.D. (1995). The length dependence of work production in rat papillary muscles in vitro. J Exp Biol \u003cem\u003e198\u003c/em\u003e, 2491-2499.\u003c/li\u003e\n \u003cli\u003eAzizi, E. (2014). Locomotor function shapes the passive mechanical properties and operating lengths of muscle. Proceedings of the Royal Society B: Biological Sciences \u003cem\u003e281\u003c/em\u003e, 20132914.\u003c/li\u003e\n \u003cli\u003eGidmark, N.J., Konow, N., LoPresti, E., and Brainerd, E.L. (2013). Bite force is limited by the force\u0026ndash;length relationship of skeletal muscle in black carp, \u003cem\u003eMylopharyngodon piceus\u003c/em\u003e. Biol Lett \u003cem\u003e9\u003c/em\u003e, 20121181.\u003c/li\u003e\n\u003c/ol\u003e\n\u003col start=\"10\"\u003e\n \u003cli\u003eMorgan, D.L., and Proske, U. (2006). Sarcomere popping requires stretch over a range where total tension decreases with length. J Physiol \u003cem\u003e574\u003c/em\u003e, 627-628.\u003c/li\u003e\n \u003cli\u003eAnapol, F., and Herring, S. (1989). Length-tension relationships of masseter and digastric muscles of miniature swine during ontogeny. J Exp Biol \u003cem\u003e143\u003c/em\u003e, 1-16.\u003c/li\u003e\n \u003cli\u003eWeijs, W.A., and van der Wielen-Drent, T.K. (1983). The relationship between sarcomere length and activation pattern in the rabbit masseter muscle. Arch Oral Biol \u003cem\u003e28\u003c/em\u003e, 307-315.\u003c/li\u003e\n \u003cli\u003eNordstrom, S.H., and Yemm, R. (1974). The relationship between jaw position and isometric active tension produced by direct stimulation of the rat masseter muscle. Arch Oral Biol \u003cem\u003e19\u003c/em\u003e, 353-359.\u003c/li\u003e\n \u003cli\u003eThexton, A.J., and Hiiemae, K.M. (1975). The twitch-contraction characteristics of opossum jaw musculature. Arch Oral Biol \u003cem\u003e20\u003c/em\u003e, 743-748.\u003c/li\u003e\n \u003cli\u003eManns, A., Miralles, R., and Palazzi, C. (1979). EMG, bite force, and elongation of the masseter muscle under isometric voluntary contractions and variations of vertical dimension. The Journal of prosthetic dentistry \u003cem\u003e42\u003c/em\u003e, 674-682.\u003c/li\u003e\n \u003cli\u003eMackenna, B.R., and T\u0026uuml;rker, K.S. (1978). Twitch tension in the jaw muscles of the cat at various degrees of mouth opening. Arch Oral Biol \u003cem\u003e23\u003c/em\u003e, 917-920.\u003c/li\u003e\n \u003cli\u003eGranatosky, M.C., and Ross, C.F. (2020). Differences in muscle mechanics underlie divergent optimality criteria between feeding and locomotor systems. Journal of Anatomy \u003cem\u003e237\u003c/em\u003e, 1072-1086.\u003c/li\u003e\n \u003cli\u003eEng, C.M., Ward, S.R., Vinyard, C.J., and Taylor, A.B. (2009). The morphology of the masticatory apparatus facilitates muscle force production at wide jaw gapes in tree-gouging common marmosets (Callithrix jacchus). J Exp Biol \u003cem\u003e212\u003c/em\u003e, 4040-4055.\u003c/li\u003e\n \u003cli\u003eRull, M., Solomon, J., and Konow, N. (2020). Elastic recoil action amplifies jaw closing speed in an aquatic feeding salamander. Proceeding of the Royal Society, series B Biology \u003cem\u003e287\u003c/em\u003e, 32429804.\u003c/li\u003e\n \u003cli\u003eKonow, N., and Roberts, T.J. (2015). The series elastic shock absorber: tendon elasticity modulates energy dissipation by muscle during burst deceleration. Proceedings of the Royal Society B: Biological Sciences \u003cem\u003e282\u003c/em\u003e, 20142800.\u003c/li\u003e\n \u003cli\u003eLaird, M.F., Granatosky, M.C., Taylor, A.B., and Ross, C.F. (2020). Muscle architecture dynamics modulate performance of the superficial anterior temporalis muscle during chewing in capuchins. Sci Rep-Uk \u003cem\u003e10\u003c/em\u003e, 6410.\u003c/li\u003e\n \u003cli\u003eAzizi, E., Brainerd, E.L., and Roberts, T.J. (2008). Variable gearing in pennate muscles. Proc Natl Acad Sci U S A \u003cem\u003e105\u003c/em\u003e, 1745-1750.\u003c/li\u003e\n \u003cli\u003eLaird, M.F., Iriarte-Diaz, J., Byron, C.D., Granatosky, M.C., Taylor, A.B., and Ross, C.F. (2023). Gape drives regional variation in temporalis architectural dynamics in tufted capuchins. \u003cem\u003e378\u003c/em\u003e, 20220550.\u003c/li\u003e\n \u003cli\u003eHiiemae, K.M. (1978). Mammalian mastication : a review of the activity of the jaw muscles and the movements they produce in chewing. In Development, Function and Evolution of Teeth, P.M.B.a.K.A. Joysey, ed., pp. 359\u0026ndash;398.\u003c/li\u003e\n \u003cli\u003eHolt, N.C., and Azizi, E. (2014). What drives activation-dependent shifts in the force-length curve? Biology letters \u003cem\u003e10\u003c/em\u003e.\u003c/li\u003e\n \u003cli\u003eSundar, S., Rimkus, B., Meemaduma, P.S., deLap, S., LaFave, N., Racca, A.W., Hettige, P., Moore, J., Gage, M., Shehaj, A., et al. (2022). Bridging the muscle genome to phenome across multiple biological scales. Journal of Experimental Biology.\u003c/li\u003e\n \u003cli\u003eTijs, C., Konow, N., and Biewener, A.A. (2021). Effect of muscle stimulation intensity on the heterogeneous function of regions within an architecturally complex muscle. Journal of Applied Physiology \u003cem\u003e130\u003c/em\u003e, 941-951.\u003c/li\u003e\n \u003cli\u003eReed, D.A., and Ross, C.F. (2010). The influence of food material properties on jaw kinematics in the primate, \u003cem\u003eCebus\u003c/em\u003e. Archives of Oral Biology \u003cem\u003e55\u003c/em\u003e, 946-962.\u003c/li\u003e\n \u003cli\u003eWilliams, S.H., Wright, B.W., Truong, V.d., Daubert, C.R., and Vinyard, C.J. (2005). Mechanical properties of foods used in experimental studies of primate masticatory function. Am J Primatol \u003cem\u003e67\u003c/em\u003e, 329-346.\u003c/li\u003e\n \u003cli\u003eAhn, A.N., Konow, N., Tijs, C., and Biewener, A.A. (2018). Different segments within vertebrate muscles can operate on different regions of their force\u0026ndash;length relationships. Integrative and Comparative Biology \u003cem\u003e58\u003c/em\u003e, 219-231.\u003c/li\u003e\n \u003cli\u003eClose, R.I. (1972). Dynamic properties of mammalian skeletal muscles. Physiol Rev \u003cem\u003e52\u003c/em\u003e, 129-197.\u003c/li\u003e\n \u003cli\u003eBuchanan, T.S. (1995). Evidence that maximum muscle stress is not a constant: differences in specific tension in elbow flexors and extensors. Med Eng Phys \u003cem\u003e17\u003c/em\u003e, 529-536.\u003c/li\u003e\n \u003cli\u003eWeijs, W.A., and Dantuma, R. (1975). Electromyography and mechanics of mastication in the albino rat. Journal of Morphology \u003cem\u003e146\u003c/em\u003e, 1-33.\u003c/li\u003e\n \u003cli\u003eHolt, N.C., and Azizi, E. (2016). The effect of activation level on muscle function during locomotion: are optimal lengths and velocities always used? Proceedings of the Royal Society of London B: Biological Sciences \u003cem\u003e283\u003c/em\u003e, 20152832.\u003c/li\u003e\n \u003cli\u003eSchwaner, M.J., Mayfield, D.L., Azizi, E., and Daley, M.A. (2024). Linking in vivo muscle dynamics to in situ force-length and force-velocity reveals that guinea fowl lateral gastrocnemius operates at shorter than optimal lengths. Journal of Experimental Biology.\u003c/li\u003e\n \u003cli\u003eRack, P.M., and Westbury, D.R. (1969). The effects of length and stimulus rate on tension in the isometric cat soleus muscle. J Physiol \u003cem\u003e204\u003c/em\u003e, 443-460.\u003c/li\u003e\n \u003cli\u003eRassier, D.E., MacIntosh, B.R., and Herzog, W. (1999). Length dependence of active force production in skeletal muscle. Journal of Applied Physiology \u003cem\u003e86\u003c/em\u003e, 1445-1457.\u003c/li\u003e\n \u003cli\u003eMayfield, D.L., and Holt, N.C. (2025). Does force depression resulting from shortening against series elasticity contribute to the activation dependence of optimum length? \u003cem\u003e328\u003c/em\u003e, C528-C540.\u003c/li\u003e\n \u003cli\u003eHolt, N.C., and Williams, C.D. (2018). Can Strain Dependent Inhibition of Cross-Bridge Binding Explain Shifts in Optimum Muscle Length? Integrative and Comparative Biology \u003cem\u003e58\u003c/em\u003e, 174-185.\u003c/li\u003e\n \u003cli\u003eHolt, N.C., and Mayfield, D.L. (2023). Muscle-tendon unit design and tuning for power enhancement, power attenuation, and reduction of metabolic cost. Journal of Biomechanics, 111585.\u003c/li\u003e\n \u003cli\u003eHENNEMAN, E. (1957). Relation between Size of Neurons and Their Susceptibility to Discharge. Science \u003cem\u003e126\u003c/em\u003e, 1345-1347.\u003c/li\u003e\n \u003cli\u003eHolt, N.C., Wakeling, J.M., and Biewener, A.A. (2014). The effect of fast and slow motor unit activation on whole-muscle mechanical performance: the size principle may not pose a mechanical paradox. Proceedings of the Royal Society B: Biological Sciences \u003cem\u003e281\u003c/em\u003e, 24695429.\u003c/li\u003e\n \u003cli\u003eLlewellyn, M.E., Thompson, K.R., Deisseroth, K., and Delp, S.L. (2010). Orderly recruitment of motor units under optical control in vivo. Nat Med \u003cem\u003e16\u003c/em\u003e, 1161-U1144.\u003c/li\u003e\n \u003cli\u003eHedrick, T.L., Tobalske, B.W., and Biewener, A.A. (2003). How cockatiels (\u003cem\u003eNymphicus hollandicus\u003c/em\u003e) modulate pectoralis power output across flight speeds. J Exp Biol \u003cem\u003e206\u003c/em\u003e, 1363-1378.\u003c/li\u003e\n \u003cli\u003eKonow, N., Collias, A., and Biewener, A.A. (2020). Skeletal muscle shape change in relation to varying foce requirements across locomotor conditions. Frontiers in Physiology \u003cem\u003e11\u003c/em\u003e, 143.\u003c/li\u003e\n \u003cli\u003eGabaldon, A.M., Nelson, F.E., and Roberts, T.J. (2008). Relative shortening velocity in locomotor muscles: turkey ankle extensors operate at low V/Vmax. Am J Physiol Regul Integr Comp Physiol \u003cem\u003e294\u003c/em\u003e, 200-210.\u003c/li\u003e\n \u003cli\u003eRome, L.C. (1998). Some advances in integrative muscle physiology. Comp Biochem Physiol B Biochem Mol Biol \u003cem\u003e120\u003c/em\u003e, 51-72.\u003c/li\u003e\n \u003cli\u003eSinha, A., and Bhargav, A. (2020). Effect of state transition, drying kinetics and moisture content on Young\u0026apos;s modulus variation during thermal drying of hygroscopic food materials. Journal of Food Engineering \u003cem\u003e279\u003c/em\u003e, 109957.\u003c/li\u003e\n \u003cli\u003eP\u0026uuml;ffel, F., Walthaus, O.K., Kang, V., and Labonte, D. (2023). Biomechanics of cutting: sharpness, wear sensitivity and the scaling of cutting forces in leaf-cutter ant mandibles. \u003cem\u003e378\u003c/em\u003e, 20220547.\u003c/li\u003e\n \u003cli\u003eAzizi, E., and Roberts, T.J. (2010). Muscle performance during frog jumping: influence of elasticity on muscle operating lengths. Proc Biol Sci \u003cem\u003e277\u003c/em\u003e, 1523-1530.\u003c/li\u003e\n \u003cli\u003eLangenbach, G.E.J., and van Eijden, T.M.G.J. (2001). Mammalian Feeding Motor Patterns. Amer. Zool. \u003cem\u003e41\u003c/em\u003e, 1338-1351.\u003c/li\u003e\n \u003cli\u003eRam, Y., and Ross, C.F. (2020). Jaw Elevator Muscle Coordination during Rhythmic Mastication in Primates: Are Triplets Units of Motor Control? Brain, Behavior and Evolution \u003cem\u003e95\u003c/em\u003e, 1-14.\u003c/li\u003e\n \u003cli\u003ePeck, C.C., and Hannam, A.G. (2007). Human jaw and muscle modelling. Archives of Oral Biology \u003cem\u003e52\u003c/em\u003e, 300-304.\u003c/li\u003e\n \u003cli\u003eShi, J., Curtis, N., Fitton, L.C., O\u0026apos;Higgins, P., and Fagan, M.J. (2012). Developing a musculoskeletal model of the primate skull: Predicting muscle activations, bite force, and joint reaction forces using multibody dynamics analysis and advanced optimisation methods. Journal of Theoretical Biology \u003cem\u003e310\u003c/em\u003e, 21-30.\u003c/li\u003e\n \u003cli\u003eCox, P.G., and Watson, P.J. (2023). Masticatory biomechanics of red and grey squirrels (\u0026lt;i\u0026gt;Sciurus vulgaris\u0026lt;/i\u0026gt; and \u0026lt;i\u0026gt;Sciurus carolinensis\u0026lt;/i\u0026gt;) modelled with multibody dynamics analysis. \u003cem\u003e10\u003c/em\u003e, 220587.\u003c/li\u003e\n \u003cli\u003eLautenschlager, S., Gill, P.G., Luo, Z.-X., Fagan, M.J., and Rayfield, E.J. (2018). The role of miniaturization in the evolution of the mammalian jaw and middle ear. Nature \u003cem\u003e561\u003c/em\u003e, 533-537.\u003c/li\u003e\n \u003cli\u003eLedogar, J.A., Senck, S., Villmoare, B.A., Smith, A.L., Weber, G.W., Richmond, B.G., Dechow, P.C., Ross, C.F., Grosse, I.R., Wright, B.W., et al. (2022). Mechanical compensation in the evolution of the early hominin feeding apparatus. \u003cem\u003e289\u003c/em\u003e, 20220711.\u003c/li\u003e\n \u003cli\u003eBasmajian, J., Stecko, G. (1962). A new bipolar electrode for electromyography. Journal of Appplied Physiology \u003cem\u003e17\u003c/em\u003e, 849.\u003c/li\u003e\n \u003cli\u003eRichards, C.T., and Biewener, A.A. (2007). Modulation of in vivo muscle power output during swimming in the African clawed frog (Xenopus laevis). J Exp Biol \u003cem\u003e210\u003c/em\u003e, 3147-3159.\u003c/li\u003e\n \u003cli\u003eAfoakwa, E.O., Paterson, A., Fowler, M., and Vieira, J. (2009). Microstructure and mechanical properties related to particle size distribution and composition in dark chocolate. International Journal of Food Science \u0026amp; Technology \u003cem\u003e44\u003c/em\u003e, 111-119.\u003c/li\u003e\n \u003cli\u003eVinyard, C.J., and Taylor, A.B. (2010). A Preliminary Analysis of the Relationship Between Jaw-Muscle Architecture and Jaw-Muscle Electromyography During Chewing Across Primates. The Anatomical Record: Advances in Integrative Anatomy and Evolutionary Biology \u003cem\u003e293\u003c/em\u003e, 572-582.\u003c/li\u003e\n \u003cli\u003eOtten, E. (1987). A Myocybernetic Model of the Jaw System of the Rat. Journal of Neuroscience Methods \u003cem\u003e21\u003c/em\u003e, 287-302.\u003c/li\u003e\n \u003cli\u003eDeLap, S., Tran, V., and Konow, N. (2024). DeepXROMM. A set of user-friendly tools for DeepLabCut-XMALab marker tracking. (https://pypi.org/project/deepxromm/).\u003c/li\u003e\n \u003cli\u003eMarsh, R.L. (2016). Speed of sound in muscle for use in sonomicrometry. Journal of Biomechanics.\u003c/li\u003e\n \u003cli\u003eEng, C.M., Konow, N., Tijs, C., Holt, N.C., and Biewener, A.A. (2019). In vivo force\u0026ndash;length and activation dynamics of two distal rat hindlimb muscles in relation to gait and grade. The Journal of Experimental Biology \u003cem\u003e222\u003c/em\u003e, jeb205559.\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Tables","content":"\u003cp\u003e\u003cstrong\u003eTable 1. Summary data from \u003cem\u003ein situ\u003c/em\u003e experiment on rat superficial masseter\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"0\" cellspacing=\"0\" cellpadding=\"0\" width=\"379\"\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003eRat\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003eStim (V)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e\u003cem\u003eF\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e (N) (submax)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e\u003cem\u003eF\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e (N) (twitch)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e (mm) (submax)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e\u003cem\u003eL\u003c/em\u003e\u003csub\u003eO\u003c/sub\u003e (mm) \u0026nbsp; \u0026nbsp; (twitch)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003eTension (submax)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.5*\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e12.8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e10.5\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.7*\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e8.7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e14.1\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e5.7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e0.4\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e6.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e12.8\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e4\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e11\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.6\u003csup\u003e#\u003c/sup\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e12.3\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e3.5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e0.7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e8.5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e9.6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e6.9\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e3.7\u003csup\u003e#\u003c/sup\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e-\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e7.3\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e1.9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e5.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e8.1\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e12\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e1.5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e3.1\u0026nbsp;\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e3.4\u0026nbsp;\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e\u0026nbsp;11.1\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.76253%;\"\u003e\n \u003cp\u003e9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 9.49868%;\"\u003e\n \u003cp\u003e8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e7.3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 14.248%;\"\u003e\n \u003cp\u003e2.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e5.8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 17.4142%;\"\u003e\n \u003cp\u003e6.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" style=\"width: 15.8311%;\"\u003e\n \u003cp\u003e9.0\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e*Twitch FL not measured. \u003csup\u003e#\u003c/sup\u003eFL curves from motor position, not fascicle length.\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":true,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Force-length operation, length-tension relationship, musculoskeletal function","lastPublishedDoi":"10.21203/rs.3.rs-6739430/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-6739430/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"The operating length ranges of mammalian jaw muscles have been estimated using twitch contractions or force measurements at the bite point, prompting a consensus that jaw muscles operate at short lengths on their force-length (FL) curve. However, since activation intensity truncates muscle optimal length (LO), we hypothesized that LO of rat superficial masseter (SM) would decrease with activation intensity, with high-force biting involving muscle shortening from long lengths on the FL curve. We measured muscle activation, strain, and force in vivo during biting on food with varying hardness and mapped the in vivo data from each muscle (N = 6) onto its FL relationship, measured in situ. Submaximal LO was approx. 12% shorter than twitch LO, and SM bite forces averaged 4.1 ± 3.9 N/cm2 (mean ± S.D.) and reached 10.6 N/cm2, corresponding to muscle activation and food hardness. Length operation ranged from 7% below LO (ascending FL plateau), to 27% beyond (descending limb). The finding that jaw muscles operate at long, potentially unstable lengths, particularly during hard food biting significantly expands our understanding of skeletal muscle function, with broad implications for craniofacial evolution, muscle mechanics and control, and healthy as well as pathological function of the jaw musculoskeletal system.","manuscriptTitle":"Rat superficial masseter operates at long, unstable lengths during biting","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-09-12 07:29:27","doi":"10.21203/rs.3.rs-6739430/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-09-22T12:14:27+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-09-22T11:25:08+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-09-18T10:54:52+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"296183554725504582459870943947317163259","date":"2025-09-18T10:36:54+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-09-15T13:35:42+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"251062299696226419620224489633500821546","date":"2025-09-15T13:29:26+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-09-12T12:58:47+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-09-11T10:17:27+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2025-09-10T17:08:26+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"22985df6-6d07-4708-8a1b-019a50fd5f77","owner":[],"postedDate":"September 12th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":54561157,"name":"Biological sciences/Evolution"},{"id":54561158,"name":"Biological sciences/Physiology"}],"tags":[],"updatedAt":"2025-11-03T16:09:17+00:00","versionOfRecord":{"articleIdentity":"rs-6739430","link":"https://doi.org/10.1038/s41598-025-21953-z","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2025-10-30 15:58:13","publishedOnDateReadable":"October 30th, 2025"},"versionCreatedAt":"2025-09-12 07:29:27","video":"","vorDoi":"10.1038/s41598-025-21953-z","vorDoiUrl":"https://doi.org/10.1038/s41598-025-21953-z","workflowStages":[]},"version":"v1","identity":"rs-6739430","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-6739430","identity":"rs-6739430","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}