Denoising spinal cord fMRI data: Approaches to acquisition and analysis

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.

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 fair 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.

Characteristic BOLD signals are detectable in white matter of the spinal cord at rest and after a stimulus

We report the reliable detection of reproducible patterns of blood-oxygenation-level-dependent (BOLD) MRI signals within the white matter (WM) of the spinal cord during a task and in a resting state. Previous functional MRI studies have shown that BOLD signals are robustly detectable not only in gray matter (GM) in the brain but also in cerebral WM as well as the GM within the spinal cord, but similar signals in WM of the spinal cord have been overlooked. In this study, we detected BOLD signals in the WM of the spinal cord in squirrel monkeys and studied their relationships with the locations and functions of ascending and descending WM tracts. Tactile sensory stimulus -evoked BOLD signal changes were detected in the ascending tracts of the spinal cord using a general-linear model. Power spectral analysis confirmed that the amplitude at the fundamental frequency of the response to a periodic stimulus was significantly higher in the ascending tracts than the descending ones. Independent component analysis of resting-state signals identified coherent fluctuations from eight WM hubs which correspond closely to the known anatomical locations of the major WM tracts. Resting-state analyses showed that the WM hubs exhibited correlated signal fluctuations across spinal cord segments in reproducible patterns that correspond well with the known neurobiological functions of WM tracts in the spinal cord. Overall, these findings provide evidence of a functional organization of intraspinal WM tracts and confirm that they produce hemodynamic responses similar to GM both at baseline and under stimulus conditions.

Functional ultrasound imaging of the human spinal cord

Utilizing the first in-human functional ultrasound imaging (fUSI) of the spinal cord, we demonstrate the integration of spinal functional responses to electrical stimulation. We record and characterize the hemodynamic responses of the spinal cord to a neuromodulatory intervention commonly used for treating pain and increasingly used for the restoration of sensorimotor and autonomic function. We found that the hemodynamic response to stimulation reflects a spatiotemporal modulation of the spinal cord circuitry not previously recognized. Our analytical capability offers a mechanism to assess blood flow changes with a new level of spatial and temporal precision in vivo and demonstrates that fUSI can decode the functional state of spinal networks in a single trial, which is of fundamental importance for developing real-time closed-loop neuromodulation systems. This work is a critical step toward developing a vital technique to study spinal cord function and effects of clinical neuromodulation.

Systematic investigation of cerebellar functional alterations and their association with surgical outcomes in patients with degenerative cervical myelopathy: a resting-state fMRI study

OBJECTIVE

The aim of this study was to systematically investigate the changes in cerebellar neural activity and cerebellar-cortical functional connectivity (FC) in patients with degenerative cervical myelopathy (DCM) using resting-state functional magnetic resonance imaging (fMRI).

METHODS

In this study, we collected clinical data and resting-state fMRI data from 54 DCM patients and 50 healthy controls (HCs). We analyzed voxel-wise regional fMRI metrics, including amplitude of low frequency fluctuation (ALFF), fractional ALFF, regional homogeneity, functional connectivity density, and voxel-mirrored homotopic connectivity. In analysis 1, we examined the differences in regional fMRI metrics within the cerebellum between the DCM patient group and the healthy control group, as well as their correlation with preoperative neurological status and prognosis. In analysis 2, we investigated cerebellar-cortical functional connectivity differences between the two groups and their correlation with preoperative neurological status and prognosis. Lastly, in analysis 3, we explored the internetwork connectivity between the 'cerebellar-SMN' (sensorimotor network) system, examined the between-group differences, and investigated its correlation with preoperative neurological status and prognosis.

RESULTS

(1) Relative to HCs, DCM patients exhibited functional alterations in wide-spread cerebellar regions; (2) DCM patients exhibited altered cerebellar-cortical FC which was associated with the preoperative neurological status and prognosis; (3) DCM patients exhibited altered internetwork connectivity between 'cerebellar-SMN' system which was associated with duration of symptom.

CONCLUSION

Wide-spread cerebellar functional alterations occur in DCM pathogenesis and the deficits in cerebellar-SMN functional connectivity may be beneficial in future studies for predicting surgical outcomes in patients with DCM.

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.

Association between activity in the ventral premotor cortex and spinal cord activation during force generation-A combined cortico-spinal fMRI study

Force generation is a crucial element of dexterity and a highly relevant skill of the human motor system. How cerebral and spinal components interact and how spinal activation is associated with the activity in the cerebral primary motor and premotor areas is poorly understood. Here, we conducted combined cortico-spinal functional magnetic resonance imaging during a simple visually guided isometric force generation task in 20 healthy young subjects. Activation was localized in the right cervical spinal cord and left primary motor and premotor areas. The main finding is that spinal activation was negatively correlated with ventral premotor cortex activation. Spinal activation was furthermore significantly correlated with primary motor cortex activation, while increasing target forces led to an increase in the amount of activation. These data indicate that human premotor areas such as the ventral premotor cortex might be functionally connected to the lower cervical spinal cord contributing to distal upper limb functions, a finding that extends our understanding of human motor function beyond the animal literature.

Spinal Cord fMRI: A New Window into the Central Nervous System

With the brain, the spinal cord forms the central nervous system. Initially considered a passive relay between the brain and the periphery, the spinal cord is now recognized as being active and plastic. Yet, it remains largely overlooked by the human neuroscience community, in stark contrast with the wealth of research investigating the brain. In this review, we argue that fMRI, traditionally used to image cerebral function, can be extended beyond the brain to help unravel spinal mechanisms involved in human behaviors. To this end, we first outline strategies that have been proposed to tackle the challenges inherent to spinal cord fMRI. Then, we discuss how they have been utilized to provide insights into the functional organization of spinal sensorimotor circuits, highlighting their potential to address fundamental and clinical questions. By summarizing guidelines and applications of spinal cord fMRI, we hope to stimulate and support further research into this promising yet underexplored field.

Denoising task-related fMRI: Balancing noise reduction against signal loss

Preprocessing fMRI data requires striking a fine balance between conserving signals of interest and removing noise. Typical steps of preprocessing include motion correction, slice timing correction, spatial smoothing, and high-pass filtering. However, these standard steps do not remove many sources of noise. Thus, noise-reduction techniques, for example, CompCor, FIX, and ICA-AROMA have been developed to further improve the ability to draw meaningful conclusions from the data. The ability of these techniques to minimize noise while conserving signals of interest has been tested almost exclusively in resting-state fMRI and, only rarely, in task-related fMRI. Application of noise-reduction techniques to task-related fMRI is particularly important given that such procedures have been shown to reduce false positive rates. Little remains known about the impact of these techniques on the retention of signal in tasks that may be associated with systemic physiological changes. In this paper, we compared two ICA-based, that is FIX and ICA-AROMA, two CompCor-based noise-reduction techniques, that is aCompCor, and tCompCor, and standard preprocessing using a large (n = 101) fMRI dataset including noxious heat and non-noxious auditory stimulation. Results show that preprocessing using FIX performs optimally for data obtained using noxious heat, conserving more signals than CompCor-based techniques and ICA-AROMA, while removing only slightly less noise. Similarly, for data obtained during non-noxious auditory stimulation, FIX noise-reduction technique before analysis with a covariate of interest outperforms the other techniques. These results indicate that FIX might be the most appropriate technique to achieve the balance between conserving signals of interest and removing noise during task-related fMRI.

Detection of resting-state functional connectivity in the lumbar spinal cord with 3T MRI

Functional MRI (fMRI) of the spinal cord is an expanding area of research with potential to investigate neuronal activity in the central nervous system. We aimed to characterize the functional connectivity features of the human lumbar spinal cord using resting-state fMRI (rs-fMRI) at 3T, using region-based and data-driven analysis approaches. A 3D multi-shot gradient echo resting-state blood oxygenation level dependent-sensitive rs-fMRI protocol was implemented in 26 healthy participants. Average temporal signal-to-noise ratio in the gray matter was 16.35 ± 4.79 after denoising. Evidence of synchronous signal fluctuations in the ventral and dorsal horns with their contralateral counterparts was observed in representative participants using interactive, exploratory seed-based correlations. Group-wise average in-slice Pearson's correlations were 0.43 ± 0.17 between ventral horns, and 0.48 ± 0.16 between dorsal horns. Group spatial independent component analysis (ICA) was used to identify areas of coherent activity¸ and revealed components within the gray matter corresponding to anatomical regions. Lower-dimensionality ICA revealed bilateral components corresponding to ventral and dorsal networks. Additional separate ICAs were run on two subsets of the participant group, yielding two sets of components that showed visual consistency and moderate spatial overlap. This work shows feasibility of rs-fMRI to probe the functional features and organization of the lumbar spinal cord.

Decoding cerebro-spinal signatures of human behavior: application to motor sequence learning

Mapping the neural patterns that drive human behavior is a key challenge in neuroscience. Even the simplest of our everyday actions stem from the dynamic and complex interplay of multiple neural structures across the central nervous system (CNS). Yet, most neuroimaging research has focused on investigating cerebral mechanisms, while the way the spinal cord accompanies the brain in shaping human behavior has been largely overlooked. Although the recent advent of functional magnetic resonance imaging (fMRI) sequences that can simultaneously target the brain and spinal cord has opened up new avenues for studying these mechanisms at multiple levels of the CNS, research to date has been limited to inferential univariate techniques that cannot fully unveil the intricacies of the underlying neural states. To address this, we propose to go beyond traditional analyses and instead use a data-driven multivariate approach leveraging the dynamic content of cerebro-spinal signals using innovation-driven coactivation patterns (iCAPs). We demonstrate the relevance of this approach in a simultaneous brain-spinal cord fMRI dataset acquired during motor sequence learning (MSL), to highlight how large-scale CNS plasticity underpins rapid improvements in early skill acquisition and slower consolidation after extended practice. Specifically, we uncovered cortical, subcortical and spinal functional networks, which were used to decode the different stages of learning with a high accuracy and, thus, delineate meaningful cerebro-spinal signatures of learning progression. Our results provide compelling evidence that the dynamics of neural signals, paired with a data-driven approach, can be used to disentangle the modular organization of the CNS. While we outline the potential of this framework to probe the neural correlates of motor learning, its versatility makes it broadly applicable to explore the functioning of cerebro-spinal networks in other experimental or pathological conditions.

Identification of synchronous BOLD signal patterns in white matter of primate spinal cord

Functional MRI studies of the brain have shown that blood-oxygenation-level-dependent (BOLD) signals are robustly detectable not only in gray matter (GM) but also in white matter (WM). Here, we report the detection and characteristics of BOLD signals in WM of spinal cord (SC) of squirrel monkeys. Tactile stimulus-evoked BOLD signal changes were detected in the ascending sensory tracts of SC using a General-Linear Model (GLM) as well as Independent Component Analysis (ICA). ICA of resting state signals identified coherent fluctuations from eight WM hubs which correspond closely with known anatomical locations of SC WM tracts. Resting state analyses showed that the WM hubs exhibited correlated signal fluctuations within and between SC segments in specific patterns that correspond well with the known neurobiological functions of WM tracts in SC. Overall, these findings suggest WM BOLD signals in SC show similar features as GM both at baseline and under stimulus conditions.

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.

Diffusion Kurtosis Imaging of Neonatal Spinal Cord in Clinical Routine

Diffusion kurtosis imaging (DKI) has undisputed advantages over the more classical diffusion magnetic resonance imaging (dMRI) as witnessed by the fast-increasing number of clinical applications and software packages widely adopted in brain imaging. However, in the neonatal setting, DKI is still largely underutilized, in particular in spinal cord (SC) imaging, because of its inherently demanding technological requirements. Due to its extreme sensitivity to non-Gaussian diffusion, DKI proves particularly suitable for detecting complex, subtle, fast microstructural changes occurring in this area at this early and critical stage of development, which are not identifiable with only DTI. Given the multiplicity of congenital anomalies of the spinal canal, their crucial effect on later developmental outcome, and the close interconnection between the SC region and the brain above, managing to apply such a method to the neonatal cohort becomes of utmost importance. This study will (i) mention current methodological challenges associated with the application of advanced dMRI methods, like DKI, in early infancy, (ii) illustrate the first semi-automated pipeline built on Spinal Cord Toolbox for handling the DKI data of neonatal SC, from acquisition setting to estimation of diffusion measures, through accurate adjustment of processing algorithms customized for adult SC, and (iii) present results of its application in a pilot clinical case study. With the proposed pipeline, we preliminarily show that DKI is more sensitive than DTI-related measures to alterations caused by brain white matter injuries in the underlying cervical SC.

Opioid analgesia alters corticospinal coupling along the descending pain system in healthy participants

Opioids are potent analgesic drugs with widespread cortical, subcortical, and spinal targets. In particular, the central pain system comprising ascending and descending pain pathways has high opioid receptor densities and is thus crucial for opioid analgesia. Here, we investigated the effects of the opioid remifentanil in a large sample (n = 78) of healthy male participants using combined corticospinal functional MRI. This approach offers the possibility to measure BOLD responses simultaneously in the brain and spinal cord, allowing us to investigate the role of corticospinal coupling in opioid analgesia. Our data show that opioids altered activity in regions involved in pain processing such as somatosensory regions, including the spinal cord and pain modulation such as prefrontal regions. Moreover, coupling strength along the descending pain system, that is, between the anterior cingulate cortex, periaqueductal gray, and spinal cord, was stronger in participants who reported stronger analgesia during opioid treatment while participants that received saline showed reduced coupling when experiencing less pain. These results indicate that coupling along the descending pain pathway is a potential mechanism of opioid analgesia and can differentiate between opioid analgesia and unspecific reductions in pain such as habituation.

Quantitative MR Markers in Non-Myelopathic Spinal Cord Compression: A Narrative Review

Degenerative spinal cord compression is a frequent pathological condition with increasing prevalence throughout aging. Initial non-myelopathic cervical spinal cord compression (NMDC) might progress over time into potentially irreversible degenerative cervical myelopathy (DCM). While quantitative MRI (qMRI) techniques demonstrated the ability to depict intrinsic tissue properties, longitudinal in-vivo biomarkers to identify NMDC patients who will eventually develop DCM are still missing. Thus, we aim to review the ability of qMRI techniques (such as diffusion MRI, diffusion tensor imaging (DTI), magnetization transfer (MT) imaging, and magnetic resonance spectroscopy (1H-MRS)) to serve as prognostic markers in NMDC. While DTI in NMDC patients consistently detected lower fractional anisotropy and higher mean diffusivity at compressed levels, caused by demyelination and axonal injury, MT and 1H-MRS, along with advanced and tract-specific diffusion MRI, recently revealed microstructural alterations, also rostrally pointing to Wallerian degeneration. Recent studies also disclosed a significant relationship between microstructural damage and functional deficits, as assessed by qMRI and electrophysiology, respectively. Thus, tract-specific qMRI, in combination with electrophysiology, critically extends our understanding of the underlying pathophysiology of degenerative spinal cord compression and may provide predictive markers of DCM development for accurate patient management. However, the prognostic value must be validated in longitudinal studies.

Activity-dependent spinal cord neuromodulation rapidly restores trunk and leg motor functions after complete paralysis

Epidural electrical stimulation (EES) targeting the dorsal roots of lumbosacral segments restores walking in people with spinal cord injury (SCI). However, EES is delivered with multielectrode paddle leads that were originally designed to target the dorsal column of the spinal cord. Here, we hypothesized that an arrangement of electrodes targeting the ensemble of dorsal roots involved in leg and trunk movements would result in superior efficacy, restoring more diverse motor activities after the most severe SCI. To test this hypothesis, we established a computational framework that informed the optimal arrangement of electrodes on a new paddle lead and guided its neurosurgical positioning. We also developed software supporting the rapid configuration of activity-specific stimulation programs that reproduced the natural activation of motor neurons underlying each activity. We tested these neurotechnologies in three individuals with complete sensorimotor paralysis as part of an ongoing clinical trial ( www.clinicaltrials.gov identifier NCT02936453). Within a single day, activity-specific stimulation programs enabled these three individuals to stand, walk, cycle, swim and control trunk movements. Neurorehabilitation mediated sufficient improvement to restore these activities in community settings, opening a realistic path to support everyday mobility with EES in people with SCI.

A visualization tool for assessment of spinal cord functional magnetic resonance imaging data quality

Functional magnetic resonance imaging (fMRI) is an extensively used neuroimaging technique to non-invasively detect neural activity. Data quality is highly variable, and fMRI analysis typically consists of a number of complex processing steps. It is crucial to visually assess images throughout analysis to ensure that data quality at each step is satisfactory. For fMRI analysis of the brain, there is a simple tool to visualize four-dimensional data on a two-dimensional plot for qualitative analysis. Despite the practicality of this method, it cannot be directly applied to fMRI data of the spinal cord, and a comparable approach does not exist for spinal cord fMRI analysis. The additional challenges encountered in spinal cord imaging, including the small size of the cord and the influence of physiological noise sources, drive the importance of developing a similar visualization technique for spinal cord fMRI. Here, we introduce a highly versatile image analysis tool to visualize spinal cord fMRI data as a simple heatmap and to co-visualize relevant traces such as physiological or motion timeseries. We present multiple variations of the plot, data features that can be identified with the heatmap, and examples of the useful qualitative analyses that can be performed using this method. The spinal cord plot can be easily integrated into an fMRI analysis pipeline and can streamline visual inspection and qualitative analysis of functional imaging data.Clinical Relevance- Implementation of this data visualization method is a simple addition to spinal cord fMRI analysis that could be used to identify normal vs. abnormal signal variation in pathologies that impact the cord, such as spinal cord injury or multiple sclerosis.

Investigating the human spinal sensorimotor pathways through functional magnetic resonance imaging

Most of our knowledge about the human spinal ascending (sensory) and descending (motor) pathways comes from non-invasive electrophysiological investigations. However, recent methodological advances in acquisition and analyses of functional magnetic resonance imaging (fMRI) data from the spinal cord, either alone or in combination with the brain, have allowed us to gain further insights into the organization of this structure. In the current review, we conducted a systematic search to produced somatotopic maps of the spinal fMRI activity observed through different somatosensory, motor and resting-state paradigms. By cross-referencing these human neuroimaging findings with knowledge acquired through neurophysiological recordings, our review demonstrates that spinal fMRI is a powerful tool for exploring, in vivo, the human spinal cord pathways. We report strong cross-validation between task-related and resting-state fMRI in accordance with well-known hemicord, postero-anterior and rostro-caudal organization of these pathways. We also highlight the specific advantages of using spinal fMRI in clinical settings to characterize better spinal-related impairments, predict disease progression, and guide the implementation of therapeutic interventions.

Functional networks in non-human primate spinal cord and the effects of injury

Spontaneous fluctuations of Blood Oxygenation-Level Dependent (BOLD) MRI signal in a resting state have previously been detected and analyzed to describe intrinsic functional networks in the spinal cord of rodents, non-human primates and human subjects. In this study we combined high resolution imaging at high field with data-driven Independent Component Analysis (ICA) to i) delineate fine-scale functional networks within and between segments of the cervical spinal cord of monkeys, and also to ii) characterize the longitudinal effects of a unilateral dorsal column injury on these networks. Seven distinct functional hubs were revealed within each spinal segment, with new hubs detected at bilateral intermediate and gray commissure regions in addition to the bilateral dorsal and ventral horns previously reported. Pair-wise correlations revealed significantly stronger connections between hubs on the dominant hand side. Unilateral dorsal-column injuries disrupted predominantly inter-segmental rather than intra-segmental functional connectivities as revealed by correlation strengths and graph-theory based community structures. The effects of injury on inter-segmental connectivity were evident along the length of the cord both below and above the lesion region. Connectivity strengths recovered over time and there was revival of inter-segmental communities as animals recovered function. BOLD signals of frequency 0.01-0.033 Hz were found to be most affected by injury. The results in this study provide new insights into the intrinsic functional architecture of spinal cord and underscore the potential of functional connectivity measures to characterize changes in networks after an injury and during recovery.

Structural and resting state functional connectivity beyond the cortex

Mapping the structural and functional connectivity of the central nervous system has become a key area within neuroimaging research. While detailed network structures across the entire brain have been probed using animal models, non-invasive neuroimaging in humans has thus far been dominated by cortical investigations. Beyond the cortex, subcortical nuclei have traditionally been less accessible due to their smaller size and greater distance from radio frequency coils. However, major neuroimaging developments now provide improved signal and the resolution required to study these structures. Here, we present an overview of the connectivity between the amygdala, brainstem, cerebellum, spinal cord and the rest of the brain. While limitations to their imaging and analyses remain, we also provide some recommendations and considerations for mapping brain connectivity beyond the cortex.

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.

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.

Effect of Physiological Noise on Thoracolumbar Spinal Cord Functional Magnetic Resonance Imaging in 3T Magnetic Field

Introduction:

Functional Magnetic Resonance Imaging (fMRI) methods have been used to study sensorimotor processing in the spinal cord. However, these techniques confront unwanted noises to the measured signal from the physiological fluctuations. In the spinal cord imaging, most of the challenges are consequences of cardiac and respiratory movement artifacts that are considered as significant sources of noise, especially in the thoracolumbar region. In this study, we investigated the effect of each source of physiological noise and their contribution to the outcome of the analysis of the blood-oxygen-level-dependent signal in the human thoracolumbar spinal cord.

Methods:

Fifteen young healthy male volunteers participated in the study, and pain stimuli were delivered on the L5 dermatome between the two malleoli. Respiratory and cardiac signals were recorded during the imaging session, and the generated respiration and cardiac regressors were included in the general linear model for quantification of the effect of each of them on the task-analysis results. The sum of active voxels of the clusters was calculated in the spinal cord in three correction states (respiration correction only, cardiac correction only, and respiration and cardiac noise corrections) and analyzed with analysis of variance statistical test and receiver operating characteristic curve.

Results:

The results illustrated that cardiac noise correction had an effective role in increasing the active voxels (Mean±SD = 23.46±9.46) compared to other noise correction methods. Cardiac effects were higher than other physiological noise sources

Conclusion:

In summary, our results indicate great respiration effects on the lumbar and thoracolumbar spinal cord fMRI, and its contribution to the heartbeat effect can be a significant variable in the individual fMRI data analysis. Displacement of the spinal cord and the effects of this noise in the thoracolumbar and lumbar spinal cord fMRI results are significant and cannot be ignored.

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 Functional Connectivity of Resting-State Spinal Cord fMRI Reveals Fine-Grained Intrinsic Architecture

The neuroimaging community has shown tremendous interest in exploring the brain's spontaneous activity using functional magnetic resonance imaging (fMRI). On the contrary, the spinal cord has been largely overlooked despite its pivotal role in processing sensorimotor signals. Only a handful of studies have probed the organization of spinal resting-state fluctuations, always using static measures of connectivity. Many innovative approaches have emerged for analyzing dynamics of brain fMRI, but they have not yet been applied to the spinal cord, although they could help disentangle its functional architecture. Here, we leverage a dynamic connectivity method based on the clustering of hemodynamic-informed transients to unravel the rich dynamic organization of spinal resting-state signals. We test this approach in 19 healthy subjects, uncovering fine-grained spinal components and highlighting their neuroanatomical and physiological nature. We provide a versatile tool, the spinal innovation-driven co-activation patterns (SpiCiCAP) framework, to characterize spinal circuits during rest and task, as well as their disruption in neurological disorders.

Assessing the spatial distribution of cervical spinal cord activity during tactile stimulation of the upper extremity in humans with functional magnetic resonance imaging

Dermatomal maps are a mainstay of clinical practice and provide information on the spatial distribution of the cutaneous innervation of spinal nerves. Dermatomal deficits can help isolate the level of spinal nerve root involvement in spinal conditions and guide clinicians in diagnosis and treatment. Dermatomal maps, however, have limitations, and the spatial distribution of spinal cord sensory activity in humans remains to be quantitatively assessed. Here we used spinal cord functional MRI to map and quantitatively compare the spatial distribution of sensory spinal cord activity during tactile stimulation of the left and right lateral shoulders (i.e. C5 dermatome) and dorsal third digits of the hands (i.e., C7 dermatome) in healthy humans (n = 24, age = 36.8 ± 11.8 years). Based on the central sites for processing of innocuous tactile sensory information, we hypothesized that the activity would be localized more to the ipsilateral dorsal spinal cord with the lateral shoulder stimulation activity being localized more superiorly than the dorsal third digit. The findings demonstrate lateralization of the activity with the left- and right-sided stimuli having more activation in the ipsilateral hemicord. Contradictory to our hypotheses, the activity for both stimulation sites was spread across the dorsal and ventral hemicords and did not demonstrate a clear superior-inferior localization. Instead, the activity for both stimuli had a broader than expected distribution, extending across the C5, C6, and C7 spinal cord segments. We highlight the complexity of the human spinal cord neuroanatomy and several sources of variability that may explain the observed patterns of activity. While the findings were not completely consistent with our a priori hypotheses, this study provides a foundation for continued work and is an important step towards developing normative quantitative spinal cord measures of sensory function, which may become useful objective MRI-based biomarkers of neurological injury and improve the management of spinal disorders.

Multiple sclerosis lesions affect intrinsic functional connectivity of the spinal cord

Patients with multiple sclerosis present with focal lesions throughout the spinal cord. There is a clinical need for non-invasive measurements of spinal cord activity and functional organization in multiple sclerosis, given the cord's critical role in the disease. Recent reports of spontaneous blood oxygenation level-dependent fluctuations in the spinal cord using functional MRI suggest that, like the brain, cord activity at rest is organized into distinct, synchronized functional networks among grey matter regions, likely related to motor and sensory systems. Previous studies looking at stimulus-evoked activity in the spinal cord of patients with multiple sclerosis have demonstrated increased levels of activation as well as a more bilateral distribution of activity compared to controls. Functional connectivity studies of brain networks in multiple sclerosis have revealed widespread alterations, which may take on a dynamic trajectory over the course of the disease, with compensatory increases in connectivity followed by decreases associated with structural damage. We build upon this literature by examining functional connectivity in the spinal cord of patients with multiple sclerosis. Using ultra-high field 7 T imaging along with processing strategies for robust spinal cord functional MRI and lesion identification, the present study assessed functional connectivity within cervical cord grey matter of patients with relapsing-remitting multiple sclerosis (n = 22) compared to a large sample of healthy controls (n = 56). Patient anatomical images were rated for lesions by three independent raters, with consensus ratings revealing 19 of 22 patients presented with lesions somewhere in the imaged volume. Linear mixed models were used to assess effects of lesion location on functional connectivity. Analysis in control subjects demonstrated a robust pattern of connectivity among ventral grey matter regions as well as a distinct network among dorsal regions. A gender effect was also observed in controls whereby females demonstrated higher ventral network connectivity. Wilcoxon rank-sum tests detected no differences in average connectivity or power of low frequency fluctuations in patients compared to controls. The presence of lesions was, however, associated with local alterations in connectivity with differential effects depending on columnar location. The patient results suggest that spinal cord functional networks are generally intact in relapsing-remitting multiple sclerosis but that lesions are associated with focal abnormalities in intrinsic connectivity. These findings are discussed in light of the current literature on spinal cord functional MRI and the potential neurological underpinnings.

Intrinsic functional architecture of the non-human primate spinal cord derived from fMRI and electrophysiology

Resting-state functional MRI (rsfMRI) has recently revealed correlated signals in the spinal cord horns of monkeys and humans. However, the interpretation of these rsfMRI correlations as indicators of functional connectivity in the spinal cord remains unclear. Here, we recorded stimulus-evoked and spontaneous spiking activity and local field potentials (LFPs) from monkey spinal cord in order to validate fMRI measures. We found that both BOLD and electrophysiological signals elicited by tactile stimulation co-localized to the ipsilateral dorsal horn. Temporal profiles of stimulus-evoked BOLD signals covaried with LFP and multiunit spiking in a similar way to those observed in the brain. Functional connectivity of dorsal horns exhibited a U-shaped profile along the dorsal-intermediate-ventral axis. Overall, these results suggest that there is an intrinsic functional architecture within the gray matter of a single spinal segment, and that rsfMRI signals at high field directly reflect this underlying spontaneous neuronal activity.

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.

Evidence for a spinal involvement in temporal pain contrast enhancement

Spatiotemporal filtering and amplification of sensory information at multiple levels during the generation of perceptual representations is a fundamental processing principle of the nervous system. While for the visual and auditory system temporal filtering of sensory signals has been noticed for a long time, respective contrast mechanisms within the nociceptive system became only recently subject of investigations, mainly in the context of offset analgesia (OA) subsequent to noxious stimulus decreases. In the present study we corroborate in a first experiment the assumption that offset analgesia involves a central component by showing that an OA-like effect accounting for 74% of a corresponding OA reference can be evoked by decomposing the stimulus offset into two separate box-car stimuli applied within the same dermatome but to separate populations of primary afferent neurons. In order to draw conclusions about the levels of the CNS at which temporal filtering of nociceptive information takes place during OA we investigate in a second experiment neuronal activity in the spinal cord during a painful thermal stimulus offset employing high-resolution fMRI in healthy volunteers. Pain-related BOLD responses in the spinal cord were significantly reduced during OA and their time course followed widely behavioral hypoalgesia, but not the thermal stimulation profile. In summary, the results suggest that temporal pain contrast enhancement during OA comprises a central mechanism and this mechanism becomes already effective at the level of the spinal cord.

Functional Neuroimaging of Nociceptive and Pain-Related Activity in the Spinal Cord and Brain: Insights From Neurovascular Coupling Studies

Spinal cord and brain processes underlie pain perception, which produces systemic cardiovascular changes. In turn, the autonomic nervous system regulates vascular function in the spinal cord and brain in order to adapt to these systemic changes, while neuronal activity induces local vascular changes. Thus, autonomic regulation and pain processes in the brain and spinal cord are tightly linked and interrelated. The objective of this topical review is to discuss work on neurovascular coupling during nociceptive processing in order to highlight supporting evidence and limitations for the use of cerebral and spinal fMRI to investigate pain mechanisms and spinal nociceptive processes. Work on functional neuroimaging of pain is presented and discussed in relation to available neurovascular coupling studies and related issues. Perspectives on future work are also discussed with an emphasis on differences between the brain and the spinal cord and on different approaches that may be useful to improve current methods, data analyses and interpretation. In summary, this review highlights the lack of data on neurovascular coupling during nociceptive stimulation and indicates that hemodynamic and BOLD responses measured with fMRI may be biased by nonspecific vascular changes. Future neuroimaging studies on nociceptive and pain-related processes would gain further understanding of neurovascular coupling in the brain and spinal cord and should take into account the effects of systemic vascular changes that may affect hemodynamic responses. Anat Rec, 301:1585-1595, 2018. © 2018 Wiley Periodicals, Inc.

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.

Altered Signaling in the Descending Pain-modulatory System after Short-Term Infusion of the μ-Opioid Agonist Remifentanil

μ-Opioid receptor agonists are widely used within the contemporary treatment of pain, but abrupt opioid suspension, even after short-term infusion, can paradoxically increase the sensitivity to noxious stimuli, a phenomenon that has been, for example, reported after application of the fast-acting μ-opioid receptor agonist remifentanil. To investigate the mechanisms underlying the effects of discontinuation of remifentanil application on pain processing in the human CNS, we analyzed neuronal responses to thermal stimuli before and after a short-term infusion of remifentanil (30 min 0.1 μg/kg body weight/min) compared with control in the brain, brainstem, and spinal cord in drug-naive male volunteers using fMRI. Subsequent to remifentanil suspension, we observed reduced heat pain thresholds and increased neuronal responses in pain-encoding as well as in key regions of the descending pain-modulatory system, such as the periaqueductal gray matter, the nucleus cuneiformis, and the rostral ventromedial medulla. Moreover, the spinal pain-related multivoxel activity pattern showed an opioid-specific change after drug suspension. Importantly, remifentanil suspension increased the functional coupling between the nucleus cuneiformis and the rostral anterior cingulate cortex, and the coupling strength between the rostral anterior cingulate cortex and the nucleus cuneiformis correlated negatively with the individual pain threshold after opioid suspension. These findings demonstrate that, already subsequent to a short-term infusion of the μ-opioid receptor agonist remifentanil, signaling in the descending pain-modulatory system is fundamentally altered and that these changes are directly related to the behavioral sensitivity to pain.SIGNIFICANCE STATEMENT Opioids are widely used in modern medicine, but, in addition to their known side effects, it is increasingly recognized that opioids can also increase sensitivity to pain subsequent to their use. Using the fast-acting μ-opioid receptor agonist remifentanil and fMRI in healthy male volunteers, this study demonstrates how signaling changes occur along the entire descending pain-modulatory pathway after opioid discontinuation and how these alterations are closely linked to increased behavioral pain sensitivity. Particularly by revealing modified responses in pain-modulatory brainstem regions that have been previously demonstrated to be causally involved in acute opioid withdrawal effects in rodents, the data provide a plausible neuronal mechanism by which the increased sensitivity to pain after opioid suspension is mediated in humans.

Thermal Stimulation Alters Cervical Spinal Cord Functional Connectivity in Humans

The spinal cord has an active role in the modulation and transmission of the neural signals traveling between the body and the brain. Recent advancements in functional magnetic resonance imaging (fMRI) have made the in vivo examination of spinal cord function in humans now possible. This technology has been recently extended to the investigation of resting state functional networks in the spinal cord, leading to the identification of distinct patterns of spinal cord functional connectivity. In this study, we expand on the previous work and further investigate resting state cervical spinal cord functional connectivity in healthy participants (n = 15) using high resolution imaging coupled with both seed-based functional connectivity analyses and graph theory-based metrics. Within spinal cord segment functional connectivity was present between the left and right ventral horns (bilateral motor network), left and right dorsal horns (bilateral sensory network), and the ipsilateral ventral and dorsal horns (unilateral sensory-motor network). Functional connectivity between the spinal cord segments was less apparent with the connectivity centered at the region of interest and spanning <20 mm along the superior-inferior axis. In a subset of participants (n = 10), the cervical spinal cord functional network was demonstrated to be state-dependent as thermal stimulation of the right ventrolateral forearm resulted in significant disruption of the bilateral sensory network, increased network global efficiency, and decreased network modularity.