Abstract

Objectives. This study is aimed at exploring the effects of virtual reality (VR) training on postural control, measured by anticipatory and compensatory postural adjustments (APAs and CPAs, respectively), in patients with chronic nonspecific low back pain (CNLBP) and the potential neuromuscular mechanism of VR training. Methods. Thirty-four patients were recruited and randomly assigned to the VR group (), the motor control exercise group (MCE, ) and the control group (CG, ). The VR group received VR training using Kinect Xbox 360 systems and magnetic therapy. Besides magnetic therapy, the participants in the MCE group performed real-time ultrasound-guided abdominal drawing-in maneuver (ADIM) and four-point kneeling exercise. The CG only received magnetic therapy. Surface muscle electromyography (sEMG) was used to record the muscle activities of transverse abdominis (TrA), multifidus (MF), lateral gastrocnemius (LG), and tibialis anterior (TA) during ball-hitting tasks. The muscle activation time and integrals of the electromyography activities (IEMGs) during the APA and CPA stages were calculated and used in the data analysis. The visual analogue scale (VAS) and Oswestry dysfunction index (ODI) scores were also recorded. Results. A significant interaction effect of was observed on the activation time of TrA () and MF (). The post-intervention activation time of the TrA was earlier in the VR group (). In contrast, the post-intervention activation time of the MF was significantly delayed in the VR group (). The IEMGs of TrA () and TA () during CPA1 significantly decreased only in the VR group after the intervention. The VAS scores of three group participants showed significant decreases after intervention (). Conclusions. Patients with CNLBP showed reciprocal muscle activation patterns of the TrA and MF muscles after VR training. VR training may be a potential intervention for enhancing the APAs of the patients with CNLBP.

1. Introduction

Chronic nonspecific low back pain (CNLBP) is one of the most common musculoskeletal disorders worldwide [1]. CNLBP may lead to a poor quality of life and increase the economic burden [2]. In recent years, altered postural control patterns have been reported as one of the most important factors that may contribute to CNLBP development [3, 4]. Postural control can be attributed to two different neuromuscular control mechanisms: anticipatory postural adjustments (APAs) and compensatory postural adjustments (CPAs). APAs are feed-forward adjustments that occur before perturbation to minimize the effects of predictable perturbations [5], whereas CPAs are reflexive adjustments made to maintain equilibrium after the onset of perturbations [6]. The central nervous system (CNS) employs APAs and CPAs to maintain stability. Previous studies have reported delayed or absent trunk muscle activation in rapid arm movement or in response to sudden loading in patients with CNLBP [3, 79]. For instance, when a resistance force at the back was suddenly released, patients with chronic low back pain (CLBP) showed delayed activation time of abdominal muscles (transverse abdominis (TrA)) in the APA phase in comparison with healthy controls [9]. These findings were supported by the pain adaptation model, which postulates that long-term pain would decrease muscle activation when the muscle is activated as an agonist. Since upper limb loading produces a flexion moment, the TrA is the agonist during sudden upper limb loading and may show decreased electrical activity in the low back pain (LBP) group [9]. The impaired APAs in patients with CNLBP could activate trunk muscles in response to sudden perturbations over time and subsequently lead to the recurrence of LBP. Thus, a training program aimed at enhancing APAs is crucial for CNLBP patients.

Motor control exercise (MCE) is a common exercise rehabilitation program for patients with CNLBP [10]. However, the influence of MCE on APAs in patients with CNLBP is inconsistent. Several studies have reported that MCE could not enhance APA capacity [1113], while Tsao and Hodges found that MCE can cause early muscle activation [14]. However, the muscle activation time after MCE was later than the time window between -100 and 50 ms, suggesting that MCE may not be a suitable rehabilitation training to enhance the APAs of CNLBP patients and that it may be more beneficial in enhancing postural stability by strengthening muscle power and endurance rather than increasing feed-forward control (e.g., APAs) by the CNS [11, 13].

Many studies have shown that virtual reality- (VR-) based training could improve participants’ balance ability by providing multisensory feedback, for example, in chronic poststroke survivors [15, 16] and elderly people [17]. Two previous studies found that VR training could enhance APA capacity [18, 19]. Ida et al. found that when healthy participants lifted one foot to avoid a real or virtual obstacle, muscle activation of the supporting leg and trunk was observed during APAs in both real and VR environments, even though muscle activation in the virtual display setting was smaller than that in the real setting. In another study, PD (Parkinson’s disease) participants showed shorter movement time and higher peak velocity of arm movements and longer APAs while catching a fast ball than they did while catching a slow ball in a VR environment. These results suggest that VR-based training can improve APAs. VR training has also been gradually applied to rehabilitation training programs in patients with CNLBP, and it has been shown to reduce pain and improve dysfunction [20]. However, the effect of VR training on APAs in patients with CNLBP is still unknown.

The present study is aimed at investigating the effect of VR-based training on APAs in patients with CNLBP by employing a ball-hitting test. We hypothesized that in comparison with MCE, VR-based training could enhance the APA capacity, including the muscle activation time and integrals of the EMG activities (IEMGs), in patients with CNLBP. The findings of the present study could help verify whether VR training is an effective treatment for enhancing postural control in patients with CNLBP.

2. Methods

2.1. Participants

Thirty-four right-handed participants with CNLBP were recruited for this study. The inclusion criteria for CNLBP participants were as follows: (1) age between 18 and 40 years, (2) persistent or periodic LBP for longer than 3 months, and (3) no referred symptoms of radiating pain below the knee or paresthesia during the straight-leg raise test. The exclusion criteria were as follows: (1) history of pelvic or spinal column surgery in the past two years; (2) diagnosis of any specific lumbar pathological condition (such as lumbar tumors, vertebral fractures, lumbar spinal stenosis, lumbar spondylolisthesis, rheumatoid arthritis, or ankylosis) and/or severe or progressive scoliosis; (3) body mass index ; (4) history of a treatment program within the past three months; (5) pregnancy; (6) history of severe dysfunction of vital organs (heart, lungs, and kidneys) and/or cognitive deficits; and (7) history of visual or hearing problems. The participants could withdraw the experiment at any time if they (1) were unwilling to participate in this experiment, (2) felt any aggravation of pain during the treatment, (3) had any other disease due to the treatment, and (4) could not complete the proposed treatment plan. Ethical approval for this study (ethics: no. [2020]476) was obtained from the First Affiliated Hospital of Sun Yat-sen University. Written informed consent was obtained from all the participants prior to the experiment.

In our pilot study, the effect sizes () for the within factor (time) and within-between interaction () of muscle activation time were achieved in TrA muscle (0.113 and 0.244, respectively). GPower (v 3.1.9.7, Germany) was employed to calculate the sample size. 30 participants were required to achieve the statistical power in TrA.

2.2. Apparatus and Data Preprocessing
2.2.1. The Ball-Hitting Test

The participants stood at the center of a platform with their feet shoulder-width apart. They were asked to keep their elbows bent at 90° while holding a metal tray in their hands. A pressure sensor at the outer center of the base of the tray was used to determine the time point (T0) when the object landed on the tray. After a sound, a load weighing 1.5 kg was suddenly released by the experimenter from the participants’ eye level above the tray (Figure 1). The participants were asked to try their best to maintain their body stability throughout the experiment. Each participant repeated five trials with a rest period of approximately 30 s between trials. Surface electromyography (sEMG) data were simultaneously recorded during the ball-hitting test. LabView 15.0.1 software (National Instruments, Austin, TX, USA) was used to simultaneously collect the data from the pressure sensor and the sEMG system at a frequency of 1000 Hz. Sufficient test-retest reliability for five trials in the ball-hitting test was observed in all the muscle activation time and IEMGs (intraclass correlation coefficient: 0.430~0.856).

2.2.2. Surface Electromyography (sEMG)

An sEMG system (Myomonitor IV; Delsys, USA) with eight channels was used to record muscle activity. The interelectrode distance of each channel was 10 mm. Before electrode attachment, the skin of the target areas was prepared by shaving, scrubbing with fine sandpaper, and rubbing with 75% alcohol to reduce impedance. Two electrodes were placed in a vertical arrangement along the muscle fibers of the transverse abdominis (TrA, 2 cm medial and inferior to the anterior superior iliac spine), multifidus (MF, at the level of the L5 spinous process on a line from the posterior superior iliac spine to the interspace between L1 and L2), lateral gastrocnemius (LG, at 1/3rd of the line from fibular head to the lateral side of the Achilles tendon insertion), and tibialis anterior (TA, at 1/3rd on the line between the tip of the fibula and the tip of the medial malleolus). The positions of the electrodes for these four muscles were based on previous studies [7, 8]. The reference electrode was placed on the patella on the dominant side. The sEMG data were sampled at a rate of 1000 Hz.

2.2.3. sEMG Data Preprocessing

The sEMG signals were processed using MATLAB software (The MathWorks Inc., Natick, MA, USA). The raw sEMG signals were rectified and bandpass-filtered (30–400 Hz). The first instance when the pressure signal acquired from the pressure sensor was equal to or greater than 5% of the peak magnitude for at least 20 ms continuously was considered as time zero (, the onset of a rise in the signal of the pressure sensor). The sEMG signals in the ball-hitting test were aligned using T0. The muscle activation time for each trial was detected in a time window from -300 ms to +200 ms in relation to T0. The previously described common methods for detecting muscle activation were not suitable in the present study because the low signal-to-noise ratio (SNR) of the sEMG signals increased the difficulty in detecting the muscle activation time. When the SNR of sEMG signals is very low, the onset time cannot be easily determined by visual inspection of the signal or by simply setting an amplitude threshold [21]. The TKE operation, which simultaneously considers the amplitude and instantaneous frequency of the surface EMG, has been previously used for detection of muscle activation time [2123]. Thus, in the present study, the TKE operation was applied to determine the onset time of muscle activity. IEMGs were calculated in four epochs, each of 150 ms duration in relation to T0 [5]. The time windows for these four epochs were as follows: (1) from -250 ms to -100 ms (anticipatory reaction, APA1); (2) -100 ms to +50 ms (anticipatory reaction, APA2); (3) +50 ms to +200 ms (early compensatory adjustment, CPA1); and (4) +200 ms to +350 ms (late compensatory adjustment, CPA2). All IEMGs were corrected by the baseline IEMGs from -600 ms to -450 ms relative to T0. The mean muscle activation times and IEMGs were used in the following data analysis.

2.3. Interventions

All participants were randomly assigned to three groups: the VR, MCE, and control groups. Participants in the control group received conventional thermal magnetic therapy only, which was performed for 20 min with a medium heat level per day. In addition to thermal magnetic therapy, the participants in the VR group received VR training, whereas those in the MCE group performed MCE. All interventions were performed for two weeks, five days per week. More details regarding VR training and MCE are described below. Three physiotherapists responsible for intervention training were blinded to participants’ outcome measures. Each of them took charge of the training program of one group. Two investigators were responsible for demographic information collection and outcome measures, respectively. These two investigators were blinded to the intervention allocation of each participant.

2.3.1. VR Training

The “Fruit Ninja” game displayed by the Kinect Xbox 360 system was employed for VR training. During the entire training session, the participant stood with feet shoulder-width apart at a distance of 1.5 m in front of the screen. The participants were asked to crush the fruit by waving their hands as much as possible, while simultaneously trying their best to avoid the “bombs” in the game (Figure 2). During training, the participants were asked to avoid trunk bending or turning. The participants needed to complete six sessions of VR training per day. Each session lasted three minutes, with a break of 2 min between sessions. It took approximately half an hour to complete the VR training per day.

2.3.2. Motor Control Exercise (MCE)

The abdominal drawing-in maneuver (ADIM), which is a key technique in MCE training, was designed to enhance coactivation of the TrA and MF to stabilize the trunk before body movement. The ADIM in the present study was based on that used in a previous study [24]. Before the training, the participant learned how to specifically activate the TrA muscle under the guidance of real-time ultrasound without obvious contraction of the internal oblique and external oblique muscles simultaneously. When the participant could perform the ADIM appropriately, MCE training began. The first step of MCE involved three sets of ultrasound-guided ADIM training per day, with a short break of approximately 2 min between sets. Each set involved 10 repetitions of the ultrasound-guided ADIM, each of which lasted for 10 s. In the second step, the participant completed a four-point kneeling exercise. In the second step, the participant completed a four-point kneeling exercise. In the first stage, the participant was instructed to lift one arm with the elbow and wrist extended for 5 s and maintain the TrA contraction at the same time. Each side of the upper limb was repeated thrice with a break of 15 s. In the second stage, the participant lifted one leg with the hip and knee extended for 5 s in a four-point kneeling position. Each side of the leg was repeated thrice with a break of 15 s. In the last stage, the participants raised one arm and the contralateral leg to a horizontal position (bird dog) and held it for 5 s in a four-point kneeling position. Each movement was repeated three times with a break of 15 s. In the last stage, the participants raised one arm and the contralateral leg to a horizontal position (bird dog) and held it for 5 s in a four-point kneeling position. Each movement was repeated three times with a break of 15 s. Each participant was instructed by an experienced physiotherapist during the MCE training. The MCE training required approximately 30 min per day after conventional thermal magnetic therapy.

2.4. Procedure

Before the training, the participants completed a demographic information questionnaire assessing gender, height, weight, age, and medical history. A visual analogue scale [25] (VAS) (anchored with “painless” at 0 and “intolerable pain” at 10) was used to measure pain intensity, and the Oswestry disability index (ODI) [26] was used to assess function disability pre- and posttraining. The ball-hitting test with sEMG recording was also conducted pre- and postintervention. VR-based training and MCE training were performed after conventional thermal magnetic therapy. The experiment flow is shown in Figure 3.

2.5. Statistical Analysis

An independent -test was used to compare between-group differences in all demographic variables except sex (Table 1). The chi-squared test was used to compare between-group differences in sex. The muscle activation times and IEMGs of the four target muscles as well as data from the clinical assessments (including VAS and ODI scores) were analyzed using the two-way mixed-design repeated-measure analysis of variance (ANOVA) with a between-subject factor of group (CG, MCE, and VR groups) and a within-subject factor of time (pre- and posttraining). Post hoc pairwise comparisons were applied when a significant effect was observed. The significance level was set at . SPSS software (version 23.0; IBM, Armonk, NY, USA) was used for data analysis.

3. Results

3.1. Demographics and Clinical Assessments

The demographic information of all participants is shown in Table 1. No between-group differences were found in sex, age, weight, height, BMI, or pain duration () (Table 1).

3.2. Muscle Activation Time

The activation times of TrA, MF, LG, and TA are presented in Figure 4. A significant main effect of time was observed on the MF activation time (, , ). A significant interaction effect of was observed on the activation time of TrA (, , ) and MF (, , ). Post hoc analysis for the significant interaction effect revealed that the activation time of TrA after VR training was significantly earlier than that before training (). However, the activation time of the MF muscle in the VR group after training was significantly delayed () in comparison with that before training. No changes in the activation times of TrA and MF were observed in the MCE (TrA: ; MF: ) and control groups (TrA: ; MF: ). Other main effects of time and group and the interaction effect of were not significant ().

3.3. IEMGs of the Four Muscles

The IEMGs in the APA1 and APA2 phases are presented in Figures 5 and 6, respectively. The ratios and values for the mixed model of the four muscles during APA1 and APA2 are shown in Table 2. A significant effect of time was observed on the IEMGs of MF (, , ) and TA (, , ) during APA1 and on the IEMGs of MF (, , ) and TA (, , ) during APA2. We also observed a significant effect of group on the IEMGs of LG (, , ) during APA1. A significant interaction effect was only observed on the TA during APA2 (, , ). Other main effects of time and group or the interaction effect of were not significant (). Post hoc analysis for the significant interaction effect showed that the IEMGs of TA during APA2 were significantly decreased after the intervention only in the VR group (), which could not be observed in the other two groups.

The IEMGs in the CPA1 and CPA2 phases are presented in Figures 7 and 8, respectively. The ratios and values for the mixed model of the four muscles during CPA1 and CPA2 are shown in Table 3. A significant effect of time was observed on the IEMGs of MF (, , ), LG (, , ), and TA (, , ) during CPA1. A interaction effect was observed on the IEMGs of TrA (, , ) and TA (, , ) during CPA1. Other main effects of time and group and the interaction effect of were not significant (). Post hoc analysis for the significant interaction effect showed that the IEMGs of TrA () and TA () during CPA1 significantly decreased after the intervention only in the VR group, which could not be observed in the other two groups.

3.4. Pain-Related Clinical Outcomes

Pain-related clinical outcomes, including VAS and ODI scores, are shown in Table 4. A significant main effect of time was observed on the VAS score (, , ). The main effect of group (, , ) and the interaction effect (, , ) were not significant for the VAS score. Time (, , ), group (, , ), and (, , ) did not show significant effects on ODI scores.

4. Discussion

The present study investigated the effect of VR training on postural control in patients with CNLBP through measurement of APAs and CPAs in an external postural perturbation task. A novel finding was that the TrA muscle showed earlier muscle activation in APAs after VR training, whereas the MF muscle showed delayed activation after VR training. These findings were not observed in the other two groups. These findings suggest that VR-based training is likely to improve the APAs of CNLBP patients.

Postural control requires the central nervous system (CNS) to receive and integrate multisensory inputs (including vision, vestibular sense, and proprioceptive and tactile information) and thereby coordinate and control the postural muscles to maintain balance and stability of the body [27]. Impaired postural control has been reported to involve APAs and CPAs of trunk muscles [28, 29]. Most previous studies reported that delayed activation of abdominal muscles, especially TrA, in patients with CNLBP was commonly observed in patients with LBP in the postural control assessment, which may increase the recurrence of LBP in patients with CNLBP [4, 30, 31]. The potential reason was that the delayed muscle activation time of TrA in the APAs of patients with CNLBP was associated with remodeling of the motor cortex [32]. Sadeghi et al. found that the prefrontal cortex is activated in patients with LBP during postural interference [33], suggesting that the prefrontal and motor cortices are involved in anticipatory processing.

In this study, the muscle onset time of TrA in CNLBP patients through two weeks of VR exercise training was significantly earlier than that before training, suggesting that VR training could improve APAs. These findings are supported by the results of previous studies [18, 19]. For example, Su et al. showed that training for a ball-catching task in a VR environment could improve PD patients’ ability to perform APAs. The possible mechanism is likely that VR training enhanced the activation of the frontoparietal and sensorimotor networks by providing visual cues and visual feedback to patients [34, 35]. The information provided by the virtual environment can increase the activation of the frontal lobe, which participates in the anticipatory process [34]. The anticipation process is a top-down cognitive process [33, 34], the preparatory state of which could be comparable to that during APAs [35]. Thus, in the present study, VR training could improve the APA capacity of CNLBP patients because the visual cues in the dynamic environment of “Fruit Ninja” elicited the participants’ anticipation of the movement in response to the objects in the virtual environment. In addition, the “Fruit Ninja” game required the participant to move the arm rapidly in different directions to cut the fruits. Hodges reported that healthy participants would activate TrA muscle earlier than other muscles to maintain postural stability, when performing the rapid shoulder flexion, abduction, extension [36]. The role of TrA muscle in the arm movement in these studies was like to increase the stiffness of the lumbar spine through raising the intra-abdominal pressure and increase the tension of thoracolumbar fascia. Thus, the participants were required to activate TrA earlier to maintain postural stability during the arm movement of the “Fruit Ninja” game.

A previous study reported opposite muscle activation patterns of TrA (delayed activation) and MF (earlier activation) in patients with CNLBP [36]. This reciprocal activation pattern in the TrA-MF pair may be related to an efficient strategy for postural control by the CNS when the perturbation is predictable [37]. The present study showed that the activation of MF was delayed after two weeks of VR training. These results suggest that after VR training, the activation pattern of the MF muscle appears to be similar to that of the control participants.

This study also found that muscle activity of the TA during APA2 was weakened after VR training. A previous study reported that patients with CLBP adopted a body-and-trunk–stiffening strategy and relied more on ankle proprioception to control their posture while standing due to the weakness of their trunk muscles [38]. In this study, the IEMGs of the TA during APA2 decreased after VR training, possibly due to the improvement in TrA muscle activation after VR training. Thus, the improvement in the TrA muscle may enable CLBP patients to achieve better coordination of the deep trunk muscle to maintain stability and rely less on the ankle strategy.

The present study also found that the muscle activities of the TrA and TA in the CPA1 stage were weakened after VR training. These results suggest that patients with CLBP tend to use the APA strategy to reduce the compensatory response of the muscle to external interference after VR training. These findings are consistent with those reported by Santos et al. [5] and Liang et al. [37]. Santos et al. found that when perturbations were predictable, stronger APAs were significantly related to smaller compensatory activities of muscles and COP displacements in response to external perturbation. The findings in Liang et al.’s study showed that after auditory training, the participants demonstrated stronger APAs and less demands on CPAs. In the present study, the IEMGs of TrA and TA of CNLBP participants decreased in CPA1 after VR training, which may result from the APA improvement due to the prediction elicited by the visual information in the virtual environment.

The muscle activation times and IEMGs showed no changes pre- and posttraining in the MCE and control groups. These results were supported by the findings reported by Vasseljen et al. [11] and Lomond et al. [12]. Vasseljen et al. found no significant difference in the activation time of abdominal muscles in patients with CLBP between pretraining and after 8 weeks of MCE training. A planned secondary analysis conducted by Lomond et al. revealed that low back stabilization or movement system impairment treatments did not ameliorate the CLBP participants’ APAs impairment, as reflected by the lack of significant effects of treatment on IEMGs. This may be because MCE can enhance postural stability mainly by strengthening muscle power and endurance [12, 13]. However, APAs rely more on brain activities, including motor planning and the excitability of the motor cortex, to predict external perturbations [35, 39]. Thus, in the present study, MCE and conventional physiotherapy did not improve muscle activation time in patients with CNLBP, because the activities in these two groups did not require the participants to perform accurate predictions during the intervention.

As for the clinical outcomes, the VAS findings suggested that LBP decreased after training in the three groups of participants, but there were no statistical differences in VAS scores among the three groups. The improvement in dysfunction revealed by ODI scores was also not significantly different after training among the three groups. Virtual walking was reported to help reduce pain and kinesiophobia and improve function in the short term in patients with chronic nonspecific LBP [40]. The potential mechanism underlying these findings is that the immediate multisensory feedback provided by VR training can improve pain processing in the CNS [41, 42]. The findings in the MCE group are supported by the results of previous studies [36, 43, 44], which showed that MCE could reduce pain and improve dysfunction in CNLBP patients [45], but the clinical improvements did not statistically differ from other treatments [36, 43, 44]. The potential reason for no between-group differences in the pain-related clinical outcomes of the present study was short-term intervention, which may be not long enough to elicit differential treatment effect. The duration of MCE intervention period of 12 weeks was reported to induce superior effect than general exercise program [46]. However, for the duration of the MCE program of 6 weeks, no significant difference was reported between MCE and generally exercise program [47]. Thus, no firm conclusion could be drawn on the comparison of the effectiveness of each intervention program on pain-related clinical outcomes. VR-based intervention may potentially be a beneficial adjunct to MCE intervention for LBP rehabilitation. Further investigation on VR-based training is at least warranted.

4.1. Limitations

There are several limitations to the present study. First, sample size was calculated only based on the core muscle of TrA rather than other muscles, which may reduce the statistical power of other muscles. The present study was a preliminary study to explore the effect of VR training on postural control. Large trial is required in the future study. Second, the age range of the sample population was only 19–30 years. Thus, the results may not be directly generalizable to participants beyond this age range. Third, the intervention period was only two weeks, which was potentially not long enough to elicit the treatment effect. Future studies should lengthen the training period to six or more weeks to further confirm the effect of VR training on APAs and pain-related clinical outcomes. In addition to assessments of behavioral outcomes, brain imaging and electrophysiological techniques could be employed to investigate the underlying neural mechanisms in future research.

4.2. Conclusions

VR-based training may be an alternative to MCE to enhance APAs by altering the muscle activation pattern of the trunk and lower limb muscles in response to perturbation. The results of this study provide a new potential treatment for APA impairment in CLBP. However, the effect of VR training on the clinical pain symptoms requires further work to verify. In addition to sEMG, future studies can use electrophysiological and brain imaging methods to investigate the underlying neural mechanisms for the effects of VR training on postural control.

Data Availability

The EMG datasets analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declared no conflict of interest.

Authors’ Contributions

Zhicheng Li and Qiuhua Yu contributed equally to this work.

Acknowledgments

We thank Mr. Ze Yuan for his help in the patients’ recruitment and data collection. In addition, Zhicheng Li gives sincere thanks to Baoyan Wang for her encouragement, patience, and care during the experiment period. This research project was supported by the National Natural Science Foundation of China (grant numbers 82002375, 81772434, and 32071316); the Guangdong Province Medical Science Technology Research Grant (grant number A2019452) and the Guangdong Basic and Applied Basic Research Foundation (No. 2020A1515011356); the Guangzhou Science and Technology Program key projects (grant numbers 201704020122 and 201907010034); and the Non-Profit Central Research Institute Fund of Chinese Academy of Medical Sciences (No. 2020-JKCS-005).