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.

A multicenter pilot study evaluating simplified central vein assessment for the diagnosis of multiple sclerosis

BACKGROUND

The central vein sign (CVS) is a proposed magnetic resonance imaging (MRI) biomarker for multiple sclerosis (MS); the optimal method for abbreviated CVS scoring is not yet established.

OBJECTIVE

The aim of this study was to evaluate the performance of a simplified approach to CVS assessment in a multicenter study of patients being evaluated for suspected MS.

METHODS

Adults referred for possible MS to 10 sites were recruited. A post-Gd 3D T2*-weighted MRI sequence (FLAIR*) was obtained in each subject. Trained raters at each site identified up to six CVS-positive lesions per FLAIR* scan. Diagnostic performance of CVS was evaluated for a diagnosis of MS which had been confirmed using the 2017 McDonald criteria at thresholds including three positive lesions (Select-3*) and six positive lesions (Select-6*). Inter-rater reliability assessments were performed.

RESULTS

Overall, 78 participants were analyzed; 37 (47%) were diagnosed with MS, and 41 (53%) were not. The mean age of participants was 45 (range: 19-64) years, and most were female (n = 55, 71%). The area under the receiver operating characteristic curve (AUROC) for the simplified counting method was 0.83 (95% CI: 0.73-0.93). Select-3* and Select-6* had sensitivity of 81% and 65% and specificity of 68% and 98%, respectively. Inter-rater agreement was 78% for Select-3* and 83% for Select-6*.

CONCLUSION

A simplified method for CVS assessment in patients referred for suspected MS demonstrated good diagnostic performance and inter-rater agreement.

Susceptibility-based imaging aids accurate distinction of pediatric-onset MS from myelin oligodendrocyte glycoprotein antibody-associated disease

BACKGROUND

Myelin oligodendrocyte glycoprotein (MOG) antibody-associated disease (MOGAD) and pediatric-onset multiple sclerosis (POMS) share clinical and magnetic resonance imaging (MRI) features but differ in prognosis and management. Early POMS diagnosis is essential to avoid disability accumulation. Central vein sign (CVS), paramagnetic rim lesions (PRLs), and central core lesions (CCLs) are susceptibility-based imaging (SbI)-related signs understudied in pediatric populations that may help discerning POMS from MOGAD.

METHODS

T2-FLAIR and SbI (three-dimensional echoplanar imaging (3D-EPI)/susceptibility-weighted imaging (SWI) or similar) were acquired on 1.5T/3T scanners. Two readers assessed CVS-positive rate (%CVS+), and their average score was used to build a receiver operator curve (ROC) assessing the ability to discriminate disease type. PRLs and CCLs were identified using a consensual approach.

RESULTS

The %CVS+ distinguished 26 POMS cases (mean age 13.7 years, 63% females, median EDSS 1.5) from 14 MOGAD cases (10.8 years, 35% females, EDSS 1.0) with ROC = 1, p < 0.0001, (cutoff 41%). PRLs were only detectable in POMS participants (mean 2.1±2.3, range 1-10), discriminating the two conditions with a sensitivity of 69% and a specificity of 100%. CCLs were more sensitive (81%) but less specific (71.43%).

CONCLUSION

The %CVS+ and PRLs are highly specific markers of POMS. After proper validation on larger multicenter cohorts, consideration should be given to including such imaging markers for diagnosing POMS at disease onset.

MWF of the corpus callosum is a robust measure of remyelination: Results from the ReBUILD trial

Myelin repair is an unrealized therapeutic goal in the treatment of multiple sclerosis (MS). Uncertainty remains about the optimal techniques for assessing therapeutic efficacy and imaging biomarkers are required to measure and corroborate myelin restoration. We analyzed myelin water fraction imaging from ReBUILD, a double-blind, randomized placebo-controlled (delayed treatment) remyelination trial, that showed a significant reduction in VEP latency in patients with MS. We focused on brain regions rich in myelin. Fifty MS subjects in two arms underwent 3T MRI at baseline and months 3 and 5. Half of the cohort was randomly assigned to receive treatment from baseline through 3 mo, whereas the other half received treatment from 3 mo to 5 mo post-baseline. We computed myelin water fraction changes occurring in normal-appearing white matter of corpus callosum, optic radiations, and corticospinal tracts. An increase in myelin water fraction was documented in the normal-appearing white matter of the corpus callosum, in correspondence with the administration of the remyelinating treatment clemastine. This study provides direct, biologically validated imaging-based evidence of medically induced myelin repair. Moreover, our work strongly suggests that significant myelin repair occurs outside of lesions. We therefore propose myelin water fraction within the normal-appearing white matter of the corpus callosum as a biomarker for clinical trials looking at remyelination.

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.

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.

Effect of GBCA Use on Detection and Diagnostic Performance of the Central Vein Sign: Evaluation Using a 3-T FLAIR* Sequence in Patients With Suspected Multiple Sclerosis

BACKGROUND. The central vein sign (CVS) is a proposed MRI biomarker of multiple sclerosis (MS). The impact of gadolinium-based contrast agent (GBCA) administration on CVS evaluation remains poorly investigated. OBJECTIVE. The purpose of this study was to assess the effect of GBCA use on CVS detection and on the diagnostic performance of the CVS for MS using a 3-T FLAIR* sequence. METHODS. This study was a secondary analysis of data from the pilot study for the prospective multicenter Central Vein Sign: A Diagnostic Biomarker in Multiple Sclerosis (CAVS-MS), which recruited adults with suspected MS from April 2018 to February 2020. Participants underwent 3-T brain MRI including FLAIR and precontrast and post-contrast echo-planar imaging T2*-weighted acquisitions. Postprocessing was used to generate combined FLAIR and T2*-weighted images (hereafter, FLAIR*). MS diagnoses were established using the 2017 McDonald criteria. Thirty participants (23 women, seven men; mean age, 45 years) were randomly selected from the CAVS-MS pilot study cohort. White matter lesions (WMLs) were marked using FLAIR* images. A single observer, blinded to clinical data and GBCA use, reviewed marked WMLs on FLAIR* images for the presence of the CVS. RESULTS. Thirteen of 30 participants had MS. Across participants, on precontrast FLAIR* imaging, 218 CVS-positive and 517 CVS-negative WMLs were identified; on post-contrast FLAIR* imaging, 269 CVS-positive and 459 CVS-negative WMLs were identified. The fraction of WMLs that were CVS-positive on precontrast and postcontrast images was 48% and 58% in participants with MS and 7% and 10% in participants without MS, respectively. The median patient-level CVS-positivity rate on precontrast and postcontrast images was 43% and 67% for participants with MS and 4% and 8% for participants without MS, respectively. In a binomial model adjusting for MS diagnoses, GBCA use was associated with an increased likelihood of at least one CVS-positive WML (odds ratio, 1.6; p < .001). At a 40% CVS-positivity threshold, the sensitivity of the CVS for MS increased from 62% on precontrast images to 92% on postcontrast images (p = .046). Specificity was not significantly different between precontrast (88%) and postcontrast (82%) images (p = .32). CONCLUSION. GBCA use increased CVS detection on FLAIR* images, thereby increasing the sensitivity of the CVS for MS diagnoses. CLINICAL IMPACT. The postcontrast FLAIR* sequence should be considered for CVS evaluation in future investigational trials and clinical practice.

Differences in Age-related Retinal and Cortical Atrophy Rates in Multiple Sclerosis

BACKGROUND AND OBJECTIVES

The timing of neurodegeneration in multiple sclerosis (MS) remains unclear. It is critical to understand the dynamics of neuroaxonal loss if we hope to prevent or forestall permanent disability in MS. We therefore used a deeply phenotyped longitudinal cohort to assess and compare rates of neurodegeneration in retina and brain throughout the MS disease course.

METHODS

We analyzed 597 patients with MS who underwent longitudinal optical coherence tomography imaging annually for 4.5 ± 2.4 years and 432 patients who underwent longitudinal MRI scans for 10 ± 3.4 years, quantifying macular ganglion cell-inner plexiform layer (GCIPL) volume and cortical gray matter (CGM) volume. The association between the slope of decline in the anatomical structure and the age of entry in the cohort (categorized by the MRI cohort's age quartiles) was assessed by hierarchical linear models.

RESULTS

The rate of CGM volume loss declined with increasing age of study entry (1.3% per year atrophy for the age of entry in the cohort younger than 35 years; 1.1% for older than 35 years and younger than 41; 0.97% for older than 41 years and younger than 49; 0.9% for older than 49 years) while the rate of GCIPL thinning was highest in patients in the youngest quartile, fell by more than 50% in the following age quartile, and then stabilized (0.7% per year thinning for the age of entry in the cohort younger than 35 years; 0.29% for age older than 35 and younger than 41 years; 0.34% for older than 41 and younger than 49 years; 0.33% for age older than 49 years).

DISCUSSION

An age-dependent reduction in retinal and cortical volume loss rates during relapsing-remitting MS suggests deceleration in neurodegeneration in the earlier period of disease and further indicates that the period of greatest adaptive immune-mediated inflammatory activity is also the period with the greatest neuroaxonal loss.

Multisite MRI reproducibility of lateral ventricular volume using the NAIMS cooperative pilot dataset

BACKGROUND AND PURPOSE

The North American Imaging in Multiple Sclerosis (NAIMS) multisite project identified interscanner reproducibility issues with T1-based whole brain volume (WBV). Lateral ventricular volume (LVV) acquired on T2-fluid-attenuated inverse recovery (FLAIR) scans has been proposed as a robust proxy measure. Therefore, we sought to determine the relative magnitude of scanner-induced T2-FLAIR-based LVV and T1-based WBV measurement errors in relation to clinically meaningful changes.

METHODS

This was a post hoc analysis of the NAIMS pilot dataset in which a relapsing-remitting MS patient with no intrastudy clinical or radiological activity was imaged twice on seven different Siemens scanners across the United States. LVV was determined using the automated NeuroSTREAM technique on T2-FLAIR and WBV was determined with SIENAX on high-resolution T1-MPRAGE. Average LVV and WBV were measured, and absolute intrascanner and interscanner coefficients of variation (CoVs) were calculated. The variabilities were compared to previously established annual pathological and clinically meaningful cutoffs of 0.40% for WBV and of 3.51% for LVV.

RESULTS

Mean LVV across all seven scan/rescan pairs was 45.87 ± 1.15 ml. Average LVV intrascanner CoV was 1.42% and interscanner CoV was 1.78%, both smaller than the reported annualized clinically meaningful cutoff of 3.51%. In contrast, intra- and interscanner CoVs for WBV (0.99% and 1.15%) were both higher than the established cutoff of 0.40%. Individually, 1/7 intrasite and 2/7 intersite pair-wise LVV comparisons were above the 3.51% cutoff, whereas 4/7 intrasite and 7/7 intersite WBV comparisons were above the 0.40% cutoff.

CONCLUSION

Fully automated LVV segmentation has higher absolute variability than WBV, but much lower relative variability compared to clinically relevant changes, and may therefore be a meaningful proxy outcome measure of neurodegeneration.

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.

Neurite Orientation Dispersion and Density Imaging for Assessing Acute Inflammation and Lesion Evolution in MS

BACKGROUND AND PURPOSE

MR imaging is essential for MS diagnosis and management, yet it has limitations in assessing axonal damage and remyelination. Gadolinium-based contrast agents add value by pinpointing acute inflammation and blood-brain barrier leakage, but with drawbacks in safety and cost. Neurite orientation dispersion and density imaging (NODDI) assesses microstructural features of neurites contributing to diffusion imaging signals. This approach may resolve the components of MS pathology, overcoming conventional MR imaging limitations.

MATERIALS AND METHODS

Twenty-one subjects with MS underwent serial enhanced MRIs (12.6 ± 9 months apart) including NODDI, whose key metrics are the neurite density and orientation dispersion index. Twenty-one age- and sex-matched healthy controls underwent unenhanced MR imaging with the same protocol. Fifty-eight gadolinium-enhancing and non-gadolinium-enhancing lesions were semiautomatically segmented at baseline and follow-up. Normal-appearing WM masks were generated by subtracting lesions and dirty-appearing WM from the whole WM.

RESULTS

The orientation dispersion index was higher in gadolinium-enhancing compared with non-gadolinium-enhancing lesions; logistic regression indicated discrimination, with an area under the curve of 0.73. At follow-up, in the 58 previously enhancing lesions, we identified 2 subgroups based on the neurite density index change across time: Type 1 lesions showed increased neurite density values, whereas type 2 lesions showed decreased values. Type 1 lesions showed greater reduction in size with time compared with type 2 lesions.

CONCLUSIONS

NODDI is a promising tool with the potential to detect acute MS inflammation. The observed heterogeneity among lesions may correspond to gradients in severity and clinical recovery after the acute phase.

Imaging correlates of visual function in multiple sclerosis

No single neuroimaging technique or sequence is capable of reflecting the functional deficits manifest in MS. Given the interest in imaging biomarkers for short- to medium-term studies, we aimed to assess which imaging metrics might best represent functional impairment for monitoring in clinical trials. Given the complexity of functional impairment in MS, however, it is useful to isolate a particular functionally relevant pathway to understand the relationship between imaging and neurological function. We therefore analyzed existing data, combining multiparametric MRI and OCT to describe MS associated visual impairment. We assessed baseline data from fifty MS patients enrolled in ReBUILD, a prospective trial assessing the effect of a remyelinating drug (clemastine). Subjects underwent 3T MRI imaging, including Neurite Orientation Dispersion and Density Imaging (NODDI), myelin content quantification, and retinal imaging, using OCT. Visual function was assessed, using low-contrast letter acuity. MRI and OCT data were studied to model visual function in MS, using a partial, least-squares, regression analysis. Measures of neurodegeneration along the entire visual pathway, described most of the observed variance in visual disability, measured by low contrast letter acuity. In those patients with an identified history of ON, however, putative myelin measures also showed correlation with visual performance. In the absence of clinically identifiable inflammatory episodes, residual disability correlates with neurodegeneration, whereas after an identifiable exacerbation, putative measures of myelin content are additionally informative.

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.

Longitudinal Disconnection Tractograms to Investigate the Functional Consequences of White Matter Damage: An Automated Pipeline

BACKGROUND AND PURPOSE

Neurosurgical resection is one of the few opportunities researchers have to image the human brain pre- and postfocal damage. A major challenge associated with brains undergoing surgical resection is that they often do not fit brain templates most image-processing methodologies are based on. Manual intervention is required to reconcile the pathology, requiring time investment and introducing reproducibility concerns, and extreme cases must be excluded.

METHODS

We propose an automatic longitudinal pipeline based on High Angular Resolution Diffusion Imaging acquisitions to facilitate a Pathway Lesion Symptom Mapping analysis relating focal white matter injury to functional deficits. This two-part approach includes (i) automatic segmentation of focal white matter injury from anisotropic power differences, and (ii) modeling disconnection using tractography on the single-subject level, which specifically identifies the disconnections associated with focal white matter damage.

RESULTS

The advantages of this approach stem from (1) objective and automatic lesion segmentation and tractogram generation, (2) objective and precise segmentation of affected tissue likely to be associated with damage to long-range white matter pathways (defined by anisotropic power), (3) good performance even in the cases of anatomical distortions by use of nonlinear tensor-based registration, which aligns images using an approach sensitive to white matter microstructure.

CONCLUSIONS

Mapping a system as variable and complex as the human brain requires sample sizes much larger than the current technology can support. This pipeline can be used to execute large-scale, sufficiently powered analyses by meeting the need for an automatic approach to objectively quantify white matter disconnection.

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

BACKGROUND AND PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Evaluation of Intra- and Interscanner Reliability of MRI Protocols for Spinal Cord Gray Matter and Total Cross-Sectional Area Measurements

BACKGROUND

In vivo quantification of spinal cord atrophy in neurological diseases using MRI has attracted increasing attention.

PURPOSE

To compare across different platforms the most promising imaging techniques to assess human spinal cord atrophy.

STUDY TYPE

Test/retest multiscanner study.

SUBJECTS

Twelve healthy volunteers.

FIELD STRENGTH/SEQUENCE

Three different 3T scanner platforms (Siemens, Philips, and GE) / optimized phase sensitive inversion recovery (PSIR), T1 -weighted (T1 -w), and T2 *-weighted (T2 *-w) protocols.

ASSESSMENT

On all images acquired, two operators assessed contrast-to-noise ratio (CNR) between gray matter (GM) and white matter (WM), and between WM and cerebrospinal fluid (CSF); one experienced operator measured total cross-sectional area (TCA) and GM area using JIM and the Spinal Cord Toolbox (SCT).

STATISTICAL TESTS

Coefficient of variation (COV); intraclass correlation coefficient (ICC); mixed effect models; analysis of variance (t-tests).

RESULTS

For all the scanners, GM/WM CNR was higher for PSIR than T2 *-w (P < 0.0001) and WM/CSF CNR for T1 -w was the highest (P < 0.0001). For TCA, using JIM, median COVs were smaller than 1.5% and ICC >0.95, while using SCT, median COVs were in the range 2.2-2.75% and ICC 0.79-0.95. For GM, despite some failures of the automatic segmentation, median COVs using SCT on T2 *-w were smaller than using JIM manual PSIR segmentations. In the mixed effect models, the subject was always the main contributor to the variance of area measurements and scanner often contributed to TCA variance (P < 0.05). Using JIM, TCA measurements on T2 *-w were different than on PSIR (P = 0.0021) and T1 -w (P = 0.0018), while using SCT, no notable differences were found between T1 -w and T2 *-w (P = 0.18). JIM and SCT-derived TCA were not different on T1 -w (P = 0.66), while they were different for T2 *-w (P < 0.0001). GM area derived using SCT/T2 *-w versus JIM/PSIR were different (P < 0.0001).

DATA CONCLUSION

The present work sets reference values for the magnitude of the contribution of different effects to cord area measurement intra- and interscanner variability.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2019;49:1078-1090.

Silent progression in disease activity-free relapsing multiple sclerosis

Objective

Rates of worsening and evolution to secondary progressive multiple sclerosis (MS) may be substantially lower in actively treated patients compared to natural history studies from the pretreatment era. Nonetheless, in our recently reported prospective cohort, more than half of patients with relapsing MS accumulated significant new disability by the 10th year of follow‐up. Notably, “no evidence of disease activity” at 2 years did not predict long‐term stability. Here, we determined to what extent clinical relapses and radiographic evidence of disease activity contribute to long‐term disability accumulation.

Methods

Disability progression was defined as an increase in Expanded Disability Status Scale (EDSS) of 1.5, 1.0, or 0.5 (or greater) from baseline EDSS = 0, 1.0–5.0, and 5.5 or higher, respectively, assessed from baseline to year 5 (±1 year) and sustained to year 10 (±1 year). Longitudinal analysis of relative brain volume loss used a linear mixed model with sex, age, disease duration, and HLA‐DRB1*15:01 as covariates.

Results

Relapses were associated with a transient increase in disability over 1‐year intervals (p = 0.012) but not with confirmed disability progression (p = 0.551). Relative brain volume declined at a greater rate among individuals with disability progression compared to those who remained stable (p < 0.05).

Interpretation

Long‐term worsening is common in relapsing MS patients, is largely independent of relapse activity, and is associated with accelerated brain atrophy. We propose the term silent progression to describe the insidious disability that accrues in many patients who satisfy traditional criteria for relapsing–remitting MS. Ann Neurol 2019;85:653–666

Differential intrinsic functional connectivity changes in semantic variant primary progressive aphasia

The semantic variant of primary progressive aphasia (svPPA) is a clinical syndrome characterized by semantic memory deficits with relatively preserved motor speech, syntax, and phonology. There is consistent evidence linking focal neurodegeneration of the anterior temporal lobes (ATL) to the semantic deficits observed in svPPA. Less is known about large-scale functional connectivity changes in this syndrome, particularly regarding the interplay between affected and spared language networks that leads to the unique cognitive dissociations typical of svPPA. Using whole-brain, seed-based connectivity on task-free Magnetic Resonance Imaging (MRI) data, we studied connectivity of networks anchored to three left-hemisphere regions crucially involved in svPPA symptomatology: ATL just posterior to the main atrophic area, opercular inferior frontal gyrus, and posterior inferior temporal lobe. First, in 32 healthy controls, these seeds isolated three networks: a ventral semantic network involving anterior middle temporal and angular gyri, a dorsal articulatory-phonological system involving inferior frontal and supramarginal regions, and a third functional connection between posterior inferior temporal and intraparietal regions likely involved in linking visual and linguistic processes. We then compared connectivity strength of these three networks between 16 svPPA patients and the 32 controls. In svPPA, decreased functional connectivity in the ventral semantic network correlated with weak semantic skills, while connectivity of the network seeded from the posterior inferior temporal lobe, though not significantly different between the two groups, correlated with pseudoword reading skills. Increased connectivity between the inferior frontal gyrus and the superior portion of the angular gyrus suggested possible adaptive changes. Our findings have two main implications. First, they support a functional subdivision of the left IPL based on its connectivity to specific language-related regions. Second, the unique neuroanatomical and linguistic profile observed in svPPA provides a compelling model for the functional interplay of these networks, being either up- or down- regulated in response to disease.

Multisite reliability and repeatability of an advanced brain MRI protocol

BACKGROUND

MRI is the imaging modality of choice for diagnosis and intervention assessment in neurological disease. Its full potential has not been realized due in part to challenges in harmonizing advanced techniques across multiple sites.

PURPOSE

To develop a method for the assessment of reliability and repeatability of advanced multisite-multisession neuroimaging studies and specifically to assess the reliability of an advanced MRI protocol, including multiband fMRI and diffusion tensor MRI, in a multisite setting.

STUDY TYPE

Prospective.

POPULATION

Twice repeated measurement of a single subject with stable relapsing-remitting multiple sclerosis (MS) at seven institutions.

FIELD STRENGTH/SEQUENCE

A 3 T MRI protocol included higher spatial resolution anatomical scans, a variable flip-angle longitudinal relaxation rate constant (R₁ ≡ 1/T₁ ) measurement, quantitative magnetization transfer imaging, diffusion tensor imaging, and a resting-state fMRI (rsFMRI) series.

ASSESSMENT

Multiple methods of assessing intrasite repeatability and intersite reliability were evaluated for imaging metrics derived from each sequence.

STATISTICAL TESTS

Student's t-test, Pearson's r, and intraclass correlation coefficient (ICC) (2,1) were employed to assess repeatability and reliability. Two new statistical metrics are introduced that frame reliability and repeatability in the respective units of the measurements themselves.

RESULTS

Intrasite repeatability was excellent for quantitative R₁ , magnetization transfer ratio (MTR), and diffusion-weighted imaging (DWI) based metrics (r > 0.95). rsFMRI metrics were less repeatable (r = 0.8). Intersite reliability was excellent for R₁ , MTR, and DWI (ICC >0.9), and moderate for rsFMRI metrics (ICC∼0.4).

DATA CONCLUSION

From most reliable to least, using a new reliability metric introduced here, MTR > R₁ > DWI > rsFMRI; for repeatability, MTR > DWI > R₁ > rsFMRI. A graphical method for at-a-glance assessment of reliability and repeatability, effect sizes, and outlier identification in multisite-multisession neuroimaging studies is introduced.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:878-888.

Measurement of spinal cord atrophy using phase sensitive inversion recovery (PSIR) imaging in motor neuron disease

BACKGROUND

The spectrum of motor neuron disease (MND) includes numerous phenotypes with various life expectancies. The degree of upper and lower motor neuron involvement can impact prognosis. Phase sensitive inversion recovery (PSIR) imaging has been shown to detect in vivo gray matter (GM) and white matter (WM) atrophy in the spinal cord of other patient populations but has not been explored in MND.

METHODS

In this study, total cord, WM and GM areas of ten patients with a diagnosis within the MND spectrum were compared to those of ten healthy controls (HC). Patients' diagnosis included amyotrophic lateral sclerosis (ALS), primary lateral sclerosis, primary muscular atrophy, facial onset sensory and motor neuronopathy and ALS-Frontotemporal dementia. Axial 2D PSIR images were acquired at four cervical disc levels (C2-C3, C3-C4, C5-C6 and C7-T1) with a short acquisition time (2 minutes) protocol. Total cross-sectional areas (TCA), GM and WM areas were measured using a combination of highly reliable manual and semi-automated methods. Cord areas in MND patients were compared with HC using linear regression analyses adjusted for age and sex. Correlation of WM and GM areas in MND patients was explored to gain insights into underlying atrophy patterns.

RESULTS

MND patients as a group had significantly smaller cervical cord GM area compared to HC at all four levels (C2-C3: p = .009; C3-C4: p = .001; C5-C6: p = .006; C7-T1: p = .002). WM area at C5-C6 level was significantly smaller (p = .001). TCA was significantly smaller at C3-C4 (p = .018) and C5-C6 (p = .002). No significant GM and WM atrophy was detected in the two patients with predominantly bulbar phenotype. Concomitant GM and WM atrophy was detected in solely upper or lower motor neuron level phenotypes. There was a significant correlation between GM and WM areas at all four levels in this diverse population of MND.

CONCLUSION

Spinal cord GM and WM atrophy can be detected in vivo in patients within the MND spectrum using a short acquisition time 2D PSIR imaging protocol. PSIR imaging shows promise as a method for quantifying spinal cord involvement and thus may be useful for diagnosis, prognosis and for monitoring disease progression.

An Automated Statistical Technique for Counting Distinct Multiple Sclerosis Lesions

BACKGROUND AND PURPOSE

Lesion load is a common biomarker in multiple sclerosis, yet it has historically shown modest association with clinical outcome. Lesion count, which encapsulates the natural history of lesion formation and is thought to provide complementary information, is difficult to assess in patients with confluent (ie, spatially overlapping) lesions. We introduce a statistical technique for cross-sectionally counting pathologically distinct lesions.

MATERIALS AND METHODS

MR imaging was used to assess the probability of a lesion at each location. The texture of this map was quantified using a novel technique, and clusters resembling the center of a lesion were counted. Validity compared with a criterion standard count was demonstrated in 60 subjects observed longitudinally, and reliability was determined using 14 scans of a clinically stable subject acquired at 7 sites.

RESULTS

The proposed count and the criterion standard count were highly correlated (r = 0.97, P < .001) and not significantly different (t₅₉ = -.83, P = .41), and the variability of the proposed count across repeat scans was equivalent to that of lesion load. After accounting for lesion load and age, lesion count was negatively associated (t₅₈ = -2.73, P < .01) with the Expanded Disability Status Scale. Average lesion size had a higher association with the Expanded Disability Status Scale (r = 0.35, P < .01) than lesion load (r = 0.10, P = .44) or lesion count (r = -.12, P = .36) alone.

CONCLUSIONS

This study introduces a novel technique for counting pathologically distinct lesions using cross-sectional data and demonstrates its ability to recover obscured longitudinal information. The proposed count allows more accurate estimation of lesion size, which correlated more closely with disability scores than either lesion load or lesion count alone.

Abnormal age-related cortical folding and neurite morphology in children with developmental dyslexia

There is increasing recognition of a relationship between regional variability in cerebral gyrification and neurodevelopment. Recent work in morphometric MRI has shown that the local gyrification index (lGI), a measure of regional brain folding, may be altered in certain neurodevelopmental disorders. Other studies report that the lGI generally decreases with age in adolescence and young adulthood; however, little is known about how these age-dependent differences in brain maturation occur in atypical neurodevelopment and mechanisms underlying gyrification, such as synaptic pruning. Organization and optimization of dendrites and axons connections across the brain might be driving gyrification and folding processes. In this study, we first assessed lGI differences in the left hemisphere in a cohort of 39 children with developmental dyslexia (DD) between the ages of 7 and 15 years in comparison to 56 typically developing controls (TDC). To better understand the microstructural basis of these changes, we next explored the relationship between lGI differences and cortical thickness and neurite morphology by applying neurite orientation dispersion and density imaging (NODDI). We identified significant differences in lGI between children with DD and TDC in left lateral temporal and middle frontal regions. Further, DD failed to show the expected age-related decreases in lGI in the same regions. Age-related differences in lGI in DD were not explained by differences in cortical thickness, but did correlate with NODDI neurite density and orientation dispersion index. Our findings suggest that gyrification changes in DD are related to abnormal neurite morphology, and are possibly an expression of differences in synaptic pruning.

The NAIMS cooperative pilot project: Design, implementation and future directions

The North American Imaging in Multiple Sclerosis (NAIMS) Cooperative represents a network of 27 academic centers focused on accelerating the pace of magnetic resonance imaging (MRI) research in multiple sclerosis (MS) through idea exchange and collaboration. Recently, NAIMS completed its first project evaluating the feasibility of implementation and reproducibility of quantitative MRI measures derived from scanning a single MS patient using a high-resolution 3T protocol at seven sites. The results showed the feasibility of utilizing advanced quantitative MRI measures in multicenter studies and demonstrated the importance of careful standardization of scanning protocols, central image processing, and strategies to account for inter-site variability.

Volumetric Analysis from a Harmonized Multisite Brain MRI Study of a Single Subject with Multiple Sclerosis

BACKGROUND AND PURPOSE

MR imaging can be used to measure structural changes in the brains of individuals with multiple sclerosis and is essential for diagnosis, longitudinal monitoring, and therapy evaluation. The North American Imaging in Multiple Sclerosis Cooperative steering committee developed a uniform high-resolution 3T MR imaging protocol relevant to the quantification of cerebral lesions and atrophy and implemented it at 7 sites across the United States. To assess intersite variability in scan data, we imaged a volunteer with relapsing-remitting MS with a scan-rescan at each site.

MATERIALS AND METHODS

All imaging was acquired on Siemens scanners (4 Skyra, 2 Tim Trio, and 1 Verio). Expert segmentations were manually obtained for T1-hypointense and T2 (FLAIR) hyperintense lesions. Several automated lesion-detection and whole-brain, cortical, and deep gray matter volumetric pipelines were applied. Statistical analyses were conducted to assess variability across sites, as well as systematic biases in the volumetric measurements that were site-related.

RESULTS

Systematic biases due to site differences in expert-traced lesion measurements were significant (P < .01 for both T1 and T2 lesion volumes), with site explaining >90% of the variation (range, 13.0-16.4 mL in T1 and 15.9-20.1 mL in T2) in lesion volumes. Site also explained >80% of the variation in most automated volumetric measurements. Output measures clustered according to scanner models, with similar results from the Skyra versus the other 2 units.

CONCLUSIONS

Even in multicenter studies with consistent scanner field strength and manufacturer after protocol harmonization, systematic differences can lead to severe biases in volumetric analyses.

Spinal cord microstructure integrating phase-sensitive inversion recovery and diffusional kurtosis imaging

PURPOSE

The aim of this prospective study was to determine the feasibility in terms of repeatability and reproducibility of diffusional kurtosis imaging (DKI) for microstructural assessment of the normal cervical spinal cord (cSC) using a phase-sensitive inversion recovery (PSIR) sequence as the anatomical reference for accurately defining white-matter (WM) and gray-matter (GM) regions of interests (ROIs).

METHODS

Thirteen young healthy subjects were enrolled to undergo DKI and PSIR sequences in the cSC. The repeatability and reproducibility of kurtosis metrics and fractional anisotropy (FA) were calculated in GM, WM, and cerebral-spinal-fluid (CSF) ROIs drawn by two independent readers on PSIR images of three different levels (C1-C4). The presence of statistically significant differences in DKI metrics for levels, ROIs (GM, WM, and CSF) repeatability, reproducibility, and inter-reader agreement was evaluated.

RESULTS

Intra-class correlation coefficients between the two readers ranged from good to excellent (0.75 to 0.90). The inferior level consistently had the highest concordance. The lower values of scan-rescan variability for all DKI parameters were found for the inferior level. Statistically significant differences in kurtosis values were not found in the lateral white-matter bundles of the spinal cord.

CONCLUSION

The integration of DKI and PSIR sequences in a clinical MR acquisition to explore the regional microstructure of the cSC in healthy subjects is feasible, and the results obtainable are reproducible. Further investigation will be required to verify the possibility to translate this method to a clinical setting to study patients with SC involvement especially in the absence of MRI abnormalities on standard sequences.

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.

Spinal cord grey matter segmentation challenge

An important image processing step in spinal cord magnetic resonance imaging is the ability to reliably and accurately segment grey and white matter for tissue specific analysis. There are several semi- or fully-automated segmentation methods for cervical cord cross-sectional area measurement with an excellent performance close or equal to the manual segmentation. However, grey matter segmentation is still challenging due to small cross-sectional size and shape, and active research is being conducted by several groups around the world in this field. Therefore a grey matter spinal cord segmentation challenge was organised to test different capabilities of various methods using the same multi-centre and multi-vendor dataset acquired with distinct 3D gradient-echo sequences. This challenge aimed to characterize the state-of-the-art in the field as well as identifying new opportunities for future improvements. Six different spinal cord grey matter segmentation methods developed independently by various research groups across the world and their performance were compared to manual segmentation outcomes, the present gold-standard. All algorithms provided good overall results for detecting the grey matter butterfly, albeit with variable performance in certain quality-of-segmentation metrics. The data have been made publicly available and the challenge web site remains open to new submissions. No modifications were introduced to any of the presented methods as a result of this challenge for the purposes of this publication.

Investigating Microstructural Abnormalities and Neurocognition in Sub-Acute and Chronic Traumatic Brain Injury Patients with Normal-Appearing White Matter: A Preliminary Diffusion Tensor Imaging Study

For a significant percentage of subjects, with chronic traumatic brain injury (TBI), who report persisting cognitive impairment and functional loss, the diagnosis is often impeded by the fact that routine neuroimaging often does not reveal any abnormalities. In this paper, we used diffusion tensor imaging (DTI) to investigate the apparently normal white matter (as assessed by routine magnetic resonance imaging) in the brains of 19 subjects with sub-acute (9) and chronic (10) TBI. We also assessed memory, executive function, and visual-motor coordination in these subjects. Using a voxel-wise approach, we investigated if parameters of diffusion were significantly different between TBI subjects and 17 healthy controls (HC), who were demographically matched to the TBI group. We also investigated if changes in DTI parameters were associated with neuropsychological performance in either group. Our results indicate significantly increased mean and axial diffusivity (MD and AD, respectively) values in widespread brain locations in TBI subjects, while controlling for age, sex, and time since injury. HC performed significantly better than the TBI subjects on tests of memory and executive function, indicating the persisting functional loss in chronic TBI. We found no correlation between diffusion parameters and performance on test of executive function in either group. We found negative correlation between FA and composite memory scores, and positive correlation between RD and visuomotor coordination test scores, in various tracts in both groups. Our study suggests that changes in MD and AD can indicate persisting micro-structure abnormalities in normal-appearing white matter in the brains of subjects with chronic TBI. Our results also suggest that FA in major white matter tracts is correlated with memory in health and in disease, alike; larger and longitudinal studies are needed to discern potential differences in these correlations in the two groups.

Long-term evolution of multiple sclerosis disability in the treatment era

Objective

To characterize the accrual of long‐term disability in a cohort of actively treated multiple sclerosis (MS) patients and to assess whether clinical and magnetic resonance imaging (MRI) data used in clinical trials have long‐term prognostic value.

Methods

This is a prospective study of 517 actively managed MS patients enrolled at a single center.

Results

More than 91% of patients were retained, with data ascertained up to 10 years after the baseline visit. At this last assessment, neurologic disability as measured by the Expanded Disability Status Scale (EDSS) was stable or improved compared to baseline in 41% of patients. Subjects with no evidence of disease activity (NEDA) by clinical and MRI criteria during the first 2 years had long‐term outcomes that were no different from those of the cohort as a whole. 25‐OH vitamin D serum levels were inversely associated with short‐term MS disease activity; however, these levels had no association with long‐term disability. At a median time of 16.8 years after disease onset, 10.7% (95% confidence interval [CI] = 7.2–14%) of patients reached an EDSS ≥ 6, and 18.1% (95% CI = 13.5–22.5%) evolved from relapsing MS to secondary progressive MS (SPMS).

Interpretation

Rates of worsening and evolution to SPMS were substantially lower when compared to earlier natural history studies. Notably, the NEDA 2‐year endpoint was not a predictor of long‐term stability. Finally, the data call into question the utility of annual MRI assessments as a treat‐to‐target approach for MS care. Ann Neurol 2016;80:499–510

Gray matter segmentation of the spinal cord with active contours in MR images

OBJECTIVE

Fully or partially automated spinal cord gray matter segmentation techniques for spinal cord gray matter segmentation will allow for pivotal spinal cord gray matter measurements in the study of various neurological disorders. The objective of this work was multi-fold: (1) to develop a gray matter segmentation technique that uses registration methods with an existing delineation of the cord edge along with Morphological Geodesic Active Contour (MGAC) models; (2) to assess the accuracy and reproducibility of the newly developed technique on 2D PSIR T1 weighted images; (3) to test how the algorithm performs on different resolutions and other contrasts; (4) to demonstrate how the algorithm can be extended to 3D scans; and (5) to show the clinical potential for multiple sclerosis patients.

METHODS

The MGAC algorithm was developed using a publicly available implementation of a morphological geodesic active contour model and the spinal cord segmentation tool of the software Jim (Xinapse Systems) for initial estimate of the cord boundary. The MGAC algorithm was demonstrated on 2D PSIR images of the C2/C3 level with two different resolutions, 2D T2* weighted images of the C2/C3 level, and a 3D PSIR image. These images were acquired from 45 healthy controls and 58 multiple sclerosis patients selected for the absence of evident lesions at the C2/C3 level. Accuracy was assessed though visual assessment, Hausdorff distances, and Dice similarity coefficients. Reproducibility was assessed through interclass correlation coefficients. Validity was assessed through comparison of segmented gray matter areas in images with different resolution for both manual and MGAC segmentations.

RESULTS

Between MGAC and manual segmentations in healthy controls, the mean Dice similarity coefficient was 0.88 (0.82-0.93) and the mean Hausdorff distance was 0.61 (0.46-0.76) mm. The interclass correlation coefficient from test and retest scans of healthy controls was 0.88. The percent change between the manual segmentations from high and low-resolution images was 25%, while the percent change between the MGAC segmentations from high and low resolution images was 13%. Between MGAC and manual segmentations in MS patients, the average Dice similarity coefficient was 0.86 (0.8-0.92) and the average Hausdorff distance was 0.83 (0.29-1.37) mm.

CONCLUSION

We demonstrate that an automatic segmentation technique, based on a morphometric geodesic active contours algorithm, can provide accurate and precise spinal cord gray matter segmentations on 2D PSIR images. We have also shown how this automated technique can potentially be extended to other imaging protocols.

Power estimation for non-standardized multisite studies

A concern for researchers planning multisite studies is that scanner and T1-weighted sequence-related biases on regional volumes could overshadow true effects, especially for studies with a heterogeneous set of scanners and sequences. Current approaches attempt to harmonize data by standardizing hardware, pulse sequences, and protocols, or by calibrating across sites using phantom-based corrections to ensure the same raw image intensities. We propose to avoid harmonization and phantom-based correction entirely. We hypothesized that the bias of estimated regional volumes is scaled between sites due to the contrast and gradient distortion differences between scanners and sequences. Given this assumption, we provide a new statistical framework and derive a power equation to define inclusion criteria for a set of sites based on the variability of their scaling factors. We estimated the scaling factors of 20 scanners with heterogeneous hardware and sequence parameters by scanning a single set of 12 subjects at sites across the United States and Europe. Regional volumes and their scaling factors were estimated for each site using Freesurfer's segmentation algorithm and ordinary least squares, respectively. The scaling factors were validated by comparing the theoretical and simulated power curves, performing a leave-one-out calibration of regional volumes, and evaluating the absolute agreement of all regional volumes between sites before and after calibration. Using our derived power equation, we were able to define the conditions under which harmonization is not necessary to achieve 80% power. This approach can inform choice of processing pipelines and outcome metrics for multisite studies based on scaling factor variability across sites, enabling collaboration between clinical and research institutions.

Identifying preoperative language tracts and predicting postoperative functional recovery using HARDI q-ball fiber tractography in patients with gliomas

OBJECT

Diffusion MRI has uniquely enabled in vivo delineation of white matter tracts, which has been applied to the segmentation of eloquent pathways for intraoperative mapping. The last decade has also seen the development from earlier diffusion tensor models to higher-order models, which take advantage of high angular resolution diffusion-weighted imaging (HARDI) techniques. However, these advanced methods have not been widely implemented for routine preoperative and intraoperative mapping.

The authors report on the application of residual bootstrap q-ball fiber tracking for routine mapping of potentially functional language pathways, the development of a system for rating tract injury to evaluate the impact on clinically assessed language function, and initial results predicting long-term language deficits following glioma resection.

METHODS

The authors have developed methods for the segmentation of 8 putative language pathways including dorsal phonological pathways and ventral semantic streams using residual bootstrap q-ball fiber tracking. Furthermore, they have implemented clinically feasible preoperative acquisition and processing of HARDI data to delineate these pathways for neurosurgical application. They have also developed a rating scale based on the altered fiber tract density to estimate the degree of pathway injury, applying these ratings to a subset of 35 patients with pre- and postoperative fiber tracking. The relationships between specific pathways and clinical language deficits were assessed to determine which pathways are predictive of long-term language deficits following surgery.

RESULTS

This tracking methodology has been routinely implemented for preoperative mapping in patients with brain gliomas who have undergone awake brain tumor resection at the University of California, San Francisco (more than 300 patients to date). In this particular study the authors investigated the white matter structure status and language correlation in a subcohort of 35 subjects both pre- and postsurgery. The rating scales developed for fiber pathway damage were found to be highly reproducible and provided significant correlations with language performance. Preservation of the left arcuate fasciculus (AF) and the temporoparietal component of the superior longitudinal fasciculus (SLF-tp) was consistent in all patients without language deficits (p < 0.001) at the long-term follow-up. Furthermore, in patients with short-term language deficits, the AF and/or SLF-tp were affected, and damage to these 2 pathways was predictive of a long-term language deficit (p = 0.005).

CONCLUSIONS

The authors demonstrated the successful application of q-ball tracking in presurgical planning for language pathways in brain tumor patients and in assessing white matter tract integrity postoperatively to predict long-term language dysfunction. These initial results predicting long-term language deficits following tumor resection indicate that postoperative injury to dorsal language pathways may be prognostic for long-term clinical language deficits.

Study results suggest the importance of dorsal stream tract preservation to reduce language deficits in patients undergoing glioma resection, as well as the potential prognostic value of assessing postoperative injury to dorsal language pathways to predict long-term clinical language deficits.

Neurite Orientation Dispersion and Density Imaging Color Maps to Characterize Brain Diffusion in Neurologic Disorders

PURPOSE

Neurite orientation dispersion and density imaging (NODDI) has recently been developed to overcome diffusion technique limitations in modeling biological systems. This manuscript reports a preliminary investigation into the use of a single color-coded map to represent NODDI-derived information.

MATERIALS AND METHODS

An optimized diffusion-weighted imaging protocol was acquired in several clinical neurological contexts including demyelinating disease, neoplastic process, stroke, and toxic/metabolic disease. The NODDI model was fitted to the diffusion datasets. NODDI is based on a three-compartment diffusion model and provides maps that quantify the contributions to the total diffusion signal in each voxel. The NODDI compartment maps were combined into a single 4-dimensional volume visualized as RGB image (red for anisotropic Gaussian diffusion, green for non-Gaussian anisotropic diffusion, and blue for isotropic Gaussian diffusion), in which the relative contributions of the different microstructural compartments can be easily appreciated.

RESULTS

The NODDI color maps better describe the heterogeneity of neoplastic as well inflammatory lesions by identifying different tissue components within areas apparently homogeneous on conventional imaging. Moreover, NODDI color maps seem to be useful for identifying vasogenic edema differently from tumor-infiltrated edema. In multiple sclerosis, the NODDI color maps enable a visual assessment of the underlying microstructural changes, possibly highlighting an increased inflammatory component, within lesions and potentially in normal-appearing white matter.

CONCLUSION

The NODDI color maps could make this technique valuable in a clinical setting, providing comprehensive and accessible information in normal and pathological brain tissues in different neurological pathologies.

Structural connectivity of the human anterior temporal lobe: A diffusion magnetic resonance imaging study

The anterior temporal lobes (ATL) have been implicated in a range of cognitive functions including auditory and visual perception, language, semantic knowledge, and social-emotional processing. However, the anatomical relationships between the ATLs and the broader cortical networks that subserve these functions have not been fully elucidated. Using diffusion tensor imaging (DTI) and probabilistic tractography, we tested the hypothesis that functional segregation of information in the ATLs is reflected by distinct patterns of structural connectivity to regions outside the ATLs. We performed a parcellation of the ATLs bilaterally based on the degree of connectivity of each voxel with eight ipsilateral target regions known to be involved in various cognitive networks. Six discrete segments within each ATL showed preferential connectivity to one of the ipsilateral target regions, via four major fiber tracts (uncinate, inferior longitudinal, middle longitudinal, and arcuate fasciculi). Two noteworthy interhemispheric differences were observed: connections between the ATL and orbito-frontal areas were stronger in the right hemisphere, while the consistency of the connection between the ATL and the inferior frontal gyrus through the arcuate fasciculus was greater in the left hemisphere. Our findings support the hypothesis that distinct regions within the ATLs have anatomical connections to different cognitive networks. Hum Brain Mapp 37:2210-2222, 2016. © 2016 Wiley Periodicals, Inc.

Secure attachment status is associated with white matter integrity in healthy young adults

The present study investigates associations between security of attachment in the mother-child relationship and patterns of brain connectivity in young adults. We hypothesized that secure attachment would relate to more efficient connectivity in white matter association fibers due to increased myelination. Attachment security was measured in 53 young adults using the Kerns Security Scale; anatomical information was acquired using diffusion tensor imaging. Higher fractional anisotropy, an index of directionality of diffusion, related to security of attachment in four left-hemisphere white matter association fibers (uncinate fasciculus, cingulum, superior longitudinal fasciculus, and inferior fronto-occipital fasciculus). As expected, this result was mainly ascribable to increased myelination, which has been independently associated with attachment security. Security of attachment may have an identifiable biological basis. Our research demonstrates the feasibility of coupling neuroimaging tools with clinical investigation.

Association Between Thoracic Spinal Cord Gray Matter Atrophy and Disability in Multiple Sclerosis

IMPORTANCE

In multiple sclerosis (MS), upper cervical cord gray matter (GM) atrophy correlates more strongly with disability than does brain or cord white matter (WM) atrophy. The corresponding relationships in the thoracic cord are unknown owing to technical difficulties in assessing GM and WM compartments by conventional magnetic resonance imaging techniques.

OBJECTIVES

To investigate the associations between MS disability and disease type with lower thoracic cord GM and WM areas using phase-sensitive inversion recovery magnetic resonance imaging at 3 T, as well as to compare these relationships with those obtained at upper cervical levels.

DESIGN, SETTING, AND PARTICIPANTS

Between July 2013 and March 2014, a total of 142 patients with MS (aged 25-75 years; 86 women) and 20 healthy control individuals were included in this cross-sectional observational study conducted at an academic university hospital.

MAIN OUTCOMES AND MEASURES

Total cord areas (TCAs), GM areas, and WM areas at the disc levels C2/C3, C3/C4, T8/9, and T9/10. Area differences between groups were assessed, with age and sex as covariates.

RESULTS

Patients with relapsing MS (RMS) had smaller thoracic cord GM areas than did age- and sex-matched control individuals (mean differences [coefficient of variation (COV)]: 0.98 mm2 [9.2%]; P = .003 at T8/T9 and 0.93 mm2 [8.0%]; P = .01 at T9/T10); however, there were no significant differences in either the WM area or TCA. Patients with progressive MS showed smaller GM areas (mean differences [COV]: 1.02 mm2 [10.6%]; P < .001 at T8/T9 and 1.37 mm2 [13.2%]; P < .001 at T9/T10) and TCAs (mean differences [COV]: 3.66 mm2 [9.0%]; P < .001 at T8/T9 and 3.04 mm2 [7.2%]; P = .004 at T9/T10) compared with patients with RMS. All measurements (GM, WM, and TCA) were inversely correlated with Expanded Disability Status Scale score. Thoracic cord GM areas were correlated with lower limb function. In multivariable models (which also included cord WM areas and T2 lesion number, brain WM volumes, brain T1 and fluid-attenuated inversion recovery lesion loads, age, sex, and disease duration), cervical cord GM areas had the strongest correlation with Expanded Disability Status Scale score followed by thoracic cord GM area and brain GM volume.

CONCLUSIONS AND RELEVANCE

Thoracic cord GM atrophy can be detected in vivo in the absence of WM atrophy in RMS. This atrophy is more pronounced in progressive MS than RMS and correlates with disability and lower limb function. Our results indicate that remarkable cord GM atrophy is present at multiple cervical and lower thoracic levels and, therefore, may reflect widespread cord GM degeneration.

2D phase-sensitive inversion recovery imaging to measure in vivo spinal cord gray and white matter areas in clinically feasible acquisition times

PURPOSE

To present and assess a procedure for measurement of spinal cord total cross-sectional areas (TCA) and gray matter (GM) areas based on phase-sensitive inversion recovery imaging (PSIR). In vivo assessment of spinal cord GM and white matter (WM) could become pivotal to study various neurological diseases, but it is challenging because of insufficient GM/WM contrast provided by conventional magnetic resonance imaging (MRI).

MATERIALS AND METHODS

We acquired 2D PSIR images at 3T at each disc level of the spinal axis in 10 healthy subjects and measured TCA, cord diameters, WM and GM areas, and GM area/TCA ratios. Second, we investigated 32 healthy subjects at four selected levels (C2-C3, C3-C4, T8-T9, T9-T10, total acquisition time <8 min) and generated normative reference values of TCA and GM areas. We assessed test-retest, intra- and interoperator reliability of the acquisition strategy, and measurement steps.

RESULTS

The measurement procedure based on 2D PSIR imaging allowed TCA and GM area assessments along the entire spinal cord axis. The tests we performed revealed high test-retest/intraoperator reliability (mean coefficient of variation [COV] at C2-C3: TCA = 0.41%, GM area = 2.75%) and interoperator reliability of the measurements (mean COV on the 4 levels: TCA = 0.44%, GM area = 4.20%; mean intraclass correlation coefficient: TCA = 0.998, GM area = 0.906).

CONCLUSION

2D PSIR allows reliable in vivo assessment of spinal cord TCA, GM, and WM areas in clinically feasible acquisition times. The area measurements presented here are in agreement with previous MRI and postmortem studies.

White matter involvement in sporadic Creutzfeldt-Jakob disease

Sporadic Creutzfeldt-Jakob disease is considered primarily a disease of grey matter, although the extent of white matter involvement has not been well described. We used diffusion tensor imaging to study the white matter in sporadic Creutzfeldt-Jakob disease compared to healthy control subjects and to correlated magnetic resonance imaging findings with histopathology. Twenty-six patients with sporadic Creutzfeldt-Jakob disease and nine age- and gender-matched healthy control subjects underwent volumetric T₁-weighted and diffusion tensor imaging. Six patients had post-mortem brain analysis available for assessment of neuropathological findings associated with prion disease. Parcellation of the subcortical white matter was performed on 3D T₁-weighted volumes using Freesurfer. Diffusion tensor imaging maps were calculated and transformed to the 3D-T₁ space; the average value for each diffusion metric was calculated in the total white matter and in regional volumes of interest. Tract-based spatial statistics analysis was also performed to investigate the deeper white matter tracts. There was a significant reduction of mean (P = 0.002), axial (P = 0.0003) and radial (P = 0.0134) diffusivities in the total white matter in sporadic Creutzfeldt-Jakob disease. Mean diffusivity was significantly lower in most white matter volumes of interest (P < 0.05, corrected for multiple comparisons), with a generally symmetric pattern of involvement in sporadic Creutzfeldt-Jakob disease. Mean diffusivity reduction reflected concomitant decrease of both axial and radial diffusivity, without appreciable changes in white matter anisotropy. Tract-based spatial statistics analysis showed significant reductions of mean diffusivity within the white matter of patients with sporadic Creutzfeldt-Jakob disease, mainly in the left hemisphere, with a strong trend (P = 0.06) towards reduced mean diffusivity in most of the white matter bilaterally. In contrast, by visual assessment there was no white matter abnormality either on T₂-weighted or diffusion-weighted images. Widespread reduction in white matter mean diffusivity, however, was apparent visibly on the quantitative attenuation coefficient maps compared to healthy control subjects. Neuropathological analysis showed diffuse astrocytic gliosis and activated microglia in the white matter, rare prion deposition and subtle subcortical microvacuolization, and patchy foci of demyelination with no evident white matter axonal degeneration. Decreased mean diffusivity on attenuation coefficient maps might be associated with astrocytic gliosis. We show for the first time significant global reduced mean diffusivity within the white matter in sporadic Creutzfeldt-Jakob disease, suggesting possible primary involvement of the white matter, rather than changes secondary to neuronal degeneration/loss.

Sporadic Creutzfeldt-Jakob disease (sCJD) is considered primarily a disease of grey matter. However, Caverzasi et al. now show a global decrease in mean diffusivity in white matter. The changes appear to be associated with reactive astrocytic gliosis and activated microglia, and suggest primary involvement of the white matter in sCJD.

Spinal cord gray matter atrophy correlates with multiple sclerosis disability

OBJECTIVE

In multiple sclerosis (MS), cerebral gray matter (GM) atrophy correlates more strongly than white matter (WM) atrophy with disability. The corresponding relationships in the spinal cord (SC) are unknown due to technical limitations in assessing SC GM atrophy. Using phase-sensitive inversion recovery (PSIR) magnetic resonance imaging, we determined the association of the SC GM and SC WM areas with MS disability and disease type.

METHODS

A total of 113 MS patients and 20 healthy controls were examined at 3T with a PSIR sequence acquired at the C2/C3 disk level. Two independent, clinically masked readers measured the cord WM and GM areas. Correlations between cord areas and Expanded Disability Status Score (EDSS) were determined. Differences in areas between groups were assessed with age and sex as covariates.

RESULTS

Relapsing MS (RMS) patients showed smaller SC GM areas than age- and sex-matched controls (p = 0.008) without significant differences in SC WM areas. Progressive MS patients showed smaller SC GM and SC WM areas compared to RMS patients (all p ≤ 0.004). SC GM, SC WM, and whole cord areas inversely correlated with EDSS (rho: -0.60, -0.32, -0.42, respectively; all p ≤ 0.001). The SC GM area was the strongest correlate of disability in multivariate models including brain GM and WM volumes, fluid-attenuated inversion recovery lesion load, T1 lesion load, SC WM area, number of SC T2 lesions, age, sex, and disease duration. Brain and spinal GM independently contributed to EDSS.

INTERPRETATION

SC GM atrophy is detectable in vivo in the absence of WM atrophy in RMS. It is more pronounced in progressive MS than RMS and contributes more to patient disability than SC WM or brain GM atrophy.

Q-ball of inferior fronto-occipital fasciculus and beyond

The inferior fronto-occipital fasciculus (IFOF) is historically described as the longest associative bundle in the human brain and it connects various parts of the occipital cortex, temporo-basal area and the superior parietal lobule to the frontal lobe through the external/extreme capsule complex. The exact functional role and the detailed anatomical definition of the IFOF are still under debate within the scientific community. In this study we present a fiber tracking dissection of the right and left IFOF by using a q-ball residual-bootstrap reconstruction of High-Angular Resolution Diffusion Imaging (HARDI) data sets in 20 healthy subjects. By defining a single seed region of interest on the coronal fractional anisotropy (FA) color map of each subject, we investigated all the pathways connecting the parietal, occipital and posterior temporal cortices to the frontal lobe through the external/extreme capsule. In line with recent post-mortem dissection studies we found more extended anterior-posterior association connections than the “classical” fronto-occipital representation of the IFOF. In particular the pathways we evidenced showed: a) diffuse projections in the frontal lobe, b) fronto-parietal lobes connections trough the external capsule in almost all the subjects and c) widespread connections in the posterior regions. Our study represents the first consistent in vivo demonstration across a large group of individuals of these novel anterior and posterior terminations of the IFOF detailed described only by post-mortem anatomical dissection. Furthermore our work establishes the feasibility of consistent in vivo mapping of this architecture with independent in vivo methodologies. In conclusion q-ball tractography dissection supports a more complex definition of IFOF, which includes several subcomponents likely underlying specific function.

Interpersonal competence in young adulthood and right laterality in white matter

The right hemisphere of the human brain is known to be involved in processes underlying emotion and social cognition. Clinical neuropsychology investigations and brain lesion studies have linked a number of personality and social disorders to abnormal white matter (WM) integrity in the right hemisphere. Here, we tested the hypothesis that interpersonal competencies are associated with integrity of WM tracts in the right hemisphere of healthy young adults. Thirty-one participants underwent diffusion tensor imaging scanning. Fractional anisotropy was used to quantify water diffusion. After the scanning session, participants completed the Adolescent Interpersonal Competence Questionnaire. Fractional anisotropy was subsequently correlated with Adolescent Interpersonal Competence Questionnaire scores using tract-based spatial statistics. Higher interpersonal competencies are related to higher WM integrity in several major tracts of the right hemisphere, in specific the uncinate fasciculus, the cingulum, the forceps minor, the infero-fronto occipital fasciculus, the inferior longitudinal fasciculus, and the superior longitudinal fasciculus. These results provide the first direct analysis of the neuroanatomical basis of interpersonal competencies and young adult self-reported skills in social contexts.

Two-step magnetic ordering in quasi-one-dimensional helimagnets: possible experimental validation of villain's conjecture about a chiral spin liquid phase

Low-temperature specific heat, magnetic susceptibility, and zero-field muon spin resonance (microSR) measurements have been performed in the quasi-one-dimensional molecular helimagnetic compound Gd(hfac)3NITEt. The specific heat presents two anomalies at T(0)=2.19+/-0.02 K and T(N)=1.88+/-0.02 K, which both disappear upon the application of a weak magnetic field. Conversely, magnetic susceptibility and muSR data show the divergence of two-spin correlation functions only at T(N)=1.88+/-0.02 K. These results suggest an experimental validation of Villain's conjecture of a two-step magnetic ordering in quasi-one-dimensional XY helimagnets; i.e., the paramagnetic phase and the helical spin solid phase are separated by a chiral spin liquid phase, where translational invariance is broken without violation of rotational invariance.

Very-low-frequency excitations in frustrated two-dimensional S = 1/2 Heisenberg antiferromagnets

Muon spin resonance and 7Li NMR relaxation measurements in frustrated two-dimensional S = 1/2 Heisenberg antiferromagnets on a square lattice are presented. It is found that, in both Li2VOSiO4 and Li2VOGeO4, spin dynamics at frequencies orders of magnitude below the Heisenberg exchange frequency are present. These dynamics are associated with the motions of walls separating coexisting collinear domains with a magnetic wave vector rotated by 90 degrees.