Multiple sclerosis lesions affect intrinsic functional connectivity of the spinal cord

Spinal Cord Imaging in Multiple Sclerosis and Related Disorders

Spinal cord MRI plays an important role in the diagnosis and prognosis of multiple sclerosis (MS) and related disorders. The ANATOMICAL, pathologic, imaging and prognostic consideriations for the spinal cord for MS and the most important other demyelinating disorders, neuromyelitis optica spectrum disorder and myelin oligodendrocyte glycoprotein-associated disease, are reviewed. Finally, differential diagnostic considerations of spinal cord MRI in MS and related disorders are discussed.

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.

Calibration-free parallel transmission of the cervical, thoracic, and lumbar spinal cord at 7T

PURPOSE

To address the limitations of spinal cord imaging at ultra-high field (UHF) due to time-consuming parallel transmit (pTx) adjustments. This study introduces calibration-free offline computed universal shim modes that can be applied seamlessly for different pTx RF coils and spinal cord target regions, substantially enhancing spinal cord imaging efficiency at UHF.

METHODS

A library of channel-wise relativeB 1 + $$ {B}_1^{+} $$ maps for the cervical spinal cord (six datasets) and thoracic and lumbar spinal cord (nine datasets) was constructed to optimize transmit homogeneity and efficiency for these regions. A tailored B0 shim was optimized for the cervical spine to enhance spatial magnetic field homogeneity further. The performance of the universal shims was validated using absolute saturation basedB 1 + $$ {B}_1^{+} $$ mapping and high-resolution 2D and 3D multi-echo gradient-recalled echo (GRE) data to assess the image quality.

RESULTS

The proposed universal shims demonstrated a 50% improvement inB 1 + $$ {B}_1^{+} $$ efficiency compared to the default (zero phase) shim mode.B 1 + $$ {B}_1^{+} $$ homogeneity was also improved by 20%. The optimized universal shims achieved performance comparable to subject-specific pTx adjustments, while eliminating the need for lengthy pTx calibration times, saving about 10 min per experiment.

CONCLUSION

The development of universal shims represents a significant advance by eliminating time-consuming subject-specific pTx adjustments. This approach is expected to make UHF spinal cord imaging more accessible and user-friendly, particularly for non-pTx experts.

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

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

Cervical spinal cord susceptibility-weighted MRI at 7T: Application to multiple sclerosis

BACKGROUND

Susceptibility-weighted imaging (SWI) has been extensively studied in the brain and in diseases of the central nervous system such as multiple sclerosis (MS) providing unique opportunities to visualize cerebral vasculature and disease-related pathology, including the central vein sign (CVS) and paramagnetic rim lesions (PRLs). However, similar studies evaluating SWI in the spinal cord of patients with MS remain severely limited.

PURPOSE

Based on our previous findings of enlarged spinal vessels in MS compared to healthy controls (HCs), we developed high-field SWI acquisition and processing methods for the cervical spinal cord with application in people with MS (pwMS) and HCs. Here, we demonstrate the vascular variability between the two cohorts and unique MS lesion features in the cervical cord.

METHODS

In this retrospective, exploratory pilot study conducted between March 2021 and March 2022, we scanned 12 HCs and 9 pwMS using an optimized non-contrast 2D T2*-weighted gradient echo sequence at 7 tesla. The overall appearance of the white and gray matter as well as tissue vasculature were compared between the two cohorts and areas of MS pathology in the patient group were assessed using both the magnitude and processed SWI images.

RESULTS

We show improved visibility of vessels and more pronounced gray and white matter contrast in the MS group compared to HCs, hypointensities surrounding the cord in the MS cohort, and identify signal changes indicative of the CVS and paramagnetic rims in 66 % of pwMS with cervical spinal lesions.

CONCLUSION

In this first study of SWI at 7T in the human spinal cord, SWI holds promise in advancing our understanding of disease processes in the cervical cord in MS.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

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.

Longitudinal multiparametric MRI of traumatic spinal cord injury in animal models

Multi-parametric MRI (mpMRI) technology enables non-invasive and quantitative assessments of the structural, molecular, and functional characteristics of various neurological diseases. Despite the recognized importance of studying spinal cord pathology, mpMRI applications in spinal cord research have been somewhat limited, partly due to technical challenges associated with spine imaging. However, advances in imaging techniques and improved image quality now allow longitudinal investigations of a comprehensive range of spinal cord pathological features by exploiting different endogenous MRI contrasts. This review summarizes the use of mpMRI techniques including blood oxygenation level-dependent (BOLD) functional MRI (fMRI), diffusion tensor imaging (DTI), quantitative magnetization transfer (qMT), and chemical exchange saturation transfer (CEST) MRI in monitoring different aspects of spinal cord pathology. These aspects include cyst formation and axonal disruption, demyelination and remyelination, changes in the excitability of spinal grey matter and the integrity of intrinsic functional circuits, and non-specific molecular changes associated with secondary injury and neuroinflammation. These approaches are illustrated with reference to a nonhuman primate (NHP) model of traumatic cervical spinal cord injuries (SCI). We highlight the benefits of using NHP SCI models to guide future studies of human spinal cord pathology, and demonstrate how mpMRI can capture distinctive features of spinal cord pathology that were previously inaccessible. Furthermore, the development of mechanism-based MRI biomarkers from mpMRI studies can provide clinically useful imaging indices for understanding the mechanisms by which injured spinal cords progress and repair. These biomarkers can assist in the diagnosis, prognosis, and evaluation of therapies for SCI patients, potentially leading to improved outcomes.

The confound of hemodynamic response function variability in human resting-state functional MRI studies

Functional magnetic resonance imaging (fMRI) is an indirect measure of neural activity with the hemodynamic response function (HRF) coupling it with unmeasured neural activity. The HRF, modulated by several non-neural factors, is variable across brain regions, individuals and populations. Yet, a majority of human resting-state fMRI connectivity studies continue to assume a non-variable HRF. In this article, with supportive prior evidence, we argue that HRF variability cannot be ignored as it substantially confounds within-subject connectivity estimates and between-subjects connectivity group differences. We also discuss its clinical relevance with connectivity impairments confounded by HRF aberrations in several disorders. We present limited data on HRF differences between women and men, which resulted in a 15.4% median error in functional connectivity estimates in a group-level comparison. We also discuss the implications of HRF variability for fMRI studies in the spinal cord. There is a need for more dialogue within the community on the HRF confound, and we hope that our article is a catalyst in the process.

A review on multiple sclerosis prognostic findings from imaging, inflammation, and mental health studies

Magnetic resonance imaging (MRI) of the brain is commonly used to detect where chronic and active lesions are in multiple sclerosis (MS). MRI is also extensively used as a tool to calculate and extrapolate brain health by way of volumetric analysis or advanced imaging techniques. In MS patients, psychiatric symptoms are common comorbidities, with depression being the main one. Even though these symptoms are a major determinant of quality of life in MS, they are often overlooked and undertreated. There has been evidence of bidirectional interactions between the course of MS and comorbid psychiatric symptoms. In order to mitigate disability progression in MS, treating psychiatric comorbidities should be investigated and optimized. New research for the prediction of disease states or phenotypes of disability have advanced, primarily due to new technologies and a better understanding of the aging brain.

Reliability of resting-state functional connectivity in the human spinal cord: assessing the impact of distinct noise sources

The investigation of spontaneous fluctuations of the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord, where it has stimulated interest from a clinical perspective. A number of resting-state functional magnetic resonance imaging (fMRI) studies have demonstrated robust functional connectivity between the time series of BOLD fluctuations in bilateral dorsal horns and between those in bilateral ventral horns, in line with the functional neuroanatomy of the spinal cord. A necessary step prior to extension to clinical studies is assessing the reliability of such resting-state signals, which we aimed to do here in a group of 45 healthy young adults at the clinically prevalent field strength of 3T. When investigating connectivity in the entire cervical spinal cord, we observed fair to good reliability for dorsal-dorsal and ventral-ventral connectivity, whereas reliability was poor for within- and between-hemicord dorsal-ventral connectivity. Considering how prone spinal cord fMRI is to noise, we extensively investigated the impact of distinct noise sources and made two crucial observations: removal of physiological noise led to a reduction in functional connectivity strength and reliability - due to the removal of stable and participant-specific noise patterns - whereas removal of thermal noise considerably increased the detectability of functional connectivity without a clear influence on reliability. Finally, we also assessed connectivity within spinal cord segments and observed that while the pattern of connectivity was similar to that of whole cervical cord, reliability at the level of single segments was consistently poor. Taken together, our results demonstrate the presence of reliable resting-state functional connectivity in the human spinal cord even after thoroughly accounting for physiological and thermal noise, but at the same time urge caution if focal changes in connectivity (e.g. due to segmental lesions) are to be studied, especially in a longitudinal manner.

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.

Functional connectivity in the dorsal network of the cervical spinal cord is correlated with diffusion tensor imaging indices in relapsing-remitting multiple sclerosis

Focal lesions may affect functional connectivity (FC) of the ventral and dorsal networks in the cervical spinal cord of people with relapsing-remitting multiple sclerosis (RRMS). Resting-state FC can be measured using functional MRI (fMRI) at 3T. This study sought to determine whether alterations in FC may be related to the degree of damage in the normal-appearing tissue. Tissue integrity and FC in the cervical spinal cord were assessed with diffusion tensor imaging (DTI) and resting-state fMRI, respectively, in a group of 26 RRMS participants with high cervical lesion load, low disability, and minimally impaired sensorimotor function, and healthy controls. Lower fractional anisotropy (FA) and higher radial diffusivity (RD) were observed in the normal-appearing white matter in the RRMS group relative to controls. Average FC in ventral and dorsal networks was similar between groups. Significant associations were found between higher FC in the dorsal sensory network and several DTI markers of pathology in the normal-appearing tissue. In the normal-appearing grey matter, dorsal FC was positively correlated with axial diffusivity (AD) (r = 0.46, p = 0.020) and mean diffusivity (MD) (r = 0.43, p = 0.032). In the normal-appearing white matter, dorsal FC was negatively correlated with FA (r = -0.43, p = 0.028) and positively correlated with RD (r = 0.49, p = 0.012), AD (r = 0.42, p = 0.037) and MD (r = 0.53, p = 0.006). These results suggest that increased connectivity, while remaining within the normal range, may represent a compensatory mechanism in response to structural damage in support of preserved sensory function in RRMS.

Lateral Corticospinal Tract and Dorsal Column Damage: Predictive Relationships With Motor and Sensory Scores at Discharge From Acute Rehabilitation After Spinal Cord Injury

OBJECTIVE

To determine if lateral corticospinal tract (LCST) integrity demonstrates a significant predictive relationship with future ipsilateral lower extremity motor function (LEMS) and if dorsal column (DC) integrity demonstrates a significant predictive relationship with future light touch (LT) sensory function post spinal cord injury (SCI) at time of discharge from inpatient rehabilitation.

DESIGN

Retrospective analyses of imaging and clinical outcomes.

SETTING

University and academic hospital.

PARTICIPANTS

A total of 151 participants (N=151) with SCI.

INTERVENTIONS

Inpatient rehabilitation.

MAIN OUTCOME MEASURES

LEMS and LT scores at discharge from inpatient rehabilitation.

RESULTS

In 151 participants, right LCST spared tissue demonstrated a significant predictive relationship with right LEMS percentage recovered (β=0.56; 95% confidence interval [CI], 0.37-0.73; R=0.43; P<.001). Left LCST spared tissue demonstrated a significant predictive relationship with left LEMS percentage recovered (β=0.66; 95% CI, 0.50-0.82; R=0.51; P<.001). DC spared tissue demonstrated a significant predictive relationship with LT percentage recovered (β=0.69; 95% CI, 0.52-0.87; R=0.55; P<.001). When subgrouping the participants into motor complete vs incomplete SCI, motor relationships were no longer significant, but the sensory relationship remained significant. Those who had no voluntary motor function but recovered some also had significantly greater LCST spared tissue than those who did not recover motor function.

CONCLUSIONS

LCST demonstrated significant moderate predictive relationships with lower extremity motor function at the time of discharge from inpatient rehabilitation, in an ipsilesional manner. DC integrity demonstrated a significant moderate predictive relationship with recovered function of LT. With further development, these neuroimaging methods might be used to predict potential deficits after SCI and to provide corresponding targeted interventions.

Preliminary Analysis of Brain Footprints in Multiple Sclerosis Females with Detrusor Sphincter Dyssynergia: A Concurrent Urodynamic and Functional Magnetic Resonance Imaging Study

Purpose

This study evaluates the grey and white brain matter characteristics in women with multiple sclerosis (MS) and detrusor sphincter dyssynergia (DSD). Grey matter is assessed via the functional connectivity (FC) of brain regions activated during voiding, using functional magnetic resonance imaging (fMRI). Two white matter tracts involved in bladder function, the anterior thalamic radiation (ATR) and superior longitudinal fasciculus (SLF), were evaluated using diffusion tensor imaging (DTI).

Methods

Twenty-seven women with MS (two groups: no-DSD (n=23) or DSD (n=4)), and eight healthy controls (HC) underwent concurrent urodynamic-fMRI evaluation with four cycles of bladder filling and emptying. A FC similarity measure (FC_sim) was calculated for each subject to express the similarity of individual FC at voiding initiation compared to all FC patterns. ATR and SLF tracts were traced and their fractional anisotropy (FA) and mean diffusivity (MD) were recorded.

Results

Mean FC_sim values were significantly different among the three groups indicating distinct FC patterns; however, no significant difference was found between DSD and no-DSD groups. DSD group showed trends of lower FA and higher MD - indicating loss of coherence - in all tracts compared to HCs, and in the left and right ATR when compared to MS women with neither DSD nor voiding dysfunction (VD), suggesting more damage in these tracts for MS women with DSD.

Conclusions

Women with MS show distinctly different FC patterns compared to HCs. There are trends showing more damage in the ATR in women with MS and DSD compared to those with neither DSD nor VD.

Clinical applications of magnetic resonance imaging based functional and structural connectivity

Advances in computational neuroimaging techniques have expanded the armamentarium of imaging tools available for clinical applications in clinical neuroscience. Non-invasive, in vivo brain MRI structural and functional network mapping has been used to identify therapeutic targets, define eloquent brain regions to preserve, and gain insight into pathological processes and treatments as well as prognostic biomarkers. These tools have the real potential to inform patient-specific treatment strategies. Nevertheless, a realistic appraisal of clinical utility is needed that balances the growing excitement and interest in the field with important limitations associated with these techniques. Quality of the raw data, minutiae of the processing methodology, and the statistical models applied can all impact on the results and their interpretation. A lack of standardization in data acquisition and processing has also resulted in issues with reproducibility. This limitation has had a direct impact on the reliability of these tools and ultimately, confidence in their clinical use. Advances in MRI technology and computational power as well as automation and standardization of processing methods, including machine learning approaches, may help address some of these issues and make these tools more reliable in clinical use. In this review, we will highlight the current clinical uses of MRI connectomics in the diagnosis and treatment of neurological disorders; balancing emerging applications and technologies with limitations of connectivity analytic approaches to present an encompassing and appropriate perspective.

Ultra-high field spinal cord MRI in multiple sclerosis: Where are we standing? A literature review

Magnetic resonance imaging (MRI) is a cornerstone in multiple sclerosis (MS) diagnostics and monitoring. Ultra-high field (UHF) MRI is being increasingly used and becoming more accessible. Due to the small diameter and mobility of the spinal cord, imaging this structure at ultra-high fields poses additional challenges compared to brain imaging. Here we review the potential benefits for the MS field by providing a literature overview of the use UHF spinal cord MRI in MS research and we elaborate on the challenges that are faced. Benefits include increased signal- and contrast-to-noise, enabling for higher spatial resolutions, which can improve MS lesion sensitivity in both the spinal white matter as well as grey matter. Additionally, these benefits can aid imaging of microstructural abnormalities in the spinal cord in MS using advanced MRI techniques like functional imaging, MR spectroscopy and diffusion-based techniques. Technical challenges include increased magnetic field inhomogeneities, distortions from physiological motion and optimalisation of sequences. Approaches including parallel imaging techniques, real time shimming and retrospective compensation of physiological motion are making it increasingly possible to unravel the potential of spinal cord UHF MRI in the context of MS research.

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.

Spontaneous neural synchrony links intrinsic spinal sensory and motor networks during unconsciousness

Non-random functional connectivity during unconsciousness is a defining feature of supraspinal networks. However, its generalizability to intrinsic spinal networks remains incompletely understood. Previously, Barry et al., 2014 used fMRI to reveal bilateral resting state functional connectivity within sensory-dominant and, separately, motor-dominant regions of the spinal cord. Here, we record spike trains from large populations of spinal interneurons in vivo in rats and demonstrate that spontaneous functional connectivity also links sensory- and motor-dominant regions during unconsciousness. The spatiotemporal patterns of connectivity could not be explained by latent afferent activity or by populations of interconnected neurons spiking randomly. We also document connection latencies compatible with mono- and disynaptic interactions and putative excitatory and inhibitory connections. The observed activity is consistent with the hypothesis that salient, experience-dependent patterns of neural transmission introduced during behavior or by injury/disease are reactivated during unconsciousness. Such a spinal replay mechanism could shape circuit-level connectivity and ultimately behavior.

Imaging Considerations in Spinal Cord Evaluation

Multiple diverse pathologies result in the clinical presentation of myelopathy. The preferred way to image the spinal cord depends on clinical history, anatomic site of interest, and patient issues limiting certain imaging modalities. This radiology-focused article discusses pertinent physiological considerations, reviews basic and newer imaging techniques, and examines several distinct disease entities in order to highlight the key role of imaging in the work-up of myelopathy.

Imaging of the Spinal Cord in Multiple Sclerosis: Past, Present, Future

Spinal cord imaging in multiple sclerosis (MS) plays a significant role in diagnosing and tracking disease progression. The spinal cord is one of four key areas of the central nervous system where documenting the dissemination in space in the McDonald criteria for diagnosing MS. Spinal cord lesion load and the severity of cord atrophy are believed to be more relevant to disability than white matter lesions in the brain in different phenotypes of MS. Axonal loss contributes to spinal cord atrophy in MS and its degree correlates with disease severity and prognosis. Therefore, measures of axonal loss are often reliable biomarkers for monitoring disease progression. With recent technical advances, more and more qualitative and quantitative MRI techniques have been investigated in an attempt to provide objective and reliable diagnostic and monitoring biomarkers in MS. In this article, we discuss the role of spinal cord imaging in the diagnosis and prognosis of MS and, additionally, we review various techniques that may improve our understanding of the disease.

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.

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.

Longitudinal assessment of recovery after spinal cord injury with behavioral measures and diffusion, quantitative magnetization transfer and functional magnetic resonance imaging

Spinal cord injuries (SCIs) are a leading cause of disability and can severely impact the quality of life. However, to date, the processes of spontaneous repair of damaged spinal cord remain incompletely understood, partly due to a lack of appropriate longitudinal tracking methods. Noninvasive, multiparametric magnetic resonance imaging (MRI) provides potential biomarkers for the comprehensive evaluation of spontaneous repair after SCI. In this study in rats, a clinically relevant contusion injury was introduced at the lumbar level that impairs both hindlimb motor and sensory functions. Quantitative MRI measurements were acquired at baseline and serially post-SCI for up to 2 wk. The progressions of injury and spontaneous recovery in both white and gray matter were tracked longitudinally using pool-size ratio (PSR) measurements derived from quantitative magnetization transfer (qMT) methods, measurements of water diffusion parameters using diffusion tensor imaging (DTI) and intrasegment functional connectivity derived from resting state functional MRI. Changes in these quantitative imaging measurements were correlated with behavioral readouts. We found (a) a progressive decrease in PSR values within 2 wk post-SCI, indicating a progressive demyelination at the center of the injury that was validated with histological staining, (b) PSR correlated closely with fractional anisotropy and transverse relaxation of free water, but did not show significant correlations with behavioral recovery, and (c) preliminary evidence that SCI induced a decrease in functional connectivity between dorsal horns below the injury site at 24 h. Findings from this study not only confirm the value of qMT and DTI methods for assessing the myelination state of injured spinal cord but indicate that they may also have further implications on whether therapies targeted towards remyelination may be appropriate. Additionally, a better understanding of changes after SCI provides valuable information to guide and assess interventions.

What role should spinal cord MRI take in the future of multiple sclerosis surveillance?

INTRODUCTION

In multiple sclerosis (MS), inflammatory, demyelinating, and neurodegenerative phenomena affect the spinal cord, with detrimental effects on patients' clinical disability. Although spinal cord imaging may be challenging, improvements in MRI technologies have contributed to better evaluate spinal cord involvement in MS.

AREAS COVERED

This review summarizes the current state-of-art of the application of conventional and advanced MRI techniques to evaluate spinal cord damage in MS. Typical features of spinal cord lesions, their role in the diagnostic work-up of suspected MS, their predictive role for subsequent disease course and clinical worsening, and their utility to define treatment response are discussed. The role of spinal cord atrophy and of other advanced MRI techniques to better evaluate the associations between spinal cord abnormalities and the accumulation of clinical disability are also evaluated. Finally, how spinal cord assessment could evolve in the future to improve monitoring of disease progression and treatment effects is examined.

EXPERT OPINION

Spinal cord MRI provides relevant additional information to brain MRI in understanding MS pathophysiology, in allowing an earlier and more accurate diagnosis of MS, and in identifying MS patients at higher risk to develop more severe disability. A future role in monitoring the effects of treatments is also foreseen.

MRI features associated with high likelihood of conversion of radiologically isolated syndrome to multiple sclerosis

Radiologically isolated syndrome (RIS) is the asymptomatic precursor to clinically isolated syndrome, relapsing-remitting multiple sclerosis (MS) or primary progressive MS. RIS is frequently diagnosed when an individual gets an MRI for an unrelated medical issue, such as headache or trauma. Treating RIS patients is controversial, but physicians may be inclined to offer prophylactic treatment for high-risk RIS patients. Identifying imaging and clinical features associated with high likelihood of early clinical conversion may prove helpful to identify a high-risk subset for potential MS therapy. The goal of this paper is to review current literatures to identify imaging and clinical features that predict early (within 5 years) conversion from RIS to MS.

Measurement of T2* in the human spinal cord at 3T

PURPOSE

To measure the transverse relaxation time T₂* in healthy human cervical spinal cord gray matter (GM) and white matter (WM) at 3T.

METHODS

Thirty healthy volunteers were recruited. Axial images were acquired using an averaged multi-echo gradient-echo (mFFE) T₂*-weighted sequence with 5 echoes. We used the signal equation for an mFFE sequence with constant dephasing gradients after each echo to jointly estimate the spin density and T₂* for each voxel.

RESULTS

No global difference in T₂* was observed between all GM (41.3 ± 5.6 ms) and all WM (39.8 ± 5.4 ms). No significant differences were observed between left (43.2 ± 6.8 ms) and right (43.4 ± 5.5 ms) ventral GM, left (38.3 ± 6.1 ms) and right (38.6 ± 6.5 ms) dorsal GM, and left (39.4 ± 5.8 ms) and right (40.3 ± 5.8 ms) lateral WM. However, significant regional differences were observed between ventral (43.4 ± 5.7 ms) and dorsal (38.4 ± 6.0 ms) GM (p < 0.05), as well as between ventral (42.9 ± 6.5 ms) and dorsal (37.9 ± 6.2 ms) WM (p < 0.05). In analyses across slices, inferior T₂* was longer than superior T₂* in GM (44.7 ms vs. 40.1 ms; p < 0.01) and in WM (41.8 ms vs. 35.9 ms; p < 0.01).

CONCLUSIONS

Significant differences in T₂* are observed between ventral and dorsal GM, ventral and dorsal WM, and superior and inferior GM and WM. There is no evidence for bilateral asymmetry in T₂* in the healthy cord. These values of T₂* in the spinal cord are notably lower than most reported values of T₂* in the cortex.

Are Magnetic Resonance Imaging Technologies Crucial to Our Understanding of Spinal Conditions?

Persistent spinal (traumatic and nontraumatic) pain is common and contributes to high societal and personal costs globally. There is an acknowledged urgency for new and interdisciplinary approaches to the condition, and soft tissues, including skeletal muscles, the spinal cord, and the brain, are rightly receiving increased attention as important biological contributors. In reaction to the recent suspicion and questioned value of imaging-based findings, this paper serves to recognize the promise that the technological evolution of imaging techniques, and particularly magnetic resonance imaging, is allowing in characterizing previously less visible morphology. We emphasize the value of quantification and data analysis of several contributors in the biopsychosocial model for understanding spinal pain. Further, we highlight emerging evidence regarding the pathobiology of changes to muscle composition (eg, atrophy, fatty infiltration), as well as advancements in neuroimaging and musculoskeletal imaging techniques (eg, fat-water imaging, functional magnetic resonance imaging, diffusion imaging, magnetization transfer imaging) for these important soft tissues. These noninvasive and objective data sources may complement known prognostic factors of poor recovery, patient self-report, diagnostic tests, and the "-omics" fields. When combined, advanced "big-data" analyses may assist in identifying associations previously not considered. Our clinical commentary is supported by empirical findings that may orient future efforts toward collaborative conversation, hypothesis generation, interdisciplinary research, and translation across a number of health fields. Our emphasis is that magnetic resonance imaging technologies and research are crucial to the advancement of our understanding of the complexities of spinal conditions. J Orthop Sports Phys Ther 2019;49(5):320-329. Epub 26 Mar 2019. doi:10.2519/jospt.2019.8793.

Advances in spinal cord imaging in multiple sclerosis

The spinal cord is frequently affected in multiple sclerosis (MS), causing motor, sensory and autonomic dysfunction. A number of pathological abnormalities, including demyelination and neuroaxonal loss, occur in the MS spinal cord and are studied in vivo with magnetic resonance imaging (MRI). The aim of this review is to summarise and discuss recent advances in spinal cord MRI. Advances in conventional spinal cord MRI include improved identification of MS lesions, recommended spinal cord MRI protocols, enhanced recognition of MRI lesion characteristics that allow MS to be distinguished from other myelopathies, evidence for the role of spinal cord lesions in predicting prognosis and monitoring disease course, and novel post-processing methods to obtain lesion probability maps. The rate of spinal cord atrophy is greater than that of brain atrophy (-1.78% versus -0.5% per year), and reflects neuroaxonal loss in an eloquent site of the central nervous system, suggesting that it can become an important outcome measure in clinical trials, especially in progressive MS. Recent developments allow the calculation of spinal cord atrophy from brain volumetric scans and evaluation of its progression over time with registration-based techniques. Fully automated analysis methods, including segmentation of grey matter and intramedullary lesions, will facilitate the use of spinal cord atrophy in trial designs and observational studies. Advances in quantitative imaging techniques to evaluate neuroaxonal integrity, myelin content, metabolic changes, and functional connectivity, have provided new insights into the mechanisms of damage in MS. Future directions of research and the possible impact of 7T scanners on spinal cord imaging will be discussed.

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.

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.

Heat shock proteins are differentially expressed in brain and spinal cord: implications for multiple sclerosis

Multiple sclerosis (MS) is a chronic neurodegenerative disease characterized by demyelination, inflammation and neurodegeneration throughout the central nervous system. Although spinal cord pathology is an important factor contributing to disease progression, few studies have examined MS lesions in the spinal cord and how they differ from brain lesions. In this study we have compared brain and spinal cord white (WM) and grey (GM) matter from MS and control tissues, focusing on small heat shock proteins (HSPB) and HSP16.2. Western blotting was used to examine protein levels of HSPB1, HSPB5, HSPB6, HSPB8 and HSP16.2 in brain and spinal cord from MS and age-matched non-neurological controls. Immunohistochemistry was used to examine expression of the HSPs in MS spinal cord lesions and controls. Expression levels were quantified using ImageJ. Western blotting revealed significantly higher levels of HSPB1, HSPB6 and HSPB8 in MS and control spinal cord compared to brain tissues. No differences in HSPB5 and HSP16.2 protein levels were observed, although HSPB5 protein levels were higher in brain WM versus GM. In MS spinal cord lesions, increased HSPB1 and HSPB5 expression was observed in astrocytes, and increased neuronal expression of HSP16.2 was observed in normal-appearing GM and type 1 GM lesions. The high constitutive expression of several HSPBs in spinal cord and increased expression of HSPBs and HSP16.2 in MS illustrate differences between brain and spinal cord in health and upon demyelination. Regional differences in HSP expression may reflect differences in astrocyte cytoskeleton composition and influence inflammation, possibly affecting the effectiveness of pharmacological agents.

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.