Pandora: 4-D White Matter Bundle Population-Based Atlases Derived from Diffusion MRI Fiber Tractography

Association of concussion history with psychiatric symptoms, limbic system structure, and kynurenine pathway metabolites in healthy, collegiate-aged athletes

Psychiatric outcomes are commonly observed in individuals with repeated concussions, though their underlying mechanism is unknown. One potential mechanism linking concussion with psychiatric symptoms is inflammation-induced activation of the kynurenine pathway, which is thought to play a role in the pathogenesis of mood disorders. Here, we investigated the association of prior concussion with multiple psychiatric-related outcomes in otherwise healthy male and female collegiate-aged athletes (N = 212) with varying histories of concussion recruited from the community. Specially, we tested the hypotheses that concussion history is associated with worse psychiatric symptoms, limbic system structural abnormalities (hippocampal volume, white matter microstructure assessed using neurite orientation dispersion and density imaging; NODDI), and elevations in kynurenine pathway (KP) metabolites (e.g., Quinolinic acid; QuinA). Given known sex-effects on concussion risk and recovery, psychiatric outcomes, and the kynurenine pathway, the moderating effect of sex was considered for all analyses. More concussions were associated with greater depression, anxiety, and anhedonia symptoms in female athletes (ps ≤ 0.005) and greater depression symptoms in male athletes (p = 0.011). More concussions were associated with smaller bilateral hippocampal tail (ps < 0.010) and left hippocampal body (p < 0.001) volumes across male and female athletes. Prior concussion was also associated with elevations in the orientation dispersion index (ODI) and lower intracellular volume fraction in several white matter tracts including the in uncinate fasciculus, cingulum-gyrus, and forceps major and minor, with evidence of female-specific associations in select regions. Regarding serum KP metabolites, more concussions were associated with elevated QuinA in females and lower tryptophan in males (ps ≤ 0.010). Finally, serum levels of QuinA were associated with elevated ODI (male and female athletes) and worse anxiety symptoms (females only), while higher ODI in female athletes and smaller hippocampal volumes in male athletes were associated with more severe anxiety and depression symptoms (ps ≤ 0.05). These data suggest that cumulative concussion is associated with psychiatric symptoms and limbic system structure in healthy athletes, with increased susceptibility to these effects in female athletes. Moreover, the associations of outcomes with serum KP metabolites highlight the KP as one potential molecular pathway underlying these observations.

A review of diffusion MRI in mood disorders: mechanisms and predictors of treatment response

By measuring the molecular diffusion of water molecules in brain tissue, diffusion MRI (dMRI) provides unique insight into the microstructure and structural connections of the brain in living subjects. Since its inception, the application of dMRI in clinical research has expanded our understanding of the possible biological bases of psychiatric disorders and successful responses to different therapeutic interventions. Here, we review the past decade of diffusion imaging-based investigations with a specific focus on studies examining the mechanisms and predictors of therapeutic response in people with mood disorders. We present a brief overview of the general application of dMRI and key methodological developments in the field that afford increasingly detailed information concerning the macro- and micro-structural properties and connectivity patterns of white matter (WM) pathways and their perturbation over time in patients followed prospectively while undergoing treatment. This is followed by a more in-depth summary of particular studies using dMRI approaches to examine mechanisms and predictors of clinical outcomes in patients with unipolar or bipolar depression receiving pharmacological, neurostimulation, or behavioral treatments. Limitations associated with dMRI research in general and with treatment studies in mood disorders specifically are discussed, as are directions for future research. Despite limitations and the associated discrepancies in findings across individual studies, evolving research strongly indicates that the field is on the precipice of identifying and validating dMRI biomarkers that could lead to more successful personalized treatment approaches and could serve as targets for evaluating the neural effects of novel treatments.

Characterizing patterns of diffusion tensor imaging variance in aging brains

Purpose

As large analyses merge data across sites, a deeper understanding of variance in statistical assessment across the sources of data becomes critical for valid analyses. Diffusion tensor imaging (DTI) exhibits spatially varying and correlated noise, so care must be taken with distributional assumptions. Here we characterize the role of physiology, subject compliance, and the interaction of subject with the scanner in the understanding of DTI variability, as modeled in spatial variance of derived metrics in homogeneous regions.

Approach

We analyze DTI data from 1035 subjects in the Baltimore Longitudinal Study of Aging (BLSA), with ages ranging from 22.4 to 103 years old. For each subject, up to 12 longitudinal sessions were conducted. We assess variance of DTI scalars within regions of interest (ROIs) defined by four segmentation methods and investigate the relationships between the variance and covariates, including baseline age, time from the baseline (referred to as "interval"), motion, sex, and whether it is the first scan or the second scan in the session.

Results

Covariate effects are heterogeneous and bilaterally symmetric across ROIs. Inter-session interval is positively related ( p ≪ 0.001 ) to FA variance in the cuneus and occipital gyrus, but negatively ( p ≪ 0.001 ) in the caudate nucleus. Males show significantly ( p ≪ 0.001 ) higher FA variance in the right putamen, thalamus, body of the corpus callosum, and cingulate gyrus. In 62 out of 176 ROIs defined by the Eve type-1 atlas, an increase in motion is associated ( p < 0.05 ) with a decrease in FA variance. Head motion increases during the rescan of DTI ( Δ μ = 0.045 millimeters per volume).

Conclusions

The effects of each covariate on DTI variance, and their relationships across ROIs are complex. Ultimately, we encourage researchers to include estimates of variance when sharing data and consider models of heteroscedasticity in analysis. This work provides a foundation for study planning to account for regional variations in metric variance.

Characterizing patterns of diffusion tensor imaging variance in aging brains

Purpose

As large analyses merge data across sites, a deeper understanding of variance in statistical assessment across the sources of data becomes critical for valid analyses. Diffusion tensor imaging (DTI) exhibits spatially varying and correlated noise, so care must be taken with distributional assumptions. Here, we characterize the role of physiology, subject compliance, and the interaction of the subject with the scanner in the understanding of DTI variability, as modeled in the spatial variance of derived metrics in homogeneous regions.

Approach

We analyze DTI data from 1035 subjects in the Baltimore Longitudinal Study of Aging, with ages ranging from 22.4 to 103 years old. For each subject, up to 12 longitudinal sessions were conducted. We assess the variance of DTI scalars within regions of interest (ROIs) defined by four segmentation methods and investigate the relationships between the variance and covariates, including baseline age, time from the baseline (referred to as "interval"), motion, sex, and whether it is the first scan or the second scan in the session.

Results

Covariate effects are heterogeneous and bilaterally symmetric across ROIs. Inter-session interval is positively related ( p ≪ 0.001 ) to FA variance in the cuneus and occipital gyrus, but negatively ( p ≪ 0.001 ) in the caudate nucleus. Males show significantly ( p ≪ 0.001 ) higher FA variance in the right putamen, thalamus, body of the corpus callosum, and cingulate gyrus. In 62 out of 176 ROIs defined by the Eve type-1 atlas, an increase in motion is associated ( p < 0.05 ) with a decrease in FA variance. Head motion increases during the rescan of DTI ( Δ μ = 0.045 mm per volume).

Conclusions

The effects of each covariate on DTI variance and their relationships across ROIs are complex. Ultimately, we encourage researchers to include estimates of variance when sharing data and consider models of heteroscedasticity in analysis. This work provides a foundation for study planning to account for regional variations in metric variance.

White matter alterations and their associations with biomarkers and behavior in subjective cognitive decline individuals: a fixel-based analysis

BACKGROUND

Subjective cognitive decline (SCD) is an early stage of dementia linked to Alzheimer's disease pathology. White matter changes were found in SCD using diffusion tensor imaging, but there are known limitations in voxel-wise tensor-based methods. Fixel-based analysis (FBA) can help understand changes in white matter fibers and how they relate to neurodegenerative proteins and multidomain behavior data in individuals with SCD.

METHODS

Healthy adults with normal cognition were recruited in the Northeastern Taiwan Community Medicine Research Cohort in 2018-2022 and divided into SCD and normal control (NC). Participants underwent evaluations to assess cognitive abilities, mental states, physical activity levels, and susceptibility to fatigue. Neurodegenerative proteins were measured using an immunomagnetic reduction technique. Multi-shell diffusion MRI data were collected and analyzed using whole-brain FBA, comparing results between groups and correlating them with multidomain assessments.

RESULTS

The final enrollment included 33 SCD and 46 NC participants, with no significant differences in age, sex, or education between the groups. SCD had a greater fiber-bundle cross-section than NC (pFWE < 0.05) at bilateral frontal superior longitudinal fasciculus II (SLFII). These white matter changes correlate negatively with plasma Aβ42 level (r = -0.38, p = 0.01) and positively with the AD8 score for subjective cognitive complaints (r = 0.42, p = 0.004) and the Hamilton Anxiety Rating Scale score for the degree of anxiety (Ham-A, r = 0.35, p = 0.019). The dimensional analysis of FBA metrics and blood biomarkers found positive correlations of plasma neurofilament light chain with fiber density at the splenium of corpus callosum (pFWE < 0.05) and with fiber-bundle cross-section at the right thalamus (pFWE < 0.05). Further examination of how SCD grouping interacts between the correlations of FBA metrics and multidomain assessments showed interactions between the fiber density at the corpus callosum with letter-number sequencing cognitive score (pFWE < 0.01) and with fatigue to leisure activities (pFWE < 0.05).

CONCLUSION

Based on FBA, our investigation suggests white matter structural alterations in SCD. The enlargement of SLFII's fiber cross-section is linked to plasma Aβ42 and neuropsychiatric symptoms, which suggests potential early axonal dystrophy associated with Alzheimer's pathology in SCD. The splenium of the corpus callosum is also a critical region of axonal degeneration and cognitive alteration for SCD.

Enhancing Hierarchical Transformers for Whole Brain Segmentation with Intracranial Measurements Integration

Whole brain segmentation with magnetic resonance imaging (MRI) enables the non-invasive measurement of brain regions, including total intracranial volume (TICV) and posterior fossa volume (PFV). Enhancing the existing whole brain segmentation methodology to incorporate intracranial measurements offers a heightened level of comprehensiveness in the analysis of brain structures. Despite its potential, the task of generalizing deep learning techniques for intracranial measurements faces data availability constraints due to limited manually annotated atlases encompassing whole brain and TICV/PFV labels. In this paper, we enhancing the hierarchical transformer UNesT for whole brain segmentation to achieve segmenting whole brain with 133 classes and TICV/PFV simultaneously. To address the problem of data scarcity, the model is first pretrained on 4859 T1-weighted (T1w) 3D volumes sourced from 8 different sites. These volumes are processed through a multi-atlas segmentation pipeline for label generation, while TICV/PFV labels are unavailable. Subsequently, the model is finetuned with 45 T1w 3D volumes from Open Access Series Imaging Studies (OASIS) where both 133 whole brain classes and TICV/PFV labels are available. We evaluate our method with Dice similarity coefficients(DSC). We show that our model is able to conduct precise TICV/PFV estimation while maintaining the 132 brain regions performance at a comparable level. Code and trained model are available at: https://github.com/MASILab/UNesT/wholebrainSeg.

The missing third dimension-Functional correlations of BOLD signals incorporating white matter

Correlations between magnetic resonance imaging (MRI) blood oxygenation level-dependent (BOLD) signals from pairs of gray matter areas are used to infer their functional connectivity, but they are unable to describe how white matter is engaged in brain networks. Recently, evidence that BOLD signals in white matter are robustly detectable and are modulated by neural activities has accumulated. We introduce a three-way correlation between BOLD signals from pairs of gray matter volumes (nodes) and white matter bundles (edges) to define the communication connectivity through each white matter bundle. Using MRI images from publicly available databases, we show, for example, that the three-way connectivity is influenced by age. By integrating functional MRI signals from white matter as a third component in network analyses, more comprehensive descriptions of brain function may be obtained.

Associations of prior concussion severity with brain microstructure using mean apparent propagator magnetic resonance imaging

Magnetic resonance imaging (MRI) diffusion studies have shown chronic microstructural tissue abnormalities in athletes with history of concussion, but with inconsistent findings. Concussions with post-traumatic amnesia (PTA) and/or loss of consciousness (LOC) have been connected to greater physiological injury. The novel mean apparent propagator (MAP) MRI is expected to be more sensitive to such tissue injury than the conventional diffusion tensor imaging. This study examined effects of prior concussion severity on microstructure with MAP-MRI. Collegiate-aged athletes (N = 111, 38 females; ≥6 months since most recent concussion, if present) completed semistructured interviews to determine the presence of prior concussion and associated injury characteristics, including PTA and LOC. MAP-MRI metrics (mean non-Gaussian diffusion [NG Mean], return-to-origin probability [RTOP], and mean square displacement [MSD]) were calculated from multi-shell diffusion data, then evaluated for associations with concussion severity through group comparisons in a primary model (athletes with/without prior concussion) and two secondary models (athletes with/without prior concussion with PTA and/or LOC, and athletes with/without prior concussion with LOC only). Bayesian multilevel modeling estimated models in regions of interest (ROI) in white matter and subcortical gray matter, separately. In gray matter, the primary model showed decreased NG Mean and RTOP in the bilateral pallidum and decreased NG Mean in the left putamen with prior concussion. In white matter, lower NG Mean with prior concussion was present in all ROI across all models and was further decreased with LOC. However, only prior concussion with LOC was associated with decreased RTOP and increased MSD across ROI. Exploratory analyses conducted separately in male and female athletes indicate associations in the primary model may differ by sex. Results suggest microstructural measures in gray matter are associated with a general history of concussion, while a severity-dependent association of prior concussion may exist in white matter.

One-shot segmentation of novel white matter tracts via extensive data augmentation and adaptive knowledge transfer

The use of convolutional neural networks (CNNs) has allowed accurate white matter (WM) tract segmentation on diffusion magnetic resonance imaging (dMRI). To train the CNN-based segmentation models, a large number of scans on which WM tracts are annotated need to be collected, and these annotated scans can be accumulated over a long period of time. However, when novel WM tracts that are different from existing annotated WM tracts are of interest, additional annotations are required for their segmentation. Due to the cost of manual annotations, methods have been developed for few-shot segmentation of novel WM tracts, where the segmentation knowledge is transferred from existing WM tracts to novel WM tracts and the amount of annotated data for novel WM tracts is reduced. Despite these developments, it is desirable to further reduce the amount of annotated data to the one-shot setting with a single annotated image. To address this problem, we develop an approach to one-shot segmentation of novel WM tracts. Our method follows the existing pretraining/fine-tuning framework that transfers segmentation knowledge from existing to novel WM tracts. First, as there is extremely scarce annotated data in the one-shot setting, we design several different data augmentation strategies so that extensive data augmentation can be performed to obtain extra synthetic training data. The data augmentation strategies are based on image masking and thus applicable to the one-shot setting. Second, to address overfitting and knowledge forgetting in the fine-tuning stage that can be more severe given limited training data, we propose an adaptive knowledge transfer strategy that selects the network weights to be updated. The data augmentation and adaptive knowledge transfer strategies are combined to train the segmentation model. Considering that the different data augmentation strategies can generate synthetic data that contain potentially conflicting information, we apply the data augmentation strategies separately, each leading to a different segmentation model. The results predicted by the different models are fused to produce the final segmentation. We validated our method on two brain dMRI datasets, including a public dataset and an in-house dataset. Different settings were considered for the validation, and the results show that the proposed method improves the one-shot segmentation of novel WM tracts.

Comparison of Anatomical Pathway Models with Tractography Estimates of the Pallidothalamic, Cerebellothalamic, and Corticospinal Tracts

Introduction: Models of structural connectivity in the human brain are typically simulated using tractographic approaches. However, the nonlinear fitting of anatomical pathway atlases to de novo subject brains represents a simpler alternative that is hypothesized to provide more anatomically realistic results. Therefore, the goal of this study was to perform a side-by-side comparison of the streamline estimates generated by either pathway atlas fits or tractographic reconstructions in the same subjects. Methods: Our analyses focused on reconstruction of the corticospinal tract (CST), cerebellothalamic (CBT), and pallidothalamic (PT) pathways using example datasets from the Human Connectome Project (HCP). We used MRtrix3 to explore whole brain, as well as manual seed-to-target, tractography approaches. In parallel, we performed nonlinear fits of an axonal pathway atlas to each HCP dataset using Advanced Normalization Tools (ANTs). Results: The different methods produced notably different estimates for each pathway in each subject. The fitted atlas pathways were highly stereotyped and exhibited low variability in their streamline trajectories. Manual tractography resulted in pathway estimates that generally corresponded with the fitted atlas pathways, but with a higher degree of variability in the individual streamlines. Pathway reconstructions derived from whole-brain tractography exhibited the highest degree of variability and struggled to create anatomically realistic representations for either the CBT or PT pathways. Conclusion: The speed, simplicity, reproducibility, and realism of anatomical pathway model fits makes them an appealing option for some forms of structural connectivity modeling in the human brain. Impact statement Axonal pathway modeling is an important component of deep brain stimulation (DBS) research studies that seek to identify the brain connections that are directly activated by stimulation. The corticospinal tract, cerebellothalamic (CBT), and pallidothalamic (PT) pathways are specifically relevant to the study of subthalamic DBS for the treatment of Parkinson's disease. Our results suggest that anatomical pathway model fits of the CBT and PT pathways to de novo subject brains represent a more anatomically realistic option than tractographic approaches when studying subthalamic DBS.

BundleWarp, streamline-based nonlinear registration of white matter tracts

Nonlinear registration plays a central role in most neuroimage analysis methods and pipelines, such as in tractography-based individual and group-level analysis methods. However, nonlinear registration is a non-trivial task, especially when dealing with tractography data that digitally represent the underlying anatomy of the brain's white matter. Furthermore, such process often changes the structure of the data, causing artifacts that can suppress the underlying anatomical and structural details. In this paper, we introduce BundleWarp, a novel and robust streamline-based nonlinear registration method for the registration of white matter tracts. BundleWarp intelligently warps two bundles while preserving the bundles' crucial topological features. BundleWarp has two main steps. The first step involves the solution of an assignment problem that matches corresponding streamlines from the two bundles (iterLAP step). The second step introduces streamline-specific point-based deformations while keeping the topology of the bundle intact (mlCPD step). We provide comparisons against streamline-based linear registration and image-based nonlinear registration methods. BundleWarp quantitatively and qualitatively outperforms both, and we show that BundleWarp can deform and, at the same time, preserve important characteristics of the original anatomical shape of the bundles. Results are shown on 1,728 pairs of bundle registrations across 27 different bundle types. In addition, we present an application of BundleWarp for quantifying bundle shape differences using the generated deformation fields.

Structural and functional connectivity of the inferior temporal numeral area

A growing body of evidence suggests that in adults, there is a spatially consistent "inferior temporal numeral area" (ITNA) in the occipitotemporal cortex that appears to preferentially process Arabic digits relative to non-numerical symbols and objects. However, very little is known about why the ITNA is spatially segregated from regions that process other orthographic stimuli such as letters, and why it is spatially consistent across individuals. In the present study, we used diffusion-weighted imaging and functional magnetic resonance imaging to contrast structural and functional connectivity between left and right hemisphere ITNAs and a left hemisphere letter-preferring region. We found that the left ITNA had stronger structural and functional connectivity than the letter region to inferior parietal regions involved in numerical magnitude representation and arithmetic. Between hemispheres, the left ITNA showed stronger structural connectivity with the left inferior frontal gyrus (Broca's area), while the right ITNA showed stronger structural connectivity to the ipsilateral inferior parietal cortex and stronger functional coupling with the bilateral IPS. Based on their relative connectivity, our results suggest that the left ITNA may be more readily involved in mapping digits to verbal number representations, while the right ITNA may support the mapping of digits to quantity representations.

Structural disconnections associated with language impairments in chronic post-stroke aphasia using disconnectome maps

Inconsistent findings have been reported about the impact of structural disconnections on language function in post-stroke aphasia. This study investigated patterns of structural disconnections associated with chronic language impairments using disconnectome maps. Seventy-six individuals with post-stroke aphasia underwent a battery of language assessments and a structural MRI scan. Support-vector regression disconnectome-symptom mapping analyses were performed to examine the correlations between disconnectome maps, representing the probability of disconnection at each white matter voxel and different language scores. To further understand whether significant disconnections were primarily representing focal damage or a more extended network of seemingly preserved but disconnected areas beyond the lesion site, results were qualitatively compared to support-vector regression lesion-symptom mapping analyses. Part of the left white matter perisylvian network was similarly disconnected in 90% of the individuals with aphasia. Surrounding this common left perisylvian disconnectome, specific structural disconnections in the left fronto-temporo-parietal network were significantly associated with aphasia severity and with lower performance in auditory comprehension, syntactic comprehension, syntactic production, repetition and naming tasks. Auditory comprehension, repetition and syntactic processing deficits were related to disconnections in areas that overlapped with and extended beyond lesion sites significant in SVR-LSM analyses. In contrast, overall language abilities as measured by aphasia severity and naming seemed to be mostly explained by focal damage at the level of the insular and central opercular cortices, given the high overlap between SVR-DSM and SVR-LSM results for these scores. While focal damage seems to be sufficient to explain broad measures of language performance, the structural disconnections between language areas provide additional information on the neural basis of specific and persistent language impairments at the chronic stage beyond lesion volume. Leveraging routinely available clinical data, disconnectome mapping furthers our understanding of anatomical connectivity constraints that may limit the recovery of some language abilities in chronic post-stroke aphasia.

Population-based tract-to-region connectome of the human brain and its hierarchical topology

Connectome maps region-to-region connectivities but does not inform which white matter pathways form the connections. Here we constructed a population-based tract-to-region connectome to fill this information gap. The constructed connectome quantifies the population probability of a white matter tract innervating a cortical region. The results show that ~85% of the tract-to-region connectome entries are consistent across individuals, whereas the remaining (~15%) have substantial individual differences requiring individualized mapping. Further hierarchical clustering on cortical regions revealed dorsal, ventral, and limbic networks based on the tract-to-region connective patterns. The clustering results on white matter bundles revealed the categorization of fiber bundle systems in the association pathways. This tract-to-region connectome provides insights into the connective topology between cortical regions and white matter bundles. The derived hierarchical relation further offers a categorization of gray and white matter structures.

The brain connectome maps region-to-region connections but often ignores the role of the connecting pathways. Here, the authors mapped the tract-to-region relations to reveal the hierarchical relation of fiber bundles and dorsal, ventral, and limbic networks.

Learning white matter subject-specific segmentation from structural MRI

PURPOSE

Mapping brain white matter (WM) is essential for building an understanding of brain anatomy and function. Tractography-based methods derived from diffusion-weighted MRI (dMRI) are the principal tools for investigating WM. These procedures rely on time-consuming dMRI acquisitions that may not always be available, especially for legacy or time-constrained studies. To address this problem, we aim to generate WM tracts from structural magnetic resonance imaging (MRI) image by deep learning.

METHODS

Following recently proposed innovations in structural anatomical segmentation, we evaluate the feasibility of training multiply spatial localized convolution neural networks to learn context from fixed spatial patches from structural MRI on standard template. We focus on six widely used dMRI tractography algorithms (TractSeg, RecoBundles, XTRACT, Tracula, automated fiber quantification (AFQ), and AFQclipped) and train 125 U-Net models to learn these techniques from 3870 T1-weighted images from the Baltimore Longitudinal Study of Aging, the Human Connectome Project S1200 release, and scans acquired at Vanderbilt University.

RESULTS

The proposed framework identifies fiber bundles with high agreement against tractography-based pathways with a median Dice coefficient from 0.62 to 0.87 on a test cohort, achieving improved subject-specific accuracy when compared to population atlas-based methods. We demonstrate the generalizability of the proposed framework on three externally available datasets.

CONCLUSIONS

We show that patch-wise convolutional neural network can achieve robust bundle segmentation from T1w. We envision the use of this framework for visualizing the expected course of WM pathways when dMRI is not available.

TractEM: Evaluation of protocols for deterministic tractography white matter atlas

Reproducible identification of white matter pathways across subjects is essential for the study of structural connectivity of the human brain. One of the key challenges is anatomical differences between subjects and human rater subjectivity in labeling. Labeling white matter regions of interest presents many challenges due to the need to integrate both local and global information. Clearly communicating the manual processes to capture this information is cumbersome, yet essential to lay a solid foundation for comprehensive atlases. Segmentation protocols must be designed so the interpretation of the requested tasks as well as locating structural landmarks is anatomically accurate, intuitive and reproducible. In this work, we quantified the reproducibility of a first iteration of an open/public multi-bundle segmentation protocol. This allowed us to establish a baseline for its reproducibility as well as to identify the limitations for future iterations. The protocol was tested/evaluated on both typical 3 T research acquisition Baltimore Longitudinal Study of Aging (BLSA) and high-acquisition quality Human Connectome Project (HCP) datasets. The results show that a rudimentary protocol can produce acceptable intra-rater and inter-rater reproducibility. However, this work highlights the difficulty in generalizing reproducible results and the importance of reaching consensus on anatomical description of white matter pathways. The protocol has been made available in open source to improve generalizability and reliability in collaboration. The goal is to improve upon the first iteration and initiate a discussion on the anatomical validity (or lack thereof) of some bundle definitions and the importance of reproducibility of tractography segmentation.