Strengths and weaknesses of state of the art fiber tractography pipelines--A comprehensive in-vivo and phantom evaluation study using Tractometer

What matters in reinforcement learning for tractography

Recently, deep reinforcement learning (RL) has been proposed to learn the tractography procedure and train agents to reconstruct the structure of the white matter without manually curated reference streamlines. While the performances reported were competitive, the proposed framework is complex, and little is still known about the role and impact of its multiple parts. In this work, we thoroughly explore the different components of the proposed framework, such as the choice of the RL algorithm, seeding strategy, the input signal and reward function, and shed light on their impact. Approximately 7,400 models were trained for this work, totalling nearly 41,000 h of GPU time. Our goal is to guide researchers eager to explore the possibilities of deep RL for tractography by exposing what works and what does not work with the category of approach. As such, we ultimately propose a series of recommendations concerning the choice of RL algorithm, the input to the agents, the reward function and more to help future work using reinforcement learning for tractography. We also release the open source codebase, trained models, and datasets for users and researchers wanting to explore reinforcement learning for tractography.

Q-ball high-resolution fiber tractography of language associated tracts: quantitative evaluation of applicability for glioma resections

BACKGROUND

To date, fiber tractography (FT) is predominantly based on diffusion tensor imaging (DTI). High angular resolution diffusion imaging (HARDI)-based reconstructions have become a focus of interest, enabling the resolution of intravoxel fiber crossing. However, experience with high resolution tractography (HRFT) for neurosurgical applications is still limited to a few reports. This prospectively designed feasibility study shares our initial experience using an analytical q-ball approach (QBI) for FT of language-associated pathways in comparison with DTI-FT, focussing on a quantitative analysis and evaluation of its applicability in clinical routine.

METHODS

Probabilistic QBI-, and DTI-FT were performed for the major components of the language-associated fiber bundles (superior longitudinal fasciculus, inferior fronto-occipital fasciculus, medial/inferior longitudinal faciculus) in 11 patients with eloquent gliomas. The data was derived from a routine DWI sequence (b=1000s/mm2, 64 gradient directions). Quantitative analysis evaluated tract volume (TV), tract length (TL) and tract density (TD). Results were correlated to tumor and edema size.

RESULTS

Quantitative analysis showed larger TV and TL of the overall fiber object using QBI-FT compared with DTI-FT (TV: 16.45±1.85 vs. 10.07±1.15cm3; P<0.0001; TL: 81.95±6.14 vs. 72.06±6.92 mm; P=0.0011). Regarding overall TD, DTI delivered significantly higher values (40.57±6.59 vs. 60.98±15.94 points/voxel; P=0.0118). Bland-Altman analysis illustrated a systematic advantage to yield lager TV and TL via QBI compared with DTI for all reconstructed pathways. The results were independent of tumor or edema volume.

CONCLUSIONS

QBI proved to be suitable for an application in the neurosurgical setting without additional expense for the patient. Quantitative analysis of FT reveals larger overall TV, longer TL with lower TD using QBI compared with DTI, suggesting the better depiction of marginal and terminal fibers according to neuroanatomical knowledge. This emphasizes the known limitation of DTI to underestimate the dimensions of a pathway. Rather than relying on DTI, sophisticated HRFT techniques should be considered for preoperative planning and intraoperative guidance in selected cases of eloquent glioma surgery.

Supervised tractogram filtering using Geometric Deep Learning

A tractogram is a virtual representation of the brain white matter. It is composed of millions of virtual fibers, encoded as 3D polylines, which approximate the white matter axonal pathways. To date, tractograms are the most accurate white matter representation and thus are used for tasks like presurgical planning and investigations of neuroplasticity, brain disorders, or brain networks. However, it is a well-known issue that a large portion of tractogram fibers is not anatomically plausible and can be considered artifacts of the tracking procedure. With Verifyber, we tackle the problem of filtering out such non-plausible fibers using a novel fully-supervised learning approach. Differently from other approaches based on signal reconstruction and/or brain topology regularization, we guide our method with the existing anatomical knowledge of the white matter. Using tractograms annotated according to anatomical principles, we train our model, Verifyber, to classify fibers as either anatomically plausible or non-plausible. The proposed Verifyber model is an original Geometric Deep Learning method that can deal with variable size fibers, while being invariant to fiber orientation. Our model considers each fiber as a graph of points, and by learning features of the edges between consecutive points via the proposed sequence Edge Convolution, it can capture the underlying anatomical properties. The output filtering results highly accurate and robust across an extensive set of experiments, and fast; with a 12GB GPU, filtering a tractogram of 1M fibers requires less than a minute.

Reconstructing the somatotopic organization of the corticospinal tract remains a challenge for modern tractography methods

The corticospinal tract (CST) is a critically important white matter fiber tract in the human brain that enables control of voluntary movements of the body. The CST exhibits a somatotopic organization, which means that the motor neurons that control specific body parts are arranged in order within the CST. Diffusion magnetic resonance imaging (MRI) tractography is increasingly used to study the anatomy of the CST. However, despite many advances in tractography algorithms over the past decade, modern, state-of-the-art methods still face challenges. In this study, we compare the performance of six widely used tractography methods for reconstructing the CST and its somatotopic organization. These methods include constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, unscented Kalman filter (UKF) tractography methods including multi-fiber (UKF2T) and single-fiber (UKF1T) models, the generalized q-sampling imaging (GQI) based deterministic tractography method, and the TractSeg method. We investigate CST somatotopy by dividing the CST into four subdivisions per hemisphere that originate in the leg, trunk, hand, and face areas of the primary motor cortex. A quantitative and visual comparison is performed using diffusion MRI data (N = 100 subjects) from the Human Connectome Project. Quantitative evaluations include the reconstruction rate of the eight anatomical subdivisions, the percentage of streamlines in each subdivision, and the coverage of the white matter-gray matter (WM-GM) interface. CST somatotopy is further evaluated by comparing the percentage of streamlines in each subdivision to the cortical volumes for the leg, trunk, hand, and face areas. Overall, UKF2T has the highest reconstruction rate and cortical coverage. It is the only method with a significant positive correlation between the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex. However, our experimental results show that all compared tractography methods are biased toward generating many trunk streamlines (ranging from 35.10% to 71.66% of total streamlines across methods). Furthermore, the coverage of the WM-GM interface in the largest motor area (face) is generally low (under 40%) for all compared tractography methods. Different tractography methods give conflicting results regarding the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex, indicating that there is generally no clear relationship, and that reconstruction of CST somatotopy is still a large challenge. Overall, we conclude that while current tractography methods have made progress toward the well-known challenge of improving the reconstruction of the lateral projections of the CST, the overall problem of performing a comprehensive CST reconstruction, including clinically important projections in the lateral (hand and face areas) and medial portions (leg area), remains an important challenge for diffusion MRI tractography.

Quantitative mapping of the brain's structural connectivity using diffusion MRI tractography: A review

Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain's white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain's structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain's structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain's white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.

A novel male Japanese quail structural connectivity atlas using ultra-high field diffusion MRI at 11.7 T

The structural connectivity of animal brains can be revealed using post-mortem diffusion-weighted magnetic resonance imaging (MRI). Despite the existence of several structural atlases of avian brains, few of them address the bird’s structural connectivity. In this study, a novel atlas of the structural connectivity is proposed for the male Japanese quail (Coturnix japonica), aiming at investigating two lines divergent on their emotionality trait: the short tonic immobility (STI) and the long tonic immobility (LTI) lines. The STI line presents a low emotionality trait, while the LTI line expresses a high emotionality trait. 21 male Japanese quail brains from both lines were scanned post-mortem for this study, using a preclinical Bruker 11.7 T MRI scanner. Diffusion-weighted MRI was performed using a 3D segmented echo planar imaging (EPI) pulsed gradient spin-echo (PGSE) sequence with a 200 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution, 75 diffusion-encoding directions and a b-value fixed at 4500 s/mm². Anatomical MRI was likewise performed using a 2D anatomical T₂-weighted spin-echo (SE) sequence with a 150 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution. This very first anatomical connectivity atlas of the male Japanese quail reveals 34 labeled fiber tracts and the existence of structural differences between the connectivity patterns characterizing the two lines. Thus, the link between the male Japanese quail’s connectivity and its underlying anatomical structures has reached a better understanding.

Tractography methods and findings in brain tumors and traumatic brain injury

White matter fiber tracking using diffusion magnetic resonance imaging (dMRI) provides a noninvasive approach to map brain connections, but improving anatomical accuracy has been a significant challenge since the birth of tractography methods. Utilizing tractography in brain studies therefore requires understanding of its technical limitations to avoid shortcomings and pitfalls. This review explores tractography limitations and how different white matter pathways pose different challenges to fiber tracking methodologies. We summarize the pros and cons of commonly-used methods, aiming to inform how tractography and its related analysis may lead to questionable results. Extending these experiences, we review the clinical utilization of tractography in patients with brain tumors and traumatic brain injury, starting from tensor-based tractography to more advanced methods. We discuss current limitations and highlight novel approaches in the context of these two conditions to inform future tractography developments.

A critical review of connectome validation studies

Diffusion MRI tractography is the most widely used macroscale method for mapping connectomes in vivo. However, tractography is prone to various errors and biases, and thus tractography-derived connectomes require careful validation. Here, we critically review studies that have developed or utilized phantoms and tracer maps to validate tractography-derived connectomes, either quantitatively or qualitatively. We identify key factors impacting connectome reconstruction accuracy, including streamline seeding, propagation and filtering methods, and consider the strengths and limitations of state-of-the-art connectome phantoms and associated validation studies. These studies demonstrate the inherent limitations of current fiber orientation models and tractography algorithms and their impact on connectome reconstruction accuracy. Reconstructing connectomes with both high sensitivity and high specificity is challenging, given that some tractography methods can generate an abundance of spurious connections, while others can overlook genuine fiber bundles. We argue that streamline filtering can minimize spurious connections and potentially improve the biological plausibility of connectomes derived from tractography. We find that algorithmic choices such as the tractography seeding methodology, angular threshold, and streamline propagation method can substantially impact connectome reconstruction accuracy. Hence, careful application of tractography is necessary to reconstruct accurate connectomes. Improvements in diffusion MRI acquisition techniques will not necessarily overcome current tractography limitations without accompanying modeling and algorithmic advances.

CHIASM, the human brain albinism and achiasma MRI dataset

We describe a collection of T1-, diffusion- and functional T2*-weighted magnetic resonance imaging data from human individuals with albinism and achiasma. This repository can be used as a test-bed to develop and validate tractography methods like diffusion-signal modeling and fiber tracking as well as to investigate the properties of the human visual system in individuals with congenital abnormalities. The MRI data is provided together with tools and files allowing for its preprocessing and analysis, along with the data derivatives such as manually curated masks and regions of interest for performing tractography.

Comparison of diffusion signal models for fiber tractography in eloquent glioma surgery - determination of accuracy under awake craniotomy conditions

OBJECTIVE

Fiber tractography(FT) has become an important non-invasive tool to ensure maximal safe tumor resection in eloquent glioma surgery. Intraoperatively applied FT is still predominantly based on Diffusion Tensor Imaging(DTI). However, reconstruction schemes of high angular resolution diffusion imaging(HARDI) data for high resolution fiber tractography(HRFT) are gaining increasing attention. The aim of this prospective study was to compare the accuracy of sophisticated HRFT-models compared with DTI-FT.

METHODS

Ten patients with eloquent gliomas underwent surgery under awake craniotomy conditions. The localization of acquisition points(AP), representing deteriorations during intraoperative electrostimulation(IOM) and neuropsychological mapping, were documented. The offsets of AP to the respective fiber bundle were calculated. Probabilistic QBI- and CSD-FT were compared to DTI-FT for the major language-associated fiber bundles (superior longitudinal fasciclus (SLF) II-IV, inferior fronto-occipital fasciculus (IFOF), inferior longitudinal fasciculus/medial longitudinal fasciculus (ILF/MLF).

RESULTS

Among 186 offset values, 46% were located closer than 10mm to the estimated fiber bundle (CSD:36%; DTI:40% and QBI:60%). Moreover, only 10 offsets were further away than 30mm (5%). Lowest mean min-offsets (SLF: 7.7±7.9mm; IFOF: 12.7±8.3mm; ILF/MLF: 17.7±6.7mm) were found for QBI, indicating a significant advantage compared with CSD or DTI (p<0.001), respectively. No significant differences were found between CSD-, and DTI-FT offsets (p=0.105), albeit for the compound SLF exclusively (p<0.001).

CONCLUSIONS

Comparing HRFT techniques QBI and CSD with DTI, QBI delivered significantly better results with lowest offsets and good correlation to IOM results. Besides, QBI-FT was feasible for neurosurgical pre- and intraoperative applications. Our findings suggest that a combined approach of QBI-FT and IOM under awake craniotomy is considerable for best preservation of neurological function in the presented setting. Overall, the implementation of selected HRFT models into neuronavigation systems seems to be a promising tool in glioma surgery.

Physical and digital phantoms for validating tractography and assessing artifacts

Fiber tractography is widely used to non-invasively map white-matter bundles in vivo using diffusion-weighted magnetic resonance imaging (dMRI). As it is the case for all scientific methods, proper validation is a key prerequisite for the successful application of fiber tractography, be it in the area of basic neuroscience or in a clinical setting. It is well-known that the indirect estimation of the fiber tracts from the local diffusion signal is highly ambiguous and extremely challenging. Furthermore, the validation of fiber tractography methods is hampered by the lack of a real ground truth, which is caused by the extremely complex brain microstructure that is not directly observable non-invasively and that is the basis of the huge network of long-range fiber connections in the brain that are the actual target of fiber tractography methods. As a substitute for in vivo data with a real ground truth that could be used for validation, a widely and successfully employed approach is the use of synthetic phantoms. In this work, we are providing an overview of the state-of-the-art in the area of physical and digital phantoms, answering the following guiding questions: "What are dMRI phantoms and what are they good for?", "What would the ideal phantom for validation fiber tractography look like?" and "What phantoms, phantom datasets and tools used for their creation are available to the research community?". We will further discuss the limitations and opportunities that come with the use of dMRI phantoms, and what future direction this field of research might take.

Tractography dissection variability: What happens when 42 groups dissect 14 white matter bundles on the same dataset?

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.

White matter microstructure alterations in cortico-striatal networks are associated with parkinsonism in schizophrenia spectrum disorders

The specific role of white matter (WM) microstructure in parkinsonism among patients with schizophrenia spectrum disorders (SSD) is largely unknown. To determine whether topographical alterations of WM microstructure contribute to parkinsonism in SSD patients, we examined healthy controls (HC, n=16) and SSD patients with and without parkinsonism, as defined by Simpson-Angus Scale total score of ≥4 (SSD-P, n=33) or <4 (SSD-nonP, n=62). We used whole brain tract-based spatial statistics (TBSS), tractometry (along tract statistics using TractSeg) and graph analytics (clustering coefficient (CCO), local betweenness centrality (BC)) to provide a framework of specific WM microstructural changes underlying parkinsonism in SSD. Using these methods, post hoc analyses showed (a) decreased fractional anisotrophy (FA), as measured via tractometry, in the corpus callosum, corticospinal tract and striato-fronto-orbital tract, and (b) increased CCO, as derived by graph analytics, in the left orbitofrontal cortex (OFC) and left superior frontal gyrus (SFG), in SSD-P patients when compared to SSD-nonP patients. Increased CCO in the left OFC and SFG was associated with SAS scores. These findings indicate the prominence of OFC alterations and aberrant connectivity with fronto-parietal regions and striatum in the pathogenesis of parkinsonism in SSD. This study further supports the notion of altered "bottom-up modulation" between basal ganglia and fronto-parietal regions in the pathobiology of parkinsonism, which may reflect an interaction between movement disorder intrinsic to SSD and antipsychotic drug-induced sensorimotor dysfunction.

Diffusion MRI tractography for neurosurgery: the basics, current state, technical reliability and challenges

Diffusion magnetic resonance imaging (dMRI) tractography is currently the only imaging technique that allows for non-invasive delineation and visualisation of white matter (WM) tractsin vivo,prompting rapid advances in related fields of brain MRI research in recent years. One of its major clinical applications is for pre-surgical planning and intraoperative image guidance in neurosurgery, where knowledge about the location of WM tracts nearby the surgical target can be helpful to guide surgical resection and optimise post-surgical outcomes. Surgical injuries to these WM tracts can lead to permanent neurological and functional deficits, making the accuracy of tractography reconstructions paramount. The quality of dMRI tractography is influenced by many modifiable factors, ranging from MRI data acquisition through to the post-processing of tractography output, with the potential of error propagation based on decisions made at each and subsequent processing steps. Research over the last 25 years has significantly improved the anatomical accuracy of tractography. An updated review about tractography methodology in the context of neurosurgery is now timely given the thriving research activities in dMRI, to ensure more appropriate applications in the clinical neurosurgical realm. This article aims to review the dMRI physics, and tractography methodologies, highlighting recent advances to provide the key concepts of tractography-informed neurosurgery, with a focus on the general considerations, the current state of practice, technical challenges, potential advances, and future demands to this field.

Trigeminal neuralgia diffusivities using Gaussian process classification and merged group tractography

ABSTRACT

Imaging of trigeminal neuralgia (TN) has demonstrated key diffusion tensor imaging-based diffusivity alterations in the trigeminal nerve; however, imaging has primarily focused on the peripheral nerve segment because of previous limitations in reliably segmenting small fiber bundles across multiple subjects. We used Selective Automated Group Integrated Tractography to study 36 subjects with TN (right-sided pain) and 36 sex-matched controls to examine the trigeminal nerve (fifth cranial nerve [CN V]), pontine decussation (TPT), and thalamocortical fibers (S1). Gaussian process classifiers were trained by scrolling a moving window over CN V, TPT, and S1 tractography centroids. Fractional anisotropy (FA), generalized FA, radial diffusivity, axial diffusivity, and mean diffusivity metrics were evaluated for both groups, analyzing TN vs control groups and affected vs unaffected sides. Classifiers that performed at greater-than-or-equal-to 70% accuracy were included. Gaussian process classifier consistently demonstrated bilateral trigeminal changes, differentiating them from controls with an accuracy of 80%. Affected and unaffected sides could be differentiated from each other with 75% accuracy. Bilateral TPT could be distinguished from controls with at least 85% accuracy. TPT left-right classification achieved 98% accuracy. Bilateral S1 could be differentiated from controls, where the affected S1 radial diffusivity classifier achieved 87% accuracy. This is the first TN study that combines group-wise merged tractography, machine learning classification, and analysis of the complete trigeminal pathways from the peripheral fibers to S1 cortex. This analysis demonstrates that TN is characterized by bilateral abnormalities throughout the trigeminal pathway compared with controls and abnormalities between affected and unaffected sides. This full pathway tractography study of TN demonstrates bilateral changes throughout the trigeminal pathway and changes between affected and unaffected sides.

Diffusion Kurtosis Imaging Fiber Tractography of Major White Matter Tracts in Neurosurgery

Diffusion tensor imaging (DTI)-based fiber tractography is routinely used in clinical applications to visualize major white matter tracts, such as the corticospinal tract (CST), optic radiation (OR), and arcuate fascicle (AF). Nevertheless, DTI is limited due to its capability of resolving intra-voxel multi-fiber populations. Sophisticated models often require long acquisition times not applicable in clinical practice. Diffusion kurtosis imaging (DKI), as an extension of DTI, combines sophisticated modeling of the diffusion process with short acquisition times but has rarely been investigated in fiber tractography. In this study, DTI- and DKI-based fiber tractography of the CST, OR, and AF was investigated in healthy volunteers and glioma patients. For the CST, significantly larger tract volumes were seen in DKI-based fiber tractography. Similar results were obtained for the OR, except for the right OR in patients. In the case of the AF, results of both models were comparable with DTI-based fiber tractography showing even significantly larger tract volumes in patients. In the case of the CST and OR, DKI-based fiber tractography contributes to advanced visualization under clinical time constraints, whereas for the AF, other models should be considered.

Multiparametric mapping of white matter microstructure in catatonia

Catatonia is characterized by motor, affective and behavioral abnormalities. To date, the specific role of white matter (WM) abnormalities in schizophrenia spectrum disorders (SSD) patients with catatonia is largely unknown. In this study, diffusion magnetic resonance imaging (dMRI) data were collected from 111 right-handed SSD patients and 28 healthy controls. Catatonic symptoms were examined on the Northoff Catatonia Rating Scale (NCRS). We used whole-brain tract-based spatial statistics (TBSS), tractometry (along tract statistics using TractSeg) and graph analytics (clustering coefficient-CCO, local betweenness centrality-BC) to provide a framework of specific WM microstructural abnormalities underlying catatonia in SSD. Following a categorical approach, post hoc analyses showed differences in fractional anisotrophy (FA) measured via tractometry in the corpus callosum, corticospinal tract and thalamo-premotor tract as well as increased CCO as derived by graph analytics of the right superior parietal cortex (SPC) and left caudate nucleus in catatonic patients (NCRS total score ≥ 3; n = 30) when compared to non-catatonic patients (NCRS total score = 0; n = 29). In catatonic patients according to DSM-IV-TR (n = 43), catatonic symptoms were associated with FA variations (tractometry) of the left corticospinal tract and CCO of the left orbitofrontal cortex, primary motor cortex, supplementary motor area and putamen. This study supports the notion that structural reorganization of WM bundles connecting orbitofrontal/parietal, thalamic and striatal regions contribute to catatonia in SSD patients.

A validation framework for neuroimaging software: The case of population receptive fields

Neuroimaging software methods are complex, making it a near certainty that some implementations will contain errors. Modern computational techniques (i.e., public code and data repositories, continuous integration, containerization) enable the reproducibility of the analyses and reduce coding errors, but they do not guarantee the scientific validity of the results. It is difficult, nay impossible, for researchers to check the accuracy of software by reading the source code; ground truth test datasets are needed. Computational reproducibility means providing software so that for the same input anyone obtains the same result, right or wrong. Computational validity means obtaining the right result for the ground-truth test data. We describe a framework for validating and sharing software implementations, and we illustrate its usage with an example application: population receptive field (pRF) methods for functional MRI data. The framework is composed of three main components implemented with containerization methods to guarantee computational reproducibility. In our example pRF application, those components are: (1) synthesis of fMRI time series from ground-truth pRF parameters, (2) implementation of four public pRF analysis tools and standardization of inputs and outputs, and (3) report creation to compare the results with the ground truth parameters. The framework was useful in identifying realistic conditions that lead to imperfect parameter recovery in all four pRF implementations, that would remain undetected using classic validation methods. We provide means to mitigate these problems in future experiments. A computational validation framework supports scientific rigor and creativity, as opposed to the oft-repeated suggestion that investigators rely upon a few agreed upon packages. We hope that the framework will be helpful to validate other critical neuroimaging algorithms, as having a validation framework helps (1) developers to build new software, (2) research scientists to verify the software's accuracy, and (3) reviewers to evaluate the methods used in publications and grants.

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.

Towards deep learning for connectome mapping: A block decomposition framework

We propose a new framework to map structural connectomes using deep learning and diffusion MRI. We show that our framework not only enables connectome mapping with a convolutional neural network (CNN), but can also be straightforwardly incorporated into conventional connectome mapping pipelines to enhance accuracy. Our framework involves decomposing the entire brain volume into overlapping blocks. Blocks are sufficiently small to ensure that a CNN can be efficiently trained to predict each block's internal connectivity architecture. We develop a block stitching algorithm to rebuild the full brain volume from these blocks and thereby map end-to-end connectivity matrices. To evaluate our block decomposition and stitching (BDS) framework independent of CNN performance, we first map each block's internal connectivity using conventional streamline tractography. Performance is evaluated using simulated diffusion MRI data generated from numerical connectome phantoms with known ground truth connectivity. Due to the redundancy achieved by allowing blocks to overlap, we find that our block decomposition and stitching steps per se can enhance the accuracy of probabilistic and deterministic tractography algorithms by up to 20-30%. Moreover, we demonstrate that our framework can improve the strength of structure-function coupling between in vivo diffusion and functional MRI data. We find that structural brain networks mapped with deep learning correlate more strongly with functional brain networks (r ​= ​0.45) than those mapped with conventional tractography (r ​= ​0.36). In conclusion, our BDS framework not only enables connectome mapping with deep learning, but its two constituent steps can be straightforwardly incorporated as part of conventional connectome mapping pipelines to enhance accuracy.

The role of diffusion tractography in refining glial tumor resection

Primary brain tumors are notoriously hard to resect surgically. Due to their infiltrative nature, finding the optimal resection boundary without damaging healthy tissue can be challenging. One potential tool to help make this decision is diffusion-weighted magnetic resonance imaging (dMRI) tractography. dMRI exploits the diffusion of water molecule along axons to generate a 3D modelization of the white matter bundles in the brain. This feature is particularly useful to visualize how a tumor affects its surrounding white matter and plan a surgical path. This paper reviews the different ways in which dMRI can be used to improve brain tumor resection, its benefits and also its limitations. We expose surgical tools that can be paired with dMRI to improve its impact on surgical outcome, such as loading the 3D tractography in the neuronavigation system and direct electrical stimulation to validate the position of the white matter bundles of interest. We also review articles validating dMRI findings using other anatomical investigation techniques, such as postmortem dissections, manganese-enhanced MRI, electrophysiological stimulations, and phantom studies with known ground truth. We will be discussing the areas of the brain where dMRI performs well and where the future challenges are. We will conclude this review with suggestions and take home messages for neurosurgeons, tractographers, and vendors for advancing the field and on how to benefit from tractography's use in clinical practice.

Presurgical simulation for neuroendoscopic procedures: Virtual study of the integrity of neurological pathways using diffusion tensor imaging tractography

Background:

White matter (WM) transgression is an unexplored concept in neuroendoscopy. Diffusion tensor image (DTI) tractography could be implemented as a planning and postoperative evaluation tool in functional disconnection procedures (FDPs), which are, currently, the subject of technological innovations. We intend to prove the usefulness of this planning method focused on the assessment of WM injury that is suitable for planning FDPs.

Methods:

Ten cranial magnetic resonance studies (20 sides) without pathological findings were processed. Fascicles were defined by two regions of interest (ROIs) using the fiber assignment method by the continuous tracking approach. Using three-dimensional (3D) simulation and DTI tractography, we created an 8-mm virtual endoscope and an uninjured inferior fronto-occipital fasciculus (IFOF) from two ROIs. The injured tract was generated using a third ROI built from the 3D model of the intersection of the oriented trajectory of the endoscope with the fascicle. Data and images were quantitatively and qualitatively analyzed.

Results:

The average percentage of the injured fibers was 32.0% (range: 12.4%–70%). The average intersected volume was 1.1 cm³ (range: 0.3–2.3 cm³). Qualitative analysis showed the inferior medial quadrant of the inferior fronto-occipital fasciculus (IFOF) as the most frequently injured region. No hemispherical asymmetry was found (P > 0.5).

Conclusion:

DTI tractography is a useful surgical planning tool that could be implemented in several endoscopic procedures. Together with a functional atlas, the presented technique provides a noninvasive method to assess the potential sequelae and thus to optimize the surgical route. The suggested method could be implemented to analyze pathological WM fascicles and to assess the surgical results of FDP such as hemispherotomy or amygdalohippocampectomy. More studies are needed to overcome the limitations of the tractography based information and to develop more anatomically and functionally reliable planning systems.

Diffusion MRI: Assessment of the Impact of Acquisition and Preprocessing Methods Using the BrainVISA-Diffuse Toolbox

Diffusion MR images are prone to severe geometric distortions induced by head movement, eddy-current and inhomogeneity of magnetic susceptibility. Various correction methods have been proposed that depend on the choice of the acquisition settings and potentially provide highly different data quality. However, the impact of this choice has not been evaluated in terms of the ratio between scan time and preprocessed data quality. This study aims at investigating the impact of six well-known preprocessing methods, each associated to specific acquisition settings, on the outcome of diffusion analyses. For this purpose, we developed a comprehensive toolbox called Diffuse which automatically guides the user to the best preprocessing pipeline according to the input data. Using MR images of 20 subjects from the HCP dataset, we compared the six pre-processing pipelines regarding the following criteria: the ability to recover brain’s true geometry, the tensor model estimation and derived indices in the white matter, and finally the spatial dispersion of six well known connectivity pathways. As expected the pipeline associated to the longer acquisition fully repeated with reversed phase-encoding (RPE) yielded the higher data quality and was used as a reference to evaluate the other pipelines. In this way, we highlighted several significant aspects of other pre-processing pipelines. Our results first established that eddy-current correction improves the tensor-fitting performance with a localized impact especially in the corpus callosum. Concerning susceptibility distortions, we showed that the use of a field map is not sufficient and involves additional smoothing, yielding to an artificial decrease of tensor-fitting error. Of most importance, our findings demonstrate that, for an equivalent scan time, the acquisition of a b0 volume with RPE ensures a better brain’s geometry reconstruction and local improvement of tensor quality, without any smoothing of the image. This was found to be the best scan time/data quality compromise. To conclude, this study highlights and attempts to quantify the strong dependence of diffusion metrics on acquisition settings and preprocessing methods.

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.

The open diffusion data derivatives, brain data upcycling via integrated publishing of derivatives and reproducible open cloud services

We describe the Open Diffusion Data Derivatives (O3D) repository: an integrated collection of preserved brain data derivatives and processing pipelines, published together using a single digital-object-identifier. The data derivatives were generated using modern diffusion-weighted magnetic resonance imaging data (dMRI) with diverse properties of resolution and signal-to-noise ratio. In addition to the data, we publish all processing pipelines (also referred to as open cloud services). The pipelines utilize modern methods for neuroimaging data processing (diffusion-signal modelling, fiber tracking, tractography evaluation, white matter segmentation, and structural connectome construction). The O3D open services can allow cognitive and clinical neuroscientists to run the connectome mapping algorithms on new, user-uploaded, data. Open source code implementing all O3D services is also provided to allow computational and computer scientists to reuse and extend the processing methods. Publishing both data-derivatives and integrated processing pipeline promotes practices for scientific reproducibility and data upcycling by providing open access to the research assets for utilization by multiple scientific communities.

Ex vivo diffusion MRI of the human brain: Technical challenges and recent advances

This review discusses ex vivo diffusion magnetic resonance imaging (dMRI) as an important research tool for neuroanatomical investigations and the validation of in vivo dMRI techniques, with a focus on the human brain. We review the challenges posed by the properties of post-mortem tissue, and discuss state-of-the-art tissue preparation methods and recent advances in pulse sequences and acquisition techniques to tackle these. We then review recent ex vivo dMRI studies of the human brain, highlighting the validation of white matter orientation estimates and the atlasing and mapping of large subcortical structures. We also give particular emphasis to the delineation of layered gray matter structure with ex vivo dMRI, as this application illustrates the strength of its mesoscale resolution over large fields of view. We end with a discussion and outlook on future and potential directions of the field.

Mapping connectomes with diffusion MRI: deterministic or probabilistic tractography?

PURPOSE

Human connectomics necessitates high-throughput, whole-brain reconstruction of multiple white matter fiber bundles. Scaling up tractography to meet these high-throughput demands yields new fiber tracking challenges, such as minimizing spurious connections and controlling for gyral biases. The aim of this study is to determine which of the two broadest classes of tractography algorithms-deterministic or probabilistic-is most suited to mapping connectomes.

METHODS

This study develops numerical connectome phantoms that feature realistic network topologies and that are matched to the fiber complexity of in vivo diffusion MRI (dMRI) data. The phantoms are utilized to evaluate the performance of tensor-based and multi-fiber implementations of deterministic and probabilistic tractography.

RESULTS

For connectome phantoms that are representative of the fiber complexity of in vivo dMRI, multi-fiber deterministic tractography yields the most accurate connectome reconstructions (F-measure = 0.35). Probabilistic algorithms are hampered by an abundance of false-positive connections, leading to lower specificity (F = 0.19). While omitting connections with the fewest number of streamlines (thresholding) improves the performance of probabilistic algorithms (F = 0.38), multi-fiber deterministic tractography remains optimal when it benefits from thresholding (F = 0.42).

CONCLUSIONS

Multi-fiber deterministic tractography is well suited to connectome mapping, while connectome thresholding is essential when using probabilistic algorithms.

When tractography meets tracer injections: a systematic study of trends and variation sources of diffusion-based connectivity

Tractography is a powerful technique capable of non-invasively reconstructing the structural connections in the brain using diffusion MRI images, but the validation of tractograms is challenging due to lack of ground truth. Owing to recent developments in mapping the mouse brain connectome, high-resolution tracer injection-based axonal projection maps have been created and quickly adopted for the validation of tractography. Previous studies using tracer injections mainly focused on investigating the match in projections and optimal tractography protocols. Being a complicated technique, however, tractography relies on multiple stages of operations and parameters. These factors introduce large variabilities in tractograms, hindering the optimization of protocols and making the interpretation of results difficult. Based on this observation, in contrast to previous studies, in this work we focused on quantifying and ranking the amount of performance variation introduced by these factors. For this purpose, we performed over a million tractography experiments and studied the variability across different subjects, injections, anatomical constraints and tractography parameters. By using N-way ANOVA analysis, we show that all tractography parameters are significant and importantly performance variations with respect to the differences in subjects are comparable to the variations due to tractography parameters, which strongly underlines the importance of fully documenting the tractography protocols in scientific experiments. We also quantitatively show that inclusion of anatomical constraints is the most significant factor for improving tractography performance. Although this critical factor helps reduce false positives, our analysis indicates that anatomy-informed tractography still fails to capture a large portion of axonal projections.

On the Origin of Individual Functional Connectivity Variability: The Role of White Matter Architecture

Fingerprint patterns derived from functional connectivity (FC) can be used to identify subjects across groups and sessions, indicating that the topology of the brain substantially differs between individuals. However, the source of FC variability inferred from resting-state functional magnetic resonance imaging remains unclear. One possibility is that these variations are related to individual differences in white matter structural connectivity (SC). However, directly comparing FC with SC is challenging given the many potential biases associated with quantifying their respective strengths. In an attempt to circumvent this, we employed a recently proposed test-retest approach that better quantifies inter-subject variability by first correcting for intra-subject nuisance variability (i.e., head motion, physiological differences in brain state, etc.) that can artificially influence FC and SC measures. Therefore, rather than directly comparing the strength of FC with SC, we asked whether brain regions with, for example, low inter-subject FC variability also exhibited low SC variability. From this, we report two main findings: First, at the whole-brain level, SC variability was significantly lower than FC variability, indicating that an individual's structural connectome is far more similar to another relative to their functional counterpart even after correcting for noise. Second, although FC and SC variability were mutually low in some brain areas (e.g., primary somatosensory cortex) and high in others (e.g., memory and language areas), the two were not significantly correlated across all cortical and sub-cortical regions. Taken together, these results indicate that even after correcting for factors that may differently affect FC and SC, the two, nonetheless, remain largely independent of one another. Further work is needed to understand the role that direct anatomical pathways play in supporting vascular-based measures of FC and to what extent these measures are dictated by anatomical connectivity.

Probing the reproducibility of quantitative estimates of structural connectivity derived from global tractography

As quantitative measures derived from fiber tractography are increasingly being used to characterize the structural connectivity of the brain, it is important to establish their reproducibility. However, no such information is as yet available for global tractography. Here we provide the first comprehensive analysis of the reproducibility of streamline counts derived from global tractography as quantitative estimates of structural connectivity. In a sample of healthy young adults scanned twice within one week, within-session and between-session test-retest reproducibility was estimated for streamline counts of connections based on regions of the AAL atlas using the intraclass correlation coefficient (ICC) for absolute agreement. We further evaluated the influence of the type of head-coil (12 versus 32 channels) and the number of reconstruction repetitions (reconstructing streamlines once or aggregated over ten repetitions). Factorial analyses demonstrated that reproducibility was significantly greater for within- than between-session reproducibility and significantly increased by aggregating streamline counts over ten reconstruction repetitions. Using a high-resolution head-coil incurred only small beneficial effects. Overall, ICC values were positively correlated with the streamline count of a connection. Additional analyses assessed the influence of different selection variants (defining fuzzy versus no fuzzy borders of the seed mask; selecting streamlines that end in versus pass through a seed) showing that an endpoint-based variant using fuzzy selection provides the best compromise between reproducibility and anatomical specificity. In sum, aggregating quantitative indices over repeated estimations and higher numbers of streamlines are important determinants of test-retest reproducibility. If these factors are taken into account, streamline counts derived from global tractography provide an adequately reproducible quantitative measure that can be used to gauge the structural connectivity of the brain in health and disease.

Network connectivity determines cortical thinning in early Parkinson's disease progression

Here we test the hypothesis that the neurodegenerative process in Parkinson’s disease (PD) moves stereotypically along neural networks, possibly reflecting the spread of toxic alpha-synuclein molecules. PD patients (n = 105) and matched controls (n = 57) underwent T1-MRI at entry and 1 year later as part of the Parkinson’s Progression Markers Initiative. Over this period, PD patients demonstrate significantly greater cortical thinning than controls in parts of the left occipital and bilateral frontal lobes and right somatomotor-sensory cortex. Cortical thinning is correlated to connectivity (measured functionally or structurally) to a “disease reservoir” evaluated by MRI at baseline. The atrophy pattern in the ventral frontal lobes resembles one described in certain cases of Alzheimer’s disease. Our findings suggest that disease propagation to the cortex in PD follows neuronal connectivity and that disease spread to the cortex may herald the onset of cognitive impairment.

In Parkinson’s disease (PD), neurodegeneration spreads from the brainstem to the cerebral cortex. Here, in a longitudinal study of PD patients, the authors found that cortical thinning followed neural connectivity from a “disease reservoir”.

The challenge of mapping the human connectome based on diffusion tractography

Tractography based on non-invasive diffusion imaging is central to the study of human brain connectivity. To date, the approach has not been systematically validated in ground truth studies. Based on a simulated human brain data set with ground truth tracts, we organized an open international tractography challenge, which resulted in 96 distinct submissions from 20 research groups. Here, we report the encouraging finding that most state-of-the-art algorithms produce tractograms containing 90% of the ground truth bundles (to at least some extent). However, the same tractograms contain many more invalid than valid bundles, and half of these invalid bundles occur systematically across research groups. Taken together, our results demonstrate and confirm fundamental ambiguities inherent in tract reconstruction based on orientation information alone, which need to be considered when interpreting tractography and connectivity results. Our approach provides a novel framework for estimating reliability of tractography and encourages innovation to address its current limitations.

Though tractography is widely used, it has not been systematically validated. Here, authors report results from 20 groups showing that many tractography algorithms produce both valid and invalid bundles.

Fiber tractography using machine learning

We present a fiber tractography approach based on a random forest classification and voting process, guiding each step of the streamline progression by directly processing raw diffusion-weighted signal intensities. For comparison to the state-of-the-art, i.e. tractography pipelines that rely on mathematical modeling, we performed a quantitative and qualitative evaluation with multiple phantom and in vivo experiments, including a comparison to the 96 submissions of the ISMRM tractography challenge 2015. The results demonstrate the vast potential of machine learning for fiber tractography.

Performance of unscented Kalman filter tractography in edema: Analysis of the two-tensor model

Diffusion MRI tractography is increasingly used in pre-operative neurosurgical planning to visualize critical fiber tracts. However, a major challenge for conventional tractography, especially in patients with brain tumors, is tracing fiber tracts that are affected by vasogenic edema, which increases water content in the tissue and lowers diffusion anisotropy. One strategy for improving fiber tracking is to use a tractography method that is more sensitive than the traditional single-tensor streamline tractography.

We performed experiments to assess the performance of two-tensor unscented Kalman filter (UKF) tractography in edema. UKF tractography fits a diffusion model to the data during fiber tracking, taking advantage of prior information from the previous step along the fiber. We studied UKF performance in a synthetic diffusion MRI digital phantom with simulated edema and in retrospective data from two neurosurgical patients with edema affecting the arcuate fasciculus and corticospinal tracts. We compared the performance of several tractography methods including traditional streamline, UKF single-tensor, and UKF two-tensor. To provide practical guidance on how the UKF method could be employed, we evaluated the impact of using various seed regions both inside and outside the edematous regions, as well as the impact of parameter settings on the tractography sensitivity. We quantified the sensitivity of different methods by measuring the percentage of the patient-specific fMRI activation that was reached by the tractography.

We expected that diffusion anisotropy threshold parameters, as well as the inclusion of a free water model, would significantly influence the reconstruction of edematous WM fiber tracts, because edema increases water content in the tissue and lowers anisotropy. Contrary to our initial expectations, varying the fractional anisotropy threshold and including a free water model did not affect the UKF two-tensor tractography output appreciably in these two patient datasets. The most effective parameter for increasing tracking sensitivity was the generalized anisotropy (GA) threshold, which increased the length of tracked fibers when reduced to 0.075. In addition, the most effective seeding strategy was seeding in the whole brain or in a large region outside of the edema.

Overall, the main contribution of this study is to provide insight into how UKF tractography can work, using a two-tensor model, to begin to address the challenge of fiber tract reconstruction in edematous regions near brain tumors.

•

Reconstruction of edematous white matter from diffusion MRI is investigated.

•

The performance of two–tensor unscented Kalman filter (UKF) tractography is assessed.

•

The two–tensor model in UKF is analyzed in phantom and patient data experiments.

•

Practical guidance on employing the UKF method in neurosurgical patients is provided

Structural connectivity differences in motor network between tremor-dominant and nontremor Parkinson's disease

Motor phenotypes of Parkinson's disease (PD) are recognized to have different prognosis and therapeutic response, but the neural basis for this clinical heterogeneity remains largely unknown. The main aim of this study was to compare differences in structural connectivity metrics of the main motor network between tremor-dominant and nontremor PD phenotypes (TD-PD and NT-PD, respectively) using probabilistic tractography-based network analysis. A total of 63 PD patients (35 TD-PD patients and 28 NT-PD patients) and 30 healthy controls underwent a 3 T MRI. Next, probabilistic tractography-based network analysis was performed to assess structural connectivity in cerebello-thalamo-basal ganglia-cortical circuits, by measuring the connectivity indices of each tract and the efficiency of each node. Furthermore, dopamine transporter single-photon emission computed tomography (DAT-SPECT) with ¹²³ I-ioflupane was used to assess dopaminergic striatal depletion in all PD patients. Both PD phenotypes showed nodal abnormalities in the substantia nigra, in agreement with DAT-SPECT evaluation. In addition, NT-PD patients displayed connectivity alterations in nigro-pallidal and fronto-striatal pathways, compared with both controls and TD-PD patients, in which the same motor connections seemed to be relatively spared. Of note, in NT-PD group, rigidity-bradykinesia score correlated with fronto-striatal connectivity abnormalities. These findings demonstrate that structural connectivity alterations occur in the cortico-basal ganglia circuit of NT-PD patients, but not in TD-PD patients, suggesting that these anatomical differences may underlie different motor phenotypes of PD. Hum Brain Mapp 38:4716-4729, 2017. © 2017 Wiley Periodicals, Inc.

White matter microstructure variations contribute to neurological soft signs in healthy adults

OBJECTIVE

Neurological soft signs (NSS) are core features of psychiatric disorders with significant neurodevelopmental origin. However, it is unclear whether NSS correlates are associated with neuropathological processes underlying the disease or if they are confounded by medication. Given that NSS are also present in healthy persons (HP), investigating HP could reveal NSS correlates, which are not biased by disease-specific processes or drug treatment. Therefore, we used a combination of diffusion MRI analysis tools to provide a framework of specific white matter (WM) microstructure variations underlying NSS in HP.

METHOD

NSS of 59 HP were examined on the Heidelberg Scale and related to diffusion associated metrics. Using tract-based spatial statistics (TBSS), we studied WM variations in fractional anisotropy (FA) as well as radial (RD), axial (AD), and mean diffusivity (MD). Using graph analytics (clustering coefficient-CC, local betweenness centrality -BC), we then explored DTI-derived structural network variations in regions identified by previous MRI studies on NSS.

RESULTS

NSS scores were negatively associated with RD, AD and MD in corpus callosum, brainstem and cerebellum (P < 0.05, corr.). NSS scores were negatively associated with CC and BC of the pallidum, the superior parietal gyrus, the precentral sulcus, the insula, and the cingulate gyrus (P < 0.05, uncorr.).

CONCLUSION

The present study supports the notion that WM microstructure variations in subcortical and cortical sensorimotor regions contribute to NSS expression in young HP. Hum Brain Mapp 38:3552-3565, 2017. © 2017 Wiley Periodicals, Inc.

Active delineation of Meyer's loop using oriented priors through MAGNEtic tractography (MAGNET)

Streamline tractography algorithms infer connectivity from diffusion MRI (dMRI) by following diffusion directions which are similarly aligned between neighboring voxels. However, not all white matter (WM) fascicles are organized in this manner. For example, Meyer's loop is a highly curved portion of the optic radiation (OR) that exhibits a narrow turn, kissing and crossing pathways, and changes in fascicle dispersion. From a neurosurgical perspective, damage to Meyer's loop carries a potential risk of inducing vision deficits to the patient, especially during temporal lobe resection surgery. To prevent such impairment, achieving an accurate delineation of Meyer's loop with tractography is thus of utmost importance. However, current algorithms tend to under-estimate the full extent of Meyer's loop, mainly attributed to the aforementioned rule for connectivity which requires a direction to be chosen across a field of orientations. In this article, it was demonstrated that MAGNEtic Tractography (MAGNET) can benefit Meyer's loop delineation by incorporating anatomical knowledge of the expected fiber orientation to overcome local ambiguities. A new ROI-mechanism was proposed which supplies additional information to streamline reconstruction algorithms by the means of oriented priors. Their results showed that MAGNET can accurately generate Meyer's loop in all of our 15 child subjects (8 males; mean age 10.2 years ± 3.1). It effectively improved streamline coverage when compared with deterministic tractography, and significantly reduced the distance between the anterior-most portion of Meyer's loop and the temporal pole by 16.7 mm on average, a crucial landmark used for preoperative planning of temporal lobe surgery. Hum Brain Mapp 38:509-527, 2017. © 2016 Wiley Periodicals, Inc.

Assessment of White Matter Transgression During Neuroendoscopic Procedures Using Diffusion Tensor Image Fiber Tracking

BACKGROUND

Presurgical planning allows anticipating intraoperative difficulties, increasing efficiency, and reducing risks. Neuroendoscopy is a minimally invasive technique whose related complications have been focused on cortical function and surface vessels injury. However, white matter disruption has been insufficiently acknowledged.

OBJECTIVE

To present a new surgical planning method based on diffusion tensor image that allows quantifying subcortical transgression and optimizing neuroendoscopic trajectories.

METHODS

Ten cranial magnetic resonance studies (20 sides) without pathologic findings were anonymized and processed. A standard transcortical approach to the frontal horn was used to study the transgression of the corpus callosum (CC) and cingulum (Ci) caused by a virtual endoscope (VE) oriented from the Kocher point to the foramen of Monro. An 8-mm VE model was created, oriented, and coregistered. VE-CC and VE-Ci intersections were segmented. The number and volume of injured fibers were measured, intersections were quantified, and the percentage of tract transgression was calculated. The areas damaged by the VE were also recorded.

RESULTS

Among the CC fibers, 16.4% were injured (range: 3.3%-37%) and 26.7% of fibers on Ci (rank: 0%-73.4%). The average intersected volumes were 19.1% (range: 4.2%-53.2%) for CC and 33.2% for Ci (range: 0%-73.7%). Qualitative analysis showed the lateral aspect of both tracts as the most frequently injured region. No hemispherical asymmetry was found (P > 0.05).

CONCLUSION

This method using tractography and oriented models of surgical instruments allows assessing white matter transgression, both qualitatively and quantitatively, for a deep brain trajectory. Thus our method permits surgeons to optimize safety and avoid transgression of eloquent tracts during surgical planning. Nevertheless, more studies are necessary.

Automated white matter fiber tract identification in patients with brain tumors

We propose a method for the automated identification of key white matter fiber tracts for neurosurgical planning, and we apply the method in a retrospective study of 18 consecutive neurosurgical patients with brain tumors. Our method is designed to be relatively robust to challenges in neurosurgical tractography, which include peritumoral edema, displacement, and mass effect caused by mass lesions. The proposed method has two parts. First, we learn a data-driven white matter parcellation or fiber cluster atlas using groupwise registration and spectral clustering of multi-fiber tractography from healthy controls. Key fiber tract clusters are identified in the atlas. Next, patient-specific fiber tracts are automatically identified using tractography-based registration to the atlas and spectral embedding of patient tractography.

Results indicate good generalization of the data-driven atlas to patients: 80% of the 800 fiber clusters were identified in all 18 patients, and 94% of the 800 fiber clusters were found in 16 or more of the 18 patients. Automated subject-specific tract identification was evaluated by quantitative comparison to subject-specific motor and language functional MRI, focusing on the arcuate fasciculus (language) and corticospinal tracts (motor), which were identified in all patients. Results indicate good colocalization: 89 of 95, or 94%, of patient-specific language and motor activations were intersected by the corresponding identified tract. All patient-specific activations were within 3mm of the corresponding language or motor tract. Overall, our results indicate the potential of an automated method for identifying fiber tracts of interest for neurosurgical planning, even in patients with mass lesions.

•

Spectral clustering machine learning approach for white matter tract identification

•

Data-driven white matter parcellation learned from healthy subjects tractography

•

White matter parcellation applied to 18 consecutive patients with brain tumors

•

Arcuate fasciculus and corticospinal tracts identified in all patients

•

All tracts within 3 mm of corresponding patient-specific functional activations

Merged Group Tractography Evaluation with Selective Automated Group Integrated Tractography

Introduction: Tractography analysis in group-based studies across large populations has been difficult to implement. We propose Selective Automated Group Integrated Tractography (SAGIT), an automated group tractography software platform that incorporates multiple diffusion magnetic resonance imaging (dMRI) practices which will allow great accessibility to group-wise dMRI. We use a merged tractography approach that permits evaluation of tractography datasets at the group level. We also introduce an image normalized overlap score (NOS) that measures the quality of the group tractography results. We deploy SAGIT to evaluate deterministic and probabilistic constrained spherical deconvolution (CST det , CST prob ) tractography, eXtended Streamline Tractography (XST), and diffusion tensor tractography (DTT) in their ability to delineate different neuroanatomy, as well as validating NOS across these different brain regions. Materials and methods: Magnetic resonance sequences were acquired from 42 healthy adults. Anatomical and group registrations were performed using Automated Normalization Tools. Cortical segmentation was performed using FreeSurfer. Four tractography algorithms were used to delineate six sets of neuroanatomy: fornix, facial/vestibular-cochlear cranial nerve complex, vagus nerve, rubral-cerebellar decussation, optic radiation, and auditory radiation. The tracts were generated both with and without region of interest filters. The generated visual reports were then evaluated by five neuroscientists. Results: At a group level, merged tractography demonstrated that different methods have different fiber distribution characteristics. CST prob is prone to false-positives, and thereby suitable in anatomy with strong priors. CST det and XST are more conservative, but have greater difficulty resolving hemispherical decussation and distant crossing projections. DTT consistently shows the worst reproducibility across the anatomies. Linear regression of rater scores against NOS shows significant (p < 0.05) correlation of the two sets of scores in filtered tractography. However, correlations are not significant (p > 0.05) for unfiltered tractography. Conclusion: The tractography results demonstrated reliable and consistent performance of SAGIT across multiple subjects and techniques. Through SAGIT, we quantifiably demonstrated that different algorithms showed different strengths and weaknesses at a group level. While no single algorithm seems to be suitable for all anatomical tasks, it is useful to consider the use of a mix of algorithms for different anatomical segments. SAGIT appears to be a promising group-wise tractography analysis approach for this purpose.