The purpose of this publication is to present our original work on a multi-muscle surface electromyographic approach to quantitatively characterize respiratory muscle activation patterns in individuals with chronic spinal cord injury using vector-based analysis.
During breathing, activation of respiratory muscles is coordinated by integrated input from the brain, brainstem, and spinal cord. When this coordination is disrupted by spinal cord injury (SCI), control of respiratory muscles innervated below the injury level is compromised1,2 leading to respiratory muscle dysfunction and pulmonary complications. These conditions are among the leading causes of death in patients with SCI3. Standard pulmonary function tests that assess respiratory motor function include spirometrical and maximum airway pressure outcomes: Forced Vital Capacity (FVC), Forced Expiratory Volume in one second (FEV1), Maximal Inspiratory Pressure (PImax) and Maximal Expiratory Pressure (PEmax)4,5. These values provide indirect measurements of respiratory muscle performance6. In clinical practice and research, a surface electromyography (sEMG) recorded from respiratory muscles can be used to assess respiratory motor function and help to diagnose neuromuscular pathology. However, variability in the sEMG amplitude inhibits efforts to develop objective and direct measures of respiratory motor function6. Based on a multi-muscle sEMG approach to characterize motor control of limb muscles7, known as the voluntary response index (VRI)8, we developed an analytical tool to characterize respiratory motor control directly from sEMG data recorded from multiple respiratory muscles during the voluntary respiratory tasks. We have termed this the Respiratory Motor Control Assessment (RMCA)9. This vector analysis method quantifies the amount and distribution of activity across muscles and presents it in the form of an index that relates the degree to which sEMG output within a test-subject resembles that from a group of healthy (non-injured) controls. The resulting index value has been shown to have high face validity, sensitivity and specificity9-11. We showed previously9 that the RMCA outcomes significantly correlate with levels of SCI and pulmonary function measures. We are presenting here the method to quantitatively compare post-spinal cord injury respiratory multi-muscle activation patterns to those of healthy individuals.
1. Settings
2. RMCA Protocol
3. Data Analysis
Figure 3 represents the electromyogram and airway pressure (on top) simultaneously recorded during MEPT from a non-injured (left) and SCI (right) individuals. Note decreased airway pressure and absence of sEMG activity in expiratory muscles in an SCI subject when compared to a non-injured individual (marked with gray ellipses). Note also that start of the task, as marked on the bottom, is associated with increased sEMG activity and raising airway pressure.
Figure 4 highlights the main steps of constructing the RV. The beginning and ending of the task (event window) was defined as data points (Step 1) according to the marker. The Root Mean Square (RMS) of sEMG within this event window represents mean sEMG activity for each muscle (Step 2). The RV is assembled using RMSs values for the specific muscle combination (Step 3).
Figure 5 illustrates calculation of the PRV and its Magnitude for a group of non-injured (healthy) individuals. Prototype response matrix was constructed using individual RVs (Step 4). Each column in the prototype response matrix includes data for each individual in the group (n=1,2,..,N) and each row represents quantified sEMG activity (RMSs) of specific muscle from the all individuals in the group. The PRV was calculated by taking the average of each row of the prototype response matrix (Step 5). The Magnitude value represents the length of the RV and was calculated according to the formula shown (Step 6).
Figure 6 shows steps for the SI calculation. The RV (Step 7) and its Magnitude (Step 8) for particular SCI individual were calculated as also shown on Figures 4 and 5. The SI was obtained by calculating the inner product of PRV and RV (Step 9).
Figures 7 and 8 illustrate the construction of the vectors and calculations of the outcomes using real data. These calculations can be made by using any appropriate software tool like Mathlab, Excel or others.
Table 1. List of specific equipment and supplies used for the Respiratory Motor Control Assessment.
Figure 1. sEMG and airway pressure recording equipment used for the Respiratory Motor Control Assessment. Click here to view larger figure.
Figure 2. T-piece Monitoring Circuit with air tube assembling used to record the airway pressure during MEPT. Note that for the MIPT, the air leaking part should be flipped to the inspiratory side.
Figure 3. sEMG activity during Maximum Expiratory Pressure Task (MEPT) in a Non-injured Individual and a subject with Spinal Cord Injury (SCI). Air pressure developed is shown on top (PRESSURE) and simultaneously recorded sEMG activity with event mark (MARKER) on the bottom. Vertical gray lines represent the 5-sec analysis windows for VRI calculations. Note decreased airway pressure and absence of activity in expiratory muscles: Right (R) and Left (L) Intercostal (IC), Rectus Abdominus (RA), and Oblique Abdominus (O) in an SCI subject compared to a non-injured individual (marked with gray ellipses). Other muscles shown: Right (R) and Left (L) Sternocleidomastoid (SC), Scalene (S); Upper Trapezius (UT), Pectoralis (P), Diaphragm (D), Lower Trapezius (LT), and Paraspinal (PS). Click here to view larger figure.
Figure 4. Steps for the calculation of the Response Vector (RV). Note that the RV during the specific task was assembled using Root Mean Square (RMS) values calculated for specific muscles.
Figure 5. Steps for the calculation of the Prototype Response Vector (PVR). Note that the RV from each healthy individual in the group was used to create the PVR and calculate its Magnitude. Click here to view larger figure.
Figure 6. Steps for the calculation of the Similarity Index (SI) and Magnitude. Note that the SI was calculated using PRV and the RV obtained from SCI subject (sRV). Note also that the Magnitude was calculated as the length of the sRV. Click here to view larger figure.
Figure 7. Example of PRV calculations using data obtained during Maximum Expiratory Task (MEPT) in 17 non-injured individuals. Click here to view larger figure.
Figure 8. Example of the VRI calculations using data obtained during Maximum Expiratory Task (MEPT) in healthy (non-injured) and SCI individuals. Note that in contrast to the non-injured subject, absence and decreased expiratory muscles activity (IC, RA, and O) in the SCI individual decreased the Similarity Index (SI) value. Also note lower overall muscle activity in an SCI individual associated with lower Magnitude value. Click here to view larger figure.
Standard clinical tests to evaluate respiratory motor function after SCI and other disorders include the pulmonary function tests and the American Spinal Injury Association Impairment Scale (AIS) evaluation14,15. However, these tools are not designed for quantitative evaluation of the trunk and respiratory motor control. In our previously published work9, we have shown that the RMCA is a valid method to quantitatively evaluate the respiratory motor function affected by SCI. We have demonstrated that this method can be used despite test-re-test and subject-to-subject EMG amplitude variability.
In order to quantify the degree of difference (SI) in the multi-muscle distribution patterns produced by a tested individual (RV) against the normative vector (PRV), the PRV and RV can be constructed for any combination of muscles. The motor tasks can vary and depend on study design as well as the list of the muscles. In contrast to the SI, Magnitude values, representing the overall sEMG activity, can be modified by subject effort, compensatory muscle activity, and physical properties of body tissues.
The SI, while providing the quantitative measure of how close the multi-muscle activation pattern is to the normative pattern, does not describe in which way the pattern may be different. For this reason, it is essential to qualitatively describe changes in the pattern and individual muscle activation. In addition, further study is needed to examine sEMG signals for additional parameters with which to characterize multi- and single-muscle activation properties.
The method presented here provides a systematic way to evaluate motor control of trunk muscles used to perform respiratory motor tasks in comparison with normative standard measures. In addition to a descriptive sEMG evaluation to investigate in which way the muscle activation pattern is altered, the VRI calculations provide the unified index values from which the severity of disruption can be compared across individuals and changes over time can be monitored. This method allows the evaluation of the respiratory motor control status and the effects of existing and novel interventions on impaired respiratory motor control in persons with SCI and other disorder.
The authors have nothing to disclose.
This work was supported by Christopher and Dana Reeve Foundation (Grant CDRF OA2-0802-2), Kentucky Spinal Cord and Head Injury Research Trust (Grant 9-10A – KSCHIRT), Craig H. Neilsen Foundation (Grant 1000056824 – HN000PCG) and National Institutes of Health: National Heart Lung and Blood Institute (Grant 1R01HL103750-01A1).
Name of the Device | Vendor | Product # | Comments |
PowerLab System 16/35 | ADInstruments | PL3516 | Number of units depends on number of channels recorded |
EMG System MA 300 | Motion Lab Systems | MA300-XVI | Number of units depends on number of channels recorded |
Low Pressure Transducer MP45 | Validyne | MP45-40-871 | |
Basic Carrier Demodulator CD15 | Validyne | CD15-A-2-A-1 | |
Air Pressure Manometer | Boehringer | 4103 | Needed for MP45 calibration |
Event Marker | Hand held switch that when pressed gives a DC voltage and sound output (including 5-sec long mark) | ||
Alcohol Wipes | Henry Schein | 1173771 | Needed for electrodes placement |
Electrode Gel | Lectron II | 36-3000-25 | Needed for electrodes placement |
Tagaderm | Henry Schein | 7779152 | Needed for electrodes placement |
Noseclip | Henry Schein | 1089460 | |
T-piece Ventilator Monitoring Circuit with One-way Valves | Alleglance (Airlife) | 1504 | |
Air Tube | UnoMedical | 400E | |
Table 1. List of specific equipment and supplies used for the Respiratory Motor Control Assessment. |