Assessment of data acquisition parameters, and analysis techniques for noise reduction in spinal cord fMRI data

Reproducible Spinal Cord Quantitative MRI Analysis with the Spinal Cord Toolbox

The spinal cord plays a pivotal role in the central nervous system, providing communication between the brain and the body and containing critical motor and sensory networks. Recent advancements in spinal cord MRI data acquisition and image analysis have shown a potential to improve the diagnostics, prognosis, and management of a variety of pathological conditions. In this review, we first discuss the significance of standardized spinal cord MRI acquisition protocol in multi-center and multi-manufacturer studies. Then, we cover open-access spinal cord MRI datasets, which are important for reproducible science and validation of new methods. Finally, we elaborate on the recent advances in spinal cord MRI data analysis techniques implemented in the open-source software package Spinal Cord Toolbox (SCT).

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.

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.

A second-order and slice-specific linear shimming technique to improve spinal cord fMRI

PURPOSE

To develop a second-order and slice-specific linear shimming technique and investigate its efficiency in the mitigation of signal loss and distortions, and the increase of temporal signal-to-noise ratio (tSNR) within the spinal cord during functional Magnetic Resonance Imaging (fMRI) of the human cervical spinal cord.

METHODS

All scans were performed on a General Electric Discovery MR750 3 T scanner, using a head, neck and spine coil and a neurovascular array. To improve B₀ homogeneity, a field map was acquired, and second-order shims (SOS) were optimized over manually defined regions of interest (ROIs). Signal loss from dephasing by susceptibility-induced gradients was reduced by optimizing slice-specific x-, y- and z-shims to maximize signal within the spinal cord. Spectral-spatial excitation pulses were used in both the slice-specific linear shimming calibration scan and fMRI acquisitions. The shimming technique's efficiency was initially tested on eight healthy volunteers by comparing tSNR between images acquired with the manufacturer's standard linear shimming and with our SOS and xyz-shimming technique. Subsequently, using an increased spatial resolution as needed for fMRI of the spinal cord, tSNR measurements were performed on resting-state fMRI images from 14 healthy participants.

RESULTS

Spinal fMRI images acquired with only the standard linear shimming suffered from severe signal loss below the C5 vertebral level. The developed shimming technique compensated for this loss especially at levels C6 and C7, while tSNR was significantly higher at all vertebral levels with SOS and xyz-shimming than without it.

CONCLUSION

A comprehensive shimming approach which includes the use of spectral-spatial excitation pulses along with both second-order and slice-specific linear shim optimization reduces regional signal loss and increases tSNR along the c-spine (C3-C7), improving the ability to record functional signals from the human spinal cord.

Distinct neural signaling characteristics between fibromyalgia and provoked vestibulodynia revealed by means of functional magnetic resonance imaging in the brainstem and spinal cord

Introduction

Fibromyalgia and provoked vestibulodynia are two chronic pain conditions that disproportionately affect women. The mechanisms underlying the pain in these conditions are still poorly understood, but there is speculation that both may be linked to altered central sensitization and autonomic regulation. Neuroimaging studies of these conditions focusing on the brainstem and spinal cord to explore changes in pain regulation and autonomic regulation are emerging, but none to date have directly compared pain and autonomic regulation in these conditions. This study compares groups of women with fibromyalgia and provoked vestibulodynia to healthy controls using a threat/safety paradigm with a predictable noxious heat stimulus.

Methods

Functional magnetic resonance imaging data were acquired at 3 tesla in the cervical spinal cord and brainstem with previously established methods. Imaging data were analyzed with structural equation modeling and ANCOVA methods during: a period of noxious stimulation, and a period before the stimulation when participants were expecting the upcoming pain.

Results

The results demonstrate several similarities and differences between brainstem/spinal cord connectivity related to autonomic and pain regulatory networks across the three groups in both time periods.

Discussion

Based on the regions and connections involved in the differences, the altered pain processing in fibromyalgia appears to be related to changes in how autonomic and pain regulation networks are integrated, whereas altered pain processing in provoked vestibulodynia is linked in part to changes in arousal or salience networks as well as changes in affective components of pain regulation.

Proof-of-concept of a novel structural equation modelling approach for the analysis of functional magnetic resonance imaging data applied to investigate individual differences in human pain responses

A novel network analysis method is demonstrated for applications with functional magnetic resonance imaging (fMRI) data. The method is based on structural equation modeling (SEM) plus modeling of physiological responses in order to explain blood oxygenation-level dependent (BOLD) responses across interconnected regions. The method, termed structural and physiological modeling (SAPM) aims to overcome a weakness of previous analysis methods by estimating both input and output signaling of every region of a network. The results also provide weighting factors (B) which describe the influence of each input signal to a region on its output signaling to another region. The SAPM method is demonstrated by applying it to fMRI data from the brainstem and spinal cord in 55 healthy participants undergoing repeated applications of a heat pain stimulation paradigm. Data are also analyzed using our established SEM method for comparison. The results with both methods indicate that individual differences in nociceptive processing are mediated by differences in descending regulation of spinal cord neurons under the influence of both the nucleus tractus solitarius and periaqueductal gray region. The SAPM results show that BOLD responses in the entire network can be explained during all periods of the stimulation paradigm based on two latent (unobserved) input signaling sources, and a model of the predicted BOLD responses to the heat stimulus. The results demonstrate the concept of our novel SAPM method and provide evidence for its validity. Additional studies are needed to further develop the method and its applications to investigations of complex neural processes across networks.

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

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

An Investigation of Descending Pain Modulation in Women With Provoked Vestibulodynia (PVD): Alterations of Spinal Cord and Brainstem Connectivity

The most common subtype of vulvodynia (idiopathic chronic vulvar pain) is provoked vestibulodynia (PVD). Previous imaging studies have shown that women with vulvodynia exhibit increased neural activity in pain-related brain regions (e.g., the secondary somatosensory cortex, insula, dorsal midcingulate, posterior cingulate, and thalamus). However, despite the recognized role of the spinal cord/brainstem in pain modulation, no previous neuroimaging studies of vulvodynia have examined the spinal cord/brainstem. Sixteen women with PVD and sixteen matched Control women underwent a spinal cord/brainstem functional magnetic resonance imaging (fMRI) session consisting of five runs with no painful thermal stimuli (No Pain), interleaved randomly with five runs with calibrated, moderately painful heat stimulation (Pain). Functional connectivity was also assessed in periods before, during, and after, pain stimulation to investigate dynamic variations in pain processing throughout the stimulation paradigm. Functional connectivity in the brainstem and spinal cord for each group was examined using structural equation modeling (SEM) for both Pain and No Pain conditions. Significant connectivity differences during stimulation were identified between PVD and Control groups within pain modulatory regions. Comparisons of Pain and No Pain conditions identified a larger number of connections in the Control group than in the PVD group, both before and during stimulation. The results suggest that women with PVD exhibit altered pain processing and indicate an insufficient response of the pain modulation system. This study is the first to examine the spinal cord/brainstem functional connectivity in women with PVD, and it demonstrates altered connectivity related to pain modulation in the spinal cord/brainstem.

Spinal cord neural activity of patients with fibromyalgia and healthy controls during temporal summation of pain: an fMRI study

The cause for the increased sensitivity of patients with fibromyalgia (FM) to painful stimuli is unclear but sensitization of dorsal horn spinal cord neurons has been suggested. There, critical changes of sensory information occur which depend on the plasticity of second-order neurons and descending pain modulation, including facilitation and inhibition. This study used repetitive stimuli that produce temporal-summation-of-second-pain (TSSP) and central sensitization, relevant mechanisms for patients with chronic pain. We examined spinal cord neural activation during TSSP in patients with FM and healthy controls (HC) and used its functional connectivity with several brainstem nuclei to model the observed blood-oxygen-level-dependent (BOLD) time-course with pain ratings. Sixteen HC and 14 FM participants received repetitive heat stimuli to the hand at 0.4 Hz to achieve TSSP during functional imaging with a 3 T-Philips Achieva MRI scanner. Stimuli were adjusted to each individual's pain sensitivity to achieve maximal pain ratings of 50 ± 10 on a numerical pain scale (0-100). Using a 16-channel neurovascular coil, multiple image series were obtained from the cervical spinal cord to the brainstem using single-shot turbo-spin echo sequences. During repetitive, sensitivity-adjusted heat stimuli, pain ratings of all subjects increased as predicted, consistent with TSSP. HC and FM participants had similar temporal patterns of spinal activation: initial BOLD increase followed by deactivation. Structural equation modeling showed that the observed spinal activity during TSSP was associated with more BOLD activity across/within the brainstem in FM subjects than HC, suggesting differences in pain modulation.NEW & NOTEWORTHY "Windup" and its behavioral correlate "temporal-summation-of-second pain" (TSSP) represent spinal cord mechanisms of pain augmentation associated with central sensitization and chronic pain. Fibromyalgia (FM) is a chronic pain disorder, where abnormal TSSP has been demonstrated. We used fMRI to study spinal cord and brainstem activation during TSSP. We characterized the time course of spinal cord and brainstem BOLD activity during TSSP which showed abnormal brainstem activity in patients with FM, possibly due to deficient pain modulation.

Investigation of the neural basis of expectation-based analgesia in the human brainstem and spinal cord by means of functional magnetic resonance imaging

•

Expectation of lower pain results in lower perceived pain in healthy humans.

•

This expectation analgesia is mediated by descending regulation of the spinal cord.

•

Connectivity analyses showed effects of expecting lower pain prior to stimulation.

•

Expectation analgesia involves regions linked to arousal and autonomic regulation.

Purpose

The expected intensity of pain resulting from a noxious stimulus has been observed to have a strong influence on the pain that is perceived. The neural basis of pain reduction, as a result of expecting lower pain, was investigated using functional magnetic resonance imaging (fMRI) in the brainstem and spinal cord.

Methods

Functional MRI studies were carried out in a region spanning the brainstem and cervical spinal cord in healthy participants. Participants were familiarized with a noxious heat stimulus and study procedures in advance, and were informed during each trial that either a heat calibrated to produce moderate pain (Base state), or a temperature 1 °C lower (Low state), would be applied to their hand. However, the Base temperature was applied in every trial.

Results

Pain ratings were significantly reduced as a result of expecting lower temperatures. FMRI results demonstrate blood oxygenation-level dependent (BOLD) signal variations in response to participants being informed of the stimulus to expect, in advance of stimulation, and in response to stimulation. Significant coordination of BOLD signals was also detected across specific brainstem and spinal cord regions, with connectivity strengths that varied significantly with the study condition, and with individual pain ratings. The results identify regions that are known to be involved with arousal and autonomic regulation.

Conclusions

Expectation-based analgesia is mediated by descending regulation of spinal cord nociceptive responses. This regulation appears to be related to arousal and autonomic regulation, consistent with the cognitive/affective dimension of pain.

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.

A comparison of the effectiveness of functional MRI analysis methods for pain research: The new normal

Studies of the neural basis of human pain processing present many challenges because of the subjective and variable nature of pain, and the inaccessibility of the central nervous system. Neuroimaging methods, such as functional magnetic resonance imaging (fMRI), have provided the ability to investigate these neural processes, and yet commonly used analysis methods may not be optimally adapted for studies of pain. Here we present a comparison of model-driven and data-driven analysis methods, specifically for the study of human pain processing. Methods are tested using data from healthy control participants in two previous studies, with separate data sets spanning the brain, and the brainstem and spinal cord. Data are analyzed by fitting time-series responses to predicted BOLD responses in order to identify significantly responding regions (model-driven), as well as with connectivity analyses (data-driven) based on temporal correlations between responses in spatially separated regions, and with connectivity analyses based on structural equation modeling, allowing for multiple source regions to explain the signal variations in each target region. The results are assessed in terms of the amount of signal variance that can be explained in each region, and in terms of the regions and connections that are identified as having BOLD responses of interest. The characteristics of BOLD responses in identified regions are also investigated. The results demonstrate that data-driven approaches are more effective than model-driven approaches for fMRI studies of pain.

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.

Coordinated Human Brainstem and Spinal Cord Networks during the Expectation of Pain Have Elements Unique from Resting-State Effects

Functional magnetic resonance imaging (fMRI) research on the human brainstem (BS) and spinal cord (SC) has identified extensive BS/SC resting-state networks (RSNs) by showing spontaneous coordinated blood oxygenation-level dependent (BOLD) signal fluctuations in the absence of a stimulus. Studies have shown that these networks can be influenced by participants’ level of arousal or attention (e.g., watching a video), and linked network function to autonomic homeostatic regulation. Here we explore how the cognitive state of expecting pain can influence connectivity in these networks. Data from two studies (a predictable pain stimulus study, and a resting-state study) were compared to show the effects of expecting pain on BS/SC networks, and how networks differed from networks associated with the resting-state. In each study, BOLD fMRI data were obtained from the cervical SC and brainstem in healthy participants at 3 tesla using a T₂-weighted single-shot fast spin-echo imaging method. Functional connectivity was investigated within the entire 3D volume by means of structural equation modeling (SEM) and analyses of covariance (ANCOVA). Results showed extensive connectivity within/across BS and SC regions during the expectation of pain, and ANCOVA analyses showed that connectivity in specific components of these networks varied with individual pain sensitivity. Comparing these results to RSN fluctuations revealed commonalities in coordination between BS and SC regions, and specific BS–BS connectivity fluctuations unique to the expectation of pain. Based on the regions involved, these results provide evidence of brainstem regulation specific to the expectation of pain.

What Has Been Learned from Magnetic Resonance Imaging Examination of the Injured Human Spinal Cord: A Canadian Perspective

Magnetic resonance imaging (MRI) has transformed the way surgeons and researchers study and treat spinal cord injury. In this narrative review, we explore the historical context of imaging the human spinal cord and describe how MRI has evolved from providing the first visualization of the human spinal cord in the 1980s to a remarkable set of imaging tools today. The article focuses in particular on the role of Canadian researchers to this field. We begin by outlining the clinical context of traumatic injury to the human spinal cord and describe why current MRI standards fall short when it comes to treating this disabling condition. Parts 2 and 3 of this work explore an exciting and dramatic shift in the use of MRI technology to aid in our understanding and treatment of traumatic injury to the spinal cord. We explore the use of functional imaging (part 2) and structural imaging (part 3) and explore how these techniques have evolved, how they are used, and the challenges that we face for continued refinement and application to patients who live with the neurological and functional deficits caused by injury to the delicate spinal cord.

Lumbar Spinal Cord Activity and Blood Biochemical Changes in Individuals With Diabetic Peripheral Neuropathy During Electrical Stimulation

It is difficult to perform an in vivo evaluation of the nerve conduction mechanism in a patient with diabetic peripheral neuropathy (DPN). We aim to explore possible activation differences to enable a further understanding of the nerve conduction mechanisms of diabetic neuropathy and to present a novel clinical method to evaluate nerve injury and recovery. DPN patients (n = 20) and healthy volunteers (n = 20) were included in this study to detect the functional activation of the lumbar spinal cord via electric stimulation. Spinal fMRI data sets were acquired via a single-shot fast spin echo (SSFSE) sequence. A task-related fMRI was performed via low-frequency electrical stimulation. After post-processing, the active voxels and the percentage of signal changes were calculated for the DPN evaluation and the correlations between the blood biochemical indexes, such as glucose, total cholesterol, and hemoglobin A1c were explored. Activation in the DPN patients was primarily observed in the T12 (10/13) vertebral level. The percentage of signal changes in DPN patients was higher than that in the control group (Z = -2.757, P < 0.05). Positive correlation between the percentage of signal changes and the total cholesterol/glucose in the DNP group was found (P < 0.05). Lumbar spinal cord fMRI, based on the SEEP effect, was determined to be feasible. The repetitive activation distribution was primarily located at the T12 vertebral level. Lumbar spinal cord fMRI might be used as a potential tool to assess and reveal the nerve conduction mechanisms in DPN.

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.

Lateralized Brainstem and Cervical Spinal Cord Responses to Aversive Sounds: A Spinal fMRI Study

Previous research has delineated the networks of brain structures involved in the perception of emotional auditory stimuli. These include the amygdala, insula, and auditory cortices, as well as frontal-lobe, basal ganglia, and cerebellar structures involved in the planning and execution of motoric behaviors. The aim of the current research was to examine whether emotional sounds also influence activity in the brainstem and cervical spinal cord. Seventeen undergraduate participants completed a spinal functional magnetic resonance imaging (fMRI) study consisting of two fMRI runs. One run consisted of three one-minute blocks of aversive sounds taken from the International Affective Digitized Sounds (IADS) stimulus set; these blocks were interleaved by 40-s rest periods. The other block consisted of emotionally neutral stimuli also drawn from the IADS. The results indicated a stark pattern of lateralization. Aversive sounds elicited greater activity than neutral sounds in the right midbrain and brainstem, and in right dorsal and ventral regions of the cervical spinal cord. Neutral stimuli, on the other hand, elicited less neural activity than aversive sounds overall; these responses were left lateralized and were found in the medial midbrain and the dorsal sensory regions of the cervical spinal cord. Together, these results demonstrate that aversive auditory stimuli elicit increased sensorimotor responses in brainstem and cervical spinal cord structures.

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.

Continuous Descending Modulation of the Spinal Cord Revealed by Functional MRI

Spontaneous variations in spinal cord activity may arise from regulation of any of a number of functions including sensory, motor, and autonomic control. Here, we use functional MRI (fMRI) of healthy participants to identify properties of blood oxygenation-level dependent (BOLD) variations in the spinal cord in response to knowledge that either a noxious stimulus is impending, or that no stimulus is to be expected. Expectation of a noxious stimulus, or no stimulus, is shown to have a significant effect on wide-spread BOLD signal variations in the spinal cord over the entire time period of the fMRI acquisition. Coordination of BOLD responses between/within spinal cord and brainstem regions are also influenced by this knowledge. We provide evidence that such signal variations are the result of continuous descending modulation of spinal cord function. BOLD signal variations in response to noxious stimulation of the hand are also shown, as in previous studies. The observation of both continuous and reactive BOLD responses to emotional/cognitive factors and noxious peripheral stimulation may have important implications, not only for our understanding of endogenous pain modulation, but also in showing that spinal cord activity is under continuous regulatory control.

Denoising spinal cord fMRI data: Approaches to acquisition and analysis

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

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

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

Changes in Pain Processing in the Spinal Cord and Brainstem after Spinal Cord Injury Characterized by Functional Magnetic Resonance Imaging

Traumatic spinal cord injury (SCI) has a number of devastating consequences, including high prevalence of chronic pain and altered pain sensitivity. The causes of altered pain states vary depending on the injury and are difficult to diagnose and treat. A better understanding of pain mechanisms after SCI is expected to lead to better diagnostic capabilities and improved treatments. We therefore applied functional magnetic resonance imaging (fMRI) of the brainstem and spinal cord in a group of participants with previous traumatic SCI to characterize changes in pain processing as a result of their injuries. The same thermal stimulus was applied to the medial palm (C8 dermatome) as a series of repeated brief noxious thermal pulses in a group of 16 participants with a cervical (n = 14) and upper thoracic (n = 2) injuries. Functional MRI of the brainstem and spinal cord was used to determine the neuronal activity evoked by the noxious stimulation, and connectivity between regions was characterized with structural equation modeling (SEM). The results show that pain ratings, the location and magnitude of blood oxygenation-level dependent fMRI results, and connectivity assessed with SEM varied widely across participants. However, the results varied in relation to the perceived pain and the level/severity of injuries, particularly in terms of hypothalamus connectivity with other regions, and descending modulation via the periaqueductal gray matter-rostral ventromedial medulla-cord pathway. The results, therefore, appear to provide sensitive indicators of each individual's pain response, and information about the mechanisms of altered pain sensitivity. The ability to characterize changes in pain processing in individuals with SCI represents a significant technological advance.

FMRI of spinal and supra-spinal correlates of temporal pain summation in fibromyalgia patients

Fibromyalgia syndrome (FM) is a debilitating chronic pain condition, which afflicts primarily females. Although the etiology of this illness is not completely understood, FM pain is thought to rely on enhanced pain sensitivity maintained by central mechanisms. One of these mechanisms is central pain amplification, which is characterized by altered temporal summation of second pain (TSSP). Here we use a TSSP paradigm and functional MRI (fMRI) of the spinal cord, brainstem, and brain to noninvasively examine the central nervous system contributions to TSSP in FM patients and normal controls (NC). Functional MRI of pain-free female adults (N = 15) and FM patients (N = 14) was conducted while brief, repetitive heat pain stimuli (0.33 Hz) were applied to the thenar eminence of the hand (C6 dermatome). The stimulus intensity was adjusted to each participant's heat pain sensitivity to achieve moderate pain. Data were analyzed by means of a General Linear Model and region-of-interest analyses. All participants demonstrated significant pain summation in the TSSP condition. FM subjects, however, required significantly lower stimulus intensities than NC to achieve similar TSSP. fMRI analyses of perceptually equal TSSP identified similar brain activity in NC and FM subjects; however, multiple areas in the brainstem (rostral ventromedial medulla and periaqueductal grey region) and spinal cord (dorsal horn) exhibited greater activity in NC subjects. Finally, increased after-sensations and enhanced dorsal horn activity was demonstrated in FM patients. In conclusion, the spinal and brainstem BOLD responses to TSSP are different between NC and FM patients, which may indicate alterations to descending pain control mechanisms suggesting contributions of these mechanisms to central sensitization and pain of FM patients.

Assessing Nociception by fMRI of the Human Spinal Cord: A Systematic Review

OBJECTIVE

To assess the use of fMRI of the spinal cord in measuring noxious stimulation.

METHODS

The Scopus, Medline, EMBASE, and Web of Science databases were searched, along with the reference lists of included articles. Two independent reviewers screened abstracts, full-text articles, and extracted data. Original research was included if fMRI of the human spinal cord was used to measure responses to noxious stimulation.

RESULTS

Of the 192 abstracts screened, 19 met the search criteria and were divided according to their focus: investigating pain responses (n = 6), methodology (n = 6), spinal cord injury (n = 2), or cognition–pain interactions (n = 5). All but one study appear to have observed activity in ipsilateral and dorsal gray matter regions in response to noxious stimuli, although contralateral or ventral activity was also widely observed.

CONCLUSIONS

Although nociception can be investigated using spinal fMRI, establishing reliability, standardizing methodology, and reporting of results will greatly advance this field.

Neural correlates of temporal summation of second pain in the human brainstem and spinal cord

Temporal summation of second pain (TSSP) occurs when painful stimuli are presented repetitively (≥ 0.33 Hz) and results from a C-fibre evoked enhancement (or "wind-up") of the dorsal horn neurons. Based on electrophysiological studies in intact animals, windup is considered a purely central phenomenon. With advancements in functional MRI (fMRI), we can now probe the central mechanisms of this pain response in humans. The aim of this study is to characterize the fMRI responses in the healthy human brainstem and spinal cord that correspond to TSSP. Functional MRI of healthy female adults (N = 15) was conducted while brief, repetitive heat pain stimuli were applied to the right thenar eminence (C6 dermatome), and TSSP (0.33 Hz) and control (0.17 Hz) heat pain paradigms were employed. The stimulus intensity was adjusted to each participant's heat pain sensitivity. Data were analyzed by means of a general linear model, and region-of-interest analyses. As predicted, participants demonstrated significant behavioural summation of pain in the TSSP condition. FMRI results identified enhanced activity in the spinal cord dorsal horn at C6 in response to the TSSP condition. Additionally, multiple areas of the brainstem (RVM and PAG) showed greater responses with the TSSP condition. These results suggest that, in humans, increased pain perception in the TSSP condition is reflected by greater responses in the dorsal horn and in regions known to play a role in the descending modulation of pain, which may modulate the spinal cord response.