BOLD signal responses to controlled hypercapnia in human spinal cord

Imaging protocols for non-traumatic spinal cord injury: current state of the art and future directions

INTRODUCTION

Non-traumatic spinal cord injury (NTSCI) is a term used to describe damage to the spinal cord from sources other than trauma. Neuroimaging techniques such as computerized tomography (CT) and magnetic resonance imaging (MRI) have improved our ability to diagnose and manage NTSCIs. Several practice guidelines utilize MRI in the diagnostic evaluation of traumatic and non-traumatic SCI to direct surgical intervention.

AREAS COVERED

The authors review practices surrounding the imaging of various causes of NTSCI as well as recent advances and future directions for the use of novel imaging modalities in this realm. The authors also present discussions around the use of simple radiographs and advanced MRI modalities in clinical settings, and briefly highlight areas of active research that seek to advance our understanding and improve patient care.

EXPERT OPINION

Although several obstacles must be overcome, it appears highly likely that novel quantitative imaging features and advancements in artificial intelligence (AI) as well as machine learning (ML) will revolutionize degenerative cervical myelopathy (DCM) care by providing earlier diagnosis, accurate localization, monitoring for deterioration and neurological recovery, outcome prediction, and standardized practice. Some intriguing findings in these areas have been published, including the identification of possible serum and cerebrospinal fluid biomarkers, which are currently in the early phases of translation.

MRI mapping of hemodynamics in the human spinal cord

Impaired spinal cord vascular function contributes to numerous neurological pathologies, making it important to be able to noninvasively characterize these changes. Here, we propose a functional magnetic resonance imaging (fMRI)-based method to map spinal cord vascular reactivity (SCVR). We used a hypercapnic breath-holding task, monitored with end-tidal CO2 (PETCO2), to evoke a systemic vasodilatory response during concurrent blood oxygenation level-dependent (BOLD) fMRI. SCVR amplitude and hemodynamic delay were mapped at the group level in 27 healthy participants as proof-of-concept of the approach, and then in two highly-sampled participants to probe feasibility/stability of individual SCVR mapping. Across the group and the highly-sampled individuals, a strong ventral SCVR amplitude was initially observed without accounting for local regional variation in the timing of the vasodilatory response. Shifted breathing traces (PETCO2) were used to account for temporal differences in the vasodilatory response across the spinal cord, producing maps of SCVR delay. These delay maps reveal an earlier ventral and later dorsal response and demonstrate distinct gray matter regions concordant with territories of arterial supply. The SCVR fMRI methods described here enable robust mapping of spatiotemporal hemodynamic properties of the human spinal cord. This noninvasive approach has exciting potential to provide early insight into pathology-driven vascular changes in the cord, which may precede and predict future irreversible tissue damage and guide the treatment of several neurological pathologies involving the spine.

Thermal stimulus task fMRI in the cervical spinal cord at 7 Tesla

Although functional magnetic resonance imaging (fMRI) is widely applied in the brain, fMRI of the spinal cord is more technically demanding. Proximity to the vertebral column and lungs results in strong spatial inhomogeneity and temporal fluctuations in B0 . Increasing field strength enables higher spatial resolution and improved sensitivity to blood oxygenation level-dependent (BOLD) signal, but amplifies the effects of B0 inhomogeneity. In this work, we present the first task fMRI in the spinal cord at 7 T. Further, we compare the performance of single-shot and multi-shot 2D echo-planar imaging (EPI) protocols, which differ in sensitivity to spatial and temporal B0 inhomogeneity. The cervical spinal cords of 11 healthy volunteers were scanned at 7 T using single-shot 2D EPI at 0.75 mm in-plane resolution and multi-shot 2D EPI at 0.75 and 0.6 mm in-plane resolutions. All protocols used 3 mm slice thickness. For each protocol, the BOLD response to 13 10-s noxious thermal stimuli applied to the right thumb was acquired in a 10-min fMRI run. Image quality, temporal signal to noise ratio (SNR), and BOLD activation (percent signal change and z-stat) at both individual- and group-level were evaluated between the protocols. Temporal SNR was highest in single-shot and multi-shot 0.75 mm protocols. In group-level analyses, activation clusters appeared in all protocols in the ipsilateral dorsal quadrant at the expected C6 neurological level. In individual-level analyses, activation clusters at the expected level were detected in some, but not all subjects and protocols. Single-shot 0.75 mm generally produced the highest mean z-statistic, while multi-shot 0.60 mm produced the best-localized activation clusters and the least geometric distortion. Larger than expected within-subject segmental variation of BOLD activation along the cord was observed. Group-level sensory task fMRI of the cervical spinal cord is feasible at 7 T with single-shot or multi-shot EPI. The best choice of protocol will likely depend on the relative importance of sensitivity to activation versus spatial localization of activation for a given experiment. PRACTITIONER POINTS: First stimulus task fMRI results in the spinal cord at 7 T. Single-shot 0.75 mm 2D EPI produced the highest mean z-statistic. Multi-shot 0.60 mm 2D EPI provided the best-localized activation and least distortion.

Recent developments and future avenues for human corticospinal neuroimaging

Non-invasive neuroimaging serves as a valuable tool for investigating the mechanisms within the central nervous system (CNS) related to somatosensory and motor processing, emotions, memory, cognition, and other functions. Despite the extensive use of brain imaging, spinal cord imaging has received relatively less attention, regardless of its potential to study peripheral communications with the brain and the descending corticospinal systems. To comprehensively understand the neural mechanisms underlying human sensory and motor functions, particularly in pathological conditions, simultaneous examination of neuronal activity in both the brain and spinal cord becomes imperative. Although technically demanding in terms of data acquisition and analysis, a growing but limited number of studies have successfully utilized specialized acquisition protocols for corticospinal imaging. These studies have effectively assessed sensorimotor, autonomic, and interneuronal signaling within the spinal cord, revealing interactions with cortical processes in the brain. In this mini-review, we aim to examine the expanding body of literature that employs cutting-edge corticospinal imaging to investigate the flow of sensorimotor information between the brain and spinal cord. Additionally, we will provide a concise overview of recent advancements in functional magnetic resonance imaging (fMRI) techniques. Furthermore, we will discuss potential future perspectives aimed at enhancing our comprehension of large-scale neuronal networks in the CNS and their disruptions in clinical disorders. This collective knowledge will aid in refining combined corticospinal fMRI methodologies, leading to the development of clinically relevant biomarkers for conditions affecting sensorimotor processing in the CNS.

Reliability of task-based fMRI in the dorsal horn of the human spinal cord

The application of functional magnetic resonance imaging (fMRI) to the human spinal cord is still a relatively small field of research and faces many challenges. Here we aimed to probe the limitations of task-based spinal fMRI at 3T by investigating the reliability of spinal cord blood oxygen level dependent (BOLD) responses to repeated nociceptive stimulation across two consecutive days in 40 healthy volunteers. We assessed the test-retest reliability of subjective ratings, autonomic responses, and spinal cord BOLD responses to short heat pain stimuli (1s duration) using the intraclass correlation coefficient (ICC). At the group level, we observed robust autonomic responses as well as spatially specific spinal cord BOLD responses at the expected location, but no spatial overlap in BOLD response patterns across days. While autonomic indicators of pain processing showed good-to-excellent reliability, both β-estimates and z-scores of task-related BOLD responses showed poor reliability across days in the target region (gray matter of the ipsilateral dorsal horn). When taking into account the sensitivity of gradient-echo echo planar imaging (GE-EPI) to draining vein signals by including the venous plexus in the analysis, we observed BOLD responses with good reliability across days. Taken together, these results demonstrate that heat pain stimuli as short as one second are able to evoke a robust and spatially specific BOLD response, which is however strongly variable within participants across time, resulting in low reliability in the dorsal horn gray matter. Further improvements in data acquisition and analysis techniques are thus necessary before event-related spinal cord fMRI as used here can be reliably employed in longitudinal designs or clinical settings.

Spatial distribution of hand-grasp motor task activity in spinal cord functional magnetic resonance imaging

Upper extremity motor paradigms during spinal cord functional magnetic resonance imaging (fMRI) can provide insight into the functional organization of the cord. Hand-grasping is an important daily function with clinical significance, but previous studies of similar squeezing movements have not reported consistent areas of activity and are limited by sample size and simplistic analysis methods. Here, we study spinal cord fMRI activation using a unimanual isometric hand-grasping task that is calibrated to participant maximum voluntary contraction (MVC). Two task modeling methods were considered: (1) a task regressor derived from an idealized block design (Ideal) and (2) a task regressor based on the recorded force trace normalized to individual MVC (%MVC). Across these two methods, group motor activity was highly lateralized to the hemicord ipsilateral to the side of the task. Activation spanned C5-C8 and was primarily localized to the C7 spinal cord segment. Specific differences in spatial distribution are also observed, such as an increase in C8 and dorsal cord activity when using the %MVC regressor. Furthermore, we explored the impact of data quantity and spatial smoothing on sensitivity to hand-grasp motor task activation. This analysis shows a large increase in number of active voxels associated with the number of fMRI runs, sample size, and spatial smoothing, demonstrating the impact of experimental design choices on motor activation.

Thermal Stimulus Task fMRI in the Cervical Spinal Cord at 7 Tesla

PURPOSE

Although functional MRI is widely applied in the brain, fMRI of the spinal cord is more technically demanding. Proximity to the vertebral column and lungs results in strong spatial inhomogeneity and temporal fluctuations in B0. Increasing field strength enables higher spatial resolution and improved sensitivity to BOLD signal, but amplifies the effects of B0 inhomogeneity. In this work, we present the first stimulus task fMRI in the spinal cord at 7 T. Further, we compare the performance of single-shot and multi-shot 2D EPI protocols, as they differ in sensitivity to spatial and temporal B0 inhomogeneity.

METHODS

The cervical spinal cords of 11 healthy volunteers were scanned at 7 T using single-shot 2D EPI at 0.75 mm in-plane resolution and multi-shot 2D EPI at 0.75 and 0.6 mm in-plane resolutions. For each protocol, the BOLD response to thirteen 10-second noxious thermal stimuli applied to the right thumb was acquired in a 10-minute fMRI run. Image quality, temporal SNR, and BOLD activation (percent signal change and z-stat) at both individual- and group-level were evaluated between the protocols.

RESULTS

Temporal SNR was highest in single-shot and multi-shot 0.75 mm protocols. In group-level analyses, activation clusters appeared in all protocols in the ipsilateral dorsal quadrant at the expected C6 neurological level. In individual-level analyses, activation clusters at the expected level were detected in some, but not all subjects and protocols. Single-shot 0.75 mm generally produced the highest mean z-statistic, while multi-shot 0.60 mm produced the best-localized activation clusters and the least geometric distortion. Larger than expected within-subject segmental variation of BOLD activation along the cord was observed.

CONCLUSION

Group-level sensory task fMRI of the cervical spinal cord is feasible at 7 T with single-shot or multi-shot EPI. The best choice of protocol will likely depend on the relative importance of sensitivity to activation versus spatial localization of activation for a given experiment.

The current state of spinal cord functional magnetic resonance imaging and its application in clinical research

Since its development, spinal cord functional magnetic resonance imaging (fMRI) has utilized various methodologies and stimulation protocols to develop a deeper understanding of a healthy human spinal cord that lays a foundation for its use in clinical research and practice. In this review, we conducted a comprehensive literature search on spinal cord fMRI studies and summarized the recent advancements and resulting scientific achievements of spinal cord fMRI in the following three aspects: the current state of spinal cord fMRI methodologies and stimulation protocols, knowledge about the healthy spinal cord's functions obtained via spinal cord fMRI, and fMRI's exemplary usage in spinal cord diseases and injuries. We conclude with a discussion that, while technical challenges exist, novel fMRI technologies for and new knowledge about the healthy human spinal cord have been established. Empowered by these developments, investigations of pathological and injury states within the spinal cord have become the next important direction of spinal cord fMRI. Recent clinical investigations into spinal cord pathologies, for example, fibromyalgia, multiple sclerosis, spinal cord injury, and cervical spondylotic myelopathy, have already provided deep insights into spinal cord impairments and the time course of impairment-caused changes. We expect that future spinal cord fMRI advancement and research development will further enhance our understanding of various spinal cord diseases and provide the foundation for evaluating existing and developing new treatment plans.

Spatial distribution of hand-grasp motor task activity in spinal cord functional magnetic resonance imaging

Upper extremity motor paradigms during spinal cord functional magnetic resonance imaging (fMRI) can provide insight into the functional organization of the cord. Hand-grasping is an important daily function with clinical significance, but previous studies of similar squeezing movements have not reported consistent areas of activity and are limited by sample size and simplistic analysis methods. Here, we study spinal cord fMRI activation using a unimanual isometric hand-grasping task that is calibrated to participant maximum voluntary contraction (MVC). Two task modeling methods were considered: (1) a task regressor derived from an idealized block design (Ideal) and (2) a task regressor based on the recorded force trace normalized to individual MVC (%MVC). Across these two methods, group motor activity was highly lateralized to the hemicord ipsilateral to the side of the task. Activation spanned C5-C8 and was primarily localized to the C7 spinal cord segment. Specific differences in spatial distribution are also observed, such as an increase in C8 and dorsal cord activity when using the %MVC regressor. Furthermore, we explored the impact of data quantity and spatial smoothing on sensitivity to hand-grasp motor task activation. This analysis shows a large increase in number of active voxels associated with the number of fMRI runs, sample size, and spatial smoothing, demonstrating the impact of experimental design choices on motor activation.

Improving T2*-weighted human cortico-spinal acquisitions with a dedicated algorithm for region-wise shimming

Cortico-spinal fMRI acquisitions aim to investigate direct interactions between brain and spinal cord, e.g. during motor output or pain processing, by covering both regions in a single measurement. Due to their large distance and location in the body, a dynamic shim update of constant and linear shim terms is required when using echo-planar imaging (EPI) to achieve reasonable image quality in both target regions. A previously presented approach with region-wise shim settings is based on a standard single-region shim algorithm and suffers from (i) non-optimal shim settings because it combines linear and second-order shim terms optimized for different volumes, and (ii) significant user interactions making it rather cumbersome, time consuming, and error-prone. Here, a dedicated ("CoSpi") shim algorithm for cortico-spinal fMRI is presented that performs joint optimization of static second-order shim terms and one set of linear and constant shim terms for each region in a single run and with minimal user interaction. Field map and T2*-weighted EPI measurements were performed on a clinical 3 T whole-body MR system in water phantoms and five healthy volunteers using the conventional region-wise and CoSpi shim settings as well as "gold standard" shim settings optimized for one of the target regions only. With CoSpi shim settings, (i) overall field inhomogeneity was reduced by about 65% / 75% (brain / spinal cord volume) compared to the conventional region-wise approach and in vivo was within 5% of the values obtained with the single-volume shim settings, (ii) geometric distortions derived from voxel displacement maps were reduced on average by about 35% / 70%, (iii) the temporal SNR determined from an EPI time series that may reflect the impact of through-slice dephasing, was increased by about 17% / 10%, and (iv) the variation of the mean field between slices, a measure targeting the predisposition to insufficient fat saturation and GRAPPA-related ghosting artifacts, was reduced by about 90% / 45%. Thus, the presented algorithm not only speeds up and simplifies the shim procedure considerably, but also provides a better field homogeneity and image quality, which both could help to significantly improve the applicability of cortico-spinal fMRI.

Advanced spinal cord MRI in multiple sclerosis: Current techniques and future directions

Spinal cord magnetic resonance imaging (MRI) has a central role in multiple sclerosis (MS) clinical practice for diagnosis and disease monitoring. Advanced MRI sequences capable of visualizing and quantifying tissue macro- and microstructure and reflecting different pathological disease processes have been used in MS research; however, the spinal cord remains under-explored, partly due to technical obstacles inherent to imaging this structure. We propose that the study of the spinal cord merits equal ambition in overcoming technical challenges, and that there is much information to be exploited to make valuable contributions to our understanding of MS. We present a narrative review on the latest progress in advanced spinal cord MRI in MS, covering in the first part structural, functional, metabolic and vascular imaging methods. We focus on recent studies of MS and those making significant technical steps, noting the challenges that remain to be addressed and what stands to be gained from such advances. Throughout we also refer to other works that presend more in-depth review on specific themes. In the second part, we present several topics that, in our view, hold particular potential. The need for better imaging of gray matter is discussed. We stress the importance of developing imaging beyond the cervical spinal cord, and explore the use of ultra-high field MRI. Finally, some recommendations are given for future research, from study design to newer developments in analysis, and the need for harmonization of sequences and methods within the field. This review is aimed at researchers and clinicians with an interest in gaining an overview of the current state of advanced MRI research in this field and what is primed to be the future of spinal cord imaging in MS research.

Brain-spinal cord interaction in long-term motor sequence learning in human: An fMRI study

The spinal cord is important for sensory guidance and execution of skilled movements. Yet its role in human motor learning is not well understood. Despite evidence revealing an active involvement of spinal circuits in the early phase of motor learning, whether long-term learning engages similar changes in spinal cord activation and functional connectivity remains unknown. Here, we investigated spinal-cerebral functional plasticity associated with learning of a specific sequence of visually-guided joystick movements (sequence task) over six days of training. On the first and last training days, we acquired high-resolution functional images of the brain and cervical cord simultaneously, while participants practiced the sequence or a random task while electromyography was recorded from wrist muscles. After six days of training, the subjects' motor performance improved in the sequence compared to the control condition. These behavioral changes were associated with decreased co-contractions and increased reciprocal activations between antagonist wrist muscles. Importantly, early learning was characterized by activation in the C8 level, whereas a more rostral activation in the C6-C7 was found during the later learning phase. Motor sequence learning was also supported by increased spinal cord functional connectivity with distinct brain networks, including the motor cortex, superior parietal lobule, and the cerebellum at the early stage, and the angular gyrus and cerebellum at a later stage of learning. Our results suggest that the early vs. late shift in spinal activation from caudal to rostral cervical segments synchronized with distinct brain networks, including parietal and cerebellar regions, is related to progressive changes reflecting the increasing fine control of wrist muscles during motor sequence learning.

Imaging and Electrophysiology for Degenerative Cervical Myelopathy [AO Spine RECODE-DCM Research Priority Number 9]

STUDY DESIGN

Narrative review.

OBJECTIVE

The current review aimed to describe the role of existing techniques and emerging methods of imaging and electrophysiology for the management of degenerative cervical myelopathy (DCM), a common and often progressive condition that causes spinal cord dysfunction and significant morbidity globally.

METHODS

A narrative review was conducted to summarize the existing literature and highlight future directions.

RESULTS

Anatomical magnetic resonance imaging (MRI) is well established in the literature as the key imaging tool to identify spinal cord compression, disc herniation/bulging, and inbuckling of the ligamentum flavum, thus facilitating surgical planning, while radiographs and computed tomography (CT) provide complimentary information. Electrophysiology techniques are primarily used to rule out competing diagnoses. However, signal change and measures of cord compression on conventional MRI have limited utility to characterize the degree of tissue injury, which may be helpful for diagnosis, prognostication, and repeated assessments to identify deterioration. Early translational studies of quantitative imaging and electrophysiology techniques show potential of these methods to more accurately reflect changes in spinal cord microstructure and function.

CONCLUSION

Currently, clinical management of DCM relies heavily on anatomical MRI, with additional contributions from radiographs, CT, and electrophysiology. Novel quantitative assessments of microstructure, perfusion, and function have the potential to transform clinical practice, but require robust validation, automation, and standardization prior to uptake.

Simultaneous brain, brainstem and spinal cord pharmacological-fMRI reveals involvement of an endogenous opioid network in attentional analgesia

Pain perception is decreased by shifting attentional focus away from a threatening event. This attentional analgesia engages parallel descending control pathways from anterior cingulate (ACC) to locus coeruleus, and ACC to periaqueductal grey (PAG) – rostral ventromedial medulla (RVM), indicating possible roles for noradrenergic or opioidergic neuromodulators. To determine which pathway modulates nociceptive activity in humans, we used simultaneous whole brain-spinal cord pharmacological-fMRI (N = 39) across three sessions. Noxious thermal forearm stimulation generated somatotopic-activation of dorsal horn (DH) whose activity correlated with pain report and mirrored attentional pain modulation. Activity in an adjacent cluster reported the interaction between task and noxious stimulus. Effective connectivity analysis revealed that ACC interacts with PAG and RVM to modulate spinal cord activity. Blocking endogenous opioids with Naltrexone impairs attentional analgesia and disrupts RVM-spinal and ACC-PAG connectivity. Noradrenergic augmentation with Reboxetine did not alter attentional analgesia. Cognitive pain modulation involves opioidergic ACC-PAG-RVM descending control which suppresses spinal nociceptive activity.

An optimized MP2RAGE sequence for studying both brain and cervical spinal cord in a single acquisition at 3T

Magnetization Prepared 2 Rapid Acquisition Gradient Echo (MP2RAGE) is a T₁ mapping technique that has been used broadly on brain and recently on cervical spinal cord (cSC). The growing interest for combined investigation of brain and SC in numerous pathologies of the central nervous system such as multiple sclerosis (MS), amyotrophic lateral sclerosis (ALS) and traumatic injuries, now brings about the need for optimization with regards to this specific investigation. This implies large spatial coverage with high spatial resolution and short acquisition time, high CNR and low B₁⁺ sensitivity, as well as high reproducibility and robust post-processing tools for T₁ quantification in different regions of brain and SC. In this work, a dedicated protocol (referred to as Pr-BSC) has been optimized for simultaneous brain and cSC T₁ MP2RAGE acquisition at 3T. After computer simulation optimization, the protocol was applied for in vivo validation experiments and compared to previously published state of the art protocols focusing on either the brain (Pr-B) or the cSC (Pr-SC). Reproducibility and in-ROI standard deviations were assessed on healthy volunteers in the perspective of future clinical use. The mean T₁ values, obtained by the Pr-BSC, in brain white, gray and deep gray matters were: (mean ± in-ROI SD) 792 ± 27 ms, 1339 ± 139 ms and 1136 ± 88 ms, respectively. In cSC, T₁ values for white matter corticospinal, posterior sensory, lateral sensory and rubro/reticulospinal tracts were 902 ± 41 ms, 920 ± 35 ms, 903 ± 46 ms, 891 ± 41 ms, respectively, and 954 ± 32 ms for anterior and intermediate gray matter. The Pr-BSC protocol showed excellent agreement with previously proposed Pr-B on brain and Pr-SC on cSC, with very high inter-scan reproducibility (coefficients of variation of 0.52 ± 0.36% and 1.12 ± 0.62% on brain and cSC, respectively). This optimized protocol covering both brain and cSC with a sub-millimetric isotropic spatial resolution in one acquisition of less than 8 min, opens up great perspectives for clinical applications focusing on degenerative tissue such as encountered in MS and ALS.

Large-scale functional ultrasound imaging of the spinal cord reveals in-depth spatiotemporal responses of spinal nociceptive circuits in both normal and inflammatory states

Despite a century of research on the physiology/pathophysiology of the spinal cord in chronic pain condition, the properties of the spinal cord were rarely studied at the large-scale level from a neurovascular point of view. This is mostly due to the limited spatial and/or temporal resolution of the available techniques. Functional ultrasound imaging (fUS) is an emerging neuroimaging approach that allows, through the measurement of cerebral blood volume, the study of brain functional connectivity or functional activations with excellent spatial (100 μm) and temporal (1 msec) resolutions and a high sensitivity. The aim of this study was to increase our understanding of the spinal cord physiology through the study of the properties of spinal hemodynamic response to the natural or electrical stimulation of afferent fibers. Using a combination of fUS and ultrasound localization microscopy, the first step of this study was the fine description of the vascular structures in the rat spinal cord. Then, using either natural or electrical stimulations of different categories of afferent fibers (Aβ, Aδ, and C fibers), we could define the characteristics of the typical hemodynamic response of the rat spinal cord experimentally. We showed that the responses are fiber-specific, located ipsilaterally in the dorsal horn, and that they follow the somatotopy of afferent fiber entries in the dorsal horn and that the C-fiber response is an N-methyl-D-aspartate receptor-dependent mechanism. Finally, fUS imaging of the mesoscopic hemodynamic response induced by natural tactile stimulations revealed a potentiated response in inflammatory condition, suggesting an enhanced response to allodynic stimulations.

Multi-shot acquisitions for stimulus-evoked spinal cord BOLD fMRI

PURPOSE

To demonstrate the feasibility of 3D multi-shot magnetic resonance imaging acquisitions for stimulus-evoked blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) in the human spinal cord in vivo.

METHODS

Two fMRI studies were performed at 3T. The first study was a hypercapnic gas challenge where data were acquired from healthy volunteers using a multi-shot 3D fast field echo (FFE) sequence as well as single-shot multi-slice echo-planar imaging (EPI). In the second study, another cohort of healthy volunteers performed an upper extremity motor task while fMRI data were acquired using a 3D multi-shot acquisition.

RESULTS

Both 2D-EPI and 3D-FFE were shown to be sensitive to BOLD signal changes in the cervical spinal cord, and had comparable contrast-to-noise ratios in gray matter. FFE exhibited much less signal drop-out and weaker geometric distortions compared to EPI. In the motor paradigm study, the mean number of active voxels was highest in the ventral gray matter horns ipsilateral to the side of the task and at the spinal level associated with innervation of finger extensors.

CONCLUSIONS

Highly multi-shot acquisition sequences such as 3D-FFE are well suited for stimulus-evoked spinal cord BOLD fMRI.

Cortico-spinal imaging to study pain

Functional magnetic resonance imaging of the brain has helped to reveal mechanisms of pain perception in health and disease. Recently, imaging approaches have been developed that allow recording neural activity simultaneously in the brain and in the spinal cord. These approaches offer the possibility to examine pain perception in the entire central pain system and in addition, to investigate cortico-spinal interactions during pain processing. Although cortico-spinal imaging is a promising technique, it bears challenges concerning data acquisition and data analysis strategies. In this review, we discuss studies that applied simultaneous imaging of the brain and spinal cord to explore central pain processing. Furthermore, we describe different MR-related acquisition techniques, summarize advantages and disadvantages of approaches that have been implemented so far and present software that has been specifically developed for the analysis of spinal fMRI data to address challenges of spinal data analysis.

Dynamic per slice shimming for simultaneous brain and spinal cord fMRI

PURPOSE

Simultaneous brain and spinal cord functional MRI is emerging as a new tool to study the central nervous system but is challenging. Poor B0 homogeneity and small size of the spinal cord are principal obstacles to this nascent technology. Here we extend a dynamic shimming approach, first posed by Finsterbusch, by shimming per slice for both the brain and spinal cord.

METHODS

We shim dynamically by a simple and fast optimization of linear field gradients and frequency offset separately for each slice in order to minimize off-resonance for both the brain and spinal cord. Simultaneous acquisition of brain and spinal cord fMRI is achieved with high spatial resolution in the spinal cord by means of an echo-planar RF pulse for reduced FOV. Brain slice acquisition is full FOV.

RESULTS

T2*-weighted images of brain and spinal cord are acquired with high clarity and minimal observable image artifacts. Fist-clenching fMRI experiments reveal task-consistent activation in motor cortices, cerebellum, and C6-T1 spinal segments.

CONCLUSIONS

High quality functional results are obtained for a sensory-motor task. Consistent activation in both the brain and spinal cord is observed at individual levels, not only at group level. Because reduced FOV excitation is applicable to any spinal cord section, future continuation of these methods holds great potential.

A practical protocol for measurements of spinal cord functional connectivity

Resting state functional magnetic resonance imaging (fMRI) has been used to study human brain function for over two decades, but only recently has this technique been successfully translated to the human spinal cord. The spinal cord is structurally and functionally unique, so resting state fMRI methods developed and optimized for the brain may not be appropriate when applied to the cord. This report therefore investigates the relative impact of different acquisition and processing choices (including run length, echo time, and bandpass filter width) on the detectability of resting state spinal cord networks at 3T. Our results suggest that frequencies beyond 0.08 Hz should be included in resting state analyses, a run length of ~8-12 mins is appropriate for reliable detection of the ventral (motor) network, and longer echo times - yet still shorter than values typically used for fMRI in the brain - may increase the detectability of the dorsal (sensory) network. Further studies are required to more fully understand and interpret the nature of resting state spinal cord networks in health and in disease, and the protocols described in this report are designed to assist such studies.

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.

Resting-state functional connectivity in the rat cervical spinal cord at 9.4 T

PURPOSE

Numerous studies have adopted resting-state functional MRI methods to infer functional connectivity between cortical regions, but very few have translated them to the spinal cord, despite its critical role in the central nervous system. Resting-state functional connectivity between gray matter horns of the spinal cord has previously been shown to be detectable in humans and nonhuman primates, but it has not been reported previously in rodents.

METHODS

Resting-state functional MRI of the cervical spinal cord of live anesthetized rats was performed at 9.4 T. The quality of the functional images acquired was assessed, and quantitative analyses of functional connectivity in C4-C7 of the spinal cord were derived.

RESULTS

Robust gray matter horn-to-horn connectivity patterns were found that were statistically significant when compared with adjacent control regions. Specifically, dorsal-dorsal and ventral-ventral connectivity measurements were most prominent, while ipsilateral dorsal-ventral connectivity was also observed but to a lesser extent. Quantitative evaluation of reproducibility also revealed moderate robustness in the bilateral sensory and motor networks that was weaker in the dorsal-ventral connections.

CONCLUSIONS

This study reports the first evidence of resting-state functional circuits within gray matter in the rat spinal cord, and verifies their detectability using resting-state functional MRI at 9.4 T. Magn Reson Med 79:2773-2783, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Spinal cord MRI at 7T

Magnetic resonance imaging (MRI) of the human spinal cord at 7T has been demonstrated by a handful of research sites worldwide, and the spinal cord remains one of the areas in which higher fields and resolution could have high impact. The small diameter of the cord (∼1 cm) necessitates high spatial resolution to minimize partial volume effects between gray and white matter, and so MRI of the cord can greatly benefit from increased signal-to-noise ratio and contrasts at ultra-high field (UHF). Herein we review the current state of UHF spinal cord imaging. Technical challenges to successful UHF spinal cord MRI include radiofrequency (B1) nonuniformities and a general lack of optimized radiofrequency coils, amplified physiological noise, and an absence of methods for robust B0 shimming along the cord to mitigate image distortions and signal losses. Numerous solutions to address these challenges have been and are continuing to be explored, and include novel approaches for signal excitation and acquisition, dynamic shimming and specialized shim coils, and acquisitions with increased coverage or optimal slice angulations.

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.

SCT: Spinal Cord Toolbox, an open-source software for processing spinal cord MRI data

For the past 25 years, the field of neuroimaging has witnessed the development of several software packages for processing multi-parametric magnetic resonance imaging (mpMRI) to study the brain. These software packages are now routinely used by researchers and clinicians, and have contributed to important breakthroughs for the understanding of brain anatomy and function. However, no software package exists to process mpMRI data of the spinal cord. Despite the numerous clinical needs for such advanced mpMRI protocols (multiple sclerosis, spinal cord injury, cervical spondylotic myelopathy, etc.), researchers have been developing specific tools that, while necessary, do not provide an integrative framework that is compatible with most usages and that is capable of reaching the community at large. This hinders cross-validation and the possibility to perform multi-center studies. In this study we introduce the Spinal Cord Toolbox (SCT), a comprehensive software dedicated to the processing of spinal cord MRI data. SCT builds on previously-validated methods and includes state-of-the-art MRI templates and atlases of the spinal cord, algorithms to segment and register new data to the templates, and motion correction methods for diffusion and functional time series. SCT is tailored towards standardization and automation of the processing pipeline, versatility, modularity, and it follows guidelines of software development and distribution. Preliminary applications of SCT cover a variety of studies, from cross-sectional area measures in large databases of patients, to the precise quantification of mpMRI metrics in specific spinal pathways. We anticipate that SCT will bring together the spinal cord neuroimaging community by establishing standard templates and analysis procedures.

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.

Simultaneous Brain-Cervical Cord fMRI Reveals Intrinsic Spinal Cord Plasticity during Motor Sequence Learning

The spinal cord participates in the execution of skilled movements by translating high-level cerebral motor representations into musculotopic commands. Yet, the extent to which motor skill acquisition relies on intrinsic spinal cord processes remains unknown. To date, attempts to address this question were limited by difficulties in separating spinal local effects from supraspinal influences through traditional electrophysiological and neuroimaging methods. Here, for the first time, we provide evidence for local learning-induced plasticity in intact human spinal cord through simultaneous functional magnetic resonance imaging of the brain and spinal cord during motor sequence learning. Specifically, we show learning-related modulation of activity in the C6–C8 spinal region, which is independent from that of related supraspinal sensorimotor structures. Moreover, a brain–spinal cord functional connectivity analysis demonstrates that the initial linear relationship between the spinal cord and sensorimotor cortex gradually fades away over the course of motor sequence learning, while the connectivity between spinal activity and cerebellum gains strength. These data suggest that the spinal cord not only constitutes an active functional component of the human motor learning network but also contributes distinctively from the brain to the learning process. The present findings open new avenues for rehabilitation of patients with spinal cord injuries, as they demonstrate that this part of the central nervous system is much more plastic than assumed before. Yet, the neurophysiological mechanisms underlying this intrinsic functional plasticity in the spinal cord warrant further investigations.

Simultaneous neuroimaging of brain and spinal cord reveals intrinsic plasticity in the spinal cord during motor sequence learning in humans, independent from that of related sensorimotor structures in the brain.

When we acquire a new motor skill—for example, learning how to play a musical instrument—new synaptic connections are induced in a distributed network of brain areas. There is ample evidence from human neuroimaging studies for this high plasticity of the brain, but what about the spinal cord, the main link between the brain and the peripheral nervous system? Literature on animal models has recently hinted that spinal cord neurons can learn during various conditioning paradigms. However, human learning models by tradition assume that the spinal cord acts as a passive relay of information from the cortex to the muscles. In this study, we simultaneously acquired functional images of both the brain and the cervical spinal cord through functional magnetic resonance imaging, and we provide evidence for local spinal cord plasticity during a well-studied motor learning task in humans. We also demonstrate a dynamic change in the interaction of the brain and spinal cord regions over the course of motor learning. The present findings have important clinical implications for rehabilitation of patients with spinal cord injuries, as they demonstrate that this part of the central nervous system is much more plastic than it was assumed before.

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.

Resting state functional connectivity in the human spinal cord

Functional magnetic resonance imaging using blood oxygenation level dependent (BOLD) contrast is well established as one of the most powerful methods for mapping human brain function. Numerous studies have measured how low-frequency BOLD signal fluctuations from the brain are correlated between voxels in a resting state, and have exploited these signals to infer functional connectivity within specific neural circuits. However, to date there have been no previous substantiated reports of resting state correlations in the spinal cord. In a cohort of healthy volunteers, we observed robust functional connectivity between left and right ventral (motor) horns, and between left and right dorsal (sensory) horns. Our results demonstrate that low-frequency BOLD fluctuations are inherent in the spinal cord as well as the brain, and by analogy to cortical circuits, we hypothesize that these correlations may offer insight into the execution and maintenance of sensory and motor functions both locally and within the cerebrum.

DOI:http://dx.doi.org/10.7554/eLife.02812.001

Brain imaging methods such as functional magnetic resonance imaging (fMRI) can provide us with a picture of what the brain is doing when a person is carrying out a specific task. For example, an fMRI scan recorded whilst someone is reading is likely to show activity in regions in the left hemisphere of the brain that are known to be involved in language comprehension. fMRI can also be used to measure patterns of neuronal activity when someone is awake but not engaged in a specific task. This approach, known as resting state fMRI, can be used to examine which regions of the resting brain are active at the same time. Researchers are interested in these patterns of brain activity because they reflect neural circuits that work together to produce different functions and behaviors.

Over 4000 papers have used resting state fMRI to study the human brain. However, to date there has been no conclusive investigation of resting state activity in the spinal cord. This is largely because the spinal cord is much smaller than the brain, and most fMRI scanners are not sensitive enough to study it in detail. Consequently, little is known about intrinsic neural circuits in the resting spinal cord.

Now Barry et al. have used advances in fMRI technology to show that resting state functional connectivity does indeed exist in the spinal cord. Correlations were found in the resting levels of activity between spatially distinct areas of the cord, specifically between the ventral horns and between the dorsal horns. The ventral horns relay motor signals to the body, whilst the dorsal horns receive sensory signals from the body.

These findings also have clinical applications. Some patients with incomplete spinal cord injuries can recover near normal function, but the mechanisms responsible for this recovery are unclear because clinicians have not been able to probe neuronal connections in the spinal cord in a non-invasive manner. The work of Barry et al. should help with efforts to understand the neuronal changes that support recovery from spinal cord injury.

DOI:http://dx.doi.org/10.7554/eLife.02812.002

Diffusion tensor imaging of the cervical spinal cord in healthy adult population: normative values and measurement reproducibility at 3T MRI

BACKGROUND

Compared to diffusion tensor imaging (DTI) of the brain, there is a paucity of reports addressing the applicability of DTI in the evaluation of the spinal cord. Most normative data of cervical spinal cord DTI consist of relatively small and arbitrarily collected populations. Comprehensive normative data are necessary for clinical decision-making.

PURPOSE

To establish normal values for cervical spinal cord DTI metrics with region of interest (ROI)- and fiber tractography (FT)-based measurements and to assess the reproducibility of both measurement methods.

MATERIAL AND METHODS

Forty healthy adults underwent cervical spinal cord 3T MRI. Sagittal and axial conventional T2 sequences and DTI in the axial plane were performed. Whole cord fractional anisotropy (FA) and apparent diffusion coefficient (ADC) values were determined at different cervical levels from C2 to C7 using the ROI method. DTI metrics (FA, axial, and radial diffusivities based on eigenvalues λ1, λ2, and λ3, and ADC) of the lateral and posterior funicles were measured at C3 level. FA and ADC of the whole cord and the lateral and posterior funicles were also measured using quantitative tractography. Intra- and inter-observer variation of the measurement methods were assessed.

RESULTS

Whole cord FA values decreased and ADC values increased in the rostral to caudal direction from C2 to C7. Between the individual white matter funicles no statistically significant difference for FA or ADC values was found. Both axial diffusivity and radial diffusivity of both lateral funicles differed significantly from those of the posterior funicle. Neither gender nor age correlated with any of the DTI metrics. Intra-observer variation of the measurements for whole cord FA and ADC showed almost perfect agreement with both ROI and tractography-based measurements. There was more variation in measurements of individual columns. Inter-observer agreement varied from moderate to strong for whole cord FA and ADC.

CONCLUSION

Both ROI- and FT-based measurements are applicable methods yielding reproducible results for cervical spinal cord DTI metrics. Normative values for both measurement methods are presented.

The current state-of-the-art of spinal cord imaging: methods

A first-ever spinal cord imaging meeting was sponsored by the International Spinal Research Trust and the Wings for Life Foundation with the aim of identifying the current state-of-the-art of spinal cord imaging, the current greatest challenges, and greatest needs for future development. This meeting was attended by a small group of invited experts spanning all aspects of spinal cord imaging from basic research to clinical practice. The greatest current challenges for spinal cord imaging were identified as arising from the imaging environment itself; difficult imaging environment created by the bone surrounding the spinal canal, physiological motion of the cord and adjacent tissues, and small cross-sectional dimensions of the spinal cord, exacerbated by metallic implants often present in injured patients. Challenges were also identified as a result of a lack of "critical mass" of researchers taking on the development of spinal cord imaging, affecting both the rate of progress in the field, and the demand for equipment and software to manufacturers to produce the necessary tools. Here we define the current state-of-the-art of spinal cord imaging, discuss the underlying theory and challenges, and present the evidence for the current and potential power of these methods. In two review papers (part I and part II), we propose that the challenges can be overcome with advances in methods, improving availability and effectiveness of methods, and linking existing researchers to create the necessary scientific and clinical network to advance the rate of progress and impact of the research.

Combined T2*-weighted measurements of the human brain and cervical spinal cord with a dynamic shim update

Important functions of the central nervous system such as sensory processing and motor execution, involve the spinal cord. Recent advances in human functional MRI have allowed to investigate spinal cord neuronal processes using the blood-oxygenation-level-dependent (BOLD) contrast. However, to assess the functional connectivity between the brain and the spinal cord, functional MRI measurements covering both regions in the same experiment are required. Unfortunately, the ideal MRI setup differs considerably for the brain and the spinal cord with respect to resolution, field-of-view, relevant receive coils, and, in particular, shim adjustments required to minimize distortion artifacts. Here, these issues are addressed for combined T2*-weighted MRI measurements of the human brain and the cervical spinal cord by using adapted parameter settings (field-of-view, in-plane resolution, slice thickness, and receiver bandwidth) for each region, a dynamic receive coil element selection where for each slice only the elements with significant signal contributions are considered, and, most importantly, the implementation of a dynamic update of the frequency and the linear shims in order to provide shim settings individually adapted to the brain and spinal cord subvolume. The feasibility of this setup for combined measurements is demonstrated in healthy volunteers at 3T. Although geometric distortions are slightly more pronounced and the temporal signal-to-noise ratio is lower as compared to measurements focusing to the brain or spinal cord only, the overall image quality can be expected to be sufficient for combined functional MRI experiments. Thus, the presented approach could help to unravel the functional coupling between the brain and the spinal cord.

CO2 blood oxygen level-dependent MR mapping of cerebrovascular reserve in a clinical population: safety, tolerability, and technical feasibility

PURPOSE

To evaluate the safety, tolerability, and technical feasibility of mapping cerebrovascular reactivity (CVR) in a clinical population by using a precise prospectively targeted CO(2) stimulus and blood oxygen level-dependent (BOLD) magnetic resonance (MR) imaging.

MATERIALS AND METHODS

A chart review was performed of all CVR studies from institutional review board-approved projects at a tertiary care hospital between January 1, 2006, and December 1, 2010. Informed consent was obtained. Records were searched for the incidence of adverse events and failed examinations. CVR maps were evaluated for diagnostic quality by two blinded observers and were categorized as good, diagnostic but suboptimal, or nondiagnostic. Outcomes were presented as raw data and descriptive statistics (means ± standard deviations). Intraclass correlation coefficient was used to determine interobserver variability.

RESULTS

Four hundred thirty-four consecutive CVR examinations from 294 patients (51.8% female patients) were studied. Patient age ranged from 9 to 88 years (mean age, 45.9 years ± 20.6). Transient symptoms, such as shortness of breath, headache, and dizziness, were reported in 48 subjects (11.1% of studies) during hypercapnic phases only. There were no neurologic ischemic events, myocardial infarctions, or other major complications. The success rate in generating CVR maps was 83.9% (364 of 434). Of the 70 (16.1%) failed examinations, 25 (35.7%) were due to discomfort; eight (11.4%), to head motion; two (2.9%), to inability to cooperate; seven (10.0%), to technical difficulties with equipment; and 28 (40.0%), to unknown or unspecified conditions. Among the 364 remaining successful examinations, good quality CVR maps were obtained in 340 (93.4%); diagnostic but suboptimal, in 12 (3.3%); and nondiagnostic, in 12 (3.3%).

CONCLUSION

CVR mapping by using a prospectively targeted CO(2) stimulus and BOLD MR imaging is safe, well tolerated, and technically feasible in a clinical patient population.

Cervical spinal cord injection of epidural corticosteroids: comprehensive longitudinal study including multiparametric magnetic resonance imaging

Despite widespread use, the efficacy of epidural corticosteroid injections (ESI) for osteoarthritis-associated neck or radicular pain remains uncertain, so even rare serious complications enter into discussions about use. However, various factors impede investigation and publication of serious adverse events. To that end, we developed new magnetic resonance imaging (MRI) techniques for spinal cord white matter quantification and used the best available physiological tests to characterize a cervical spinal cord lesion caused by inadvertent intramedullary injection of Depo-Medrol. A 29-year-old woman with mild cervical osteoarthritis had 2 years of headache and neck pain (concussion and whiplash) after 2 minor motor vehicle accidents. During C5-6 ESI, she developed new left-sided motor and sensory symptoms, and MRI demonstrated a new left dorsal spinal cord cavity. Mild left-sided motor and sensory symptoms have persisted for more than 2.5 years, during which time we performed serial neurological examinations, standard electrodiagnostics, somatosensory evoked potentials, and transcranial measurement of corticospinal central motor conduction time (CMCT). We used 3-Tesla MRI with a 32-channel coil developed for high-resolution cervical spinal cord structural imaging, diffusion tensor imaging (DTI), and magnetization transfer (MT). T(2)(∗)-weighted signal and DTI and MT metrics showed delayed spread of the lesion across 4 vertebral levels rostrally, consistent with Wallerian degeneration within the ascending left dorsal columns. However, only CMCT metrics detected objective correlates of her left hemiparesis and bilateral hyperreflexia. DTI and MT metrics may better distinguish between post-traumatic demyelination and axonal degeneration than conventional MRI. These tests should be considered to better characterize similar spinal cord injuries.

Reduction of physiological noise with independent component analysis improves the detection of nociceptive responses with fMRI of the human spinal cord

The evaluation of spinal cord neuronal activity in humans with functional magnetic resonance imaging (fMRI) is technically challenging. Major difficulties arise from cardiac and respiratory movement artifacts that constitute significant sources of noise. In this paper we assessed the Correction of Structured noise using spatial Independent Component Analysis (CORSICA). FMRI data of the cervical spinal cord were acquired in 14 healthy subjects using gradient-echo EPI. Nociceptive electrical stimuli were applied to the thumb. Additional data with short TR (250 ms, to prevent aliasing) were acquired to generate a spatial map of physiological noise derived from Independent Component Analysis (ICA). Physiological noise was subsequently removed from the long-TR data after selecting independent components based on the generated noise map. Stimulus-evoked responses were analyzed using the general linear model, with and without CORSICA and with a regressor generated from the cerebrospinal fluid region. Results showed higher sensitivity to detect stimulus-related activation in the targeted dorsal segment of the cord after CORSICA. Furthermore, fewer voxels showed stimulus-related signal changes in the CSF and outside the spinal region, suggesting an increase in specificity. ICA can be used to effectively reduce physiological noise in spinal cord fMRI time series.

Stimulus site and modality dependence of functional activity within the human spinal cord

Chronic pain is thought to arise because of maladaptive changes occurring within the peripheral nervous system and CNS. The transition from acute to chronic pain is known to involve the spinal cord (Woolf and Salter, 2000). Therefore, to investigate altered human spinal cord function and translate results obtained from other species, a noninvasive neuroimaging technique is desirable. We have investigated the functional response in the cervical spinal cord of 18 healthy human subjects (aged 22-40 years) to noxious thermal and non-noxious tactile stimulation of the left and right forearms. Physiological noise, which is a significant source of signal variability in the spinal cord, was accounted for in the general linear model. Group analysis, performed using a mixed-effects model, revealed distinct regions of activity that were dependent on both the side and the type of stimulation. In particular, thermal stimulation on the medial aspect of the wrist produced activity within the C6/C5 segment ipsilateral to the side of stimulation. Similar to data recorded in animals (Fitzgerald, 1982), painful thermal stimuli produced increased ipsilateral and decreased contralateral blood flow, which may reflect, respectively, excitatory and inhibitory processes. Nonpainful punctate stimulation of the thenar eminence provoked more diffuse activity but was still ipsilateral to the side of stimulation. These results present the first noninvasive evidence for a lateralized response to noxious and non-noxious stimuli in the human spinal cord. The development of these techniques opens the path to understanding, at a subject-specific level, central sensitization processes that contribute to chronic pain states.

Effect of cerebrovascular changes on brain DTI quantitation: a hypercapnia study

Quantitative diffusion tensor imaging (DTI) offers a valuable tool to probe the microstructural changes in neural tissues in vivo, where absolute quantitation accuracy and reproducibility are essential. It has been long recognized that measurement of apparent diffusion coefficient (ADC) using DTI could be influenced by the presence of water molecules in cerebrovasculature. However, little is known about to what extent such blood signal affects DTI quantitation. In this study, we quantitatively examined the effect of cerebral hemodynamic change on DTI indices by using a standard multislice echo planar imaging (EPI) spin echo (SE) DTI acquisition protocol and a rat model of hypercapnia. In response to 5% CO(2) challenge, mean, radial and axial diffusivities measured with diffusion factor (b-value) of b=1.0 ms/μm(2) were found to increase in whole brain (1.52%±0.22%, 1.66%±0.16% and 1.35%±0.37%, respectively), gray matter (1.56%±0.23%, 1.63%±0.14% and 1.47%±0.45%, respectively) and white matter regions (1.45%±0.28%, 1.88%±0.33% and 1.10%±0.26%, respectively). Fractional anisotropy (FA) was found to decrease by 1.67%±0.38%, 1.91%±0.59% and 1.46%±0.30% in whole brain, gray matter and white matter regions, respectively. In addition, these diffusivity increases and FA decreases became more pronounced at a lower b-value (b=0.3 ms/μm(2)). The results indicated that in vivo DTI quantitation in brain can be contaminated by vascular factors on the order of few percentages. Consequently, alterations in cerebrovasculature and hemodynamics can affect the DTI quantitation and its efficacy in characterizing the neural tissue microstructures in normal and diseased states. Caution should be taken in designing and interpreting quantitative DTI studies as all DTI indices can be potentially confounded by physiologic conditions and by cerebrovascular and hemodynamic characteristics.

Measurement and characterization of the human spinal cord SEEP response using event-related spinal fMRI

Although event-related fMRI is able to reliably detect brief changes in brain activity and is now widely used throughout systems and cognitive neuroscience, there have been no previous reports of event-related spinal cord fMRI. This is likely attributable to the various technical challenges associated with spinal fMRI (e.g., imaging a suitable length of the cord, reducing image artifacts from the vertebrae and intervertebral discs, and dealing with physiological noise from spinal cord motion). However, with many of these issues now resolved, the largest remaining impediment for event-related spinal fMRI is a deprived understanding of the spinal cord fMRI signal time course. Therefore, in this study, we used a proton density-weighted HASTE sequence, with functional contrast based on signal enhancement by extravascular water protons (SEEP), and a motion-compensating GLM analysis to (i) characterize the SEEP response function in the human cervical spinal cord and (ii) demonstrate the feasibility of event-related spinal fMRI. This was achieved by applying very brief (1 s) epochs of 22°C thermal stimulation to the palm of the hand and measuring the impulse response function. Our results suggest that the spinal cord SEEP response (time to peak ≈8 s; FWHM ≈4 s; and probably lacking pre- and poststimulus undershoots) is slower than previous estimates of SEEP or BOLD responses in the brain, but faster than previously reported spinal cord BOLD responses. Finally, by detecting and mapping consistent signal-intensity changes within and across subjects, and validating these regions with a block-designed experiment, this study represents the first successful demonstration of event-related spinal fMRI.

32-channel RF coil optimized for brain and cervical spinal cord at 3 T

Diffusion and functional magnetic resonance imaging of the spinal cord remain challenging due to the small cross-sectional size of the cord and susceptibility-related distortions. Although partially addressable through parallel imaging, few highly parallel array coils have been implemented for the cervical cord. Here, we developed a 32-channel coil that fully covers the brain and c-spine and characterized its performance in comparison with a commercially available head/neck/spine array. Image and temporal signal-to-noise ratio were, respectively, increased by 2× and 1.8× in the cervical cord. Averaged g-factors at 4× acceleration were lowered by 22% in the brain and by 39% in the spinal cord, enabling 1-mm isotropic R = 4 multi-echo magnetization prepared gradient echo of the full brain and c-spine in 3:20 min. Diffusion imaging of the cord at 0.6 × 0.6 × 5 mm(3) resolution and tractography of the full brain and c-spine at 1.7-mm isotropic resolution were feasible without noticeable distortion. Improvements of this nature potentially enhance numerous basic and clinical research studies focused on spinal and supraspinal regions.

Identification and characterisation of midbrain nuclei using optimised functional magnetic resonance imaging

Localising activity in the human midbrain with conventional functional MRI (fMRI) is challenging because the midbrain nuclei are small and located in an area that is prone to physiological artefacts. Here we present a replicable and automated method to improve the detection and localisation of midbrain fMRI signals. We designed a visual fMRI task that was predicted would activate the superior colliculi (SC) bilaterally. A limited number of coronal slices were scanned, orientated along the long axis of the brainstem, whilst simultaneously recording cardiac and respiratory traces. A novel anatomical registration pathway was used to optimise the localisation of the small midbrain nuclei in stereotactic space. Two additional structural scans were used to improve registration between functional and structural T1-weighted images: an echo-planar image (EPI) that matched the functional data but had whole-brain coverage, and a whole-brain T2-weighted image. This pathway was compared to conventional registration pathways, and was shown to significantly improve midbrain registration. To reduce the physiological artefacts in the functional data, we estimated and removed structured noise using a modified version of a previously described physiological noise model (PNM). Whereas a conventional analysis revealed only unilateral SC activity, the PNM analysis revealed the predicted bilateral activity. We demonstrate that these methods improve the measurement of a biologically plausible fMRI signal. Moreover they could be used to investigate the function of other midbrain nuclei.

Demyelination and degeneration in the injured human spinal cord detected with diffusion and magnetization transfer MRI

Characterizing demyelination/degeneration of spinal pathways in traumatic spinal cord injured (SCI) patients is crucial for assessing the prognosis of functional rehabilitation. Novel techniques based on diffusion-weighted (DW) magnetic resonance imaging (MRI) and magnetization transfer (MT) imaging provide sensitive and specific markers of white matter pathology. In this paper we combined for the first time high angular resolution diffusion-weighted imaging (HARDI), MT imaging and atrophy measurements to evaluate the cervical spinal cord of fourteen SCI patients and age-matched controls. We used high in-plane resolution to delineate dorsal and ventrolateral pathways. Significant differences were detected between patients and controls in the normal-appearing white matter for fractional anisotropy (FA, p<0.0001), axial diffusivity (p<0.05), radial diffusivity (p<0.05), generalized fractional anisotropy (GFA, p<0.0001), magnetization transfer ratio (MTR, p<0.0001) and cord area (p<0.05). No significant difference was detected in mean diffusivity (p=0.41), T1-weighted (p=0.76) and T2-weighted (p=0.09) signals. MRI metrics were remarkably well correlated with clinical disability (Pearson's correlations, FA: p<0.01, GFA: p<0.01, radial diffusivity: p=0.01, MTR: p=0.04 and atrophy: p<0.01). Stepwise linear regressions showed that measures of MTR in the dorsal spinal cord predicted the sensory disability whereas measures of MTR in the ventro-lateral spinal cord predicted the motor disability (ASIA score). However, diffusion metrics were not specific to the sensorimotor scores. Due to the specificity of axial and radial diffusivity and MT measurements, results suggest the detection of demyelination and degeneration in SCI patients. Combining HARDI with MT imaging is a promising approach to gain specificity in characterizing spinal cord pathways in traumatic injury.

Characterization of the hemodynamic response in the rat lumbar spinal cord using intrinsic optical imaging and laser speckle

Quantifying spinal cord functions is crucial for understanding neurophysiological mechanisms governing the intact and the injured spinal cord. Intrinsic optical imaging (IOI) and laser speckle provides measures of deoxyhemoglobin (HbR) and oxyhemoglobin (HbO(2)) concentrations, blood volume (BV) and blood flow (BF) at high spatial and temporal resolution. In this study we used IOI and laser speckle to characterize the hemodynamic response to neuronal activation in the lumbar spinal cord of anaesthetized rats (N=9). We report consistent temporal variations of HbR, HbO(2), BV and BF located ipsilaterally at L3-L5. Responses were significantly higher when stimulation intensity was increased. Vascular changes extended several millimetres from the epicenter, supporting the venous drainage observed in functional magnetic resonance imaging studies.