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

Automated Detection and Measurement of Dural Sack Cross-Sectional Area in Lumbar Spine MRI Using Deep Learning

Lumbar spine magnetic resonance imaging (MRI) is a critical diagnostic tool for the assessment of various spinal pathologies, including degenerative disc disease, spinal stenosis, and spondylolisthesis. The accurate identification and quantification of the dural sack cross-sectional area are essential for the evaluation of these conditions. Current manual measurement methods are time-consuming and prone to inter-observer variability. Our study developed and validated deep learning models, specifically U-Net, Attention U-Net, and MultiResUNet, for the automated detection and measurement of the dural sack area in lumbar spine MRI, using a dataset of 515 patients with symptomatic back pain and externally validating the results based on 50 patient scans. The U-Net model achieved an accuracy of 0.9990 and 0.9987 on the initial and external validation datasets, respectively. The Attention U-Net model reported an accuracy of 0.9992 and 0.9989, while the MultiResUNet model displayed a remarkable accuracy of 0.9996 and 0.9995, respectively. All models showed promising precision, recall, and F1-score metrics, along with reduced mean absolute errors compared to the ground truth manual method. In conclusion, our study demonstrates the potential of these deep learning models for the automated detection and measurement of the dural sack cross-sectional area in lumbar spine MRI. The proposed models achieve high-performance metrics in both the initial and external validation datasets, indicating their potential utility as valuable clinical tools for the evaluation of lumbar spine pathologies. Future studies with larger sample sizes and multicenter data are warranted to validate the generalizability of the model further and to explore the potential integration of this approach into routine clinical practice.

Deep Learning-Based Auto-Segmentation of Spinal Cord Internal Structure of Diffusion Tensor Imaging in Cervical Spondylotic Myelopathy

Cervical spondylotic myelopathy (CSM) is a chronic disorder of the spinal cord. ROI-based features on diffusion tensor imaging (DTI) provide additional information about spinal cord status, which would benefit the diagnosis and prognosis of CSM. However, the manual extraction of the DTI-related features on multiple ROIs is time-consuming and laborious. In total, 1159 slices at cervical levels from 89 CSM patients were analyzed, and corresponding fractional anisotropy (FA) maps were calculated. Eight ROIs were drawn, covering both sides of lateral, dorsal, ventral, and gray matter. The UNet model was trained with the proposed heatmap distance loss for auto-segmentation. Mean Dice coefficients on the test dataset for dorsal, lateral, and ventral column and gray matter were 0.69, 0.67, 0.57, 0.54 on the left side and 0.68, 0.67, 0.59, 0.55 on the right side. The ROI-based mean FA value based on segmentation model strongly correlated with the value based on manual drawing. The percentages of the mean absolute error between the two values of multiple ROIs were 0.07, 0.07, 0.11, and 0.08 on the left side and 0.07, 0.1, 0.1, 0.11, and 0.07 on the right side. The proposed segmentation model has the potential to offer a more detailed spinal cord segmentation and would be beneficial for quantifying a more detailed status of the cervical spinal cord.

A pipeline to quantify spinal cord atrophy with deep learning: Application to differentiation of MS and NMOSD patients

PURPOSE

Quantitative measurement of various anatomical regions of the brain and spinal cord (SC) in MRI images are used as unique biomarkers to consider progress and effects of demyelinating diseases of the central nervous system. This paper presents a fully-automated image processing pipeline which quantifies the SC volume of MRI images.

METHODS

In the proposed pipeline, after conducting some pre-processing tasks, a deep convolutional network is utilized to segment the spinal cord cross-sectional area (SCCSA) of each slice. After full segmentation, certain extra slices interpolate between each two adjacent slices using the shape-based interpolation method. Then, a 3D model of the SC is reconstructed, and, by counting the voxels of it, the SC volume is calculated. The performance of the proposed method for the SCCSA segmentation is evaluated on 140 MRI images. Subsequently, to demonstrate the application of the proposed pipeline, we study the differentiations of SC atrophy between 38 Multiple Sclerosis (MS) and 25 Neuromyelitis Optica Spectrum Disorder (NMOSD) patients.

RESULTS

The experimental results of the SCCSA segmentation indicate that the proposed method, adapted by Mask R-CNN, presented the most satisfactory result with the average Dice coefficient of 0.96. For this method, statistical metrics including sensitivity, specificity, accuracy, and precision are 97.51%, 99.98%, 99.92%, and 98.04% respectively. Moreover, the t-test result (p-value = 0.00089) verified a significant difference between the SC atrophy of MS and NMOSD patients.

CONCLUSION

The pipeline efficiently quantifies the SC volume of MRI images and can be utilized as an affordable computer-aided tool for diagnostic purposes.

Voxelwise analysis of diffusion MRI of cervical spinal cord using tract-based spatial statistics

Robust voxelwise analysis using tract-based spatial statistics (TBSS) together with permutation statistical method is standardly used in analyzing diffusion tensor imaging (DTI) of brain. A similar analytical method could be useful when studying DTI of cervical spinal cord. Based on anatomical data of sixty-four healthy volunteers, white (WM) and gray matter (GM) masks were created and subsequently registered into DTI space. Using TBSS, two skeleton types were created (single line and dilated for WM as well as GM). From anatomical data, percentage rates of overlap were calculated for all skeletons in relation to WM and GM masks. Voxelwise analysis of fractional anisotropy values depending on age and sex was conducted. Correlation of fraction anisotropy values with age of subjects was also evaluated. The two WM skeleton types showed a high overlap rate with WM masks (~94%); GM skeletons showed lower rates (56% and 42%, respectively, for single line and dilated). WM and GM areas where fraction anisotropy values differ between sexes were identified (p < .05). Furthermore, using voxelwise analysis such WM voxels were identified where fraction anisotropy values differ depending on age (p < .05) and in these voxels linear dependence of fraction anisotropy and age (r = -0.57, p < .001) was confirmed by regression analysis. This dependence was not proven when using WM anatomical masks (r = -0.21, p = .10). The analytical approach presented shown to be useful for group analysis of DTI data for cervical spinal cord.

Prospective Follow-up of Intramedullary Slitlike Cavities: A Consecutive Series of 48 Patients

Object: Predicting whether intramedullary slitlike cavity (SC) will worsen over time or remain stable is an outstanding clinical challenge. The aim of this study was to identify early features of SC (clinical and magnetic resonance imaging [MRI] findings).

Methods: We prospectively included all patients referred to our institution following the discovery of a SC and divided them in two groups: typical SC (defined as a cavity spanning fewer than three vertebrae, not enlarging the spinal cord, and located at the midline between the anterior third and posterior two-thirds of the spinal cord) or atypical SC (all others). Clinical evolution and changes in MRI features were evaluated during follow-up. In some patients, diffusion tensor imaging was performed and cervical cord cross-sectional area was analyzed.

Results: A total of 48 consecutive patients were included in the study. The mean follow-up was 58 months. Of the seven patients presenting with deficits at first consultation, two worsened and five remained stable. Of the 41 patients without deficits, seven worsened and 34 remained stable. None of the patients developed severe motor deficits or experienced enlargement of the cavity; 7% of patients who presented with typical SC worsened compared with 35% with atypical SC. The negative predictive value was 0.93 (P = 0.02).

Conclusion: Most patients remained stable and a subset of patients developed minor motor deficits. For clinical management, we propose surveillance of patients with a typical SC and close follow-up of those with an atypical SC and/or presenting with deficits.

Fully automatic 3D segmentation of the thoracolumbar spinal cord and the vertebral canal from T2-weighted MRI using K-means clustering algorithm

STUDY DESIGN

Method development.

OBJECTIVES

To develop a reliable protocol for automatic segmentation of Thoracolumbar spinal cord using MRI based on K-means clustering algorithm in 3D images.

SETTING

University-based laboratory, Tehran, Iran.

METHODS

T2 structural volumes acquired from the spinal cord of 20 uninjured volunteers on a 3T MR scanner. We proposed an automatic method for spinal cord segmentation based on the K-means clustering algorithm in 3D images and compare our results with two available segmentation methods (PropSeg, DeepSeg) implemented in the Spinal Cord Toolbox. Dice and Hausdorff were used to compare the results of our method (K-Seg) with the manual segmentation, PropSeg, and DeepSeg.

RESULTS

The accuracy of our automatic segmentation method for T2-weighted images was significantly better or similar to the SCT methods, in terms of 3D DC (p < 0.001). The 3D DCs were respectively (0.81 ± 0.04) and Hausdorff Distance (12.3 ± 2.48) by the K-Seg method in contrary to other SCT methods for T2-weighted images.

CONCLUSIONS

The output with similar protocols showed that K-Seg results match the manual segmentation better than the other methods especially on the thoracolumbar levels in the spinal cord due to the low image contrast as a result of poor SNR in these areas.

Spinal Cord Imaging in Amyotrophic Lateral Sclerosis: Historical Concepts-Novel Techniques

Amyotrophic lateral sclerosis (ALS) is the most common adult onset motor neuron disease with no effective disease modifying therapies at present. Spinal cord degeneration is a hallmark feature of ALS, highlighted in the earliest descriptions of the disease by Lockhart Clarke and Jean-Martin Charcot. The anterior horns and corticospinal tracts are invariably affected in ALS, but up to recently it has been notoriously challenging to detect and characterize spinal pathology in vivo. With recent technological advances, spinal imaging now offers unique opportunities to appraise lower motor neuron degeneration, sensory involvement, metabolic alterations, and interneuron pathology in ALS. Quantitative spinal imaging in ALS has now been used in cross-sectional and longitudinal study designs, applied to presymptomatic mutation carriers, and utilized in machine learning applications. Despite its enormous clinical and academic potential, a number of physiological, technological, and methodological challenges limit the routine use of computational spinal imaging in ALS. In this review, we provide a comprehensive overview of emerging spinal cord imaging methods and discuss their advantages, drawbacks, and biomarker potential in clinical applications, clinical trial settings, monitoring, and prognostic roles.

Tracking a Fast-Moving Disease: Longitudinal Markers, Monitoring, and Clinical Trial Endpoints in ALS

Amyotrophic lateral sclerosis (ALS) encompasses a heterogeneous group of phenotypes with different progression rates, varying degree of extra-motor involvement and divergent progression patterns. The natural history of ALS is increasingly evaluated by large, multi-time point longitudinal studies, many of which now incorporate presymptomatic and post-mortem assessments. These studies not only have the potential to characterize patterns of anatomical propagation, molecular mechanisms of disease spread, but also to identify pragmatic monitoring markers. Sensitive markers of progressive neurodegenerative change are indispensable for clinical trials and individualized patient care. Biofluid markers, neuroimaging indices, electrophysiological markers, rating scales, questionnaires, and other disease-specific instruments have divergent sensitivity profiles. The discussion of candidate monitoring markers in ALS has a dual academic and clinical relevance, and is particularly timely given the increasing number of pharmacological trials. The objective of this paper is to provide a comprehensive and critical review of longitudinal studies in ALS, focusing on the sensitivity profile of established and emerging monitoring markers.

Atlas-Free Cervical Spinal Cord Segmentation on Midsagittal T2-Weighted Magnetic Resonance Images

An automatic atlas-free method for segmenting the cervical spinal cord on midsagittal T2-weighted magnetic resonance images (MRI) is presented. Pertinent anatomical knowledge is transformed into constraints employed at different stages of the algorithm. After picking up the midsagittal image, the spinal cord is detected using expectation maximization and dynamic programming (DP). Using DP, the anterior and posterior edges of the spinal canal and the vertebral column are detected. The vertebral bodies and the intervertebral disks are then segmented using region growing. Then, the anterior and posterior edges of the spinal cord are detected using median filtering followed by DP. We applied this method to 79 noncontrast MRI studies over a 3-month period. The spinal cords were detected in all cases, and the vertebral bodies were successfully labeled in 67 (85%) of them. Our algorithm had very good performance. Compared to manual segmentation results, the Jaccard indices ranged from 0.937 to 1, with a mean of 0.980 ± 0.014. The Hausdorff distances between the automatically detected and manually delineated anterior and posterior spinal cord edges were both 1.0 ± 0.5 mm. Used alone or in combination, our method lays a foundation for computer-aided diagnosis of spinal diseases, particularly cervical spondylotic myelopathy.

Virtual brain biopsies in amyotrophic lateral sclerosis: Diagnostic classification based on in vivo pathological patterns

Background

Diagnostic uncertainty in ALS has serious management implications and delays recruitment into clinical trials. Emerging evidence of presymptomatic disease-burden provides the rationale to develop diagnostic applications based on the evaluation of in-vivo pathological patterns early in the disease.

Objectives

To outline and test a diagnostic classification approach based on an array of complementary imaging metrics in key disease-associated anatomical structures.

Methods

Data from 75 ALS patients and 75 healthy controls were randomly allocated in a ‘training’ and ‘validation’ cohort. Spatial masks were created for anatomical foci which best discriminate patients from controls in the ‘training sample’. In a virtual ‘brain biopsy’, data was then retrieved from these key disease-associated brain regions. White matter diffusivity indices, grey matter T1-signal intensity values and basal ganglia volumes were evaluated as predictor variables in a canonical discriminant function.

Results

Following predictor variable selection, a classification specificity of 85.5% and sensitivity of 89.1% was achieved in the training sample and 90% specificity and 90% sensitivity in the validation sample.

Discussion

This study evaluates disease-associated imaging measures in a dummy diagnostic application. Although larger samples will be required for robust validation, the study confirms the potential of multimodal quantitative imaging in future clinical applications.

•

Reliable diagnostic, monitoring and prognostic biomarkers are urgently in ALS.

•

Accurate diagnostic classification may be achieved based on MRI metrics.

•

Basal ganglia, grey and white matter indices were integrated in a diagnostic model.

•

85.5% specificity and 89.1% sensitivity were achieved in the training sample.

•

90% specificity and 90% sensitivity were achieved in the validation sample.

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.

Cervical Spinal Cord Atrophy Profile in Adult SMN1-Linked SMA

PURPOSE

The mechanisms underlying the topography of motor deficits in spinal muscular atrophy (SMA) remain unknown. We investigated the profile of spinal cord atrophy (SCA) in SMN1-linked SMA, and its correlation with the topography of muscle weakness.

MATERIALS AND METHODS

Eighteen SMN1-linked SMA patients type III/V and 18 age/gender-matched healthy volunteers were included. Patients were scored on manual muscle testing and functional scales. Spinal cord was imaged using 3T MRI system. Radial distance (RD) and cord cross-sectional area (CSA) measurements in SMA patients were compared to those in controls and correlated with strength and disability scores.

RESULTS

CSA measurements revealed a significant cord atrophy gradient mainly located between C3 and C6 vertebral levels with a SCA rate ranging from 5.4% to 23% in SMA patients compared to controls. RD was significantly lower in SMA patients compared to controls in the anterior-posterior direction with a maximum along C4 and C5 vertebral levels (p-values < 10-5). There were no correlations between atrophy measurements, strength and disability scores.

CONCLUSIONS

Spinal cord atrophy in adult SMN1-linked SMA predominates in the segments innervating the proximal muscles. Additional factors such as neuromuscular junction or intrinsic skeletal muscle defects may play a role in more complex mechanisms underlying weakness in these patients.

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.

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.