Test-retest reproducibility of white matter parcellation using diffusion MRI tractography fiber clustering

A diffusion MRI tractography atlas for concurrent white matter mapping across Eastern and Western populations

The study of brain differences across Eastern and Western populations provides vital insights for understanding potential cultural and genetic influences on cognition and mental health. Diffusion MRI (dMRI) tractography is an important tool in assessing white matter (WM) connectivity and brain tissue microstructure across different populations. However, a comprehensive investigation into WM fiber tracts between Eastern and Western populations is challenged due to the lack of a cross-population WM atlas and the large site-specific variability of dMRI data. This study presents a dMRI tractography atlas, namely the East-West WM Atlas, for concurrent WM mapping between Eastern and Western populations and creates a large, harmonized dMRI dataset (n=306) based on the Human Connectome Project and the Chinese Human Connectome Project. The curated WM atlas, as well as subject-specific data including the harmonized dMRI data, the whole brain tractography data, and parcellated WM fiber tracts and their diffusion measures, are publicly released. This resource is a valuable addition to facilitating the exploration of brain commonalities and differences across diverse cultural backgrounds.

TractoSCR: a novel supervised contrastive regression framework for prediction of neurocognitive measures using multi-site harmonized diffusion MRI tractography

Neuroimaging-based prediction of neurocognitive measures is valuable for studying how the brain's structure relates to cognitive function. However, the accuracy of prediction using popular linear regression models is relatively low. We propose a novel deep regression method, namely TractoSCR, that allows full supervision for contrastive learning in regression tasks using diffusion MRI tractography. TractoSCR performs supervised contrastive learning by using the absolute difference between continuous regression labels (i.e., neurocognitive scores) to determine positive and negative pairs. We apply TractoSCR to analyze a large-scale dataset including multi-site harmonized diffusion MRI and neurocognitive data from 8,735 participants in the Adolescent Brain Cognitive Development (ABCD) Study. We extract white matter microstructural measures using a fine parcellation of white matter tractography into fiber clusters. Using these measures, we predict three scores related to domains of higher-order cognition (general cognitive ability, executive function, and learning/memory). To identify important fiber clusters for prediction of these neurocognitive scores, we propose a permutation feature importance method for high-dimensional data. We find that TractoSCR obtains significantly higher accuracy of neurocognitive score prediction compared to other state-of-the-art methods. We find that the most predictive fiber clusters are predominantly located within the superficial white matter and projection tracts, particularly the superficial frontal white matter and striato-frontal connections. Overall, our results demonstrate the utility of contrastive representation learning methods for regression, and in particular for improving neuroimaging-based prediction of higher-order cognitive abilities. Our code will be available at: https://github.com/SlicerDMRI/TractoSCR.

Short fiber bundle filtering and test-retest reproducibility of the Superficial White Matter

In recent years, there has been a growing interest in studying the Superficial White Matter (SWM). The SWM consists of short association fibers connecting near giry of the cortex, with a complex organization due to their close relationship with the cortical folding patterns. Therefore, their segmentation from dMRI tractography datasets requires dedicated methodologies to identify the main fiber bundle shape and deal with spurious fibers. This paper presents an enhanced short fiber bundle segmentation based on a SWM bundle atlas and the filtering of noisy fibers. The method was tuned and evaluated over HCP test-retest probabilistic tractography datasets (44 subjects). We propose four fiber bundle filters to remove spurious fibers. Furthermore, we include the identification of the main fiber fascicle to obtain well-defined fiber bundles. First, we identified four main bundle shapes in the SWM atlas, and performed a filter tuning in a subset of 28 subjects. The filter based on the Convex Hull provided the highest similarity between corresponding test-retest fiber bundles. Subsequently, we applied the best filter in the 16 remaining subjects for all atlas bundles, showing that filtered fiber bundles significantly improve test-retest reproducibility indices when removing between ten and twenty percent of the fibers. Additionally, we applied the bundle segmentation with and without filtering to the ABIDE-II database. The fiber bundle filtering allowed us to obtain a higher number of bundles with significant differences in fractional anisotropy, mean diffusivity, and radial diffusivity of Autism Spectrum Disorder patients relative to controls.

Brains of endurance athletes differ in the association areas but not in the primary areas

Regular participation in sports results in a series of physiological adaptations. However, little is known about the brain adaptations to physical activity. Here we aimed to investigate whether young endurance athletes and non-athletes differ in the gray and white matter of the brain and whether cardiorespiratory fitness (CRF) is associated with these differences. We assessed the CRF, volumes of the gray and white matter of the brain using structural magnetic resonance imaging (sMRI), and brain white matter connections using diffusion magnetic resonance imaging (dMRI) in 20 young male endurance athletes and 21 healthy non-athletes. While total brain volume was similar in both groups, the white matter volume was larger and the gray matter volume was smaller in the athletes compared to non-athletes. The reduction of gray matter was located in the association areas of the brain that are specialized in processing of sensory stimuli. In the microstructure analysis, significant group differences were found only in the association tracts, for example, the inferior occipito-frontal fascicle (IOFF) showing higher fractional anisotropy and lower radial diffusivity, indicating stronger myelination in this tract. Additionally, gray and white matter brain volumes, as well as association tracts correlated with CRF. No changes were observed in other brain areas or tracts. In summary, the brain signature of the endurance athlete is characterized by changes in the integration of sensory and motor information in the association areas.

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) Study® has collected data from over 10,000 children across 21 sites, providing insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a dataset of harmonized and processed ABCD dMRI data (from release 3) has been created, comprising quality-controlled imaging data from 9,345 subjects, focusing exclusively on the baseline session, i.e., the first time point of the study. This resource required substantial computational time (approx. 50,000 CPU hours) for harmonization, whole-brain tractography, and white matter parcellation. The dataset includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts in full and low resolution, and 804 different dMRI-derived measures per subject (72.3 TB total size). Accessible via the NIMH Data Archive, it offers a large-scale dMRI dataset for studying structural connectivity in child and adolescent neurodevelopment. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Tractography-based automated identification of the retinogeniculate visual pathway with novel microstructure-informed supervised contrastive learning

The retinogeniculate visual pathway (RGVP) is responsible for carrying visual information from the retina to the lateral geniculate nucleus. Identification and visualization of the RGVP are important in studying the anatomy of the visual system and can inform the treatment of related brain diseases. Diffusion MRI (dMRI) tractography is an advanced imaging method that uniquely enables in vivo mapping of the 3D trajectory of the RGVP. Currently, identification of the RGVP from tractography data relies on expert (manual) selection of tractography streamlines, which is time-consuming, has high clinical and expert labor costs, and is affected by inter-observer variability. In this paper, we present a novel deep learning framework, DeepRGVP , to enable fast and accurate identification of the RGVP from dMRI tractography data. We design a novel microstructure-informed supervised contrastive learning method that leverages both streamline label and tissue microstructure information to determine positive and negative pairs. We propose a simple and successful streamline-level data augmentation method to address highly imbalanced training data, where the number of RGVP streamlines is much lower than that of non-RGVP streamlines. We perform comparisons with several state-of-the-art deep learning methods that were designed for tractography parcellation, and we show superior RGVP identification results using DeepRGVP. In addition, we demonstrate a good generalizability of DeepRGVP to dMRI tractography data from neurosurgical patients with pituitary tumors and we show DeepRGVP can successfully identify RGVPs despite the effect of lesions affecting the RGVPs. Overall, our study shows the high potential of using deep learning to automatically identify the RGVP.

Multivariate Associations Among White Matter, Neurocognition, and Social Cognition Across Individuals With Schizophrenia Spectrum Disorders and Healthy Controls

BACKGROUND AND HYPOTHESIS

Neurocognitive and social cognitive abilities are important contributors to functional outcomes in schizophrenia spectrum disorders (SSDs). An unanswered question of considerable interest is whether neurocognitive and social cognitive deficits arise from overlapping or distinct white matter impairment(s).

STUDY DESIGN

We sought to fill this gap, by harnessing a large sample of individuals from the multi-center Social Processes Initiative in the Neurobiology of the Schizophrenia(s) (SPINS) dataset, unique in its collection of advanced diffusion imaging and an extensive battery of cognitive assessments. We applied canonical correlation analysis to estimates of white matter microstructure, and cognitive performance, across people with and without an SSD.

STUDY RESULTS

Our results established that white matter circuitry is dimensionally and strongly related to both neurocognition and social cognition, and that microstructure of the uncinate fasciculus and the rostral body of the corpus callosum may assume a "privileged role" subserving both. Further, we found that participant-wise estimates of white matter microstructure, weighted by cognitive performance, were largely consistent with participants' categorical diagnosis, and predictive of (cross-sectional) functional outcomes.

CONCLUSIONS

The demonstrated strength of the relationship between white matter circuitry and neurocognition and social cognition underscores the potential for using relationships among these variables to identify biomarkers of functioning, with potential prognostic and therapeutic implications.

High-Resolution Tractography Protocol to Investigate the Pathways between Human Mediodorsal Thalamic Nucleus and Prefrontal Cortex

Animal studies have established that the mediodorsal nucleus (MD) of the thalamus is heavily and reciprocally connected with all areas of the prefrontal cortex (PFC). In humans, however, these connections are difficult to investigate. High-resolution imaging protocols capable of reliably tracing the axonal tracts linking the human MD with each of the PFC areas may thus be key to advance our understanding of the variation, development, and plastic changes of these important circuits, in health and disease. Here, we tested in adult female and male humans the reliability of a new reconstruction protocol based on in vivo diffusion MRI to trace, measure, and characterize the fiber tracts interconnecting the MD with 39 human PFC areas per hemisphere. Our protocol comprised the following three components: (1) defining regions of interest; (2) preprocessing diffusion data; and, (3) modeling white matter tracts and tractometry. This analysis revealed largely separate PFC territories of reciprocal MD-PFC tracts bearing striking resemblance with the topographic layout observed in macaque connection-tracing studies. We then examined whether our protocol could reliably reconstruct each of these MD-PFC tracts and their profiles across test and retest sessions. Results revealed that this protocol was able to trace and measure, in both left and right hemispheres, the trajectories of these 39 area-specific axon bundles with good-to-excellent test-retest reproducibility. This protocol, which has been made publicly available, may be relevant for cognitive neuroscience and clinical studies of normal and abnormal PFC function, development, and plasticity.SIGNIFICANCE STATEMENT Reciprocal MD-PFC interactions are critical for complex human cognition and learning. Reliably tracing, measuring and characterizing MD-PFC white matter tracts using high-resolution noninvasive methods is key to assess individual variation of these systems in humans. Here, we propose a high-resolution tractography protocol that reliably reconstructs 39 area-specific MD-PFC white matter tracts per hemisphere and quantifies structural information from diffusion MRI data. This protocol revealed a detailed mapping of thalamocortical and corticothalamic MD-PFC tracts in four different PFC territories (dorsal, medial, orbital/frontal pole, inferior frontal) showing structural connections resembling those observed in tracing studies with macaques. Furthermore, our automated protocol revealed high test-retest reproducibility and is made publicly available, constituting a step forward in mapping human MD-PFC circuits in clinical and academic research.

Exploring the impact of hippocampal sclerosis on white matter tracts and memory in individuals with mesial temporal lobe epilepsy

OBJECTIVE

To investigate how the presence/side of hippocampal sclerosis (HS) are related to the white matter structure of cingulum bundle (CB), arcuate fasciculus (AF), and inferior longitudinal fasciculus (ILF) in mesial temporal lobe epilepsy (MTLE).

METHODS

We acquired diffusion-weighted magnetic resonance imaging (MRI) from 86 healthy and 71 individuals with MTLE (22 righ-HS; right-HS, 34 left-HS; left-HS, and 15 nonlesional MTLE). We utilized two-tensor tractography and fiber clustering to compare fractional anisotropy (FA) of each side/tract between groups. Additionally, we examined the association between FA and nonverbal (WMS-R) and verbal (WMS-R, RAVLT codification) memory performance for MTLE individuals.

RESULTS

White matter abnormalities depended on the side and presence of HS. The left-HS demonstrated widespread abnormalities for all tracts, the right-HS showed lower FA for ipsilateral tracts and the nonlesional MTLE group did not differ from healthy individuals. Results indicate no differences in verbal/nonverbal memory performance between the groups, but trend-level associations between higher FA of visual memory and the left CB (r = 0.286, P = 0.018), verbal memory (RAVLT) and -left CB (r = 0.335, P = 0.005), -right CB (r = 0.286, P = 0.016), and -left AF (r = 0.287, P = 0.017).

SIGNIFICANCE

Our results highlight that the presence and side of HS are crucial to understand the pathophysiology of MTLE. Specifically, left-sided HS seems to be related to widespread bilateral white matter abnormalities. Future longitudinal studies should focus on developing diagnostic and treatment strategies dependent on HS's presence/side.

PPA: Principal parcellation analysis for brain connectomes and multiple traits

Our understanding of the structure of the brain and its relationships with human traits is largely determined by how we represent the structural connectome. Standard practice divides the brain into regions of interest (ROIs) and represents the connectome as an adjacency matrix having cells measuring connectivity between pairs of ROIs. Statistical analyses are then heavily driven by the (largely arbitrary) choice of ROIs. In this article, we propose a human trait prediction framework utilizing a tractography-based representation of the brain connectome, which clusters fiber endpoints to define a data-driven white matter parcellation targeted to explain variation among individuals and predict human traits. This leads to Principal Parcellation Analysis (PPA), representing individual brain connectomes by compositional vectors building on a basis system of fiber bundles that captures the connectivity at the population level. PPA eliminates the need to choose atlases and ROIs a priori, and provides a simpler, vector-valued representation that facilitates easier statistical analysis compared to the complex graph structures encountered in classical connectome analyses. We illustrate the proposed approach through applications to data from the Human Connectome Project (HCP) and show that PPA connectomes improve power in predicting human traits over state-of-the-art methods based on classical connectomes, while dramatically improving parsimony and maintaining interpretability. Our PPA package is publicly available on GitHub, and can be implemented routinely for diffusion image data.

Deep Learning Methods for Identification of White Matter Fiber Tracts: Review of State-of-the-Art and Future Prospective

Quantitative analysis of white matter fiber tracts from diffusion Magnetic Resonance Imaging (dMRI) data is of great significance in health and disease. For example, analysis of fiber tracts related to anatomically meaningful fiber bundles is highly demanded in pre-surgical and treatment planning, and the surgery outcome depends on accurate segmentation of the desired tracts. Currently, this process is mainly done through time-consuming manual identification performed by neuro-anatomical experts. However, there is a broad interest in automating the pipeline such that it is fast, accurate, and easy to apply in clinical settings and also eliminates the intra-reader variabilities. Following the advancements in medical image analysis using deep learning techniques, there has been a growing interest in using these techniques for the task of tract identification as well. Recent reports on this application show that deep learning-based tract identification approaches outperform existing state-of-the-art methods. This paper presents a review of current tract identification approaches based on deep neural networks. First, we review the recent deep learning methods for tract identification. Next, we compare them with respect to their performance, training process, and network properties. Finally, we end with a critical discussion of open challenges and possible directions for future works.

The organization of frontostriatal brain wiring in non-affective early psychosis compared with healthy subjects using a novel diffusion imaging fiber cluster analysis

BACKGROUND

Alterations in brain connectivity may underlie neuropsychiatric conditions such as schizophrenia. We here assessed the degree of convergence of frontostriatal fiber projections in 56 young adult healthy controls (HCs) and 108 matched Early Psychosis-Non-Affective patients (EP-NAs) using our novel fiber cluster analysis of whole brain diffusion magnetic resonance imaging tractography.

METHODS

Using whole brain tractography and our fiber clustering methodology on harmonized diffusion magnetic resonance imaging data from the Human Connectome Project for Early Psychosis we identified 17 white matter fiber clusters that connect frontal cortex (FCtx) and caudate (Cd) per hemisphere in each group. To quantify the degree of convergence and, hence, topographical relationship of these fiber clusters, we measured the inter-cluster mean distances between the endpoints of the fiber clusters at the level of the FCtx and of the Cd, respectively.

RESULTS

We found (1) in both groups, bilaterally, a non-linear relationship, yielding convex curves, between FCtx and Cd distances for FCtx-Cd connecting fiber clusters, driven by a cluster projecting from inferior frontal gyrus; however, in the right hemisphere, the convex curve was more flattened in EP-NAs; (2) that cluster pairs in the right (p = 0.03), but not left (p = 0.13), hemisphere were significantly more convergent in HCs vs EP-NAs; (3) in both groups, bilaterally, similar clusters projected significantly convergently to the Cd; and, (4) a significant group by fiber cluster pair interaction for 2 right hemisphere fiber clusters (numbers 5, 11; p = .00023; p = .00023) originating in selective PFC subregions.

CONCLUSIONS

In both groups, we found the FCtx-Cd wiring pattern deviated from a strictly topographic relationship and that similar clusters projected significantly more convergently to the Cd. Interestingly, we also found a significantly more convergent pattern of connectivity in HCs in the right hemisphere and that 2 clusters from PFC subregions in the right hemisphere significantly differed in their pattern of connectivity between groups.

A Proposed Human Structural Brain Connectivity Matrix in the Center for Morphometric Analysis Harvard-Oxford Atlas Framework: A Historical Perspective and Future Direction for Enhancing the Precision of Human Structural Connectivity with a Novel Neuroanatomical Typology

A complete structural definition of the human nervous system must include delineation of its wiring diagram (e.g., Swanson LW. Brain architecture: understanding the basic plan, 2012). The complete formulation of the human brain circuit diagram (BCD [Front Neuroanat. 2020;14:18]) has been hampered by an inability to determine connections in their entirety (i.e., not only pathway stems but also origins and terminations). From a structural point of view, a neuroanatomic formulation of the BCD should include the origins and terminations of each fiber tract as well as the topographic course of the fiber tract in three dimensions. Classic neuroanatomical studies have provided trajectory information for pathway stems and their speculative origins and terminations [Dejerine J and Dejerine-Klumpke A. Anatomie des Centres Nerveux, 1901; Dejerine J and Dejerine-Klumpke A. Anatomie des Centres Nerveux: Méthodes générales d'étude-embryologie-histogénèse et histologie. Anatomie du cerveau, 1895; Ludwig E and Klingler J. Atlas cerebri humani, 1956; Makris N. Delineation of human association fiber pathways using histologic and magnetic resonance methodologies; 1999; Neuroimage. 1999 Jan;9(1):18-45]. We have summarized these studies previously [Neuroimage. 1999 Jan;9(1):18-45] and present them here in a macroscale-level human cerebral structural connectivity matrix. A matrix in the present context is an organizational construct that embodies anatomical knowledge about cortical areas and their connections. This is represented in relation to parcellation units according to the Harvard-Oxford Atlas neuroanatomical framework established by the Center for Morphometric Analysis at Massachusetts General Hospital in the early 2000s, which is based on the MRI volumetrics paradigm of Dr. Verne Caviness and colleagues [Brain Dev. 1999 Jul;21(5):289-95]. This is a classic connectional matrix based mainly on data predating the advent of DTI tractography, which we refer to as the "pre-DTI era" human structural connectivity matrix. In addition, we present representative examples that incorporate validated structural connectivity information from nonhuman primates and more recent information on human structural connectivity emerging from DTI tractography studies. We refer to this as the "DTI era" human structural connectivity matrix. This newer matrix represents a work in progress and is necessarily incomplete due to the lack of validated human connectivity findings on origins and terminations as well as pathway stems. Importantly, we use a neuroanatomical typology to characterize different types of connections in the human brain, which is critical for organizing the matrices and the prospective database. Although substantial in detail, the present matrices may be assumed to be only partially complete because the sources of data relating to human fiber system organization are limited largely to inferences from gross dissections of anatomic specimens or extrapolations of pathway tracing information from nonhuman primate experiments [Front Neuroanat. 2020;14:18, Front Neuroanat. 2022;16:1035420, and Brain Imaging Behav. 2021;15(3):1589-1621]. These matrices, which embody a systematic description of cerebral connectivity, can be used in cognitive and clinical studies in neuroscience and, importantly, to guide research efforts for further elucidating, validating, and completing the human BCD [Front Neuroanat. 2020;14:18].

Deep Diffusion MRI Registration (DDMReg): A Deep Learning Method for Diffusion MRI Registration

In this paper, we present a deep learning method, DDMReg, for accurate registration between diffusion MRI (dMRI) datasets. In dMRI registration, the goal is to spatially align brain anatomical structures while ensuring that local fiber orientations remain consistent with the underlying white matter fiber tract anatomy. DDMReg is a novel method that uses joint whole-brain and tract-specific information for dMRI registration. Based on the successful VoxelMorph framework for image registration, we propose a novel registration architecture that leverages not only whole brain information but also tract-specific fiber orientation information. DDMReg is an unsupervised method for deformable registration between pairs of dMRI datasets: it does not require nonlinearly pre-registered training data or the corresponding deformation fields as ground truth. We perform comparisons with four state-of-the-art registration methods on multiple independently acquired datasets from different populations (including teenagers, young and elderly adults) and different imaging protocols and scanners. We evaluate the registration performance by assessing the ability to align anatomically corresponding brain structures and ensure fiber spatial agreement between different subjects after registration. Experimental results show that DDMReg obtains significantly improved registration performance compared to the state-of-the-art methods. Importantly, we demonstrate successful generalization of DDMReg to dMRI data from different populations with varying ages and acquired using different acquisition protocols and different scanners.

Brainstem Diffusion Tensor Tractography and Clinical Applications in Pain

The brainstem is one of the most vulnerable brain structures in many neurological conditions, such as pain, sleep problems, autonomic dysfunctions, and neurodegenerative disorders. Diffusion tensor imaging and tractography provide structural details and quantitative measures of brainstem fiber pathways. Until recently, diffusion tensor tractographic studies have mainly focused on whole-brain MRI acquisition. Due to the brainstem's spatial localization, size, and tissue characteristics, and limits of imaging techniques, brainstem diffusion MRI poses particular challenges in tractography. We provide a brief overview on recent advances in diffusion tensor tractography in revealing human pathways connecting the brainstem to the subcortical regions (e.g., basal ganglia, mesolimbic, basal forebrain), and cortical regions. Each of these pathways contains different distributions of fiber tracts from known neurotransmitter-specific nuclei in the brainstem. We compare the brainstem tractographic approaches in literature and our in-lab developed automated brainstem tractography in terms of atlas building, technical advantages, and neuroanatomical implications on neurotransmitter systems. Lastly, we summarize recent investigations of using brainstem tractography as a promising tool in association with pain.

Connectivity-based Cortical Parcellation via Contrastive Learning on Spatial-Graph Convolution

Objective. Objective of this work is the development and evaluation of a cortical parcellation framework based on tractography-derived brain structural connectivity. Impact Statement. The proposed framework utilizes novel spatial-graph representation learning methods for solving the task of cortical parcellation, an important medical image analysis and neuroscientific problem. Introduction. The concept of "connectional fingerprint" has motivated many investigations on the connectivity-based cortical parcellation, especially with the technical advancement of diffusion imaging. Previous studies on multiple brain regions have been conducted with promising results. However, performance and applicability of these models are limited by the relatively simple computational scheme and the lack of effective representation of brain imaging data. Methods. We propose the Spatial-graph Convolution Parcellation (SGCP) framework, a two-stage deep learning-based modeling for the graph representation brain imaging. In the first stage, SGCP learns an effective embedding of the input data through a self-supervised contrastive learning scheme with the backbone encoder of a spatial-graph convolution network. In the second stage, SGCP learns a supervised classifier to perform voxel-wise classification for parcellating the desired brain region. Results. SGCP is evaluated on the parcellation task for 5 brain regions in a 15-subject DWI dataset. Performance comparisons between SGCP, traditional parcellation methods, and other deep learning-based methods show that SGCP can achieve superior performance in all the cases. Conclusion. Consistent good performance of the proposed SGCP framework indicates its potential to be used as a general solution for investigating the regional/subregional composition of human brain based on one or more connectivity measurements.

Assessing the Reliability of Template-Based Clustering for Tractography in Healthy Human Adults

Tractography is a non-invasive technique to investigate the brain’s structural pathways (also referred to as tracts) that connect different brain regions. A commonly used approach for identifying tracts is with template-based clustering, where unsupervised clustering is first performed on a template in order to label corresponding tracts in unseen data. However, the reliability of this approach has not been extensively studied. Here, an investigation into template-based clustering reliability was performed, assessing the output from two datasets: Human Connectome Project (HCP) and MyConnectome project. The effect of intersubject variability on template-based clustering reliability was investigated, as well as the reliability of both deep and superficial white matter tracts. Identified tracts were evaluated by assessing Euclidean distances from a dataset-specific tract average centroid, the volumetric overlap across corresponding tracts, and along-tract agreement of quantitative values. Further, two template-based techniques were employed to evaluate the reliability of different clustering approaches. Reliability assessment can increase the confidence of a tract identifying technique in future applications to study pathways of interest. The two different template-based approaches exhibited similar reliability for identifying both deep white matter tracts and the superficial white matter.

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.

Automated facial-vestibulocochlear nerve complex identification based on data-driven tractography clustering

Small size and intricate anatomical environment are the main difficulties facing tractography of the facial-vestibulocochlear nerve complex (FVN), and lead to challenges in fiber orientation distribution (FOD) modeling, fiber tracking, region-of-interest selection, and fiber filtering. Experts need rich experience in anatomy and tractography, as well as substantial labor costs, to identify the FVN. Thus, we present a pipeline to identify the FVN automatically, in what we believe is the first study of the automated identification of the FVN. First, we created an FVN template. Forty high-resolution multishell data were used to perform data-driven fiber clustering based on the multishell multitissue constraint spherical deconvolution FOD model and deterministic tractography. We selected the brainstem and cerebellum (BS-CB) region as the seed region and removed the fibers that reach other brain regions. We then performed spectral fiber clustering twice. The first clustering was to create a BS-CB atlas and separate the fibers that pass through the cerebellopontine angle, and the other one was to extract the FVN. Second, we registered the subject-specific fibers in the space of the FVN template and assigned each fiber to the closest cluster to identify the FVN automatically by spectral embedding. We applied the proposed method to different acquirement sites, including two different healthy datasets and two tumor patient datasets. Experimental results showed that our automatic identification results have ideal colocalization with expert manual identification in terms of spatial overlap and visualization. Importantly, we successfully applied our method to tumor patient data. The FVNs identified by the proposed method were in agreement with intraoperative findings.

The Organization of Frontostriatal Brain Wiring in Healthy Subjects Using a Novel Diffusion Imaging Fiber Cluster Analysis

To assess normal organization of frontostriatal brain wiring, we analyzed diffusion magnetic resonance imaging (dMRI) scans in 100 young adult healthy subjects (HSs). We identified fiber clusters intersecting the frontal cortex and caudate, a core component of associative striatum, and quantified their degree of deviation from a strictly topographic pattern. Using whole brain dMRI tractography and an automated tract parcellation clustering method, we extracted 17 white matter fiber clusters per hemisphere connecting the frontal cortex and caudate. In a novel approach to quantify the geometric relationship among clusters, we measured intercluster endpoint distances between corresponding cluster pairs in the frontal cortex and caudate. We show first, the overall frontal cortex wiring pattern of the caudate deviates from a strictly topographic organization due to significantly greater convergence in regionally specific clusters; second, these significantly convergent clusters originate in subregions of ventrolateral, dorsolateral, and orbitofrontal prefrontal cortex (PFC); and, third, a similar organization in both hemispheres. Using a novel tractography method, we find PFC-caudate brain wiring in HSs deviates from a strictly topographic organization due to a regionally specific pattern of cluster convergence. We conjecture cortical subregions projecting to the caudate with greater convergence subserve functions that benefit from greater circuit integration.

Comparison of multiple tractography methods for reconstruction of the retinogeniculate visual pathway using diffusion MRI

The retinogeniculate visual pathway (RGVP) conveys visual information from the retina to the lateral geniculate nucleus. The RGVP has four subdivisions, including two decussating and two nondecussating pathways that cannot be identified on conventional structural magnetic resonance imaging (MRI). Diffusion MRI tractography has the potential to trace these subdivisions and is increasingly used to study the RGVP. However, it is not yet known which fiber tracking strategy is most suitable for RGVP reconstruction. In this study, four tractography methods are compared, including constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, and multi-fiber (UKF-2T) and single-fiber (UKF-1T) unscented Kalman filter (UKF) methods. Experiments use diffusion MRI data from 57 subjects in the Human Connectome Project. The RGVP is identified using regions of interest created by two clinical experts. Quantitative anatomical measurements and expert anatomical judgment are used to assess the advantages and limitations of the four tractography methods. Overall, we conclude that UKF-2T and iFOD1 produce the best RGVP reconstruction results. The iFOD1 method can better quantitatively estimate the percentage of decussating fibers, while the UKF-2T method produces reconstructed RGVPs that are judged to better correspond to the known anatomy and have the highest spatial overlap across subjects. Overall, we find that it is challenging for current tractography methods to both accurately track RGVP fibers that correspond to known anatomy and produce an approximately correct percentage of decussating fibers. We suggest that future algorithm development for RGVP tractography should take consideration of both of these two points.

MK-Curve improves sensitivity to identify white matter alterations in clinical high risk for psychosis

Diffusion kurtosis imaging (DKI) is a diffusion MRI approach that enables the measurement of brain microstructural properties, reflecting molecular restrictions and tissue heterogeneity. DKI parameters such as mean kurtosis (MK) provide additional subtle information to that provided by popular diffusion tensor imaging (DTI) parameters, and thus have been considered useful to detect white matter abnormalities, especially in populations that are not expected to show severe brain pathologies. However, DKI parameters often yield artifactual output values that are outside of the biologically plausible range, which diminish sensitivity to identify true microstructural changes. Recently we have proposed the mean-kurtosis-curve (MK-Curve) method to correct voxels with implausible DKI parameters, and demonstrated its improved performance against other approaches that correct artifacts in DKI. In this work, we aimed to evaluate the utility of the MK-Curve method to improve the identification of white matter abnormalities in group comparisons. To do so, we compared group differences, with and without the MK-Curve correction, between 115 individuals at clinical high risk for psychosis (CHR) and 93 healthy controls (HCs). We also compared the correlation of the corrected and uncorrected DKI parameters with clinical characteristics. Following the MK-curve correction, the group differences had larger effect sizes and higher statistical significance (i.e., lower p-values), demonstrating increased sensitivity to detect group differences, in particular in MK. Furthermore, the MK-curve-corrected DKI parameters displayed stronger correlations with clinical variables in CHR individuals, demonstrating the clinical relevance of the corrected parameters. Overall, following the MK-curve correction our analyses found widespread lower MK in CHR that overlapped with lower fractional anisotropy (FA), and both measures were significantly correlated with a decline in functioning and with more severe symptoms. These observations further characterize white matter alterations in the CHR stage, demonstrating that MK and FA abnormalities are widespread, and mostly overlap. The improvement in group differences and stronger correlation with clinical variables suggest that applying MK-curve would be beneficial for the detection and characterization of subtle group differences in other experiments as well.

Sex-Related Differences in White Matter Asymmetry and Its Implications for Verbal Working Memory in Psychosis High-Risk State

Objective: Sexual dimorphism has been investigated in schizophrenia, although sex-specific differences among individuals who are at clinical high-risk (CHR) for developing psychosis have been inconclusive. This study aims to characterize sexual dimorphism of language areas in the brain by investigating the asymmetry of four white matter tracts relevant to verbal working memory in CHR patients compared to healthy controls (HC). HC typically show a leftward asymmetry of these tracts. Moreover, structural abnormalities in asymmetry and verbal working memory dysfunctions have been associated with neurodevelopmental abnormalities and are considered core features of schizophrenia. Methods: Twenty-nine subjects with CHR (17 female/12 male) for developing psychosis and twenty-one HC (11 female/10 male) matched for age, sex, and education were included in the study. Two-tensor unscented Kalman filter tractography, followed by an automated, atlas-guided fiber clustering approach, were used to identify four fiber tracts related to verbal working memory: the superior longitudinal fasciculi (SLF) I, II and III, and the superior occipitofrontal fasciculus (SOFF). Using fractional anisotropy (FA) of tissue as the primary measure, we calculated the laterality index for each tract. Results: There was a significantly greater right>left asymmetry of the SLF-III in CHR females compared to HC females, but no hemispheric difference between CHR vs. HC males. Moreover, the laterality index of SLF-III for CHR females correlated negatively with Backward Digit Span performance, suggesting a greater rightward asymmetry was associated with poorer working memory functioning. Conclusion: This study suggests increased rightward asymmetry of the SLF-III in CHR females. This finding of sexual dimorphism in white matter asymmetry in a language-related area of the brain in CHR highlights the need for a deeper understanding of the role of sex in the high-risk state. Future work investigating early sex-specific pathophysiological mechanisms, may lead to the development of novel personalized treatment strategies aimed at preventing transition to a more chronic and difficult-to-treat disorder.

Microstructural alterations in medial forebrain bundle are associated with interindividual pain sensitivity

The perception of pain to noxious stimuli, also known as pain sensitivity, varies among individuals. The comprised brain structures and their white matter pathways are complex and elusive. Here, we aimed to investigate whether variation of microstructure of the medial forebrain bundle (MFB), a tract connecting the basal forebrain with the brain stem, is associated with interindividual pain sensitivity. We assessed interindividual pain sensitivity as a rating of pain intensity to heat stimuli (45, 47, and 48.9°C) in 38 healthy men (age: 27.05 ± 5.7 years). We also reconstructed the MFB using multitensor tractography from diffusion magnetic resonance imaging (dMRI) and calculated free‐water corrected dMRI measures of fractional anisotropy (FA_t), radial diffusivity (RD_t), and axial diffusivity (AD_t). Lower ratings of interindividual pain intensity correlated with higher FA_t and lower RD_t of the MFB. As changes in FA_t and RD_t may reflect abnormalities in myelination, the results might be interpreted as that a lower pain rating is associated with higher degree of myelination of the MFB and could represent an inhibitory pathway of pain. Our results suggest that alteration of microstructure in the MFB contributes to the interindividual variation of pain perception.

Creation of a novel trigeminal tractography atlas for automated trigeminal nerve identification

Diffusion MRI (dMRI) tractography has been successfully used to study the trigeminal nerves (TGNs) in many clinical and research applications. Currently, identification of the TGN in tractography data requires expert nerve selection using manually drawn regions of interest (ROIs), which is prone to inter-observer variability, time-consuming and carries high clinical and labor costs. To overcome these issues, we propose to create a novel anatomically curated TGN tractography atlas that enables automated identification of the TGN from dMRI tractography. In this paper, we first illustrate the creation of a trigeminal tractography atlas. Leveraging a well-established computational pipeline and expert neuroanatomical knowledge, we generate a data-driven TGN fiber clustering atlas using tractography data from 50 subjects from the Human Connectome Project. Then, we demonstrate the application of the proposed atlas for automated TGN identification in new subjects, without relying on expert ROI placement. Quantitative and visual experiments are performed with comparison to expert TGN identification using dMRI data from two different acquisition sites. We show highly comparable results between the automatically and manually identified TGNs in terms of spatial overlap and visualization, while our proposed method has several advantages. First, our method performs automated TGN identification, and thus it provides an efficient tool to reduce expert labor costs and inter-operator bias relative to expert manual selection. Second, our method is robust to potential imaging artifacts and/or noise that can prevent successful manual ROI placement for TGN selection and hence yields a higher successful TGN identification rate.

Deep white matter analysis (DeepWMA): Fast and consistent tractography segmentation

White matter tract segmentation, i.e. identifying tractography fibers (streamline trajectories) belonging to anatomically meaningful fiber tracts, is an essential step to enable tract quantification and visualization. In this study, we present a deep learning tractography segmentation method (DeepWMA) that allows fast and consistent identification of 54 major deep white matter fiber tracts from the whole brain. We create a large-scale training tractography dataset of 1 million labeled fiber samples, and we propose a novel 2D multi-channel feature descriptor (FiberMap) that encodes spatial coordinates of points along each fiber. We learn a convolutional neural network (CNN) fiber classification model based on FiberMap and obtain a high fiber classification accuracy of 90.99% on the training tractography data with ground truth fiber labels. Then, the method is evaluated on a test dataset of 597 diffusion MRI scans from six independently acquired populations across genders, the lifespan (1 day - 82 years), and different health conditions (healthy control, neuropsychiatric disorders, and brain tumor patients). We perform comparisons with two state-of-the-art tract segmentation methods. Experimental results show that our method obtains a highly consistent tract segmentation result, where on average over 99% of the fiber tracts are successfully identified across all subjects under study, most importantly, including neonates and patients with space-occupying brain tumors. We also demonstrate good generalization of the method to tractography data from multiple different fiber tracking methods. The proposed method leverages deep learning techniques and provides a fast and efficient tool for brain white matter segmentation in large diffusion MRI tractography datasets.

GeoSP: A parallel method for a cortical surface parcellation based on geodesic distance

We present GeoSP, a parallel method that creates a parcellation of the cortical mesh based on a geodesic distance, in order to consider gyri and sulci topology. The method represents the mesh with a graph and performs a K-means clustering in parallel. It has two modes of use, by default, it performs the geodesic cortical parcellation based on the boundaries of the anatomical parcels provided by the Desikan-Killiany atlas. The other mode performs the complete parcellation of the cortex. Results for both modes and with different values for the total number of sub-parcels show homogeneous sub-parcels. Furthermore, the execution time is 82s for the whole cortex mode and 18s for the Desikan-Killiany atlas subdivision, for a parcellation into 350 sub-parcels. The proposed method will be available to the community to perform the evaluation of data-driven cortical parcellations. As an example, we compared GeoSP parcellation with Desikan-Killiany and Destrieux atlases in 50 subjects, obtaining more homogeneous parcels for GeoSP and minor differences in structural connectivity reproducibility across subjects.

Test-retest reproducibility of a multi-atlas automated segmentation tool on multimodality brain MRI

INTRODUCTION

The increasing use of large sample sizes for population and personalized medicine requires high-throughput tools for imaging processing that can handle large amounts of data with diverse image modalities, perform a biologically meaningful information reduction, and result in comprehensive quantification. Exploring the reproducibility of these tools reveals the specific strengths and weaknesses that heavily influence the interpretation of results, contributing to transparence in science.

METHODS

We tested-retested the reproducibility of MRICloud, a free automated method for whole-brain, multimodal MRI segmentation and quantification, on two public, independent datasets of healthy adults.

RESULTS

The reproducibility was extremely high for T1-volumetric analysis, high for diffusion tensor images (DTI) (however, regionally variable), and low for resting-state fMRI.

CONCLUSION

In general, the reproducibility of the different modalities was slightly superior to that of widely used software. This analysis serves as a normative reference for planning samples and for the interpretation of structure-based MRI studies.

Test-retest reproducibility of white matter parcellation using diffusion MRI tractography fiber clustering

There are two popular approaches for automated white matter parcellation using diffusion MRI tractography, including fiber clustering strategies that group white matter fibers according to their geometric trajectories and cortical-parcellation-based strategies that focus on the structural connectivity among different brain regions of interest. While multiple studies have assessed test-retest reproducibility of automated white matter parcellations using cortical-parcellation-based strategies, there are no existing studies of test-retest reproducibility of fiber clustering parcellation. In this work, we perform what we believe is the first study of fiber clustering white matter parcellation test-retest reproducibility. The assessment is performed on three test-retest diffusion MRI datasets including a total of 255 subjects across genders, a broad age range (5-82 years), health conditions (autism, Parkinson's disease and healthy subjects), and imaging acquisition protocols (three different sites). A comprehensive evaluation is conducted for a fiber clustering method that leverages an anatomically curated fiber clustering white matter atlas, with comparison to a popular cortical-parcellation-based method. The two methods are compared for the two main white matter parcellation applications of dividing the entire white matter into parcels (i.e., whole brain white matter parcellation) and identifying particular anatomical fiber tracts (i.e., anatomical fiber tract parcellation). Test-retest reproducibility is measured using both geometric and diffusion features, including volumetric overlap (wDice) and relative difference of fractional anisotropy. Our experimental results in general indicate that the fiber clustering method produced more reproducible white matter parcellations than the cortical-parcellation-based method.

Automated fiber tract reconstruction for surgery planning: Extensive validation in language-related white matter tracts

Diffusion MRI and tractography hold great potential for surgery planning, especially to preserve eloquent white matter during resections. However, fiber tract reconstruction requires an expert with detailed understanding of neuroanatomy. Several automated approaches have been proposed, using different strategies to reconstruct the white matter tracts in a supervised fashion. However, validation is often limited to comparison with manual delineation by overlap-based measures, which is limited in characterizing morphological and topological differences. In this work, we set up a fully automated pipeline based on anatomical criteria that does not require manual intervention, taking advantage of atlas-based criteria and advanced acquisition protocols available on clinical-grade MRI scanners. Then, we extensively validated it on epilepsy patients with specific focus on language-related bundles. The validation procedure encompasses different approaches, including simple overlap with manual segmentations from two experts, feasibility ratings from external multiple clinical raters and relation with task-based functional MRI. Overall, our results demonstrate good quantitative agreement between automated and manual segmentation, in most cases better performances of the proposed method in qualitative terms, and meaningful relationships with task-based fMRI. In addition, we observed significant differences between experts in terms of both manual segmentation and external ratings. These results offer important insights on how different levels of validation complement each other, supporting the idea that overlap-based measures, although quantitative, do not offer a full perspective on the similarities and differences between automated and manual methods.