Anatomical accuracy of standard-practice tractography algorithms in the motor system - A histological validation in the squirrel monkey brain

Advancements in Diffusion MRI Tractography for Neurosurgery

ABSTRACT

Diffusion magnetic resonance imaging tractography is a noninvasive technique that enables the visualization and quantification of white matter tracts within the brain. It is extensively used in preoperative planning for brain tumors, epilepsy, and functional neurosurgical procedures such as deep brain stimulation. Over the past 25 years, significant advancements have been made in imaging acquisition, fiber direction estimation, and tracking methods, resulting in considerable improvements in tractography accuracy. The technique enables the mapping of functionally critical pathways around surgical sites to avoid permanent functional disability. When the limitations are adequately acknowledged and considered, tractography can serve as a valuable tool to safeguard critical white matter tracts and provides insight regarding changes in normal white matter and structural connectivity of the whole brain beyond local lesions. In functional neurosurgical procedures such as deep brain stimulation, it plays a significant role in optimizing stimulation sites and parameters to maximize therapeutic efficacy and can be used as a direct target for therapy. These insights can aid in patient risk stratification and prognosis. This article aims to discuss state-of-the-art tractography methodologies and their applications in preoperative planning and highlight the challenges and new prospects for the use of tractography in daily clinical practice.

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.

The spatial extent of anatomical connections within the thalamus varies across the cortical hierarchy in humans and macaques

Each cortical area has a distinct pattern of anatomical connections within the thalamus, a central subcortical structure composed of functionally and structurally distinct nuclei. Previous studies have suggested that certain cortical areas may have more extensive anatomical connections that target multiple thalamic nuclei, which potentially allows them to modulate distributed information flow. However, there is a lack of quantitative investigations into anatomical connectivity patterns within the thalamus. Consequently, it remains unknown if cortical areas exhibit systematic differences in the extent of their anatomical connections within the thalamus. To address this knowledge gap, we used diffusion magnetic resonance imaging (dMRI) to perform brain-wide probabilistic tractography for 828 healthy adults from the Human Connectome Project. We then developed a framework to quantify the spatial extent of each cortical area's anatomical connections within the thalamus. Additionally, we leveraged resting-state functional MRI, cortical myelin, and human neural gene expression data to test if the extent of anatomical connections within the thalamus varied along the cortical hierarchy. Our results revealed two distinct corticothalamic tractography motifs: 1) a sensorimotor cortical motif characterized by focal thalamic connections targeting posterolateral thalamus, associated with fast, feed-forward information flow; and 2) an associative cortical motif characterized by diffuse thalamic connections targeting anteromedial thalamus, associated with slow, feed-back information flow. These findings were consistent across human subjects and were also observed in macaques, indicating cross-species generalizability. Overall, our study demonstrates that sensorimotor and association cortical areas exhibit differences in the spatial extent of their anatomical connections within the thalamus, which may support functionally-distinct cortico-thalamic information flow.

Self-supervised segmentation and characterization of fiber bundles in anatomic tracing data

Anatomic tracing is the gold standard tool for delineating brain connections and for validating more recently developed imaging approaches such as diffusion MRI tractography. A key step in the analysis of data from tracer experiments is the careful, manual charting of fiber trajectories on histological sections. This is a very time-consuming process, which limits the amount of annotated tracer data that are available for validation studies. Thus, there is a need to accelerate this process by developing a method for computer-assisted segmentation. Such a method must be robust to the common artifacts in tracer data, including variations in the intensity of stained axons and background, as well as spatial distortions introduced by sectioning and mounting the tissue. The method should also achieve satisfactory performance using limited manually charted data for training. Here we propose the first deeplearning method, with a self-supervised loss function, for segmentation of fiber bundles on histological sections from macaque brains that have received tracer injections. We address the limited availability of manual labels with a semi-supervised training technique that takes advantage of unlabeled data to improve performance. We also introduce anatomic and across-section continuity constraints to improve accuracy. We show that our method can be trained on manually charted sections from a single case and segment unseen sections from different cases, with a true positive rate of ~0.80. We further demonstrate the utility of our method by quantifying the density of fiber bundles as they travel through different white-matter pathways. We show that fiber bundles originating in the same injection site have different levels of density when they travel through different pathways, a finding that can have implications for microstructure-informed tractography methods. The code for our method is available at https://github.com/v-sundaresan/fiberbundle_seg_tracing.

Tractography passes the test: Results from the diffusion-simulated connectivity (disco) challenge

Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods.

Using light and X-ray scattering to untangle complex neuronal orientations and validate diffusion MRI

Disentangling human brain connectivity requires an accurate description of nerve fiber trajectories, unveiled via detailed mapping of axonal orientations. However, this is challenging because axons can cross one another on a micrometer scale. Diffusion magnetic resonance imaging (dMRI) can be used to infer axonal connectivity because it is sensitive to axonal alignment, but it has limited spatial resolution and specificity. Scattered light imaging (SLI) and small-angle X-ray scattering (SAXS) reveal axonal orientations with microscopic resolution and high specificity, respectively. Here, we apply both scattering techniques on the same samples and cross-validate them, laying the groundwork for ground-truth axonal orientation imaging and validating dMRI. We evaluate brain regions that include unidirectional and crossing fibers in human and vervet monkey brain sections. SLI and SAXS quantitatively agree regarding in-plane fiber orientations including crossings, while dMRI agrees in the majority of voxels with small discrepancies. We further use SAXS and dMRI to confirm theoretical predictions regarding SLI determination of through-plane fiber orientations. Scattered light and X-ray imaging can provide quantitative micrometer 3D fiber orientations with high resolution and specificity, facilitating detailed investigations of complex fiber architecture in the animal and human brain.

In vivo probabilistic atlas of white matter tracts of the human subthalamic area combining track density imaging and optimized diffusion tractography

The human subthalamic area is a region of high anatomical complexity, tightly packed with tiny fiber bundles. Some of them, including the pallidothalamic, cerebello-thalamic, and mammillothalamic tracts, are relevant targets in functional neurosurgery for various brain diseases. Diffusion-weighted imaging-based tractography has been suggested as a useful tool to map white matter pathways in the human brain in vivo and non-invasively, though the reconstruction of these specific fiber bundles is challenging due to their small dimensions and complex anatomy. To the best of our knowledge, a population-based, in vivo probabilistic atlas of subthalamic white matter tracts is still missing. In the present work, we devised an optimized tractography protocol for reproducible reconstruction of the tracts of subthalamic area in a large data sample from the Human Connectome Project repository. First, we leveraged the super-resolution properties and high anatomical detail provided by short tracks track-density imaging (stTDI) to identify the white matter bundles of the subthalamic area on a group-level template. Tracts identification on the stTDI template was also aided by visualization of histological sections of human specimens. Then, we employed this anatomical information to drive tractography at the subject-level, optimizing tracking parameters to maximize between-subject and within-subject similarities as well as anatomical accuracy. Finally, we gathered subject level tracts reconstructed with optimized tractography into a large-scale, normative population atlas. We suggest that this atlas could be useful in both clinical anatomy and functional neurosurgery settings, to improve our understanding of the complex morphology of this important brain region.

Insights from the IronTract challenge: Optimal methods for mapping brain pathways from multi-shell diffusion MRI

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

Evaluating reproducibility and subject-specificity of microstructure-informed connectivity

Tractography enables identifying and evaluating the healthy and diseased brain's white matter pathways from diffusion-weighted magnetic resonance imaging data. As previous evaluation studies have reported significant false-positive estimation biases, recent microstructure-informed tractography algorithms have been introduced to improve the trade-off between specificity and sensitivity. However, a major limitation for characterizing the performance of these techniques is the lack of ground truth brain data. In this study, we compared the performance of two relevant microstructure-informed tractography methods, SIFT2 and COMMIT, by assessing the subject specificity and reproducibility of their derived white matter pathways. Specifically, twenty healthy young subjects were scanned at eight different time points at two different sites. Subject specificity and reproducibility were evaluated using the whole-brain connectomes and a subset of 29 white matter bundles. Our results indicate that although the raw tractograms are more vulnerable to the presence of false-positive connections, they are highly reproducible, suggesting that the estimation bias is subject-specific. This high reproducibility was preserved when microstructure-informed tractography algorithms were used to filter the raw tractograms. Moreover, the resulting track-density images depicted a more uniform coverage of streamlines throughout the white matter, suggesting that these techniques could increase the biological meaning of the estimated fascicles. Notably, we observed an increased subject specificity by employing connectivity pre-processing techniques to reduce the underlaying noise and the data dimensionality (using principal component analysis), highlighting the importance of these tools for future studies. Finally, no strong bias from the scanner site or time between measurements was found. The largest intraindividual variance originated from the sole repetition of data measurements (inter-run).

The Case for Optimized Edge-Centric Tractography at Scale

The anatomic validity of structural connectomes remains a significant uncertainty in neuroimaging. Edge-centric tractography reconstructs streamlines in bundles between each pair of cortical or subcortical regions. Although edge bundles provides a stronger anatomic embedding than traditional connectomes, calculating them for each region-pair requires exponentially greater computation. We observe that major speedup can be achieved by reducing the number of streamlines used by probabilistic tractography algorithms. To ensure this does not degrade connectome quality, we calculate the identifiability of edge-centric connectomes between test and re-test sessions as a proxy for information content. We find that running PROBTRACKX2 with as few as 1 streamline per voxel per region-pair has no significant impact on identifiability. Variation in identifiability caused by streamline count is overshadowed by variation due to subject demographics. This finding even holds true in an entirely different tractography algorithm using MRTrix. Incidentally, we observe that Jaccard similarity is more effective than Pearson correlation in calculating identifiability for our subject population.

Investigating the Anatomy and Microstructure of the Dentato-rubro-thalamic and Subthalamo-ponto-cerebellar Tracts in Parkinson's Disease

Cerebellar-thalamic connections play a central role in deep brain stimulation-based treatment of tremor syndromes. Here, we used diffusion Magnetic Resonance Imaging (MRI) tractography to delineate the main cerebellar peduncles as well as two main white matter tracts that connect the cerebellum with the thalamus, the dentato-rubro-thalamic tract (DRTT) and the subthalamo-ponto-cerebellar tract (SPCT). We first developed a reconstruction protocol in young healthy adults with high-resolution diffusion imaging data and then demonstrate feasibility of transferring this protocol to clinical studies using standard diffusion MRI data from a cohort of patients with Parkinson's disease (PD) and their matched healthy controls. The tracts obtained closely corresponded to the previously described anatomical pathways and features of the DRTT and the SPCT. Second, we investigated the microstructure of these tracts with fractional anisotropy (FA), radial diffusivity (RD), and hindrance modulated orientational anisotropy (HMOA) in patients with PD and healthy controls. By reducing dimensionality of both the microstructural metrics and the investigated cerebellar and cerebellar-thalamic tracts using principal component analyses, we found global differences between patients with PD and controls, suggestive of higher fractional anisotropy, lower radial diffusivity, and higher hindrance modulated orientational anisotropy in patients. However, separate analyses for each of the tracts did not yield any significant differences. Our findings contribute to the characterization of the distinct anatomical connections between the cerebellum and the diencephalon. Microstructural differences between patients and controls in the cerebellar pathways suggest involvement of these structures in PD, complementing previous functional and diffusion imaging studies.

Mapping the human corticoreticular pathway with multimodal delineation of the gigantocellular reticular nucleus and high-resolution diffusion tractography

The corticoreticular pathway (CRP) is a major motor tract that transmits cortical input to the reticular formation motor nuclei and may be an important mediator of motor recovery after central nervous system damage. However, its cortical origins, trajectory and laterality are incompletely understood in humans. This study aimed to map the human CRP and generate an average CRP template in standard MRI space. Following recently established guidelines, we manually delineated the primary reticular formation motor nucleus (gigantocellular reticular nucleus [GRN]) using several group-mean MRI contrasts from the Human Connectome Project (HCP). CRP tractography was then performed with HCP diffusion-weighted MRI data (N = 1065) by selecting diffusion streamlines that reached both the cortex and GRN. Corticospinal tract (CST) tractography was also performed for comparison. Results suggest that the human CRP has widespread origins, which overlap with the CST across most of the motor cortex and include additional exclusive inputs from the medial and anterior prefrontal cortices. The estimated CRP projected through the anterior and posterior limbs of the internal capsule before partially decussating in the midbrain tegmentum and converging bilaterally on the pontomedullary reticular formation. Thus, the CRP trajectory appears to partially overlap the CST, while being more distributed and anteromedial to the CST in the cerebrum before moving posterior to the CST in the brainstem. These findings have important implications for neurophysiologic testing, cortical stimulation and movement recovery after brain lesions. We expect that our GRN and tract maps will also facilitate future CRP research.

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.

Diffusion MRI and anatomic tracing in the same brain reveal common failure modes of tractography

Anatomic tracing is recognized as a critical source of knowledge on brain circuitry that can be used to assess the accuracy of diffusion MRI (dMRI) tractography. However, most prior studies that have performed such assessments have used dMRI and tracer data from different brains and/or have been limited in the scope of dMRI analysis methods allowed by the data. In this work, we perform a quantitative, voxel-wise comparison of dMRI tractography and anatomic tracing data in the same macaque brain. An ex vivo dMRI acquisition with high angular resolution and high maximum b-value allows us to compare a range of q-space sampling, orientation reconstruction, and tractography strategies. The availability of tracing in the same brain allows us to localize the sources of tractography errors and to identify axonal configurations that lead to such errors consistently, across dMRI acquisition and analysis strategies. We find that these common failure modes involve geometries such as branching or turning, which cannot be modeled well by crossing fibers. We also find that the default thresholds that are commonly used in tractography correspond to rather conservative, low-sensitivity operating points. While deterministic tractography tends to have higher sensitivity than probabilistic tractography in that very conservative threshold regime, the latter outperforms the former as the threshold is relaxed to avoid missing true anatomical connections. On the other hand, the q-space sampling scheme and maximum b-value have less of an impact on accuracy. Finally, using scans from a set of additional macaque brains, we show that there is enough inter-individual variability to warrant caution when dMRI and tracer data come from different animals, as is often the case in the tractography validation literature. Taken together, our results provide insights on the limitations of current tractography methods and on the critical role that anatomic tracing can play in identifying potential avenues for improvement.

The Complex Hodological Architecture of the Macaque Dorsal Intraparietal Areas as Emerging from Neural Tracers and DW-MRI Tractography

In macaque monkeys, dorsal intraparietal areas are involved in several daily visuomotor actions. However, their border and sources of cortical afferents remain loosely defined. Combining retrograde histologic tracing and MRI diffusion-based tractography, we found a complex hodology of the dorsal bank of the intraparietal sulcus (db-IPS), which can be subdivided into a rostral intraparietal area PEip, projecting to the spinal cord, and a caudal medial intraparietal area MIP lacking such projections. Both include an anterior and a posterior sector, emerging from their ipsilateral, gradient-like connectivity profiles. As tractography estimations, we used the cross-sectional area of the white matter bundles connecting each area with other parietal and frontal regions, after selecting regions of interest (ROIs) corresponding to the injection sites of neural tracers. For most connections, we found a significant correlation between the proportions of cells projecting to all sectors of PEip and MIP along the continuum of the db-IPS and tractography. The latter also revealed “false positive” but plausible connections awaiting histologic validation.

Brain connections derived from diffusion MRI tractography can be highly anatomically accurate-if we know where white matter pathways start, where they end, and where they do not go

MR Tractography, which is based on MRI measures of water diffusivity, is currently the only method available for noninvasive reconstruction of fiber pathways in the brain. However, it has several fundamental limitations that call into question its accuracy in many applications. Therefore, there has been intense interest in defining and mitigating the intrinsic limitations of the method. Recent studies have reported that tractography is inherently limited in its ability to accurately reconstruct the connections of the brain, when based on voxel-averaged estimates of local fiber orientation alone. Several validation studies have confirmed that tractography techniques are plagued by both false-positive and false-negative connections. However, these validation studies which quantify sensitivity and specificity, particularly in animal models, have not utilized prior anatomical knowledge, as is done in the human literature, for virtual dissection of white matter pathways, instead assessing tractography implemented in a relatively unconstrained manner. Thus, they represent a worse-case scenario for bundle-segmentation techniques and may not be indicative of the anatomical accuracy in the process of bundle segmentation, where streamline filtering using inclusion and exclusion regions-of-interest is common. With this in mind, the aim of the current study is to investigate and quantify the upper bounds of accuracy using current tractography methods. Making use of the same dataset utilized in two seminal validation papers, we show that prior anatomical knowledge in the form of manually placed or template-driven constraints can significantly improve the anatomical accuracy of estimated brain connections. Thus, we show that it is possible to achieve a high sensitivity and high specificity simultaneously, and conclude that current tractography algorithms, in combination with anatomically driven constraints, can result in reconstructions which very accurately reflect the ground truth white matter connections.

On the cortical connectivity in the macaque brain: A comparison of diffusion tractography and histological tracing data

Diffusion-weighted magnetic resonance imaging (DW-MRI) tractography is a non-invasive tool to probe neural connections and the structure of the white matter. It has been applied successfully in studies of neurological disorders and normal connectivity. Recent work has revealed that tractography produces a high incidence of false-positive connections, often from "bottleneck" white matter configurations. The rich literature in histological connectivity analysis studies in the macaque monkey enables quantitative evaluation of the performance of tractography algorithms. In this study, we use the intricate connections of frontal, cingulate, and parietal areas, well established by the anatomical literature, to derive a symmetrical histological connectivity matrix composed of 59 cortical areas. We evaluate the performance of fifteen diffusion tractography algorithms, including global, deterministic, and probabilistic state-of-the-art methods for the connectivity predictions of 1711 distinct pairs of areas, among which 680 are reported connected by the literature. The diffusion connectivity analysis was performed on a different ex-vivo macaque brain, acquired using multi-shell DW-MRI protocol, at high spatial and angular resolutions. Across all tested algorithms, the true-positive and true-negative connections were dominant over false-positive and false-negative connections, respectively. Moreover, three-quarters of streamlines had endpoints location in agreement with histological data, on average. Furthermore, probabilistic streamline tractography algorithms show the best performances in predicting which areas are connected. Altogether, we propose a method for quantitative evaluation of tractography algorithms, which aims at improving the sensitivity and the specificity of diffusion-based connectivity analysis. Overall, those results confirm the usefulness of tractography in predicting connectivity, although errors are produced. Many of the errors result from bottleneck white matter configurations near the cortical grey matter and should be the target of future implementation of methods.

Open access resource for cellular-resolution analyses of corticocortical connectivity in the marmoset monkey

Understanding the principles of neuronal connectivity requires tools for efficient quantification and visualization of large datasets. The primate cortex is particularly challenging due to its complex mosaic of areas, which in many cases lack clear boundaries. Here, we introduce a resource that allows exploration of results of 143 retrograde tracer injections in the marmoset neocortex. Data obtained in different animals are registered to a common stereotaxic space using an algorithm guided by expert delineation of histological borders, allowing accurate assignment of connections to areas despite interindividual variability. The resource incorporates tools for analyses relative to cytoarchitectural areas, including statistical properties such as the fraction of labeled neurons and the percentage of supragranular neurons. It also provides purely spatial (parcellation-free) data, based on the stereotaxic coordinates of 2 million labeled neurons. This resource helps bridge the gap between high-density cellular connectivity studies in rodents and imaging-based analyses of human brains.

Understanding principles of neuronal connectivity requires tools for quantification and visualization of large datasets. Here, the authors introduce an online resource encompassing the coordinates of two million neurons labelled by tracer injections in the marmoset cortex, and analysis tools.

Comparison between diffusion MRI tractography and histological tract-tracing of cortico-cortical structural connectivity in the ferret brain

The anatomical wiring of the brain is a central focus in network neuroscience. Diffusion MRI tractography offers the unique opportunity to investigate the brain fiber architecture in vivo and noninvasively. However, its reliability is still highly debated. Here, we explored the ability of diffusion MRI tractography to match invasive anatomical tract-tracing connectivity data of the ferret brain. We also investigated the influence of several state-of-the-art tractography algorithms on this match to ground truth connectivity data. Tract-tracing connectivity data were obtained from retrograde tracer injections into the occipital, parietal, and temporal cortices of adult ferrets. We found that the relative densities of projections identified from the anatomical experiments were highly correlated with the estimates from all the studied diffusion tractography algorithms (Spearman's rho ranging from 0.67 to 0.91), while only small, nonsignificant variations appeared across the tractography algorithms. These results are comparable to findings reported in mouse and monkey, increasing the confidence in diffusion MRI tractography results. Moreover, our results provide insights into the variations of sensitivity and specificity of the tractography algorithms, and hence into the influence of choosing one algorithm over another.

Neurite orientation dispersion and density imaging quantifies corticospinal tract microstructural organization in children with unilateral cerebral palsy

Children with unilateral cerebral palsy (UCP) due to early brain injury exhibit disrupted connectivity of corticospinal tracts (CSTs), which can be quantified using diffusion-weighted magnetic resonance imaging (DWI). Diffusion tensor imaging (DTI) is commonly used to quantify white matter organization, however, this model lacks the biological specificity to accurately describe underlying microstructural properties. Newer approaches, such as neurite orientation dispersion and density imaging (NODDI), may provide more biologically accurate information regarding CST microstructure. In this study, we directly compared metrics of CST microstructure using NODDI and DTI models to characterize the microstructural organization of corticospinal pathways. Twenty participants with UCP participating in a neuromodulation/rehabilitation intervention underwent imaging including multi-shell DWI; 10 participants' datasets were adequately completed for neuroimaging analysis. Task fMRI-guided probabilistic tractography from motor cortex to brainstem was performed at baseline and follow-up to reconstruct the CSTs. Diffusion metrics were compared between hemispheres at baseline, and between baseline and follow-up to test for intervention effects. Correlation analyses were used to compare baseline metrics to changes in hand function following the intervention. DTI results showed that mean fractional anisotropy in lesioned and nonlesioned CSTs did not significantly differ, but mean, axial, and radial diffusivity were greater in the lesioned CST. For NODDI, intracellular volume fraction (ICVF) and orientation dispersion index (ODI) were lower in the lesioned CST. Unimanual function was strongly correlated with ICVF, but not FA. NODDI may reveal distinct properties of CST microstructure that are linked to motor function, indicating their potential in characterizing brain structure and development.

Histologically derived fiber response functions for diffusion MRI vary across white matter fibers-An ex vivo validation study in the squirrel monkey brain

Understanding the relationship between the diffusion-weighted MRI signal and the arrangement of white matter fibers is fundamental for accurate voxel-wise reconstruction of the fiber orientation distribution (FOD) and subsequent fiber tractography. Spherical deconvolution reconstruction techniques model the diffusion signal as the convolution of the FOD with a response function that represents the signal profile of a single fiber orientation. Thus, given the signal and a fiber response function, the FOD can be estimated in every imaging voxel by deconvolution. However, the selection of the appropriate response function remains relatively under-studied, and requires further validation. In this work, using 3D histologically defined FODs and the corresponding diffusion signal from three ex vivo squirrel monkey brains, we derive the ground truth response functions. We find that the histologically derived response functions differ from those conventionally used. Next, we find that response functions statistically vary across brain regions, which suggests that the practice of using the same kernel throughout the brain is not optimal. We show that different kernels lead to different FOD reconstructions, which in turn can lead to different tractography results depending on algorithmic parameters, with large variations in the accuracy of resulting reconstructions. Together, these results suggest there is room for improvement in estimating and understanding the relationship between the diffusion signal and the underlying FOD.

Challenges in diffusion MRI tractography - Lessons learned from international benchmark competitions

Diffusion MRI (dMRI) fiber tractography has become a pillar of the neuroimaging community due to its ability to noninvasively map the structural connectivity of the brain. Despite widespread use in clinical and research domains, these methods suffer from several potential drawbacks or limitations. Thus, validating the accuracy and reproducibility of techniques is critical for sound scientific conclusions and effective clinical outcomes. Towards this end, a number of international benchmark competitions, or "challenges", has been organized by the diffusion MRI community in order to investigate the reliability of the tractography process by providing a platform to compare algorithms and results in a fair manner, and evaluate common and emerging algorithms in an effort to advance the state of the field. In this paper, we summarize the lessons from a decade of challenges in tractography, and give perspective on the past, present, and future "challenges" that the field of diffusion tractography faces.

LEARNING 3D WHITE MATTER MICROSTRUCTURE FROM 2D HISTOLOGY

Histological analysis is typically the gold standard for validating measures of tissue microstructure derived from magnetic resonance imaging (MRI) contrasts. However, most histological investigations are inherently 2-dimensional (2D), due to increased field-of-view, higher in-plane resolutions, ease of acquisition, decreased costs, and a large number of available contrasts compared to 3-dimensional (3D) analysis. Because of this, it would be of great interest to be able to learn the 3D tissue microstructure from 2D histology. In this study, we use diffusion MRI (dMRI) of a squirrel monkey brain and corresponding myelin stained sections in combination with a convolution neural network to learn the relationship between the 3D diffusion estimated axonal fiber orientation distributions and the 2D myelin stain. We find that we are able to estimate the 3D fiber distribution with moderate to high angular agreement with the ground truth (median angular correlation coefficients of 0.48 across the unseen slices). This network could be used to validate dMRI neuronal structural measurements in 3D, even if only 2D histology is available for validation. Generalization is possible to transfer this network to human stained sections to infer the 3D fiber distribution at resolutions currently unachievable with dMRI, which would allow diffusion fiber tractography at unprecedented resolutions. We envision the use of similar networks to learn other 3D microstructural measures from an array of potential common 2D histology contrasts.

Limits to anatomical accuracy of diffusion tractography using modern approaches

Diffusion MRI fiber tractography is widely used to probe the structural connectivity of the brain, with a range of applications in both clinical and basic neuroscience. Despite widespread use, tractography has well-known pitfalls that limits the anatomical accuracy of this technique. Numerous modern methods have been developed to address these shortcomings through advances in acquisition, modeling, and computation. To test whether these advances improve tractography accuracy, we organized the 3-D Validation of Tractography with Experimental MRI (3D-VoTEM) challenge at the ISBI 2018 conference. We made available three unique independent tractography validation datasets - a physical phantom and two ex vivo brain specimens - resulting in 176 distinct submissions from 9 research groups. By comparing results over a wide range of fiber complexities and algorithmic strategies, this challenge provides a more comprehensive assessment of tractography's inherent limitations than has been reported previously. The central results were consistent across all sub-challenges in that, despite advances in tractography methods, the anatomical accuracy of tractography has not dramatically improved in recent years. Taken together, our results independently confirm findings from decades of tractography validation studies, demonstrate inherent limitations in reconstructing white matter pathways using diffusion MRI data alone, and highlight the need for alternative or combinatorial strategies to accurately map the fiber pathways of the brain.