Abstract
So far, surface electromyography (sEMG) has been the method of choice to detect and evaluate muscle fatigue. However, recent advancements in non-cryogenic quantum sensors, such as optically pumped magnetometers (OPMs), enable interesting possibilities to flexibly record biomagnetic signals. Yet, a magnetomyographic investigation of muscular fatigue is still missing. Here, we simultaneously used sEMG (4 surface electrode) and OPM-based magnetomyography (OPM-MMG, 4 sensors) to detect muscle fatigue during a 3 × 1-min isometric contractions of the left rectus femoris muscle in 7 healthy participants. Both signals exhibited the characteristic spectral compression distinctive for muscle fatigue. OPM-MMG and sEMG slope values, used to quantify the spectral compression of the signals, were positively correlated, displaying similarity between the techniques. Additionally, the analysis of the different components of the magnetic field vector enabled speculations regarding the propagation of the muscle action potentials (MAPs). Altogether these results show the feasibility of the magnetomyographic approach with OPMs and propose a potential alternative to sEMG for the study of muscle fatigue.
Introduction
Muscle fatigue is defined as a decrease in maximal force production in response to sustained or repetitive muscular activity (Enoka and Duchateau, 2008; Wan et al., 2017). The mechanisms underlying this process are manifold. One major contributor seems to be the accumulation of lactate in the muscles as a byproduct of anaerobic glycolysis during condition of high energy demand (Tesch et al., 1978; Sahlin, 1986; Ishii and Nishida, 2013); however, there is no consensus and existing evidence questions this assumption (Brooks, 2001; Gladden, 2004; Sola-Penna, 2008). Therefore, a global mechanism cannot be addressed (Allen et al., 2008), leading to a variety of available protocols and methods for the study of muscular fatigue (Vøllestad, 1997; Cifrek et al., 2009; Ma et al., 2009; Muramatsu and Kobayashi, 2014).
Surface electromyography (sEMG) is a commonly used and well-established means to continuously and non-invasively monitor muscle fatigue (Jurell, 1998; Cifrek et al., 2009; Marco et al., 2017; Liu et al., 2019) by recording the muscle action potentials (MAP) produced during muscle contraction. An early investigation of the frequency components of the electromyography (EMG) signal reported a shift of the power spectrum toward lower frequencies as an index of muscle fatigue (Kogi and Hakamada, 1962). A phenomenon also referred to as spectral compression (De Luca, 1993; Lowery et al., 2000). This result has been replicated several times (Viitasalo and Komi, 1977; Petrofsky et al., 1982; De Luca, 1984; Duchêne and Goubel, 1990; Doud and Walsh, 1995; Lowery et al., 2000; Beck et al., 2014) confirming the validity of this measure for the evaluation of muscle fatigue (Cifrek et al., 2009). Nevertheless, in spite of the recognized reliability of sEMG, the last few decades have seen the progressive development and miniaturization of magnetic sensors (Zuo et al., 2020), for the study of muscle activity.
Magnetomyography (MMG), as defined in the pioneering work by Cohen and Givler (1972), seeks to investigate muscle contraction physiology by recording the magnetic fields which arise from the same ionic current that generate the EMG signal (Cohen and Givler, 1972; Parker and Wikswo, 1997). The working principle of MMG is based on Biot-Savart law. The ionic current flowing within a muscle fiber generates concentric and circular magnetic fields, which can be recorded with a sufficiently sensitive magnetometer. Even though MMG and EMG signals originate from the same ionic currents, and present comparable temporal and spectral profile (Cohen and Givler, 1972; Parker and Wikswo, 1997), the potential superiority of MMG for the investigation of muscle contractions has been suggested (Garcia and Baffa, 2015). The magnetic permeability of human tissue is almost the same as that of empty space (Schenck, 1996); therefore, MMG signals are much less distorted by the different tissue layer interposed between the source and the skin surface. Besides, noise voltages originating at the electrode-skin interface during sEMG recording might interfere with the signal (Huigen et al., 2002; De Luca et al., 2010) and could be avoided using contactless MMG recording. Indeed, to date, MMG has shown a broad variety of applications (Fenici et al., 2005; Broser et al., 2018; Escalona-Vargas et al., 2019; Strand et al., 2019; Barchiesi et al., 2020; Blanks and Eswaran, 2020; Elzenheimer et al., 2020; Llinás et al., 2020; Marquetand et al., 2021; Broser et al., 2021a,b).
Here, we use MMG with Optically Pumped Magnetometers (OPMs) to monitor muscle fatigue of the rectus femoris muscle during an isometric contraction. The small size of the sensors and the unnecessity of cryogenic cooling (for a detailed technical description of OPM sensors see: Shah and Wakai, 2013; Boto et al., 2017; Tierney et al., 2019) allow a simultaneous Optically Pumped Magnetometer Magnetomyography (OPM-MMG) and sEMG recording. We compare the spectral profile of the signals over time using time-frequency decomposition and quantify the spectral compression by calculating the decrease of the spectral center of gravity over time. Additionally, we correlate OPM-MMG and sEMG slope values, used to quantify the spectral compression of the signals, to estimate the similarity between the techniques. Finally, we analyze the signal magnitude of the different components of the magnetic field vector to infer coarse spatial information regarding the main direction of propagation of the MAPs during isometric contraction.
Materials and Methods
Participants
Twelve healthy subjects participated in the study (6 males; mean age: 26.5 ± 2.9 SD years; and mean body-mass-index: 22.3 ± 1.5 SD kg/m2). The experiment was conducted at the MEG Center of the University of Tübingen (Germany) in March 2021, in accordance with the standards of the World Medical Association. All the participants are also authors of this publication and gave informed consent for their data to be published.
Experimental Setup
Prior to the measurement, ultrasound imaging (Mindray TE7, 14Mhz-linear probe) of the left rectus femoris muscle was performed for each subject by a certified ultrasound user (JM) to determine the longitudinal axis of the muscle. Following, participants sat down on a comfortable chair inside a magnetically shielded room (Ak3b, VAC Vacuumschmelze, Hanau, Germany), extending the legs on the chair at a controlled knee angle of 150°.
Four OPM sensors (QZFM-gen-1.5, QuSpin Inc., Louisville, CO, United States) were placed linearly in a distal to proximal order parallel to the longitudinal extent of the left rectus femoris muscle. The OPMs were held by a custom-built plastic holder attached to an aluminum structure. The sensors were placed at 6 cm distance from each other, with the first sensor placed 6 cm proximal to the patella, and 1–3 cm (mean distance 2.03 ± 0.65 SD cm) above the skin surface. In an alternating pattern, with respect to the OPM sensors, 4 paramagnetic surface electrodes (Conmed, Cleartrace2 MR-ECG-electrodes) were place along the rectus femoris muscle for the recording of the electromyographic signal (Figures 1B,D). Additionally, one ground and one reference electrode were placed, respectively, on the right shoulder and on the ipsilateral knee. As to the best of our knowledge, no paramagnetic disposable EMG electrodes are commercially available, we used the MR-ECG-electrodes and twirled a paramagnetic cable wire around the electrode button, as previously done in Marquetand et al. (2021).
Figure 1
During the measurement, the left ankle was pressed against the chair by a second person in the magnetically shielded room, and the participants were asked to push against it with the strongest, most continuous force possible. Three isometric contractions of the rectus femoris muscle were performed by the participants in 60 s blocks, with 30 s rest between each block (Figure 1A). Since the utilized OPMs could record only two orthogonal components of the magnetic field vector at a time (y and z, OPM sensor coordinates, Figure 1C), we repeated the measurements after rotating the sensors by 90° around their z-axis to record the magnetic field in all three spatial directions (X, Y, and Z, leg coordinates, Figures 1B,D). For this reason, the Z direction was recorded twice; therefore, we will refer to ZX and ZY to indicate the Z component orthogonal, respectively, to the X and Y directions. To avoid unbalanced fatigue effects, the starting orientation of the OPM sensors was alternated across the subjects. Following the same logic, we labeled the EMG recordings according to the OPMs orientation, thus respecting the same alternating pattern. For clarity, from now on whenever we refer to OPM vector components or spatial directions, we refer to the leg coordinates system.
Data Acquisition
MEG system (CTF Omega 275, Coquitlam, BC, Canada) electronics were used to record the analog output signal of the OPMs. The EEG channels of the same system were used to record the signals from the sEMG. Both OPM and EMG recordings were acquired with a sampling rate of 2343.8 Hz. The employed OPMs provide a magnetic field sensitivity in the order of 15 fT/ in a bandwidth of 3–135 Hz, an operating range below 200 nT, and a dynamic range of a few nanoteslas. To adapt to a non-zero magnetic background field, the sensors are equipped with internal compensation coils (Osborne et al., 2018) that can cancel magnetic background fields of up to 200 nT (operating range). For all measurements, the internal output gain factor of 3 was applied in the user interface of the sensors, which corresponds to a conversion factor for the analog output of the OPMs to magnetic flux of 1,11 nT/V.
Data Analysis
Data analysis was performed in MATLAB (Version R2018b, MathWorks Inc., Natick, Massachusetts United States) using custom-built scripts and the open-source FieldTrip toolbox (Oostenveld et al., 2010). Continuously recorded EMG and OPM data were segmented in 3 trials of 60 s length, accordingly with the 3 isometric contraction sessions. In all the trials, 1 s of data-padding remained at the beginning and end of each trial to avoid edge artifacts due to filtering processes. The signals were demeaned and filtered using a 10 Hz high-pass, zero-phase, sixth order Butterworth infinite impulse response (IIR) filter. To suppress the powerline noise, a band-stop (frequency ranges 48–52; 98–102; 148–152), zero-phase fourth order Butterworth IIR filter was applied. Time-frequency analysis was performed using Morlet wavelet spectral power decomposition, with a Gaussian width of 15 (number of cycles) to maximize the frequency resolution. Time window was set to 100 ms, and the frequencies of interest ranged from 20 to 90 Hz. The magnitude of the signal was computed taking the square root of the power and averaged across the subjects. For the comparison with sEMG, the orthogonal vector components Y−Z and X−Z at each OPMs were summed by applying the Pythagorean theorem so that as:
where represents the magnitude of the frequency ; x, y, and z are the 3 different directions for the OPM sensor. Accordingly, represents the resultant magnitude of directions X and ZX; represents the resultant magnitude of directions Y and ZY. The computed resultant components were further averaged across subjects.
To quantify the spectral compression of the EMG and MMG signals, the spectral center of gravity (Cg) was used. The spectral center of gravity, which has also been defined mean frequency, centroid, or central frequency (Farina and Merletti, 2000; Georgakis et al., 2003; Du and Vuskovic, 2004; Oskoei and Hu, 2008; Phinyomark et al., 2012), represents here the weighted average of the frequency, calculated as the sum of product of the EMG/MMG spectrum magnitude and the frequency divided by the total sum of the spectrum magnitude. The definition of Cg is given by:
where indicates the frequency bin, its respective magnitude, and is the length of the frequency bins. The time-frequency decomposition of the signal was done using a sliding time window of 100 ms length. The Cg was calculated for each time window to quantify the shift of the spectral magnitude over time.
Finally, an individual estimation of the average signal strength was computed. The magnitude of the 4 vector components (separately considering ZX and ZY) was averaged for each OPM sensor to estimate the mean strength of the signal at different sensor position regardless of the component. In the same way, the signal magnitude of each vector component was averaged across the 4 OPMs to obtain a mean estimation of the strength of each component independently from the sensor position.
Statistics
To test the hypothesis of a shift over time of the frequency spectrum as an indicator of muscle fatigue, the relation between time and Cg, in both sEMG and OPM signals, was estimated. After checking Cg values for normality using the Kolmogorov-Smirnov test, correlation (Pearson, two tailed) between Cg and time was calculated. A least-square regression line was fitted to the Cg data points to derive intercept and slope value as fatigue indices (Merletti et al., 1990; Farina and Merletti, 2000; Georgakis et al., 2003). Slope values of all subjects for sEMG and OPMs were further correlated (Pearson, one tailed) to quantify the similarity of the results between the two techniques. Finally, a comparison of all the individual averages signal magnitude, recorded at different sensors position and for different vector components, was performed. Repeated-measures ANOVA were used to assess the effect of sensors (within-subject factor, four levels: S1, S2, S3, and S4) and components (within-subject factor, four levels: X, ZX, Y, and ZY) related to the signal strength. Greenhouse–Geisser correction was applied to adjust for violation of sphericity assumption. Two-tailed, paired sample t-tests were used to post-hoc compare the results.
Results
Following the preprocessing, data were visually inspected. Five of the original twelve subjects had to be discharged due to a significant data loss caused by the partial detachment of the EMG reference electrode. The detachment was characterized by a clear flat line in the electromyographic signal which did not allow any further analysis. Hence, seven subjects were finally included in the analysis.
Time-Frequency Analysis and Spectral Center of Gravity in OPM-MMG
Spectral decomposition was performed on each component (X, Y, and Z directions) of the magnetic field vector. In all three spatial directions, the spectral shift of the signal magnitude toward lower frequencies was already visible in the time-frequency analysis (Figures 2A,C,E,G). Correlation (Pearson, two tailed) between spectral Cg and time highlights a negative relationship in all the spatial direction. However, differences between the vector components were appreciable, with the X direction displaying a stronger decay (Figure 2B, r = −0.69, p < 0.0001, slope = −0.0626 Hz/s, intercept = 54.6090 Hz) of the spectral Cg over time, compared with the Y direction (Figure 2F, r = −0.34, p < 0.0001, slope = −0.0220 Hz/s, intercept = 54.4718 Hz). Interestingly, the same pattern was founded also for the Z components, where the ZX component displayed a stronger decay (Figure 2D, r = −0.35, p < 0.0001, slope = −0.0242 Hz/s, intercept = 55.4986 Hz) compared to the ZY component (Figure 2H, r = −0.16, p = 0.0002, slope = −0.0101 Hz/s, intercept = 54.2747 Hz).
Figure 2
Comparison between sEMG and OPM-MMG
After combing the orthogonal components of the magnetic field vectors X−ZX and Y−ZY by applying the Pythagorean theorem, the computed resultant components named here X + ZX and Y + ZY were compared with the sEMG signals. Again the spectral shift was visible in the time-frequency analysis of the 2 signals (Figures 3A,C,E,G), and confirmed by the analysis of the spectral Cg. However, while the spectral compression depicted by the sEMG was similar between measurement 1 (Figure 3B, r = −0.52, p < 0.0001, slope = −0.0337 Hz/s, intercept = 53.7368 Hz) and measurement 2 (Figure 3D, r = −0.53, p < 0.0001, slope = −0.0359 Hz/s, intercept = 53.4894 Hz), OPM-MMG spectral compression varies between components: X + ZX (Figure 3F, r = −0.67, p < 0.0001, slope = −0.0505 Hz/s, intercept = 54.9551 Hz); Y + ZY (Figure 3H, r = − 0.24, p < 0.0001, slope = −0.0131 Hz/s, intercept = 54.3622 Hz).
Figure 3
To quantify the similarity of the results between the two techniques, OPM slope values of each subject, for each component (Supplementary Material), were correlated (Pearson, one tailed) with the slope values of the respective EMG measurement session (Figure 4). The values reported for the X component revealed a high degree of similarity with the corresponding EMG values (Figure 4A, r = 0.86, p = 0.0069), while for all the other components only a tendency toward a positive correlation was found (ZXFigure 4B, r = 0.66, p = 0.0529; Y Figure 4C, r = 0.67, p = 0.0512; and ZYFigure 4D, r = 0.66, p = 0.0543).
Figure 4
Comparison of the Signal Strength Between OPM Sensors Position and Vector Components
Considering the differences between vector components highlighted by the previous results, an analysis of the average signal strength measured at different sensors position for different components was performed. ANOVA on signal magnitude revealed a significant main effect of both sensors position F(3,18, ε = 0.471) = 6.241, p = 0.027 and vector components F(3,18, ε = 0.598) = 4.7, p = 0.037. No significant sensors position x components interaction was found F(9,54, ε = 0.164) = 1.317, p = 0.303. Post-hoc comparisons showed significant differences in signal magnitude between both component X (p = 0.021, two-tailed paired t-test) and ZY (p = 0.036, two-tailed paired t-test) compared to component Y (Figure 5A, and between sensor 1 and 4 p = 0.027, two-tailed paired t-test; Figure 5B).
Figure 5
Discussion
In this first magnetomyographic investigation of muscle fatigue with OPM sensors, we showed a spectral compression of the MMG signal, highlighted by the decrease of the spectral center of gravity over time. This compression was comparable with the one of the electromyographic signal that was simultaneously recorded through sEMG, and which is regarded as a distinctive marker of muscle fatigue (Kogi and Hakamada, 1962; Viitasalo and Komi, 1977; Petrofsky et al., 1982; De Luca, 1984; Duchêne and Goubel, 1990; Doud and Walsh, 1995; Lowery et al., 2000; Beck et al., 2014).
Interestingly, while for the EMG signal, the spectral compression was closely matched between measurements, MMG signal compression varied considerably depending on the different recorded components of the magnetic field vector (refer to leg coordinates system Figures 1B,D). The clearest spectral compression was found for the X component. This finding was also supported by a significant correlation between the single subjects’ slope values for the X component, with the values of the corresponding EMG measurement, which highlight the high similarity of the results reported by the two techniques. For all the other vector components (ZX, Y, and ZY), the picture is more blurred. Even though all the different signals displayed a significant spectral compression over time, the entity of this compression was lower compared to the one showed by the X component. Moreover, the correlation of the single subjects’ slope values with the corresponding EMG measurement reported only a tendency rather than a significant positive correlation. This was true not only for the components Y and ZY which were recorded on a different measurement after the rotation of the sensors, but also for the ZX component that is orthogonal to X. However, since the sample size sensibly influences the significance of the results even in presence of a strong correlation, the latter finding can be explained by the small number of subjects included in the experiment and can therefore be regarded as a limitation of the study. Yet, it is evident that different components of the MMG signal relate differentially to the EMG signal, enabling speculations on the major direction of the magnetic field.
As depicted in the model shown in Figure 5C, the strengths of the recorded vector components vary as a function of the position and orientation of the OPM sensors. Assuming a motor action potential traveling along the Y axis, a circular magnetic field in the X−Z plane is generated. In this situation, a sensor placed in position I would record only the X component of the magnetic field vector, a sensor placed in position II would capture a combination of the X and Z components, while a sensor placed in position III would measure exclusively the Z component. Since the magnetic field vector in the Y direction would be zero, no magnetic activity could be recorded in this direction regardless the orientation of the sensors. In our data (Figure 5A), the lowest signal strength was observed for the Y component, while components X and Z were similar. By comparing this finding with the model, we can speculate that the propagation of the MAPs occurs mainly along the Y axis. Considering the comparable signal strength of the X and Z components, we can reasonably exclude the possibilities described by position I and III where a low signal would be expected in two of the recorded components (Z and Y for position I; X and Y for position III). Hence, a situation like the one depicted by position II seems to be the one that better fit our data, in which the only one component displaying a low signal correspond to the propagation direction of the muscular action potentials. The interpretation of this result must be made cautiously since we can only speculate on the propagation direction of the MAPs, by looking at the major direction of the magnetic field. Any spatial information regarding the propagation of the signal at the source level requires the implementation of source localization algorithms that involve volume conduction models which are currently not yet developed for the rectus femoris muscle. Also, the comparison of the signal magnitude recorded at different OPM sensors (Figure 5B), highlighted a progressive increment of the strength of the signal between the sensors, with a significant difference between sensor 1 and 4. The understanding of this finding needs further investigations. One possible explanation might be related to the structure of the rectus femoris. The muscle originates from the anterior inferior iliac spine and the portion of the ilium above the acetabulum and develops with a fusiform shape until the insertion into the base of the patella. Due to its fusiform shape, the superficial surface of the rectus femoris increases from distal to proximal to decrease again at the origin. Considering the arrangement of the sensors in our study (Figure 5D), sensor 1 was in a distal position compared to sensor 4, and therefore, the higher signal magnitude recorded in the latter could reflect the overall muscle surface size. Another explanation for this finding, however, is the varying distance of the sensors from the source. Depending on the sensor position, the distance of the OPM from the skin surface, hence from the muscle, changed during the isometric contraction. Because of the decay of the magnetic field strength proportional to 1/r3, where r represents the distance from the source, the magnitude of the signal increases with the proximity to the source. Therefore, the lower signal strength recorded at sensor 1 might also reflect the larger distance of the sensor from the source, which decrease from distal to proximal resulting in a significant higher signal strength at sensor 4. Even though this was beyond the scope of this study, further investigations need to account for this factor. A possible solution would be to integrate a light distance detector to each sensor to real-time monitor the position of each OPM with respect to the skin surface.
Finally, some considerations regarding the potential of OPM-MMG to become a valid alternative to sEMG. Typically, sEMG used disposable silver/silver chloride (Ag/AgCl) electrodes which require a careful skin preparation (Hermens et al., 1999), and the application of a conductive gel to reduce the electrode-skin impedance and enhance conductivity. From a clinical perspective, this process is time consuming, it might be unpleasant for the patients since it requires a mild skin abrasion (Tam and Webster, 1977). And in some rare cases may even lead to adverse dermatological responses to electrolyte gels (Searle and Kirkup, 2000). The presence of skin rashes, scars, scabs, or areas with lesioned skin (e.g., burns or cuts) pose additional limitations to the use of sEMG. Furthermore, intrinsic noise sources, such as electro-chemical noise, originating at the electrode-skin interface may contaminate the sEMG signal (Huigen et al., 2002; De Luca et al., 2010). For all these limitations, contactless OPM-MMG measurement can offer an advantage. Yet, OPM technology is currently under development; therefore, the high cost and the requirement of magnetic shielding still limit its use to research context. However, in view of the promising advancement of wearable OPM systems for MEG application (Boto et al., 2019; Roberts et al., 2019; Boto et al., 2021; Seymour et al., 2021), and the progress in development of OPM-optimized magnetic shielding solutions (Holmes et al., 2018; Iivanainen et al., 2019; Borna et al., 2020; Zhang et al., 2020), it is foreseeable a broader diffusion of the technology in the future with a consequential decrease of the costs. This would paved the way toward the development of wearable OPM-MMG systems (Zuo et al., 2020), for the measurement and assessment of muscle activity for research and diagnostic purposes.
Conclusion
It can be summarized that is possible to detect muscle fatigue through magnetomyography with OPMs. OPM-MMG and sEMG slope values, which quantified the spectral compression of the signals over time as index of fatigue, were positively correlated exhibiting similarity between the techniques. Additionally, the capability of OPMs to record different components of the magnetic field vector enables speculations regarding the propagation direction of the MAPs during isometric contraction. Yet, these speculations need to be verified through the implementation of source reconstruction algorithms, which account for sensors position, muscle geometry, and volume conduction properties, allowing to estimate muscle activity at a source level. Overall, this first magnetomyographic investigation of muscular fatigue suggests that OPM-based MMG offers a promising new approach of study in muscle physiology.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The data supporting the findings of this study are available on request to the corresponding author.
Ethics statement
The study was performed in accordance to the Declaration of Helsinki (World Medical Association, 2001). The participants are also authors of this publication and gave informed consent for their data to be published.
Author contributions
JM and DSo: study design, data collection, data analysis, writing and revision. LS: writing and revision. SB: data analysis and revision. HC, GR, CK, MK, AR, JH, DSt, and JO: data collection and revision. TM: instrumentation, setup and revision. JD: setup, study design and data collection. CB and PB: data analysis and revision. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.724755/full#supplementary-material
Supplementary Material 1 | Single subjects’ spectral center of gravity, Pearson’s r, slope values, and point of intercept for the X and ZX components.
Supplementary Material 2 | Single subjects’ spectral center of gravity, Pearson’s r, slope values, and point of intercept for the Y and ZY components.
Supplementary Material 3 | Single subjects’ spectral center of gravity, Pearson’s r, slope values, and point of intercept for the two EMG measurements.
- OPM
Optically Pumped Magnetometer
- EMG
Electromyography
- sEMG
Surface Electromyography
- MMG
Magnetomyography
- OPM-MMG
Optically Pumped Magnetometer Magnetomyography
- MEG
Magnetoencephalography
- MAP
Muscular Action Potential
Abbreviations
References
1
AllenD. G.LambG. D.WesterbladH. (2008). Skeletal muscle fatigue: cellular mechanisms. Physiol. Rev.88, 287–332. doi: 10.1152/physrev.00015.2007
2
BarchiesiG.DemarchiG.WilhelmF. H.HauswaldA.SanchezG.WeiszN. (2020). Head magnetomyography (hMMG): A novel approach to monitor face and whole head muscular activity. Psychophysiology57:e13507. doi: 10.1111/psyp.13507
3
BeckT. W.StockM. S.DefreitasJ. M. (2014). Shifts in EMG spectral power during fatiguing dynamic contractions. Muscle Nerve50, 95–102. doi: 10.1002/mus.24098
4
BlanksA. M.EswaranH. (2020). Measurement of uterine electrophysiological activity. Curr. Opin. Physiol.13, 38–42. doi: 10.1016/j.cophys.2019.09.009
5
BornaA.CarterT. R.ColomboA. P.JauY.-Y.McKayJ.WeisendM.et al. (2020). Non-invasive functional-brain-imaging with an OPM-based magnetoencephalography system. PLoS One15:e0227684. doi: 10.1371/journal.pone.0227684
6
BotoE.HillR. M.ReaM.HolmesN.SeedatZ. A.LeggettJ.et al. (2021). Measuring functional connectivity with wearable MEG. NeuroImage230:117815. doi: 10.1016/j.neuroimage.2021.117815
7
BotoE.MeyerS. S.ShahV.AlemO.KnappeS.KrugerP.et al. (2017). A new generation of magnetoencephalography: room temperature measurements using optically-pumped magnetometers. NeuroImage149, 404–414. doi: 10.1016/j.neuroimage.2017.01.034
8
BotoE.SeedatZ. A.HolmesN.LeggettJ.HillR. M.RobertsG.et al. (2019). Wearable neuroimaging: combining and contrasting magnetoencephalography and electroencephalography. NeuroImage201:116099. doi: 10.1016/j.neuroimage.2019.116099
9
BrooksG. A. (2001). Lactate doesn’t necessarily cause fatigue: why are we surprised?J. Physiol.536:1. doi: 10.1111/j.1469-7793.2001.t01-1-00001.x
10
BroserP. J.KnappeS.KajalD.-S.NouryN.AlemO.ShahV.et al. (2018). Optically pumped magnetometers for magneto-Myography to study the innervation of the hand. IEEE Trans. Neural Syst. Rehabil. Eng.26, 2226–2230. doi: 10.1109/TNSRE.2018.2871947
11
BroserP. J.MarquetandJ.MiddelmannT.SomettiD.BraunC. (2021a). Investigation of the temporal and spatial dynamics of muscular action potentials through optically pumped magnetometers. J. Electromyogr. Kinesiol.59:102571. doi: 10.1016/j.jelekin.2021.102571
12
BroserP. J.MiddelmannT.SomettiD.BraunC. (2021b). Optically pumped magnetometers disclose magnetic field components of the muscular action potential. J. Electromyogr. Kinesiol.56:102490. doi: 10.1016/j.jelekin.2020.102490
13
CifrekM.MedvedV.TonkovićS.OstojićS. (2009). Surface EMG based muscle fatigue evaluation in biomechanics. Clin. Biomech.24, 327–340. doi: 10.1016/j.clinbiomech.2009.01.010
14
CohenD.GivlerE. (1972). Magnetomyography: magnetic fields around the human body produced by skeletal muscles. Appl. Phys. Lett.21, 114–116. doi: 10.1063/1.1654294
15
De LucaC. J. (1984). Myoelectrical manifestations of localized muscular fatigue in humans. Crit. Rev. Biomed. Eng.11, 251–279. PMID:
16
De LucaC. J. (1993). Spectral compression of the EMG signal as an index of muscle fatigue. In Electrophysiological Kinesiology (pp. 57–57) IOS Press.
17
De LucaC. J.Donald GilmoreL.KuznetsovM.RoyS. H. (2010). Filtering the surface EMG signal: movement artifact and baseline noise contamination. J. Biomech.43, 1573–1579. doi: 10.1016/j.jbiomech.2010.01.027
18
DoudJ. R.WalshJ. M. (1995). Muscle fatigue and muscle length interaction: effect on the EMG frequency components. Electromyogr. Clin. Neurophysiol.35, 331–339. PMID:
19
DuS.VuskovicM. (2004). “Temporal vs. spectral approach to feature extraction from prehensile EMG signals.” In Proceedings of the 2004 IEEE International Conference on Information Reuse and Integration, 2004. IRI 2004. (pp. 344–350) IEEE,
20
DuchêneJ.GoubelF. (1990). EMG spectral shift as an indicator of fatigability in an heterogeneous muscle group. Eur. J. Appl. Physiol.61, 81–87. doi: 10.1007/BF00236698
21
ElzenheimerE.LaufsH.Schulte-MattlerW.SchmidtG. (2020). Magnetic measurement of electrically evoked muscle responses With optically pumped magnetometers. IEEE Trans. Neural Syst. Rehabil. Eng.28, 756–765. doi: 10.1109/TNSRE.2020.2968148
22
EnokaR. M.DuchateauJ. (2008). Muscle fatigue: what, why and how it influences muscle function. J. Physiol.586, 11–23. doi: 10.1113/jphysiol.2007.139477
23
Escalona-VargasD.OliphantS.SiegelE. R.EswaranH. (2019). Characterizing pelvic floor muscles activities using magnetomyography. Neurourol. Urodyn.38, 151–157. doi: 10.1002/nau.23870
24
FarinaD.MerlettiR. (2000). Comparison of algorithms for estimation of EMG variables during voluntary isometric contractions. J. Electromyogr. Kinesiol.10, 337–349. doi: 10.1016/S1050-6411(00)00025-0
25
FeniciR.BrisindaD.MeloniA. M. (2005). Clinical application of magnetocardiography. Expert. Rev. Mol. Diagn.5, 291–313. doi: 10.1586/14737159.5.3.291
26
GarciaM. A.BaffaO. (2015). Magnetic fields from skeletal muscles: a valuable physiological measurement?Front. Physiol.6:228. doi: 10.3389/fphys.2015.00228
27
GeorgakisA.StergioulasL. K.GiakasG. (2003). Fatigue analysis of the surface EMG signal in isometric constant force contractions using the averaged instantaneous frequency. I.E.E.E. Trans. Biomed. Eng.50, 262–265. doi: 10.1109/TBME.2002.807641
28
GladdenL. B. (2004). Lactate metabolism: a new paradigm for the third millennium. J. Physiol.558, 5–30. doi: 10.1113/jphysiol.2003.058701
29
HermensH. J.FreriksB.MerlettiR.StegemanD.BlokJ.RauG.et al. (1999). European recommendations for surface electromyography. Roessingh Res. Dev.8, 13–54.
30
HolmesN.LeggettJ.BotoE.RobertsG.HillR. M.TierneyT. M.et al. (2018). A bi-planar coil system for nulling background magnetic fields in scalp mounted magnetoencephalography. NeuroImage181, 760–774. doi: 10.1016/j.neuroimage.2018.07.028
31
HuigenE.PeperA.GrimbergenC. A. (2002). Investigation into the origin of the noise of surface electrodes. Med. Biol. Eng. Comput.40, 332–338. doi: 10.1007/BF02344216
32
IivanainenJ.ZetterR.GrönM.HakkarainenK.ParkkonenL. (2019). On-scalp MEG system utilizing an actively shielded array of optically-pumped magnetometers. NeuroImage194, 244–258. doi: 10.1016/j.neuroimage.2019.03.022
33
IshiiH.NishidaY. (2013). Effect of lactate accumulation during exercise-induced muscle fatigue on the sensorimotor cortex. J. Phys. Ther. Sci.25, 1637–1642. doi: 10.1589/jpts.25.1637
34
JurellK. C. (1998). Surface Emg and fatigue. Phys. Med. Rehabil. Clin.9, 933–947. doi: 10.1016/S1047-9651(18)30242-0
35
KogiK.HakamadaT. (1962). Frequency analysis of the surface electromyogram in muscle fatigue. Rodo Kagaku J. Sci. Labour38, 519–528.
36
LiuS.-H.LinC.-B.ChenY.ChenW.HuangT.-S.HsuC.-Y. (2019). An EMG patch for the real-time monitoring of muscle-fatigue conditions During exercise. Sensors19:3108. doi: 10.3390/s19143108
37
LlinásR. R.UstininM.RykunovS.WaltonK. D.RabelloG. M.GarciaJ.et al. (2020). Noninvasive muscle activity imaging using magnetography. Proc. Natl. Acad. Sci.117, 4942–4947. doi: 10.1073/pnas.1913135117
38
LoweryM. M.VaughanC. L.NolanP. J.O’MalleyM. J. (2000). Spectral compression of the electromyographic signal due to decreasing muscle fiber conduction velocity. IEEE Trans. Rehabil. Eng.8, 353–361. doi: 10.1109/86.867877
39
MaL.ChablatD.BennisF.ZhangW. (2009). A new simple dynamic muscle fatigue model and its validation. Int. J. Ind. Ergon.39, 211–220. doi: 10.1016/j.ergon.2008.04.004
40
MarcoG.AlbertoB.TaianV. (2017). Surface EMG and muscle fatigue: multi-channel approaches to the study of myoelectric manifestations of muscle fatigue. Physiol. Meas.38, R27–R60. doi: 10.1088/1361-6579/aa60b9
41
MarquetandJ.MiddelmannT.DaxJ.BaekS.SomettiD.GrimmA.et al. (2021). Optically pumped magnetometers reveal fasciculations non-invasively. Clin. Neurophysiol.132, 2681–2684. doi: 10.1016/j.clinph.2021.06.009
42
MerlettiR.KnaflitzM.De LucaC. J. (1990). Myoelectric manifestations of fatigue in voluntary and electrically elicited contractions. J. Appl. Physiol.69, 1810–1820. doi: 10.1152/jappl.1990.69.5.1810
43
MuramatsuY.KobayashiH. (2014). Assessment of local muscle fatigue by NIRS–development and evaluation of muscle suit. Robomech J.1:19. doi: 10.1186/s40648-014-0019-2
44
OostenveldR.FriesP.MarisE.SchoffelenJ.-M. (2010). FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput. Intell. Neurosci.2011:e156869. doi: 10.1155/2011/156869
45
OsborneJ.OrtonJ.AlemO.ShahV. (2018). “Fully integrated standalone zero field optically pumped magnetometer for biomagnetism,” in Steep Dispersion Engineering and Opto-Atomic Precision Metrology XI (USA: SPIE), 89–95.
46
OskoeiM. A.HuH. (2008). Support vector machine-based classification scheme for myoelectric control applied to upper limb. I.E.E.E. Trans. Biomed. Eng.55, 1956–1965. doi: 10.1109/TBME.2008.919734
47
ParkerK. K.WikswoJ. P. (1997). A model of the magnetic fields created by single motor unit compound action potentials in skeletal muscle. I.E.E.E. Trans. Biomed. Eng.44, 948–957. doi: 10.1109/10.634647
48
PetrofskyJ. S.GlaserR. M.PhillipsC. A.LindA. R.WilliamsC. (1982). Evaluation of the amplitude and frequency components of the surface EMG as an index of muscle fatigue. Ergonomics25, 213–223. doi: 10.1080/00140138208924942
49
PhinyomarkA.ThongpanjaS.HuH.PhukpattaranontP.LimsakulC. (2012). “The usefulness of mean and median frequencies in electromyography analysis” in Computational Intelligence in Electromyography Analysis - A Perspective on Current Applications and Future Challenges. ed. NaikG. R. (Croatia: InTech), 195–220.
50
RobertsG.HolmesN.AlexanderN.BotoE.LeggettJ.HillR. M.et al. (2019). Towards OPM-MEG in a virtual reality environment. NeuroImage199, 408–417. doi: 10.1016/j.neuroimage.2019.06.010
51
SahlinK. (1986). Muscle fatigue and lactic acid accumulation. Acta Physiol. Scand. Suppl.556, 83–91. PMID:
52
SchenckJ. F. (1996). The role of magnetic susceptibility in magnetic resonance imaging: MRI magnetic compatibility of the first and second kinds. Med. Phys.23, 815–850. doi: 10.1118/1.597854
53
SearleA.KirkupL. (2000). A direct comparison of wet, dry and insulating bioelectric recording electrodes. Physiol. Meas.21, 271–283. doi: 10.1088/0967-3334/21/2/307
54
SeymourR. A.AlexanderN.MellorS.O’NeillG. C.TierneyT. M.BarnesG. R.et al. (2021). Using OPMs to measure neural activity in standing, mobile participants. NeuroImage244:118604. doi: 10.1016/j.neuroimage.2021.118604
55
ShahV. K.WakaiR. T. (2013). A compact, high performance atomic magnetometer for biomedical applications. Phys. Med. Biol.58, 8153–8161. doi: 10.1088/0031-9155/58/22/8153
56
Sola-PennaM. (2008). Metabolic regulation by lactate. IUBMB Life60, 605–608. doi: 10.1002/iub.97
57
StrandS.LutterW.StrasburgerJ. F.ShahV.BaffaO.WakaiR. T. (2019). Low-cost fetal Magnetocardiography: A comparison of superconducting quantum interference device and optically pumped magnetometers. J. Am. Heart Assoc.8:e013436. doi: 10.1161/JAHA.119.013436
58
TamH.WebsterJ. G. (1977). Minimizing electrode motion artifact by skin abrasion. I.E.E.E. Trans. Biomed. Eng.BME-24, 134–139. doi: 10.1109/TBME.1977.326117
59
TeschP.SjödinB.ThorstenssonA.KarlssonJ. (1978). Muscle fatigue and its relation to lactate accumulation and LDH activity in man. Acta Physiol. Scand.103, 413–420. doi: 10.1111/j.1748-1716.1978.tb06235.x
60
TierneyT. M.HolmesN.MellorS.LópezJ. D.RobertsG.HillR. M.et al. (2019). Optically pumped magnetometers: from quantum origins to multi-channel magnetoencephalography. NeuroImage199, 598–608. doi: 10.1016/j.neuroimage.2019.05.063
61
ViitasaloJ. H. T.KomiP. V. (1977). Signal characteristics of EMG during fatigue. Eur. J. Appl. Physiol.37, 111–121. doi: 10.1007/BF00421697
62
VøllestadN. K. (1997). Measurement of human muscle fatigue. J. Neurosci. Methods74, 219–227. doi: 10.1016/S0165-0270(97)02251-6
63
WanJ.QinZ.WangP.SunY.LiuX. (2017). Muscle fatigue: general understanding and treatment. Exp. Mol. Med.49, e384–e384. doi: 10.1038/emm.2017.194
64
ZhangX.ChenC.ZhangM.MaC.ZhangY.WangH.et al. (2020). Detection and analysis of MEG signals in occipital region with double-channel OPM sensors. J. Neurosci. Methods346:108948. doi: 10.1016/j.jneumeth.2020.108948
65
ZuoS.HeidariH.FarinaD.NazarpourK. (2020). Miniaturized magnetic sensors for implantable Magnetomyography. Adv. Mater. Technol.5:2000185. doi: 10.1002/admt.202000185
Summary
Keywords
OPM, sEMG, magnetomyography, muscle fatigue, quantum sensors
Citation
Sometti D, Semeia L, Baek S, Chen H, Righetti G, Dax J, Kronlage C, Kirchgässner M, Romano A, Heilos J, Staber D, Oppold J, Middelmann T, Braun C, Broser P and Marquetand J (2021) Muscle Fatigue Revisited – Insights From Optically Pumped Magnetometers. Front. Physiol. 12:724755. doi: 10.3389/fphys.2021.724755
Received
14 June 2021
Accepted
25 November 2021
Published
17 December 2021
Volume
12 - 2021
Edited by
Martina Krüger, Heinrich Heine University of Düsseldorf, Germany
Reviewed by
Mark Elisabeth Willems, University of Chichester, United Kingdom; Diana Escalona-Vargas, University of Arkansas for Medical Sciences, United States
Updates
Copyright
© 2021 Sometti, Semeia, Baek, Chen, Righetti, Dax, Kronlage, Kirchgässner, Romano, Heilos, Staber, Oppold, Middelmann, Braun, Broser and Marquetand.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Justus Marquetand, Justus.marquetand@med.uni-tuebingen.de
This article was submitted to Striated Muscle Physiology, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.