The present protocol describes concise experimental details on the evaluation and interpretation of in vivo torque data obtained via electrical stimulation of the common peroneal nerve in anesthetized pigs.
Reliable assessment of skeletal muscle strength is arguably the most important outcome measure in neuromuscular and musculoskeletal disease and injury studies, particularly when evaluating regenerative therapies' efficacy. Additionally, a critical aspect of translating many regenerative therapies is the demonstration of scalability and effectiveness in a large animal model. Various physiological preparations have been established to evaluate intrinsic muscle function properties in basic science studies, primarily in small animal models. The practices may be categorized as: in vitro (isolated fibers, fiber bundles, or whole muscle), in situ (muscle with intact vascularization and innervation but distal tendon attached to a force transducer), and in vivo (structures of the muscle or muscle unit remain intact). There are strengths and weaknesses to each of these preparations; however, a clear advantage of in vivo strength testing is the ability to perform repeated measurements in the same animal. Herein, the materials and methods to reliably assess isometric torque produced by the hindlimb dorsiflexor muscles in vivo in response to standard peroneal electrical stimulation in anesthetized pigs are presented.
The primary function of skeletal muscle is to produce force, which ultimately makes activities such as breathing, eating, and ambulating possible. Conditions that reduce skeletal muscle functional capacity can lead to diminished performance (occupational or sport), disability, or death. For example, the maintenance of muscle mass and function in aging populations is positively associated with quality of life and the capacity to perform basic and instrumental activities of daily living1,2. And, declining muscle strength in Duchenne muscular dystrophy patients results in the inability to ambulate and respiratory failure, ultimately contributing to premature mortality3,4,5. Thus, muscle strength measurement is a critical outcome measure in studies involving neuromuscular disease or injury.
Maximal voluntary isometric or isokinetic torque (and/or fatigue index) is often used as an index of functional capacity in clinical studies6. In animal studies, analogous measurements can be made in vivo using electrical nerve stimulation while under anesthesia. Notably, in vivo preparations are minimally invasive with musculature, tendons, vasculature, and innervation remaining intact and therefore permit repeated functional assessments7,8,9,10,11. This preparation is commonly used in small rodent models and to a lesser extent in larger animal models such as rabbits12, dogs13,14, sheep15, and pigs16,17. The general use of such methodology could be impactful to many translational research studies, such as in genetically engineered porcine (pig) models of spinal muscular atrophy (SMA)18. Herein, methods to assess nerve stimulation-induced maximal isometric torque of the porcine dorsiflexor muscle group in vivo are presented. The techniques presented were initially adapted from those developed originally to assess mouse anterior crural muscle torque19,20 and subsequently refined through experience investigating torque producing capacity following injury17,21,22,23,24,25,26,27,28 and during development16 in various porcine models.
This protocol highlights in vivo isometric torque measurement using methodology that requires a computer integrated with a load cell and electrical stimulator. The methods presented here use a commercially available integrated Swine Isometric Footplate Test Apparatus, platform apparatus, and corresponding software (see Table of Materials). However, the methodology can be adapted to use other commercially available or custom-made software, data-acquisition devices, and stimulators. These methods are intended for use in a dedicated large animal surgical suite replete with standard equipment such as: locking surgical table, second locking table of equal height for the testing platform, ventilator and monitoring devices, and heating mat or other devices to maintain body temperature.
The following team members are needed to conduct these methods: one skilled anesthesia technician and two study personnel to perform the functional testing. These people will work together for the initial stabilization of the limb on the platform apparatus. Then, one of the two personnel will be responsible for the electrode placement/positioning and the other for the computer applications during the testing.
All animal experiments were conducted in compliance with the Animal Welfare Act, the implementing Animal Welfare Regulations, and the principles of the Guide for the Care and Use of Laboratory Animals. Prior testing has demonstrated that these methods are reliable26 and have no adverse effects on the health or limb function of the pig. Testing has been conducted as often as weekly without any adverse events23. Additionally, testing pre- and post-surgical interventions during the same day can be performed without placing untoward stress on the animal or inducing neuromuscular dysfunction.
1. Computer set up
2. Anesthesia preparation and maintenance
3. Evaluation of in vivo isometric torque
Figure 1: Pictures from various vantage points show the pig attachment to the foot-plate and anatomical alignment onto the frame. Anatomical landmarks for the anterior compartment muscles, fibular head, knee, tibial plateau, and femur are noted. Note the placement of subdermal electrode pairs on the lateral side of the leg. Please click here to view a larger version of this figure.
4. Protocol for torque-joint angle analysis
5. Protocol for torque-frequency analysis
6. Data analysis
Reliability and optimization of the in vivo testing parameters of the anterior compartment of the pig have previously been reported26. Comparative data across rodents and pigs for torque-frequency has also been reported27.
During the in vivo assessment, visualization of the torque waveform is needed in real-time to ensure appropriate anterior compartment activation. The waveforms should only reflect dorsiflexion. The waveforms should have a smooth, rounded appearance and an apparent tetanic plateau (Figure 2A). Inconsistencies or perturbations of the waveform indicate various experimental limitations, such as inadequate stimulation, improper electrode placement, or inadequate depth of anesthesia (Figure 2B).
Figure 3A is a twitch torque-time tracing with an arrow indicating 50% max torque. Time-to-peak contraction should start at the initiation of the stimulator and end when maximal twitch torque is achieved (representative time bars are shown below the tracing). Half-relaxation for a twitch should start at the maximal twitch torque and end at 50% maximal twitch torque (representative time bars are shown below the tracing). Figure 3B is a tetanic torque-time tracing with an arrow indicating 50% max torque. Unlike twitches that are ideal in terms of a definitive and timely maximal torque, tetanic contractions have greater variability in the timing of maximal torque concerning when the stimulator starts and ends, requiring a more nuanced approach to contractile property analysis. Time-to-peak contraction should begin with the initiation of the stimulator and stop somewhere between 90%-100% of maximal torque. The time bars in Figure 3B show a cutoff of 95% maximal torque. This is helpful in cases such as the selected representative data because the maximal torque is not reached until late in the plateau phase. A complementary analysis to time-to-peak is the average rate of contraction. The dashed bars on the ascending limb of the torque-time tracing represent a range of 30%-70% maximal torque. The average rate of contraction should be started at the start of stimulation and capture the average rate change between 30%-70% maximal torque. These are recommended ranges, and individual research groups can determine the ideal range around 50% (e.g., ±10%). The important part is to be consistent within and across studies. In contrast to the twitch, tetanic contraction half-relaxation should start at the end of stimulation instead of maximal torque for the same reason mentioned above with time-to-peak. The time bars in Figure 3B represent the time between the end of the stimulation and reaching 50% relaxation. A complementary analysis to half-relaxation is the average rate of relaxation. The dashed bars on the descending limb of the torque-tracing represent the same 30%-70% maximal torque range as the ascending limb. The average relaxation rate should start at the end of stimulation and capture the average rate of change between 30%-70% maximal torque. Again, these are recommended ranges. One critical note: do not confuse the average rate of contraction/relaxation with the maximal rate of contraction/relaxation. The maximal rate represents the single most remarkable rate change between two adjacent data points and can be widely variable.
Several twitch and contractile properties can be analyzed to gain insight into fiber type and excitation-contraction coupling attributes of the skeletal muscles10,32. Over-interpretation of the twitch and contractile properties are cautioned; they represent suggestions and rationale for further cellular-level interrogation and are not necessarily indicative. In general, rates of contractility can reflect sarcoplasmic reticulum calcium release and myosin heavy chain isoform enzymatic rate. In contrast, relaxation rates can reflect sarco(endo)plasmic reticulum calcium ATPase enzyme rate and isoform. These properties can be influenced by fatigue, muscle damage, exercise training, and numerous pathologies (e.g., disuse atrophy).
Figure 4 depicts representative values for torque-frequency and torque-joint angle relationships for uninjured limbs. These data are representative of a wide range of pig sizes.
A representative, experimental analysis of surface EMG was conducted during in vivo muscle analysis (Figure 5) to demonstrate experimental control of rate coding and total muscle activity. Adhesive EMG electrodes were placed at the mid-belly of the peroneus tertius. A ground electrode was placed on the knee to minimize stimulation artifact, and stimulation electrode needles were placed around the peroneal nerve proximal to the muscle location. Simultaneous torque and EMG recordings were made at stimulation frequencies of 20, 60, and 100 Hz. The number of stimulator pulses (red bars in Figure 5) reflects the quotient of stimulation duration and time between pulses. For example, a 20 Hz stimulation frequency means a pulse every 50 ms; therefore, 400 ms stimulation duration divided by 50 ms between pulses equals eight pulses delivered (Figure 5A). The stimulator pulses are delivered to the nerve axon via percutaneous needle electrode placement and produce a similar number of electrical muscle pulses (i.e., 20 Hz equals 8 EMG recordings), demonstrating the experimental control of action potential frequency of the muscle group of interest. The raw EMG recordings can be converted via root-mean-square analysis (EMG RMS) to visualize the total muscle activity with increasing stimulation frequency. The area under the curve (AUC) analysis is one way to quantify the EMG RMS to determine changes in whole muscle activity. Representative AUCs for each EMG RMS stimulation frequency is provided (Figure 5A–C).
Figure 2: Representative high and low-quality waveforms. (A) Isometric waveforms present in a square-wave appearance, with a notable fluid plateau. (B) Low-quality waveforms can be due to inadequate stimulation or improper electrode placement. In these instances, repositioning of the electrodes is needed. For both A and B, the stimulator pulses (red bars) are indicated. Please click here to view a larger version of this figure.
Figure 3: Twitch and tetanic contractile property analysis. (A) Representative twitch (1 Hz) and (B) tetanic (100 Hz) torque-time tracings are modified to detail contractile properties. The red arrow on each graph shows 50% maximal torque. Blue and black bars beneath the tracings show time-to-peak and half-relaxation time durations, respectively. Dashed bars on the ascending and descending limbs of the tetanic torque-time tracing represent a range of 30%-70% maximal torque that can be used to determine the average rate of contraction or relaxation. Please click here to view a larger version of this figure.
Figure 4: Example data of torque-joint angle and torque-frequency. Data provided is from a range of female Yorkshire Cross pigs at 2.9-6.3 months; 39.4-75.4 kg body mass; all considered healthy control at the evaluation time. During all testing, core body temperature was maintained at 37 °C. (A) Torque normalized to body mass is evaluated at ankle joints of 0-50° of plantarflexion; note the peak torque is at determined at 30°. (B) Torque normalized to body mass is evaluated at various stimulation frequencies from 10-100 Hz; note these evaluations were conducted with the ankle joint at 30° of plantarflexion. (C) Individual torque tracings for each of the stimulation frequencies were evaluated. Please click here to view a larger version of this figure.
Figure 5: Concurrent in vivo isometric torque and EMG measurements. Simultaneous EMG and torque recordings at representative stimulation frequencies of (A) 20, (B) 60, and (C) 100 Hz collected from a female Yorkshire pig (~90 kg body mass). Stimulator pulses (red bars) were delivered according to the set stimulation frequency. Raw EMG recordings were converted to root-mean-square (EMG RMS) to visualize total muscle activity with increasing stimulation frequency. Representative EMG RMS curves were analyzed for the area under the curve (AUC), and AUCs are provided for each stimulation frequency. Please click here to view a larger version of this figure.
Critical steps, modifications, and troubleshooting
For minimizing data variability and maximizing the approach's success, the following critical steps are highlighted.
Optimal nerve stimulation
This experimental approach starts with nerve axon depolarization and relies on correct electrode placement and optimized electrical stimulation. A post-mortem analysis of nerve anatomy related to boney landmarks can help visualize proper electrode placement during testing. Acquiring maximal twitch torque helps determine appropriate current (in milliamperes; mA) delivered to the nerve axon. There are two values to keep in mind when optimizing nerve stimulation at the onset of testing: (1) the twitch-to-tetanic ratio is ~1:5, e.g., ~2 N·m twitch torque corresponds to a 10 N·m tetanic torque (Figure 3); and (2) the typical torque to body mass is ~0.3 N·m per kg body mass (Figure 4). If the peak twitch torques appear low, remove the electrodes and attempt another placement. Be sure to check stimulator settings, BNC connections, and electrode connections. Electrode re-placement may be needed in between contractions if there is too much movement during positioning of the limb between joint angles, as noted above (Figure 2). Please note experimental and interventional approaches could impact these values.
Proper biomechanical alignment
Starting muscle length influences muscle contractile force (the length-tension relationship), and muscle length can change based on hip, knee, and ankle joint alignment. The joint angles must be standardized between limbs and among pigs. A 90° ankle joint angle is strongly recommended for the hip and knee. A slightly plantarflexed ankle position (~30° from the neutral 0° ankle joint angle) is optimal for peak strength. It reflects the natural anatomical position of the ankle joint in both pigs and dogs while standing. All joints should also be parallel with the foot-pedal and torque transducers to avoid loss of measurable torque due to the contribution of a perpendicular torque vector. Inspecting the hip-knee-ankle joint angles and foot-pedal-joint alignment is strongly recommended after securing the foot to the foot-pedal and securing the knee joint with the limb clamping bars (Figure 1). If there is misalignment, unlock and remove the bars and reposition the pig on the surgical table. While standardizing joint angles across studies is critical for minimizing data variance, there are limitations to biomechanical alignment that are notable, discussed below.
Significance with respect to existing or alternative methods
Alterative examples of clinically relevant and non-invasive assessments of muscle function that could be used for porcine models include treadmill walking distance, EMG, and active muscle shear wave electrography. As the 6 min walk test in humans, a treadmill walking test can evaluate disease progression and intervention success in large animals33,34,35. Typically, after an acclimation period, animals are walked until the end of compliance at different treadmill velocities and/or incline levels. Food rewards are often necessary to achieve maximal motivation. However, treadmill walking outcomes offer only indirect interpretations of muscle contractile function due to limitations such as subject motivation, non-maximal motor unit recruitment, and inherent co-dependence on other body systems such as the cardiovascular, skeletal, and respiratory systems.
On the other hand, EMG offers a slightly better direct assessment of the skeletal muscle system, as EMG electrodes are placed directly on the muscle group of interest36,37,38. EMG electrodes then measure the collective muscle activity (depolarized muscle fibers). This muscle activity is based on motor unit recruitment and rate coding (the frequency of action potentials sent to recruited motor units). However, separating the relative contributions of motor unit recruitment versus rate coding is impossible with surface EMG. Further, EMG relies on subject willingness to generate maximal contractions, and this level of cooperation is unlikely in large-animal models. While it can be informative to assess changes in EMG during the gait cycle, these data do not represent a maximal functional ability of the skeletal muscle group of interest. Ultrasound-based imaging utilizing B-mode and shear-wave elastography is another non-invasive modality used to evaluate muscle function. There is a good correlation between Young's modulus measured by elastography and increasing muscle loads39,40. Shear-wave elastography has been validated and used as a quantitative measure of passive tissue stiffness41,42,43,44,45, including in a porcine volumetric muscle loss injury model23. It may also be used as an indirect measurement of active muscle force production39. However, limitations akin to EMG for subject willingness and cooperation to perform contractions are still present.
The in vivo protocol described here, in contrast to treadmill walking distance and EMG, provides a reliable, reproducible, and maximal assessment of muscle function. This protocol evokes muscle contractions in a controlled, quantifiable manner that is independent of motivation. Specifically, percutaneous electrodes are used to stimulate nerves axons bypassing the central nervous system. Depolarization of the nerve axons engages all motor units eliminating variability associated with motor unit recruitment. Additionally, the investigator controls rate coding (stimulation frequency). The resulting neuromuscular physiology that applies to this approach starts with voltage-gated sodium channel activation at the nodes of Ranvier. All subsequent (or downstream) physiology is engaged, including excitation-contraction coupling and cross-bridge cycling. A significant advantage of the in vivo non-invasive muscle analysis is that contractile muscle function can be measured repeatedly, for example, weekly, to monitor muscle strength after injury, intervention, or over a disease progression.
Limitations of the method
The in vivo equipment described in this protocol permits passive and active isometric torque as a function of joint angle and stimulation frequency. The testing apparatus used does not support the measurement of dynamic contractions (e.g., isokinetic eccentric or concentric contractions). The apparatus allows a 105° range of motion to characterize the torque-joint angle relationship and uses a load cell with a maximum torque range of ~50 N·m. Specific experimental questions may require performance characteristics outside of these specifications. Notably, the load cell on this described apparatus may be exchanged for greater torque ranges if needed.
The protocol described herein to measure maximal neuromuscular strength in vivo has notable limitations. First, this method requires anesthesia, which may be conducted differently per animal facility protocols and resources. Anesthetics are known to have varying effects on neuromuscular function and have been shown to alter mouse in vivo dorsiflexor torque production in an anesthetic-type and -dose-dependent manner29. The differential effects of anesthetics on the large-animal in vivo torque are unclear; therefore, control and experimental groups must have the same anesthesia agents (e.g., all groups administered ketamine) to control this variability. Second, reliance on in vivo diffusion patterns limits exploration of cellular mechanisms of contractile dysfunction and acute drug toxicities. For example, caffeine can be used during in vitro organ bath testing of an isolated muscle to stimulate sarcoplasmic reticulum calcium release, bypassing excitation-contraction coupling46 directly. The amount of caffeine to induce this effect (mM) is lethal in an in vivo setting. Drug influences on the whole body (e.g., kidney/liver stress) and subsequent factors secreted into circulation will need to be considered if this approach is used for drug screening on acute muscle strength23. Third, the use of maximal electrical nerve stimulation deviates from voluntary recruitment strategies, as discussed above, and therefore does not reflect changes in strength that may be due to neuromuscular recruitment adaptations.
In vivo torque measurements may also be limited with regard to establishing a specific mechanism for experimental observations. For instance, torque about the ankle joint depends not only on muscle force production but also on the tendon and joint and connective tissue properties. Moreover, force is generated by groups of muscles, specifically the plantar flexors (gastrocnemius, soleus, and plantaris muscles) and the dorsiflexors (peroneus tertius, tibialis, and digitorum muscles) in pigs. Therefore, interpretations of maximal in vivo torque data require consideration of potential musculotendinous and anatomical changes and are limited to muscle groups, not individual muscles. Relatedly, muscle groups are often made up of a mixture of predominately fast- and slow-muscle fibers, such as the gastrocnemius and soleus muscle, respectively, of the plantar flexors. Contractile properties such as rate of contraction and relaxation (or time-to-peak contraction and half-relaxation time) are not reliable indicators of fiber type physiology using in vivo versus isolated muscle preparations, such as in vitro or in situ testing protocols47. Isolated muscle preparations are also superior in understanding the influence of biomechanical parameters on muscle function because properties such as muscle length can be precisely controlled; it is important to stress that the joint angle-torque relationship is not directly equivalent to the muscle length-force relationship, as tendon (e.g., slack), muscle (e.g., pennation angle, sarcomere overlap), and joint (e.g., moment arm) properties that contribute to torque production are dependent on the joint angle. To that end, large animal in situ functional testing48 could be a valuable addition to in vivo testing, bearing in mind that in situ testing is a terminal experiment. Other advancements to the current protocol that may be explored in the future to improve the mechanistic insight of experimental findings include using ultrasound B-mode imaging to measure muscle and tendon architectural properties and implantation of a tendon force transducer to measure muscle force during voluntary and electrically stimulated contractions49.
Importance and potential applications of the method
This protocol evaluates in vivo torque-producing capacity of the porcine dorsiflexor muscle group, demonstrating a non-invasive method to assess gain or loss of muscle function in a physiological setting. Because the methodology is non-terminal to the pig, it can also be used to evaluate muscle function in the same subjects longitudinally during the progression of a disease, or before, during, and following a treatment strategy. As such, a repeated measures experimental design may allow for robust statistical comparisons with greater power and fewer animals compared to independent measures. Additionally, skeletal muscle dysfunction is a salient component of various disease processes and conditions, such as chronic disease-associated muscle wasting (e.g., heart failure, kidney failure, AIDS, cancer, etc.), muscular dystrophy, neurodegenerative diseases (e.g., SMA or amyotrophic lateral sclerosis; ALS), aging (i.e., sarcopenia), and drug toxicities. Skeletal muscle functional capacity is a critical primary outcome measure for interventions such as exercise, nutrition, and drug and regenerative medicine therapies. Thus, the protocol described herein to reliably evaluate porcine torque producing capacity in vivo may be used across numerous study applications. It may be instrumental in acquiring extensive animal data for the translation of developing therapies.
The authors have nothing to disclose.
Work and data presented were supported broadly by the US Army Medical Research and Material Command to BTC and SMG (#MR140099; #C_003_2015_USAISR; #C_001_2018_USAISR); and the Department of Veterans Affairs, Veterans Health Administration, Office of Research and Development (I21 RX003188) to JAC and Dr. Luke Brewster. The authors gratefully acknowledge the USAISR Veterinary Service and Comparative Pathology Branches and UMN Advanced Preclinical Imaging Center for technical assistance in completing these studies.
615A Dynamic Muscle Control LabBook and Analysis Software Suite | Aurora Scientific Inc. | 615A | Compatible Win Vista/7/10 |
892A Swine Isometric Footplate Test Apparatus | Aurora Scientific Inc. | 892A | Includes Isometric Load Cell, Pig Footplate, Goniometer stage and positioners |
Calibration Weights | Ohaus or similar | 80850116 | |
Computer | Aurora Scientific or any vendor | 601A | Computer must include data acquisition card and interface for software |
Gauze pad | Various vendors | 4 by 4 squares or similar | |
Monopolar Needle Electrodes | Chalgren, Electrode Store, or similar vendor | 242-550-24TP, or DTM-2.00SAF | |
Non-adhesive Flexiable Tape | 3M, Coflex, or similar | 4 inch by 5 yard role | |
Stimulator | Aurora Scientific or comparable | 701C | Must include constant current stimulation mode |