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

Evaluation of the structural integrity of different spinal cord tracts with magnetization transfer ratio in degenerative cervical myelopathy

PURPOSE

Degenerative cervical myelopathy (DCM) is a common cause of spinal cord dysfunction. In this study, we explored the potential of magnetization transfer ratio (MTR) for evaluating the structural integrity of spinal cord tracts in patients with clinically significant DCM.

METHODS

Fifty-three patients with DCM and 41 patients with cervical radiculopathy were evaluated using high-resolution cervical spinal cord magnetic resonance imaging (MRI), which included the magnetization transfer technique. MRI data were analyzed with the Spinal Cord Toolbox (v5.5); MTR values in each spinal tract were calculated and compared between groups after correction for patient age and sex. Correlations between MTR values and patients' clinical disability rate were also evaluated.

RESULTS

A statistically significant reduction in the average MTR of the spinal cord white matter, as well as the MTR of the ventral columns and lateral funiculi, was revealed in the DCM group (adjusted p < 0.01 for all comparisons). Furthermore, reductions in MTR values in the fasciculus cuneatus, spinocerebellar, rubrospinal, and reticulospinal tracts were found in patients with DCM (adjusted p < 0.01 for all comparisons). Positive correlations between the JOA score and the MTR within the ventral columns of the spinal cord (R = 0.38, adjusted p < 0.05) and the ventral spinocerebellar tract (R = 0.41, adjusted p < 0.05) were revealed.

CONCLUSION

The findings of our study indicate that demyelination in patients with DCM primarily affects the spinal tracts of the extrapyramidal system, and the extent of these changes is related to the severity of the condition.

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.

Brainstem and subcortical regions volume loss in patients with degenerative cervical myelopathy and its association with spinal cord compression severity

BACKGROUND

In recent years, structural and functional reorganization of the brain and changes in brainstem structural connectivity have been shown in patients with degenerative cervical myelopathy (DCM). We hypothesized that volume loss in the basal ganglia, thalami, and brainstem structures exists and is associated with spinal cord compression severity in patients with DCM.

METHODS

Forty-seven patients with DCM and 25 patients with cervical radiculopathy were evaluated using cervical spinal cord and brain magnetic resonance imaging (MRI). Brainstem structures, basal ganglia, and thalami volumes were evaluated with FreeSurfer and compared between groups with correction for individual intracranial volume, as well as patient age and sex. Additionally, spinal cord MRI data were analysed with the Spinal Cord Toolbox, and cross-sectional area (CSA) and fractional anisotropy (FA) values were calculated. Correlations between MR-morphometry data and spinal cord structural changes, as well as disease duration, were also evaluated in patients with DCM.

RESULTS

A statistically significant reduction in the volume of the whole brainstem was revealed in the DCM group compared to the radiculopathy group (p < 0.01, FDR-corrected). Additionally, reductions in medulla oblongata, pons and midbrain volumes were found in patients with DCM (p < 0.01, p < 0.01 and p < 0.05, respectively, FDR-corrected). Additionally, a trend in the loss of volume of the left putamen was found (p = 0.087, FDR-corrected). Furthermore, medulla oblongata volume was correlated with spinal cord compression severity (R = 0.54, adjusted p < 0.001) and white matter damage (R = 0.46, adjusted p < 0.05) in patients with DCM. Negative correlations between the duration of the disease and the severity of spinal cord compression (R = -0.42, adjusted p < 0.05) and white matter damage (R = -0.49, adjusted p < 0.05) were also revealed, as well as a trend toward a negative association between the duration of the disease and the volume of the medulla oblongata (R = -0.35; adjusted p < 0.1).

CONCLUSIONS

We revealed a reduction in the volume of brainstem structures in patients with DCM compared to patients with radiculopathy. Moreover, we found that these changes are associated with cord compression severity.

Reliability of multi-parameter mapping (MPM) in the cervical cord: A multi-center multi-vendor quantitative MRI study

MRI based multicenter studies which target neurological pathologies affecting the spinal cord and brain - including spinal cord injury (SCI) - require standardized acquisition protocols and image processing methods. We have optimized and applied a multi-parameter mapping (MPM) protocol that simultaneously covers the brain and the cervical cord within a traveling heads study across six clinical centers (Leutritz et al., 2020). The MPM protocol includes quantitative maps (magnetization transfer saturation (MT), proton density (PD), longitudinal (R1), and effective transverse (R2*) relaxation rates) sensitive to myelination, water content, iron concentration, and morphometric measures, such as cross-sectional cord area. Previously, we assessed the repeatability and reproducibility of the brain MPM data acquired in the five healthy participants who underwent two scan-rescans (Leutritz et al., 2020). This study focuses on the cervical cord MPM data derived from the same acquisitions to determine its repeatability and reproducibility in the cervical cord. MPM matrices of the cervical cord were generated and processed using the hMRI and the spinal cord toolbox. To determine reliability of the cervical MPM data, the intra-site (i.e., scan-rescan) coefficient of variation (CoV), inter-site CoV, and bias within region of interests (C1, C2 and C3 levels) were determined. The range of the mean intra- and inter-site CoV of MT, R1 and PD was between 2.5% and 12%, and between 1.1% and 4.0% for the morphometric measures. In conclusion, the cervical MPM data showed a high repeatability and reproducibility for key imaging biomarkers and hence can be employed as a standardized tool in multi-center studies, including clinical trials.

MAGNIMS recommendations for harmonization of MRI data in MS multicenter studies

•

Sharing data from cooperative studies is essential to develop new biomarkers in MS.

•

Differences in MRI acquisition, analysis, storage represent a substantial constraint.

•

We review the state of the art and developments in the harmonization of MRI.

•

We provide recommendations to harmonize large MRI datasets in the MS field.

There is an increasing need of sharing harmonized data from large, cooperative studies as this is essential to develop new diagnostic and prognostic biomarkers. In the field of multiple sclerosis (MS), the issue has become of paramount importance due to the need to translate into the clinical setting some of the most recent MRI achievements. However, differences in MRI acquisition parameters, image analysis and data storage across sites, with their potential bias, represent a substantial constraint. This review focuses on the state of the art, recent technical advances, and desirable future developments of the harmonization of acquisition, analysis and storage of large-scale multicentre MRI data of MS cohorts. Huge efforts are currently being made to achieve all the requirements needed to provide harmonized MRI datasets in the MS field, as proper management of large imaging datasets is one of our greatest opportunities and challenges in the coming years. Recommendations based on these achievements will be provided here. Despite the advances that have been made, the complexity of these tasks requires further research by specialized academical centres, with dedicated technical and human resources. Such collective efforts involving different professional figures are of crucial importance to offer to MS patients a personalised management while minimizing consumption of resources.

Minimum detectable spinal cord atrophy with automatic segmentation: Investigations using an open-access dataset of healthy participants

•

Evaluate the robustness of an automated analysis pipeline for detecting SC atrophy.

•

Simulate spinal cord atrophy and scan-rescan variability.

•

Fully automated analysis method available on an open access database.

•

Evaluation of sample size and inter/intra-subject variability for T1w and T2w images.

Spinal cord atrophy is a well-known biomarker in multiple sclerosis (MS) and other diseases. It is measured by segmenting the spinal cord on an MRI image and computing the average cross-sectional area (CSA) over a few slices. Introduced about 25 years ago, this procedure is highly sensitive to the quality of the segmentation and is prone to rater-bias. Recently, fully-automated spinal cord segmentation methods, which remove the rater-bias and enable the automated analysis of large populations, have been introduced. A lingering question related to these automated methods is: How reliable are they at detecting atrophy? In this study, we evaluated the precision and accuracy of automated atrophy measurements by simulating scan-rescan experiments.

Spinal cord MRI data from the open-access spine-generic project were used. The dataset aggregates 42 sites worldwide and consists of 260 healthy subjects and includes T1w and T2w contrasts. To simulate atrophy, each volume was globally rescaled at various scaling factors. Moreover, to simulate patient repositioning, random rigid transformations were applied. Using the DeepSeg algorithm from the Spinal Cord Toolbox, the spinal cord was segmented and vertebral levels were identified. Then, the average CSA between C3-C5 vertebral levels was computed for each Monte Carlo sample, allowing us to derive measures of atrophy, intra/inter-subject variability, and sample-size calculations.

The minimum sample size required to detect an atrophy of 2% between unpaired study arms, commonly seen in MS studies, was 467 +/− 13.9 using T1w and 467 +/− 3.2 using T2w images. The minimum sample size to detect a longitudinal atrophy (between paired study arms) of 0.8% was 60 +/− 25.1 using T1w and 10 +/− 1.2 using T2w images. At the intra-subject level, the estimated CSA, observed in this study, showed good precision compared to other studies with COVs (across Monte Carlo transformations) of 0.8% for T1w and 0.6% for T2w images.

While these sample sizes seem small, we would like to stress that these results correspond to a “best case” scenario, in that the dataset used here was of particularly good quality and the model for simulating atrophy does not encompass all the variability met in real-life datasets. The simulated atrophy and scan-rescan variability may over-simplify the biological reality. The proposed framework is open-source and available at https://csa-atrophy.readthedocs.io/.