Time-frequency patterns of somatosensory evoked potentials in predicting the location of spinal cord injury

Identifying Intraoperative Spinal Cord Injury Location from Somatosensory Evoked Potentials' Time-Frequency Components

Excessive distraction in corrective spine surgery can lead to iatrogenic distraction spinal cord injury. Diagnosis of the location of the spinal cord injury helps in early removal of the injury source. The time-frequency components of the somatosensory evoked potential have been reported to provide information on the location of spinal cord injury, but most studies have focused on contusion injuries of the cervical spine. In this study, we established 19 rat models of distraction spinal cord injury at different levels and collected the somatosensory evoked potentials of the hindlimb and extracted their time-frequency components. Subsequently, we used k-medoid clustering and naive Bayes to classify spinal cord injury at the C5 and C6 level, as well as spinal cord injury at the cervical, thoracic, and lumbar spine, respectively. The results showed that there was a significant delay in the latency of the time-frequency components distributed between 15 and 30 ms and 50 and 150 Hz in all spinal cord injury groups. The overall classification accuracy was 88.28% and 84.87%. The results demonstrate that the k-medoid clustering and naive Bayes methods are capable of extracting the time-frequency component information depending on the spinal cord injury location and suggest that the somatosensory evoked potential has the potential to diagnose the location of a spinal cord injury.

Identification of injury type using somatosensory and motor evoked potentials in a rat spinal cord injury model

The spinal cord is at risk of injury during spinal surgery. If intraoperative spinal cord injury is identified early, irreversible impairment or loss of neurological function can be prevented. Different types of spinal cord injury result in damage to different spinal cord regions, which may cause different somatosensory and motor evoked potential signal responses. In this study, we examined electrophysiological and histopathological changes between contusion, distraction, and dislocation spinal cord injuries in a rat model. We found that contusion led to the most severe dorsal white matter injury and caused considerable attenuation of both somatosensory and motor evoked potentials. Dislocation resulted in loss of myelinated axons in the lateral region of the injured spinal cord along the rostrocaudal axis. The amplitude of attenuation in motor evoked potential responses caused by dislocation was greater than that caused by contusion. After distraction injury, extracellular spaces were slightly but not significantly enlarged; somatosensory evoked potential responses slightly decreased and motor evoked potential responses were lost. Correlation analysis showed that histological and electrophysiological findings were significantly correlated and related to injury type. Intraoperative monitoring of both somatosensory and motor evoked potentials has the potential to identify iatrogenic spinal cord injury type during surgery.

Different Time-frequency Distribution Patterns of Somatosensory Evoked Potentials in Dual- and Single-level Spinal Cord Compression

Among patients with cervical myelopathy, the most common level of stenosis at spinal cord of all ages was reported to be between cervical levels C5-6. Previous studies found that time-frequency components (TFCs) of somatosensory evoked potentials (SEPs) possess location information of spinal cord injury (SCI) in single-level deficits in the spinal cord. However, the clinical reality is that there are multiple compressions at multiple spinal cord segments. This study proposed a new algorithm to differentiate distribution patterns of SEP TFCs between the dual-level compression and the corresponding single-level compression, which is potentials in providing precise diagnosis of cervical myelopathy. In the present animal study, a group of rats with dual-level compressive (C5+6) injury to cervical spinal cord was investigated. SEPs were collected at 2 weeks after surgery, while SEP TFCs were calculated. The SEP TFCs under dual-level compression were compared to an existent dataset with one sham control group and three single level compression groups at C4, C5, C6. Behavioral evaluation showed very similar scale of injury severity between individual rats, while histology evaluation confirmed the precise location of injury. According to time-frequency distribution patterns, it showed that the middle-energy components of dual-level showed similar patterns as that of each single-level group. In addition, the low-energy components of the dual-level C5+6 group had the highest correlation with C5 (R=0.3423, p<0.01) and C6 (R=0.4000, p<0.01) groups, but much lower with C4 group (R=0.1071, p=0.012). These results indicated that SEP TFCs components possess information regarding the location of neurological lesion after spinal cord compression. It preliminarily demonstrated that SEP TFCs are likely a useful measure to provide location information of neurological lesions after compression SCI.

Interhemispheric differences in time-frequency representation of motor evoked potentials in brain tumor patients

OBJECTIVE

Conventional time-series parameters are unreliable descriptors of motor-evoked potentials (MEPs) in brain tumor patients. Frequency domain analysis is suggested to provide additional information about the status of the cortico-spinal motor system. Aim of the present study was to describe the time-frequency representation of MEPs and its relation to the motor performance.

METHODS

This study enrolled 17 consecutive brain tumor patients with impaired dexterity. After brain mapping of the affected (AH) and non-affected (NAH) hemisphere, TMS was applied to the hotspots of the abductor pollicis brevis muscles (APB). Using a Morlet wavelet approach, event-related spectral perturbation (ERSP) and inter-trial coherence (ITC) of the MEPs were calculated and compared to the Grooved Pegboard Test (GPT). Additionally, inter- and intra-subject reliability was assessed by the intraclass correlation coefficient (ICC).

RESULTS

MEPs were projecting to a frequency band between 30 and 400 Hz with a local maximum between 100 and 150 Hz. There was a significant ERSP and ITC reduction of the AH in comparison to the NAH. In contrast, no interhemispheric differences were depicted in the conventional time-series analysis. ERSP and ITC values correlated significantly with GPT results (r = 0.35 and r = 0.50). Time-frequency MEP description had good inter-and intra-subject reliability (ICC = 0.63).

CONCLUSIONS

Brain tumors affect corticospinal transmission resulting in a reduction of temporal and spectral MEP synchronization correlating with the dexterity performance.

SIGNIFICANCE

Time-frequency representation of MEPs provide additional information beyond conventional time-domain features.

Time-Frequency Representation of Motor Evoked Potentials in Brain Tumor Patients

Background: The integrity of the motor system can be examined by applying navigated transcranial magnetic stimulation (nTMS) to the cortex. The corresponding motor-evoked potentials (MEPs) in the target muscles are mirroring the status of the human motor system, far beyond corticospinal integrity. Commonly used time domain features of MEPs (e.g., peak-to-peak amplitudes and onset latencies) exert a high inter-subject and intra-subject variability. Frequency domain analysis might help to resolve or quantify disease-related MEP changes, e.g., in brain tumor patients. The aim of the present study was to describe the time-frequency representation of MEPs in brain tumor patients, its relation to clinical and imaging findings, and the differences to healthy subject.

Methods: This prospective study compared 12 healthy subjects with 12 consecutive brain tumor patients (with and without a paresis) applying nTMS mapping. Resulting MEPs were evaluated in the time series domain (i.e., amplitudes and latencies). After transformation into the frequency domain using a Morlet wavelet approach, event-related spectral perturbation (ERSP), and inter-trial coherence (ITC) were calculated and compared to diffusion tensor imaging (DTI) results.

Results: There were no significant differences in the time series characteristics between groups. MEPs were projecting to a frequency band between 30 and 300 Hz with a local maximum around 100 Hz for both healthy subjects and patients. However, there was ERSP reduction for higher frequencies (>100 Hz) in patients in contrast to healthy subjects. This deceleration was mirrored in an increase of the inter-peak MEP latencies. Patients with a paresis showed an additional disturbance in ITC in these frequencies. There was no correlation between the CST integrity (as measured by DTI) and the MEP parameters.

Conclusion: Time-frequency analysis may provide additional information above and beyond classical MEP time domain features and the status of the corticospinal system in brain tumor patients. This first evaluation indicates that brain tumors might affect cortical physiology and the responsiveness of the cortex to TMS resulting in a temporal dispersion of the corticospinal transmission.

A translational study of somatosensory evoked potential time-frequency components in rats, goats, and humans

Somatosensory evoked potentials (SEPs) have been widely used to assess neurological function in clinical practice. A good understanding of the association between SEP signals and neurological function is helpful for precise diagnosis of impairment location. Previous studies on SEPs have been reported in animal models. However, few studies have reported the relationships between SEP waveforms in animals and those in humans. In this study, we collected normal SEP waveforms and decomposed them into specific time-frequency components (TFCs). Our results showed three stable TFC distribution regions in intact goats and rats and in humans. After we induced spinal cord injury in the animal models, a greater number of small TFC distribution regions were observed in the injured goat and rat groups than in the normal group. Moreover, there were significant correlations (P < 0.05) and linear relationships between the main SEP TFCs of the human group and those of the goat and rat groups. A stable TFC distribution of SEP components was observed in the human, goat and rat groups, and the TFC distribution modes were similar between the three groups. Results in various animal models in this study could be translated to future clinical studies based on SEP TFC analysis. Human studies were approved by the Institutional Review Board of the University of Hong Kong/Hospital Authority Hong Kong West Cluster (approval No. UM 05-312 T/975) on December 5, 2005. Rat experiments were approved by the Committee on the Use of Live Animals in Teaching and Research of Li Ka Shing Faculty of Medicine of the University of Hong Kong (approval No. CULART 2912-12) on January 28, 2013. Goat experiments were approved by the Animal Ethics Committee of Affiliated Hospital of Guangdong Medical University (approval No. GDY2002132) on March 5, 2018.

Exploration of Cervical Myelopathy Location From Somatosensory Evoked Potentials Using Random Forests Classification

Studies using time-frequency analysis have reported that somatosensory evoked potentials provide information regarding the location of spinal cord injury. However, a better understanding of the time-frequency components derived from somatosensory evoked potentials is essential for developing more reliable algorithms that can diagnosis level (location) of cervical injury. In the present study, we proposed a random forests machine learning approach, for separating somatosensory evoked potentials depending on spinal cord state. For data acquisition, we established rat models of compression spinal cord injury at the C4, C5, and C6 levels to induce cervical myelopathy. After making the compression injury, we collected somatosensory evoked potentials and extracted their time-frequency components. We then used the random forests classification system to analyze the evoked potential dataset that was obtained from the three groups of model rats. Evaluation of the classifier performance revealed an overall classification accuracy of 84.72%, confirming that the random forests method was able to separate the time-frequency components of somatosensory evoked potentials from rats under different conditions. Features of the time-frequency components contained information that could identify the location of the cervical spinal cord injury, demonstrating the potential benefits of using time-frequency components of somatosensory evoked potentials to diagnose the level of cervical injury in cervical myelopathy.

Sorafenib promotes sensory conduction function recovery via miR-142-3p/AC9/cAMP axis post dorsal column injury

Spinal cord injury results in sensation dysfunction. This study explored miR-142-3p, which acts a critical role in sciatic nerve conditioning injury (SNCI) promoting the repair of the dorsal column injury and validated its function on primary sensory neuron(DRG). miR-142-3p expression increased greatly in the spinal cord dorsal column lesion (SDCL) group and increased slightly in the SNCI group. Subsequently, the expression of adenylate cyclase 9 (AC9), the target gene of miR-142-3p, declined sharply in the SDCL group and declined limitedly in the SNCI group. The expression trend of cAMP was opposite to that of miR-142-3p. MiR-142-3p inhibitor improved the axon length, upregulated the expression of AC9, cAMP, p-CREB, IL-6, and GAP43, and downregulated the expression of GTP-RhoA. miR-142-3p inhibitor combined with AC9 siRNA showed shorter axon length, the expression of AC9, cAMP, p-CREB, IL-6, and GAP43 was decreased, and the expression of GTP-RhoA was increased. H89 and AG490, inhibitors of cAMP/PKA pathway and IL6/STAT3/GAP43 axis, respectively, declined the enhanced axonal growth by miR-142-3p inhibitor and altered the expression level of the corresponding proteins. Thus, a substitution therapy using Sorafenib that downregulates the miR-142-3p expression for SNCI was investigated. The results showed the effect of Sorafenib was similar to that of miR-142-3p inhibitor and SNCI on both axon growth in vitro and sensory conduction function recovery in vivo. In conclusion, miR-142-3p acts a pivotal role in SNCI promoting the repair of dorsal column injury. Sorafenib mimics the treatment effect of SNCI via downregulation of miR-142-3p, subsequently, promoting sensory conduction function recovery post dorsal column injury.

MiR-20a Plays a Key Regulatory Role in the Repair of Spinal Cord Dorsal Column Lesion via PDZ-RhoGEF/RhoA/GAP43 Axis in Rat

Spinal cord injury (SCI) causes sensory dysfunctions such as paresthesia, dysesthesia, and chronic neuropathic pain. MiR-20a facilitates the axonal outgrowth of the cortical neurons. However, the role of miR-20a in the axonal outgrowth of primary sensory neurons and spinal cord dorsal column lesion (SDCL) is yet unknown. Therefore, the role of miR-20a post-SDCL was investigated in rat. The NF-200 immunofluorescence staining was applied to observe whether axonal outgrowth of dorsal root ganglion (DRG) neurons could be altered by miR-20a or PDZ-RhoGEF modulation in vitro. The expression of miR-20a was quantized with RT-PCR. Western blotting analyzed the expression of PDZ-RhoGEF/RhoA/GAP43 axis after miR-20a or PDZ-RhoGEF was modulated. The spinal cord sensory conduction function was assessed by somatosensory-evoked potentials and tape removal test. The results demonstrated that the expression of miR-20a decreased in a time-dependent manner post-SDCL. The regulation of miR-20a modulated the axonal growth and the expression of PDZ-RhoGEF/RhoA/GAP43 axis in vitro. The in vivo regulation of miR-20a altered the expression of miR-20a-PDZ-RhoGEF/RhoA/GAP43 axis and promoted the recovery of ascending sensory function post-SDCL. The results indicated that miR-20a/PDZ-RhoGEF/RhoA/GAP43 axis is associated with the pathophysiological process of SDCL. Thus, targeting the miR-20a/PDZ-RhoGEF /RhoA/GAP43 axis served as a novel strategy in promoting the sensory function recovery post-SCI.

Component analysis of somatosensory evoked potentials for identifying spinal cord injury location

This study aims to determine whether the time-frequency components (TFCs) of somatosensory evoked potentials (SEPs) can be used to identify the specific location of a compressive spinal cord injury using a classification technique. Waveforms of SEPs after compressive injuries at various locations (C4, C5 and C6) in rat spinal cords were decomposed into a series of TFCs using a high-resolution time-frequency analysis method. A classification method based on support vector machine (SVM) was applied to the distributions of these TFCs among different pathological locations. The difference among injury locations manifests itself in different categories of SEP TFCs. High-energy TFCs of normal-state SEPs have significantly higher power and frequency than those of injury-state SEPs. The location of C5 is characterized by a unique distribution pattern of middle-energy TFCs. The difference between C4 and C6 is evidenced by the distribution pattern of low-energy TFCs. The proposed classification method based on SEP TFCs offers a discrimination accuracy of 80.2%. In this study, meaningful information contained in various SEP components was investigated and used to propose a new application of SEPs for identification of the location of pathological changes in the cervical spinal cord.

Topography of Synchronization of Somatosensory Evoked Potentials Elicited by Stimulation of the Sciatic Nerve in Rat

Purpose: Traditionally, the topography of somatosensory evoked potentials (SEPs) is generated based on amplitude and latency. However, this operation focuses on the physical morphology and field potential-power, so it suffers from difficulties in performing identification in an objective manner. In this study, measurement of the synchronization of SEPs is proposed as a method to explore brain functional networks as well as the plasticity after peripheral nerve injury.

Method: SEPs elicited by unilateral sciatic nerve stimulation in twelve adult male Sprague-Dawley (SD) rats in the normal group were compared with SEPs evoked after unilateral sciatic nerve hemisection in four peripheral nerve injured SD rats. The characterization of synchronized networks from SEPs was conducted using equal-time correlation, correlation matrix analysis, and comparison to randomized surrogate data. Eigenvalues of the correlation matrix were used to identify the clusters of functionally synchronized neuronal activity, and the participation index (PI) was calculated to indicate the involvement of each channel in the cluster. The PI value at the knee point of the PI histogram was used as a threshold to demarcate the cortical boundary.

Results: Ten out of the twelve normal rats showed only one synchronized brain network. The remaining two normal rats showed one strong and one weak network. In the peripheral nerve injured group, only one synchronized brain network was found in each rat. In the normal group, all network shapes appear regular and the network is largely contained in the posterior cortex. In the injured group, the network shapes appear irregular, the network extends anteriorly and posteriorly, and the network area is significantly larger. There are considerable individual variations in the shape and location of the network after peripheral nerve injury.

Conclusion: The proposed method can detect functional brain networks. Compared to the results of the traditional SEP-morphology-based analysis method, the synchronized functional network area is much larger. Furthermore, the proposed method can also characterize the rapid cortical plasticity after a peripheral nerve is acutely injured.