Resting-state brain and spinal cord networks in humans are functionally integrated

Integrating brainstem and cortical functional architectures

The brainstem is a fundamental component of the central nervous system, yet it is typically excluded from in vivo human brain mapping efforts, precluding a complete understanding of how the brainstem influences cortical function. In this study, we used high-resolution 7-Tesla functional magnetic resonance imaging to derive a functional connectome encompassing cortex and 58 brainstem nuclei spanning the midbrain, pons and medulla. We identified a compact set of integrative hubs in the brainstem with widespread connectivity with cerebral cortex. Patterns of connectivity between brainstem and cerebral cortex manifest as neurophysiological oscillatory rhythms, patterns of cognitive functional specialization and the unimodal-transmodal functional hierarchy. This persistent alignment between cortical functional topographies and brainstem nuclei is shaped by the spatial arrangement of multiple neurotransmitter receptors and transporters. We replicated all findings using 3-Tesla data from the same participants. Collectively, this work demonstrates that multiple organizational features of cortical activity can be traced back to the brainstem.

Dynamic Functional Network Connectivity in Acute Incomplete Cervical Cord Injury Patients and Its Associations With Sensorimotor Dysfunction Measures

BACKGROUND

Dynamic functional network connectivity (dFNC) captures temporal variations in functional connectivity during magnetic resonance imaging acquisition. However, the neural mechanisms driving dFNC alterations in the brain networks of patients with acute incomplete cervical cord injury (AICCI) remain unclear.

METHODS

This study included 16 AICCI patients and 16 healthy controls. Initially, independent component analysis was employed to extract whole-brain independent components from resting-state functional magnetic resonance imaging data. Subsequently, a sliding time window approach, combined with k-means clustering, was used to estimate dFNC states for each participant. Finally, a correlation analysis was conducted to examine the association between sensorimotor dysfunction scores in AICCI patients and the temporal characteristics of dFNC.

RESULTS

Independent component analysis was employed to extract 26 whole-brain independent components. Subsequent dynamic analysis identified 4 distinct connectivity states across the entire cohort. Notably, AICCI patients demonstrated a significant preference for State 3 compared to healthy controls, as evidenced by a higher frequency and longer duration spent in this state. Conversely, State 4 exhibited a reduced frequency and shorter dwell time in AICCI patients. Moreover, correlation analysis revealed a positive association between sensorimotor dysfunction and both the mean dwell time and the fraction of time spent in State 3.

CONCLUSIONS

Patients with AICCI demonstrate abnormal connectivity within dFNC states, and the temporal characteristics of dFNC are associated with sensorimotor dysfunction scores. These findings highlight the potential of dFNC as a sensitive biomarker for detecting network functional changes in AICCI patients, providing valuable insights into the dynamic alterations in brain connectivity related to sensorimotor dysfunction in this population.

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.

Concurrent spinal and brain imaging with optically pumped magnetometers

BACKGROUND

The spinal cord and its interactions with the brain are fundamental for movement control and somatosensation. However, brain and spinal electrophysiology in humans have largely been treated as distinct enterprises, in part due to the relative inaccessibility of the spinal cord. Consequently, there is a dearth of knowledge on human spinal electrophysiology, including the multiple pathologies that affect the spinal cord as well as the brain.

NEW METHOD

Here we exploit recent advances in the development of wearable optically pumped magnetometers (OPMs) which can be flexibly arranged to provide coverage of both the spinal cord and the brain in relatively unconstrained environments. This system for magnetospinoencephalography (MSEG) measures both spinal and cortical signals simultaneously by employing custom-made scanning casts.

RESULTS

We evidence the utility of such a system by recording spinal and cortical evoked responses to median nerve stimulation at the wrist. MSEG revealed early (10 - 15 ms) and late (>20 ms) responses at the spinal cord, in addition to typical cortical evoked responses (i.e., N20).

COMPARISON WITH EXISTING METHODS

Early spinal evoked responses detected were in line with conventional somatosensory evoked potential recordings.

CONCLUSION

This MSEG system demonstrates the novel ability for concurrent non-invasive millisecond imaging of brain and spinal cord.

Functional ultrasound imaging of the human spinal cord

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

Toward a brighter constellation: multiorgan neuroimaging of neural and vascular dynamics in the spinal cord and brain

Significance

Pain comprises a complex interaction between motor action and somatosensation that is dependent on dynamic interactions between the brain and spinal cord. This makes understanding pain particularly challenging as it involves rich interactions between many circuits (e.g., neural and vascular) and signaling cascades throughout the body. As such, experimentation on a single region may lead to an incomplete and potentially incorrect understanding of crucial underlying mechanisms.

Aim

We aimed to develop and validate tools to enable detailed and extended observation of neural and vascular activity in the brain and spinal cord. The first key set of innovations was targeted to developing novel imaging hardware that addresses the many challenges of multisite imaging. The second key set of innovations was targeted to enabling bioluminescent (BL) imaging, as this approach can address limitations of fluorescent microscopy including photobleaching, phototoxicity, and decreased resolution due to scattering of excitation signals.

Approach

We designed 3D-printed brain and spinal cord implants to enable effective surgical implantations and optical access with wearable miniscopes or an open window (e.g., for one- or two-photon microscopy or optogenetic stimulation). We also tested the viability for BL imaging and developed a novel modified miniscope optimized for these signals (BLmini).

Results

We describe "universal" implants for acute and chronic simultaneous brain-spinal cord imaging and optical stimulation. We further describe successful imaging of BL signals in both foci and a new miniscope, the "BLmini," which has reduced weight, cost, and form-factor relative to standard wearable miniscopes.

Conclusions

The combination of 3D-printed implants, advanced imaging tools, and bioluminescence imaging techniques offers a coalition of methods for understanding spinal cord-brain interactions. Our work has the potential for use in future research into neuropathic pain and other sensory disorders and motor behavior.

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.

Spinal fMRI demonstrates segmental organisation of functionally connected networks in the cervical spinal cord: A test-retest reliability study

Resting functional magnetic resonance imaging (fMRI) studies have identified intrinsic spinal cord activity, which forms organised motor (ventral) and sensory (dorsal) resting-state networks. However, to facilitate the use of spinal fMRI in, for example, clinical studies, it is crucial to first assess the reliability of the method, particularly given the unique anatomical, physiological, and methodological challenges associated with acquiring the data. Here, we characterise functional connectivity relationships in the cervical cord and assess their between-session test-retest reliability in 23 young healthy volunteers. Resting-state networks were estimated in two ways (1) by estimating seed-to-voxel connectivity maps and (2) by calculating seed-to-seed correlations. Seed regions corresponded to the four grey matter horns (ventral/dorsal and left/right) of C5-C8 segmental levels. Test-retest reliability was assessed using the intraclass correlation coefficient. Spatial overlap of clusters derived from seed-to-voxel analysis between sessions was examined using Dice coefficients. Following seed-to-voxel analysis, we observed distinct unilateral dorsal and ventral organisation of cervical spinal resting-state networks that was largely confined in the rostro-caudal extent to each spinal segmental level, with more sparse connections observed between segments. Additionally, strongest correlations were observed between within-segment ipsilateral dorsal-ventral connections, followed by within-segment dorso-dorsal and ventro-ventral connections. Test-retest reliability of these networks was mixed. Reliability was poor when assessed on a voxelwise level, with more promising indications of reliability when examining the average signal within clusters. Reliability of correlation strength between seeds was highly variable, with the highest reliability achieved in ipsilateral dorsal-ventral and dorso-dorsal/ventro-ventral connectivity. However, the spatial overlap of networks between sessions was excellent. We demonstrate that while test-retest reliability of cervical spinal resting-state networks is mixed, their spatial extent is similar across sessions, suggesting that these networks are characterised by a consistent spatial representation over time.

FASB: an integrated processing pipeline for Functional Analysis of simultaneous Spinal cord-Brain fMRI

Simultaneous functional magnetic resonance imaging (fMRI) of the spinal cord and brain represents a powerful method for examining both ascending sensory and descending motor pathways in humans in vivo . However, its image acquisition protocols, and processing pipeline are less well established. This limitation is mainly due to technical difficulties related to spinal cord fMRI, and problems with the logistics stemming from a large field of view covering both brain and cervical cord. Here, we propose an acquisition protocol optimized for both anatomical and functional images, as well as an optimized integrated image processing pipeline, which consists of a novel approach for automatic modeling and mitigating the negative impact of spinal voxels with low temporal signal to noise ratio (tSNR). We validate our integrated pipeline, named FASB, using simultaneous fMRI data acquired during the performance of a motor task, as well as during resting-state conditions. We demonstrate that FASB outperforms the current spinal fMRI processing methods in three domains, including motion correction, registration to the spinal cord template, and improved detection power of the group-level analysis by removing the effects of participant-specific low tSNR voxels, typically observed at the disk level. Using FASB, we identify significant task-based activations in the expected sensorimotor network associated with a unilateral handgrip force production task across the entire central nervous system, including the contralateral sensorimotor cortex, thalamus, striatum, cerebellum, brainstem, as well as ipsilateral ventral horn at C5-C8 cervical levels. Additionally, our results show significant task-based functional connectivity between the key sensory and motor brain areas and the dorsal and ventral horns of the cervical cord. Overall, our proposed acquisition protocol and processing pipeline provide a robust method for characterizing the activation and functional connectivity of distinct cortical, subcortical, brainstem and spinal cord regions in humans.

Brain MRI changes in degenerative cervical myelopathy: a systematic review

BACKGROUND

Degenerative cervical myelopathy (DCM) is the most common cause of adult spinal cord dysfunction globally. Associated neurological symptoms and signs have historically been explained by pathobiology within the cervical spine. However, recent advances in imaging have shed light on numerous brain changes in patients with DCM, and it is hypothesised that these changes contribute to DCM pathogenesis. The aetiology, significance, and distribution of these supraspinal changes is currently unknown. The objective was therefore to synthesise all current evidence on brain changes in DCM.

METHODS

A systematic review was performed. Cross-sectional and longitudinal studies with magnetic resonance imaging on a cohort of patients with DCM were eligible. PRISMA guidelines were followed. MEDLINE and Embase were searched to 28th August 2023. Duplicate title/abstract screening, data extraction and risk of bias assessments were conducted. A qualitative synthesis of the literature is presented as per the Synthesis Without Meta-Analysis (SWiM) reporting guideline. The review was registered with PROSPERO (ID: CRD42022298538).

FINDINGS

Of the 2014 studies that were screened, 47 studies were identified that used MRI to investigate brain changes in DCM. In total, 1500 patients with DCM were included in the synthesis, with a mean age of 53 years. Brain alterations on MRI were associated with DCM both before and after surgery, particularly within the sensorimotor network, visual network, default mode network, thalamus and cerebellum. Associations were commonly reported between brain MRI alterations and clinical measures, particularly the Japanese orthopaedic association (JOA) score. Risk of bias of included studies was low to moderate.

INTERPRETATION

The rapidly expanding literature provides mounting evidence for brain changes in DCM. We have identified key structures and pathways that are altered, although there remains uncertainty regarding the directionality and clinical significance of these changes. Future studies with greater sample sizes, more detailed phenotyping and longer follow-up are now needed.

FUNDING

ODM is supported by an Academic Clinical Fellowship at the University of Cambridge. BMD is supported by an NIHR Clinical Doctoral Fellowship at the University of Cambridge (NIHR300696). VFJN is supported by an NIHR Rosetrees Trust Advanced Fellowship (NIHR302544). This project was supported by an award from the Rosetrees Foundation with the Storygate Trust (A2844).

Towards a Brighter Constellation: Multi-Organ Neuroimaging of Neural and Vascular Dynamics in the Spinal Cord and Brain

Significance

Pain is comprised of a complex interaction between motor action and somatosensation that is dependent on dynamic interactions between the brain and spinal cord. This makes understanding pain particularly challenging as it involves rich interactions between many circuits (e.g., neural and vascular) and signaling cascades throughout the body. As such, experimentation on a single region may lead to an incomplete and potentially incorrect understanding of crucial underlying mechanisms.

Aim

Here, we aimed to develop and validate new tools to enable detailed and extended observation of neural and vascular activity in the brain and spinal cord. The first key set of innovations were targeted to developing novel imaging hardware that addresses the many challenges of multi-site imaging. The second key set of innovations were targeted to enabling bioluminescent imaging, as this approach can address limitations of fluorescent microscopy including photobleaching, phototoxicity and decreased resolution due to scattering of excitation signals.

Approach

We designed 3D-printed brain and spinal cord implants to enable effective surgical implantations and optical access with wearable miniscopes or an open window (e.g., for one- or two-photon microscopy or optogenetic stimulation). We also tested the viability for bioluminescent imaging, and developed a novel modified miniscope optimized for these signals (BLmini).

Results

Here, we describe novel 'universal' implants for acute and chronic simultaneous brain-spinal cord imaging and optical stimulation. We further describe successful imaging of bioluminescent signals in both foci, and a new miniscope, the 'BLmini,' which has reduced weight, cost and form-factor relative to standard wearable miniscopes.

Conclusions

The combination of 3D printed implants, advanced imaging tools, and bioluminescence imaging techniques offers a new coalition of methods for understanding spinal cord-brain interactions. This work has the potential for use in future research into neuropathic pain and other sensory disorders and motor behavior.

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

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

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

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

Integrating brainstem and cortical functional architectures

The brainstem is a fundamental component of the central nervous system yet it is typically excluded from in vivo human brain mapping efforts, precluding a complete understanding of how the brainstem influences cortical function. Here we use high-resolution 7 Tesla fMRI to derive a functional connectome encompassing cortex as well as 58 brainstem nuclei spanning the midbrain, pons and medulla. We identify a compact set of integrative hubs in the brainstem with widespread connectivity with cerebral cortex. Patterns of connectivity between brainstem and cerebral cortex manifest as multiple emergent phenomena including neurophysiological oscillatory rhythms, patterns of cognitive functional specialization, and the unimodal-transmodal functional hierarchy. This persistent alignment between cortical functional topographies and brainstem nuclei is shaped by the spatial arrangement of multiple neurotransmitter receptors and transporters. We replicate all findings using 3 Tesla data from the same participants. Collectively, we find that multiple organizational features of cortical activity can be traced back to the brainstem.

Integrating brainstem and cortical functional architectures

The brainstem is a fundamental component of the central nervous system yet it is typically excluded from in vivo human brain mapping efforts, precluding a complete understanding of how the brainstem influences cortical function. Here we use high-resolution 7 Tesla fMRI to derive a functional connectome encompassing cortex as well as 58 brainstem nuclei spanning the midbrain, pons and medulla. We identify a compact set of integrative hubs in the brainstem with widespread connectivity with cerebral cortex. Patterns of connectivity between brainstem and cerebral cortex manifest as multiple emergent phenomena including neurophysiological oscillatory rhythms, patterns of cognitive functional specialization, and the unimodal-transmodal functional hierarchy. This persistent alignment between cortical functional topographies and brainstem nuclei is shaped by the spatial arrangement of multiple neurotransmitter receptors and transporters. We replicate all findings using 3 Tesla data from the same participants. Collectively, we find that multiple organizational features of cortical activity can be traced back to the brainstem.

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.

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.

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.

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

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

Simultaneous multislice imaging with slice-specific z-shim

PURPOSE

To implement slice-specific z-shim in simultaneous multislice (SMS) imaging in order to minimize signal losses in slice-accelerated T₂ *-weighted acquisitions, such as for spinal cord functional neuroimaging.

METHODS

The RF envelopes of the individual slice bands are temporally shifted on the plateau of the slice-selection gradient pulse before being combined to the multiband RF envelope. Thus, optimum z-shims can be realized for each slice of an SMS excitation, which is in contrast to conventional z-shimming. EPI with 2-fold SMS acceleration was performed on a 3T whole-body MR system in phantoms and the cervical spinal cord of healthy volunteers (i) without z-shim, (ii) with conventional z-shim using the average value of the slices of the SMS excitation, and (iii) with optimal, slice-specific z-shims for each slice using envelope shifts.

RESULTS

Phantom experiments demonstrate the equivalence of the envelope shift and conventional z-shimming for non-SMS excitations. With SMS, the best image quality is obtained with "mixed" z-shim, where only the z-shim differences of the slices of an SMS excitation are implemented by an envelope shift while their mean z-shim is applied conventionally with a gradient pulse after the echoes acquired for N/2 ghost correction. In phantoms and in vivo, this setup outperforms the approaches without slice-specific z-shim with respect to signal amplitude and temporal SNR at the expense of slight TE differences (<1 ms) between the slices.

CONCLUSION

With RF envelope shifts, slice-specific z-shims can be combined with SMS imaging, which could improve slice-accelerated functional neuroimaging in the spinal cord.

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

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

Altered Spinal Cord Functional Connectivity Associated with Parkinson's Disease Progression

BACKGROUND

Parkinson's disease (PD) has traditionally been viewed as an α-synucleinopathy brain pathology. Yet evidence based on postmortem human and animal experimental models indicates that the spinal cord may also be affected.

OBJECTIVE

Functional magnetic resonance imaging (fMRI) seems to be a promising candidate to better characterize spinal cord functional organization in PD patients.

METHODS

Resting-state spinal fMRI was performed in 70 PD patients and 24 age-matched healthy controls, the patients being divided into three groups based on their motor symptom severity: PD_low (n = 24), PD_med (n = 22), and PD_adv (n = 24) groups. A combination of independent component analysis (ICA) and a seed-based approach was applied.

RESULTS

When pooling all participants, the ICA revealed distinct ventral and dorsal components distributed along the rostro-caudal axis. This organization was highly reproducible within subgroups of patients and controls. PD severity, assessed by Unified Parkinson's Disease Rating Scale (UPDRS) scores, was associated with a decrease in spinal functional connectivity (FC). Notably, we observed a reduced intersegmental correlation in PD as compared to controls, the latter being negatively associated with patients' upper-limb UPDRS scores (P = 0.0085). This negative association between FC and upper-limb UPDRS scores was significant between adjacent C4-C5 (P = 0.015) and C5-C6 (P = 0.20) cervical segments, levels associated with upper-limb functions.

CONCLUSIONS

The present study provides the first evidence of spinal cord FC changes in PD and opens new avenues for the effective diagnosis and therapeutic strategies in PD. This underscores how spinal cord fMRI can serve as a powerful tool to characterize, in vivo, spinal circuits for a variety of neurological diseases. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Remodeling of the brain correlates with gait instability in cervical spondylotic myelopathy

Introduction

Cervical spondylotic myelopathy (CSM) is a common form of non-traumatic spinal cord injury (SCI) and usually leads to remodeling of the brain and spinal cord. In CSM with gait instability, the remodeling of the brain and cervical spinal cord is unclear. We attempted to explore the remodeling of these patients' brains and spinal cords, as well as the relationship between the remodeling of the brain and spinal cord and gait instability.

Methods

According to the CSM patients' gait, we divided patients into two groups: normal gait patients (nPT) and abnormal gait patients (aPT). Voxel-wise z-score transformation amplitude of low-frequency fluctuations (zALFF) and resting-state functional connectivity (rs-FC) were performed for estimating brain changes. Cross-sectional area (CSA) and fractional anisotropy (FA) of the spinal cord were computed by Spinal cord toolbox. Correlations of these measures and the modified Japanese Orthopedic Association (mJOA) score were analyzed.

Results

We found that the zALFF of caudate nucleus in aPT was higher than that in healthy controls (HC) and lower than that in nPT. The zALFF of the right postcentral gyrus and paracentral lobule in HC was higher than those of aPT and nPT. Compared with the nPT, the aPT showed increased functional connectivity between the caudate nucleus and left angular gyrus, bilateral precuneus and bilateral posterior cingulate cortex (PCC), which constitute a vital section of the default mode network (DMN). No significantly different FA values or CSA of spinal tracts at the C2 level were observed between the HC, nPT and aPT groups. In CSM, the right paracentral lobule's zALFF was negatively correlated with the FA value of fasciculus gracilis (FCG), and the right caudate zALFF was positively correlated with the FA value of the fasciculus cuneatus (FCC). The results showed that the functional connectivity between the right caudate nucleus and DMN was negatively correlated with the CSA of the lateral corticospinal tract (CST).

Discussion

The activation of the caudate nucleus and the strengthening functional connectivity between the caudate nucleus and DMN were associated with gait instability in CSM patients. Correlations between spinal cord and brain function might be related to the clinical symptoms in CSM.

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.

Automated slice-specific z-shimming for functional magnetic resonance imaging of the human spinal cord

Functional magnetic resonance imaging (fMRI) of the human spinal cord faces many challenges, such as signal loss due to local magnetic field inhomogeneities. This issue can be addressed with slice-specific z-shimming, which compensates for the dephasing effect of the inhomogeneities using a slice-specific gradient pulse. Here, we aim to address outstanding issues regarding this technique by evaluating its effects on several aspects that are directly relevant for spinal fMRI and by developing two automated procedures in order to improve upon the time-consuming and subjective nature of manual selection of z-shims: one procedure finds the z-shim that maximizes signal intensity in each slice of an EPI reference-scan and the other finds the through-slice field inhomogeneity for each EPI-slice in field map data and calculates the required compensation gradient moment. We demonstrate that the beneficial effects of z-shimming are apparent across different echo times, hold true for both the dorsal and ventral horn, and are also apparent in the temporal signal-to-noise ratio (tSNR) of EPI time-series data. Both of our automated approaches were faster than the manual approach, lead to significant improvements in gray matter tSNR compared to no z-shimming and resulted in beneficial effects that were stable across time. While the field-map-based approach performed slightly worse than the manual approach, the EPI-based approach performed as well as the manual one and was furthermore validated on an external corticospinal data-set (N > 100). Together, automated z-shimming may improve the data quality of future spinal fMRI studies and lead to increased reproducibility in longitudinal studies.

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

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

The Motor Engram of Functional Connectivity Generated by Acute Whole-Body Dynamic Balance Training

PURPOSE

Whole-body dynamic balance is necessary for both athletic activities and activities of daily living. This study aimed to investigate the effect of acute dynamic balance training on neural networks.

METHODS

We evaluated resting-state functional connectivity (rs-FC), white matter fiber density, fiber-bundle cross-section, and gray matter volume in 28 healthy young adults (14 women) before and after 30 min of slackline training using a randomized, counterbalanced crossover design.

RESULTS

The rs-FC between the left lateral prefrontal cortex (PFC) and the foot area of the primary sensorimotor (SM1) cortex increased significantly after slackline training compared with that after a control condition involving ergometer-based aerobic exercise. In addition, changes in rs-FC between the left lateral PFC and the primary sensorimotor were correlated with performance changes after training (i.e., offline process) rather than online learning. We also observed a main effect of time between the hippocampus and the cingulate cortex, including the anterior areas, and between the bilateral lateral PFC. Although we observed no structural changes, fiber density in the commissural fiber pathway before the first balance assessment was correlated with initial balance capability.

CONCLUSIONS

Our findings demonstrate that acute whole-body dynamic balance training alters specific rs-FC, and that this change is associated with performance changes after training. In addition, rs-FC changes in cognitive regions were modulated by both acute dynamic balance training and aerobic exercise. These findings have the potential to influence various fields (e.g., sports neuroscience, neurorehabilitation) and may aid in the development of methods that can improve motor and cognitive performance.

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

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

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.

Connectomic mapping of brain-spinal cord neural networks: future directions in assessing spinal cord injury at rest

Following spinal cord injury (SCI), the central nervous system undergoes significant reconstruction. The dynamic change in the interaction of the brain-spinal cord axis as well as in structure-function relations plays a vital role in the determination of neurological functions, which might have important clinical implications for the treatment and its efficacy evaluation of patients with SCI. Brain connectomes based on neuroimaging data is a relatively new field of research that maps the brain's large-scale structural and functional networks at rest. Importantly, increasing evidence shows that such resting-state signals can also be seen in the spinal cord. In the present review, we focus on the reconstruction of multi-level neural circuits after SCI. We also describe how the connectome concept could further our understanding of neuroplasticity after SCI. We propose that mapping the cortical-subcortical-spinal cord networks can provide novel insights into the pathologies of SCI.

Generic acquisition protocol for quantitative MRI of the spinal cord

Quantitative spinal cord (SC) magnetic resonance imaging (MRI) presents many challenges, including a lack of standardized imaging protocols. Here we present a prospectively harmonized quantitative MRI protocol, which we refer to as the spine generic protocol, for users of 3T MRI systems from the three main manufacturers: GE, Philips and Siemens. The protocol provides guidance for assessing SC macrostructural and microstructural integrity: T1-weighted and T2-weighted imaging for SC cross-sectional area computation, multi-echo gradient echo for gray matter cross-sectional area, and magnetization transfer and diffusion weighted imaging for assessing white matter microstructure. In a companion paper from the same authors, the spine generic protocol was used to acquire data across 42 centers in 260 healthy subjects. The key details of the spine generic protocol are also available in an open-access document that can be found at https://github.com/spine-generic/protocols . The protocol will serve as a starting point for researchers and clinicians implementing new SC imaging initiatives so that, in the future, inclusion of the SC in neuroimaging protocols will be more common. The protocol could be implemented by any trained MR technician or by a researcher/clinician familiar with MRI acquisition.

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.

Observation of others' painful heat stimulation involves responses in the spinal cord

Observing others' aversive experiences is central to know what is dangerous for ourselves. Hence, observation often elicits behavioral and physiological responses comparable to first-hand aversive experiences and engages overlapping brain activation. While brain activation to first-hand aversive experiences relies on connections to the spinal cord, it is unresolved whether merely observing aversive stimulation also involves responses in the spinal cord. Here, we show that observation of others receiving painful heat stimulation involves neural responses in the spinal cord, located in the same cervical segment as first-hand heat pain. However, while first-hand painful experiences are coded within dorsolateral regions of the spinal cord, observation of others' painful heat stimulation involves medial regions. Dorsolateral areas that process first-hand pain exhibit negative responses when observing pain in others. Our results suggest a distinct processing between self and others' pain in the spinal cord when integrating social information.

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.