Skip to main content

STUDY PROTOCOL article

Front. Neurol., 23 August 2022
Sec. Neurorehabilitation

Immediate and short-term effects of continuous theta burst transcranial magnetic stimulation over contralesional premotor area on post-stroke spasticity in patients with severe hemiplegia: Study protocol for a randomized controlled trial

\nXiupan Wei,&#x;Xiupan Wei1,2Nan Xia,&#x;Nan Xia1,2Yang-An Li,Yang-An Li1,2Minghui Gu,Minghui Gu1,2Tongming Zhang,Tongming Zhang1,2Wei Gao
Wei Gao3*Yali Liu,
Yali Liu1,2*
  • 1Department of Rehabilitation Medicine, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
  • 2World Health Organization Collaborating Centre for Training and Research in Rehabilitation, Wuhan, China
  • 3Department of Traumatic Surgery, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Background: Post-stroke spasticity is an important complication that greatly affects survivors' functional prognosis and daily activities. Increasing evidence points to aberrant contralesional neuromodulation compensation after brain injury as a possible culprit for increased spasticity in patients with severe stroke. Hyperactivity of the contralesional premotor area (cPMA) was supposed to be highly correlated with this progression. This study aims to demonstrate the immediate and short-term efficacy of continuous theta-burst stimulation (cTBS) targeting cPMA on upper limb spasticity in severe subacute stroke patients.

Methods: This trial is a single-center, prospective, three-group randomized controlled trial. Forty-five eligible patients will be recruited and randomized into three groups: the sham-cTBS group (sham cTBS targeting contralesional PMA), the cTBS-cM1 group (cTBS targeting contralesional M1), and the cTBS-cPMA group (cTBS targeting contralesional PMA). All subjects will undergo comprehensive rehabilitation and the corresponding cTBS interventions once a day, five times a week for 4 weeks. Clinical scales, neurophysiological examinations, and neuroimaging will be used as evaluation tools in this study. As the primary outcome, clinical performance on muscle spasticity of elbow/wrist flexor/extensors and upper-limb motor function will be evaluated with the modified Ashworth scale and the Fugl-Meyer Assessment of Upper Extremity Scale, respectively. These scale scores will be collected at baseline, after 4 weeks of treatment, and at follow-up. The secondary outcomes were neurophysiological examinations and Neuroimaging. In neurophysiological examinations, motor evoked potentials, startle reflex, and H reflexes will be used to assess the excitability of the subject's motor cortex, reticulospinal pathway, and spinal motor neurons, respectively. Results of them will be recorded before and after the first cTBS treatment, at post-intervention (at 4 weeks), and at follow-up (at 8 weeks). Neuroimaging tests with diffusion tensor imaging for all participants will be evaluated at baseline and after the 4-week treatment.

Discussion: Based on the latest research progress on post-stroke spasticity, we innovatively propose a new neuromodulation target for improving post-stroke spasticity via cTBS. We expected that cTBS targeting cPMA would have significant immediate and short-term effects on spasticity and related neural pathways. The effect of cTBS-cPMA may be better than that of cTBS via conventional cM1. The results of our study will provide robust support for the application of cTBS neuromodulation in post-stroke spasticity after a severe stroke.

Clinical trial registration: This trial was registered with chictr.org.cn on June 13, 2022 (protocol version). http://www.chictr.org.cn/showproj.aspx?proj=171759.

Introduction

Post-stroke spasticity is defined as a velocity-dependent increase in exaggerated stretch reflex responses due to the hyperexcitable descending excitatory brainstem pathways and results in the hyperexcitability of alpha motor neurons in the spinal cord (1). Spasticity is a common complication after stroke, with the prevalence ranging from 30% to 80% of stroke survivors (2). The appearance of spasticity after stroke can significantly affect the patient's recovery of motor function and cause long-term problems such as joint contractures, abnormal pain, and deformity (3). The spasticity of most patients with mild to moderate stroke will gradually improve during the movement recovery process, but that of some patients with severe injury will progress in the opposite way (4). Reorganization of the descending cortical spinal tract (CST) and reticular spinal tract (RST) from both hemispheres may be a plausible explanation for this distinct development (5).

Both CST and RST are the two important descending pathways of postural and motor control in primates. RST is usually considered to control proximal and axial muscles, and to be involved mainly in movement and posture adjustment, while CST is thought to be involved in fine motor control (6). The two tracts are competitive but in a wonderful state of balance. Significant brain lesion after stroke causes structural damage to these descending pathways, especially in the compact pyramidal cell populations such as CST (7). However, the RST dominated by ipsilateral innervation has relatively small damage, and its participation in motor control is relatively increased and manifested in the patient's movement, such as the clumsy and synergic movement of the upper limbs (8, 9). Spasticity after a stroke occurs and diminishes in synchrony with the appearance and disappearance of these pathological movement patterns (10, 11). Although it is still not clear how the regulation of muscle tone by the two major systems is carried out, it is speculated that the loss of inhibitory regulation caused by the imbalance of CST and RST triggers the abnormal hyperactivation of the medial reticulospinal tract (mRST) (9, 1113). This imbalance is especially pronounced in patients with severe brain lesions. Recent imaging evidence suggests that ipsilateral PMA is the advanced cortex of RST in motor control, which plays a key role in ipsilateral movement (14, 15). After cortical damage, the contralesional PMA area is hyperactivated and may trigger a weakened inhibition of the ipsilateral mRST pathway (12). Due to the persistence of severe injury, the injured CST is difficult to recover, and this persistent weakening of inhibition can be inappropriately reinforced during the neural reorganization of the contralesional RST and related cortices (1). With appropriate treatment, attenuating the hyperexcitability of this pathway via contralesional PMA may be a therapeutic idea to stop spasticity from exacerbating.

Limited interventions were proven effective currently in improving spasticity. The only botulinum toxin injection backed by a lot of evidence is expensive and may require multiple iterations (16). Some central modulating strategies have also been attempted to treat post-stroke spasticity. Transcranial magnetic stimulation (TMS) is one of the options. However, ambiguity in research results leads to the disapproval of this technology on post-stroke spasticity (17, 18). The reason for this may be that previous study designs did not take into account the dynamic reorganization of CST/RST, and mainly focused on the promotion of motor recovery (1). Almost all studies with positive results were interventions targeting the bilateral M1 areas (part of CST) (1921), improvement in spasticity accompanied by significant motor recovery. This can be interpreted as a weakening of RST regulation by enhancing CST control, which may not the most direct way to regulate spasticity-related pathways. Moreover, the traditional inhibitory TMS therapy developed based on the theory of hemispheric inhibition (22) was insufficient to reverse cortical excitability in all subjects or produced short-lived or variable results (23). Continuous theta burst magnetic stimulation (cTBS) is a novel patterned protocol of TMS, it can induce long-term depression (LTD) effect on cortical excitability and lead to long-term therapeutic benefits in facilitating induction of neuronal plasticity mechanisms (24). TBS has the advantages of shorter stimulation duration, low stimulation pulse intensity, longer post-acting time, and probably improved efficiency compared with ordinary rTMS, which is regarded as a very promising therapy (25). In addition, how to assess post-stroke spasticity objectively also set a barrier for related research. Most previous studies of TBS on post-stroke spasticity mainly use the modified Ashworth scale (MAS) as their primary outcome (18), which only roughly reflects changes in spasticity, and does not benefit mechanistic exploration.

Combined with the recent mechanism research on post-stroke spasticity, we propose a novel central regulatory target for cTBS treatment of post-stroke spasticity, namely cPMA. At the same time, the contralateral M1 area, which was extensively used in previous studies, was also introduced into our study design. Based on a sham control, we will explore the effects of “direct” intervention (cTBS via cPMA) and “indirect” intervention (cTBS via cM1) on post-stroke spasticity and verify our hypothesis on the abnormal neural regulation of spasticity after stroke.

In addition to clinical assessment scales, we will introduce neurophysiological examinations and diffusion tensor imaging (DTI) for this study. Motor evoked potentials (MEP), startle reflexes (SR), and H reflex will be used to assess the excitability of the subject's motor cortex, reticulospinal pathway, and spinal motor neurons, respectively. DTI will be performed to show the contralesional neural reorganization.

We hypothesize that cTBS in cPMA ameliorates spasticity in patients with severe subacute stroke by inhibiting the overactivated mRST and reducing the excitability of spinal motor neurons, thereby improving motor outcomes of the upper extremity. This study aims to evaluate the efficacy of cTBS therapy in the contralesional PMA for severe post-stroke spasticity and confirm the above hypothesis.

Methods and analysis

Study design

The study will be a non-blinded pilot randomized controlled trial with three parallel groups. The study will be performed in accordance with the Declaration of Helsinki. This protocol follows the Consolidated Standards of Reporting Trials (CONSORT) Statement on randomized trials and it will be conducted according to the Recommendations for Interventional Trials (SPIRIT). The trial was registered in the Chinese Clinical Trial Registry (No. 2200060885). It will be carried out in the Department of Rehabilitation Medicine of Tongji Hospital in Wuhan, China. A flowchart overview of the study is presented in Figure 1. The Standard Protocol Items: SPIRIT table for enrolment, interventions, and assessments is presented in Table 1. Forty-five eligible subjects will be randomly allocated into three groups equally: (1) sham-cTBS group, (2) cTBS-cM1 group (cTBS targeting contralesional M1), (3) cTBS-cPMA group (cTBS targeting contralesional PMA).

FIGURE 1
www.frontiersin.org

Figure 1. The flow diagram for this study. MEP, motor evoked potential; MAS, modified Ashworth scale; FMA-UE, Fugl-Meyer Assessment of Upper Extremity Scale; SR, startle reflex; DTI, diffusion tensor imaging; Hslp/Mslp, the ratio of the slope of H and M-waves; Hmax/Mmax, the ratio of the maximum of H and M-waves.

TABLE 1
www.frontiersin.org

Table 1. Tabulation of enrolment, interventions, and assessments throughout the trial.

Participants

The inclusion criteria are as follows: (1) age 18–80 years old. (2) first unilateral stroke (cerebral infarction, primary intracerebral hemorrhage, subarachnoid hemorrhage) confirmed by MRI or CT within 6 months before enrolment. (3) Perceptible muscle spasticity with MAS score >1. (4) Severe motor impairment with the Fugl-Meyer Assessment of Upper Extremity (FMA-UE) score 0–35 (9, 26). (5) Failure to induce MEP in the upper extremity of the hemiplegic side (27, 28). (6) Mini-Mental State Examination (MMSE) score ≥ 22 and compliance with the interventions. (7) sign the informed consent.

The exclusion criteria are as follows: (1) more than one stroke (patients with previous transient ischemic attack could participate), bilateral cerebral hemisphere lesion. (2) history of seizures or seizures after stroke, severe skull fracture. (3) conditions unsuitable for TBS (29), e.g., intracranial implants, cardiac pacemakers or pregnancy, implanted drug pumps, etc. (4) Medications (such as lignocaine, penicillins, cephalosporins, amphotericin, tricyclics, selective serotonin reuptake inhibitors, azatadine, aminophylline, bupropion, imipramine, clozapine, olanzapine, etc.) that may lower the seizure threshold have been used continuously for nearly 3 months. (5) serious impairments of heart, lung, liver, kidney, and other organs, cannot tolerate training. (6) participation in other spasticity rehabilitation studies.

Randomization and blinding

The physician will screen out the patients who meet the eligibility criteria, explain the whole content of the trial and ask the consent to sign the informed consent. Their demographic information, medical history, and medication details will then be collected. We will use an electronic random sequence generator (www.random.org) to randomly assign all of the individuals to one of the three groups in a 1:1:1 ratio. An independent researcher will conceal the allocation using sealed and opaque envelopes, numbered consecutively. In this trial, subjects, assessors, and data analysts are blind to the allocation results. Because this study involves locating therapeutic targets in the cerebral cortex, the intervenor cannot be blind. The treatment and outcome evaluation regimen are depicted in Table 1.

Procedure

The study will consist of four phases: (1) A baseline assessment involving clinical scales (MAS and FMA-UE) and DTI before the stimulation sessions; (2) Comprehensive rehabilitation, and the corresponding cTBS interventions will be provided once a day, 5 times a week for 4 weeks. Neurophysiological assessments such as H reflex, MEP and SR will be collected before and after the first cTBS treatment at this phase; (3) Post-intervention (4 weeks) clinical, neurophysiological and DTI assessment; (4) Follow-up (8 weeks) clinical and neurophysiological assessment. Figure 1 shows the flow diagram of this study.

All stroke patients will receive a routine treatment programme including rehabilitation, medication, multidisciplinary therapy. The patient's rehabilitation included various forms of physical and occupational therapy, ~3 to 4 hours a day, 5 days a week, for 4 weeks. The intensity of active training will be appropriately adjusted according to the patient's tolerance. At the same time, an additional cTBS treatment will be provided for those participants for 4 weeks, once a day. Each patient's cTBS intervention will be performed by the same therapist.

A TMS-navigation system (Localite TMS navigator, Germany) will be used to locate the therapeutic targets, record data, and delivery the stimulation. Each patient's MRI data will be collected and imported into the navigation system to build a personalized brain model. Participants will be asked to sit in a comfortable chair with adhesive surface electrodes placed over the muscle belly and tendon of the first dorsal interosseous (FDI). The cortical target eliciting the largest responses and minimum latency of the motor evoked potential (MEP) in the first dorsal interosseous (FDI) muscle at resting status will be determined as the hotspot. The resting motor threshold (rMT) of contralesional hemisphere will be recorded accordingly, which is defined as the lowest stimulus intensity producing a MEP (peak to peak amplitude ≥ 50 μV in at least 5 of 10 trials) in the relaxed condition (30). T.M.S Motor Threshold Assessment Tool (MTAT) 2.0 (http://www.clinicalresearcher.org/software.html) will be used to assess the motor threshold (MT) (31). The location of the FDI hotspot will be set as the cTBS target for cTBS-cM1 group. As to the treatment target for cTBS-cPMA group, the dorsal portion of the superior precentral sulcus (above the superior frontal sulcus) in the contralesional hemisphere will be selected as the landmark for targeting the PMA (32). The sham group will receive cTBS with the same parameters (intensity, time) as the cTBS group on PMA of the unaffected hemisphere but with the coil rotated 90° away from the scalp so that the induced current is minimized. A transcranial magnetic stimulation (TMS MagVenture® MagPro R30, Denmark) with a cool-B65 A/P figure-of-eight-shaped coil will be applied as the stimulator, and the stimulus intensity is set to 70% of contralesional hemisphere's rMT for the three groups. The participants will receive an uninterrupted 40 s train of sham cTBS stimulation over contralesional PMA or an uninterrupted 40 s train of cTBS over contralesional M1 or PMA, respectively. The basic element is a burst ofthree stimuli at 50 Hz, which is repeated at intervals of 200 ms (i.e., 5 Hz) (24), for a total of 600 pulses. Figure 2 shows the cTBS paradigm.

FIGURE 2
www.frontiersin.org

Figure 2. Continuous theta burst stimulation (cTBS) paradigm and neurophysiological effects. The basic element of cTBS is a burst of 3 stimuli at 50 Hz, which is repeated every 200 ms. In cTBS paradigm, an uninterrupted 40 s train of TBS was delivered (600 pulses).

Tolerability and safety

Participants will be monitored during the treatment to identify any signals or adverse reactions, such as pain on stimulation site, headache, nausea, and very unlikely seizures (0.08/1,000 sessions) (33), that would require their exclusion from the study. Furthermore, any discomfort experienced by the participants will be recorded and reported.

Outcome measurements

The primary outcomes are MAS and FMA-UE. Changes before and after the 4 weeks of treatment will be used as the primary measure to assess improvement in spasticity and abnormal movement synergies (9). And the secondary outcomes are results from neurophysiological examinations including H reflex, M response, MEP and SR, and changes of neuroimaging with DTI.

Clinical assessments

Modified Ashworth Scale

The Modified Ashworth Scale (MAS) is widely used to evaluate muscle spasticity by measuring resistance during passive soft-tissue stretching. The MAS has a six-point scale (ranging from 0 to 4, respectively 0, 1, 1+, 2, 3, 4) with a score of 0 representing no spasticity and a score of 4 representing severe spasticity (34, 35). An experienced clinician will use MAS to evaluate the muscle tone of the elbow flexor, elbow extensor, and wrist flexor and extensor on the hemiplegic side of the patient and record them, respectively.

Fugl-Meyer Assessment of Upper Extremity Scale

We will also assess the mean change in the Fugl-Meyer Assessment of the Upper Extremity (FMA-UE) Scale. This is an assessment table for the upper extremity, upper arm, and wrist/hand, which examines motion, coordination, and reflexes. The FMA-UE scale has three points (respectively, 0, 1, 2 score, 0 is the worst, 2 is the best) for each item with a maximum of 66 points. It is widely used to measure post-stroke motor recovery of the upper limb because of its convenience, effectiveness, and reliability, and the higher the score, the greater the motor function (36, 37). Some of the recent refinement of flexor or extensor synergy movements and out-of-synergy movements based on the movements assessed in FMA-UE will help reflect changes in motor innervation pathways through score changes (9).

Those clinical assessments will be performed at baseline, after 4 weeks of treatment, and at follow-up.

Neurophysiological examinations

CST excitability: Motor evoked potential

TMS evoked neurophysiological parameters are useful measures for monitoring post-stroke patients and predicting recovery of motor function (38). MEP is used to investigate CST activities frequently. In this study, we will collect the contralesional hemisphere's resting motor threshold (rMT) and motor evoked potentials (MEP) from the first dorsal interosseous (FDI) muscles. All subjects will be seated in a comfortable chair with their hands completely relaxed during the evaluations. The transcranial magnetic stimulation coil will be held tangentially in a posterior anterior plane at a 45° angle from the midline over the hotspot. As mentioned before, we will primarily find the hot spots of FDI, and measure rMT. After that, MEPs will be sampled during TMS at 110% rMT and will be kept constant during the following evaluations. We will record the MEP as zero if it cannot be evoked (39).

RST excitability: Startle reflex

The subject will be instructed to sit in a chair with a backrest, with the hands resting on the thighs and remaining relaxed. Two pairs of surface electromyography electrodes were applied to the left and right sternocleidomastoid muscles (SEM) belly record muscle response. After wearing a headphone (Sennheiser HD25-I; Wedemark, Germany), the subjects will receive a total of 30 sound stimuli. Half of them are 40 ms beeps at 80 dB and half are 40 ms white noises at 120 dB. The stimulus firing sequence was randomly generated, with a random interval of 6–10 s between stimuli to ensure subjects were fully relaxed. The Psychtoolbox-3 toolkit within MATLAB (2017b, MathWorks, Natick, MA, USA) will be used for the design and execution of this trial. The sEMG signals will be collected by the Ultimu EMG system (Noraxon USA Inc., Scottsdale, AZ, USA) at a sampling rate of 2,000 Hz. Raw sEMG data will be output and processed under the same software. It will be segmented, bandpass filtered at 30–300 Hz, and notch filtered at 50 Hz. A threshold detection method (mean with 3SD) based on the amplitudes at time window 2,500–500 ms before each stimulus will be set to determine the onset of muscles (40). A positive startle response will be screened out if the onset of either SEM at the time widow of 30–130 ms after stimuli. The positive rate, as well as SEM onset delays (41) will be recorded for analysis.

Motor neuron pool excitability: H reflex and M response

Hslp/Mslp (the ratio of the slope of H and M-waves) denotes the relationship between stimulus intensity and motor neuron reflex recruitment (42), and previous researches have proposed Hslp/Mslp and Hmax/Mmax (the ratio of the maximum of H and M-waves) as sensitive and objective measures for the motor neuron pool excitability (42, 43). Poststroke patients have poor autonomic activity, which is thought to contribute to increased reflex excitability manifested by greater Hslp and Hmax on the paretic side (4446).

A physician will examine all patients with Dantec, two channel EMG. All subjects will be seated comfortably in a chair with their hands resting on thighs and elbows flexed 20°-30°. Paired surface electrodes will be attached to the skin on the belly of the flexor carpi radialis (FCR) muscle. We will stimulate (1 ms duration, once every 3 s) the median nerve at medial to the bicep tendon in the antecubital fossa to elicit H reflexes and M-responses from the flexor carpi radialis muscle. The intensity of these pulses will increase gradually from below the threshold of the H-reflex to supramaximal and saturated condition for M response. At each stimulus intensity, we will deliver 3–5 pulses successively. The H-reflexes and the M responses evoked by these pulses will be amplified in the bandwidth 5 Hz−3 kHz. All amplification procedures are controlled by software running on a signal processor. All signals will be checked and derived to a hard disk for offline analyses. All data regarding the size and series stimulus intensities of both the H-reflex and M response will be used in the construction of the respective recruitment curves and to calculate the average of their maximum values.

Electrophysiological examinations will be performed before and after the first cTBS treatment, at post-intervention (at 4 weeks), and at follow-up (at 8 weeks).

Neuroimaging: diffusion tensor imaging

Diffusion Tensor Imaging (DTI) will be evaluated at baseline and after the 4-week treatment. And it will be conducted as follows. All individuals will be subjected to 3.0 T MRI scanner (Discovery LS MR 750; GE Healthcare, Chicago, IL, USA) equipped with a 16-channel head & neck coil to acquire 3D T1-weighted and DTI images before and after 20 times therapy.

DTI is used in diffusion-weighted gradients along 60 non-collinear directions and two non-diffusion-weighted volumes. Scanning parameters will include matrix size = 128 × 128, field of view = 230 × 230 mm2, TR = 5,800 ms, TE = 71 ms, flip angle = 180°, number of averages = 1, b value = 1,000 s/mm2, slice thickness = 4 mm, and voxel size = 1.7 × 1.7 × 4 mm3 (47).

DTI data processing will be carried out as follows. (1) Normalization: All data will be spatially normalized using parameters derived from 3D T1WI processing. We will use FMRIB's Diffusion Toolbox and TBSS (Tract-Based Spatial Statistics) (both from the FMRIB Software Library 6.0, FSL, https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FSL) to analyze DTI data. (2) Registration: We will use FSL tools to get the transformation matrix and then correct the eddy current distortions in the DTI datasets. After that, each subject's raw diffusion-weighted pictures will be linearly aligned to a non-diffusion weighted image (b0) with the removal of the non-brain tissues using a brain extraction tool (BET). (3) Local fitting: The extracted brain will be used for diffusion tensors. The diffusion tensor will be estimated at each voxel to obtain the fractional anisotropy (FA) pictures (47). We will select the region of interest (ROI) method to measure tract integrity. The red nucleus (RN) is an appropriate part to draw ROI (7). We will draw ROI on the posterior limb of internal capsule and the anterior pons on a color-coded map for CST and on the RN for RST individually. In each ROI, we will extract fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD) and axial diffusivity (AD). Three dimensional reconstructions of the CST and RST will be performed for each patient.

DTI has been commonly used to quantitatively measure the integrity of tissues. Previous studies reported the DTI biomarkers (e.g., FA) measured after stroke have emerged as potential predictors of motor recovery (48, 49).

Sample size calculation

Sample size is estimated through software G*Power 3.1 (50). We used MAS as the primary outcome for the calculation of the effect size. It is determined to be 12 per group based on a power of 0.8 with α level at 0.05. According to a previous study (21), effect size d = 0.53 was calculated based on the difference in MAS score from their main research results. Through conversion, the effect size f used in our study is about 0.25. Repeated measures correlation in our study is set at 0.5 to detect a moderate within-between effect of Conhen's f = 0.25. Considering a drop-out rate of 20%, we aim to recruit totally 45 stroke subjects.

Data monitoring and management

The data monitoring committee (DMC) 's responsibilities include safety monitoring, test data monitoring and test design adjustment recommendations with no conflicts of interest. Rehabilitation Department of Tongji Hospital will be in charge of informed consent quality assurance, eligible participant recruiting, implementation of interventions, and data administration. Designated person will be responsible for collecting case report forms (CRFs) as well as data transfer and analysis. The investigator and director are responsible for retaining all records, and the data center will keep anonymized case report form data. Electronic data will be preserved on password-protected computers, while all paper data will be kept in secure filing cabinets.

Statistical analysis

The SPSS 21.0 program will be used for all statistical analysis and the significance level set for all the analysis will be P ≤ 0.05. The Shapiro-Wilk test will be used (P > 0.05) to check normality. In this study, Continuous variables will be expressed as mean with standard deviation and median for non-normal distributions and categorical variables will be presented as percentage or frequencies. Between the three groups, data will be first compared using Student's t-test for continuous variables and the chi-square or Fisher's exact tests for categorical variables. Appropriate metrics will be processed further using repeated measures ANOVA. Bonferroni correction method and False Discovery Rate <5% will be considered for multiple comparisons. All analyses will be done with the intention to treat principle, and missing data will be estimated using multiple imputation. The mean comparisons for related outcomes in two different moments will be studied with the Wilcoxon-signed rank test. The major strategy for dealing with missing data will be a sensitivity analysis and weighted estimating equations.

Discussion

The motor recovery after stroke follows a basically predictable pattern, from flaccid, spastic to recovered (4, 10). In patients with severe stroke, motor function may be arrested at a certain stage limited by spasticity for a long time in the process of motor recovery (4), which is a thorny problem in the rehabilitation process. Motor rehabilitation relies on a complex combination of spontaneous recovery and motor learning to promote neural plasticity (11, 51). The spinal plastic change relates to the severity of cortical damage and has a time-limited window, almost entirely in the first weeks (52, 53). The anatomical structure determines that in patients with severe stroke, CST is more vulnerable to damage, ultimately leading to the facilitatory medial RST being unopposed and independent of cortical control (13), which plays a significant role in spasticity. Spasticity reflects a phenomenon of maladaptive neuroplasticity in motor control pathways (12), and can easily develop if there is not appropriate intervention to regulate in the course of recovery (54).

Non-invasive magnetic stimulation is a promising therapy option for modulating neural plasticity, but the undetected MEP (27) and the extensive recombination of the cerebral cortex (55) after severe stroke lead to doubtful targeting in the cortex for treatment. Although there is weak-grade evidence to support that the contralesional TMS targeting cM1 improves post-stroke spasticity, the effect is inconclusive (54). Exploring the spasticity mechanism in severe stroke patients may be more helpful to obtain optimized intervention targets and clarify the effect of this non-invasive treatment.

As mentioned previously, spasticity is thought to be a top-down manifestation resulting from hyperexcitability of the advance central PMA and its corresponding descending pathway of medial RST (12). This study will compare the acute and short-term effects of cTBS over cPMA and M1 in severe stroke patients. Moreover, we aim to demonstrate that the hyperexcitable cPMA after stroke is the cortical manifestation causing spasticity, and reducing the hyperexcitability of cPMA may help to down regulating the excitability mRST of hemiplegic side and to facilitate adaptive neuroplasticity of RST in spasticity relief and motor recovery.

Several characteristics distinguish this study from previous reports. Firstly, it is the first clinical study that TBS improves spasticity by regulating the neuroplasticity of cPMA and mRST. The results will verify the relationship between spasticity, cPMA and mRST, and provide theoretical and practical support for effective improvement of spasticity. Otherwise, the changes of neuroplasticity and spasticity after stroke will be evaluated objectively and comprehensively from clinical, electrophysiological and neuroimaging perspectives. In this study, we use the FMA-UE scale to monitor the recovery course of the upper extremity's motor function (56). Moreover, subsections and items of FMA-UE such as synergy or out-of-synergy movements can roughly reveal the pathways (CST/RST) that patients currently dominate their motor and spastic recovery (9). Some personalized rehabilitation decisions based on this feature may help develop more effective treatment plans. DTI biomarkers of CST have been considered as predictors of motor recovery (48). In this study, DTI will be used to evaluate the neuroplasticity changes of both RST (47) and CST (48) in post-stroke spasticity. The results will be used to further clarify the dynamic competitive relationship between the both (9). We expect that the results of this study can provide powerful support for the application of TBS neuromodulation in improving spasticity, and we may supplement targeted training [motor training, auditory stimulation training (57), etc.] in different stages of motor recovery in the near future, so that stroke patients can reasonably utilize CST and RST pathways to facilitate motor function recovery.

Ethics statement

The studies involving human participants were reviewed and approved by Institutional Ethical Committee of Tongji Hospital. The patients/participants provided their written informed consent to participate in this study.

Author contributions

YL, WG, NX, and XW contributed to the conception and design of the study. XW and NX wrote the first draft of the manuscript and generated the figures. All authors contributed to the refinement of the study protocol and approved the final manuscript.

Acknowledgments

The authors gratefully acknowledge the contribution of medical staffs in the Department of Rehabilitation Medicine of Tongji Hospital for their support and assistance for the study initiation and implementation.

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.

Abbreviations

Abbreviations: CST, corticospinal tract; RST, reticulospinal tract; cPMA, contralesional premotor area; cM1, contralesional primary motor cortex; mRST, the medial reticulospinal tract; LTD, long-term depression; FMA-UE, Fugl-Meyer Assessment of Upper Extremity; MEP, motor evoked potential; SR, startle reflex; MAS, Modified Ashworth scale; DTI, Diffusion tensor imaging; Hslp/Mslp, the ratio of the slope of H and M-waves; Hmax/Mmax, the ratio of the maximum of H and M-waves; RN, red nucleus.

References

1. Li S, Francisco GE, Rymer WZ. A new definition of poststroke spasticity and the interference of spasticity with motor recovery from acute to chronic stages. Neurorehabil Neural Repair. (2021) 35:601–10. doi: 10.1177/15459683211011214

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Kuo C-L, Hu G-C. Post-stroke spasticity: a review of epidemiology, pathophysiology, and treatments. Int J Gerontol. (2018) 12:280–4. doi: 10.1016/j.ijge.2018.05.005

CrossRef Full Text | Google Scholar

3. Bavikatte G, Bavikatte TG, Gaber T. Approach to spasticity after stroke in general practice. Br J Med Pract. (2009) 2:29–34. Available online at: https://www.bjmp.org/content/approach-spasticity-general-practice

PubMed Abstract | Google Scholar

4. Thomas ET. The restoration of motor function following hemiplegia in brain. Brain. (1951) 74:443–80. doi: 10.1093/brain/74.4.443

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Nudo RJ. Mechanisms for recovery of motor function following cortical damage. Curr Opin Neurobiol. (2006) 16:638–44. doi: 10.1016/j.conb.2006.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Baker SN. The primate reticulospinal tract, hand function and functional recovery. J Physiol. (2011) 589(Pt 23):5603–12. doi: 10.1113/jphysiol.2011.215160

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Kim H, Lee H, Jung KI, Ohn SH, Yoo WK. Changes in diffusion metrics of the red nucleus in chronic stroke patients with severe corticospinal tract injury: a preliminary study. Ann Rehabil Med. (2018) 42:396–405. doi: 10.5535/arm.2018.42.3.396

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Karbasforoushan H, Cohen-Adad J, Dewald JPA. Brainstem and spinal cord MRI identifies altered sensorimotor pathways post-stroke. Nat Commun. (2019) 10:3524. doi: 10.1038/s41467-019-11244-3

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Senesh MR, Barragan K, Reinkensmeyer DJ. Rudimentary dexterity corresponds with reduced ability to move in synergy after stroke: evidence of competition between corticoreticsulospinal and corticospinal tracts? Neurorehabil Neural Repair. (2020) 34:904–14. doi: 10.1177/1545968320943582

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Brunnstrom S. Motor testing procedures in hemiplegia: based on sequential recovery stages. Phys Ther. (1966) 46:357–75. doi: 10.1093/ptj/46.4.357

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Li S. Spasticity, motor recovery, and neural plasticity after stroke. Front Neurol. (2017) 8:120. doi: 10.3389/fneur.2017.00120

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Li S, Chen YT, Francisco GE, Zhou P, Rymer WZ. A Unifying pathophysiological account for post-stroke spasticity and disordered motor control. Front Neurol. (2019) 10:468. doi: 10.3389/fneur.2019.00468

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Li S, Francisco GE. New insights into the pathophysiology of post-stroke spasticity. Front Hum Neurosci. (2015) 9:192. doi: 10.3389/fnhum.2015.00192

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Pundik S, Falchook AD, McCabe J, Litinas K, Daly JJ. Functional brain correlates of upper limb spasticity and its mitigation following rehabilitation in chronic stroke survivors. Stroke Res Treat. (2014) 2014:306325. doi: 10.1155/2014/306325

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Fregosi M, Contestabile A, Hamadjida A, Rouiller EM. Corticobulbar projections from distinct motor cortical areas to the reticular formation in Macaque monkeys. Eur J Neurosci. (2017) 45:1379–95. doi: 10.1111/ejn.13576

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kinnear BZ, Lannin NA, Cusick A, Harvey LA, Rawicki B. Rehabilitation therapies after botulinum toxin-a injection to manage limb spasticity a systematic review. Phys Ther. (2014) 94:1569–81. doi: 10.2522/ptj.20130408

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Etoh S, Noma T, Ikeda K, Jonoshita Y, Ogata A, Matsumoto S, et al. Effects of repetitive trascranial magnetic stimulation on repetitive facilitation exercises of the hemiplegic hand in chronic stroke patients. J Rehabil Med. (2013) 45:843–7. doi: 10.2340/16501977-1175

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Xu P, Huang Y, Wang J, An X, Zhang T, Li Y, et al. Repetitive transcranial magnetic stimulation as an alternative therapy for stroke with spasticity: a systematic review and meta-analysis. J Neurol. (2021) 268:4013–22. doi: 10.1007/s00415-020-10058-4

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Dos Santos RBC, Galvao SCB, Frederico LMP, Amaral NSL, Carneiro MIS, de Moura Filho AG, et al. Cortical and spinal excitability changes after repetitive transcranial magnetic stimulation combined to physiotherapy in stroke spastic patients. Neurol Sci. (2019) 40:1199–207. doi: 10.1007/s10072-019-03765-y

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Liu Y, Li H, Zhang J, Zhao QQ, Mei HN, Ma J, et al. Meta-analysis: whether repetitive transcranial magnetic stimulation improves dysfunction caused by stroke with lower limb spasticity. Evid Based Complement Alternat Med. (2021) 2021:7219293. doi: 10.1155/2021/7219293

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Barros, Galvao SC. Borba Costa dos Santos R, Borba dos Santos P, Cabral ME, Monte-Silva K. Efficacy of coupling repetitive transcranial magnetic stimulation and physical therapy to reduce upper-limb spasticity in patients with stroke: a randomized controlled trial. Arch Phys Med Rehabil. (2014) 95:222–9. doi: 10.1016/j.apmr.2013.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Daskalakis ZJ, Christensen BK, Fitzgerald PB, Roshan L, Chen R. The mechanisms of interhemispheric inhibition in the human motor cortex. J Physiol. (2002) 543(Pt 1):317–26. doi: 10.1113/jphysiol.2002.017673

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Rounis E, Huang YZ. Theta burst stimulation in humans: a need for better understanding effects of brain stimulation in health and disease. Exp Brain Res. (2020) 238:1707–14. doi: 10.1007/s00221-020-05880-1

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Huang YZ, Edwards MJ, Rounis E, Bhatia KP, Rothwell JC. Theta burst stimulation of the human motor cortex. Neuron. (2005) 45:201–6. doi: 10.1016/j.neuron.2004.12.033

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Chu HT, Cheng CM, Liang CS, Chang WH, Juan CH, Huang YZ, et al. Efficacy and tolerability of theta-burst stimulation for major depression: a systematic review and meta-analysis. Prog Neuropsychopharmacol Biol Psychiatry. (2021) 106:110168. doi: 10.1016/j.pnpbp.2020.110168

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Woytowicz EJ, Rietschel JC, Goodman RN, Conroy SS, Sorkin JD, Whitall J, et al. Determining levels of upper extremity movement impairment by applying a cluster analysis to the fugl-meyer assessment of the upper extremity in chronic stroke. Arch Phys Med Rehabil. (2017) 98:456–62. doi: 10.1016/j.apmr.2016.06.023

PubMed Abstract | CrossRef Full Text | Google Scholar

27. van Kuijk AA, Pasman JW, Hendricks HT, Zwarts MJ, Geurts AC. Predicting hand motor recovery in severe stroke,the role of motor evoked potentials in relation to early clinical assessment. Neurorehabil Neural Repair. (2009) 23:45–51. doi: 10.1177/1545968308317578

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lai C-J, Wang C-P, Tsai P-Y, Chan R-C, Lin S-H, Lin F-G, et al. Corticospinal integrity and motor impairment predict outcomes after excitatory repetitive transcranial magnetic stimulation: a preliminary study. Arch Phys Med Rehabil. (2015) 96:69–75. doi: 10.1016/j.apmr.2014.08.014

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Eric MW. Risk and safety of repetitive transcranial magnetic stimulation:report and suggested guidelines from the International Workshop on the Safety of Repetitive Transcranial Magnetic Stimulation, June 5–7, 1996. Electroencephalogr Clin Neurophysiol. (1998) 108:1–16. doi: 10.1016/s0168-5597(97)00096-8

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Matsuura A, Karita T, Nakaso N, Kondo Y, Mori F. Mirror movement-like muscle hypertonia stroke model based on ipsilateral weight load. Neuroreport. (2020) 31:251–5. doi: 10.1097/WNR.0000000000001403

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Awiszus F. Tms and threshold hunting. Suppl Clin Neurophysiol. (2003) 56:13–23. doi: 10.1016/S1567-424X(09)70205-3

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Ahdab R, Ayache SS, Farhat WH, Mylius V, Schmidt S, Brugieres P, et al. Reappraisal of the anatomical landmarks of motor and premotor cortical regions for image-guided brain navigation in tms practice. Hum Brain Mapp. (2014) 35:2435–47. doi: 10.1002/hbm.22339

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lerner AJ, Wassermann EM, Tamir DI. Seizures from transcranial magnetic stimulation 2012-2016: results of a survey of active laboratories and clinics. Clin Neurophysiol. (2019) 130:1409–16. doi: 10.1016/j.clinph.2019.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Naghdi S, Ansari NN, Mansouri K, Olyaei GR, Asgari A, Kazemnejad A. The correlation between modified ashworth scale scores and the new index of alpha motoneurones excitability in post-stroke patients. Electromyogr Clin Neurophysiol. (2008) 48:109–15.

PubMed Abstract | Google Scholar

35. Ewoldt JK, Lazzaro EC, Roth EJ, Suresh NL. Quantification of a single score (1+) in the Modified Ashworth Scale (Mas), a Clinical Assessment of Spasticity. Annu Int Conf IEEE Eng Med Biol Soc. (2016) 2016:1737–40. doi: 10.1109/EMBC.2016.7591052

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kim H, Her J, Ko J, Park D, Woo JH, You Y, et al. Reliability, concurrent validity, and responsiveness of the Fugl-Meyer Assessment (Fma) for Hemiplegic Patients. J Phys Ther Sci. (2012) 24:839–99. doi: 10.1589/jpts.24.893

CrossRef Full Text | Google Scholar

37. See J, Dodakian L, Chou C, Chan V, McKenzie A, Reinkensmeyer DJ, et al. A standardized approach to the Fugl-Meyer assessment and its implications for clinical trials. Neurorehabil Neural Repair. (2013) 27:732–41. doi: 10.1177/1545968313491000

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Cakar E, Akyuz G, Durmus O, Bayman L, Yagci I, Karadag-Saygi E, et al. The relationships of motor-evoked potentials to hand dexterity, motor function, and spasticity in chronic stroke patients: a transcranial magnetic stimulation study. Acta Neurol Belg. (2016) 116:481–7. doi: 10.1007/s13760-016-0633-2

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Stinear C. Prediction of recovery of motor function after stroke. Lancet Neurol. (2010) 9:1228–32. doi: 10.1016/S1474-4422(10)70247-7

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Xia N, He C, Li YA, Gu M, Chen Z, Wei X, et al. Startle increases the incidence of anticipatory muscle activations but does not change the task-specific muscle onset for patients after subacute stroke. Front Neurol. (2021) 12:789176. doi: 10.3389/fneur.2021.789176

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Rangarajan V, Schreiber JJ, Barragan B, Schaefer SY, Honeycutt CF. Delays in the reticulospinal system are associated with a reduced capacity to learn a simulated feeding task in older adults. Front Neural Circuits. (2021) 15:681706. doi: 10.3389/fncir.2021.681706

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Funase K, Higashi T, Yoshimura T, Imanaka K, Nishihira Y. Evident difference in the excitability of the motoneuron pool between normal subjects and patients with spasticity assessed by a new method using H reflex and M response. Neurosci Lett. (1996) 203:127–30. doi: 10.1016/0304-3940(95)12284-2

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Walton C, Kalmar J, Cafarelli E. Caffeine increases spinal excitability in humans. Muscle Nerve. (2003) 28:359–64. doi: 10.1002/mus.10457

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Lundbye-Jensen J, Nielsen JB. Central nervous adaptations following 1 wk of wrist and hand immobilization. J Appl Physiol (1985). (2008) 105:139–51. doi: 10.1152/japplphysiol.00687.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Clark BC, Manini TM, Bolanowski SJ, Ploutz-Snyder LL. Adaptations in human neuromuscular function following prolonged unweighting: ii. neurological properties and motor imagery efficacy. J Appl Physiol (1985). (2006) 101:264–72. doi: 10.1152/japplphysiol.01404.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Phadke CP, Robertson CT, Condliffe EG, Patten C. Upper-extremity H-reflex measurement post-stroke: reliability and inter-limb differences. Clin Neurophysiol. (2012) 123:1606–15. doi: 10.1016/j.clinph.2011.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Ko SH, Kim T, Min JH, Kim M, Ko HY, Shin YI. Corticoreticular pathway in post-stroke spasticity: a diffusion tensor imaging study. J Pers Med. (2021) 11:1151. doi: 10.3390/jpm11111151

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Moura LM, Luccas R, Paiva JPQ, Amaro E, Leemans A, Leite CDC, et al. Diffusion tensor imaging biomarkers to predict motor outcomes in stroke: a narrative review. Front Neurol. (2019) 10:445. doi: 10.3389/fneur.2019.00445

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Soares JM, Marques P, Alves V, Sousa N. A Hitchhiker's guide to diffusion tensor imaging. Front Neurosci. (2013) 7:31. doi: 10.3389/fnins.2013.00031

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Kang H. Sample size determination and power analysis using the G*Power software. J Educ Eval Health Prof. (2021) 18:17. doi: 10.3352/jeehp.2021.18.17

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Langhorne P, Bernhardt J, Kwakkel G. Stroke rehabilitation. Lancet. (2011) 377:1693–702. doi: 10.1016/S0140-6736(11)60325-5

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Sist B, Fouad K, Winship IR. Plasticity beyond peri-infarct cortex: spinal up regulation of structural plasticity, neurotrophins, and inflammatory cytokines during recovery from cortical stroke. Exp Neurol. (2014) 252:47–56. doi: 10.1016/j.expneurol.2013.11.019

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Wiersma AM, Fouad K, Winship IR. Enhancing spinal plasticity amplifies the benefits of rehabilitative training and improves recovery from stroke. J Neurosci. (2017) 37:10983–97. doi: 10.1523/JNEUROSCI.0770-17.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Hatem SM, Saussez G, Della Faille M, Prist V, Zhang X, Dispa D, et al. Rehabilitation of motor function after stroke: a multiple systematic review focused on techniques to stimulate upper extremity recovery. Front Hum Neurosci. (2016) 10:442. doi: 10.3389/fnhum.2016.00442

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Grefkes C, Ward NS. Cortical reorganization after stroke: how much and how functional? Neuroscientist. (2014) 20:56–70. doi: 10.1177/1073858413491147

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Gladstone DJ, Danells CJ, Black SE. The fugl-meyer assessment of motor recovery after stroke: a critical review of its measurement properties. Neurorehabil Neural Repair. (2002) 16:232–40. doi: 10.1177/154596802401105171

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Aluru V, Lu Y, Leung A, Verghese J, Raghavan P. Effect of auditory constraints on motor performance depends on stage of recovery post-stroke. Front Neurol. (2014) 5:106. doi: 10.3389/fneur.2014.00106

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: spasticity, stroke, premotor area, theta burst stimulation (TBS), H reflex, randomized controlled trial, diffusion tensor imaging (DTI)

Citation: Wei X, Xia N, Li Y-A, Gu M, Zhang T, Gao W and Liu Y (2022) Immediate and short-term effects of continuous theta burst transcranial magnetic stimulation over contralesional premotor area on post-stroke spasticity in patients with severe hemiplegia: Study protocol for a randomized controlled trial. Front. Neurol. 13:895580. doi: 10.3389/fneur.2022.895580

Received: 14 March 2022; Accepted: 10 August 2022;
Published: 23 August 2022.

Edited by:

Antonio Suppa, Sapienza University of Rome, Italy

Reviewed by:

Kenneth N. K. Fong, Hong Kong Polytechnic University, Hong Kong SAR, China
Tsung-Hsun Hsieh, Chang Gung University, Taiwan

Copyright © 2022 Wei, Xia, Li, Gu, Zhang, Gao and Liu. 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: Yali Liu, bGl1eWFsaTc0JiN4MDAwNDA7MTI2LmNvbQ==; Wei Gao, Z2FvYnVsbCYjeDAwMDQwOzEyNi5jb20=

These authors have contributed equally to this work and share first authorship

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.