Investigation of human cervical and upper thoracic spinal cord motion: implications for imaging spinal cord structure and function

DeepRetroMoCo: deep neural network-based retrospective motion correction algorithm for spinal cord functional MRI

Background and purpose

There are distinct challenges in the preprocessing of spinal cord fMRI data, particularly concerning the mitigation of voluntary or involuntary movement artifacts during image acquisition. Despite the notable progress in data processing techniques for movement detection and correction, applying motion correction algorithms developed for the brain cortex to the brainstem and spinal cord remains a challenging endeavor.

Methods

In this study, we employed a deep learning-based convolutional neural network (CNN) named DeepRetroMoCo, trained using an unsupervised learning algorithm. Our goal was to detect and rectify motion artifacts in axial T2*-weighted spinal cord data. The training dataset consisted of spinal cord fMRI data from 27 participants, comprising 135 runs for training and 81 runs for testing.

Results

To evaluate the efficacy of DeepRetroMoCo, we compared its performance against the sct_fmri_moco method implemented in the spinal cord toolbox. We assessed the motion-corrected images using two metrics: the average temporal signal-to-noise ratio (tSNR) and Delta Variation Signal (DVARS) for both raw and motion-corrected data. Notably, the average tSNR in the cervical cord was significantly higher when DeepRetroMoCo was utilized for motion correction, compared to the sct_fmri_moco method. Additionally, the average DVARS values were lower in images corrected by DeepRetroMoCo, indicating a superior reduction in motion artifacts. Moreover, DeepRetroMoCo exhibited a significantly shorter processing time compared to sct_fmri_moco.

Conclusion

Our findings strongly support the notion that DeepRetroMoCo represents a substantial improvement in motion correction procedures for fMRI data acquired from the cervical spinal cord. This novel deep learning-based approach showcases enhanced performance, offering a promising solution to address the challenges posed by motion artifacts in spinal cord fMRI data.

Mechanisms of Action of Dorsal Root Ganglion Stimulation

The dorsal root ganglion (DRG) serves as a pivotal site for managing chronic pain through dorsal root ganglion stimulation (DRG-S). In recent years, the DRG-S has emerged as an attractive modality in the armamentarium of neuromodulation therapy due to its accessibility and efficacy in alleviating chronic pain refractory to conventional treatments. Despite its therapeutic advantages, the precise mechanisms underlying DRG-S-induced analgesia remain elusive, attributed in part to the diverse sensory neuron population within the DRG and its modulation of both peripheral and central sensory processing pathways. Emerging evidence suggests that DRG-S may alleviate pain by several mechanisms, including the reduction of nociceptive signals at the T-junction of sensory neurons, modulation of pain gating pathways within the dorsal horn, and regulation of neuronal excitability within the DRG itself. However, elucidating the full extent of DRG-S mechanisms necessitates further exploration, particularly regarding its supraspinal effects and its interactions with cognitive and affective networks. Understanding these mechanisms is crucial for optimizing neurostimulation technologies and improving clinical outcomes of DRG-S for chronic pain management. This review provides a comprehensive overview of the DRG anatomy, mechanisms of action of the DRG-S, and its significance in neuromodulation therapy for chronic pain.

Closing the loop and raising the bar: Automated control systems in neuromodulation

INTRODUCTION

Neuromodulation has emerged as a promising therapy for the management of chronic pain, movement disorders, and other neurological conditions. Spinal cord stimulation (SCS) is a widely used form of neuromodulation that involves the delivery of electrical impulses to the spinal cord to modulate the transmission of pain signals to the brain. In recent years, there has been increasing interest in the use of automation systems to improve the efficacy and safety of SCS. This narrative review summarizes the status of Food and Drug Administration-approved autonomous neuromodulation devices including closed loop, feedforward, and feedback systems. The review discusses the advantages and disadvantages of each system and focuses specifically on the use of these systems for SCS. It is important for clinicians to understand the expanding role of automation in neuromodulation in order to select appropriate therapies founded on automation systems to the specific needs of the patient and the underlying condition.

CONCLUSION

The review also provides insights into the current state of the art in neuromodulation automation systems and discusses potential future directions for research in this field.

Increased cranio-caudal spinal cord oscillations are the cardinal pathophysiological change in degenerative cervical myelopathy

Introduction

Degenerative cervical myelopathy (DCM) is the most common cause of non-traumatic incomplete spinal cord injury, but its pathophysiology is poorly understood. As spinal cord compression observed in standard MRI often fails to explain a patient's status, new diagnostic techniques to assess DCM are one of the research priorities. Minor cardiac-related cranio-caudal oscillations of the cervical spinal cord are observed by phase-contrast MRI (PC-MRI) in healthy controls (HCs), while they become pathologically increased in patients suffering from degenerative cervical myelopathy. Whether transversal oscillations (i.e., anterior-posterior and right-left) also change in DCM patients is not known.

Methods

We assessed spinal cord motion simultaneously in all three spatial directions (i.e., cranio-caudal, anterior-posterior, and right-left) using sagittal PC-MRI and compared physiological oscillations in 18 HCs to pathological changes in 72 DCM patients with spinal canal stenosis. The parameter of interest was the amplitude of the velocity signal (i.e., maximum positive to maximum negative peak) during the cardiac cycle.

Results

Most patients suffered from mild DCM (mJOA score 16 (14-18) points), and the majority (68.1%) presented with multisegmental stenosis. The spinal canal was considerably constricted in DCM patients in all segments compared to HCs. Under physiological conditions in HCs, the cervical spinal cord oscillates in the cranio-caudal and anterior-posterior directions, while right-left motion was marginal [e.g., segment C5 amplitudes: cranio-caudal: 0.40 (0.27-0.48) cm/s; anterior-posterior: 0.18 (0.16-0.29) cm/s; right-left: 0.10 (0.08-0.13) cm/s]. Compared to HCs, DCM patients presented with considerably increased cranio-caudal oscillations due to the cardinal pathophysiologic change in non-stenotic [e.g., segment C5 amplitudes: 0.79 (0.49-1.32) cm/s] and stenotic segments [.g., segment C5 amplitudes: 0.99 (0.69-1.42) cm/s]). In contrast, right-left [e.g., segment C5 amplitudes: non-stenotic segment: 0.20 (0.13-0.32) cm/s; stenotic segment: 0.11 (0.09-0.18) cm/s] and anterior-posterior oscillations [e.g., segment C5 amplitudes: non-stenotic segment: 0.26 (0.15-0.45) cm/s; stenotic segment: 0.11 (0.09-0.18) cm/s] remained on low magnitudes comparable to HCs.

Conclusion

Increased cranio-caudal oscillations of the cervical cord are the cardinal pathophysiologic change and can be quantified using PC-MRI in DCM patients. This study addresses spinal cord oscillations as a relevant biomarker reflecting dynamic mechanical cord stress in DCM patients, potentially contributing to a loss of function.

A systematic review of repeatability and reproducibility studies of diffusion tensor imaging of cervical spinal cord

OBJECTIVES

Diffusion tensor imaging (DTI) techniques are being studied as a possible diagnostic and predictive tool for the evaluation of cervical spinal cord disease. This systematic review aims to evaluate the previous DTI studies that specifically investigated the repeatability and reproducibility of DTI in the cervical spinal cord.

METHODS AND MATERIALS

A search in the PubMed, Scopus, Web of Science and Ovid electronic databases was conducted for articles published between January 1990 and February 2022 that related to the repeatability and reproducibility of DTI in evaluating the cervical spinal cord using one of the following measurements: the intraclass correlation coefficient (ICC) and/or the coefficient of variation (CV), and/or Bland-Altman (BA) differences analysis methods. DTI studies that presented full statistical analysis of repeatability and/or reproducibility tests of the cervical spinal cord in peer-reviewed full-text publications published in journals were included. Articles that included at least one of the keywords within the titles or abstracts were identified. Additional full-text papers were found by searching the citations and reference lists of related articles. This review has followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidance. Risk of bias was evaluated with 13 criteria weighted toward methodological quality of reported studies using the QuADS assessment criteria. This assessment only included full-text articles written in English.

RESULTS

A total of 11 studies were included and assessed for different characteristics, including sample size,(3-34) re-test time interval (<1 h to >3 months), test-retest reproducibility scores and acquisition method. Six studies used ICC which ranged from poor (ICC<0.37) to excellent reproducibility (ICC 0.91-0.99). Four studies reported an overall CV lower than 40% for all DTI metrics. Three studies reported the Bland-Altman (BA) differences and reported a minimum percentage showing no strong differences between repeated measurements. Quantitative analysis was not undertaken due to heterogeneity of methods. Repeatability and reproducibility measures were generally found to be good.

CONCLUSION

This study revealed that the application of DTI and its related measures in a clinical setting in the assessment of cervical spinal cord changes is feasible and reproducible. However, cervical spinal cord DTI suffers from some existing limitations that prevent it from being routinely used in research and clinical settings.

ADVANCES IN KNOWLEDGE

DTI with its parametric maps provide broad evaluation of the tissue structure of axonal white matter and are being studied as a possible diagnostic and predictive tool for the assessment of cervical spinal cord (CSC) disease.

Utility of MRI in Quantifying Tissue Injury in Cervical Spondylotic Myelopathy

Cervical spondylotic myelopathy (CSM) is a progressive disease that worsens over time if untreated. However, the rate of progression can vary among individuals and may be influenced by various factors, such as the age of the patients, underlying conditions, and the severity and location of the spinal cord compression. Early diagnosis and prompt treatment can help slow the progression of CSM and improve symptoms. There has been an increased use of magnetic resonance imaging (MRI) methods in diagnosing and managing CSM. MRI methods provide detailed images and quantitative structural and functional data of the cervical spinal cord and brain, allowing for an accurate evaluation of the extent and location of tissue injury. This review aims to provide an understanding of the use of MRI methods in interrogating functional and structural changes in the central nervous system in CSM. Further, we identified several challenges hindering the clinical utility of these neuroimaging methods.

Variation in Diffusion Tensor Imaging Parameters in the Cervical and Thoracic Spinal Cord (C1-C5 and C6-T2) Segments of Normal Beagle Dogs

This study aimed to determine the characteristics and reference values of each vertebra in the cervicothoracic region by performing diffusion tensor imaging (DTI) scans and analyzing DTI parameters in normal Beagle dogs. In five adult Beagles under anesthetic maintenance, DTI was performed using a 1.5-T magnetic resonance imaging (MRI) scanner. Axial DTI was performed using three overlapping slabs to cover the cervical and thoracic spinal cords. After post-processing, DTI parameters were calculated along the entire spinal cord. Among DTI parameters, fractional anisotropy, relative anisotropy, and axonal diffusivity significantly decreased in the caudal direction. However, the apparent diffusion coefficient, radial diffusivity, and mean diffusivity values were not significantly correlated with vertebral levels. We provide evidence for the existence of segment-dependent DTI parameters in the canine cervical spinal cord. Therefore, comparisons of DTI parameters between lesions at different vertebral levels should be avoided unless normative data are available. Furthermore, the DTI data obtained in this study may contribute to the development of a clinical reference for spinal cord evaluation in dogs using DTI parameters.

Effects of variability in manually contoured spinal cord masks on fMRI co-registration and interpretation

Functional magnetic resonance imaging (fMRI) of the human spinal cord (SC) is a unique non-invasive method for characterizing neurovascular responses to stimuli. Group-analysis of SC fMRI data involves co-registration of subject-level data to standard space, which requires manual masking of the cord and may result in bias of group-level SC fMRI results. To test this, we examined variability in SC masks drawn in fMRI data from 21 healthy participants from a completed study mapping responses to sensory stimuli of the C7 dermatome. Masks were drawn on temporal mean functional image by eight raters with varying levels of neuroimaging experience, and the rater from the original study acted as a reference. Spatial agreement between rater and reference masks was measured using the Dice Similarity Coefficient, and the influence of rater and dataset was examined using ANOVA. Each rater's masks were used to register functional data to the PAM50 template. Gray matter-white matter signal contrast of registered functional data was used to evaluate the spatial normalization accuracy across raters. Subject- and group-level analyses of activation during left- and right-sided sensory stimuli were performed for each rater's co-registered data. Agreement with the reference SC mask was associated with both rater (F(7, 140) = 32.12, P < 2 × 10-16, η2 = 0.29) and dataset (F(20, 140) = 20.58, P < 2 × 10-16, η2 = 0.53). Dataset variations may reflect image quality metrics: the ratio between the signal intensity of spinal cord voxels and surrounding cerebrospinal fluid was correlated with DSC results (p < 0.001). As predicted, variability in the manually-drawn masks influenced spatial normalization, and GM:WM contrast in the registered data showed significant effects of rater and dataset (rater: F(8, 160) = 23.57, P < 2 × 10-16, η2 = 0.24; dataset: F(20, 160) = 22.00, P < 2 × 10-16, η2 = 0.56). Registration differences propagated into subject-level activation maps which showed rater-dependent agreement with the reference. Although group-level activation maps differed between raters, no systematic bias was identified. Increasing consistency in manual contouring of spinal cord fMRI data improved co-registration and inter-rater agreement in activation mapping, however our results suggest that improvements in image acquisition and post-processing are also critical to address.

Comparison of axial and sagittal spinal cord motion measurements in degenerative cervical myelopathy

BACKGROUND AND PURPOSE

The timing of decision-making for a surgical intervention in patients with mild degenerative cervical myelopathy (DCM) is challenging. Spinal cord motion phase contrast MRI (PC-MRI) measurements can reveal the extent of dynamic mechanical strain on the spinal cord to potentially identify high-risk patients. This study aims to determine the comparability of axial and sagittal PC-MRI measurements of spinal cord motion with the prospect of improving the clinical workup.

METHODS

Sixty-four DCM patients underwent a PC-MRI scan assessing spinal cord motion. The agreement of axial and sagittal measurements was determined by means of intraclass correlation coefficients (ICCs) and Bland-Altman analyses.

RESULTS

The comparability of axial and sagittal PC-MRI measurements was good to excellent at all cervical levels (ICCs motion amplitude: .810-.940; p &lt; .001). Significant differences between axial and sagittal amplitude values could be found at segments C3 and C4, while its magnitude was low (C3: 0.07 ± 0.19 cm/second; C4: -0.12 ± 0.30 cm/second). Bland-Altman analysis showed a good agreement between axial and sagittal PC-MRI scans (coefficients of repeatability: minimum -0.23 cm/second at C2; maximum -0.58 cm/second at C4). Subgroup analysis regarding anatomic conditions (stenotic vs. nonstenotic segments) and different velocity encoding (2 vs. 3 cm/second) showed comparable results.

CONCLUSIONS

This study demonstrates good comparability between axial and sagittal spinal cord motion measurements in DCM patients. To this end, axial and sagittal PC-MRI are both accurate and sensitive in detecting pathologic cord motion. Therefore, such measures could identify high-risk patients and improve clinical decision-making (ie, timing of decompression).

Morphology of Herniated Disc as a Predictor for Outcomes of Posterior Percutaneous Full-endoscopic Cervical Discectomy in Treating Cervical Spondylotic Radiculopathy

Objective

To quantitively characterize the morphology of cervical disc herniation (CDH) causing cervical spondylotic radiculopathy (CSR) and investigate whether the morphological features of CDH are associated with clinical outcomes in CSR patients treated by posterior percutaneous full‐endoscopic cervical discectomy (PPECD).

Methods

This is a single‐center retrospective study. Eighty‐seven PPECD‐treated patients meeting the inclusion criteria were included between May 2017 and May 2019. Based on preoperative T2‐weighted magnetic resonance imaging (MRI), we designed and measured six morphological parameters of CDH for all patients to reflect its relative position to cervical spinal cord and protruding degree: DC‐SC distance from the center of disc (DC) and the center of spinal cord (SC); DC‐DP distance from the center of cervical disc (DC) to the peak of herniation (DP); internal diameter of the disc; axial length of CDH; central angle of CDH formed by central axes of CDH and spinal cord; the modified index of CDH. We recorded general information, neck disability index (NDI) scores, visual analog scale (VAS) scores of neck and arm of all patients preoperatively and postoperatively at 1‐year follow‐up. The association of preoperative general variables and morphological parameters with clinical outcomes were explored by utilizing logistic regression and receiver operating characteristic curve (ROC) analysis.

Results

The preoperative neck‐VAS, arm‐VAS, and NDI were significantly decreased after PPECD and remained at a low value at follow‐up. In regards to the morphological parameters of CDH, the mean value of DC‐SC distance, DC‐DP distance, internal diameter of the disc, axial length of CDH, central angle of CDH, and modified index of CDH were 1.61 ± 0.30 cm, 1.66 ± 0.32cm, 1.04 ± 0.21 cm, 0.63 ± 0.19cm, 39.38° ± 11.94°, and 0.39 ± 0.24, respectively. For patients grouped by difference in the recovery rate of NDI and arm‐VAS (excellent improved group, EI; and limited improved group, LI), there were no differences in the age, gender, surgical segments, and morphological parameters, except for the central angle of CDH. According to binary logistic regression analysis, only the preoperative central angle of CDH was significantly associated with postoperative NDI recovery (odds ratio: 0.873; 95% confidence interval: 0.819–0.931, P = 0.002). ROC analysis showed the optimal cut‐off value of the central angle of CDH for predicting the postoperative improvement of functional outcomes is 33.788°.

Conclusion

Preoperative morphology of CDH is related to the outcomes of CSR patients after PPECD. Patients with a large central angle of CDH (>33.788°) have more likelihood of ameliorating neurological symptoms of CSR. There is the potential to select the central angle of CDH as a predictor for outcomes of PPECD in treating CSR.

Spinal Cord Motion in Degenerative Cervical Myelopathy: The Level of the Stenotic Segment and Gender Cause Altered Pathodynamics

In degenerative cervical myelopathy (DCM), focally increased spinal cord motion has been observed for C5/C6, but whether stenoses at other cervical segments lead to similar pathodynamics and how severity of stenosis, age, and gender affect them is still unclear. We report a prospective matched-pair controlled trial on 65 DCM patients. A high-resolution 3D T2 sampling perfection with application-optimized contrasts using different flip angle evolution (SPACE) and a phase-contrast magnetic resonance imaging (MRI) sequence were performed and automatically segmented. Anatomical and spinal cord motion data were assessed per segment from C2/C3 to C7/T1. Spinal cord motion was focally increased at a level of stenosis among patients with stenosis at C4/C5 (n = 14), C5/C6 (n = 33), and C6/C7 (n = 10) (p < 0.033). Patients with stenosis at C2/C3 (n = 2) and C3/C4 (n = 6) presented a similar pattern, not reaching significance. Gender was a significant predictor of higher spinal cord dynamics among men with stenosis at C5/C6 (p = 0.048) and C6/C7 (p = 0.033). Age and severity of stenosis did not relate to spinal cord motion. Thus, the data demonstrates focally increased spinal cord motion depending on the specific level of stenosis. Gender-related effects lead to dynamic alterations among men with stenosis at C5/C6 and C6/C7. The missing relation of motion to severity of stenosis underlines a possible additive diagnostic value of spinal cord motion analysis in DCM.

Spinal cord fMRI with MB-SWIFT for assessing epidural spinal cord stimulation in rats

Purpose

Electrical epidural spinal cord stimulation (SCS) is used as a treatment for chronic pain as well as to partially restore motor function after a spinal cord injury. Monitoring the spinal cord activity during SCS with fMRI could provide important and objective measures of integrative responses to treatment. Unfortunately, spinal cord fMRI is severely challenged by motion and susceptibility artifacts induced by the implanted electrode and bones. This pilot study introduces multi‐band sweep imaging with Fourier transformation (MB‐SWIFT) technique for spinal cord fMRI during SCS in rats. Given the close to zero acquisition delay and high bandwidth in 3 dimensions, MB‐SWIFT is demonstrated to be highly tolerant to motion and susceptibility‐induced artifacts and thus holds promise for fMRI during SCS.

Methods

MB‐SWIFT with 0.78 × 0.78 × 1.50 mm³ spatial resolution and 3‐s temporal resolution was used at 9.4 Tesla in rats undergoing epidural SCS at different frequencies. Its performance was compared with spin echo EPI. The origin of the functional contrast was also explored using suppression bands.

Results

MB‐SWIFT was tolerant to electrode‐induced artifacts and respiratory motion, leading to substantially higher fMRI sensitivity than spin echo fMRI. Clear stimulation frequency‐dependent responses to SCS were detected in the rat spinal cord close to the stimulation site. The origin of MB‐SWIFT fMRI signals was consistent with dominant inflow effects.

Conclusion

fMRI of the rat spinal cord during SCS can be consistently achieved with MB‐SWIFT, thus providing a valuable experimental framework for assessing the effects of SCS on the central nervous system.

Focal cervical spinal stenosis causes mechanical strain on the entire cervical spinal cord tissue - A prospective controlled, matched-pair analysis based on phase-contrast MRI

BACKGROUND

Focally increased spinal cord motion at the level of cervical spinal stenosis has been revealed by phase-contrast MRI (PC-MRI).

OBJECTIVE

To investigate spinal cord motion among patients suffering of degenerative cervical myelopathy (DCM) across the entire cervical spine applying automated segmentation and standardized PC-MRI post-processing protocols.

METHODS

Prospective, matched-pair controlled trial on 29 patients with stenosis at C5/C6. MRI-protocol covering all cervical segments: 3D T2-SPACE, prospectively ECG-triggered sagittal PC-MRI. Segmentation by trained 3D hierarchical deep convolutional neural network and data processing were conducted via in-house software pipeline. Parameters per segment: maximum velocity, peak-to-peak (PTP)-amplitude, total displacement, PTP-amplitude^HB (PTP-amplitude per duration of heartbeat), and, for characterization of intraindividual alterations, the PTP-amplitude index between the cervical segments C3/C4-C7/T1 and C2/C3.

RESULTS

Spinal cord motion was increased at C4/C5, C5/C6 and C6/C7 among patients (all parameters, p < 0.001-0.025). The PTP-amplitude index revealed an increase from C3/C4 to C4/C5 (p = 0.002), C4/C5 to C5/C6 (p = 0.037) and a decrease from C5/C6 to C6/C7 and C6/C7 to C7/T1 (p < 0.001, each). This implied an up-building stretch on spinal cord tissue cranial and a mechanical compression caudal of the stenotic level. Furthermore, significant far range effects across the entire cervical spinal cord were observed (e.g. PTP-amplitude C2/C3 vs. C6/C7, p = 0.026) in contrast to controls (p = 1.00).

CONCLUSION

This study revealed the nature and extends of mechanical stress on the entire cervical spinal cord tissue due to focal stenosis. These pathophysiological alterations of spinal cord motion can be expected to be clinically relevant.

Effect of distortion corrections on the tractography quality in spinal cord diffusion-weighted imaging

PURPOSE

To assess the impact of a different distortion correction (DC) method and patient geometry (sagittal balance) on the quality of spinal cord tractography rendering according to different tractography approaches.

METHODS

Forty-four adults free of spinal cord diseases underwent cervical diffusion-weighted imaging. The phase-encoding direction was head→foot. Sequence with opposed polarities (foot→head) was acquired to perform DC. Eddy-current, motion effects, and susceptibility artifact correction methods were used for DC, and two deterministic and one probabilistic tractography approaches were evaluated using MRtrix and DSI Studio tractography software. Fiber length and number of fibers were extracted to evaluate the quality of the tractography rendering. For each subject, cervical lordosis was measured to assess patient geometry. The angle between the main direction of the spinal cord and the orientation of the acquisition box were computed at each spine level to assess acquisition geometry and define an angle threshold for which a tractography of good quality is no longer possible.

RESULTS

There was a significant improvement in tractography quality after performing DC with susceptibility artifact correction using a deterministic approach based on tensor. Before DC, the angle threshold was defined at C6 (15.2°) compared with C7 (21.9°) after corrections, demonstrating the importance of spinal cord angulation for DC.

CONCLUSION

The impact of DC on tractography quality is greatly impacted by acquisition geometry. To obtain a good-quality tractography, we propose as a future perspective to adapt the acquisition geometry to that of the patient by automatically adjusting the acquisition box.

The Dancing Cord: Inherent Spinal Cord Motion and Its Effect on Cord Dose in Spine Stereotactic Body Radiation Therapy

BACKGROUND

Spinal cord dose limits are critically important for the safe practice of spine stereotactic body radiotherapy (SBRT). However, the effect of inherent spinal cord motion on cord dose in SBRT is unknown.

OBJECTIVE

To assess the effects of cord motion on spinal cord dose in SBRT.

METHODS

Dynamic balanced fast field echo (BFFE) magnetic resonance imaging (MRI) was obtained in 21 spine metastasis patients treated with SBRT. Planning computed tomography (CT), conventional static T2-weighted MRI, BFFE MRI, and dose planning data were coregistered. Spinal cord from the dynamic BFFE images (corddyn) was compared with the T2-weighted MRI (cordstat) to analyze motion of corddyn beyond the cordstat (Dice coefficient, Jaccard index), and beyond cordstat with added planning organ at risk volume (PRV) margins. Cord dose was compared between cordstat, and corddyn (Wilcoxon signed-rank test).

RESULTS

Dice coefficient (0.70-0.95, median 0.87) and Jaccard index (0.54-0.90, median 0.77) demonstrated motion of corddyn beyond cordstat. In 62% of the patients (13/21), the dose to corddyn exceeded that of cordstat by 0.6% to 13.8% (median 4.3%). The corddyn spatially excursed outside the 1-mm PRV margin of cordstat in 9 patients (43%); among these dose to corddyn exceeded dose to cordstat >+ 1-mm PRV margin in 78% of the patients (7/9). Corddyn did not excurse outside the 1.5-mm or 2-mm PRV cord cordstat margin.

CONCLUSION

Spinal cord motion may contribute to increases in radiation dose to the cord from SBRT for spine metastasis. A PRV margin of at least 1.5 to 2 mm surrounding the cord should be strongly considered to account for inherent spinal cord motion.

Feasibility of human spinal cord perfusion mapping using dynamic susceptibility contrast imaging at 7T: Preliminary results and identified guidelines

PURPOSE

To explore the feasibility of dynamic susceptibility contrast MRI at 7 Tesla for human spinal cord perfusion mapping and fill the gap between brain and spinal cord perfusion mapping techniques.

METHODS

Acquisition protocols for high-resolution single shot EPI in the spinal cord were optimized for both spin-echo and gradient-echo preparations, including cardiac gating, acquisition times and breathing cycle recording. Breathing-induced MRI signal fluctuations were investigated in healthy volunteers. A specific image- and signal-processing pipeline was implemented to address them. Dynamic susceptibility contrast was then evaluated in 3 healthy volunteers and 5 patients. Bolus depiction on slice-wise signal within cord was investigated, and maps of relative perfusion indices were computed.

RESULTS

Signal fluctuations were increased by 1.9 and 2.3 in free-breathing compared to apnea with spin-echo and gradient-echo, respectively. The ratio between signal fluctuations and bolus peak in healthy volunteers was 5.0% for spin-echo and 3.8% for gradient-echo, allowing clear depiction of the bolus on every slice and yielding relative blood flow and volume maps exhibiting the expected higher perfusion of gray matter. However, signal fluctuations in patients were increased by 4 in average (using spin-echo), compromising the depiction of the bolus in slice-wise signal. Moreover, 3 of 18 slices had to be discarded because of fat-aliasing artifacts.

CONCLUSION

Dynamic susceptibility contrast MRI at 7 Tesla showed great potential for spinal cord perfusion mapping with a reliability never achieved thus far for single subject and single slice measurements. Signal stability needs to be improved in acquisition conditions associated with patients; guidelines to achieve that have been identified and shared.

Ultrahigh-resolution quantitative spinal cord MRI at 9.4T

PURPOSE

To present the results of the first human spinal cord in vivo MRI scans at 9.4T.

METHODS

A human brain coil was used to image the human spinal cord at 9.4T. All anatomical images were acquired with a T₂ *-weighted gradient-echo sequence. A comparison of the influence of four different B₀ shimming routines on the image quality was performed. Intrinsic signal-to-noise-ratio maps were determined using a pseudo-multiple replica approach. Measurements with different echo times were compared and processed to one multiecho data image combination image. Based on the multiecho acquisitions, T₂ *-relaxation time maps were calculated. Algorithmic spinal cord detection and gray matter/white matter segmentation were tested.

RESULTS

An echo time between 9 and 13.8 ms compromised best between gray matter/white matter contrast and image quality. A maximum in-plane resolution of 0.15 × 0.15 mm² was achieved for anatomical images. These images offered excellent image quality and made small structures of the spinal cord visible. The scanner vendor implemented B₀ shimming routine performed best during this work. Intrinsic signal-to-noise-ratio values of between 6600 and 8060 at the upper cervical spinal cord were achieved. Detection and segmentation worked reliably. An average T₂ *-time of 24.88 ms ± 6.68 ms for gray matter and 19.37 ms ± 8.66 ms for white matter was calculated.

CONCLUSION

The proposed human brain coil can be used to image the spinal cord. The maximum in-plane resolution in this work was higher compared with the 7T results from the literature. The 9.4T acquisitions made the small structures of the spinal cord clearly visible.

Unique Characteristics of the Dorsal Root Ganglion as a Target for Neuromodulation

Objective

The dorsal root ganglion (DRG) is a novel target for neuromodulation, and DRG stimulation is proving to be a viable option in the treatment of chronic intractable neuropathic pain. Although the overall principle of conventional spinal cord stimulation (SCS) and DRG stimulation—in which an electric field is applied to a neural target with the intent of affecting neural pathways to decrease pain perception—is similar, there are significant differences in the anatomy and physiology of the DRG that make it an ideal target for neuromodulation and may account for the superior outcomes observed in the treatment of certain chronic neuropathic pain states. This review highlights the anatomy of the DRG, its function in maintaining homeostasis and its role in neuropathic pain, and the unique value of DRG as a target in neuromodulation for pain.

Methods

A narrative literature review was performed.

Results

Overall, the DRG is a critical structure in sensory transduction and modulation, including pain transmission and the maintenance of persistent neuropathic pain states. Unique characteristics including selective somatic organization, specialized membrane characteristics, and accessible and consistent location make the DRG an ideal target for neuromodulation. Because DRG stimulation directly recruits the somata of primary sensory neurons and harnesses the filtering capacity of the pseudounipolar neural architecture, it is differentiated from SCS, peripheral nerve stimulation, and other neuromodulation options.

Conclusions

There are several advantages to targeting the DRG, including lower energy usage, more focused and posture-independent stimulation, reduced paresthesia, and improved clinical outcomes.

Spinal cord stimulation programming: a crash course

The aim of this comprehensive review is to provide an instructional guide for providers regarding the parameters and programming of spinal cord stimulation (SCS) devices. Knowing these fundamentals will aid in providing superior pain relief to patients. SCS has four programmable parameters: contact (electrode) selection, amplitude, pulse width, and frequency. Each parameter needs to be accounted for when assessing which program works for which patient. Traditional open-loop systems allow for different "programs," or combinations of these four parameters, to be pre-set by the provider and medical device representative. These allow for flexibility in the type of stimulation delivered to the patient depending on activity. Patients are also given control over programs and changing the amplitudes of these programs. However, some open-loop systems place the burden of toggling between programs to manage pain control on patients, though this tends to be less in subparesthesia programs. Newer closed-loop systems make it possible for stimulation settings to automatically adjust in response to accelerometry and evoked compound action potential feedback, and therefore have the potential to streamline the patient experience. This article provides practitioners with the basic knowledge of SCS parameters and programming systems. Understanding their use is essential to providing optimal pain relief to patients.

Assessment of spinal cord motion as a new diagnostic MRI-parameter in cervical spinal canal stenosis: study protocol on a prospective longitudinal trial

Background

Increased spinal cord motion has been proven to be a relevant finding within spinal canal stenosis disclosed by phase-contrast MRI (PC-MRI). Adapted PC-MRI is a suitable and reliable method within the well deliberated setting. As the decision between conservative and operative treatment can be challenging in some cases, further diagnostic marker would facilitate the diagnostic process. We hypothesize that increased spinal cord motion will correlate to clinical course and functional impairment and will contribute as a new diagnostic marker.

Methods

A monocentric, prospective longitudinal observational trial on cervical spinal canal stenosis will be conducted at the University Medical Center Freiburg. Patients (n = 130) with relevant cervical spinal canal stenosis, being defined by at least contact to the spinal cord, will be included. Also, we will examine a control group of healthy volunteers (n = 20) as proof-of-principle. We will observe two openly assigned branches of participants undergoing conservative and surgical decompressive treatment (based on current German Guidelines) over a time course of 12 month, including a total of 4 visits. We will conduct a broad assessment of clinical parameters, standard scores and gradings, electrophysiological measurements, standard MRI, and adapted functional PC-MRI of spinal cord motion. Primary endpoint is the evaluation of an expected negative correlation of absolute spinal cord displacement to clinical impairment. Secondary endpoints are the evaluation of positive correlation of increased absolute spinal cord displacement to prolonged evoked potentials, prediction of clinical course by absolute spinal cord displacement, and demonstration of normalized spinal cord motion after decompressive surgery.

Discussion

With the use of adapted, non-invasive PC-MRI as a quantitative method for assessment of spinal cord motion, further objective diagnostic information can be gained, that might improve the therapeutic decision-making process. This study will offer the needed data in order to establish PC-MRI on spinal cord motion within the diagnostic work-up of patients suffering from spinal canal stenosis.

Trial registration

German Clinical Trials Register, ID: DRKS00012962, Register date 2018/01/17

Automatic Spinal Cord Gray Matter Quantification: A Novel Approach

BACKGROUND AND PURPOSE

Currently, accurate and reproducible spinal cord GM segmentation remains challenging and a noninvasive broadly accepted reference standard for spinal cord GM measurements is still a matter of ongoing discussion. Our aim was to assess the reproducibility and accuracy of cervical spinal cord GM and WM cross-sectional area measurements using averaged magnetization inversion recovery acquisitions images and a fully-automatic postprocessing segmentation algorithm.

MATERIALS AND METHODS

The cervical spinal cord of 24 healthy subjects (14 women; mean age, 40 ± 11 years) was scanned in a test-retest fashion on a 3T MR imaging system. Twelve axial averaged magnetization inversion recovery acquisitions slices were acquired over a 48-mm cord segment. GM and WM were both manually segmented by 2 experienced readers and compared with an automatic variational segmentation algorithm with a shape prior modified for 3D data with a slice similarity prior. Precision and accuracy of the automatic method were evaluated using coefficients of variation and Dice similarity coefficients.

RESULTS

The mean GM area was 17.20 ± 2.28 mm² and the mean WM area was 72.71 ± 7.55 mm² using the automatic method. Reproducibility was high for both methods, while being better for the automatic approach (all mean automatic coefficients of variation, ≤4.77%; all differences, P < .001). The accuracy of the automatic method compared with the manual reference standard was excellent (mean Dice similarity coefficients: 0.86 ± 0.04 for GM and 0.90 ± 0.03 for WM). The automatic approach demonstrated similar coefficients of variation between intra- and intersession reproducibility as well as among all acquired spinal cord slices.

CONCLUSIONS

Our novel approach including the averaged magnetization inversion recovery acquisitions sequence and a fully-automated postprocessing segmentation algorithm demonstrated an accurate and reproducible spinal cord GM and WM segmentation. This pipeline is promising for both the exploration of longitudinal structural GM changes and application in clinical settings in disorders affecting the spinal cord.

Quantification of DTI in the Pediatric Spinal Cord: Application to Clinical Evaluation in a Healthy Patient Population

BACKGROUND AND PURPOSE

The purpose of the study is to characterize diffusion tensor imaging indices in the developing spinal cord, evaluating differences based on age and cord region. Describing the progression of DTI indices in the pediatric cord increases our understanding of spinal cord development.

MATERIALS AND METHODS

A retrospective analysis was performed on DTI acquired in 121 pediatric patients (mean, 8.6 years; range, 0.3-18.0 years) at Monroe Carell Jr. Children's Hospital at Vanderbilt from 2017 to 2018. Diffusion-weighted images (15 directions; b = 750 s/mm²; slice thickness, 5 mm; in-plane resolution, 1.0 × 1.0 mm²) were acquired on a 3T scanner in the cervicothoracic and/or thoracolumbar cord. Manual whole-cord segmentation was performed. Images were masked and further segmented into cervical, upper thoracic, thoracolumbar, and conus regions. Analyses of covariance were performed for each DTI-derived index to investigate how age affects diffusion across cord regions, and 95% confidence intervals were calculated across age for each derived index and region. Post hoc testing was performed to analyze regional differences.

RESULTS

Analyses of covariance revealed significant correlations of age with axial diffusivity, mean diffusivity, and fractional anisotropy (all, P < .001). There were also significant differences among cord regions for axial diffusivity, radial diffusivity, mean diffusivity, and fractional anisotropy (all, P < .001).

CONCLUSIONS

This research demonstrates that diffusion evolves in the pediatric spinal cord during development, dependent on both cord region and the diffusion index of interest. Future research could investigate how diffusion may be affected by common pediatric spinal pathologies.

A method for correcting breathing-induced field fluctuations in T2*-weighted spinal cord imaging using a respiratory trace

PURPOSE

Spinal cord MRI at ultrahigh field is hampered by time-varying magnetic fields associated with the breathing cycle, giving rise to ghosting artifacts in multi-shot acquisitions. Here, we suggest a correction approach based on linking the signal from a respiratory bellows to field changes inside the spinal cord. The information is used to correct the data at the image reconstruction level.

METHODS

The correction was demonstrated in the context of multi-shot T2*-weighted imaging of the cervical spinal cord at 7T. A respiratory trace was acquired during a high-resolution multi-echo gradient-echo sequence, used for structural imaging and quantitative T2* mapping, and a multi-shot EPI time series, as would be suitable for fMRI. The coupling between the trace and the breathing-induced fields was determined by a short calibration scan in each individual. Images were reconstructed with and without trace-based correction.

RESULTS

In the multi-echo acquisition, breathing-induced fields caused severe ghosting in images with long TE, which led to a systematic underestimation of T2* in the spinal cord. The trace-based correction reduced the ghosting and increased the estimated T2* values. Breathing-related ghosting was also observed in the multi-shot EPI images. The correction largely removed the ghosting, thereby improving the temporal signal-to-noise ratio of the time series.

CONCLUSIONS

Trace-based retrospective correction of breathing-induced field variations can reduce ghosting and improve quantitative metrics in multi-shot structural and functional T2*-weighted imaging of the spinal cord. The method is straightforward to implement and does not rely on sequence modifications or additional hardware beyond a respiratory bellows.

Physiological motion of the optic chiasm and its impact on stereotactic radiosurgery dose

OBJECTIVE

Avoidance of radiation-induced optic neuropathy (RION) from stereotactic radiosurgery (SRS) requires precise anatomical localization; however, no prior studies have characterized the physiologic motion of the optic chiasm. We measured the extent of chiasm motion and its impact on SRS dose.

METHODS

In this cross-sectional study, serial MRI was performed in multiple planes in 11 human subjects without optic pathway abnormalities to determine chiasm motion across time. Subsequently, the measured displacement was applied to the hypothetical chiasm dose received in 11 patients treated with SRS to a perichiasmatic lesion.

RESULTS

On sagittal images, the average anteroposterior chiasm displacement was 0.51 mm [95% confidence interval (CI) 0.27 - 0.75 mm], and the average superior-inferior displacement was 0.48 mm (95% CI 0.22 - 0.74 mm). On coronal images, the average superior-inferior displacement was 0.42 mm (95% CI 0.13 - 0.71 mm), and the average lateral displacement was 0.75 mm (95% CI 0.42 - 1.08 mm). In 11 patients who underwent SRS to a perichiasmatic lesion, the average displacements increased the maximum chiasm dose (Dmax) by a mean of 14 % (range 6-23 %; p < 0.001).

CONCLUSION

Average motion of the optic chiasm was approximately 0.50-0.75 mm, which increased chiasm Dmax by a mean of 14%. In the occasional patient with higher-than-average chiasm motion in a region of steep dose gradient, the increase in chiasm Dmax and risk of RION could be even larger. Similarly, previously reported chiasm dose constraints may underestimate the true dose received during radiosurgery.

ADVANCES IN KNOWLEDGE

To limit the risk of RION, clinicians may consider adding a 0.50-0.75 mm expansion to the chiasm avoidance structure.

Segmental differences of cervical spinal cord motion: advancing from confounders to a diagnostic tool

Increased cranio-caudal spinal cord motion is associated with clinical impairment in degenerative cervical myelopathy. However, whether spinal cord motion holds potential as a neuroimaging biomarker requires further validation. Different confounders (i.e. subject characteristics, methodological problems such as phase drift, etc.) on spinal cord motion readouts have to be considered. Twenty-two healthy subjects underwent phase contrast MRI, a subset of subjects (N = 9) had repeated scans. Parameters of interest included amplitude of velocity signal, maximum cranial respectively maximum caudal velocity, displacement (=area under curve of the velocity signal). The cervical spinal cord showed pulse synchronic oscillatory motions with significant differences in all readouts across cervical segments, with a maximum at C5. The Inter-rater reliability was excellent for all readouts. The test-retest reliability was excellent for all parameters at C2 to C6, but not for maximum cranial velocity at C6 and all readouts at C7. Spinal cord motion was correlated with spinal canal size, heart rate and body size. This is the first study to propose a standardized MRI measurement of spinal cord motion for further clinical implementation based on satisfactory phase drift correction and excellent reliability. Understanding the influence of confounders (e.g. structural conditions of the spine) is essential for introducing cord motion into the diagnostic work up.

Spinal Cord Imaging in Amyotrophic Lateral Sclerosis: Historical Concepts-Novel Techniques

Amyotrophic lateral sclerosis (ALS) is the most common adult onset motor neuron disease with no effective disease modifying therapies at present. Spinal cord degeneration is a hallmark feature of ALS, highlighted in the earliest descriptions of the disease by Lockhart Clarke and Jean-Martin Charcot. The anterior horns and corticospinal tracts are invariably affected in ALS, but up to recently it has been notoriously challenging to detect and characterize spinal pathology in vivo. With recent technological advances, spinal imaging now offers unique opportunities to appraise lower motor neuron degeneration, sensory involvement, metabolic alterations, and interneuron pathology in ALS. Quantitative spinal imaging in ALS has now been used in cross-sectional and longitudinal study designs, applied to presymptomatic mutation carriers, and utilized in machine learning applications. Despite its enormous clinical and academic potential, a number of physiological, technological, and methodological challenges limit the routine use of computational spinal imaging in ALS. In this review, we provide a comprehensive overview of emerging spinal cord imaging methods and discuss their advantages, drawbacks, and biomarker potential in clinical applications, clinical trial settings, monitoring, and prognostic roles.

Diffusion tensor imaging shows mechanism-specific differences in injury pattern and progression in rat models of acute spinal cord injury

We investigate the ability of diffusion tensor imaging (DTI) to distinguish between three experimental rat models of spinal cord injury mechanism - contusion, dislocation, and distraction. Ex vivo DTI scans were performed on cord specimens that were preserved at different time points of the acute injury (3 hr, 24 hr, and 7 days post-injury) across all three injury mechanisms. White matter was classified as abnormal if their DTI metric was substantially different from regional values measured from a set of uninjured controls, thus allowing generation of binary "white matter damage maps" which categorizes each pixel in the DTI image as "normal" or "damaged". Damage classification was most robust using thresholds in the longitudinal diffusivity, which supports previous studies that show that longitudinal diffusivity is the most robust DTI metric in depicting damage in SCI. Furthermore, the spatial damage patterns from all subjects in the same group were consolidated into a "damage occurrence ratio map", which illustrates an average damage shape that characterizes the injury mechanism. Our analysis has yielded a dataset which highlights the differences in injury pattern due to the initial mode of mechanical injury. For example, contusion produced an initial injury that emanated radially outward from the central canal, with subsequent damage along the caudal corticospinal tract and rostral gracile fasciculus; dislocation injuries showed a high level of involvement in the lateral and ventral white matter which became less apparent by 7 days post-injury, and distraction injuries were found to be less focal and more distributed rostrocaudally. This work represents a first step in adopting the use of the primary injury mechanism as a clinical prognostic factor in SCI, which may help to inform the trialing of existing neuroprotective treatment candidates, the development of new therapies as well as personalize the management of SCI for the individual patient.

Ten Key Insights into the Use of Spinal Cord fMRI

A comprehensive review of the literature-to-date on functional magnetic resonance imaging (fMRI) of the spinal cord is presented. Spinal fMRI has been shown, over more than two decades of work, to be a reliable tool for detecting neural activity. We discuss 10 key points regarding the history, development, methods, and applications of spinal fMRI. Animal models have served a key purpose for the development of spinal fMRI protocols and for experimental spinal cord injury studies. Applications of spinal fMRI span from animal models across healthy and patient populations in humans using both task-based and resting-state paradigms. The literature also demonstrates clear trends in study design and acquisition methods, as the majority of studies follow a task-based, block design paradigm, and utilize variations of single-shot fast spin-echo imaging methods. We, therefore, discuss the similarities and differences of these to resting-state fMRI and gradient-echo EPI protocols. Although it is newly emerging, complex connectivity and network analysis is not only possible, but has also been shown to be reliable and reproducible in the spinal cord for both task-based and resting-state studies. Despite the technical challenges associated with spinal fMRI, this review identifies reliable solutions that have been developed to overcome these challenges.

Relevance of time-dependence for clinically viable diffusion imaging of the spinal cord

Purpose

Time‐dependence is a key feature of the diffusion‐weighted (DW) signal, knowledge of which informs biophysical modelling. Here, we study time‐dependence in the human spinal cord, as its axonal structure is specific and different from the brain.

Methods

We run Monte Carlo simulations using a synthetic model of spinal cord white matter (WM) (large axons), and of brain WM (smaller axons). Furthermore, we study clinically feasible multi‐shell DW scans of the cervical spinal cord (b = 0; b = 711 s mm⁻²; b = 2855 s mm⁻²), obtained using three diffusion times (Δ of 29, 52 and 76 ms) from three volunteers.

Results

Both intra‐/extra‐axonal perpendicular diffusivities and kurtosis excess show time‐dependence in our synthetic spinal cord model. This time‐dependence is reflected mostly in the intra‐axonal perpendicular DW signal, which also exhibits strong decay, unlike our brain model. Time‐dependence of the total DW signal appears detectable in the presence of noise in our synthetic spinal cord model, but not in the brain. In WM in vivo, we observe time‐dependent macroscopic and microscopic diffusivities and diffusion kurtosis, NODDI and two‐compartment SMT metrics. Accounting for large axon calibers improves fitting of multi‐compartment models to a minor extent.

Conclusions

Time‐dependence of clinically viable DW MRI metrics can be detected in vivo in spinal cord WM, thus providing new opportunities for the non‐invasive estimation of microstructural properties. The time‐dependence of the perpendicular DW signal may feature strong intra‐axonal contributions due to large spinal axon caliber. Hence, a popular model known as “stick” (zero‐radius cylinder) may be sub‐optimal to describe signals from the largest spinal axons.

Effect of cardiac-related translational motion in diffusion MRI of the spinal cord

Cardiac-related spinal cord motion affects diffusion-weighted (DWI) signal. The goal of this study was to further quantify the specific detrimental effect of cord translational motion on the DWI signal in order to make better informed decisions about the cost-benefit of cardiac gating. We designed an MRI-compatible phantom mimicking the spinal cord translational motion. Cardiac-gated DWI data were acquired by varying the trigger delay and the b-values. Evaluation of the effect of motion on the DWI signal was done by computing the apparent diffusion coefficient (ADC) along (z-direction) and orthogonal (y- and x-directions) to the phantom. The computed ADCs of the phantom moving along Z were similar for the three orthogonal diffusion-encoding directions, with an average value of 1.65·10^-9 , 1.66·10^-9 and 1.65·10^-9 m²/s along X, Y and Z respectively. DW phase images on the other hand showed the expected linear relationship with phantom velocity. Pure translational motion has minor effect on the diffusion-weighted magnitude signal. The sudden signal drop typically observed in in vivo spinal cord DWI is likely not caused by translational motion of the spinal cord, and possibly originates from non-rigid compression/stretching of the cord and/or from intra-voxel incoherent motion (IVIM).

Age related diffusion and tractography changes in typically developing pediatric cervical and thoracic spinal cord

Background and objective

Diffusion tensor imaging (DTI) and diffusion tensor tractography (DTT) are two techniques that can measure white matter integrity of the spinal cord. Recently, DTI indices have been shown to change with age. The purpose of this study is (a) to evaluate the maturational states of the entire pediatric spinal cord using DTI and DTT indices including fractional anisotropy (FA), mean diffusivity (MD), mean length of white matter fiber tracts and tract density and (b) to analyze the DTI and DTT parameters along the entire spinal cord as a function of spinal cord levels and age.

Method

A total of 23 typically developing (TD) pediatric subjects ranging in age from 6 to 16 years old (11.94 ± 3.26 (mean ± standard deviation), 13 females and 10 males) were recruited, and scanned using 3.0 T MR scanner. Reduced FOV diffusion tensor images were acquired axially in the same anatomical location prescribed for the T2-weighted images to cover the entire spinal cord (C1-mid L1 levels). To mitigate motion induced artifacts, diffusion directional images were aligned with the reference image (b0) using a rigid body registration algorithm performed by in-house software developed in Matlab (MathWorks, Natick, Massachusetts). Diffusion tensor maps (FA and MD) and streamline deterministic tractography were then generated from the motion corrected DTI dataset. DTI and DTT parameters were calculated by using ROIs drawn to encapsulate the whole cord along the entire spinal cord by an independent board certified neuroradiologist. These indices then were compared between two age groups (age group A = 6-11 years (n = 11) and age group B = 12-16 years (n = 12)) based on similar standards and age definitions used for reporting spinal cord injury in the pediatric population. Standard least squared linear regression based on a restricted maximum likelihood (REML) method was used to evaluate the relationship between age and DTI and DTT parameters.

Results

An increase in FA (group A = 0.42 ± 0.097, group B = 0.49 ± 0.116), white matter tract density (group A = 368.01 ± 236.88, group B = 440.13 ± 245.24) and mean length of fiber tracts (group A = 48.16 ± 20.48 mm, group B = 60.28 ± 23.87 mm) and a decrease in MD (group A = 1.06 ± 0.23 × 10-3 mm2/s, group B = 0.82 ± 0.24 × 10-3 mm2/s) were observed with age along the entire spinal cord. Statistically significant increases have been shown in FA (p = 0.004, R2 = 0.57), tract density (p = 0.0004, R2 = 0.58), mean length of fiber tracts (p < 0.001, R2 = 0.5) and a significant decrease has been shown in MD (p = 0.002, R2 = 0.59) between group A and group B. Also, it has been shown DTI and DTT parameters vary along the spinal cord as a function of intervertebral disk and mid-vertebral body level.

Conclusion

This study provides an initial understanding of age related changes of DTI values as well as DTT metrics of the spinal cord. The results show significant differences in DTI and DTT parameters which may result from decreasing water content, myelination of fiber tracts, and the thickening diameter of fiber tracts during the maturation process. Consequently, when quantitative DTI and DTT of the spinal cord is undertaken in the pediatric population an age and level matched normative dataset should be used to accurately interpret the quantitative results.

Reduced Field of View Diffusion Tensor Imaging and Fiber Tractography of the Pediatric Cervical and Thoracic Spinal Cord Injury

The aim of this study is to assess the utility and effectiveness of diffusion tensor imaging (DTI) and diffusion tensor tractography (DTT) of the entire pediatric cervical and thoracic spinal cord toward discrimination of typically developing (TD) controls and subjects with spinal cord injury (SCI). A total of 43 pediatric subjects, including 23 TD subjects ranging in age from 6 to 16 years old and 20 subjects with SCI ranging in age from 7 to 16 years, were recruited and scanned using a 3.0 Tesla magnetic resonance scanner. Reduced field of view diffusion tensor images were acquired axially to cover the entire spinal cord across two slabs. For DTI analysis, motion correction was performed by coregistration of the diffusion-weighted images to the reference image (b0). Streamline deterministic tractography results were generated from the preprocessed data. DTI and DTT parameters of the whole cord, including fractional anisotropy (FA), mean diffusivity (MD), tract length, and tract density, were calculated, averaged across the whole spinal cord, and compared between the TD and SCI groups. Statistically significant decreases have been shown in FA (TD = 0.46 ± 0.11; SCI = 0.37 ± 0.09; p < 0.0001) and tract density (TD = 405.93 ± 243.84; SCI = 268.90 ± 270.34; p < 0.0001). However, the mean length of tracts and MD did not show significant differences. When investigating differences in DTI and DTT parameters above and below the injury site, it was shown that the FA and tract density in patients with cervical SCI decreased significantly in the thoracic region. An identical trend was observed in the cervical region for patients with thoracic SCI as well. When comparing TD and SCI subjects, FA and tract density were the most sensitive parameters in detecting functional changes of the spinal cord in chronic pediatric SCI. The results show that both DTI and DTT have the potential to be imaging biomarkers in the diagnosis of SCI.

Effective Relief of Pain and Associated Symptoms With Closed-Loop Spinal Cord Stimulation System: Preliminary Results of the Avalon Study

OBJECTIVES

Conventional spinal cord stimulation (SCS) delivers a fixed-input of energy into the dorsal column. Physiologic effects such as heartbeat, respiration, spinal cord movement, and history of stimulation can cause both the perceived intensity and recruitment of stimulation to increase or decrease, with clinical consequences. A new SCS system controls stimulation dose by measuring the recruitment of fibers in the dorsal column and by using the amplitude of the evoked compound action potentials (ECAPs) to maintain stimulation within an individualized therapeutic range. Safety and efficacy of this closed-loop system was evaluated through six-month postimplantation.

MATERIALS AND METHODS

Chronic pain subjects with back and/or leg pain who were successfully trialed received a permanent system (Evoke; Saluda Medical, Sydney, Australia). Ratings of pain (100-mm visual analogue scale [VAS] and Brief Pain Instrument [BPI]), quality of life (EuroQol instrument [EQ-5D-5L]), function (Oswestry Disability Index [ODI]), and sleep (Pittsburgh Sleep Quality Index [PSQI]) were collected at baseline and repeated three and six months after implantation.

RESULTS

Fifty-one subjects underwent a trial procedure; permanent implants were placed in 36 subjects. The proportion of subjects with ≥50% relief was 92.6% (back) and 91.3% (leg) at three months, and 85.7% (back) and 82.6% (leg) at six months. The proportion with ≥80% pain relief was 70.4% (back) and 56.5% (leg) at three months, and 64.3% (back) and 60.9% (leg) at six months. Statistically significant improvements in mean BPI, EQ-5D-5L, ODI, and PSQI were also observed at both time points.

CONCLUSIONS

The majority of subjects experienced profound pain relief at three and six months, providing preliminary evidence for the effectiveness of the closed-loop SCS system. The exact mechanism of action for these outcomes is still being explored, although one likely hypothesis holds that ECAP feedback control may minimize recruitment of Aβ nociceptors and Aδ fibers during daily use of SCS.

Clinical utility for diffusion MRI sequence in emergency and inpatient spine protocols

Diffusion imaging of the spine has the potential to change clinical management, but is challenging due to the small size of the cord and susceptibility artifacts from adjacent structures. Reduced field-of-view (rFOV) diffusion can improve image quality by decreasing the echo train length. Over the past 2 years, we have acquired a rFOV diffusion sequence for MRI spine protocols on most inpatients and emergency room patients. We provide selected imaging diagnoses to illustrate the utility of including diffusion spine MRI in clinical practice. Our experiences support using diffusion MRI to improve diagnostic certainty and facilitate prompt treatment or clinical management.

Investigating resting-state functional connectivity in the cervical spinal cord at 3T

The study of spontaneous fluctuations in the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord. Two ultra-high field functional magnetic resonance imaging (fMRI) studies in humans have provided evidence for reproducible resting-state connectivity between the dorsal horns as well as between the ventral horns, and a study in non-human primates has shown that these resting-state signals are impacted by spinal cord injury. As these studies were carried out at ultra-high field strengths using region-of-interest (ROI) based analyses, we investigated whether such resting-state signals could also be observed at the clinically more prevalent field strength of 3 T. In a reanalysis of a sample of 20 healthy human participants who underwent a resting-state fMRI acquisition of the cervical spinal cord, we were able to observe significant dorsal horn connectivity as well as ventral horn connectivity, but no consistent effects for connectivity between dorsal and ventral horns, thus replicating the human 7 T results. These effects were not only observable when averaging along the acquired length of the spinal cord, but also when we examined each of the acquired spinal segments separately, which showed similar patterns of connectivity. Finally, we investigated the robustness of these resting-state signals against variations in the analysis pipeline by varying the type of ROI creation, temporal filtering, nuisance regression and connectivity metric. We observed that – apart from the effects of band-pass filtering – ventral horn connectivity showed excellent robustness, whereas dorsal horn connectivity showed moderate robustness. Together, our results provide evidence that spinal cord resting-state connectivity is a robust and spatially consistent phenomenon that could be a valuable tool for investigating the effects of pathology, disease progression, and treatment response in neurological conditions with a spinal component, such as spinal cord injury.

Denoising spinal cord fMRI data: Approaches to acquisition and analysis

Functional magnetic resonance imaging (fMRI) of the human spinal cord is a difficult endeavour due to the cord's small cross-sectional diameter, signal drop-out as well as image distortion due to magnetic field inhomogeneity, and the confounding influence of physiological noise from cardiac and respiratory sources. Nevertheless, there is great interest in spinal fMRI due to the spinal cord's role as the principal sensorimotor interface between the brain and the body and its involvement in a variety of sensory and motor pathologies. In this review, we give an overview of the various methods that have been used to address the technical challenges in spinal fMRI, with a focus on reducing the impact of physiological noise. We start out by describing acquisition methods that have been tailored to the special needs of spinal fMRI and aim to increase the signal-to-noise ratio and reduce distortion in obtained images. Following this, we concentrate on image processing and analysis approaches that address the detrimental effects of noise. While these include variations of standard pre-processing methods such as motion correction and spatial filtering, the main focus lies on denoising techniques that can be applied to task-based as well as resting-state data sets. We review both model-based approaches that rely on externally acquired respiratory and cardiac signals as well as data-driven approaches that estimate and correct for noise using the data themselves. We conclude with an outlook on techniques that have been successfully applied for noise reduction in brain imaging and whose use might be beneficial for fMRI of the human spinal cord.

Organization of the intrinsic functional network in the cervical spinal cord: A resting state functional MRI study

Resting state functional magnetic resonance imaging (rsfMRI) has been extensively applied to investigate the organization of functional networks in the brain. As an essential part of the central nervous system (CNS), the spinal cord has not been well explored about its intrinsic functional network. In this study, we aim to thoroughly investigate the characteristics of the intrinsic functional network in the spinal cord using rsfMRI. Functional connectivity and graph theory analysis were employed to evaluate the organization of the functional network, including its topology and network communication properties. Furthermore, the reproducibility of rsfMRI analysis on the spinal cord was also examined by intra-class correlation (ICC). Comprehensive evaluation of the intrinsic functional organization presented a non-uniform distribution of topological characteristics of the functional network, in which the upper levels (C2 and C3 vertebral levels) of the cervical spinal cord showed high levels of connectivity. The present results revealed the significance of the upper cervical cord in the intrinsic functional network of the human cervical spinal cord. In addition, this study demonstrated the efficiency of the cervical spinal cord functional network and the reproducibility of rsfMRI analysis on the spinal cord was also confirmed. As knowledge expansion of intrinsic functional network from the brain to the spinal cord, this study shed light on the organization of the spinal cord functional network in both normal development and clinical disorders.

Motion correction and frequency stabilization for MRS of the human spinal cord

Subject motion is challenging for MRS, because it can falsify results. For spinal cord MRS in particular, subject movement is critical, since even a small movement > 1 mm) can lead to a voxel shift out of the desired measurement region. Therefore, the identification of motion corrupted MRS scans is essential. In this investigation, MR navigators acquired simultaneously with the MRS data are used to identify a displacement of the spinal cord due to subject motion. It is shown that navigators are able to recognize substantial subject motion (>1 mm) without impairing the MRS measurement. In addition, navigators are easy to apply to the measurement, because no additional hardware and just a minor additional user effort are needed. Moreover, no additional scan time is required, because navigators can be applied in the deadtime of the MRS sequence. Furthermore, in this work, retrospective motion correction combined with frequency stabilization is presented by combining navigators with non-water-suppressed (1)H-MRS, resulting in an improved spectral quality of the spinal cord measurements.

fMRI Localization of Spinal Cord Processing Underlying Female Sexual Arousal

Using functional magnetic resonance imaging, the authors aimed to determine the roles of the human spinal cord in mediating sexual responses in women. Functional magnetic resonance imaging of the entire lower thoracic, lumbar, and sacral spinal cord was performed using a sexual stimulation paradigm designed to elicit psychological and physical components of sexual arousal. Responses were measured in 9 healthy adult women during 3 consecutive conditions: (a) erotic audiovisual, (b) manual clitoral, and (c) audiovisual plus manual stimulation. Functional magnetic resonance imaging results in healthy subjects demonstrate that this method is sensitive for mapping sexual function in the spinal cord, and identify several key regions involved in human sexual response, including the intermediolateral cell column, the dorsal commissural nucleus, and the sacral parasympathetic nucleus. Using spinal functional magnetic resonance imaging, this study identified many of the spinal cord regions involved in female sexual responses. Results from audiovisual and manual clitoral stimulation correspond with previous data regarding lumbar and sacral neurologic changes during sexual arousal. This study provides the first characterization of neural activity in the human spinal cord underlying healthy female sexual responses and sets a foundation for future studies aimed at mapping changes that result from sexual dysfunction, spinal cord trauma or disease.

Diffusion tensor imaging of the cervical spinal cord in children

PURPOSE

Obtaining fast, reliable, high-resolution diffusion tensor imaging (DTI) of the pediatric cervical spinal cord (CSC) is challenging, given the multitude of technical limitations involved. Overcoming these limitations may further potentiate DTI as a valuable quantitative tool in evaluating the pediatric CSC.

METHODS

Sixteen patients (9 girls and 7 boys) with hypoxic brain injury, craniocervical junction malformations, and head trauma were included in this retrospective study. Region of interests were placed from C1-C2 through C7-T1 consecutively at the cervical intervertebral disc levels. DTI metrics were compared with a pediatric DTI database of healthy controls. Clinical background and outcomes were tabulated.

RESULTS

Patients with hypoxic brain injury, Chiari I and II malformations, and head trauma demonstrated lower fractional anisotropy values than that of healthy controls at certain cervical intervertebral disc levels.

CONCLUSIONS

DTI may be a promising modality for providing additional information beyond that of conventional magnetic resonance imaging in pediatric central nervous system disorders.

Poro-elastic modeling of Syringomyelia - a systematic study of the effects of pia mater, central canal, median fissure, white and gray matter on pressure wave propagation and fluid movement within the cervical spinal cord

Syringomyelia, fluid-filled cavities within the spinal cord, occurs frequently in association with a Chiari I malformation and produces some of its most severe neurological symptoms. The exact mechanism causing syringomyelia remains unknown. Since syringomyelia occurs frequently in association with obstructed cerebrospinal fluid (CSF) flow, it has been hypothesized that syrinx formation is mechanically driven. In this study we model the spinal cord tissue either as a poro-elastic medium or as a solid linear elastic medium, and simulate the propagation of pressure waves through an anatomically plausible 3D geometry, with boundary conditions based on in vivo CSF pressure measurements. Then various anatomic and tissue properties are modified, resulting in a total of 11 variations of the model that are compared. The results show that an open segment of the central canal and a stiff pia (relative to the cord) both increase the radial pressure gradients and enhance interstitial fluid flow in the central canal. The anterior median fissure, anisotropic permeability of the white matter, and Poisson ratio play minor roles.

Magnetic resonance imaging assessment of spinal cord and cauda equina motion in supine patients with spinal metastases planned for spine stereotactic body radiation therapy

PURPOSE

To assess motion of the spinal cord and cauda equina, which are critical neural tissues (CNT), which is important when evaluating the planning organ-at-risk margin required for stereotactic body radiation therapy.

METHODS AND MATERIALS

We analyzed CNT motion in 65 patients with spinal metastases (11 cervical, 39 thoracic, and 24 lumbar spinal segments) in the supine position using dynamic axial and sagittal magnetic resonance imaging (dMRI, 3T Verio, Siemens) over a 137-second interval. Motion was segregated according to physiologic cardiorespiratory oscillatory motion (characterized by the average root mean square deviation) and random bulk shifts associated with gross patient motion (characterized by the range). Displacement was evaluated in the anteroposterior (AP), lateral (LR), and superior-inferior (SI) directions by use of a correlation coefficient template matching algorithm, with quantification of random motion measure error over 3 separate trials. Statistical significance was defined according to P<.05.

RESULTS

In the AP, LR, and SI directions, significant oscillatory motion was observed in 39.2%, 35.1%, and 10.8% of spinal segments, respectively, and significant bulk motions in all cases. The median oscillatory CNT motions in the AP, LR, and SI directions were 0.16 mm, 0.17 mm, and 0.44 mm, respectively, and the maximal statistically significant oscillatory motions were 0.39 mm, 0.41 mm, and 0.77 mm, respectively. The median bulk displacements in the AP, LR, and SI directions were 0.51 mm, 0.59 mm, and 0.66 mm, and the maximal statistically significant displacements were 2.21 mm, 2.87 mm, and 3.90 mm, respectively. In the AP, LR, and SI directions, bulk displacements were greater than 1.5 mm in 5.4%, 9.0%, and 14.9% of spinal segments, respectively. No significant differences in axial motion were observed according to cord level or cauda equina.

CONCLUSIONS

Oscillatory CNT motion was observed to be relatively minor. Our results support the importance of controlling bulk patient motion and the practice of applying a planning organ-at-risk margin.

Intrinsically organized resting state networks in the human spinal cord

Spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signals of the brain have repeatedly been observed when no task or external stimulation is present. These fluctuations likely reflect baseline neuronal activity of the brain and correspond to functionally relevant resting-state networks (RSN). It is not known however, whether intrinsically organized and spatially circumscribed RSNs also exist in the spinal cord, the brain's principal sensorimotor interface with the body. Here, we use recent advances in spinal fMRI methodology and independent component analysis to answer this question in healthy human volunteers. We identified spatially distinct RSNs in the human spinal cord that were clearly separated into dorsal and ventral components, mirroring the functional neuroanatomy of the spinal cord and likely reflecting sensory and motor processing. Interestingly, dorsal (sensory) RSNs were separated into right and left components, presumably related to ongoing hemibody processing of somatosensory information, whereas ventral (motor) RSNs were bilateral, possibly related to commissural interneuronal networks involved in central pattern generation. Importantly, all of these RSNs showed a restricted spatial extent along the spinal cord and likely conform to the spinal cord's functionally relevant segmental organization. Although the spatial and temporal properties of the dorsal and ventral RSNs were found to be significantly different, these networks showed significant interactions with each other at the segmental level. Together, our data demonstrate that intrinsically highly organized resting-state fluctuations exist in the human spinal cord and are thus a hallmark of the entire central nervous system.