Kangyong Zheng
Changcheng Chen
Suyong Yang
Xueqiang Wang- 1Department of Sport Rehabilitation, Shanghai University of Sport, Shanghai, China
- 2Department of Rehabilitation Medicine, Qingtian People’s Hospital, Zhejiang, China
- 3School of Psychology, Shanghai University of Sport, Shanghai, China
- 4Department of Rehabilitation Medicine, Shanghai Shangti Orthopaedic Hospital, Shanghai, China
In this study, electroencephalography (EEG) was utilized to explore the neurophysiological mechanisms of aerobic exercise-induced hypoalgesia (EIH) and provide a theoretical basis for the application of aerobic exercise in pain assessment and treatment. Forty-five healthy subjects were randomly divided into moderate-intensity aerobic exercise [70% heart rate reserve (HRR)], low-intensity aerobic exercise (50% HRR), or control groups (sitting). Aerobic exercise was performed with cycling. Pressure pain threshold (PPT), heat pain threshold (HPT), event-related potential (ERP) induced by contact heat stimulus and pain scoring were measured before and after the intervention. We found that moderate-intensity aerobic exercise can increase the PPT (rectus femoris: t = −2.71, p = 0.017; tibialis anterior muscle: t = −2.36, p = 0.033) and HPT (tibialis anterior muscle: t = −2.219, p = 0.044) of proximal intervention sites rather than distal sites, and decreased pain scorings of contact heat stimulus. After moderate-intensity aerobic exercise, alpha oscillation power reflecting the central descending inhibitory function was enhanced (t = −2.31, p < 0.05). Low-intensity aerobic exercise mainly reduced the pain unpleasantness rating (Block 1: t = 2.415, p = 0.030; Block 2: t = 3.287, p = 0.005; Block 4: t = 2.646, p = 0.019; Block 5: t = 2.567, p = 0.022). Aerobic exercise had an overall EIH effect. Its hypoalgesic effect was related to exercise intensity and affected by the site and type of pain stimulus. Moderate-intensity aerobic exercise effectively reduced the sensitivity to various painful stimuli, and low-intensity aerobic exercise selectively inhibited the negative emotional pain response. The hypoalgesic mechanism of aerobic exercise involves the enhancement of the central descending inhibitory function.
Introduction
As a non-drug treatment, exercise has been widely used for chronic pain management (Brosseau et al., 2008; van Middelkoop et al., 2010; Chen et al., 2012; Loew et al., 2012; Steele et al., 2020). Exercise reduces the sensitivity to noxious stimuli through a process called exercise-induced hypoalgesia (EIH), which is an endogenous pain regulation (Thorén et al., 1990). EIH is manifested as an increased pain perception threshold and pain tolerance, and a change in pain ratings during or after exercise (Koltyn, 2000).
Existing research has shown that three types of exercise produce EIH effects, and these types are aerobic, dynamic resistance, and isometric contraction exercise (Naugle et al., 2012). Compared with the two other types of exercise, aerobic exercise is more diversified and is thus conducive to personalized exercise programs. However, the factors that affect aerobic EIH are complex; in particular, the influence of exercise intensity remains unclear (Polaski et al., 2019). Early studies have reported that exercise intensity must exceed 70% maximal oxygen uptake [VO (2) max], and as the intensity increases, the effect of aerobic EIH increases further (Kemppainen et al., 1990). The dependence of aerobic EIH effects on high intensity may limit the clinical application of aerobic exercise (Naugle et al., 2012). Given that pain is a dynamic process of physiology and psychology, accurately assessing the pain sensitivity is difficult. Therefore, using accurate sensory assessment methods is conducive to identifying the aerobic EIH effect. Recently, some researchers provided evidence that aerobic exercise intensities of 70% HRR and 50% HRR reduce pain ratings, especially the former for a greater dose–response effect (i.e., the more intense exercise produces larger effects) (Naugle et al., 2014).
Researchers have recently discovered that the aerobic EIH effect also has a clinical evaluation value. In some patients with chronic musculoskeletal pain, the effect of aerobic EIH is weakened and even reversed (i.e., increased pain sensitivity) for potential mechanisms of abnormal descending inhibition or excessive activation of muscle nociceptive afferents (Vierck et al., 2001; Staud et al., 2005; Naugle et al., 2012). High pain sensitivity is an important feature and risk factor of chronic pain (Meints et al., 2019; Nahman-Averbuch et al., 2019). In addition, the EIH efficiency before total knee arthroplasty is positively correlated with pain relief after surgery (Vaegter et al., 2017a). Therefore, evaluations of EIH are helpful in monitoring and evaluating the endogenous pain regulation system function (Gomolka et al., 2019).
Currently, the exercise programs concerning aerobic EIH have contained various exercise types (e.g., cycling, running), pain induction techniques (e.g., pressure and heat stimuli) and measurements [e.g., functional magnetic resonance imaging, event-related potential (ERP)] (Janal et al., 1984; Scheef et al., 2012; Jones et al., 2016; Vaegter et al., 2017b). Methodological differences among studies make it difficult to identify the effect of exercise on different pain aspects (Naugle et al., 2012, 2014; Jones et al., 2019). For example, it was reported that parameters of the exercise, like types, durations, and intensities of exercise, may determine which system is activated (Naugle et al., 2012).
Given that the exact neural mechanisms are still unclear, some studies have been proposed to elucidate EIH (Micalos et al., 2014; Jones et al., 2016). A recent study evaluated the impact of cycling on brain activation to pain in fibromyalgia and found that cycling seemed to activate brain areas, like the left dorsalateral prefrontal lobe of the anterior insula, which involved in descending pain inhibition, decreasing pain sensitivity (Ellingson et al., 2016). A novel ERP study found that the effect of exercise on the amplitude of somatosensory and laser evoked potentials was negligible when compared to that of controlled condition, yet failing to identify the straightforward mechanism of EIH (Jones et al., 2016). Perhaps a potential mechanism is the activation of endogenous opioid system, which links to changes in pain sensitivity (Naugle et al., 2012). Additionally, a plausible basis for EIH is by a functional restoration of the descending pain-inhibition pathways and/or desensitization (Naugle and Riley, 2014; Micalos et al., 2015). However, human research has not provided consistent evidences (Naugle et al., 2012). Further research should be conducted on the hypoalgesia effect and central nervous system mechanisms of moderate- and low-intensity aerobic exercise to promote the application of aerobic exercise in pain management.
This study evaluated the hypoalgesic effect of aerobic exercise among healthy subjects. The main research objectives are to (i) explore the influence of low- and moderate-intensity aerobic exercise on pain sensitivity and (ii) examine the neurophysiological mechanisms of aerobic EIH.
Materials and Methods
Participants
Forty-five healthy subjects (23 males and 22 females) were recruited online and via posters. The average age of the subjects was 24.47 years old. The inclusion criteria were as follows: 18–35 years old; right-handed; without non-persistent or intermittent pain in the last 3 months; no smoking and drinking habits; good health; non-professional athlete. The exclusion criteria were as follows: has participated in the same type of experiment in the past week; women who are menstruating; has a cardiovascular disease, such as hypertension, arrhythmia, and cardiomyopathy; has a history of neurological diseases, such as epilepsy, cerebral palsy, and spinal cord injury; has a disease that affects the musculoskeletal system, such as arthritis, tendonitis, and disc herniation; long-term use of drugs that affect the nervous and musculoskeletal systems; and presence of other acute symptoms, such as colds, fever, and cough, on the day of the experiment.
This study was approved by the Ethics Committee of Shanghai University of Sport. All subjects signed an informed consent form before the test. Meanwhile, demographics and clinical characteristics of the subjects were collected, including age, height, weight, grip strength, body mass index (BMI), sex ratio, years of education, and physical activity levels. Physical activity levels were assessed with the International Physical Activity Questionnaire–Short (IPSQ-S) which is a widely used self-reporting questionnaire for assessing physical activity (van Poppel et al., 2010; Lee et al., 2011). For this questionnaire, subjects were asked about the weekly frequency of walking, moderate intensity and high intensity activities in the past 7 days and the cumulative time per day, gaining an overall estimate of physical activity (Craig et al., 2003).
Interventions
This study was conducted in the psychological experiment center of Shanghai University of Sport. The experimental process is shown in Figure 1. Forty-five healthy subjects were recruited and randomly (1:1:1) divided into three groups via computer generated sequence, namely, moderate-intensity exercise group (Group 1), low-intensity exercise group (Group 2), and control group (Group 3). Aerobic exercise was individually performed in the form of cycling (Monark, Switzerland), and an H10 Polar heart rate belt (Polar, Finland) was used to measure the rest and exercise heart rates. Groups 1 and 2 received a power bicycle exercise intervention for 25 min (5 min of warm-up and 20 min of training). The intervention time for Group 3 was also 25 min (2 min of warm up and 23 min of rest). The exercise intensity for Groups 1–3 was 70% HRR, 50% HRR, and sitting for rest, respectively (American College of Sports Medicine Position Stand, 1998; Duncan et al., 2013). The subjects were instructed to report their rating of perceived exertion (RPE) scores, and their heart rates were recorded every 5 min. A Borg 6–20 RPE scale, a subjective rating of the intensity of exercise based on the feeling of exertion, ranging from 6 “no exertion at all” to 20 “maximal exertion,” was clearly visible to the subjects during the exercise (Williams, 2017).
Figure 1. Experimental flow chart. PPT, pressure pain threshold; HPT, heat pain threshold; CHEPs, contact heat-evoked potentials; HRR, heart rate reserve.
Test Methods
The three groups of subjects were tested for pressure pain threshold (PPT), heat pain threshold (HPT), and ERP induced by contact heat stimulus before and after the intervention.
Pressure Pain Threshold Test Method
PPT was tested using a pressure algometer (Wagner FPX, United States) with a probe diameter of 1.0 cm that exerts pressure up to 20 kg at a rate of 0.5 kg/s (Gerhardt et al., 2017). Four test sites on the right side were tibialis anterior muscle (one-half of the distance between the tibial tuberosity and lateral malleolus), rectus femoris (10 cm above the upper margin of the patella), sacrospinalis (2 cm away from the spinous process of lumbar 3), and thenar muscle (the thenar eminence of the hand). Tibialis anterior, one of the most important muscles during bicycling, did not change moving range dramatically during bicycling, allowing stable and reliable pain assessments (Blake et al., 2012; Jobson et al., 2013; Flouris et al., 2015; Jongerius et al., 2021). Three marks were placed approximately 1 inch apart over test muscle to ensure similar test point before and after inventions. As soon as the subject felt pain and said “pain”, the tester stopped pressing immediately and the value displayed on the algometer was PPT. The test was performed three times for each site with a time gap of 10 seconds, and the average value was obtained.
Heat Pain Threshold Test Method
HPT test was performed on three sites on the right side by using the PATHWAY sensory evaluation system (Medoc Ltd., Israel). The three sites were the forearm (5 cm above the volar wrist transverse striation), rectus femoris (10 cm above the upper margin of the patella), and tibialis anterior muscle (one-half of the distance between the tibial tuberosity and lateral malleolus). Three marks were placed approximately 1 inch apart over each test muscle to ensure similar test point before and after inventions. A slow-heating probe (AST stimulator) was used at a baseline temperature of 36 °C, heating rate of 1 °C/s and cooling rate of 3 °C/s. HPT was determined as the temperature when subjects initially felt pain from the rising temperature and pressed a button to terminate the rise. The test was performed three times for each site with a time gap of 5 s, and the average value was obtained.
Contact Heat Stimulus Parameter Setting and Pain Scoring
A PATHWAY sensory evaluation system with a fast-heating probe was used. The stimulation intensity was 50°C the heating speed was 70°C/s, and the cooling speed was 40°C/s (Meng et al., 2013). The stimulation sites were tibialis anterior muscles (one-half of the distance between the tibial tuberosity and lateral malleolus) of both legs. The interval between two stimuli was 18–20 s, and 5 stimuli formed a block. After a block ended, switched to the opposite side, alternating left and right. Each subject received contact heat stimuli in 6 blocks and 30 stimuli in total. About 3 s after each stimulation, the subjects were instructed to orally report pain intensity (no pain—worst pain imaginable) and unpleasantness (no unpleasantness—most unpleasantness) on numerical rating scale (NRS) that ranged from 0 to 10. The pain intensity and pain unpleasantness of each block were averaged, respectively, and the average values were obtained (Li et al., 2020).
Electroencephalography Data Collection and Processing
Electroencephalography Data Collection
The Brain Products system (BP, Germany) with a band pass filter of 0.01–100 Hz and sampling frequency of 1,000 Hz was used for data collection. EEG data were recorded by 64 Ag–AgCl scalp electrodes placed in accordance with the International 10–20 System. FCz was the reference (Lei and Liao, 2017), and the ground electrode was AFz. Eye blinks and ocular movement signals were recorded with a vertical electrooculographic electrode placed 1 cm below the lower eyelid. The electrode impedances did not exceed 10 kΩ.
Electroencephalography Data Preprocessing
The BP Analyzer 2.1 software was used to preprocess the EEG data. First, the data were re-referenced to TP9/TP10. S, 1–30 Hz band-pass filtering was performed. Third, eye blink and movement signals were removed through independent component analysis. Fourth, the EEG data were segmented into epochs with a time window of 1,500 ms (ranging from 500 ms pre-stimulus to 1,000 ms post-stimulus), and baseline correction was performed with the pre-stimulus interval (–500 ms to 0 ms) (Hu et al., 2014). Lastly, Epochs exceeding ± 80 μV were rejected as considered contaminated by artifacts.
Time-Domain Analysis
For ERP, the peak latencies and amplitudes of N2 and P2 waves were detected from the average waveform recorded by Cz electrode (Wydenkeller et al., 2008; Jones et al., 2016, 2019). The time for contact heat stimulus to reach the target temperature was about 200 ms. N2 and P2 waves were defined as the most negative and positive deflections between 350 and 700 ms after stimulus onset, respectively (Tu et al., 2016; Jin et al., 2018).
Frequency-Domain Analysis of Pre-stimulus Event-Related Potential Signal Oscillation
Fast Fourier transform (FFT) was applied to the ERP time domain signals before the stimulus to explore the influence of exercise on nerve oscillation at different frequencies. The procedure yielded an ERP spectrum ranging from 1 Hz to 30 Hz. Individual- and group-level ERP spectra were obtained after superposing and averaging.
Statistical Analysis
One-way ANOVA was applied for inter-group comparison of PPT, HPT and ERP results of pretest, posttest and pre–post differences (values before minus after intervention), respectively; and paired T-test was used for intra-group comparison. For the pain intensity and pain unpleasantness ratings, test time (pretest, posttest) × block (Blocks 1–6) × intervention condition (moderate-intensity exercise, low-intensity exercise, and control) were analyzed through repeated-measure ANOVA. Greenhouse-Geisser correction was applied when the spherical test was not met (Lo et al., 2019). Spearman correlation analysis was conducted on the three types of exercise intensity index (RPE score, exercise heart rate, and %HRR) with a difference in the various pain indicators before and after intervention. SPSS 22.0 was applied for statistical analysis, and p < 0.05 was considered statistically significant.
Results
Demographics and Clinical Characteristics and Exercise Intensity Monitoring
The forty-five subjects were randomly divided into three groups at a ratio of 1:1:1. Demographics and clinical characteristics, such as age, height, weight, resting heart rate, and physical activity level, showed no statistical difference among the groups (p > 0.05; Table 1). The differences in the exercise intensity monitoring results among groups, including heart rate,%HRR and RPE score during exercise, were all statistically significant (p < 0.001); that is, the highest values of heart rate,%HRR and RPE score were found in the moderate-intensity exercise group, and the lowest ones existed in control group.
Table 1. Demographics and clinical characteristics and exercise intensity monitoring of the subjects.
Behavioral Results
Pressure Pain Threshold Results
The test results of PPT at different sites before and after intervention for the three groups are shown in Figure 2. No statistical difference in PPT existed among the groups before exercise. After intervention, a statistically significant difference in thenar PPT was observed among the three groups (p = 0.011). The difference between Groups 1 and 3 was the most statistically significant (p = 0.004). The difference in PPT before and after intervention was compared at the group level and found to be not statistically significant (p > 0.05).
Figure 2. PPT in different sites before and after intervention in the three groups. (A) Thenar muscle PPT. (B) Sacrospinalis PPT. (C) Rectus femoris PPT. (D) Tibialis anterior PPT. On the left of the vertical dotted line are the PPT values before and after intervention of individuals (dots: purple for Group 1, green for Group 2, and red for Group 3) and groups (black line: X ± SD). The data on the right of the vertical dotted line are the individual (dot) and group (black line: X ± SD) values of the difference in PPT after intervention (△ = posttest–pretest). PPT: pressure pain threshold; Pre: pretest; Post: posttest. *p < 0.05; **p < 0.01.
The paired t-test showed that after intervention, the PPT of multiple sites in Group 1 increased significantly, including the PPT in rectus femoris (p = 0.017) and tibialis anterior muscle (p = 0.033).
Heat Pain Threshold Results
The results of HPT at different sites before and after intervention are shown in Figure 3. One-way ANOVA showed that no statistical difference in HPT in the forearm, rectus femoris and tibialis anterior muscle existed between groups before intervention, after intervention and their differences. The paired t-test revealed that the HPT of tibialis anterior muscle increased significantly after moderate-intensity aerobic exercise (p = 0.044).
Figure 3. HPT in different sites in the three groups before and after intervention. (A) Forearm HPT. (B) Rectus femoris HPT. (C) Tibialis anterior HPT. The data on the left side of the vertical dotted line are the HPT values of individuals (dots: purple for Group 1, green for Group 2, and red for Group 3) and groups (black line: X ± SD) before and after intervention. The data on the right side of the vertical dotted line are the values of the difference in HPT after intervention (△ = posttest–pretest) for individuals (dot) and groups (black line: X ± SD). HPT, heat pain threshold; Pre, pretest; Post, posttest; *p < 0.05.
Contact Heat Stimulus Results
(1) Pain intensity ratings (Figure 4). At the intra-group level, the main effect of test time (pretest and posttest) was significant (p = 0.002), and the main effect of the block (Blocks 1–6) was significant (p < 0.001). At the intergroup level, the main effect of intervention conditions (moderate-intensity exercise, low-intensity exercise, and control) was not significant (p = 0.259). The interaction effect between test time and intervention conditions was significant (p = 0.024).
Figure 4. Pain intensity ratings (X ± SD) for RPHS before and after intervention in the three groups. (A) Group 1. (B) Group 2. (C) Group 3. NRS, numeric rating scales; Pre, pretest; Post posttest. *p < 0.05; **p < 0.01.
The paired t-test indicated that the pain intensity ratings of contact heat stimuli in Blocks 2, 3, 4, 5, and 6 decreased significantly in Group 1 (Block 2: p = 0.031; Block 3: p < 0.001; Block 4: p = 0.014; Block 5: p = 0.003; Block 6: p < 0.001). The difference between the posttest and pretest exhibited an increasing trend, suggesting that the hypoalgesic effect of moderate-intensity aerobic exercise still existed after the end of contact heat stimulus.
(2) Pain unpleasantness ratings (Figure 5). At the intra-group level, the main effects of test time and block were significant (p < 0.001 for test time and p < 0.001 for block). At the inter-group level, the main effect of intervention conditions was not significant (p = 0.856). The interaction effect between test time and intervention conditions was not significant (p = 0.874).
Figure 5. Pain unpleasantness ratings (X ± SD) of RPHS before and after intervention in the three groups. (A) Group 1. (B) Group 2. (C) Group 3. NRS, numeric rating scales; Pre, pretest; Post, posttest. *p < 0.05; **p < 0.01.
The paired t-test revealed that the pain unpleasantness ratings in Blocks 3, 4, 5, and 6 decreased significantly in Group 1 (Block 3: p < 0.001; Block 4: p = 0.005; Block 5: p = 0.003; Block 6: p = 0.010). The pain unpleasantness ratings in Blocks 1, 2, 4, and 5 decreased significantly in Group 2 (Block 1: p = 0.030; Block 2: p = 0.005; Block 4: p = 0.019; Block 5: p = 0.022). The difference between the posttest and pretest showed no increasing trend, suggesting that moderate-intensity aerobic exercise had different durations in regulating the pain emotional response and pain perception. Groups 1 and 2 had different blocks of regulating effects, suggesting that the regulation of the pain emotional response after the two aerobic exercises had a time difference.
Event-Related Potential Results
Time-Domain Results
The ERP waveforms and topographic maps for the three groups before and after intervention are shown in Figure 6. The N2 and P2 amplitudes are presented in Table 2. No statistical significance existed in N2 and P2 amplitudes among the groups before and after intervention. At the intra-group level, the P2 amplitude of Group 1 after exercise was significantly lower than that before exercise (p < 0.05).
Figure 6. Contact heat-evoked potential waveforms and N2 and P2 topographic maps of the central parietal electrode (Cz) in the three groups before and after intervention. (A) Group 1. (B) Group 2. (C) Group 3. Pre, pretest; Post, posttest.
Table 2. N2 and P2 amplitudes before and after interventions [X ± SD (μV)].
Spectrum Results
The change in alpha oscillation power is shown in Table 3. The power density distribution and alpha band topographic maps before and after contact heat stimuli are shown in Figure 7. At the intra-group and inter-group levels, no statistical difference in the power density of the alpha band existed before and after intervention. The paired t-test showed that the power density of the alpha band in Group 1 increased significantly after intervention (p = 0.046).
Table 3. Comparison of power density of alpha oscillation before and after intervention in three groups [X ± SD (μV2/Hz)].
Figure 7. Power density distribution and alpha band topographic maps of the three groups before contact heat pain stimulation. (A) Group 1. (B) Group 2. (C) Group 3. (D) The topographic maps of alpha band oscillation power density before and after intervention. The alpha band (8–12Hz) is marked with a rectangular border. The alpha oscillation power density is largest in the frontal region. Pre, pretest; Post, posttest.
Correlation Results
As shown in Table 4, a Spearman correlation analysis was conducted on the three types of exercise intensity index (RPE score, exercise heart rate, and %HRR) with a difference in the various indicators before and after intervention. The difference between pain intensity rating and %HRR showed a significant negative correlation (r = –0.438, p = 0.007). Alpha power density was significantly and positively correlated with RPE score and exercise heart rate (r = 0.372, p = 0.020; r = 0.374, p = 0.019). The difference between the PPT of the thenar muscle and exercise heart rate had a significant positive correlation (r = 0.403, p = 0.012). A significant positive correlation also existed between the PPT difference of the tibialis anterior muscle and %HRR (r = 0.386, p = 0.016).
Table 4. The correlation between the three exercise intensity indexes with the difference of various indicators before and after intervention.
Discussion
Behavioral Evidence for Aerobic Exercise-Induced Hypoalgesia
The PPT results were consistent with those in previous studies. PPT and HPT increased after moderate-intensity aerobic exercise (Meeus et al., 2010; Kodesh and Weissman-Fogel, 2014), whereas low-intensity aerobic exercise had a limited effect on pressure pain sensitivity (Hoffman et al., 2004). The correlation analysis in this study indicated that the hypoalgesic effect of aerobic exercise was related to exercise intensity. The hypoalgesic effect may also be related to RPE score; thus, appropriate improvement of subjective fatigue within a certain range can enhance aerobic EIH (Peterson et al., 2019).
In terms of HPT, only the tibialis anterior muscle increased significantly after moderate-intensity aerobic exercise. No significant changes were observed in rectus femoris. Compared with mechanical stimulation, aerobic exercise exerted a limited influence on heat pain sensitivity (Jones et al., 2019). Previous research has shown that peripheral afference can partly explain EIH differences in different pain stimuli (de Souza et al., 2013). These differences are also related to the increase in skin temperature after aerobic exercise (Pertovaara et al., 1996). Likewise, an increase in local temperature can lead to heat pain sensitization. Cross-sectional studies have revealed that long-term vigorous physical exercise is associated with low pressure pain sensitivity and low heat pain sensitivity (Andrzejewski et al., 2010; Ellingson et al., 2012), which is considered to be the accumulation of the EIH effects of long-term exercise. Another possibility is that exercise exerts a greater impact on nociceptors located in the muscles than those in the skin. Pressure pain is likely to originate from deep nociceptors in the muscles, while heat pain is likely to originate from superficial nociceptors in the skin (Jones et al., 2019). Recent studies have found that blocking the blood flow of the limb during exercise can significantly reduce the EIH effect on the blocked limb (Jones et al., 2017). This condition may be related to the blocked transmission of hypoalgesic substances in the blood (Thorén et al., 1990). Further studies have shown that mechanical and heat pain are regulated by δ and μ opioid receptors, respectively (Scherrer et al., 2009). Meanwhile, cannabinoid receptors can regulate the input of mechanical and heat stimulation but have a greater impact on mechanical stimulation (Agarwal et al., 2007). Notably, a difference in test time increases the difference in test results. These explanations indicate that the differences in the EIH effects of different stimulus types still require further study.
In the present study, moderate-intensity aerobic exercise reduced the pain intensity and pain unpleasantness ratings in contact heat stimuli, whereas low-intensity aerobic exercise only influenced the pain unpleasantness ratings. These results indicate that the emotional response to pain is more affected by aerobic exercise than by pain perception. Some researchers used positron emission tomography to prove that pain stimulation, like hot water immersion stimulation, can selectively change the brain’s unpleasant response without changing the perception of pain intensity (Rainville et al., 1997). They also found that anterior cingulate gyrus neuron activity, encoding pain unpleasantness, is significantly changed, whereas the activation of the primary somatosensory cortex, encoding discriminative properties of somatosensory stimuli (i.e., intensity of noxious stimulation), remains unchanged; there seems to be at least a partial segregation of function between pain affect and sensation. These points were consistent with the findings of present study, for moderate-intensity aerobic exercise effectively reducing pain intensity and pain unpleasantness of contact heat stimulus and low-intensity aerobic exercise just selectively inhibiting the pain unpleasantness, for example.
The behavioral results of this study provided evidence that the EIH effect of moderate-intensity aerobic exercise is better than that of low-intensity aerobic exercise.
Event-Related Potential Evidence for Aerobic Exercise-Induced Hypoalgesia
The hypoalgesic effect of aerobic exercise was further supported by the ERP results, which exhibited a significant difference in P2 amplitude before and after exercise. P2 is one of the main components of pain-evoked potential, and its amplitude is closely related to stimulus intensity and the subject’s pain perception (Iannetti et al., 2006; Huang et al., 2013). P2 mainly originates from the anterior and medial cingulate gyrus and reflect the cognitive and emotional processing of pain perception (Garcia-Larrea et al., 2003; Valentini et al., 2012).
The spectrum analysis results showed that after moderate-intensity aerobic exercise, alpha oscillation power was enhanced before contact heat stimulus. The correlation analysis indicated that alpha oscillation power was positively correlated with exercise intensity. Alpha oscillations before stimulation can regulate pain perception after stimulation (Haegens et al., 2011; Tu et al., 2016). Brain oscillation before stimulation reflects the central preparation for external stimulation and predicts the neural response caused by perception and subsequent stimulation (Laufs et al., 2003; Boly et al., 2007; Gilbert and Sigman, 2007). Studies have shown that alpha oscillation is an indicator of sensory cortical excitability and attention resource allocation of the somatosensory system (Haegens et al., 2011). In general, alpha oscillation before stimulation partly reflects the activation of the resting state sensory-motor neural network (Anderson and Ding, 2011; Weisz et al., 2014). In addition, alpha oscillation reflects the central descending inhibitory function (Mathewson et al., 2011; Klimesch, 2012). Therefore, this study indicated that the hypoalgesic mechanisms of aerobic exercise may involve the enhancement of the central descending inhibitory function for a significant increase in the power density of alpha oscillation after moderate intensity aerobic exercise.
This study has its limitations. First, the blood lactate threshold of the subjects was not measured in this study, and the subjects had different tolerance to exercise, which may have influenced the aerobic EIH efficiency. Second, the pain test indicators were limited, and the experimental pain caused by contact heat stimulus could not fully reflect the neurophysiological characteristics of daily musculoskeletal pain. Lastly, the EIH efficiency of low-intensity aerobic exercise in this work was insufficient and still needs to be evaluated comprehensively by increasing the number of exercise interventions. The correlation between pain sensitivity and psychological factors has drawn increasing attention to whether psychological factors mediate the effects of EIH, and this area requires further research.
In summary, aerobic exercise exerts an overall EIH effect. Its hypoalgesic effect is related to exercise intensity and is affected by the site and type of pain stimulus. Moderate-intensity aerobic exercise effectively reduces the pain sensitivity to various painful stimuli, and low-intensity aerobic exercise selectively inhibits the negative emotional pain response. An increase in alpha oscillation power before a stimulus indicates that the hypoalgesic mechanisms of aerobic exercise may involve the enhancement of the central descending inhibitory function. This study focuses on the neurophysiological mechanisms of the aerobic EIH effect of moderate- and low-intensity aerobic exercise. The findings can provide theoretical guidance for optimizing pain exercise prescription in clinical practice.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Shanghai University of Sport. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
KZ, CC, SY, and XW conceptualized this study and contributed to revising and approving the final version of the manuscript. KZ and CC contributed to collecting data, analyzed the data, and drafted the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Science and Technology Innovation Program of Shanghai Science and Technology Committee (21S31902400), the Fok Ying-Tong Education Foundation of China (161092), the Scientific and Technological Research Program of the Shanghai Science and Technology Committee (19080503100), and the Shanghai Key Lab of Human Performance (Shanghai University of Sport) (11DZ2261100).
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.
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.
Acknowledgments
We thank all the healthy subjects who participated in this experiment and Shanghai University of Sport for providing the psychological laboratory.
References
Agarwal, N., Pacher, P., Tegeder, I., Amaya, F., Constantin, C. E., Brenner, G. J., et al. (2007). Cannabinoids mediate analgesia largely via peripheral type 1 cannabinoid receptors in nociceptors. Nat. Neurosci. 10, 870–879. doi: 10.1038/nn1916
American College of Sports Medicine Position Stand (1998). The recommended quantity and quality of exercise for developing and maintaining cardiorespiratory and muscular fitness, and flexibility in healthy adults. Med. Sci. Sports Exerc. 30, 975–991. doi: 10.1097/00005768-199806000-00032
Anderson, K. L., and Ding, M. (2011). Attentional modulation of the somatosensory mu rhythm. Neuroscience 180, 165–180. doi: 10.1016/j.neuroscience.2011.02.004
Andrzejewski, W., Kassolik, K., Brzozowski, M., and Cymer, K. (2010). The influence of age and physical activity on the pressure sensitivity of soft tissues of the musculoskeletal system. J. Bodyw. Mov. Ther. 14, 382–390. doi: 10.1016/j.jbmt.2009.07.004
Blake, O. M., Champoux, Y., and Wakeling, J. M. (2012). Muscle coordination patterns for efficient cycling. Med. Sci. Sports Exerc. 44, 926–938. doi: 10.1249/MSS.0b013e3182404d4b
Boly, M., Balteau, E., Schnakers, C., Degueldre, C., Moonen, G., Luxen, A., et al. (2007). Baseline brain activity fluctuations predict somatosensory perception in humans. Proc. Natl. Acad. Sci. U.S.A. 104, 12187–12192. doi: 10.1073/pnas.0611404104
Brosseau, L., Wells, G. A., Tugwell, P., Egan, M., Wilson, K. G., Dubouloz, C. J., et al. (2008). Ottawa Panel evidence-based clinical practice guidelines for aerobic fitness exercises in the management of fibromyalgia: part 1. Phys. Ther. 88, 857–871. doi: 10.2522/ptj.20070200
Chen, B. L., Zhong, Y., Huang, Y. L., Zeng, L. W., Li, Y. Q., Yang, X. X., et al. (2012). Systematic back muscle exercise after percutaneous vertebroplasty for spinal osteoporotic compression fracture patients: a randomized controlled trial. Clin. Rehabil. 26, 483–492. doi: 10.1177/0269215511423557
Craig, C. L., Marshall, A. L., Sjöström, M., Bauman, A. E., Booth, M. L., Ainsworth, B. E., et al. (2003). International physical activity questionnaire: 12-country reliability and validity. Med. Sci. Sports Exerc. 35, 1381–1395. doi: 10.1249/01.Mss.0000078924.61453.Fb
de Souza, G. G., Duarte, I. D., and de Castro Perez, A. (2013). Differential involvement of central and peripheral α2 adrenoreceptors in the antinociception induced by aerobic and resistance exercise. Anesth. Analg. 116, 703–711. doi: 10.1213/ANE.0b013e31827ab6e4
Duncan, M., Smith, M., and Lyons, M. (2013). The effect of exercise intensity on coincidence anticipation performance at different stimulus speeds. Eur. J. Sport Sci. 13, 559–566. doi: 10.1080/17461391.2012.752039
Ellingson, L. D., Colbert, L. H., and Cook, D. B. (2012). Physical activity is related to pain sensitivity in healthy women. Med. Sci. Sports Exerc. 44, 1401–1406. doi: 10.1249/MSS.0b013e318248f648
Ellingson, L. D., Stegner, A. J., Schwabacher, I. J., Koltyn, K. F., and Cook, D. B. (2016). Exercise strengthens central nervous system modulation of pain in fibromyalgia. Brain Sci. 6:8. doi: 10.3390/brainsci6010008
Flouris, A. D., Dinas, P. C., Tsitoglou, K., Patramani, I., Koutedakis, Y., and Kenny, G. P. (2015). Non-invasive measurement of tibialis anterior muscle temperature during rest, cycling exercise and post-exercise recovery. Physiol. Meas. 36, N103–N113. doi: 10.1088/0967-3334/36/7/n103
Garcia-Larrea, L., Frot, M., and Valeriani, M. (2003). Brain generators of laser-evoked potentials: from dipoles to functional significance. Neurophysiol. Clin. 33, 279–292. doi: 10.1016/j.neucli.2003.10.008
Gerhardt, A., Eich, W., Treede, R. D., and Tesarz, J. (2017). Conditioned pain modulation in patients with nonspecific chronic back pain with chronic local pain, chronic widespread pain, and fibromyalgia. Pain 158, 430–439. doi: 10.1097/j.pain.0000000000000777
Gilbert, C. D., and Sigman, M. (2007). Brain states: top-down influences in sensory processing. Neuron 54, 677–696. doi: 10.1016/j.neuron.2007.05.019
Gomolka, S., Vaegter, H. B., Nijs, J., Meeus, M., Gajsar, H., Hasenbring, M. I., et al. (2019). Assessing endogenous pain inhibition: test-retest reliability of exercise-induced hypoalgesia in local and remote body parts after aerobic cycling. Pain Med. 20, 2272–2282. doi: 10.1093/pm/pnz131
Haegens, S., Händel, B. F., and Jensen, O. (2011). Top-down controlled alpha band activity in somatosensory areas determines behavioral performance in a discrimination task. J. Neurosci. 31, 5197–5204. doi: 10.1523/jneurosci.5199-10.2011
Hoffman, M. D., Shepanski, M. A., Ruble, S. B., Valic, Z., Buckwalter, J. B., and Clifford, P. S. (2004). Intensity and duration threshold for aerobic exercise-induced analgesia to pressure pain. Arch. Phys. Med. Rehabil. 85, 1183–1187. doi: 10.1016/j.apmr.2003.09.010
Hu, L., Xiao, P., Zhang, Z. G., Mouraux, A., and Iannetti, G. D. (2014). Single-trial time-frequency analysis of electrocortical signals: baseline correction and beyond. Neuroimage 84, 876–887. doi: 10.1016/j.neuroimage.2013.09.055
Huang, G., Xiao, P., Hung, Y. S., Iannetti, G. D., Zhang, Z. G., and Hu, L. (2013). A novel approach to predict subjective pain perception from single-trial laser-evoked potentials. Neuroimage 81, 283–293. doi: 10.1016/j.neuroimage.2013.05.017
Iannetti, G. D., Zambreanu, L., and Tracey, I. (2006). Similar nociceptive afferents mediate psychophysical and electrophysiological responses to heat stimulation of glabrous and hairy skin in humans. J. Physiol. 577, 235–248. doi: 10.1113/jphysiol.2006.115675
Janal, M. N., Colt, E. W. D., Clark, C. W., and Glusman, M. (1984). Pain sensitivity, mood and plasma endocrine levels in man following long-distance running: effects of naloxone. Pain 19, 13–25. doi: 10.1016/0304-3959(84)90061-7
Jin, Q. Q., Wu, G. Q., Peng, W. W., Xia, X. L., Hu, L., and Iannetti, G. D. (2018). Somatotopic representation of second pain in the primary somatosensory cortex of humans and rodents. J. Neurosci. 38, 5538–5550. doi: 10.1523/jneurosci.3654-17.2018
Jobson, S. A., Hopker, J., Arkesteijn, M., and Passfield, L. (2013). Inter- and intra-session reliability of muscle activity patterns during cycling. J. Electromyogr. Kinesiol. 23, 230–237. doi: 10.1016/j.jelekin.2012.08.013
Jones, M. D., Nuzzo, J. L., Taylor, J. L., and Barry, B. K. (2019). Aerobic exercise reduces pressure more than heat pain sensitivity in healthy adults. Pain Med. 20, 1534–1546. doi: 10.1093/pm/pny289
Jones, M. D., Taylor, J. L., and Barry, B. K. (2017). Occlusion of blood flow attenuates exercise-induced hypoalgesia in the occluded limb of healthy adults. J. Appl. Physiol. 122, 1284–1291. doi: 10.1152/japplphysiol.01004.2016
Jones, M. D., Taylor, J. L., Booth, J., and Barry, B. K. (2016). Exploring the mechanisms of exercise-induced hypoalgesia using somatosensory and laser evoked potentials. Front. Physiol. 7:581. doi: 10.3389/fphys.2016.00581
Jongerius, N., Wainwright, B., Wheat, J., and Bissas, A. (2021). Prevalence and functional implications of soleus and tibialis anterior activation strategies during cycling. J. Sports Sci 17, 1–8. doi: 10.1080/02640414.2021.1939981
Kemppainen, P., Paalasmaa, P., Pertovaara, A., Alila, A., and Johansson, G. (1990). Dexamethasone attenuates exercise-induced dental analgesia in man. Brain Res. 519, 329–332. doi: 10.1016/0006-8993(90)90096-t
Klimesch, W. (2012). α-band oscillations, attention, and controlled access to stored information. Trends Cogn. Sci. 16, 606–617. doi: 10.1016/j.tics.2012.10.007
Kodesh, E., and Weissman-Fogel, I. (2014). Exercise-induced hypoalgesia–interval versus continuous mode. Appl. Physiol. Nutr. Metab. 39, 829–834. doi: 10.1139/apnm-2013-0481
Koltyn, K. F. (2000). Analgesia following exercise: a review. Sports Med. 29, 85–98. doi: 10.2165/00007256-200029020-00002
Laufs, H., Krakow, K., Sterzer, P., Eger, E., Beyerle, A., Salek-Haddadi, A., et al. (2003). Electroencephalographic signatures of attentional and cognitive default modes in spontaneous brain activity fluctuations at rest. Proc. Natl. Acad. Sci. U.S.A. 100, 11053–11058. doi: 10.1073/pnas.1831638100
Lee, P. H., Macfarlane, D. J., Lam, T. H., and Stewart, S. M. (2011). Validity of the international physical activity questionnaire short form (IPAQ-SF): a systematic review. Int. J. Behav. Nutr. Phys. Act. 8:115. doi: 10.1186/1479-5868-8-115
Lei, X., and Liao, K. (2017). Understanding the influences of EEG reference: a large-scale brain network perspective. Front. Neurosci. 11:205. doi: 10.3389/fnins.2017.00205
Li, X., Lo, W. L. A., Lu, S. W., Liu, H., Lin, K. Y., Lai, J. Y., et al. (2020). Trunk muscle activity during pressure feedback monitoring among individuals with and without chronic low Back pain. BMC Musculoskelet. Disord. 21:569. doi: 10.1186/s12891-020-03565-y
Lo, W. L. A., Liang, Z., Li, W., Luo, S., Zou, Z., Chen, S., et al. (2019). The Effect of judo training on set-shifting in school children. Biomed. Res. Int. 2019:2572016. doi: 10.1155/2019/2572016
Loew, L., Brosseau, L., Wells, G. A., Tugwell, P., Kenny, G. P., Reid, R., et al. (2012). Ottawa panel evidence-based clinical practice guidelines for aerobic walking programs in the management of osteoarthritis. Arch. Phys. Med. Rehabil. 93, 1269–1285. doi: 10.1016/j.apmr.2012.01.024
Mathewson, K. E., Lleras, A., Beck, D. M., Fabiani, M., Ro, T., and Gratton, G. (2011). Pulsed out of awareness: EEG alpha oscillations represent a pulsed-inhibition of ongoing cortical processing. Front. Psychol. 2:99. doi: 10.3389/fpsyg.2011.00099
Meeus, M., Roussel, N. A., Truijen, S., and Nijs, J. (2010). Reduced pressure pain thresholds in response to exercise in chronic fatigue syndrome but not in chronic low back pain: an experimental study. J. Rehabil. Med. 42, 884–890. doi: 10.2340/16501977-0595
Meints, S. M., Mawla, I., Napadow, V., Kong, J., Gerber, J., Chan, S. T., et al. (2019). The relationship between catastrophizing and altered pain sensitivity in patients with chronic low-back pain. Pain 160, 833–843. doi: 10.1097/j.pain.0000000000001461
Meng, J., Jackson, T., Chen, H., Hu, L., Yang, Z., Su, Y., et al. (2013). Pain perception in the self and observation of others: an ERP investigation. Neuroimage 72, 164–173. doi: 10.1016/j.neuroimage.2013.01.024
Micalos, P. S., Harris, J., Drinkwater, E. J., Cannon, J., and Marino, F. E. (2015). Perceptual and cerebro-spinal responses to graded innocuous and noxious stimuli following aerobic exercise. J. Sports Med. Phys. Fitness 55, 1407–1415.
Micalos, P. S., Korgaonkar, M. S., Drinkwater, E. J., Cannon, J., and Marino, F. E. (2014). Cerebral responses to innocuous somatic pressure stimulation following aerobic exercise rehabilitation in chronic pain patients: a functional magnetic resonance imaging study. Int. J. Gen. Med. 7, 425–432. doi: 10.2147/ijgm.S55169
Nahman-Averbuch, H., Leon, E., Hunter, B. M., Ding, L., Hershey, A. D., Powers, S. W., et al. (2019). Increased pain sensitivity but normal pain modulation in adolescents with migraine. Pain 160, 1019–1028. doi: 10.1097/j.pain.0000000000001477
Naugle, K. M., and Riley, J. L. III (2014). Self-reported physical activity predicts pain inhibitory and facilitatory function. Med. Sci. Sports Exerc. 46, 622–629. doi: 10.1249/MSS.0b013e3182a69cf1
Naugle, K. M., Fillingim, R. B., and Riley, J. L. III (2012). A meta-analytic review of the hypoalgesic effects of exercise. J. Pain 13, 1139–1150. doi: 10.1016/j.jpain.2012.09.006
Naugle, K. M., Naugle, K. E., Fillingim, R. B., Samuels, B., and Riley, J. L. III (2014). Intensity thresholds for aerobic exercise-induced hypoalgesia. Med. Sci. Sports Exerc. 46, 817–825. doi: 10.1249/mss.0000000000000143
Pertovaara, A., Kauppila, T., and Hämäläinen, M. M. (1996). Influence of skin temperature on heat pain threshold in humans. Exp. Brain Res. 107, 497–503. doi: 10.1007/bf00230429
Peterson, J. A., Schubert, D. J., Campbell, J., Bemben, M. G., and Black, C. D. (2019). Endogenous pain inhibitory function: endurance-trained athletes vs active controls. Pain Med. 20, 1822–1830. doi: 10.1093/pm/pnz014
Polaski, A. M., Phelps, A. L., Kostek, M. C., Szucs, K. A., and Kolber, B. J. (2019). Exercise-induced hypoalgesia: a meta-analysis of exercise dosing for the treatment of chronic pain. PLoS One 14:e0210418. doi: 10.1371/journal.pone.0210418
Rainville, P., Duncan, G. H., Price, D. D., Carrier, B., and Bushnell, M. C. (1997). Pain affect encoded in human anterior cingulate but not somatosensory cortex. Science 277, 968–971. doi: 10.1126/science.277.5328.968
Scheef, L., Jankowski, J., Daamen, M., Weyer, G., Klingenberg, M., Renner, J., et al. (2012). An fMRI study on the acute effects of exercise on pain processing in trained athletes. Pain 153, 1702–1714. doi: 10.1016/j.pain.2012.05.008
Scherrer, G., Imamachi, N., Cao, Y. Q., Contet, C., Mennicken, F., O’Donnell, D., et al. (2009). Dissociation of the opioid receptor mechanisms that control mechanical and heat pain. Cell 137, 1148–1159. doi: 10.1016/j.cell.2009.04.019
Staud, R., Robinson, M. E., and Price, D. D. (2005). Isometric exercise has opposite effects on central pain mechanisms in fibromyalgia patients compared to normal controls. Pain 118, 176–184. doi: 10.1016/j.pain.2005.08.007
Steele, J. R., Coltman, C. E., and McGhee, D. E. (2020). Effects of obesity on breast size, thoracic spine structure and function, upper torso musculoskeletal pain and physical activity in women. J. Sport Health Sci. 9, 140–148. doi: 10.1016/j.jshs.2019.05.003
Thorén, P., Floras, J. S., Hoffmann, P., and Seals, D. R. (1990). Endorphins and exercise: physiological mechanisms and clinical implications. Med. Sci. Sports Exerc. 22, 417–428.
Tu, Y., Zhang, Z., Tan, A., Peng, W., Hung, Y. S., Moayedi, M., et al. (2016). Alpha and gamma oscillation amplitudes synergistically predict the perception of forthcoming nociceptive stimuli. Hum. Brain Mapp. 37, 501–514. doi: 10.1002/hbm.23048
Vaegter, H. B., Handberg, G., Emmeluth, C., and Graven-Nielsen, T. (2017a). Preoperative hypoalgesia after cold pressor test and aerobic exercise is associated with pain relief 6 months after total knee replacement. Clin. J. Pain 33, 475–484. doi: 10.1097/ajp.0000000000000428
Vaegter, H. B., Hoeger Bement, M., Madsen, A. B., Fridriksson, J., Dasa, M., and Graven-Nielsen, T. (2017b). Exercise increases pressure pain tolerance but not pressure and heat pain thresholds in healthy young men. Eur. J. Pain 21, 73–81. doi: 10.1002/ejp.901
Valentini, E., Hu, L., Chakrabarti, B., Hu, Y., Aglioti, S. M., and Iannetti, G. D. (2012). The primary somatosensory cortex largely contributes to the early part of the cortical response elicited by nociceptive stimuli. Neuroimage 59, 1571–1581. doi: 10.1016/j.neuroimage.2011.08.069
van Middelkoop, M., Rubinstein, S. M., Verhagen, A. P., Ostelo, R. W., Koes, B. W., and van Tulder, M. W. (2010). Exercise therapy for chronic nonspecific low-back pain. Best Pract. Res. Clin. Rheumatol. 24, 193–204. doi: 10.1016/j.berh.2010.01.002
van Poppel, M. N., Chinapaw, M. J., Mokkink, L. B., van Mechelen, W., and Terwee, C. B. (2010). Physical activity questionnaires for adults: a systematic review of measurement properties. Sports Med. 40, 565–600. doi: 10.2165/11531930-000000000-00000
Vierck, C. J. Jr., Staud, R., Price, D. D., Cannon, R. L., Mauderli, A. P., and Martin, A. D. (2001). The effect of maximal exercise on temporal summation of second pain (windup) in patients with fibromyalgia syndrome. J. Pain 2, 334–344. doi: 10.1054/jpai.2001.25533
Weisz, N., Wühle, A., Monittola, G., Demarchi, G., Frey, J., Popov, T., et al. (2014). Prestimulus oscillatory power and connectivity patterns predispose conscious somatosensory perception. Proc. Natl. Acad. Sci. U.S.A. 111, E417–E425. doi: 10.1073/pnas.1317267111
Williams, N. (2017). The borg rating of perceived exertion (RPE) scale. Occup. Med. 67, 404–405. doi: 10.1093/occmed/kqx063
Keywords: aerobic exercise, hypoalgesia, event-related potential, oscillation, exercise intensity
Citation: Zheng KY, Chen CC, Yang SY and Wang XQ (2021) Aerobic Exercise Attenuates Pain Sensitivity: An Event-Related Potential Study. Front. Neurosci. 15:735470. doi: 10.3389/fnins.2021.735470
Received: 02 July 2021; Accepted: 27 August 2021;
Published: 21 September 2021.
Edited by:
Julian Keil, University of Kiel, GermanyReviewed by:
Daya Shankar Gupta, Husson University, United StatesZhenyong Lyu, Yangzhou University, China
Stacy Christine Stolzman, Concordia University (Wisconsin), United States
Copyright © 2021 Zheng, Chen, Yang and Wang. 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: Suyong Yang, ysy77@163.com; Xueqiang Wang, wangxueqiang@sus.edu.cn
†These authors have contributed equally to this work