Tractostorm: The what, why, and how of tractography dissection reproducibility.

Investigative studies of white matter (WM) brain structures using diffusion MRI (dMRI) tractography frequently require manual WM bundle segmentation, often called "virtual dissection." Human errors and personal decisions make these manual segmentations hard to reproduce, which have not yet been quantified by the dMRI community. It is our opinion that if the field of dMRI tractography wants to be taken seriously as a widespread clinical tool, it is imperative to harmonize WM bundle segmentations and develop protocols aimed to be used in clinical settings. The EADC-ADNI Harmonized Hippocampal Protocol achieved such standardization through a series of steps that must be reproduced for every WM bundle. This article is an observation of the problematic. A specific bundle segmentation protocol was used in order to provide a real-life example, but the contribution of this article is to discuss the need for reproducibility and standardized protocol, as for any measurement tool. This study required the participation of 11 experts and 13 nonexperts in neuroanatomy and "virtual dissection" across various laboratories and hospitals. Intra-rater agreement (Dice score) was approximately 0.77, while inter-rater was approximately 0.65. The protocol provided to participants was not necessarily optimal, but its design mimics, in essence, what will be required in future protocols. Reporting tractometry results such as average fractional anisotropy, volume or streamline count of a particular bundle without a sufficient reproducibility score could make the analysis and interpretations more difficult. Coordinated efforts by the diffusion MRI tractography community are needed to quantify and account for reproducibility of WM bundle extraction protocols in this era of open and collaborative science.

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.

Evaluating syntactic comprehension during awake intraoperative cortical stimulation mapping.

OBJECTIVE: Electrocortical stimulation mapping (ECS) is widely used to identify essential language areas, but sentence-level processing has rarely been investigated. METHODS: While undergoing awake surgery in the dominant left hemisphere, 6 subjects were asked to comprehend sentences varying in their demands on syntactic processing. RESULTS: In all 6 subjects, stimulation of the inferior frontal gyrus disrupted comprehension of passive sentences, which critically depend on syntactic processing to correctly assign grammatical roles, without disrupting comprehension of simpler tasks. In 4 of the 6 subjects, these sites were localized to the pars opercularis. Sentence comprehension was also disrupted by stimulation of other perisylvian sites, but in a more variable manner. CONCLUSIONS: These findings suggest that there may be language regions that differentially contribute to sentence processing and which therefore are best identified using sentence-level tasks. The functional consequences of resecting these sites remain to be investigated.

Cortically constrained shape recognition: Automated white matter tract segmentation validated in the pediatric brain.

BACKGROUND AND PURPOSE: Manual segmentation of white matter (WM) bundles requires extensive training and is prohibitively labor-intensive for large-scale studies. Automated segmentation methods are necessary in order to eliminate operator dependency and to enable reproducible studies. Significant changes in the WM landscape throughout childhood require flexible methods to capture the variance across the span of brain development. METHODS: Here, we describe a novel automated segmentation tool called Cortically Constrained Shape Recognition (CCSR), which combines two complementary approaches: (1) anatomical connectivity priors based on FreeSurfer-derived regions of interest and (2) shape priors based on 3-dimensional streamline bundle atlases applied using RecoBundles. We tested the performance and repeatability of this approach by comparing volume and diffusion metrics of the main language WM tracts that were both manually and automatically segmented in a pediatric cohort acquired at the UCSF Dyslexia Center (n = 59; 25 females; average age: 11 ± 2; range: 7-14). RESULTS: The CCSR approach showed high agreement with the expert manual segmentations: across all tracts, the spatial overlap between tract volumes showed an average Dice Similarity Coefficient (DSC) of 0.76, and the fractional anisotropy (FA) on average had a Lin's Concordance Correlation Coefficient (CCC) of 0.81. The CCSR's repeatability in a subset of this cohort achieved a DSC of 0.92 on average across all tracts. CONCLUSION: This novel automated segmentation approach is a promising tool for reproducible large-scale tractography analyses in pediatric populations and particularly for the quantitative assessment of structural connections underlying various clinical presentations in neurodevelopmental disorders.

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.

fMRI-Targeted High-Angular Resolution Diffusion MR Tractography to Identify Functional Language Tracts in Healthy Controls and Glioma Patients.

BACKGROUND: MR Tractography enables non-invasive preoperative depiction of language subcortical tracts, which is crucial for the presurgical work-up of brain tumors; however, it cannot evaluate the exact function of the fibers. PURPOSE: A systematic pipeline was developed to combine tractography reconstruction of language fiber bundles, based on anatomical landmarks (Anatomical-T), with language fMRI cortical activations. A fMRI-targeted Tractography (fMRI-T) was thus obtained, depicting the subsets of the anatomical tracts whose endpoints are located inside a fMRI activation. We hypothesized that fMRI-T could provide additional functional information regarding the subcortical structures, better reflecting the eloquent white matter structures identified intraoperatively. METHODS: Both Anatomical-T and fMRI-T of language fiber tracts were performed on 16 controls and preoperatively on 16 patients with left-hemisphere brain tumors, using a q-ball residual bootstrap algorithm based on High Angular Resolution Diffusion Imaging (HARDI) datasets (b = 3000 s/mm2; 60 directions); fMRI ROIs were obtained using picture naming, verbal fluency, and auditory verb generation tasks. In healthy controls, normalized MNI atlases of fMRI-T and Anatomical-T were obtained. In patients, the surgical resection of the tumor was pursued by identifying eloquent structures with intraoperative direct electrical stimulation mapping and extending surgery to the functional boundaries. Post-surgical MRI allowed to identify Anatomical-T and fMRI-T non-eloquent portions removed during the procedure. RESULTS: MNI Atlases showed that fMRI-T is a subset of Anatomical-T, and that different task-specific fMRI-T involve both shared subsets and task-specific subsets - e.g., verbal fluency fMRI-T strongly involves dorsal frontal tracts, consistently with the phonogical-articulatory features of this task. A quantitative analysis in patients revealed that Anatomical-T removed portions of AF-SLF and IFOF were significantly greater than verbal fluency fMRI-T ones, suggesting that fMRI-T is a more specific approach. In addition, qualitative analyses showed that fMRI-T AF-SLF and IFOF predict the exact functional limits of resection with increased specificity when compared to Anatomical-T counterparts, especially the superior frontal portion of IFOF, in a subcohort of patients. CONCLUSION: These results suggest that performing fMRI-T in addition to the 'classic' Anatomical-T may be useful in a preoperative setting to identify the 'high-risk subsets' that should be spared during the surgical procedure.

Cluster Confidence Index: A Streamline-Wise Pathway Reproducibility Metric for Diffusion-Weighted MRI Tractography.

BACKGROUND: Diffusion-weighted magnetic resonance imaging tractography can be used to create models of white matter fascicles. Anatomical and pathological variability between subjects can drastically alter the tractography output, so standardizing results across a cohort is nontrivial. Furthermore, tractography methods have inherently low reproducibility due to stochasticity (for probabilistic methods) and subjective decisions, since the final fascicle model often requires a manual intervention step performed by an expert human operator to control both outliers and systematic false-positive pathways, as defined by prior knowledge of anatomy. METHODS: We present an approach that computationally assigns a cluster confidence index (CCI) reflecting the reproducibility of that pathway in the context of a streamline dataset. This metric is a tractography algorithm-agnostic tool that can be applied to any dataset of streamlines. RESULTS: Applications of this metric include systematic elimination of outlier streamlines using a CCI threshold and interactive filtering by CCI to facilitate manual segmentation of fascicle models. CONCLUSIONS: This method is intended to replace the application of a streamline density threshold so that outliers are eliminated based on low pathway density instead of voxel-wise density.

Subcortical stimulation mapping of descending motor pathways for perirolandic gliomas: assessment of morbidity and functional outcome in 702 cases.

OBJECTIVE: Herein, the authors report their experience with intraoperative stimulation mapping to locate the descending subcortical motor pathways in patients undergoing surgery for hemispheric gliomas within or adjacent to the rolandic cortex, with particular description of the morbidity and functional outcomes associated with this technique. METHODS: This is a retrospective analysis of patients who, in the period between 1997 and 2016, had undergone resection of hemispheric perirolandic gliomas within or adjacent to descending motor pathways. Data regarding intraoperative stimulation mapping and patient postoperative neurological status were collected. RESULTS: Of 702 patients, stimulation mapping identified the descending motor pathways in 300 cases (43%). A new or worsened motor deficit was seen postoperatively in 210 cases (30%). Among these 210 cases, there was improvement in motor function to baseline levels by 3 months postoperatively in 161 cases (77%), whereas the deficit remained in 49 cases (23%). The majority (65%) of long-term deficits (persisting beyond 3 months) were mild or moderate (antigravity strength or better). On multivariate analysis, patients in whom the subcortical motor pathways had been identified with stimulation mapping during surgery were more likely to develop an additional and/or worsened motor deficit postoperatively than were those in whom the subcortical pathways had not been found (45% vs 19%, respectively, p < 0.001). This difference remained when considering the likelihood of a long-term deficit (i.e., persisting > 3 months; 12% vs 3.2%, p < 0.001). A higher tumor grade and the presence of a preoperative motor deficit were also associated with higher rates of motor deficits persisting long-term. A region of restricted diffusion adjacent to the resection cavity was seen in 20 patients with long-term deficits (41%) and was more common in cases in which the motor pathways were not identified (69%). Long-term deficits that occur in settings in which the subcortical motor pathways are not identified seem in large part due to ischemic injury to descending tracts. CONCLUSIONS: Stimulation mapping allows surgeons to identify the descending motor pathways during resection of tumors in perirolandic regions and to attain an acceptable rate of morbidity in these high-risk cases.

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.

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.