Validation of structural brain connectivity networks: The impact of scanning parameters

Macaque Brainnetome Atlas: A multifaceted brain map with parcellation, connection, and histology

The rhesus macaque (Macaca mulatta) is a crucial experimental animal that shares many genetic, brain organizational, and behavioral characteristics with humans. A macaque brain atlas is fundamental to biomedical and evolutionary research. However, even though connectivity is vital for understanding brain functions, a connectivity-based whole-brain atlas of the macaque has not previously been made. In this study, we created a new whole-brain map, the Macaque Brainnetome Atlas (MacBNA), based on the anatomical connectivity profiles provided by high angular and spatial resolution ex vivo diffusion MRI data. The new atlas consists of 248 cortical and 56 subcortical regions as well as their structural and functional connections. The parcellation and the diffusion-based tractography were evaluated with invasive neuronal-tracing and Nissl-stained images. As a demonstrative application, the structural connectivity divergence between macaque and human brains was mapped using the Brainnetome atlases of those two species to uncover the genetic underpinnings of the evolutionary changes in brain structure. The resulting resource includes: (1) the thoroughly delineated Macaque Brainnetome Atlas (MacBNA), (2) regional connectivity profiles, (3) the postmortem high-resolution macaque diffusion and T2-weighted MRI dataset (Brainnetome-8), and (4) multi-contrast MRI, neuronal-tracing, and histological images collected from a single macaque. MacBNA can serve as a common reference frame for mapping multifaceted features across modalities and spatial scales and for integrative investigation and characterization of brain organization and function. Therefore, it will enrich the collaborative resource platform for nonhuman primates and facilitate translational and comparative neuroscience research.

Evolutionary continuity and divergence of auditory dorsal and ventral pathways in primates revealed by ultra-high field diffusion MRI

Auditory dorsal and ventral pathways in the human brain play important roles in supporting speech and language processing. However, the evolutionary root of the dual auditory pathways in the primate brain is unclear. By parcellating the auditory cortex of marmosets (a New World monkey species), macaques (an Old World monkey species), and humans using the same individual-based analysis method and tracking the pathways from the auditory cortex based on multi-shell diffusion-weighted MRI (dMRI), homologous auditory dorsal and ventral fiber tracks were identified in these primate species. The ventral pathway was found to be well conserved in all three primate species analyzed but extend to more anterior temporal regions in humans. In contrast, the dorsal pathway showed a divergence between monkey and human brains. First, frontal regions in the human brain have stronger connections to the higher-level auditory regions than to the lower-level auditory regions along the dorsal pathway, while frontal regions in the monkey brain show opposite connection patterns along the dorsal pathway. Second, the left lateralization of the dorsal pathway is only found in humans. Moreover, the connectivity strength of the dorsal pathway in marmosets is more similar to that of humans than macaques. These results demonstrate the continuity and divergence of the dual auditory pathways in the primate brains along the evolutionary path, suggesting that the putative neural networks supporting human speech and language processing might have emerged early in primate evolution.

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.

Towards reliable reconstruction of the mouse brain corticothalamic connectivity using diffusion MRI

Diffusion magnetic resonance imaging (dMRI) tractography has yielded intriguing insights into brain circuits and their relationship to behavior in response to gene mutations or neurological diseases across a number of species. Still, existing tractography approaches suffer from limited sensitivity and specificity, leading to uncertain interpretation of the reconstructed connections. Hence, in this study, we aimed to optimize the imaging and computational pipeline to achieve the best possible spatial overlaps between the tractography and tracer-based axonal projection maps within the mouse brain corticothalamic network. We developed a dMRI-based atlas of the mouse forebrain with structural labels imported from the Allen Mouse Brain Atlas (AMBA). Using the atlas and dMRI tractography, we first reconstructed detailed node-to-node mouse brain corticothalamic structural connectivity matrices using different imaging and tractography parameters. We then investigated the effects of each condition for accurate reconstruction of the corticothalamic projections by quantifying the similarities between the tractography and the tracer data from the Allen Mouse Brain Connectivity Atlas (AMBCA). Our results suggest that these parameters significantly affect tractography outcomes and our atlas can be used to investigate macroscopic structural connectivity in the mouse brain. Furthermore, tractography in mouse brain gray matter still face challenges and need improved imaging and tractography methods.

Optimization of quasi-diffusion magnetic resonance imaging for quantitative accuracy and time-efficient acquisition

PURPOSE

Quasi-diffusion MRI (QDI) is a novel quantitative technique based on the continuous time random walk model of diffusion dynamics. QDI provides estimates of the diffusion coefficient, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> in mm²  s^-1 and a fractional exponent, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> , defining the non-Gaussianity of the diffusion signal decay. Here, the b-value selection for rapid clinical acquisition of QDI tensor imaging (QDTI) data is optimized.

METHODS

Clinically appropriate QDTI acquisitions were optimized in healthy volunteers with respect to a multi-b-value reference (MbR) dataset comprising 29 diffusion-sensitized images arrayed between <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>b</mml:mi> <mml:mo>=</mml:mo> <mml:mn>0</mml:mn></mml:mrow> <mml:annotation>$$ b=0 $$</mml:annotation></mml:semantics> </mml:math> and 5000 s mm^-2 . The effects of varying maximum b-value ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> ), number of b-value shells, and the effects of Rician noise were investigated.

RESULTS

QDTI measures showed <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> dependence, most significantly for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, which monotonically decreased with higher <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> leading to improved tissue contrast. Optimized 2 b-value shell acquisitions showed small systematic differences in QDTI measures relative to MbR values, with overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ \kern0.50em {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and underestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, and overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> anisotropies in gray and white matter. Additional shells improved the accuracy, precision, and reliability of QDTI estimates with 3 and 4 shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>5000</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=5000 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , and 4 b-value shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>3960</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=3960 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , providing minimal bias in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> compared to the MbR.

CONCLUSION

A highly detailed optimization of non-Gaussian dMRI for in vivo brain imaging was performed. QDI provided robust parameterization of non-Gaussian diffusion signal decay in clinically feasible imaging times with high reliability, accuracy, and precision of QDTI measures.

Improving EEG-Based Driver Distraction Classification Using Brain Connectivity Estimators

This paper discusses a novel approach to an EEG (electroencephalogram)-based driver distraction classification by using brain connectivity estimators as features. Ten healthy volunteers with more than one year of driving experience and an average age of 24.3 participated in a virtual reality environment with two conditions, a simple math problem-solving task and a lane-keeping task to mimic the distracted driving task and a non-distracted driving task, respectively. Independent component analysis (ICA) was conducted on the selected epochs of six selected components relevant to the frontal, central, parietal, occipital, left motor, and right motor areas. Granger-Geweke causality (GGC), directed transfer function (DTF), partial directed coherence (PDC), and generalized partial directed coherence (GPDC) brain connectivity estimators were used to calculate the connectivity matrixes. These connectivity matrixes were used as features to train the support vector machine (SVM) with the radial basis function (RBF) and classify the distracted and non-distracted driving tasks. GGC, DTF, PDC, and GPDC connectivity estimators yielded the classification accuracies of 82.27%, 70.02%, 86.19%, and 80.95%, respectively. Further analysis of the PDC connectivity estimator was conducted to determine the best window to differentiate between the distracted and non-distracted driving tasks. This study suggests that the PDC connectivity estimator can yield better classification accuracy for driver distractions.

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.

Intra-Areal Visual Topography in Primate Brains Mapped with Probabilistic Tractography of Diffusion-Weighted Imaging

Noninvasive diffusion-weighted magnetic resonance imaging (dMRI) can be used to map the neural connectivity between distinct areas in the intact brain, but the standard resolution achieved fundamentally limits the sensitivity of such maps. We investigated the sensitivity and specificity of high-resolution postmortem dMRI and probabilistic tractography in rhesus macaque brains to produce retinotopic maps of the lateral geniculate nucleus (LGN) and extrastriate cortical visual area V5/MT based on their topographic connections with the previously established functional retinotopic map of primary visual cortex (V1). We also replicated the differential connectivity of magnocellular and parvocellular LGN compartments with V1 across visual field positions. Predicted topographic maps based on dMRI data largely matched the established retinotopy of both LGN and V5/MT. Furthermore, tractography based on in vivo dMRI data from the same macaque brains acquired at standard field strength (3T) yielded comparable topographic maps in many cases. We conclude that tractography based on dMRI is sensitive enough to reveal the intrinsic organization of ordered connections between topographically organized neural structures and their resultant functional organization.

Uncovering Cortical Units of Processing From Multi-Layered Connectomes

Modern diffusion and functional magnetic resonance imaging (dMRI/fMRI) provide non-invasive high-resolution images from which multi-layered networks of whole-brain structural and functional connectivity can be derived. Unfortunately, the lack of observed correspondence between the connectivity profiles of the two modalities challenges the understanding of the relationship between the functional and structural connectome. Rather than focusing on correspondence at the level of connections we presently investigate correspondence in terms of modular organization according to shared canonical processing units. We use a stochastic block-model (SBM) as a data-driven approach for clustering high-resolution multi-layer whole-brain connectivity networks and use prediction to quantify the extent to which a given clustering accounts for the connectome within a modality. The employed SBM assumes a single underlying parcellation exists across modalities whilst permitting each modality to possess an independent connectivity structure between parcels thereby imposing concurrent functional and structural units but different structural and functional connectivity profiles. We contrast the joint processing units to their modality specific counterparts and find that even though data-driven structural and functional parcellations exhibit substantial differences, attributed to modality specific biases, the joint model is able to achieve a consensus representation that well accounts for both the functional and structural connectome providing improved representations of functional connectivity compared to using functional data alone. This implies that a representation persists in the consensus model that is shared by the individual modalities. We find additional support for this viewpoint when the anatomical correspondence between modalities is removed from the joint modeling. The resultant drop in predictive performance is in general substantial, confirming that the anatomical correspondence of processing units is indeed present between the two modalities. Our findings illustrate how multi-modal integration admits consensus representations well-characterizing each individual modality despite their biases and points to the importance of multi-layered connectomes as providing supplementary information regarding the brain's canonical processing units.

Directed functional and structural connectivity in a large-scale model for the mouse cortex

Inferring the structural connectivity from electrophysiological measurements is a fundamental challenge in systems neuroscience. Directed functional connectivity measures, such as the generalized partial directed coherence (GPDC), provide estimates of the causal influence between areas. However, the relation between causality estimates and structural connectivity is still not clear. We analyzed this problem by evaluating the effectiveness of GPDC to estimate the connectivity of a ground-truth, data-constrained computational model of a large-scale network model of the mouse cortex. The model contains 19 cortical areas composed of spiking neurons, with areas connected by long-range projections with weights obtained from a tract-tracing cortical connectome. We show that GPDC values provide a reasonable estimate of structural connectivity, with an average Pearson correlation over simulations of 0.74. Moreover, even in a typical electrophysiological recording scenario containing five areas, the mean correlation was above 0.6. These results suggest that it may be possible to empirically estimate structural connectivity from functional connectivity even when detailed whole-brain recordings are not achievable.

We analyzed the relationship between structural and directed functional connectivity by evaluating the effectiveness of generalized partial directed coherence (GPDC) to estimate the connectivity of a ground-truth, data-constrained computational model of a large-scale network model of the mouse cortex. We show that GPDC values provide a reasonable estimate of structural connectivity even in a typical electrophysiological recording scenario containing few areas. These results suggest that it may be possible to empirically estimate structural connectivity from functional connectivity even when detailed whole-brain recordings are not achievable.

Prevalence of white matter pathways coming into a single white matter voxel orientation: The bottleneck issue in tractography

Characterizing and understanding the limitations of diffusion MRI fiber tractography is a prerequisite for methodological advances and innovations which will allow these techniques to accurately map the connections of the human brain. The so-called "crossing fiber problem" has received tremendous attention and has continuously triggered the community to develop novel approaches for disentangling distinctly oriented fiber populations. Perhaps an even greater challenge occurs when multiple white matter bundles converge within a single voxel, or throughout a single brain region, and share the same parallel orientation, before diverging and continuing towards their final cortical or sub-cortical terminations. These so-called "bottleneck" regions contribute to the ill-posed nature of the tractography process, and lead to both false positive and false negative estimated connections. Yet, as opposed to the extent of crossing fibers, a thorough characterization of bottleneck regions has not been performed. The aim of this study is to quantify the prevalence of bottleneck regions. To do this, we use diffusion tractography to segment known white matter bundles of the brain, and assign each bundle to voxels they pass through and to specific orientations within those voxels (i.e. fixels). We demonstrate that bottlenecks occur in greater than 50-70% of fixels in the white matter of the human brain. We find that all projection, association, and commissural fibers contribute to, and are affected by, this phenomenon, and show that even regions traditionally considered "single fiber voxels" often contain multiple fiber populations. Together, this study shows that a majority of white matter presents bottlenecks for tractography which may lead to incorrect or erroneous estimates of brain connectivity or quantitative tractography (i.e., tractometry), and underscores the need for a paradigm shift in the process of tractography and bundle segmentation for studying the fiber pathways of the human brain.

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.

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.

Current and Future Challenges of Functional MRI and Diffusion Tractography in the Surgical Setting: From Eloquent Brain Mapping to Neural Plasticity

Decades ago, Spetzler (1986) and Sawaya (1998) provided a rough brain segmentation of the eloquent areas of the brain, aimed to help surgical decisions in cases of vascular malformations and tumors, respectively. Currently in clinical use, their criteria are in need of revision. Defining functions (eg, sensorimotor, language and visual) that should be preserved during surgery seems a straightforward task. In practice, locating the specific areas that could cause a permanent vs transient deficit is not an easy task. This is particularly true for the associative cortex and cognitive domains such as language. The old model, with Broca's and Wernicke's areas at the forefront, has been superseded by a dual-stream model of parallel language processing; named ventral and dorsal pathways. This complicated network of cortical hubs and subcortical white matter pathways needing preservation during surgery is a work in progress. Preserving not only cortical regions but most importantly preserving the connections, or white matter fiber bundles, of core regions in the brain is the new paradigm. For instance, the arcuate fascicululs and inferior fronto-occipital fasciculus are key components of the dorsal and ventral language pathways, respectively; and their damage result in permanent language deficits. Interestedly, the damage of the temporal portions of these bundles -where there is a crossroad with other multiple bundles-, appears to be more important (permanent) than the damage of the frontal portions - where plasticity and contralateral activation could help. Although intraoperative direct cortical and subcortical stimulation have contributed largely, advanced MR techniques such as functional MRI (fMRI) and diffusion tractography (DT), are at the epi-center of our current understanding. Nevertheless, these techniques posse important challenges: such as neurovascular uncoupling or venous bias on fMRI; and appropriate anatomical validation or accurate representation of crossing fibers on DT. These limitations should be well understood and taken into account in clinical practice. Unifying multidisciplinary research and clinical efforts is desirable, so these techniques could contribute more efficiently not only to locate eloquent areas but to improve outcomes and our understanding of neural plasticity. Finally, although there are constant anatomical and functional regions at the individual level, there is a known variability at the inter-individual level. This concept should strengthen the importance of a personalized approach when evaluating these regions on fMRI and DT. It should strengthen the importance of personalized treatments as well, aimed to meet tailored needs and expectations.

Ventralis intermedius nucleus anatomical variability assessment by MRI structural connectivity

The ventralis intermedius nucleus (Vim) is centrally placed in the dentato-thalamo-cortical pathway (DTCp) and is a key surgical target in the treatment of severe medically refractory tremor. It is not visible on conventional MRI sequences; consequently, stereotactic targeting currently relies on atlas-based coordinates. This fails to capture individual anatomical variability, which may lead to poor long-term clinical efficacy. Probabilistic tractography, combined with known anatomical connectivity, enables localisation of thalamic nuclei at an individual subject level. There are, however, a number of confounds associated with this technique that may influence results. Here we focused on an established method, using probabilistic tractography to reconstruct the DTCp, to identify the connectivity-defined Vim (cd-Vim) in vivo. Using 100 healthy individuals from the Human Connectome Project, our aim was to quantify cd-Vim variability across this population, measure the discrepancy with atlas-defined Vim (ad-Vim), and assess the influence of potential methodological confounds. We found no significant effect of any of the confounds. The mean cd-Vim coordinate was located within 1.88 mm (left) and 2.12 mm (right) of the average midpoint and 3.98 mm (left) and 5.41 mm (right) from the ad-Vim coordinates. cd-Vim location was more variable on the right, which reflects hemispheric asymmetries in the probabilistic DTC reconstructed. The method was reproducible, with no significant cd-Vim location differences in a separate test-retest cohort. The superior cerebellar peduncle was identified as a potential source of artificial variance. This work demonstrates significant individual anatomical variability of the cd-Vim that atlas-based coordinate targeting fails to capture. This variability was not related to any methodological confound tested. Lateralisation of cerebellar functions, such as speech, may contribute to the observed asymmetry. Tractography-based methods seem sensitive to individual anatomical variability that is missed by conventional neurosurgical targeting; these findings may form the basis for translational tools to improve efficacy and reduce side-effects of thalamic surgery for tremor.

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.

Using connectomics for predictive assessment of brain parcellations

The organization of the human brain remains elusive, yet is of great importance to the mechanisms of integrative brain function. At the macroscale, its structural and functional interpretation is conventionally assessed at the level of cortical units. However, the definition and validation of such cortical parcellations are problematic due to the absence of a true gold standard. We propose a framework for quantitative evaluation of brain parcellations via statistical prediction of connectomics data. Specifically, we evaluate the extent in which the network representation at the level of cortical units (defined as parcels) accounts for high-resolution brain connectivity. Herein, we assess the pertinence and comparative ranking of ten existing parcellation atlases to account for functional (FC) and structural connectivity (SC) data based on data from the Human Connectome Project (HCP), and compare them to data-driven as well as spatially-homogeneous geometric parcellations including geodesic parcellations with similar size distributions as the atlases. We find substantial discrepancy in parcellation structures that well characterize FC and SC and differences in what well represents an individual's functional connectome when compared against the FC structure that is preserved across individuals. Surprisingly, simple spatial homogenous parcellations generally provide good representations of both FC and SC, but are inferior when their within-parcellation distribution of individual parcel sizes is matched to that of a valid atlas. This suggests that the choice of fine grained and coarse representations used by existing atlases are important. However, we find that resolution is more critical than the exact border location of parcels.

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.