Accurate quantification methods to evaluate cervical cord atrophy in multiple sclerosis patients

Intersubject Variability and Normalization Strategies for Spinal Cord Total Cross-Sectional and Gray Matter Areas

BACKGROUND AND PURPOSE

The quantification of spinal cord (SC) atrophy by MRI has assumed an important role in assessment of neuroinflammatory/neurodegenerative diseases and traumatic SC injury. Recent technical advances make possible the quantification of gray matter (GM) and white matter tissues in clinical settings. However, the goal of a reliable diagnostic, prognostic or predictive marker is still elusive, in part due to large intersubject variability of SC areas. Here, we investigated the sources of this variability and explored effective strategies to reduce it.

METHODS

One hundred twenty-nine healthy subjects (mean age: 41.0 ± 15.9) underwent MRI on a Siemens 3T Skyra scanner. Two-dimensional PSIR at the C2-C3 vertebral level and a sagittal 1 mm³ 3D T1-weighted brain acquisition extended to the upper cervical cord were acquired. Total cross-sectional area and GM area were measured at C2-C3, as well as measures of the vertebra, spinal canal and the skull. Correlations between the different metrics were explored using Pearson product-moment coefficients. The most promising metrics were used to normalize cord areas using multiple regression analyses.

RESULTS

The most effective normalization metrics were the V-scale (from SienaX) and the product of the C2-C3 spinal canal diameters. Normalization methods based on these metrics reduced the intersubject variability of cord areas of up to 17.74%. The measured cord areas had a statistically significant sex difference, while the effect of age was moderate.

CONCLUSIONS

The present work explored in a large cohort of healthy subjects the source of intersubject variability of SC areas and proposes effective normalization methods for its reduction.

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.

Reliable and fast volumetry of the lumbar spinal cord using cord image analyser (Cordial)

OBJECTIVE

To validate the precision and accuracy of the semi-automated cord image analyser (Cordial) for lumbar spinal cord (SC) volumetry in 3D T1w MRI data of healthy controls (HC).

MATERIALS AND METHODS

40 3D T1w images of 10 HC (w/m: 6/4; age range: 18-41 years) were acquired at one 3T-scanner in two MRI sessions (time interval 14.9±6.1 days). Each subject was scanned twice per session, allowing determination of test-retest reliability both in back-to-back (intra-session) and scan-rescan images (inter-session). Cordial was applied for lumbar cord segmentation twice per image by two raters, allowing for assessment of intra- and inter-rater reliability, and compared to a manual gold standard.

RESULTS

While manually segmented volumes were larger (mean: 2028±245 mm³ vs. Cordial: 1636±300 mm³, p<0.001), accuracy assessments between manually and semi-automatically segmented images showed a mean Dice-coefficient of 0.88±0.05. Calculation of within-subject coefficients of variation (COV) demonstrated high intra-session (1.22-1.86%), inter-session (1.26-1.84%), as well as intra-rater (1.73-1.83%) reproducibility. No significant difference was shown between intra- and inter-session reproducibility or between intra-rater reliabilities. Although inter-rater reproducibility (COV: 2.87%) was slightly lower compared to all other reproducibility measures, between rater consistency was very strong (intraclass correlation coefficient: 0.974).

CONCLUSION

While under-estimating the lumbar SCV, Cordial still provides excellent inter- and intra-session reproducibility showing high potential for application in longitudinal trials.

KEY POINTS

• Lumbar spinal cord segmentation using the semi-automated cord image analyser (Cordial) is feasible. • Lumbar spinal cord is 40-mm cord segment 60 mm above conus medullaris. • Cordial provides excellent inter- and intra-session reproducibility in lumbar spinal cord region. • Cordial shows high potential for application in longitudinal trials.

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.

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.

Fast and accurate semi-automated segmentation method of spinal cord MR images at 3T applied to the construction of a cervical spinal cord template

Objective

To design a fast and accurate semi-automated segmentation method for spinal cord 3T MR images and to construct a template of the cervical spinal cord.

Materials and Methods

A semi-automated double threshold-based method (DTbM) was proposed enabling both cross-sectional and volumetric measures from 3D T₂-weighted turbo spin echo MR scans of the spinal cord at 3T. Eighty-two healthy subjects, 10 patients with amyotrophic lateral sclerosis, 10 with spinal muscular atrophy and 10 with spinal cord injuries were studied. DTbM was compared with active surface method (ASM), threshold-based method (TbM) and manual outlining (ground truth). Accuracy of segmentations was scored visually by a radiologist in cervical and thoracic cord regions. Accuracy was also quantified at the cervical and thoracic levels as well as at C2 vertebral level. To construct a cervical template from healthy subjects’ images (n=59), a standardization pipeline was designed leading to well-centered straight spinal cord images and accurate probability tissue map.

Results

Visual scoring showed better performance for DTbM than for ASM. Mean Dice similarity coefficient (DSC) was 95.71% for DTbM and 90.78% for ASM at the cervical level and 94.27% for DTbM and 89.93% for ASM at the thoracic level. Finally, at C2 vertebral level, mean DSC was 97.98% for DTbM compared with 98.02% for TbM and 96.76% for ASM. DTbM showed similar accuracy compared with TbM, but with the advantage of limited manual interaction.

Conclusion

A semi-automated segmentation method with limited manual intervention was introduced and validated on 3T images, enabling the construction of a cervical spinal cord template.

Validation of a semiautomated spinal cord segmentation method

PURPOSE

To validate semiautomated spinal cord segmentation in healthy subjects and patients with neurodegenerative diseases and trauma.

MATERIALS AND METHODS

Forty-nine healthy subjects, as well as 29 patients with amyotrophic lateral sclerosis, 19 with spinal muscular atrophy, and 14 with spinal cord injuries were studied. Cord area was measured from T2 -weighted 3D turbo spin echo images (cord levels from C2 to T9) using the semiautomated segmentation method of Losseff et al (Brain [1996] 119(Pt 3):701-708), compared with manual segmentation. Reproducibility was evaluated using the inter- and intraobserver coefficient of variation (CoV). Accuracy was assessed using the Dice similarity coefficient (DSC). Robustness to initialization was assessed by simulating modifications to the contours drawn manually prior to segmentation.

RESULTS

Mean interobserver CoV was 4.00% for manual segmentation (1.90% for Losseff's method) in the cervical region and 5.62% (respectively 2.19%) in the thoracic region. Mean intraobserver CoV was 2.34% for manual segmentation (1.08% for Losseff's method) in the cervical region and 2.35% (respectively 1.34%) in the thoracic region. DSC was high (0.96) in both cervical and thoracic regions. DSC remained higher than 0.8 even when modifying initial contours by 50%.

CONCLUSION

The semiautomated segmentation method showed high reproducibility and accuracy in measuring spinal cord area.

Normalized upper cervical spinal cord atrophy in multiple sclerosis

BACKGROUND AND PURPOSE

To find an optimal normalizing factor for upper cervical spinal cord area (UCCA) and to establish whether, in a cross-sectional study, the normalized UCCA correlates better with the neurological disability than the absolute measurement in multiple sclerosis patients.

METHODS

UCCA and three potential normalizing factors were estimated from magnetic resonance imaging data in 51 control subjects. Their reliability was assessed and the linear relationships between UCCA and three potential correction factors were investigated. UCCA was then normalized by these factors respectively. On the basis of these results, an optimal factor was selected and applied to 29 MS subjects.

RESULTS

An extremely strong correlation between UCCA and LECA was seen (r= .88, P < .01). The coefficient of variation (COV) of UCCA was reduced to 4.4% from 9.3% after correction by LECA. The normalized measurement of UCCA correlated better with the expanded disability status scale (EDSS) than the absolute measurement especially in relapsing-remitting multiple sclerosis patients. Moreover, more spinal cord atrophy was identified in corrected data than uncorrected data.

CONCLUSION

Our findings suggest that LECA is an optimal correction factor for UCCA and normalized UCCA is preferable to absolute measurement in cross-sectional study.

MRI in multiple sclerosis: what's inside the toolbox?

Magnetic resonance imaging (MRI) has played a central role in the diagnosis and management of multiple sclerosis (MS). In addition, MRI metrics have become key supportive outcome measures to explore drug efficacy in clinical trials. Conventional MRI measures have contributed to the understanding of MS pathophysiology at the macroscopic level yet have failed to provide a complete picture of underlying MS pathology. They also show relatively weak relationships to clinical status such as predictive strength for clinical progression. Advanced quantitative MRI measures such as magnetization transfer, spectroscopy, diffusion imaging, and relaxometry techniques are somewhat more specific and sensitive for underlying pathology. These measures are particularly useful in revealing diffuse damage in cerebral white and gray matter and therefore may help resolve the dissociation between clinical and conventional MRI findings. In this article, we provide an overview of the array of tools available with brain and spinal cord MRI technology as it is applied to MS. We review the most recent data regarding the role of conventional and advanced MRI techniques in the assessment of MS. We focus on the most relevant pathologic and clinical correlation studies relevant to these measures.