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

Cervical spinal cord morphometrics in degenerative cervical myelopathy: quantification using semi-automated normalized technique and correlation with neurological dysfunctions

BACKGROUND CONTEXT

Degenerative cervical myelopathy (DCM) is characterized by spinal cord atrophy. Accurate estimation of spinal cord atrophy is key to the understanding of neurological diseases, including DCM. However, its clinical application is hampered by difficulties in its precise and consistent estimation due to significant variability in spinal cord morphometry along the cervical spine, both within and between individuals.

PURPOSE

To characterize morphometrics of the compressed spinal cord in DCM patients. We employed our semiautomated analysis framework that incorporates the Spinal Cord Toolbox (SCT) and a normalization approach to effectively address the challenges posed by cord compression in these patients. Additionally, we examined the clinical relevance of these morphometric measures to enhance our understanding of DCM pathophysiology.

STUDY DESIGN

Prospective study.

PATIENT SAMPLE

This study investigated 36 DCM patients and 31 healthy controls (HCs).

OUTCOME MEASURES

Clinical scores including 9-hole peg test for hand dexterity, hand grip strength, balance, gait speed, modified Japanese Orthopaedic Association (mJOA) score, and imaging-based spinal cord morphometrics.

METHOD

Using the generic spine acquisition protocol and our semiautomated analysis pipeline, spinal cord morphometrics, including cross-sectional area (CSA), anterior-posterior (AP) and transverse (RL) diameters, eccentricity, and solidity, were estimated from sagittal T2w magnetic resonance imaging (MRI) images using the Spinal Cord Toolbox (SCT). Normalized metrics were extracted from the C1 to C7 vertebral levels and compared between DCM patients and HC. Morphometric data at regions of maximum spinal cord compression (MSCC) were correlated with the clinical scores. A subset of participants underwent follow-up scans at 6 months to monitor longitudinal changes in spinal cord atrophy.

RESULTS

Spinal cord morphometric data were normalized against the healthy population morphometry (PAM50 database) and extracted for all participants. DCM patients showed a notable reduction in CSA, AP, and RL diameter across all vertebral levels compared to HC. MSCC metrics correlated significantly with clinical scores like dexterity, grip strength, and mJOA scores. Longitudinal analysis indicated a decrease in CSA and worsening clinical scores in DCM patients.

CONCLUSION

Our processing pipeline offers a reliable method for assessing spinal cord compression in DCM patients. Normalized spinal cord morphometrics, particularly the CSA could have potential for monitoring DCM disease severity and progression, guiding treatment decisions. Furthermore, to our knowledge our study is the first to apply the generic spinal cord acquisition protocol, ensuring consistent imaging across different MRI scanners and settings. Coupled with our semiautomated analysis pipeline, this protocol is key for the detailed morphometric characterization of compressed spinal cords in patients with DCM, a disease that is both complex and heterogenous. This study was funded by the National Institute of Neurological Disorders and Stroke (NINDS) (K23:NS091430) and (R01: NS129852-01A1).

Biological variation in cervical spinal cord MRI morphometry in healthy individuals and people with multiple sclerosis

BACKGROUND AND PURPOSE

Conclusions from prior literature regarding the impact of sex, age, and height on spinal cord (SC) MRI morphometrics are conflicting, while the effect of body weight on SC morphometrics has been found to be nonsignificant. The purpose of this case-control study is to assess the associations between cervical SC MRI morphometric parameters and age, sex, height, and weight to establish their potential role as confounding variables in a clinical study of people with multiple sclerosis (MS) compared to a cohort of healthy volunteers.

METHODS

Sixty-nine healthy volunteers and 31 people with MS underwent cervical SC MRI at 3 Tesla field strength. Images were centered at the C3/C4 intervertebral disc and processed using Spinal Cord Toolbox v.4.0.2. Mixed-effects linear regression models were used to evaluate the effects of biological variables and disease status on morphometric parameters.

RESULTS

Sex, age, and height had significant effects on cord and gray matter (GM) cross-sectional area (CSA) as well as the GM:cord CSA ratio. There were no significant effects of body weight on morphometric parameters. The effect of MS disease duration on cord CSA in the C4 level was significant when controlling for all other variables.

CONCLUSIONS

Studies of disease-related changes in SC morphometry should control for sex, age, and height to account for physiological variation.

Cervical and thoracic spinal cord gray matter atrophy is associated with disability in patients with amyotrophic lateral sclerosis

BACKGROUND AND PURPOSE

In amyotrophic lateral sclerosis (ALS), there is an unmet need for more precise patient characterization through quantitative, ideally operator-independent, assessments of disease extent and severity. Radially sampled averaged magnetization inversion recovery acquisitions (rAMIRA) magnetic resonance imaging enables gray matter (GM) and white matter (WM) area quantitation in the cervical and thoracic spinal cord (SC) with optimized contrast. We aimed to investigate rAMIRA-derived SC GM and SC WM areas and their association with clinical phenotype and disability in ALS.

METHODS

A total of 36 patients with ALS (mean [SD] age 61.7 [12.6] years, 14 women) and 36 healthy, age- and sex-matched controls (HCs; mean [SD] age 63.1 [12.1] years, 14 women) underwent two-dimensional axial rAMIRA imaging at the inter-vertebral disc levels C2/3-C5/C6 and the lumbar enlargement level Tmax. ALS Functional Rating Scale-revised (ALSFRS-R) score, muscle strength, and sniff nasal inspiratory pressure (SNIP) were assessed.

RESULTS

Compared to HCs, GM and WM areas were reduced in patients at all cervical levels (p < 0.0001). GM area (p = 0.0001), but not WM area, was reduced at Tmax. Patients with King's Stage 3 showed significant GM atrophy at all levels, while patients with King's Stage 1 showed significant GM atrophy selectively at Tmax. SC GM area was significantly associated with muscle force at corresponding myotomes. GM area at C3/C4 was associated with ALSFRS-R (p < 0.001) and SNIP (p = 0.0016).

CONCLUSION

Patients with ALS assessed by rAMIRA imaging show significant cervical and thoracic SC GM and SC WM atrophy. SC GM area correlates with muscle strength and clinical disability. GM area reduction at Tmax may be an early disease sign. Longitudinal studies are warranted.

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.

Improved inter-subject alignment of the lumbosacral cord for group-level in vivo gray and white matter assessments: A scan-rescan MRI study at 3T

INTRODUCTION

Magnetic resonance imaging (MRI) enables the investigation of pathological changes in gray and white matter at the lumbosacral enlargement (LSE) and conus medullaris (CM). However, conducting group-level analyses of MRI metrics in the lumbosacral spinal cord is challenging due to variability in CM length, lack of established image-based landmarks, and unknown scan-rescan reliability. This study aimed to improve inter-subject alignment of the lumbosacral cord to facilitate group-level analyses of MRI metrics. Additionally, we evaluated the scan-rescan reliability of MRI-based cross-sectional area (CSA) measurements and diffusion tensor imaging (DTI) metrics.

METHODS

Fifteen participants (10 healthy volunteers and 5 patients with spinal cord injury) underwent axial T2*-weighted and diffusion MRI at 3T. We assessed the reliability of spinal cord and gray matter-based landmarks for inter-subject alignment of the lumbosacral cord, the inter-subject variability of MRI metrics before and after adjusting for the CM length, the intra- and inter-rater reliability of CSA measurements, and the scan-rescan reliability of CSA measurements and DTI metrics.

RESULTS

The slice with the largest gray matter CSA as an LSE landmark exhibited the highest reliability, both within and across raters. Adjusting for the CM length greatly reduced the inter-subject variability of MRI metrics. The intra-rater, inter-rater, and scan-rescan reliability of MRI metrics were the highest at and around the LSE (scan-rescan coefficient of variation <3% for CSA measurements and <7% for DTI metrics within the white matter) and decreased considerably caudal to it.

CONCLUSIONS

To facilitate group-level analyses, we recommend using the slice with the largest gray matter CSA as a reliable LSE landmark, along with an adjustment for the CM length. We also stress the significance of the anatomical location within the lumbosacral cord in relation to the reliability of MRI metrics. The scan-rescan reliability values serve as valuable guides for power and sample size calculations in future longitudinal studies.

Pontomedullary junction as a reference for spinal cord cross-sectional area: validation across neck positions

Spinal cord cross-sectional area (CSA) is an important MRI biomarker to assess spinal cord atrophy in various neurodegenerative and traumatic spinal cord diseases. However, the conventional method of computing CSA based on vertebral levels is inherently flawed, as the prediction of spinal levels from vertebral levels lacks reliability, leading to considerable variability in CSA measurements. Computing CSA from an intrinsic neuroanatomical reference, the pontomedullary junction (PMJ), has been proposed in previous work to overcome limitations associated with using a vertebral reference. However, the validation of this alternative approach, along with its variability across and within participants under variable neck extensions, remains unexplored. The goal of this study was to determine if the variability of CSA across neck flexions/extensions is reduced when using the PMJ, compared to vertebral levels. Ten participants underwent a 3T MRI T2w isotropic scan at 0.6 mm3 for 3 neck positions: extension, neutral and flexion. Spinal cord segmentation, vertebral labeling, PMJ labeling, and CSA were computed automatically while spinal segments were labeled manually. Mean coefficient of variation for CSA across neck positions was 3.99 ± 2.96% for the PMJ method vs. 4.02 ± 3.01% for manual spinal segment method vs. 4.46 ± 3.10% for the disc method. These differences were not statistically significant. The PMJ method was slightly more reliable than the disc-based method to compute CSA at specific spinal segments, although the difference was not statistically significant. This suggests that the PMJ can serve as a valuable alternative and reliable method for estimating CSA when a disc-based approach is challenging or not feasible, such as in cases involving fused discs in individuals with spinal cord injuries.

Reliability of spinal cord measures based on synthetic T1-weighted MRI derived from multiparametric mapping (MPM)

Short MRI acquisition time, high signal-to-noise ratio, and high reliability are crucial for image quality when scanning healthy volunteers and patients. Cross-sectional cervical cord area (CSA) has been suggested as a marker of neurodegeneration and potential outcome measure in clinical trials and is conventionally measured on T₁-weigthed 3D Magnetization Prepared Rapid Acquisition Gradient-Echo (MPRAGE) images. This study aims to reduce the acquisition time for the comprehensive assessment of the spinal cord, which is typically based on MPRAGE for morphometry and multi-parameter mapping (MPM) for microstructure. The MPRAGE is replaced by a synthetic T₁-w MRI (synT₁-w) estimated from the MPM, in order to measure CSA. SynT₁-w images were reconstructed using the MPRAGE signal equation based on quantitative maps of proton density (PD), longitudinal (R₁) and effective transverse (R₂*) relaxation rates. The reliability of CSA measurements from synT₁-w images was determined within a multi-center test-retest study format and validated against acquired MPRAGE scans by assessing the agreement between both methods. The response to pathological changes was tested by longitudinally measuring spinal cord atrophy following spinal cord injury (SCI) for synT₁-w and MPRAGE using linear mixed effect models. CSA measurements based on the synT₁-w MRI showed high intra-site (Coefficient of variation [CoV]: 1.43% to 2.71%) and inter-site repeatability (CoV: 2.90% to 5.76%), and only a minor deviation of -1.65 mm² compared to MPRAGE. Crucially, by assessing atrophy rates and by comparing SCI patients with healthy controls longitudinally, differences between synT₁-w and MPRAGE were negligible. These results demonstrate that reliable estimates of CSA can be obtained from synT₁-w images, thereby reducing scan time significantly.

Association of daily physical activity with brain volumes and cervical spinal cord areas in multiple sclerosis

BACKGROUND

Remote activity monitoring has the potential to evaluate real-world, motor function, and disability at home. The relationships of daily physical activity with spinal cord white matter and gray matter (GM) areas, multiple sclerosis (MS) disability and leg function, are unknown.

OBJECTIVE

Evaluate the association of structural central nervous system pathology with ambulatory disability.

METHODS

Fifty adults with progressive or relapsing MS with motor disability who could walk >2 minutes were assessed using clinician-evaluated, patient-reported outcomes, and quantitative brain and spinal cord magnetic resonance imaging (MRI) measures. Fitbit Flex2, worn on the non-dominant wrist, remotely assessed activity over 30 days. Univariate and multivariate analyses were performed to assess correlations between physical activity and other disability metrics.

RESULTS

Mean age was 53.3 years and median Expanded Disability Status Scale (EDSS) was 4.0. Average daily step counts (STEPS) were highly correlated with EDSS and walking measures. Greater STEPS were significantly correlated with greater C2-C3 spinal cord GM areas (ρ = 0.39, p = 0.04), total cord area (TCA; ρ = 0.35, p = 0.04), and cortical GM volume (ρ = 0.32, p = 0.04).

CONCLUSION

These results provide preliminary evidence that spinal cord GM area is a neuroanatomical substrate associated with STEPS. STEPS could serve as a proxy to alert clinicians and researchers to possible changes in structural nervous system pathology.

Association of volumetric MRI measures and disability in MS patients of the same age: Descriptions from a birth year cohort

BACKGROUND AND OBJECTIVES

Although MRI-based markers of neuroinflammation have proven crucial for the diagnosis of multiple sclerosis (MS), predicting clinical progression with inflammation remains difficult. Neurodegenerative markers such as brain volume loss show stronger clinical (predictive) correlations, but also harbor age-related variation that must be disentangled from disease duration. In this study we investigated how clinical disability is related to volumetric MRI measures in a cohort of MS patients and healthy controls (HC) of the same age: Project Y.

METHODS

This study included 234 MS patients born in 1966 and 112 HC born between 1965 and 1967 in the Netherlands. Disability was quantified using the expanded disability status scale (EDSS), nine hole peg test (9HPT), and timed 25 foot walking test (T25FWT). Volumes were quantified on 3T MRI as normalized whole brain (NBV) and regional gray matter (GM) volumes using the same scanner and MRI protocol: cortical (normalized cortical gray matter volume; NCGMV), deep (NDGMV), thalamic (NThalV), and cerebellar (NCbV) GM volumes. In addition, mean upper cervical cord area (MUCCA), white matter lesion volume (LV), and spinal cord lesions were assessed. These measures were compared between patients and HC, and related to disability measures using linear regression.

RESULTS

Mean age of people with MS (PwMS) was 52.8 years (SD 0.9) and median disease duration 15.8 years (IQR 8.7-24.8). All global and regional brain measures were lower in MS patients compared to HC. Univariate regression models showed that NDGMV (β = -0.20) and MUCCA (β = -0.38) were most strongly related to the EDSS in all PwMS. After subtype stratification, MUCCA was most strongly related to the EDSS (β = -0.60) and 9HPT (β = -0.55) in secondary progressive PwMS. Multivariate regression models demonstrated that in all PwMS, the EDSS was best explained by lower MUCCA, longer disease durations and a progressive disease course (adjusted-R (Sastre-Garriga et al., 2017) = 0.26, p < 0.001). MUCCA was a consistent correlate in separate models of the EDSS for all PwMS, relapsing and progressive onset PwMS. The 9HPT (adjusted-R (Sastre-Garriga et al., 2017) = 0.20, p < 0.001) was best explained by lower MUCCA, higher LV and pack years, while lower limb disability (adjusted-R (Sastre-Garriga et al., 2017) = 0.11, p < 0.001) was best explained by lower MUCCA, progressive onset MS and female sex.

DISCUSSION

Our results indicate that in a cohort unbiased by age differences, spinal cord and deep gray matter volumes best related to physical disability. Our results support the use of these measures in clinical practice and trials.

Longitudinal assessment of cervical spinal cord compartments in multiple sclerosis

BACKGROUND

Although cervical spinal cord (cSC) area is an established biomarker in MS, there is currently a lack of longitudinal assessments of cSC gray and white matter areas.

OBJECTIVE

We conducted an explorative analysis of longitudinal changes of cSC gray and white matter areas in MS patients.

METHODS

65 MS patients (33 relapsing-remitting; 20 secondary progressive and 12 primary progressive) and 20 healthy controls (HC) received clinical and upper cSC MRI assessments over 1.10±0.28 years. cSC compartments were quantified on MRI using the novel averaged magnetization inversion recovery acquisitions sequence (in-plane resolution=0.67 × 0.67mm²), and in-house developed post-processing methods. Patients were stratified regarding clinical progression.

RESULTS

Patients with clinical progression showed faster reduction of cSC areas over time at the level of cSC enlargement (approximate vertebral level C4-C5) compared to stable patients (p<0.05). In addition, when compared to the rostral-cSC (approximate vertebral level C2-C3), a preferential reduction of cSC and white matter areas over time at the level of cSC enlargement (p<0.05 and p<0.01, respectively) was demonstrated only in patients with clinical progression, but not in stable MS patients and HC. Compared to HC, MS patients showed comparable changes over time in all cSC compartments.

CONCLUSIONS

MS patients with clinical disease progression demonstrate subtle signs of a more pronounced tissue loss at the level of cSC enlargement. Future studies should consider larger sample sizes and more extended observation periods.

Fully Automatic Method for Reliable Spinal Cord Compartment Segmentation in Multiple Sclerosis

BACKGROUND AND PURPOSE

Fully automatic quantification methods of spinal cord compartments are needed to study pathologic changes of the spinal cord GM and WM in MS in vivo. We propose a novel method for automatic spinal cord compartment segmentation (SCORE) in patients with MS.

MATERIALS AND METHODS

The cervical spinal cords of 24 patients with MS and 24 sex- and age-matched healthy controls were scanned on a 3T MR imaging system, including an averaged magnetization inversion recovery acquisition sequence. Three experienced raters manually segmented the spinal cord GM and WM, anterior and posterior horns, gray commissure, and MS lesions. Subsequently, manual segmentations were used to train neural segmentation networks of spinal cord compartments with multidimensional gated recurrent units in a 3-fold cross-validation fashion. Total intracranial volumes were quantified using FreeSurfer.

RESULTS

The intra- and intersession reproducibility of SCORE was high in all spinal cord compartments (eg, mean relative SD of GM and WM: ≤ 3.50% and ≤1.47%, respectively) and was better than manual segmentations (all P < .001). The accuracy of SCORE compared with manual segmentations was excellent, both in healthy controls and in patients with MS (Dice similarity coefficients of GM and WM: ≥ 0.84 and ≥0.92, respectively). Patients with MS had lower total WM areas (P < .05), and total anterior horn areas (P < .01 respectively), as measured with SCORE.

CONCLUSIONS

We demonstrate a novel, reliable quantification method for spinal cord tissue segmentation in healthy controls and patients with MS and other neurologic disorders affecting the spinal cord. Patients with MS have reduced areas in specific spinal cord tissue compartments, which may be used as MS biomarkers.

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.

Cervical spinal cord atrophy in amyotrophic lateral sclerosis across disease stages

OBJECTIVE

Spinal cord degeneration is a hallmark of amyotrophic lateral sclerosis. The assessment of gray matter and white matter cervical spinal cord atrophy across clinical stages defined using the King's staging system could advance the understanding of amyotrophic lateral sclerosis progression.

METHODS

We assessed the in vivo spatial pattern of gray and white matter atrophy along cervical spinal cord (C2 to C6 segments) using 2D phase-sensitive inversion recovery imaging in a cohort of 44 amyotrophic lateral sclerosis patients, evaluating its change across the King's stages and the correlation with disability scored by the amyotrophic lateral sclerosis functional rating scale revised (ALSFRS-R) and disease duration. A mathematical model inferring the potential onset of cervical gray matter atrophy was developed.

RESULTS

In amyotrophic lateral sclerosis patients at King's stage 1, significant cervical spinal cord alterations were mainly identified in gray matter, whereas they involved both gray and white matter in patients at King's stage ≥ 2. Gray and white matter areas correlated with clinical disability at all cervical segments. C3-C4 level was the segment showing early gray matter atrophy starting about 7 to 20 months before symptom onset according to our model.

INTERPRETATION

Our findings suggest that cervical spinal cord atrophy spreads from gray to white matter across King's stages in amyotrophic lateral sclerosis, making spinal cord magnetic resonance imaging an in vivo assessment tool to measure the progression of the disease.

Dynamics of progressive degeneration of major spinal pathways following spinal cord injury: A longitudinal study

BACKGROUND

Following spinal cord injury (SCI), disease processes spread gradually along the spinal cord forming a spatial gradient with most pronounced changes located at the lesion site. However, the dynamics of this gradient in SCI patients is not established.

OBJECTIVE

This study tracks the spatiotemporal dynamics of remote anterograde and retrograde spinal tract degeneration in the upper cervical cord following SCI over two years utilizing quantitative MRI.

METHODS

Twenty-three acute SCI patients (11 paraplegics, 12 tetraplegics) and 21 healthy controls were scanned with a T1-weighted sequence for volumetry and a FLASH sequence for myelin-sensitive magnetization transfer saturation (MTsat) of the upper cervical cord. We estimated myelin content from MTsat maps within the corticospinal tracts (CST) and dorsal columns (DC) and measured spinal cord atrophy by means of left-right width (LRW) and anterior-posterior width (APW) on the T1-weighted images across cervical levels C1-C3. MTsat in the CST and LRW were considered proxies for retrograde degeneration, while MTsat in the DC and APW provided evidence for anterograde degeneration, respectively. Using regression models, we compared the temporal and spatial trajectories of these MRI readouts between tetraplegics, paraplegics, and controls over a 2-year period and assessed their associations with clinical improvement.

RESULTS

Linear rates and absolute differences in myelin-sensitive MTsat indicated retrograde and anterograde neurodegeneration in the CST and DC, respectively. Changes in MTsat within the CST and in LRW progressively developed over time forming a gradient towards lower cervical levels by 2 years after injury, especially in tetraplegics (change per cervical level in MTsat: -0.247 p.u./level, p = 0.034; in LRW: -0.323 mm/level, p = 0.024). MTsat within the DC was already decreased at cervical levels C1-C3 at baseline (1.5 months after injury) in both tetra- and paraplegics, while linear decreases in APW over time were similar across C1-C3, preserving the spatial gradient. The relative improvement in light touch score was associated with MTsat within the DC at baseline (r_s = 0.575, p = 0.014).

CONCLUSION

Rostral and remote to the injury, the CST and DC show ongoing structural changes, indicative of myelin reductions and atrophy within 2 years after SCI. While anterograde degeneration in the DC was already detectable uniformly at C1-C3 early following SCI, retrograde degeneration in the CST developed over time revealing specific spatial and temporal neurodegenerative gradients. Disentangling and quantifying such dynamic pathological processes may provide biomarkers for regenerative and remyelinating therapies along entire spinal pathways.

Progressive Spinal Cord Degeneration in Friedreich's Ataxia: Results from ENIGMA-Ataxia

BACKGROUND

Spinal cord damage is a hallmark of Friedreich's ataxia (FRDA), but its progression and clinical correlates remain unclear.

OBJECTIVE

The objective of this study was to perform a characterization of cervical spinal cord structural damage in a large multisite FRDA cohort.

METHODS

We performed a cross-sectional analysis of cervical spinal cord (C1-C4) cross-sectional area (CSA) and eccentricity using magnetic resonance imaging data from eight sites within the ENIGMA-Ataxia initiative, including 256 individuals with FRDA and 223 age- and sex-matched control subjects. Correlations and subgroup analyses within the FRDA cohort were undertaken based on disease duration, ataxia severity, and onset age.

RESULTS

Individuals with FRDA, relative to control subjects, had significantly reduced CSA at all examined levels, with large effect sizes (d > 2.1) and significant correlations with disease severity (r < -0.4). Similarly, we found significantly increased eccentricity (d > 1.2), but without significant clinical correlations. Subgroup analyses showed that CSA and eccentricity are abnormal at all disease stages. However, although CSA appears to decrease progressively, eccentricity remains stable over time.

CONCLUSIONS

Previous research has shown that increased eccentricity reflects dorsal column (DC) damage, while decreased CSA reflects either DC or corticospinal tract (CST) damage, or both. Hence our data support the hypothesis that damage to the DC and damage to CST follow distinct courses in FRDA: developmental abnormalities likely define the DC, while CST alterations may be both developmental and degenerative. These results provide new insights about FRDA pathogenesis and indicate that CSA of the cervical spinal cord should be investigated further as a potential biomarker of disease progression. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Anterior horn atrophy in the cervical spinal cord: A new biomarker in progressive multiple sclerosis

BACKGROUND

Spinal cord (SC) gray and white matter pathology plays a central role in multiple sclerosis (MS).

OBJECTIVE

We aimed to investigate the extent, pattern, and clinical relevance of SC gray and white matter atrophy in vivo.

METHODS

39 relapsing-remitting patients (RRMS), 40 progressive MS patients (PMS), and 24 healthy controls (HC) were imaged at 3T using the averaged magnetization inversion recovery acquisitions sequence. Total and lesional cervical gray and white matter, and posterior (SCPH) and anterior horn (SCAH) areas were automatically quantified. Clinical assessment included the expanded disability status scale, timed 25-foot walk test, nine-hole peg test, and the 12-item MS walking scale.

RESULTS

PMS patients had significantly reduced cervical SCAH - but not SCPH - areas compared with HC and RRMS (both p < 0.001). In RRMS and PMS, the cervical SCAH areas increased significantly less in the region of cervical SC enlargement compared with HC (all p < 0.001). This reduction was more pronounced in PMS compared with RRMS (both p < 0.001). In PMS, a lower cervical SCAH area was the most important magnetic resonance imaging (MRI)-variable for higher disability scores.

CONCLUSION

MS patients show clinically relevant cervical SCAH atrophy, which is more pronounced in PMS and at the level of cervical SC enlargement.

Automatic measure and normalization of spinal cord cross-sectional area using the pontomedullary junction

Spinal cord cross-sectional area (CSA) is a relevant biomarker to assess spinal cord atrophy in neurodegenerative diseases. However, the considerable inter-subject variability among healthy participants currently limits its usage. Previous studies explored factors contributing to the variability, yet the normalization models required manual intervention and used vertebral levels as a reference, which is an imprecise prediction of the spinal levels. In this study we implemented a method to measure CSA automatically from a spatial reference based on the central nervous system (the pontomedullary junction, PMJ), we investigated factors to explain variability, and developed normalization strategies on a large cohort (N = 804). Following automatic spinal cord segmentation, vertebral labeling and PMJ labeling, the spinal cord CSA was computed on T1w MRI scans from the UK Biobank database. The CSA was computed using two methods. For the first method, the CSA was computed at the level of the C2-C3 intervertebral disc. For the second method, the CSA was computed at 64 mm caudally from the PMJ, this distance corresponding to the average distance between the PMJ and the C2-C3 disc across all participants. The effect of various demographic and anatomical factors was explored, and a stepwise regression found significant predictors; the coefficients of the best fit model were used to normalize CSA. CSA measured at C2-C3 disc and using the PMJ differed significantly (paired t-test, p-value = 0.0002). The best normalization model included thalamus, brain volume, sex and the interaction between brain volume and sex. The coefficient of variation went down for PMJ CSA from 10.09 (without normalization) to 8.59%, a reduction of 14.85%. For CSA at C2-C3, it went down from 9.96 to 8.42%, a reduction of 15.13 %. This study introduces an end-to-end automatic pipeline to measure and normalize cord CSA from a neurological reference. This approach requires further validation to assess atrophy in longitudinal studies. The inter-subject variability of CSA can be partly accounted for by demographics and anatomical factors.

MRI diffusion tensor imaging scalar values in dogs with intervertebral disc herniation: A comparison between manual and semiautomated region of interest methods

Magnetic resonance imaging (MRI) diffusion tensor imaging (DTI) measures have been described as methods for quantifying spinal cord injury and predicting outcome in dogs with intervertebral disc herniation (IVDH); however, studies comparing methods for selecting regions of interest (ROIs) are currently lacking. The aims of this retrospective, methods comparison, observational study were to compare DTI measurements acquired using manual (mROI) versus semiautomated ROI (sROI) methods and to compare DTI measurements with patient outcomes. Magnetic resonance imaging scans that included DTI pulse sequences were retrieved for 65 dogs with confirmed IVDH. Regions of interest were placed at one vertebral length cranial and caudal to the region of spinal cord compression (RSCC) using the mROI and sROI methods. Scalar values based on the mROI and sROI methods were compared. There was a significant difference for all DTI measures (P < 0.0001), where fractional anisotropy was higher (95% confidence interval [CI]: 0.15, 0.19) and mean diffusivity (MD; CI: -0.41, -0.35), axial diffusivity (AD; CI: -0.47, -0.36) and radial diffusivity (RD; CI: -0.36, -0.27) were lower for the mROI than for the sROI. For both the mROI and sROI, MD, AD, and RD were significantly lower (p < 0.05) at the RSCC in paraplegic dogs that did not regain motor function. The findings indicated that DTI methods for quantifying SCI using open source software and ROI were feasible for use in dogs with IVDH; however, values based on sROI methods differed from values based on mROI methods. Some DTI measures based on both the mROI and sROI methods were predictive of poor patient outcome.

Simultaneous assessment of regional distributions of atrophy across the neuraxis in MS patients

BACKGROUND

The ability to assess brain and cord atrophy simultaneously would improve the efficiency of MRI to track disease evolution.

OBJECTIVE

To test a promising tool to simultaneously map the regional distribution of atrophy in multiple sclerosis (MS) patients across the brain and cord.

METHODS

Voxel-based morphometry combined with a statistical parametric mapping probabilistic brain-spinal cord (SPM-BSC) template was applied to standard T1-weighted magnetic resonance imaging (MRI) scans covering the brain and cervical cord from 37 MS patients and 20 healthy controls (HC). We also measured the cord area at C2-C3 with a semi-automatic segmentation method using (i) the same T1-weighted acquisitions used for the new voxel-based analysis and (ii) dedicated spinal cord phase sensitive inversion recovery (PSIR) acquisitions. Cervical cord findings derived from the three approaches were compared to each other and the goodness to fit to clinical scores was assessed by regression analyses.

RESULTS

The SPM-BSC approach revealed a severity-dependent pattern of atrophy across the cervical cord and thalamus in MS patients when compared to HCs. The magnitude of cord atrophy was confirmed by the semi-automatic extraction approach at C2-C3 using both standard brain T1-weighted and advanced cord dedicated acquisitions. Associations between atrophy of cord and thalamus with disability and cognition were demonstrated.

CONCLUSION

Atrophy in the brain and cervical cord of MS patients can be identified simultaneously and rapidly at the voxel-level. The SPM-BSC approach yields similar results as available standard processing tools with the added advantage of performing the analysis simultaneously and faster.

A New Framework for Investigating the Biological Basis of Degenerative Cervical Myelopathy [AO Spine RECODE-DCM Research Priority Number 5]: Mechanical Stress, Vulnerability and Time

STUDY DESIGN

Literature Review (Narrative).

OBJECTIVE

To propose a new framework, to support the investigation and understanding of the pathobiology of DCM, AO Spine RECODE-DCM research priority number 5.

METHODS

Degenerative cervical myelopathy is a common and disabling spinal cord disorder. In this perspective, we review key knowledge gaps between the clinical phenotype and our biological models. We then propose a reappraisal of the key driving forces behind DCM and an individual's susceptibility, including the proposal of a new framework.

RESULTS

Present pathobiological and mechanistic knowledge does not adequately explain the disease phenotype; why only a subset of patients with visualized cord compression show clinical myelopathy, and the amount of cord compression only weakly correlates with disability. We propose that DCM is better represented as a function of several interacting mechanical forces, such as shear, tension and compression, alongside an individual's vulnerability to spinal cord injury, influenced by factors such as age, genetics, their cardiovascular, gastrointestinal and nervous system status, and time.

CONCLUSION

Understanding the disease pathobiology is a fundamental research priority. We believe a framework of mechanical stress, vulnerability, and time may better represent the disease as a whole. Whilst this remains theoretical, we hope that at the very least it will inspire new avenues of research that better encapsulate the full spectrum of disease.

Extent of cord pathology in the lumbosacral enlargement in non-traumatic versus traumatic spinal cord injury

This study compares remote neurodegenerative changes caudal to a cervical injury in degenerative cervical myelopathy (DCM) (i.e., non-traumatic) and incomplete traumatic spinal cord injury (tSCI) patients, using MRI-based tissue area measurements and diffusion tensor imaging (DTI). Eighteen mild to moderate DCM patients with sensory impairments (mJOA score: 16.2±1.9), 14 incomplete tetraplegic tSCI patients (AIS C&D), and 20 healthy controls were recruited. All participants received DTI and T2*-weighted scans in the lumbosacral enlargement (caudal to injury) and at C2/C3 (rostral to injury). MRI readouts included DTI metrics in the white matter (WM) columns and cross-sectional WM and gray matter area. One-way ANOVA with Tukey's post-hoc comparison (p<0.05) was used to assess group differences. In the lumbosacral enlargement, compared to DCM, tSCI patients exhibited decreased fractional anisotropy in the lateral (tSCI vs. DCM, -11.9%, p=0.007) and ventral WM column (-8.0%, p=0.021), and showed trend toward lower values in the dorsal column (-8.9%, p=0.068). At C2/C3, compared to controls, fractional anisotropy was lower in both groups in the dorsal (DCM vs. controls, -7.9%, p=0.024; tSCI vs. controls, -10.0%, p=0.007) and in the lateral column (DCM: -6.2%, p=0.039; tSCI: -13.3%, p<0.001), while tSCI patients had lower fractional anisotropy than DCM patients in the lateral column (-7.6%, p=0.029). WM areas were not different between patient groups but were lower compared to controls in the lumbosacral enlargement (DCM: -16.9%, p<0.001; tSCI, -10.5%, p=0.043) and at C2/C3 (DCM: -16.0%, p<0.001; tSCI: -18.1%, p<0.001). In conclusion, mild to moderate DCM and incomplete tSCI lead to similar degree of degeneration of the dorsal and lateral columns at C2/C3, but tSCI results in more widespread white matter damage in the lumbosacral enlargement. These remote changes are likely to contribute to the patients' impairment and recovery. DTI is a sensitive tool to assess remote pathological changes in DCM and tSCI patients.

Spinal cord atrophy predicts progressive disease in relapsing multiple sclerosis

Objective

A major challenge in multiple sclerosis (MS) research is the understanding of silent progression and Progressive MS. Using a novel method to accurately capture upper cervical cord area from legacy brain MRI scans we aimed to study the role of spinal cord and brain atrophy for silent progression and conversion to secondary progressive disease (SPMS).

Methods

From a single‐center observational study, all RRMS (n = 360) and SPMS (n = 47) patients and 80 matched controls were evaluated. RRMS patient subsets who converted to SPMS (n = 54) or silently progressed (n = 159), respectively, during the 12‐year observation period were compared to clinically matched RRMS patients remaining RRMS (n = 54) or stable (n = 147), respectively. From brain MRI, we assessed the value of brain and spinal cord measures to predict silent progression and SPMS conversion.

Results

Patients who developed SPMS showed faster cord atrophy rates (−2.19%/yr) at least 4 years before conversion compared to their RRMS matches (−0.88%/yr, p < 0.001). Spinal cord atrophy rates decelerated after conversion (−1.63%/yr, p = 0.010) towards those of SPMS patients from study entry (−1.04%). Each 1% faster spinal cord atrophy rate was associated with 69% (p < 0.0001) and 53% (p < 0.0001) shorter time to silent progression and SPMS conversion, respectively.

Interpretation

Silent progression and conversion to secondary progressive disease are predominantly related to cervical cord atrophy. This atrophy is often present from the earliest disease stages and predicts the speed of silent progression and conversion to Progressive MS. Diagnosis of SPMS is rather a late recognition of this neurodegenerative process than a distinct disease phase. ANN NEUROL 2022;91:268–281

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

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

Normalization of Spinal Cord Total Cross-Sectional and Gray Matter Areas as Quantified With Radially Sampled Averaged Magnetization Inversion Recovery Acquisitions

Background: MR imaging of the spinal cord (SC) gray matter (GM) at the cervical and lumbar enlargements' level may be particularly informative in lower motor neuron disorders, e. g., spinal muscular atrophy, but also in other neurodegenerative or autoimmune diseases affecting the SC. Radially sampled averaged magnetization inversion recovery acquisition (rAMIRA) is a novel approach to perform SC imaging in clinical settings with favorable contrast and is well-suited for SC GM quantitation. However, before applying rAMIRA in clinical studies, it is important to understand (i) the sources of inter-subject variability of total SC cross-sectional areas (TCA) and GM area (GMA) measurements in healthy subjects and (ii) their relation to age and sex to facilitate the detection of pathology-associated changes. In this study, we aimed to develop normalization strategies for rAMIRA-derived SC metrics using skull and spine-based metrics to reduce anatomical variability.

Methods: Sixty-one healthy subjects (age range 11–93 years, 37.7% women) were investigated with axial two-dimensional rAMIRA imaging at 3T MRI. Cervical and thoracic levels including the level of the cervical (C4/C5) and lumbar enlargements (T_max) were examined. SC T2-weighted sagittal images and high-resolution 3D whole-brain T1-weighted images were acquired. TCA and GMAs were quantified. Anatomical variables with associations of |r| > 0.30 in univariate association with SC areas, and age and sex were used to construct normalization models using backward selection with TCA_C4/C5 as outcome. The effect of the normalization was assessed by % relative standard deviation (RSD) reductions.

Results: Mean inter-individual variability and the SD of the SC area metrics were considerable: TCA_C4/5: 8.1%/9.0; TCA_Tmax: 8.9%/6.5; GMA_C4/C5: 8.6%/2.2; GMA_Tmax: 12.2%/3.8. Normalization based on sex, brain WM volume, and spinal canal area resulted in RSD reductions of 23.7% for TCAs and 12.0% for GM areas at C4/C5. Normalizations based on the area of spinal canal alone resulted in RSD reductions of 10.2% for TCAs and 9.6% for GM areas at C4/C5, respectively.

Discussion: Anatomic inter-individual variability of SC areas is substantial. This study identified effective normalization models for inter-subject variability reduction in TCA and SC GMA in healthy subjects based on rAMIRA imaging.

Improving Accuracy of Brainstem MRI Volumetry: Effects of Age and Sex, and Normalization Strategies

Background: Brainstem-mediated functions are impaired in neurodegenerative diseases and aging. Atrophy can be visualized by MRI. This study investigates extrinsic sources of brainstem volume variability, intrinsic sources of anatomical variability, and the influence of age and sex on the brainstem volumes in healthy subjects. We aimed to develop efficient normalization strategies to reduce the effects of intrinsic anatomic variability on brainstem volumetry.

Methods: Brainstem segmentation was performed from MPRAGE data using our deep-learning-based brainstem segmentation algorithm MD-GRU. The extrinsic variability of brainstem volume assessments across scanners and protocols was investigated in two groups comprising 11 (median age 33.3 years, 7 women) and 22 healthy subjects (median age 27.6 years, 50% women) scanned twice and compared using Dice scores. Intrinsic anatomical inter-individual variability and age and sex effects on brainstem volumes were assessed in segmentations of 110 healthy subjects (median age 30.9 years, range 18–72 years, 53.6% women) acquired on 1.5T (45%) and 3T (55%) scanners. The association between brainstem volumes and predefined anatomical covariates was studied using Pearson correlations. Anatomical variables with associations of |r| > 0.30 as well as the variables age and sex were used to construct normalization models using backward selection. The effect of the resulting normalization models was assessed by % relative standard deviation reduction and by comparing the inter-individual variability of the normalized brainstem volumes to the non-normalized values using paired t- tests with Bonferroni correction.

Results: The extrinsic variability of brainstem volumetry across different field strengths and imaging protocols was low (Dice scores > 0.94). Mean inter-individual variability/SD of total brainstem volumes was 9.8%/7.36. A normalization based on either total intracranial volume (TICV), TICV and age, or v-scale significantly reduced the inter-individual variability of total brainstem volumes compared to non-normalized volumes and similarly reduced the relative standard deviation by about 35%.

Conclusion: The extrinsic variability of the novel brainstem segmentation method MD-GRU across different scanners and imaging protocols is very low. Anatomic inter-individual variability of brainstem volumes is substantial. This study presents efficient normalization models for variability reduction in brainstem volumetry in healthy subjects.

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/.

Simultaneous voxel-wise analysis of brain and spinal cord morphometry and microstructure within the SPM framework

To validate a simultaneous analysis tool for the brain and cervical cord embedded in the statistical parametric mapping (SPM) framework, we compared trauma-induced macro- and microstructural changes in spinal cord injury (SCI) patients to controls. The findings were compared with results obtained from existing processing tools that assess the brain and spinal cord separately. A probabilistic brain-spinal cord template (BSC) was generated using a generative semi-supervised modelling approach. The template was incorporated into the pre-processing pipeline of voxel-based morphometry and voxel-based quantification analyses in SPM. This approach was validated on T1-weighted scans and multiparameter maps, by assessing trauma-induced changes in SCI patients relative to controls and comparing the findings with the outcome from existing analytical tools. Consistency of the MRI measures was assessed using intraclass correlation coefficients (ICC). The SPM approach using the BSC template revealed trauma-induced changes across the sensorimotor system in the cord and brain in SCI patients. These changes were confirmed with established approaches covering brain or cord, separately. The ICC in the brain was high within regions of interest, such as the sensorimotor cortices, corticospinal tracts and thalamus. The simultaneous voxel-wise analysis of brain and cervical spinal cord was performed in a unique SPM-based framework incorporating pre-processing and statistical analysis in the same environment. Validation based on a SCI cohort demonstrated that the new processing approach based on the brain and cord is comparable to available processing tools, while offering the advantage of performing the analysis simultaneously across the neuraxis.

MRI Measurement of Upper Cervical Spinal Cord Cross-Sectional Area in Children

BACKGROUND AND PURPOSE

Neurological and neurodegenerative diseases can affect the spinal cord (SC) of pediatric patients. Magnetic resonance imaging (MRI) allows for in vivo quantification of SC atrophy via cross-sectional area (CSA). The study of CSA values in the general population is important to disentangle disease-related changes from intersubject variability. This study aimed at providing normative values for cervical CSA in children, extending our previous work performed with adults.

METHODS

Seventy-eight children (age 7-17 years) were selected from a Developmental Dyslexia study. All subjects underwent a 3T brain MRI session and any incidental findings were reported on the scans. A sagittal 1 mm³ 3-dimensional T₁ -weighted brain acquisition extended to the upper cervical cord was used to measure CSA at C2-C3, as well as spinal canal area and skull volume (V-scale). These three metrics were linearly fitted as a function of age to extract trends and percentage annual changes. Sex differences of CSA were assessed using least squares regression analyses, adjusting for age. We tested normalization strategies proven to be effective in reducing the intersubject variability of adults' CSA.

RESULTS

CSA changed as a function of age at a faster rate when compared with skull volume (CSA: 1.82% increase, V-scale: .60% reduction). Sex had a statistically significant effect on CSA. Normalization methods based on canal area and skull volume reduced the CSA intersubject variability up to 16.84%.

CONCLUSIONS

We present CSA normative values in a large cohort of children, reporting on sources of intersubject variability and how to reduce them applying normalization methods previously developed.