Multicenter Validation of Mean Upper Cervical Cord Area Measurements from Head 3D T1-Weighted MR Imaging in Patients with Multiple Sclerosis

Exploring subcortical pathology and processing speed in neuromyelitis optica spectrum disorder with myelin water imaging

BACKGROUND AND PURPOSE

Neuromyelitis optica spectrum disorder (NMOSD) affects the optic nerves and spinal cord but can also cause focal brain inflammation. Subcortical pathology may contribute to the etiology of cognitive deficits in NMOSD. Using myelin water imaging, we investigated cerebral normal-appearing white matter (NAWM) and thalamic metrics and their association with cognition in NMOSD participants compared to healthy controls (HC).

METHODS

Seventeen NMOSD participants and 21 HC were scanned on a 3.0-Tesla MRI scanner using a multicomponent driven-equilibrium single-pulse observation of T1 and T2 protocol. Tissue compartment and thalamic volumes (normalized to intracranial volume), T1 relaxation time, and myelin water fraction (MWF) were reported. Eleven NMOSD participants underwent the Symbol Digit Modalities Test (SDMT) for cognitive evaluation. Group comparisons were performed using Student's t-test. The association between thalamic metrics and SDMT score was assessed using multiple regression analysis with age as a covariate.

RESULTS

Compared to HC, NMOSD participants had reduced white matter volume (-14.2%, p < .0001), increased T1 relaxation time (+2.29%, p = .022), and lower MWF (-3.64%, p = .024) in NAWM. NMOSD group had a trend for smaller thalamic volumes than HC (-5.52%, p = .082) and no differences in thalamic MWF (p = .258) or T1 (p = .714). Thalamic T1 predicted SDMT score (adjusted R2 = .51, p = .04) when controlling for age.

CONCLUSIONS

NAWM in NMOSD demonstrates diffuse abnormalities with increased water content and demyelination, suggesting a diffuse disease process overlooked by focal inflammation measures. Increased water content, as a biomarker for diffuse thalamic pathology, may partially explain cognitive impairment in NMOSD.

Association of Mean Upper Cervical Spinal Cord Cross-Sectional Area With Cerebral Small Vessel Disease: A Community-Based Cohort Study

BACKGROUND

The purpose of this study was to investigate the association between the mean upper cervical spinal cord cross-sectional area (MUCCA) and the risk and severity of cerebral small vessel disease (CSVD).

METHODS

Community-dwelling residents in Lishui City, China, from the cross-sectional survey in the PRECISE cohort study (Polyvascular Evaluation for Cognitive Impairment and Vascular Events) conducted from 2017 to 2019. We included 1644 of 3067 community-dwelling adults in the PRECISE study after excluding those with incorrect, incomplete, insufficient, or missing clinical or imaging data. Total and modified total CSVD scores, as well as magnetic resonance imaging features, including white matter hyperintensity, lacunes, cerebral microbleeds, enlarged perivascular spaces, and brain atrophy, were assessed at the baseline. The Spinal Cord Toolbox was used to measure the upper cervical spinal cord cross-sectional area of the C1 to C3 segments of the spinal cord and its average value was taken as MUCCA. Participants were divided into 4 groups according to quartiles of MUCCA. Associations were analyzed using linear regression models adjusted for age, sex, current smoking and drinking, medical history, intracranial volume, and total cortical volume.

RESULTS

The means±SD age of the participants was 61.4±6.5 years, and 635 of 1644 participants (38.6%) were men. The MUCCA was smaller in patients with CSVD than those without CSVD. Using the total CSVD score as a criterion, the MUCCA was 61.78±6.12 cm2 in 504 of 1644 participants with CSVD and 62.74±5.94 cm2 in 1140 of 1644 participants without CSVD. Using the modified total CSVD score, the MUCCA was 61.81±6.04 cm2 in 699 of 1644 participants with CSVD and 62.91±5.94 cm2 in 945 of 1644 without CSVD. There were statistical differences between the 2 groups after adjusting for covariates in 3 models. The MUCCA was negatively associated with the total and modified total CSVD scores (adjusted β value, -0.009 [95% CI, -0.01 to -0.003] and -0.007 [95% CI, -0.01 to -0.0006]) after adjustment for covariates. Furthermore, the MUCCA was negatively associated with the white matter hyperintensity burden (adjusted β value, -0.01 [95% CI, -0.02 to -0.003]), enlarged perivascular spaces in the basal ganglia (adjusted β value, -0.005 [95% CI, -0.009 to -0.001]), lacunes (adjusted β value, -0.004 [95% CI, -0.007 to -0.0007]), and brain atrophy (adjusted β value, -0.009 [95% CI, -0.01 to -0.004]).

CONCLUSIONS

The MUCCA and CSVD were correlated. Spinal cord atrophy may serve as an imaging marker for CSVD; thus, small vessel disease may involve the spinal cord in addition to being intracranial.

MRI assessment of cervical spinal cord cross-sectional area in patients with multiple sclerosis

Objectives

Spinal cord abnormalities including cervical cord atrophy are common in multiple sclerosis (MS). This study aimed to assess the cervical spinal cord cross-sectional area (CSA) using magnetic resonance imaging (MRI) in MS patients.

Materials and Methods

Sixty participants were enrolled in this study (16 male and 44 female), 30 patients with MS, diagnosed according to the revised McDonald criteria, and 30 apparently healthy individuals as the control group. CSA of the spinal cord was measured on axial T2-weighted images of the cervical MRI studies from C2 to C7 vertebral levels.

Results

There was a significant difference between MS patients and the control group in mean CSA at a different level. The mean CSA at C2, in MS cases, was significantly lower than controls (67.7 ± 9.4 mm2 vs. 81.3 ± 4.6 mm2). Similarly, the mean CSA at C7 (64.4 ± 9.9 mm2) and average C2-7 (68 ± 9.1 mm2) of MS cases were significantly lower than the control. There was a strong inverse correlation between mean cervical cord CSA and duration of the disease and disability score. The reduction in cervical cord CSA was more prominent in patients with secondary progressive MS. There was no significant difference regarding age, gender, type of treatment, or the number of cervical cord lesions.

Conclusion

The mean CSA was significantly lower in patients with MS than in the control group and was lesser in progressive types. Patients with a longer duration of MS and a high disability score tend to have smaller CSA.

Validation of a New Semiautomated Segmentation Pipeline Based on the Spinal Cord Toolbox DeepSeg Algorithm to Estimate the Cervical Canal Area

BACKGROUND AND PURPOSE

As in the brain reserve concept, a larger cervical canal area may also protect against disability. In this context, a semiautomated pipeline has been developed to obtain quantitative estimations of the cervical canal area. The aim of the study was to validate the pipeline, to evaluate the consistency of the cervical canal area measurements during a 1-year period, and to compare cervical canal area estimations obtained from brain and cervical MRI acquisitions.

MATERIALS AND METHODS

Eight healthy controls and 18 patients with MS underwent baseline and follow-up 3T brain and cervical spine sagittal 3D MPRAGE. The cervical canal area was measured in all acquisitions, and estimations obtained with the proposed pipeline were compared with manual segmentations performed by 1 evaluator using the Dice similarity coefficient. The cervical canal area estimations obtained on baseline and follow-up T1WI were compared; brain and cervical cord acquisitions were also compared using the individual and average intraclass correlation coefficients.

RESULTS

The agreement between the manual cervical canal area masks and the masks provided by the proposed pipeline was excellent, with a mean Dice similarity coefficient mean of 0.90 (range, 0.73-0.97). The cervical canal area estimations obtained from baseline and follow-up scans showed a good level of concordance (intraclass correlation coefficient = 0.76; 95% CI, 0.44-0.88); estimations obtained from brain and cervical MRIs also had good agreement (intraclass correlation coefficient = 0.77; 95% CI, 0.45-0.90).

CONCLUSIONS

The proposed pipeline is a reliable tool to estimate the cervical canal area. The cervical canal area is a stable measure across time; moreover, when cervical sequences are not available, the cervical canal area could be estimated using brain T1WI.

Spinal cord reserve in multiple sclerosis

BACKGROUND

The spinal cord (SC) is a preferential target of multiple sclerosis (MS) damage highly relevant towards disability. Differential impact of such damage could be due to the initial amount of SC tissue, as described for the brain parenchyma (brain reserve concept). We aimed to test the existence of SC reserve by using spinal canal area (SCaA) as a proxy.

METHODS

Brain sagittal three-dimensional T1-weighted scans covering down to C5 level were acquired in 2930 people with MS and 43 healthy controls (HCs) in a cross-sectional, multicentre study. SC area (SCA) and SCaA were obtained with the Spinal Cord Toolbox. Demographical data and patient-derived disability scores were obtained. SC parameters were compared between groups with age-adjusted and sex-adjusted linear regression models. The main outcome of the study, the existence of an association between SCaA and Patient Determined Disease Steps, was tested with scaled linear models.

RESULTS

1747 persons with MS (mean age: 46.35 years; 73.2% female) and 42 HCs (mean age: 45.56 years; 78.6% female) were analysed after exclusion of post-processing errors and application of quality criteria. SCA (60.41 mm² vs 65.02 mm², p<0.001) was lower in people with MS compared with HC; no differences in SCaA were observed (213.24 mm² vs 212.61 mm², p=0.125). Adjusted scaled linear models showed that a larger SCaA was significantly associated with lower scores on Patient Determined Disease Steps (beta coefficient: -0.12, p=0.0124) independently of spinal cord atrophy, brain T2 lesion volume, age and sex.

CONCLUSIONS

A larger SCaA may be protective against disability in MS, possibly supporting the existence of SC reserve.

Neural correlates of digital measures shown by structural MRI: a post-hoc analysis of a smartphone-based remote assessment feasibility study in multiple sclerosis

BACKGROUND

A study was undertaken to evaluate remote monitoring via smartphone sensor-based tests in people with multiple sclerosis (PwMS). This analysis aimed to explore regional neural correlates of digital measures derived from these tests.

METHODS

In a 24-week, non-randomized, interventional, feasibility study (NCT02952911), sensor-based tests on the Floodlight Proof-of-Concept app were used to assess cognition (smartphone-based electronic Symbol Digit Modalities Test), upper extremity function (Draw a Shape Test, Pinching Test), and gait and balance (Static Balance Test, Two-Minute Walk Test, U-Turn Test). In this post-hoc analysis, digital measures and standard clinical measures (e.g., Nine-Hole Peg Test [9HPT]) were correlated against regional structural magnetic resonance imaging outcomes. Seventy-six PwMS aged 18-55 years with an Expanded Disability Status Scale score of 0.0-5.5 were enrolled from two different sites (USA and Spain). Sixty-two PwMS were included in this analysis.

RESULTS

Worse performance on digital and clinical measures was associated with smaller regional brain volumes and larger ventricular volumes. Whereas digital and clinical measures had many neural correlates in common (e.g., putamen, globus pallidus, caudate nucleus, lateral occipital cortex), some were observed only for digital measures. For example, Draw a Shape Test and Pinching Test measures, but not 9HPT score, correlated with volume of the hippocampus (r = 0.37 [drawing accuracy over time on the Draw a Shape Test]/ - 0.45 [touching asynchrony on the Pinching Test]), thalamus (r = 0.38/ - 0.41), and pons (r = 0.35/ - 0.35).

CONCLUSIONS

Multiple neural correlates were identified for the digital measures in a cohort of people with early MS. Digital measures showed associations with brain regions that clinical measures were unable to demonstrate, thus providing potential novel information on functional ability compared with standard clinical assessments.

Advanced spinal cord MRI in multiple sclerosis: Current techniques and future directions

Spinal cord magnetic resonance imaging (MRI) has a central role in multiple sclerosis (MS) clinical practice for diagnosis and disease monitoring. Advanced MRI sequences capable of visualizing and quantifying tissue macro- and microstructure and reflecting different pathological disease processes have been used in MS research; however, the spinal cord remains under-explored, partly due to technical obstacles inherent to imaging this structure. We propose that the study of the spinal cord merits equal ambition in overcoming technical challenges, and that there is much information to be exploited to make valuable contributions to our understanding of MS. We present a narrative review on the latest progress in advanced spinal cord MRI in MS, covering in the first part structural, functional, metabolic and vascular imaging methods. We focus on recent studies of MS and those making significant technical steps, noting the challenges that remain to be addressed and what stands to be gained from such advances. Throughout we also refer to other works that presend more in-depth review on specific themes. In the second part, we present several topics that, in our view, hold particular potential. The need for better imaging of gray matter is discussed. We stress the importance of developing imaging beyond the cervical spinal cord, and explore the use of ultra-high field MRI. Finally, some recommendations are given for future research, from study design to newer developments in analysis, and the need for harmonization of sequences and methods within the field. This review is aimed at researchers and clinicians with an interest in gaining an overview of the current state of advanced MRI research in this field and what is primed to be the future of spinal cord imaging in MS research.

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.

Cervical Spinal Cord Atrophy can be Accurately Quantified Using Head Images

Background

Spinal cord atrophy provides a clinically relevant metric for monitoring MS. However, the spinal cord is imaged far less frequently than brain due to artefacts and acquisition time, whereas MRI of the brain is routinely performed.

Objective

To validate spinal cord cross-sectional area measurements from routine 3DT1 whole-brain MRI versus those from dedicated cord MRI in healthy controls and people with MS.

Methods

We calculated cross-sectional area at C1 and C2/3 using T2*-weighted spinal cord images and 3DT1 brain images, for 28 healthy controls and 73 people with MS. Correlations for both groups were assessed between: (1) C1 and C2/3 using cord images; (2) C1 from brain and C1 from cord; and (3) C1 from brain and C2/3 from cord.

Results and Conclusion

C1 and C2/3 from cord were strongly correlated in controls (r = 0.94, p<0.0001) and MS (r = 0.85, p<0.0001). There was strong agreement between C1 from brain and C2/3 from cord in controls (r = 0.84, p<0.0001) and MS (r = 0.81, p<0.0001). This supports the use of C1 cross-sectional area calculated from brain imaging as a surrogate for the traditional C2/3 cross-sectional area measure for spinal cord atrophy.

The effect of gadolinium-based contrast-agents on automated brain atrophy measurements by FreeSurfer in patients with multiple sclerosis

OBJECTIVE

To determine whether reliable brain atrophy measures can be obtained from post-contrast 3D T1-weighted images in patients with multiple sclerosis (MS) using FreeSurfer.

METHODS

Twenty-two patients with MS were included, in which 3D T1-weighted MR images were obtained during the same scanner visit, with the same acquisition protocol, before and after administration of gadolinium-based contrast agents (GBCAs). Two FreeSurfer versions (v.6.0.1 and v.7.1.1.) were applied to calculate grey matter (GM) and white matter (WM) volumes and global and regional cortical thickness. The consistency between measures obtained in pre- and post-contrast images was assessed by intra-class correlation coefficient (ICC), the difference was investigated by paired t-tests, and the mean percentage increase or decrease was calculated for total WM and GM matter volume, total deep GM and thalamus volume, and mean cortical thickness.

RESULTS

Good to excellent reliability was found between all investigated measures, with ICC ranging from 0.926 to 0.996, all p values < 0.001. GM volumes and cortical thickness measurements were significantly higher in post-contrast images by 3.1 to 17.4%, while total WM volume decreased significantly by 1.7% (all p values < 0.001).

CONCLUSION

The consistency between values obtained from pre- and post-contrast images was excellent, suggesting it may be possible to extract reliable brain atrophy measurements from T1-weighted images acquired after administration of GBCAs, using FreeSurfer. However, absolute values were systematically different between pre- and post-contrast images, meaning that such images should not be compared directly. Potential systematic effects, possibly dependent on GBCA dose or the delay time after contrast injection, should be investigated.

TRIAL REGISTRATION

Clinical trials.gov. identifier: NCT00360906.

KEY POINTS

• The influence of gadolinium-based contrast agents (GBCAs) on atrophy measurements is still largely unknown and challenges the use of a considerable source of historical and prospective real-world data. • In 22 patients with multiple sclerosis, the consistency between brain atrophy measurements obtained from pre- and post-contrast images was excellent, suggesting it may be possible to extract reliable atrophy measurements in T1-weighted images acquired after administration of GBCAs, using FreeSurfer. • Absolute values were systematically different between pre- and post-contrast images, meaning that such images should not be compared directly, and measurements extracted from certain regions (e.g., the temporal pole) should be interpreted with caution.

Quantification of Cervical Cord Cross-Sectional Area: Which Acquisition, Vertebra Level, and Analysis Software? A Multicenter Repeatability Study on a Traveling Healthy Volunteer

Background: Considerable spinal cord (SC) atrophy occurs in multiple sclerosis (MS). While MRI-based techniques for SC cross-sectional area (CSA) quantification have improved over time, there is no common agreement on whether to measure at single vertebral levels or across larger regions and whether upper SC CSA can be reliably measured from brain images. Aim: To compare in a multicenter setting three CSA measurement methods in terms of repeatability at different anatomical levels. To analyze the agreement between measurements performed on the cervical cord and on brain MRI. Method: One healthy volunteer was scanned three times on the same day in six sites (three scanner vendors) using a 3T MRI protocol including sagittal 3D T1-weighted imaging of the brain (covering the upper cervical cord) and of the SC. Images were analyzed using two semiautomated methods [NeuroQLab (NQL) and the Active Surface Model (ASM)] and the fully automated Spinal Cord Toolbox (SCT) on different vertebral levels (C1-C2; C2/3) on SC and brain images and the entire cervical cord (C1-C7) on SC images only. Results: CSA estimates were significantly smaller using SCT compared to NQL and ASM (p < 0.001), regardless of the cord level. Inter-scanner repeatability was best in C1-C7: coefficients of variation for NQL, ASM, and SCT: 0.4, 0.6, and 1.0%, respectively. CSAs estimated in brain MRI were slightly lower than in SC MRI (all p ≤ 0.006 at the C1-C2 level). Despite protocol harmonization between the centers with regard to image resolution and use of high-contrast 3D T1-weighted sequences, the variability of CSA was partly scanner dependent probably due to differences in scanner geometry, coil design, and details of the MRI parameter settings. Conclusion: For CSA quantification, dedicated isotropic SC MRI should be acquired, which yielded best repeatability in the entire cervical cord. In the upper part of the cervical cord, use of brain MRI scans entailed only a minor loss of CSA repeatability compared to SC MRI. Due to systematic differences between scanners and the CSA quantification software, both should be kept constant within a study. The MRI dataset of this study is available publicly to test new analysis approaches.

A matter of atrophy: differential impact of brain and spine damage on disability worsening in multiple sclerosis

As atrophy represents the most relevant driver of progression in multiple sclerosis (MS), we investigated the impact of different patterns of brain and spinal cord atrophy on disability worsening in MS. We acquired clinical and MRI data from 90 patients with relapsing–remitting MS and 24 healthy controls (HC). Clinical progression at follow-up (mean 3.7 years) was defined according to the Expanded Disability Status Scale-Plus. Brain and spinal cord volumes were computed on MRI brain scans. After normalizing each participants’ brain and spine volume to the mean of the HC, z-score cut-offs were applied to separate pathologically atrophic from normal brain and spine volumes (accepting a 2.5% error probability). Accordingly, MS patients were classified into four groups (Group I: no brain or spinal cord atrophy N = 40, Group II: brain atrophy/no spinal cord atrophy N = 11, Group III: no brain atrophy/ spinal cord atrophy N = 32, Group IV: both brain and spinal cord atrophy N = 7). All patients’ groups showed significantly lower brain volume than HC (p < 0.0001). Group III and IV showed lower spine volume than HC (p < 0.0001 for both). Higher brain lesion load was identified in Group II (p = 0.049) and Group IV (p = 0.023) vs Group I, and in Group IV (p = 0.048) vs Group III. Spinal cord atrophy (OR = 3.75, p = 0.018) and brain + spinal cord atrophy (OR = 5.71, p = 0.046) were significant predictors of disability progression. The presence of concomitant brain and spinal cord atrophy is the strongest correlate of progression over time. Isolated spinal cord atrophy exerts a similar effect, confirming the leading role of spinal cord atrophy in the determination of motor disability.

Cortical volume reductions as a sign of secondary cerebral and cerebellar impairment in patients with degenerative cervical myelopathy

•

Degenerative cervical myelopathy is the most common cause of chronic impairment of the spinal cord.

•

MRI-based anatomical assessment of cerebral and cerebellar areas revealed significant tissue volume reduction in DCM patients compared to healthy controls.

•

Disease severity correlated with cerebral and cerebellar atrophy in the primary motor cortex, primary somatosensory cortex and cerebellar areas.

•

Chronic injury to the spinal cord seems to have impact on remote anatomical structures in the brain.

This study investigated supra- and infratentorial structural gray and white matter (GM, WM) alterations in patients with degenerative cervical myelopathy (DCM) as an indicator of secondary harm due to chronic cervical cord compression and micro trauma. With MRI-based anatomical assessment and subsequent voxel-based morphometry analyses, pre- and postoperative volume alterations in the primary motor cortex (MI), the primary somatosensory cortex (SI), the supplementary motor area (SMA), and the cerebellum were analyzed in 43 DCM patients and 20 controls. We assessed disease-related symptom severity by the modified Japanese Orthopaedic Association scale (mJOA). The study also explored symptom severity-based brain volume alterations as well as their association with clinical status. Patients had lower mJOA scores (p = .000) and lower GM volume than controls in SI (p = .016) and cerebellar regions (p = .001). Symptom severity-based subgroup analyses revealed volume reductions in almost all investigated GM ROIs (MI: p = .001; CB: p = .040; SMA: p = .007) in patients with severe clinical symptoms as well as atrophy already present in patients with moderate symptom severity. Clinical symptoms in DCM were associated with cortical and cerebellar volume reduction. GM volume alterations may serve as an indicator of both disease severity and ongoing disease progression in DCM, and should be considered in further patient care and treatment monitoring.

The sequence of structural, functional and cognitive changes in multiple sclerosis

BACKGROUND

As disease progression remains poorly understood in multiple sclerosis (MS), we aim to investigate the sequence in which different disease milestones occur using a novel data-driven approach.

METHODS

We analysed a cohort of 295 relapse-onset MS patients and 96 healthy controls, and considered 28 features, capturing information on T2-lesion load, regional brain and spinal cord volumes, resting-state functional centrality ("hubness"), microstructural tissue integrity of major white matter (WM) tracts and performance on multiple cognitive tests. We used a discriminative event-based model to estimate the sequence of biomarker abnormality in MS progression in general, as well as specific models for worsening physical disability and cognitive impairment.

RESULTS

We demonstrated that grey matter (GM) atrophy of the cerebellum, thalamus, and changes in corticospinal tracts are early events in MS pathology, whereas other WM tracts as well as the cognitive domains of working memory, attention, and executive function are consistently late events. The models for disability and cognition show early functional changes of the default-mode network and earlier changes in spinal cord volume compared to the general MS population. Overall, GM atrophy seems crucial due to its early involvement in the disease course, whereas WM tract integrity appears to be affected relatively late despite the early onset of WM lesions.

CONCLUSION

Data-driven modelling revealed the relative occurrence of both imaging and non-imaging events as MS progresses, providing insights into disease propagation mechanisms, and allowing fine-grained staging of patients for monitoring purposes.

Spinal cord atrophy in a primary progressive multiple sclerosis trial: Improved sample size using GBSI

•

The GBSI provided clinically meaningful measurements of spinal cord atrophy, with low sample size.

•

Deriving spinal cord atrophy from brain MRI using the GBSI is easier than spinal cord MRI.

•

Spinal cord atrophy on GBSI could be used as a secondary outcome measure.

Background

We aimed to evaluate the implications for clinical trial design of the generalised boundary-shift integral (GBSI) for spinal cord atrophy measurement.

Methods

We included 220 primary-progressive multiple sclerosis patients from a phase 2 clinical trial, with baseline and week-48 3DT1-weighted MRI of the brain and spinal cord (1 × 1 × 1 mm³), acquired separately. We obtained segmentation-based cross-sectional spinal cord area (CSA) at C1-2 (from both brain and spinal cord MRI) and C2-5 levels (from spinal cord MRI) using DeepSeg, and, then, we computed corresponding GBSI.

Results

Depending on the spinal cord segment, we included 67.4–98.1% patients for CSA measurements, and 66.9–84.2% for GBSI. Spinal cord atrophy measurements obtained with GBSI had lower measurement variability, than corresponding CSA. Looking at the image noise floor, the lowest median standard deviation of the MRI signal within the cerebrospinal fluid surrounding the spinal cord was found on brain MRI at the C1-2 level. Spinal cord atrophy derived from brain MRI was related to the corresponding measures from dedicated spinal cord MRI, more strongly for GBSI than CSA. Spinal cord atrophy measurements using GBSI, but not CSA, were associated with upper and lower limb motor progression.

Discussion

Notwithstanding the reduced measurement variability, the clinical correlates, and the possibility of using brain acquisitions, spinal cord atrophy using GBSI should remain a secondary outcome measure in MS studies, until further advancements increase the quality of acquisition and reliability of processing.

Considerations for Mean Upper Cervical Cord Area Implementation in a Longitudinal MRI Setting: Methods, Interrater Reliability, and MRI Quality Control

BACKGROUND AND PURPOSE

Spinal cord atrophy is commonly measured from cerebral MRIs, including the upper cervical cord. However, rescan intraparticipant measures have not been investigated in healthy cohorts. This study investigated technical and rescan variability in the mean upper cervical cord area calculated from T1-weighted cerebral MRIs.

MATERIALS AND METHODS

In this retrospective study, 8 healthy participants were scanned and rescanned with non-distortion- and distortion-corrected MPRAGE sequences (11-50 sessions in 6-8 months), and 50 participants were scanned once with distortion-corrected MPRAGE sequences in the Day2day daily variability study. From another real-world observational cohort, we collected non-distortion-corrected MPRAGE scans from 27 healthy participants (annually for 2-4 years) and cross-sectionally from 77 participants. Statistical analyses included coefficient of variation, smallest real difference, intraclass correlation coefficient, Bland-Altman limits of agreement, and paired t tests.

RESULTS

Distortion- versus non-distortion-corrected MPRAGE-derived mean upper cervical cord areas were similar; however, a paired t test showed incomparability (t = 11.0, P = <.001). Higher variability was found in the mean upper cervical cord areas calculated from an automatic segmentation method. Interrater analysis yielded incomparable measures in the same participant scans (t = 4.5, P = <.001). Non-distortion-corrected mean upper cervical cord area measures were shown to be robust in real-world data (t = -1.04, P = .31). The main sources of variability were found to be artifacts from movement, head/neck positioning, and/or metal implants.

CONCLUSIONS

Technical variability in cord measures decreased using non-distortion-corrected MRIs, a semiautomatic segmentation approach, and 1 rater. Rescan variability was within ±4.4% for group mean upper cervical cord area when MR imaging quality criteria were met.

Clinically relevant cranio-caudal patterns of cervical cord atrophy evolution in MS

OBJECTIVE

To characterize the distribution and regional evolution of cervical cord atrophy in patients with multiple sclerosis (MS) in a multicenter dataset.

METHODS

MRI and clinical evaluations were acquired from 179 controls and 435 patients (35 clinically isolated syndromes [CIS], 259 relapsing-remitting multiple sclerosis [RRMS], 99 secondary progressive multiple sclerosis [SPMS], and 42 primary progressive multiple sclerosis [PPMS]). Sixty-nine controls and 178 patients underwent a 1-year MRI and clinical follow-up. Patients were classified as clinically stable/worsened according to their disability change. Longitudinal changes of cord atrophy were investigated with linear mixed-effect models. Sample size calculations were performed using age-, sex- and site-adjusted annualized percentage normalized cord cross-sectional area (CSAn) changes.

RESULTS

Baseline CSAn was lower in patients with MS vs controls (p < 0.001), but not different between controls and patients with CIS or between patients with early RRMS (disease duration ≤5 years) and patients with CIS. Patients with late RRMS (disease duration >5 years) showed significant cord atrophy vs patients with early RRMS (p = 0.02). Patients with progressive MS had decreased CSAn (p < 0.001) vs patients with RRMS. Atrophy was located between C1/C2 and C5 in patients with RRMS vs patients with CIS, and widespread along the cord in patients with progressive MS vs patients with RRMS, with an additional C5/C6 involvement in patients with SPMS vs patients with PPMS. At follow-up, CSAn decreased in all phenotypes (p < 0.001), except CIS. Cord atrophy rates were highest in patients with early RRMS and clinically worsened patients, who had a more widespread cord involvement than stable patients. The sample size per arm required to detect a 50% treatment effect was 118 for patients with early RRMS.

CONCLUSIONS

Cord atrophy increased in MS during 1 year, except for CIS. Faster atrophy contributed to explain clinical worsening.

Validation of mean upper cervical cord area (MUCCA) measurement techniques in multiple sclerosis (MS): High reproducibility and robustness to lesions, but large software and scanner effects

Introduction

Atrophy of the spinal cord is known to occur in multiple sclerosis (MS). The mean upper cervical cord area (MUCCA) can be used to measure this atrophy. Currently, several (semi-)automated methods for MUCCA measurement exist, but validation in clinical magnetic resonance (MR) images is lacking.

Methods

Five methods to measure MUCCA (SCT-PropSeg, SCT-DeepSeg, NeuroQLab, Xinapse JIM and ITK-SNAP) were investigated in a predefined upper cervical cord region. First, within-scanner reproducibility and between-scanner robustness were assessed using intra-class correlation coefficient (ICC) and Dice's similarity index (SI) in scan-rescan 3DT1-weighted images (brain, including cervical spine using a head coil) performed on three 3 T MR machines (GE MR750, Philips Ingenuity, Toshiba Vantage Titan) in 21 subjects with MS and 6 healthy controls (dataset A). Second, sensitivity of MUCCA measurement to lesions in the upper cervical cord was assessed with cervical 3D T1-weighted images (3 T GE HDxT using a head-neck-spine coil) in 7 subjects with MS without and 14 subjects with MS with cervical lesions (dataset B), using ICC and SI with manual reference segmentations.

Results

In dataset A, MUCCA differed between MR machines (p < 0.001) and methods (p < 0.001) used, but not between scan sessions. With respect to MUCCA values, Xinapse JIM showed the highest within-scanner reproducibility (ICC absolute agreement = 0.995) while Xinapse JIM and SCT-PropSeg showed the highest between-scanner robustness (ICC consistency = 0.981 and 0.976, respectively). Reproducibility of segmentations between scan sessions was highest in Xinapse JIM and SCT-PropSeg segmentations (median SI ≥ 0.921), with a significant main effect of method (p < 0.001), but not of MR machine or subject group. In dataset B, SI with manual outlines did not differ between patients with or without cervical lesions for any of the segmentation methods (p > 0.176). However, there was an effect of method for both volumetric and voxel wise agreement of the segmentations (both p < 0.001). Highest volumetric and voxel wise agreement was obtained with Xinapse JIM (ICC absolute agreement = 0.940 and median SI = 0.962).

Conclusion

Although MUCCA is highly reproducible within a scanner for each individual measurement method, MUCCA differs between scanners and between methods. Cervical cord lesions do not affect MUCCA measurement performance.

•

Mean upper cervical cord area (MUCCA) was obtained with five different methods.

•

MUCCA was determined in a unique scan-rescan multi-vendor MR study.

•

Reproducibility: MUCCA did not differ between scan-rescan images for any method.

•

Robustness: MUCCA differed between methods and between scanners.

•

Performance of MUCCA methods was not affected by the presence of lesions.

Influence of CNS T2-focal lesions on cervical cord atrophy and disability in multiple sclerosis

Background:

Mechanisms associated with cervical spinal cord (CSC) and upper thoracic spinal cord (TSC) atrophy in multiple sclerosis (MS) are poorly understood.

Objective:

To assess the influence of brain, CSC and TSC T2-hyperintense lesions on cord atrophy and disability in MS.

Methods:

Thirty-four MS patients underwent 3T brain, cervical and thoracic cord magnetic resonance imaging (MRI) and Expanded Disability Status Scale (EDSS) score assessment. CSC/TSC lesion number and volume (LV), whole-brain and cortico-spinal tract (CST) LVs were obtained. Normalized whole CSC and upper TSC cross-sectional areas (CSAn) were also derived. Age- and sex-adjusted regression models assessed associations of brain/cord lesions with CSAn and EDSS and identified variables independently associated with CSAn and EDSS with a stepwise variable selection.

Results:

CSC CSAn (β = −0.36, p = 0.03) and TSC CSAn (β = −0.60, p < 0.001) were associated with CSC T2 LV. EDSS (median = 3.0) was correlated with CSC T2 LV (β = 0.42, p = 0.01), brain (β = 0.34, p = 0.04) and CST LV (β = 0.35, p = 0.03). The multivariate analysis retained CSC LV as significant predictor of CSC CSAn ( R2 = 0.20, p = 0.023) and TSC CSAn ( R2 = 0.51, p < 0.001) and retained CSC and CST LVs as significant predictors of EDSS ( R2 = 0.55, p = 0.001).

Conclusions:

CSC LV is an independent predictor of cord atrophy. When neurological impairment is relatively mild, central nervous system (CNS) lesion burden is a better correlate of disability than atrophy.

Advances in spinal cord imaging in multiple sclerosis

The spinal cord is frequently affected in multiple sclerosis (MS), causing motor, sensory and autonomic dysfunction. A number of pathological abnormalities, including demyelination and neuroaxonal loss, occur in the MS spinal cord and are studied in vivo with magnetic resonance imaging (MRI). The aim of this review is to summarise and discuss recent advances in spinal cord MRI. Advances in conventional spinal cord MRI include improved identification of MS lesions, recommended spinal cord MRI protocols, enhanced recognition of MRI lesion characteristics that allow MS to be distinguished from other myelopathies, evidence for the role of spinal cord lesions in predicting prognosis and monitoring disease course, and novel post-processing methods to obtain lesion probability maps. The rate of spinal cord atrophy is greater than that of brain atrophy (-1.78% versus -0.5% per year), and reflects neuroaxonal loss in an eloquent site of the central nervous system, suggesting that it can become an important outcome measure in clinical trials, especially in progressive MS. Recent developments allow the calculation of spinal cord atrophy from brain volumetric scans and evaluation of its progression over time with registration-based techniques. Fully automated analysis methods, including segmentation of grey matter and intramedullary lesions, will facilitate the use of spinal cord atrophy in trial designs and observational studies. Advances in quantitative imaging techniques to evaluate neuroaxonal integrity, myelin content, metabolic changes, and functional connectivity, have provided new insights into the mechanisms of damage in MS. Future directions of research and the possible impact of 7T scanners on spinal cord imaging will be discussed.

MRI-Based Methods for Spinal Cord Atrophy Evaluation: A Comparison of Cervical Cord Cross-Sectional Area, Cervical Cord Volume, and Full Spinal Cord Volume in Patients with Aquaporin-4 Antibody Seropositive Neuromyelitis Optica Spectrum Disorders

BACKGROUND AND PURPOSE

Measures for spinal cord atrophy have become increasingly important as imaging biomarkers in the assessment of neuroinflammatory diseases, especially in neuromyelitis optica spectrum disorders. The most commonly used method, mean upper cervical cord area, is relatively easy to measure and can be performed on brain MRIs that capture cervical myelon. Measures of spinal cord volume (eg, cervical cord volume or total cord volume) require longer scanning and more complex analysis but are potentially better suited as spinal cord atrophy measures. This study investigated spinal cord atrophy measures in a cohort of healthy subjects and patients with aquaporin-4 antibody seropositive neuromyelitis optica spectrum disorders and evaluated the discriminatory performance of mean upper cervical cord cross-sectional area compared with cervical cord volume and total cord volume.

MATERIALS AND METHODS

Mean upper cervical cord area, cervical cord volume, and total cord volume were measured using 3T MRIs from healthy subjects (n = 19) and patients with neuromyelitis optica spectrum disorders (n = 30). Group comparison and receiver operating characteristic analyses between healthy controls and patients with neuromyelitis optica spectrum disorders were performed.

RESULTS

Mean upper cervical cord area, cervical cord volume, and total cord volume measures showed similar and highly significant group differences between healthy control subjects and patients with neuromyelitis optica spectrum disorders (P < .01 for all). All 3 measures showed similar receiver operating characteristic-area under the curve values (mean upper cervical cord area = 0.70, cervical cord volume = 0.75, total cord volume = 0.77) with no significant difference between them. No associations among mean upper cervical cord cross-sectional area, cervical cord volume, or total cord volume with disability measures were found.

CONCLUSIONS

All 3 measures showed similar discriminatory power between healthy control and neuromyelitis optica spectrum disorders groups. Mean upper cervical cord area is easier to obtain compared with cervical cord volume and total cord volume and can be regarded as an efficient representative measure of spinal cord atrophy in the neuromyelitis optica spectrum disorders context.

Future Brain and Spinal Cord Volumetric Imaging in the Clinic for Monitoring Treatment Response in MS

PURPOSE OF REVIEW

Volumetric analysis of brain imaging has emerged as a standard approach used in clinical research, e.g., in the field of multiple sclerosis (MS), but its application in individual disease course monitoring is still hampered by biological and technical limitations. This review summarizes novel developments in volumetric imaging on the road towards clinical application to eventually monitor treatment response in patients with MS.

RECENT FINDINGS

In addition to the assessment of whole-brain volume changes, recent work was focused on the volumetry of specific compartments and substructures of the central nervous system (CNS) in MS. This included volumetric imaging of the deep brain structures and of the spinal cord white and gray matter. Volume changes of the latter indeed independently correlate with clinical outcome measures especially in progressive MS. Ultrahigh field MRI and quantitative MRI added to this trend by providing a better visualization of small compartments on highly resolving MR images as well as microstructural information. New developments in volumetric imaging have the potential to improve sensitivity as well as specificity in detecting and hence monitoring disease-related CNS volume changes in MS.

Cervical spinal cord atrophy: An early marker of progressive MS onset

Objective

To assess whether cervical spinal cord atrophy heralds the onset of progressive MS.

Methods

We studied 34 individuals with radiologically isolated syndrome (RIS) and 31 patients with relapsing-remitting MS (RRMS) age matched to 25 patients within a year of onset of secondary progressive MS (SPMS). Two raters independently measured (twice per rater) the cervical spinal cord average segmental area (CASA) (mm²) of axial T2-weighted images between C2 and C7 landmarks. The midsagittal T2-weighted image from the end of C2 to the end of C7 vertebra was used to measure the cervical spine (c-spine) length (mm). Sex, age at cervical MRI, number and location of cervical spinal cord lesions, c-spine length, and diagnoses were analyzed against the outcome measures of CASA and C2 and C7 slice segmental areas.

Results

Intrarater and interrater agreement was excellent (intraclass correlation coefficient >0.97). The CASA area (p = 0.03) and C7 area (p = 0.002) were smaller in SPMS compared with RRMS. The C2 area (p = 0.027), CASA (p = 0.004), and C7 area (p = 0.003) were smaller in SPMS compared with RIS. The C2 area did not differ between SPMS and RRMS (p = 0.09). The C2 area (p = 0.349), CASA (p = 0.136), and C7 area (p = 0.228) did not differ between RIS and MS (SPMS and RRMS combined). In the multivariable model, ≥2 cervical spinal cord lesions were associated with the C2 area (p = 0.008), CASA (p = 0.009), and C7 area independent of disease course (p = 0.017). Progressive disease course was associated with the C7 area independent of the cervical spinal cord lesion number (p = 0.004).

Conclusion

Cervical spinal cord atrophy is evident at the onset of progressive MS and seems partially independent of the number of cervical spinal cord lesions.

Classification of evidence

This study provides Class III evidence that MRI cervical spinal cord atrophy distinguishes patients at the onset of progressive MS from those with RIS and RRMS.

Gradient nonlinearity effects on upper cervical spinal cord area measurement from 3D T1 -weighted brain MRI acquisitions

PURPOSE

To explore (i) the variability of upper cervical cord area (UCCA) measurements from volumetric brain 3D T₁ -weighted scans related to gradient nonlinearity (GNL) and subject positioning; (ii) the effect of vendor-implemented GNL corrections; and (iii) easily applicable methods that can be used to retrospectively correct data.

METHODS

A multiple sclerosis patient was scanned at seven sites using 3T MRI scanners with the same 3D T₁ -weighted protocol without GNL-distortion correction. Two healthy subjects and a phantom were additionally scanned at a single site with varying table positions. The 2D and 3D vendor-implemented GNL-correction algorithms and retrospective methods based on (i) phantom data fit, (ii) normalization with C2 vertebral body diameters, and (iii) the Jacobian determinant of nonlinear registrations to a template were tested.

RESULTS

Depending on the positioning of the subject, GNL introduced up to 15% variability in UCCA measurements from volumetric brain T₁ -weighted scans when no distortion corrections were used. The 3D vendor-implemented correction methods and the three proposed methods reduced this variability to less than 3%.

CONCLUSIONS

Our results raise awareness of the significant impact that GNL can have on quantitative UCCA studies, and point the way to prospectively and retrospectively managing GNL distortions in a variety of settings, including clinical environments. Magn Reson Med 79:1595-1601, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Elevated CSF neurofilament proteins predict brain atrophy: A 15-year follow-up study

Background:

Body fluid and structural imaging biomarkers give information on neurodegeneration. The relationship over time is not known in multiple sclerosis.

Objective:

To investigate the temporal relationship of elevated cerebrospinal fluid (CSF) neurofilament (Nf) protein levels, a biomarker for axonal loss, with magnetic resonance imaging (MRI) atrophy measures.

Methods:

In patients with multiple sclerosis, CSF Nf heavy chain (NfH) phosphoform levels were quantified at baseline and dichotomised into ‘normal’ and ‘high’. Atrophy was assessed by MRI at baseline and 15-year follow-up using SIENAX and FreeSurfer software.

Results:

High baseline CSF NfH SMI35 levels predicted pronounced atrophy at 15-year follow-up (odds ratio (OR): 36, p < 0.01), in the absence of baseline brain atrophy (OR: 28, p < 0.05), for the averaged MRI normalised brain volume (1.44 L vs 1.33 L, p < 0.05), normalised grey matter volume (0.77 L vs 0.69 L, p < 0.01) and putamen (12.7 mL vs 10.7 mL, p < 0.05). Region-specific calculations including the spinal cord showed that a power of >80% is reached with 14–50 patients.

Conclusion:

These data suggest that high CSF NfH levels are an early predictor of later brain and spinal cord atrophy using structural imaging biomarkers and can be investigated in reasonably sized patient cohorts.