Body size interacts with the structure of the central nervous system: A multi-center in vivo neuroimaging study

Clinical research emphasizes the implementation of rigorous and reproducible study designs that rely on between-group matching or controlling for sources of biological variation such as subject's sex and age. However, corrections for body size (i.e. height and weight) are mostly lacking in clinical neuroimaging designs. This study investigates the importance of body size parameters in their relationship with spinal cord (SC) and brain magnetic resonance imaging (MRI) metrics. Data were derived from a cosmopolitan population of 267 healthy human adults (age 30.1±6.6 years old, 125 females). We show that body height correlated strongly or moderately with brain gray matter (GM) volume, cortical GM volume, total cerebellar volume, brainstem volume, and cross-sectional area (CSA) of cervical SC white matter (CSA-WM; 0.44≤r≤0.62). In comparison, age correlated weakly with cortical GM volume, precentral GM volume, and cortical thickness (-0.21≥r≥-0.27). Body weight correlated weakly with magnetization transfer ratio in the SC WM, dorsal columns, and lateral corticospinal tracts (-0.20≥r≥-0.23). Body weight further correlated weakly with the mean diffusivity derived from diffusion tensor imaging (DTI) in SC WM (r=-0.20) and dorsal columns (-0.21), but only in males. CSA-WM correlated strongly or moderately with brain volumes (0.39≤r≤0.64), and weakly with precentral gyrus thickness and DTI-based fractional anisotropy in SC dorsal columns and SC lateral corticospinal tracts (-0.22≥r≥-0.25). Linear mixture of sex and age explained 26±10% of data variance in brain volumetry and SC CSA. The amount of explained variance increased at 33±11% when body height was added into the mixture model. Age itself explained only 2±2% of such variance. In conclusion, body size is a significant biological variable. Along with sex and age, body size should therefore be included as a mandatory variable in the design of clinical neuroimaging studies examining SC and brain structure.

Investigating the associations between lumbar paraspinal muscle health and age, BMI, sex, physical activity, and back pain using an automated computer-vision model: A UK Biobank study

BACKGROUND CONTEXT

The role of lumbar paraspinal muscle health in back pain (BP) is not straightforward. Challenges in this field have included the lack of tools and large, heterogenous datasets to interrogate the association between muscle health and BP. Computer-vision models have been transformative in this space, enabling the automated quantification of muscle health and the processing of large datasets.

PURPOSE

To investigate the associations between lumbar paraspinal muscle health and age, sex, BMI, physical activity, and BP in a large, heterogenous dataset using an automated computer-vision model.

DESIGN

Cross-sectional study.

PATIENT SAMPLE

Participants from the UK Biobank with abdominal Dixon fat-water MRI (N=9,564) were included (41.8% women, mean [SD] age: 63.5 [7.6] years, BMI: 26.4 [4.1] kg/m2) of whom 6,953 reported no pain, 930 acute BP, and 1,681 chronic BP.

OUTCOME MEASURES

Intramuscular fat (IMF) and average cross-sectional area (aCSA) were automatically derived using a computer-vision model for the left and right lumbar multifidus (LM), erector spinae (ES), and psoas major (PM) from the L1 to L5 vertebral levels.

METHODS

Two-tailed partial Pearson correlations were generated for each muscle to assess the relationships between the muscle measures (IMF and aCSA) and age (controlling for BMI, sex, and physical activity), BMI (controlling for age, sex, and physical activity), and physical activity (controlling for age, sex, and BMI). One-way ANCOVA was used to identify sex differences in IMF and aCSA for each muscle while controlling for age, BMI, and physical activity. Similarly, one-way ANCOVA was used to identify between-group differences (no pain, acute BP, and chronic BP) for each muscle and along the superior-inferior expanse of the lumbar spine while controlling for age, BMI, sex, and physical activity (α=0.05).

RESULTS

Females had higher IMF (LM mean difference [MD]=11.1%, ES MD=10.2%, PM MD=0.3%, p<.001) and lower aCSA (LM MD=47.6 mm2, ES MD=350.0 mm2, PM MD=321.5 mm2, p<.001) for all muscles. Higher age was associated with higher IMF and lower aCSA for all muscles (r≥0.232, p<.001) except for LM and aCSA (r≤0.013, p≥.267). Higher BMI was associated with higher IMF and aCSA for all muscles (r≥0.174, p<.001). Higher physical activity was associated with lower IMF and higher aCSA for all muscles (r≥0.036, p≤.002) except for LM and aCSA (r≤0.010, p≥.405). People with chronic BP had higher IMF and lower aCSA than people with no pain (IMF MD≤1.6%, aCSA MD≤27.4 mm2, p<.001) and higher IMF compared to acute BP (IMF MD≤1.1%, p≤.044). The differences between people with BP and people with no pain were not spatially localized to the inferior lumbar levels but broadly distributed across the lumbar spine.

CONCLUSIONS

Paraspinal muscle health is associated with age, BMI, sex, and physical activity with the exception of the association between LM aCSA and age and physical activity. People with BP (chronic>acute) have higher IMF and lower aCSA than people reporting no pain. The differences were not localized but broadly distributed across the lumbar spine. When interpreting measures of paraspinal muscle health in the research or clinical setting, the associations with age, BMI, sex, and physical activity should be considered.

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.

Severe central nervous system demyelination in Sanfilippo disease

Introduction

Chronic progressive neuroinflammation is a hallmark of neurological lysosomal storage diseases, including mucopolysaccharidosis III (MPS III or Sanfilippo disease). Since neuroinflammation is linked to white matter tract pathology, we analyzed axonal myelination and white matter density in the mouse model of MPS IIIC HgsnatP304L and post-mortem brain samples of MPS III patients.

Methods

Brain and spinal cord tissues of human MPS III patients, 6-month-old HgsnatP304L mice and age- and sex-matching wild type mice were analyzed by immunofluorescence to assess levels of myelin-associated proteins, primary and secondary storage materials, and levels of microgliosis. Corpus callosum (CC) region was studied by transmission electron microscopy to analyze axon myelination and morphology of oligodendrocytes and microglia. Mouse brains were analyzed ex vivo by high-filed MRI using Diffusion Basis Spectrum Imaging in Python-Diffusion tensor imaging algorithms.

Results

Analyses of CC and spinal cord tissues by immunohistochemistry revealed substantially reduced levels of myelin-associated proteins including Myelin Basic Protein, Myelin Associated Glycoprotein, and Myelin Oligodendrocyte Glycoprotein. Furthermore, ultrastructural analyses revealed disruption of myelin sheath organization and reduced myelin thickness in the brains of MPS IIIC mice and human MPS IIIC patients compared to healthy controls. Oligodendrocytes (OLs) in the CC of MPS IIIC mice were scarce, while examination of the remaining cells revealed numerous enlarged lysosomes containing heparan sulfate, GM3 ganglioside or "zebra bodies" consistent with accumulation of lipids and myelin fragments. In addition, OLs contained swollen mitochondria with largely dissolved cristae, resembling those previously identified in the dysfunctional neurons of MPS IIIC mice. Ex vivo Diffusion Basis Spectrum Imaging revealed compelling signs of demyelination (26% increase in radial diffusivity) and tissue loss (76% increase in hindered diffusivity) in CC of MPS IIIC mice.

Discussion

Our findings demonstrate an important role for white matter injury in the pathophysiology of MPS III. This study also defines specific parameters and brain regions for MRI analysis and suggests that it may become a crucial non-invasive method to evaluate disease progression and therapeutic response.

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.

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.

Centering inclusivity in the design of online conferences-An OHBM-Open Science perspective

As the global health crisis unfolded, many academic conferences moved online in 2020. This move has been hailed as a positive step towards inclusivity in its attenuation of economic, physical, and legal barriers and effectively enabled many individuals from groups that have traditionally been underrepresented to join and participate. A number of studies have outlined how moving online made it possible to gather a more global community and has increased opportunities for individuals with various constraints, e.g., caregiving responsibilities. Yet, the mere existence of online conferences is no guarantee that everyone can attend and participate meaningfully. In fact, many elements of an online conference are still significant barriers to truly diverse participation: the tools used can be inaccessible for some individuals; the scheduling choices can favour some geographical locations; the set-up of the conference can provide more visibility to well-established researchers and reduce opportunities for early-career researchers. While acknowledging the benefits of an online setting, especially for individuals who have traditionally been underrepresented or excluded, we recognize that fostering social justice requires inclusivity to actively be centered in every aspect of online conference design. Here, we draw from the literature and from our own experiences to identify practices that purposefully encourage a diverse community to attend, participate in, and lead online conferences. Reflecting on how to design more inclusive online events is especially important as multiple scientific organizations have announced that they will continue offering an online version of their event when in-person conferences can resume.

Open-access quantitative MRI data of the spinal cord and reproducibility across participants, sites and manufacturers

In a companion paper by Cohen-Adad et al. we introduce the spine generic quantitative MRI protocol that provides valuable metrics for assessing spinal cord macrostructural and microstructural integrity. This protocol was used to acquire a single subject dataset across 19 centers and a multi-subject dataset across 42 centers (for a total of 260 participants), spanning the three main MRI manufacturers: GE, Philips and Siemens. Both datasets are publicly available via git-annex. Data were analysed using the Spinal Cord Toolbox to produce normative values as well as inter/intra-site and inter/intra-manufacturer statistics. Reproducibility for the spine generic protocol was high across sites and manufacturers, with an average inter-site coefficient of variation of less than 5% for all the metrics. Full documentation and results can be found at https://spine-generic.rtfd.io/ . The datasets and analysis pipeline will help pave the way towards accessible and reproducible quantitative MRI in the spinal cord.

7 T imaging reveals a gradient in spinal cord lesion distribution in multiple sclerosis

We used 7 T MRI to: (i) characterize the grey and white matter pathology in the cervical spinal cord of patients with early relapsing-remitting and secondary progressive multiple sclerosis; (ii) assess the spinal cord lesion spatial distribution and the hypothesis of an outside-in pathological process possibly driven by CSF-mediated immune cytotoxic factors; and (iii) evaluate the association of spinal cord pathology with brain burden and its contribution to neurological disability. We prospectively recruited 20 relapsing-remitting, 15 secondary progressive multiple sclerosis participants and 11 age-matched healthy control subjects to undergo 7 T imaging of the cervical spinal cord and brain as well as conventional 3 T brain acquisition. Cervical spinal cord imaging at 7 T was used to segment grey and white matter, including lesions therein. Brain imaging at 7 T was used to segment cortical and white matter lesions and 3 T imaging for cortical thickness estimation. Cervical spinal cord lesions were mapped voxel-wise as a function of distance from the inner central canal CSF pool to the outer subpial surface. Similarly, brain white matter lesions were mapped voxel-wise as a function of distance from the ventricular system. Subjects with relapsing-remitting multiple sclerosis showed a greater predominance of spinal cord lesions nearer the outer subpial surface compared to secondary progressive cases. Inversely, secondary progressive participants presented with more centrally located lesions. Within the brain, there was a strong gradient of lesion formation nearest the ventricular system that was most evident in participants with secondary progressive multiple sclerosis. Lesion fractions within the spinal cord grey and white matter were related to the lesion fraction in cerebral white matter. Cortical thinning was the primary determinant of the Expanded Disability Status Scale, white matter lesion fractions in the spinal cord and brain of the 9-Hole Peg Test and cortical thickness and spinal cord grey matter cross-sectional area of the Timed 25-Foot Walk. Spinal cord lesions were localized nearest the subpial surfaces for those with relapsing-remitting and the central canal CSF surface in progressive disease, possibly implying CSF-mediated pathogenic mechanisms in lesion development that may differ between multiple sclerosis subtypes. These findings show that spinal cord lesions involve both grey and white matter from the early multiple sclerosis stages and occur mostly independent from brain pathology. Despite the prevalence of cervical spinal cord lesions and atrophy, brain pathology seems more strongly related to physical disability as measured by the Expanded Disability Status Scale.

Spatial distribution of multiple sclerosis lesions in the cervical spinal cord

Spinal cord lesions detected on MRI hold important diagnostic and prognostic value for multiple sclerosis. Previous attempts to correlate lesion burden with clinical status have had limited success, however, suggesting that lesion location may be a contributor. Our aim was to explore the spatial distribution of multiple sclerosis lesions in the cervical spinal cord, with respect to clinical status. We included 642 suspected or confirmed multiple sclerosis patients (31 clinically isolated syndrome, and 416 relapsing-remitting, 84 secondary progressive, and 73 primary progressive multiple sclerosis) from 13 clinical sites. Cervical spine lesions were manually delineated on T2- and T2*-weighted axial and sagittal MRI scans acquired at 3 or 7 T. With an automatic publicly-available analysis pipeline we produced voxelwise lesion frequency maps to identify predilection sites in various patient groups characterized by clinical subtype, Expanded Disability Status Scale score and disease duration. We also measured absolute and normalized lesion volumes in several regions of interest using an atlas-based approach, and evaluated differences within and between groups. The lateral funiculi were more frequently affected by lesions in progressive subtypes than in relapsing in voxelwise analysis (P < 0.001), which was further confirmed by absolute and normalized lesion volumes (P < 0.01). The central cord area was more often affected by lesions in primary progressive than relapse-remitting patients (P < 0.001). Between white and grey matter, the absolute lesion volume in the white matter was greater than in the grey matter in all phenotypes (P < 0.001); however when normalizing by each region, normalized lesion volumes were comparable between white and grey matter in primary progressive patients. Lesions appearing in the lateral funiculi and central cord area were significantly correlated with Expanded Disability Status Scale score (P < 0.001). High lesion frequencies were observed in patients with a more aggressive disease course, rather than long disease duration. Lesions located in the lateral funiculi and central cord area of the cervical spine may influence clinical status in multiple sclerosis. This work shows the added value of cervical spine lesions, and provides an avenue for evaluating the distribution of spinal cord lesions in various patient groups.

Automatic segmentation of the spinal cord and intramedullary multiple sclerosis lesions with convolutional neural networks

The spinal cord is frequently affected by atrophy and/or lesions in multiple sclerosis (MS) patients. Segmentation of the spinal cord and lesions from MRI data provides measures of damage, which are key criteria for the diagnosis, prognosis, and longitudinal monitoring in MS. Automating this operation eliminates inter-rater variability and increases the efficiency of large-throughput analysis pipelines. Robust and reliable segmentation across multi-site spinal cord data is challenging because of the large variability related to acquisition parameters and image artifacts. In particular, a precise delineation of lesions is hindered by a broad heterogeneity of lesion contrast, size, location, and shape. The goal of this study was to develop a fully-automatic framework - robust to variability in both image parameters and clinical condition - for segmentation of the spinal cord and intramedullary MS lesions from conventional MRI data of MS and non-MS cases. Scans of 1042 subjects (459 healthy controls, 471 MS patients, and 112 with other spinal pathologies) were included in this multi-site study (n = 30). Data spanned three contrasts (T1-, T2-, and T2∗-weighted) for a total of 1943 vol and featured large heterogeneity in terms of resolution, orientation, coverage, and clinical conditions. The proposed cord and lesion automatic segmentation approach is based on a sequence of two Convolutional Neural Networks (CNNs). To deal with the very small proportion of spinal cord and/or lesion voxels compared to the rest of the volume, a first CNN with 2D dilated convolutions detects the spinal cord centerline, followed by a second CNN with 3D convolutions that segments the spinal cord and/or lesions. CNNs were trained independently with the Dice loss. When compared against manual segmentation, our CNN-based approach showed a median Dice of 95% vs. 88% for PropSeg (p ≤ 0.05), a state-of-the-art spinal cord segmentation method. Regarding lesion segmentation on MS data, our framework provided a Dice of 60%, a relative volume difference of -15%, and a lesion-wise detection sensitivity and precision of 83% and 77%, respectively. In this study, we introduce a robust method to segment the spinal cord and intramedullary MS lesions on a variety of MRI contrasts. The proposed framework is open-source and readily available in the Spinal Cord Toolbox.

Monitoring for myelopathic progression with multiparametric quantitative MRI

Background

Patients with mild degenerative cervical myelopathy (DCM) are often managed non-operatively, and surgery is recommended if neurological progression occurs. However, detection of progression is often subjective. Quantitative MRI (qMRI) directly measures spinal cord (SC) tissue changes, detecting axonal injury, demyelination, and atrophy. This longitudinal study compared multiparametric qMRI with clinical measures of progression in non-operative DCM patients.

Methods

26 DCM patients were followed. Clinical data included modified Japanese Orthopedic Association (mJOA) and additional assessments. 3T qMRI data included cross sectional area, diffusion fractional anisotropy, magnetization transfer ratio, and T2*-weighted white/grey matter signal ratio, extracted from the compressed SC and above/below. Progression was defined as 1) patients’ subjective impression, 2) 2-point mJOA decrease, 3) ≥3 clinical measures worsening ≥5%, 4) increased compression on MRI, or 5) ≥1 of 10 qMRI measures or composite score worsening (p < 0.004, corrected).

Results

Follow-up (13.5 ± 4.9 months) included mJOA in all 26 patients, MRI in 25, and clinical/qMRI in 22. 42.3% reported subjective worsening, compared with mJOA (11.5%), MRI (20%), comprehensive assessments (54.6%), and qMRI (68.2%). Relative to subjective worsening, qMRI showed 100% sensitivity and 53.3% specificity compared with comprehensive assessments (75%, 60%), mJOA (27.3%, 100%), and MRI (18.2%, 81.3%). A decision-making algorithm incorporating qMRI identified progression and recommended surgery for 11 subjects (42.3%).

Conclusions

Quantitative MRI shows high sensitivity to detect myelopathic progression. Our results suggest that neuroplasticity and behavioural adaptation may mask progressive SC tissue injury. qMRI appears to be a useful method to confirm subtle myelopathic progression in individual patients, representing an advance toward clinical translation of qMRI.

Can microstructural MRI detect subclinical tissue injury in subjects with asymptomatic cervical spinal cord compression? A prospective cohort study

OBJECTIVES

Degenerative cervical myelopathy (DCM) involves extrinsic spinal cord compression causing tissue injury and neurological dysfunction. Asymptomatic spinal cord compression (ASCC) is more common, but its significance is poorly defined. This study investigates if: (1) ASCC can be automatically diagnosed using spinal cord shape analysis; (2) multiparametric quantitative MRI can detect similar spinal cord tissue injury as previously observed in DCM.

DESIGN

Prospective observational longitudinal cohort study.

SETTING

Single centre, tertiary care and research institution.

PARTICIPANTS

40 neurologically intact subjects (19 female, 21 male) divided into groups with and without ASCC.

INTERVENTIONS

None.

OUTCOME MEASURES

Clinical assessments: modified Japanese Orthopaedic Association score and physical examination. 3T MRI assessments: automated morphometric analysis compared with consensus ratings of spinal cord compression, and measures of tissue injury: cross-sectional area, diffusion fractional anisotropy, magnetisation transfer ratio and T2*-weighted imaging white to grey matter signal intensity ratio (T2*WI WM/GM) extracted from rostral (C1-3), caudal (C6-7) and maximally compressed levels.

RESULTS

ASCC was present in 20/40 subjects. Diagnosis with automated shape analysis showed area under the curve >97%. Five MRI metrics showed differences suggestive of tissue injury in ASCC compared with uncompressed subjects (p<0.05), while a composite of all 10 measures (average of z scores) showed highly significant differences (p=0.002). At follow-up (median 21 months), two ASCC subjects developed DCM.

CONCLUSIONS

ASCC appears to be common and can be accurately and objectively diagnosed with automated morphometric analysis. Quantitative MRI appears to detect subclinical tissue injury in ASCC prior to the onset of neurological symptoms and signs. These findings require further validation, but offer the intriguing possibility of presymptomatic diagnosis and treatment of DCM and other spinal pathologies.

Automatic spinal cord localization, robust to MRI contrasts using global curve optimization

During the last two decades, MRI has been increasingly used for providing valuable quantitative information about spinal cord morphometry, such as quantification of the spinal cord atrophy in various diseases. However, despite the significant improvement of MR sequences adapted to the spinal cord, automatic image processing tools for spinal cord MRI data are not yet as developed as for the brain. There is nonetheless great interest in fully automatic and fast processing methods to be able to propose quantitative analysis pipelines on large datasets without user bias. The first step of most of these analysis pipelines is to detect the spinal cord, which is challenging to achieve automatically across the broad range of MRI contrasts, field of view, resolutions and pathologies. In this paper, a fully automated, robust and fast method for detecting the spinal cord centerline on MRI volumes is introduced. The algorithm uses a global optimization scheme that attempts to strike a balance between a probabilistic localization map of the spinal cord center point and the overall spatial consistency of the spinal cord centerline (i.e. the rostro-caudal continuity of the spinal cord). Additionally, a new post-processing feature, which aims to automatically split brain and spine regions is introduced, to be able to detect a consistent spinal cord centerline, independently from the field of view. We present data on the validation of the proposed algorithm, known as "OptiC", from a large dataset involving 20 centers, 4 contrasts (T₂-weighted n = 287, T₁-weighted n = 120, T₂^∗-weighted n = 307, diffusion-weighted n = 90), 501 subjects including 173 patients with a variety of neurologic diseases. Validation involved the gold-standard centerline coverage, the mean square error between the true and predicted centerlines and the ability to accurately separate brain and spine regions. Overall, OptiC was able to cover 98.77% of the gold-standard centerline, with a mean square error of 1.02 mm. OptiC achieved superior results compared to a state-of-the-art spinal cord localization technique based on the Hough transform, especially on pathological cases with an averaged mean square error of 1.08 mm vs. 13.16 mm (Wilcoxon signed-rank test p-value < .01). Images containing brain regions were identified with a 99% precision, on which brain and spine regions were separated with a distance error of 9.37 mm compared to ground-truth. Validation results on a challenging dataset suggest that OptiC could reliably be used for subsequent quantitative analyses tasks, opening the door to more robust analysis on pathological cases.

Spinal cord grey matter segmentation challenge

An important image processing step in spinal cord magnetic resonance imaging is the ability to reliably and accurately segment grey and white matter for tissue specific analysis. There are several semi- or fully-automated segmentation methods for cervical cord cross-sectional area measurement with an excellent performance close or equal to the manual segmentation. However, grey matter segmentation is still challenging due to small cross-sectional size and shape, and active research is being conducted by several groups around the world in this field. Therefore a grey matter spinal cord segmentation challenge was organised to test different capabilities of various methods using the same multi-centre and multi-vendor dataset acquired with distinct 3D gradient-echo sequences. This challenge aimed to characterize the state-of-the-art in the field as well as identifying new opportunities for future improvements. Six different spinal cord grey matter segmentation methods developed independently by various research groups across the world and their performance were compared to manual segmentation outcomes, the present gold-standard. All algorithms provided good overall results for detecting the grey matter butterfly, albeit with variable performance in certain quality-of-segmentation metrics. The data have been made publicly available and the challenge web site remains open to new submissions. No modifications were introduced to any of the presented methods as a result of this challenge for the purposes of this publication.

Topologically preserving straightening of spinal cord MRI

PURPOSE

To propose a robust and accurate method for straightening magnetic resonance (MR) images of the spinal cord, based on spinal cord segmentation, that preserves spinal cord topology and that works for any MRI contrast, in a context of spinal cord template-based analysis.

MATERIALS AND METHODS

The spinal cord curvature was computed using an iterative Non-Uniform Rational B-Spline (NURBS) approximation. Forward and inverse deformation fields for straightening were computed by solving analytically the straightening equations for each image voxel. Computational speed-up was accomplished by solving all voxel equation systems as one single system. Straightening accuracy (mean and maximum distance from straight line), computational time, and robustness to spinal cord length was evaluated using the proposed and the standard straightening method (label-based spline deformation) on 3T T₂ - and T₁ -weighted images from 57 healthy subjects and 33 patients with spinal cord compression due to degenerative cervical myelopathy (DCM).

RESULTS

The proposed algorithm was more accurate, more robust, and faster than the standard method (mean distance = 0.80 vs. 0.83 mm, maximum distance = 1.49 vs. 1.78 mm, time = 71 vs. 174 sec for the healthy population and mean distance = 0.65 vs. 0.68 mm, maximum distance = 1.28 vs. 1.55 mm, time = 32 vs. 60 sec for the DCM population).

CONCLUSION

A novel image straightening method that enables template-based analysis of quantitative spinal cord MRI data is introduced. This algorithm works for any MRI contrast and was validated on healthy and patient populations. The presented method is implemented in the Spinal Cord Toolbox, an open-source software for processing spinal cord MRI data.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2017;46:1209-1219.

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

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

Fully-integrated framework for the segmentation and registration of the spinal cord white and gray matter

The spinal cord white and gray matter can be affected by various pathologies such as multiple sclerosis, amyotrophic lateral sclerosis or trauma. Being able to precisely segment the white and gray matter could help with MR image analysis and hence be useful in further understanding these pathologies, and helping with diagnosis/prognosis and drug development. Up to date, white/gray matter segmentation has mostly been done manually, which is time consuming, induces a bias related to the rater and prevents large-scale multi-center studies. Recently, few methods have been proposed to automatically segment the spinal cord white and gray matter. However, no single method exists that combines the following criteria: (i) fully automatic, (ii) works on various MRI contrasts, (iii) robust towards pathology and (iv) freely available and open source. In this study we propose a multi-atlas based method for the segmentation of the spinal cord white and gray matter that addresses the previous limitations. Moreover, to study the spinal cord morphology, atlas-based approaches are increasingly used. These approaches rely on the registration of a spinal cord template to an MR image, however the registration usually doesn't take into account the spinal cord internal structure and thus lacks accuracy. In this study, we propose a new template registration framework that integrates the white and gray matter segmentation to account for the specific gray matter shape of each individual subject. Validation of segmentation was performed in 24 healthy subjects using T₂^*-weighted images, in 8 healthy subjects using diffusion weighted images (exhibiting inverted white-to-gray matter contrast compared to T₂*-weighted), and in 5 patients with spinal cord injury. The template registration was validated in 24 subjects using T₂*-weighted data. Results of automatic segmentation on T₂*-weighted images was in close correspondence with the manual segmentation (Dice coefficient in the white/gray matter of 0.91/0.71 respectively). Similarly, good results were obtained in data with inverted contrast (diffusion-weighted image) and in patients. When compared to the classical template registration framework, the proposed framework that accounts for gray matter shape significantly improved the quality of the registration (comparing Dice coefficient in gray matter: p=9.5×10^-6). While further validation is needed to show the benefits of the new registration framework in large cohorts and in a variety of patients, this study provides a fully-integrated tool for quantitative assessment of white/gray matter morphometry and template-based analysis. All the proposed methods are implemented in the Spinal Cord Toolbox (SCT), an open-source software for processing spinal cord multi-parametric MRI data.

Brainhack: a collaborative workshop for the open neuroscience community

Brainhack events offer a novel workshop format with participant-generated content that caters to the rapidly growing open neuroscience community. Including components from hackathons and unconferences, as well as parallel educational sessions, Brainhack fosters novel collaborations around the interests of its attendees. Here we provide an overview of its structure, past events, and example projects. Additionally, we outline current innovations such as regional events and post-conference publications. Through introducing Brainhack to the wider neuroscience community, we hope to provide a unique conference format that promotes the features of collaborative, open science.

Segmentation of the human spinal cord

Segmenting the spinal cord contour is a necessary step for quantifying spinal cord atrophy in various diseases. Delineating gray matter (GM) and white matter (WM) is also useful for quantifying GM atrophy or for extracting multiparametric MRI metrics into specific WM tracts. Spinal cord segmentation in clinical research is not as developed as brain segmentation, however with the substantial improvement of MR sequences adapted to spinal cord MR investigations, the field of spinal cord MR segmentation has advanced greatly within the last decade. Segmentation techniques with variable accuracy and degree of complexity have been developed and reported in the literature. In this paper, we review some of the existing methods for cord and WM/GM segmentation, including intensity-based, surface-based, and image-based methods. We also provide recommendations for validating spinal cord segmentation techniques, as it is important to understand the intrinsic characteristics of the methods and to evaluate their performance and limitations. Lastly, we illustrate some applications in the healthy and pathological spinal cord. One conclusion of this review is that robust and automatic segmentation is clinically relevant, as it would allow for longitudinal and group studies free from user bias as well as reproducible multicentric studies in large populations, thereby helping to further our understanding of the spinal cord pathophysiology and to develop new criteria for early detection of subclinical evolution for prognosis prediction and for patient management. Another conclusion is that at the present time, no single method adequately segments the cord and its substructure in all the cases encountered (abnormal intensities, loss of contrast, deformation of the cord, etc.). A combination of different approaches is thus advised for future developments, along with the introduction of probabilistic shape models. Maturation of standardized frameworks, multiplatform availability, inclusion in large suite and data sharing would also ultimately benefit to the community.

Fully automated segmentation of the cervical cord from T1-weighted MRI using PropSeg: Application to multiple sclerosis

Spinal cord (SC) atrophy, i.e. a reduction in the SC cross-sectional area (CSA) over time, can be measured by means of image segmentation using magnetic resonance imaging (MRI). However, segmentation methods have been limited by factors relating to reproducibility or sensitivity to change. The purpose of this study was to evaluate a fully automated SC segmentation method (PropSeg), and compare this to a semi-automated active surface (AS) method, in healthy controls (HC) and people with multiple sclerosis (MS). MRI data from 120 people were retrospectively analysed; 26 HC, 21 with clinically isolated syndrome, 26 relapsing remitting MS, 26 primary and 21 secondary progressive MS. MRI data from 40 people returning after one year were also analysed. CSA measurements were obtained within the cervical SC. Reproducibility of the measurements was assessed using the intraclass correlation coefficient (ICC). A comparison between mean CSA changes obtained with the two methods over time was performed using multivariate structural equation regression models. Associations between CSA measures and clinical scores were investigated using linear regression models. Compared to the AS method, the reproducibility of CSA measurements obtained with PropSeg was high, both in patients and in HC, with ICC > 0.98 in all cases. There was no significant difference between PropSeg and AS in terms of detecting change over time. Furthermore, PropSeg provided measures that correlated with physical disability, similar to the AS method. PropSeg is a time-efficient and reliable segmentation method, which requires no manual intervention, and may facilitate large multi-centre neuroprotective trials in progressive MS.

•

PropSeg is a fully automated segmentation method to measure cord atrophy in MS.

•

PropSeg offers comparable results to a widely used semi-automated method.

•

PropSeg is available for immediate clinical utility.

Automatic Segmentation of the Spinal Cord and Spinal Canal Coupled With Vertebral Labeling

Quantifying spinal cord (SC) atrophy in neurodegenerative and traumatic diseases brings important diagnosis and prognosis information for the clinician. We recently developed the PropSeg method, which allows for fast, accurate and automatic segmentation of the SC on different types of MRI contrast (e.g., T1-, T2- and T2(∗) -weighted sequences) and any field of view. However, comparing measurements from the SC between subjects is hindered by the lack of a generic coordinate system for the SC. In this paper, we present a new framework combining PropSeg and a vertebral level identification method, thereby enabling direct inter- and intra-subject comparison of SC measurements for large cohort studies as well as for longitudinal studies. Our segmentation method is based on the multi-resolution propagation of tubular deformable models. Coupled with an automatic intervertebral disk identification method, our segmentation pipeline provides quantitative metrics of the SC and spinal canal such as cross-sectional areas and volumes in a generic coordinate system based on vertebral levels. This framework was validated on 17 healthy subjects and on one patient with SC injury against manual segmentation. Results have been compared with an existing active surface method and show high local and global accuracy for both SC and spinal canal (Dice coefficients =0.91 ± 0.02) segmentation. Having a robust and automatic framework for SC segmentation and vertebral-based normalization opens the door to bias-free measurement of SC atrophy in large cohorts.