Validation of tractography: Comparison with manganese tracing

Beyond the surface: how ex-vivo diffusion-weighted imaging reveals large animal brain microstructure and connectivity

Diffusion-weighted Imaging (DWI) is an effective and state-of-the-art neuroimaging method that non-invasively reveals the microstructure and connectivity of tissues. Recently, novel applications of the DWI technique in studying large brains through ex-vivo imaging enabled researchers to gain insights into the complex neural architecture in different species such as those of Perissodactyla (e.g., horses and rhinos), Artiodactyla (e.g., bovids, swines, and cetaceans), and Carnivora (e.g., felids, canids, and pinnipeds). Classical in-vivo tract-tracing methods are usually considered unsuitable for ethical and practical reasons, in large animals or protected species. Ex-vivo DWI-based tractography offers the chance to examine the microstructure and connectivity of formalin-fixed tissues with scan times and precision that is not feasible in-vivo. This paper explores DWI's application to ex-vivo brains of large animals, highlighting the unique insights it offers into the structure of sometimes phylogenetically different neural networks, the connectivity of white matter tracts, and comparative evolutionary adaptations. Here, we also summarize the challenges, concerns, and perspectives of ex-vivo DWI that will shape the future of the field in large brains.

Alterations of White Matter Connectivity in Adults with Essential Hypertension

Purpose

To explore the topology of the white matter network in individuals with essential hypertension by graph theory.

Patients and Methods

T1-weighted image and diffusion tensor imaging (DTI) data from 43 patients diagnosed with essential hypertension (EHT) and 33 individuals with normotension (healthy controls, HCs) were incorporated in this cross-sectional study. Furthermore, structural networks were constructed by graph theory to calculate whole brain network characteristics and intracerebral node characteristics.

Results

Both EHT and HC groups displayed small-worldness in their structural networks. The area under the curve (AUC) of the small-worldness coefficient (σ) was higher in the EHT group compared to the HC group, whereas the AUC of assortativity was lower in the EHT group in contrast to the HC group. The nodal clustering coefficient (CP) and local efficiency (Eloc) of the EHT group decreased in the right dorsolateral superior frontal gyrus and the left medial superior frontal gyrus. These values increased in the left anterior cingulate and paracingulate gyrus. Furthermore, weight and body mass index (BMI) were positively correlated with σ.

Conclusion

The EHT group showed brain network separation and integration dysfunction. Weight and BMI were positively correlated with σ. The data acquired in this investigation implied that altered structural connectivity in the prefrontal region may be a potential neuroimaging marker in EHT patients.

Whole-brain, gray, and white matter time-locked functional signal changes with simple tasks and model-free analysis

Recent studies have revealed the production of time-locked blood oxygenation level-dependent (BOLD) functional MRI (fMRI) signals throughout the entire brain in response to tasks, challenging the existence of sparse and localized brain functions and highlighting the pervasiveness of potential false negative fMRI findings. "Whole-brain" actually refers to gray matter, the only tissue traditionally studied with fMRI. However, several reports have demonstrated reliable detection of BOLD signals in white matter, which have previously been largely ignored. Using simple tasks and analyses, we demonstrate BOLD signal changes across the whole brain, in both white and gray matters, in similar manner to previous reports of whole brain studies. We investigated whether white matter displays time-locked BOLD signals across multiple structural pathways in response to a stimulus in a similar manner to the cortex. We find that both white and gray matter show time-locked activations across the whole brain, with a majority of both tissue types showing statistically significant signal changes for all task stimuli investigated. We observed a wide range of signal responses to tasks, with different regions showing different BOLD signal changes to the same task. Moreover, we find that each region may display different BOLD responses to different stimuli. Overall, we present compelling evidence that, just like all gray matter, essentially all white matter in the brain shows time-locked BOLD signal changes in response to multiple stimuli, challenging the idea of sparse functional localization and the prevailing wisdom of treating white matter BOLD signals as artifacts to be removed.

PPA: Principal parcellation analysis for brain connectomes and multiple traits

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

Variational Autoencoders for Generating Synthetic Tractography-Based Bundle Templates in a Low-Data Setting

White matter tracts generated from whole brain tractography are often processed using automatic segmentation methods with standard atlases. Atlases are generated from hundreds of subjects, which becomes time-consuming to create and difficult to apply to all populations. In this study, we extended our prior work on using a deep generative model - a Convolutional Variational Autoencoder - to map complex and data-intensive streamlines to a low-dimensional latent space given a limited sample size of 50 subjects from the ADNI3 dataset, to generate synthetic population-specific bundle templates using Kernel Density Estimation (KDE) on streamline embeddings. We conducted a quantitative shape analysis by calculating bundle shape metrics, and found that our bundle templates better capture the shape distribution of the bundles than the atlas data used in the original segmentation derived from young healthy adults. We further demonstrated the use of our framework for direct bundle segmentation from whole-brain tractograms.

Variational Autoencoders for Generating Synthetic Tractography-Based Bundle Templates in a Low-Data Setting

White matter tracts generated from whole brain tractography are often processed using automatic segmentation methods with standard atlases. Atlases are generated from hundreds of subjects, which becomes time-consuming to create and difficult to apply to all populations. In this study, we extended our prior work on using a deep generative model a Convolutional Variational Autoencoder - to map complex and data-intensive streamlines to a low-dimensional latent space given a limited sample size of 50 subjects from the ADNI3 dataset, to generate synthetic population-specific bundle templates using Kernel Density Estimation (KDE) on streamline embeddings. We conducted a quantitative shape analysis by calculating bundle shape metrics, and found that our bundle templates better capture the shape distribution of the bundles than the atlas data used in the original segmentation derived from young healthy adults. We further demonstrated the use of our framework for direct bundle segmentation from whole-brain tractograms.

Whole-brain, gray and white matter time-locked functional signal changes with simple tasks and model-free analysis

Recent studies have revealed the production of time-locked blood oxygenation-level dependent (BOLD) functional MRI (fMRI) signals throughout the entire brain in response to a task, challenging the idea of sparse and localized brain functions, and highlighting the pervasiveness of potential false negative fMRI findings. In these studies, 'whole-brain' refers to gray matter regions only, which is the only tissue traditionally studied with fMRI. However, recent reports have also demonstrated reliable detection and analyses of BOLD signals in white matter which have been largely ignored in previous reports. Here, using model-free analysis and simple tasks, we investigate BOLD signal changes in both white and gray matters. We aimed to evaluate whether white matter also displays time-locked BOLD signals across all structural pathways in response to a stimulus. We find that both white and gray matter show time-locked activations across the whole-brain, with a majority of both tissue types showing statistically significant signal changes for all task stimuli investigated. We observed a wide range of signal responses to tasks, with different regions showing very different BOLD signal changes to the same task. Moreover, we find that each region may display different BOLD responses to different stimuli. Overall, we present compelling evidence that the whole brain, including both white and gray matter, show time-locked activation to multiple stimuli, not only challenging the idea of sparse functional localization, but also the prevailing wisdom of treating white matter BOLD signals as artefacts to be removed.

Fixel-based analysis of the diffusion properties of the patients with brain injury and chronic health symptoms

The diffusion properties from diffusion tensor imaging (DTI) are sensitive to white matter (WM) abnormalities and could serve as indicators of diffuse axonal damages incurred during a traumatic brain injury (TBI). Analyses of diffusion metrics in the regions of interest (ROIs) were used to compare the differences in the 18 major fiber tracts in 46 participants, between TBI participants with (n = 17) or without (n = 16) chronic symptoms (CS) and a control group (CG, n = 13). In addition to the widely used diffusion metrics, such as fractional anisotropy (FA), mean (MD), axial (AD) and radial (RD) diffusivities, apparent fiber density (AFD), complexity (CX) and fixel number (FN) derived from Mrtrix3 software package were used to characterize WM tracts and compare between participant groups in the ROIs defined by the fixel numbers. Significant differences were found in FA, AFD, MD, RD and CX in ROIs with different FNs in the corpus callosum forceps minor, left and right inferior longitudinal fasciculus, and left and right uncinate fasciculus for both TBI groups compared to controls. Diffusion properties in ROIs with different FNs can serve as detailed biomarkers of WM abnormalities, especially for individuals with chronic TBI related symptoms.

Diffusion Tensor Imaging as a Prognostic Tool for Recovery in Acute and Hyperacute Stroke

Stroke represents a major cause of mortality and long-term disability among adult populations, leaving a devastating socioeconomic impact globally. Clinical manifestation of stroke is characterized by great diversity, ranging from minor disability to considerable neurological impairment interfering with activities of daily living and even death. Prognostic ambiguity has stimulated the interest for implementing stroke recovery biomarkers, including those provided by structural neuroimaging techniques, i.e., diffusion tensor imaging (DTI) and tractography for the study of white matter (WM) integrity. Considering the necessity of prompt and accurate prognosis in stroke survivors along with the potential capacity of DTI as a relevant imaging biomarker, the purpose of our study was to review the pertinent literature published within the last decade regarding DTI as a prognostic tool for recovery in acute and hyperacute stroke. We conducted a thorough literature search in two databases (MEDLINE and Science Direct) in order to trace all relevant studies published between 1 January 2012 and 16 March 2022 using predefined terms as key words. Only full-text human studies published in the English language were included. Forty-four studies were identified and are included in this review. We present main findings and by describing several methodological issues, we highlight shortcomings and gaps in the current literature so that research priorities for future research can be outlined. Our review suggests that DTI can track longitudinal changes and identify prognostic correlates in acute and hyperacute stroke patients.

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.

Cingulum bundle connectivity in treatment-refractory compared to treatment-responsive patients with bipolar disorder and healthy controls: a tractography and surgical targeting analysis

OBJECTIVE

The clinical response of patients with bipolar disorder to medical treatment is variable. A better understanding of the underlying neural circuitry involved in bipolar treatment responsivity subtypes may provide insight into treatment resistance and aid in identifying an effective surgical target for deep brain stimulation (DBS) specific to the disorder. Despite considerable imaging research related to the disease, a paucity of comparative imaging analyses of treatment responsiveness exists. There are also no DBS targets designed expressly for patients with bipolar disorder. Therefore, the authors analyzed cingulum bundle axonal connectivity in relation to cortico-striatal-thalamo-cortical (CSTC) loops implicated in bipolar disorder across subjects who are responsive to treatment (RSP) and those who are refractory to therapy (REF), compared to healthy controls (HCs).

METHODS

Twenty-five subjects with bipolar disorder (13 RSP and 12 REF), diagnosed using the Mini International Neuropsychiatric Interview and classified with standardized rating scales, and 14 HCs underwent MRI with diffusion sequences for probabilistic diffusion-weighted tractography analysis. Image processing and tractography were performed using MRTrix. Region of interest (ROI) masks were created manually for 10 anterior cingulum bundle subregions, including surgical targets previously evaluated for the treatment of bipolar disorder (cingulotomy and subgenual cingulate DBS targets). Cortical and subcortical ROIs of brain areas thought to be associated with bipolar disorder and described in animal tract-tracing models were created via FreeSurfer. The number of axonal projections from the cingulum bundle subregion ROIs to cortical/subcortical ROIs for each group was compared.

RESULTS

Significant differences were found across groups involving cingulum bundle and CSTC loops. Subjects in the RSP group had increased connections from rostral cingulum bundle to medial orbitofrontal cortex, which is part of the limbic CSTC loop, whereas subjects in the REF group had increased connectivity from rostral cingulum bundle to thalamus. Additionally, compared to HCs, both RSP and REF subjects had decreased cingulum bundle dorsal connectivity (dorsal anterior/posterior cingulate, dorsomedial/lateral frontal cortex) and increased cingulum bundle ventral connectivity (subgenual cingulate, frontal pole, lateral orbitofrontal cortex) involving limbic and associative CSTC loops.

CONCLUSIONS

Findings demonstrate that bipolar treatment responsivity may be associated with significant differences in cingulum bundle connectivity in relation to CSTC loops, which may help identify a surgical target for bipolar disorder treatment via DBS in the future.

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.

Investigating secondary white matter degeneration following ischemic stroke by modelling affected fiber tracts

•

Secondary white matter degeneration was studied in 11 ischemic stroke patients.

•

We used a custom-developed approach to model damaged fibers associated with a lesion.

•

This approach tackles the inter-subject variability in lesion size and location.

•

Findings suggest that secondary degeneration spreads along an entire fiber’s length.

Secondary white matter degeneration is a common occurrence after ischemic stroke, as identified by Diffusion Tensor Imaging (DTI). However, despite recent advances, the time course of the process is not completely understood. The primary aim of this study was to assess secondary degeneration using an approach whereby we create a patient-specific model of damaged fibers based on the volumetric characteristics of lesions. We also examined the effects of secondary degeneration along the modelled streamlines at different distances from the primary infarction using DTI. Eleven patients who presented with upper limb motor deficits at the time of a first-ever ischemic stroke were included. They underwent scanning at weeks 6 and 29 post-stroke. The fractional anisotropy (FA), mean diffusivity (MD), primary eigenvalue (λ₁), and transverse eigenvalue (λ₂₃) were measured. Using regions of interest based on the simulation output, the differences between the modelled fibers and matched contralateral areas were analyzed. The longitudinal change between the two time points and across five distances from the primary lesion was also assessed using the ratios of diffusion quantities (rFA, rMD, rλ₁, and rλ₂₃) between the ipsilesional and contralesional hemisphere. At week 6 post-stroke, significantly decreased λ₁ was found along the ipsilesional corticospinal tract (CST) with a trend towards lower FA, reduced MD and λ₂₃. At week 29 post-stroke, significantly decreased FA was shown relative to the non-lesioned side, with a trend towards lower λ₁, unchanged MD, and higher λ₂₃. Along the ipsilesional tract, the rFA diminished, whereas the rMD, rλ₁, and rλ₂₃ significantly increased over time. No significant variations in the time progressive effect with distance were demonstrated. The findings support previously described mechanisms of secondary degeneration and suggest that it spreads along the entire length of a damaged tract. Future investigations using higher-order tractography techniques can further explain the intravoxel alterations caused by ischemic injury.

Chronic recording of cortical activity underlying vocalization in awake minipigs

BACKGROUND

Investigating brain dynamics underlying vocal production in animals is a powerful way to inform on the neural bases of human speech. In particular, brain networks underlying vocal production in non-human primates show striking similarities with the human speech production network. However, despite increasing findings also in birds and more recently in rodents, the extent to which the primate vocal cortical network model generalizes to other non-primate mammals remains unclear. Especially, no domestic species has yet been proposed to investigate vocal brain activity using electrophysiological approaches.

NEW METHOD

In the present study, we introduce a novel experimental paradigm to identify the cortical dynamics underlying vocal production in behaving minipigs. A key problem to chronically implant cortical probes in pigs is the presence and growth of frontal sinuses extending caudally to the parietal bone and preventing safe access to neural structures with conventional craniotomy in adult animals.

RESULTS

Here we first show that implantations of soft ECoG grids can be done safely using conventional craniotomy in minipigs younger than 5 months, a period when sinuses are not yet well developed. Using wireless recordings in behaving animals, we further show activation of the motor and premotor cortex around the onset of vocal production of grunts, the most common vocalization of pigs.

CONCLUSION

These results suggest that minipigs, which are very loquacious and social animals, can be a good experimental large animal model to study the cortical bases of vocal production.

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.

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.

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.

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.

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.

Connectome analysis of male world-class gymnasts using probabilistic multishell, multitissue constrained spherical deconvolution tracking

In athletes, long‐term intensive training has been shown to increase unparalleled athletic ability and might induce brain plasticity. We evaluated the structural connectome of world‐class gymnasts (WCGs), as mapped by diffusion‐weighted magnetic resonance imaging probabilistic tractography and a multishell, multitissue constrained spherical deconvolution method to increase the precision of tractography at the tissue interfaces. The connectome was mapped in 10 Japanese male WCGs and in 10 age‐matched male controls. Network‐based statistic identified subnetworks with increased connectivity density in WCGs, involving the sensorimotor, default mode, attentional, visual, and limbic areas. It also revealed a significant association between the structural connectivity of some brain structures with functions closely related to the gymnastic skills and the D‐score, which is used as an index of the gymnasts' specific physical abilities for each apparatus. Furthermore, graph theory analysis demonstrated the characteristics of brain anatomical topology in the WCGs. They displayed significantly increased global connectivity strength with decreased characteristic path length at the global level and higher nodal strength and degree in the sensorimotor, default mode, attention, and limbic/subcortical areas at the local level as compared with controls. Together, these findings extend the current understanding of neural mechanisms that distinguish WCGs from controls and suggest brain anatomical network plasticity in WCGs resulting from long‐term intensive training. Future studies should assess the contribution of genetic or early‐life environmental factors in the brain network organization of WCGs. Furthermore, the indices of brain topology (i.e., connection density and graph theory indices) could become markers for the objective evaluation of gymnastic performance.

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.

The Topography of the Frontal Terminations of the Uncinate Fasciculus Revisited Through Focused Fiber Dissections: Shedding Light on a Current Controversy and Introducing the Insular Apex as a Key Anatomoclinical Area

BACKGROUND

Recent studies advocate a connectivity pattern wider than previously believed of the uncinate fasciculus that extends to the ventrolateral and dorsolateral prefrontal cortices. These new percepts on the connectivity of the tract suggest a more expansive role for the uncinate fasciculus. Our aim was to shed light on this controversy through fiber dissections.

METHODS

Twenty normal adult human formalin-fixed cerebral hemispheres were used. Focused dissections on the insular, orbitofrontal, ventromedial, ventrolateral, and dorsolateral prefrontal areas were performed to record the topography of the frontal terminations of the uncinate fasciculus.

RESULTS

Three discrete fiber layers were consistently disclosed: the first layer was recorded to terminate at the posterior orbital gyrus and pars orbitalis, the second layer at the posterior two thirds of the gyrus rectus, and the last layer at the posterior one third of the paraolfactory gyrus. The insular apex was documented as a crucial landmark regarding the topographic differentiation of the uncinate and occipitofrontal fasciculi (i.e., fibers that travel ventrally belong to the uncinate fasciculus whereas those traveling dorsally are occipitofrontal fibers).

CONCLUSIONS

The frontal terminations of the uncinate fasciculus were consistently documented to project to the posterior orbitofrontal area. The area of the insular apex is introduced for the first time as a crucial surface landmark to effectively distinguish the stems of the uncinate and occipitofrontal fasciculi. This finding could refine the spatial resolution of awake subcortical mapping, especially for insular lesions, and improve the accuracy of in vivo diffusion tensor imaging protocols.

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.

Structural brain disconnectivity mapping of post-stroke fatigue

•

We tested for associations between post stroke fatigue (PSF) and both lesion characteristics and brain structural disconnectome in 84 S patients.

•

Results provided no evidence supporting a simple association between PSF severity and lesion characteristics or disconnectivity.

•

PSF was strongly correlated with depression.

•

Further studies including patients with more severe symptoms are needed to generalize the findings across a wider clinical spectrum.

Stroke patients commonly suffer from post stroke fatigue (PSF). Despite a general consensus that brain perturbations constitute a precipitating event in the multifactorial etiology of PSF, the specific predictive value of conventional lesion characteristics such as size and localization remains unclear. The current study represents a novel approach to assess the neural correlates of PSF in chronic stroke patients. While previous research has focused primarily on lesion location or size, with mixed or inconclusive results, we targeted the extended structural network implicated by the lesion, and evaluated the added explanatory value of a structural disconnectivity approach with regards to the brain correlates of PSF. To this end, we estimated individual structural brain disconnectome maps in 84 S survivors in the chronic phase (≥3 months post stroke) using information about lesion location and normative white matter pathways obtained from 170 healthy individuals. PSF was measured by the Fatigue Severity Scale (FSS). Voxel wise analyses using non-parametric permutation-based inference were conducted on disconnectome maps to estimate regional effects of disconnectivity. Associations between PSF and global disconnectivity and clinical lesion characteristics were tested by linear models, and we estimated Bayes factor to quantify the evidence for the null and alternative hypotheses, respectively. The results revealed no significant associations between PSF and disconnectome measures or lesion characteristics, with moderate evidence in favor of the null hypothesis. These results suggest that symptoms of post-stroke fatigue among chronic stroke patients are not simply explained by lesion characteristics or the extent and distribution of structural brain disconnectome, and are discussed in light of methodological considerations.

Comparison of diffusion MRI and CLARITY fiber orientation estimates in both gray and white matter regions of human and primate brain

Diffusion MRI (dMRI) represents one of the few methods for mapping brain fiber orientations non-invasively. Unfortunately, dMRI fiber mapping is an indirect method that relies on inference from measured diffusion patterns. Comparing dMRI results with other modalities is a way to improve the interpretation of dMRI data and help advance dMRI technologies. Here, we present methods for comparing dMRI fiber orientation estimates with optical imaging of fluorescently labeled neurofilaments and vasculature in 3D human and primate brain tissue cuboids cleared using CLARITY. The recent advancements in tissue clearing provide a new opportunity to histologically map fibers projecting in 3D, which represents a captivating complement to dMRI measurements. In this work, we demonstrate the capability to directly compare dMRI and CLARITY in the same human brain tissue and assess multiple approaches for extracting fiber orientation estimates from CLARITY data. We estimate the three-dimensional neuronal fiber and vasculature orientations from neurofilament and vasculature stained CLARITY images by calculating the tertiary eigenvector of structure tensors. We then extend CLARITY orientation estimates to an orientation distribution function (ODF) formalism by summing multiple sub-voxel structure tensor orientation estimates. In a sample containing part of the human thalamus, there is a mean angular difference of 19o±15o between the primary eigenvectors of the dMRI tensors and the tertiary eigenvectors from the CLARITY neurofilament stain. We also demonstrate evidence that vascular compartments do not affect the dMRI orientation estimates by showing an apparent lack of correspondence (mean angular difference = 49o±23o) between the orientation of the dMRI tensors and the structure tensors in the vasculature stained CLARITY images. In a macaque brain dataset, we examine how the CLARITY feature extraction depends on the chosen feature extraction parameters. By varying the volume of tissue over which the structure tensor estimates are derived, we show that orientation estimates are noisier with more spurious ODF peaks for sub-voxels below 30 µm3 and that, for our data, the optimal gray matter sub-voxel size is between 62.5 µm3 and 125 µm3. The example experiments presented here represent an important advancement towards robust multi-modal MRI-CLARITY comparisons.

A comparison of diffusion tractography techniques in simulating the generalized Ising model to predict the intrinsic activity of the brain

Diffusion tractography is a non-invasive technique that is being used to estimate the location and direction of white matter tracts in the brain. Identifying the characteristics of white matter plays an important role in research as well as in clinical practice that relies on finding the relationship between the structure and function of the brain. An Ising model implemented on a structural connectivity (SC) has proven to explain the spontaneous fluctuations in the brain at criticality using brain's structure depicted by white matter tracts. Since the SC is the only input of the model, identifying the tractography technique which provides a SC that delivers the highest prediction of the brain's intrinsic activity via the generalized Ising model (GIM) is essential. Hence an Ising model is simulated on SCs generated using two different acquisition schemes (single and multi-shell) and two different tractography approaches (deterministic and probabilistic) and analyzed at criticality across 69 healthy subjects. Results showed that by introducing the GIM, predictability of the empirical correlation matrix increases on average from 0.2 to 0.6 compared to the predictability using the empirical connectivity matrix directly. It is also observed that the SC generated using deterministic tractography without fractional anisotropy resulted in the highest correlation coefficient value of 0.65 between the simulated and empirical correlation matrices. Additionally, calculated dimensionalities per simulation illustrated that the dimensionality depends upon the method of tractography that has been used to extract the SC.

Variability of white matter anatomy in the subcallosal cingulate area

The subcallosal cingulate (SCC) area is a putative hub in the brain network underlying depression. Deep brain stimulation (DBS) targeting a particular subregion of SCC, identified as the intersection of forceps minor (FM), uncinate fasciculus (UCF), cingulum and fronto-striatal fiber bundles, may be critical to a therapeutic response in patients with severe, treatment-resistant forms of major depressive disorder (MDD). The pattern and variability of the white matter anatomy and organization within SCC has not been extensively characterized across individuals. The goal of this study is to investigate the variability of white matter bundles within the SCC that structurally connect this region with critical nodes in the depression network. Structural and diffusion data from 100 healthy subjects from the Human Connectome Project database were analyzed. Anatomically defined SCC regions were used as seeds to perform probabilistic tractography and to estimate the connectivity from the SCC to subject-specific target areas believed to be involved in the pathology of MDD including ventral striatum (VS), UCF, anterior cingulate cortex (ACC), and medial prefrontal cortex (mPFC). Four distinct areas of connectivity were identified within SCC across subjects: (a) postero-lateral SCC connectivity to medial temporal regions via UCF, (b) postero-medial connectivity to VS, (c) superior-medial connectivity to ACC via cingulum bundle, and (d) antero-lateral connectivity to mPFC regions via forceps minor. Assuming white matter connectivity is critical to therapeutic response, the improved anatomic understanding of SCC as well as an appreciation of the intersubject variability are critical to developing optimized therapeutic targeting for SCC DBS.

Two cortical representations of voice control are differentially involved in speech fluency

Recent studies have identified two distinct cortical representations of voice control in humans, the ventral and the dorsal laryngeal motor cortex. Strikingly, while persistent developmental stuttering has been linked to a white-matter deficit in the ventral laryngeal motor cortex, intensive fluency-shaping intervention modulated the functional connectivity of the dorsal laryngeal motor cortical network. Currently, it is unknown whether the underlying structural network organization of these two laryngeal representations is distinct or differently shaped by stuttering intervention. Using probabilistic diffusion tractography in 22 individuals who stutter and participated in a fluency shaping intervention, in 18 individuals who stutter and did not participate in the intervention and in 28 control participants, we here compare structural networks of the dorsal laryngeal motor cortex and the ventral laryngeal motor cortex and test intervention-related white-matter changes. We show (i) that all participants have weaker ventral laryngeal motor cortex connections compared to the dorsal laryngeal motor cortex network, regardless of speech fluency, (ii) connections of the ventral laryngeal motor cortex were stronger in fluent speakers, (iii) the connectivity profile of the ventral laryngeal motor cortex predicted stuttering severity (iv) but the ventral laryngeal motor cortex network is resistant to a fluency shaping intervention. Our findings substantiate a weaker structural organization of the ventral laryngeal motor cortical network in developmental stuttering and imply that assisted recovery supports neural compensation rather than normalization. Moreover, the resulting dissociation provides evidence for functionally segregated roles of the ventral laryngeal motor cortical and dorsal laryngeal motor cortical networks.

Clinical Evaluation of Cingulum Bundle Connectivity for Neurosurgical Hypothesis Development

BACKGROUND

The cingulum bundle (CB) has long been a target for psychiatric neurosurgical procedures, but with limited understanding of the brain networks being impacted. Recent advances in human tractography could provide a foundation to better understand the effects of neurosurgical interventions on the CB; however, the reliability of tractography remains in question.

OBJECTIVE

To evaluate the ability of different tractography techniques, derived from typical, human diffusion-weighted imaging (DWI) data, to characterize CB connectivity described in animal models. This will help validate the clinical applicability of tractography, and generate insight on current and future neurosurgical targets for psychiatric disorders.

METHODS

Connectivity of the CB in 15 healthy human subjects was evaluated using DWI-based tractography, and compared to tract-tracing findings from nonhuman primates. Brain regions of interest were defined to coincide with the animal model. Tractography was performed using 3 techniques (FSL probabilistic, Camino probabilistic, and Camino deterministic). Differences in connectivity were assessed, and the CB segment with the greatest connectivity was determined.

RESULTS

Each tractography technique successfully reproduced the animal tracing model with a mean accuracy of 72% (68-75%, P < .05). Additionally, one region of the CB, the rostral dorsal segment, had significantly greater connectivity to associated brain structures than all other CB segments (P < .05).

CONCLUSION

Noninvasive, in vivo human analysis of the CB, using clinically available DWI for tractography, consistently reproduced the results of an animal tract-tracing model. This suggests that tractography of the CB can be used for clinical applications, which may aid in neurosurgical targeting for psychiatric disorders.

Reproducing macaque lateral grasping and oculomotor networks using resting state functional connectivity and diffusion tractography

Cortico-cortical networks involved in motor control have been well defined in the macaque using a range of invasive techniques. The advent of neuroimaging has enabled non-invasive study of these large-scale functionally specialized networks in the human brain; however, assessing its accuracy in reproducing genuine anatomy is more challenging. We set out to assess the similarities and differences between connections of macaque motor control networks defined using axonal tracing and those reproduced using structural and functional connectivity techniques. We processed a cohort of macaques scanned in vivo that were made available by the open access PRIME-DE resource, to evaluate connectivity using diffusion imaging tractography and resting state functional connectivity (rs-FC). Sectors of the lateral grasping and exploratory oculomotor networks were defined anatomically on structural images, and connections were reproduced using different structural and functional approaches (probabilistic and deterministic whole-brain and seed-based tractography; group template and native space functional connectivity analysis). The results showed that parieto-frontal connections were best reproduced using both structural and functional connectivity techniques. Tractography showed lower sensitivity but better specificity in reproducing connections identified by tracer data. Functional connectivity analysis performed in native space had higher sensitivity but lower specificity and was better at identifying connections between intrasulcal ROIs than group-level analysis. Connections of AIP were most consistently reproduced, although those connected with prefrontal sectors were not identified. We finally compared diffusion MR modelling with histology based on an injection in AIP and speculate on anatomical bases for the observed false negatives. Our results highlight the utility of precise ex vivo techniques to support the accuracy of neuroimaging in reproducing connections, which is relevant also for human studies.

Diffusion Tensor Imaging-Based Studies at the Group-Level Applied to Animal Models of Neurodegenerative Diseases

The understanding of human and non-human microstructural brain alterations in the course of neurodegenerative diseases has substantially improved by the non-invasive magnetic resonance imaging (MRI) technique of diffusion tensor imaging (DTI). Animal models (including disease or knockout models) allow for a variety of experimental manipulations, which are not applicable to humans. Thus, the DTI approach provides a promising tool for cross-species cross-sectional and longitudinal investigations of the neurobiological targets and mechanisms of neurodegeneration. This overview with a systematic review focuses on the principles of DTI analysis as used in studies at the group level in living preclinical models of neurodegeneration. The translational aspect from in-vivo animal models toward (clinical) applications in humans is covered as well as the DTI-based research of the non-human brains' microstructure, the methodological aspects in data processing and analysis, and data interpretation at different abstraction levels. The aim of integrating DTI in multiparametric or multimodal imaging protocols will allow the interrogation of DTI data in terms of directional flow of information and may identify the microstructural underpinnings of neurodegeneration-related patterns.

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.

Mapping Short Association Fibers in the Early Cortical Visual Processing Stream Using In Vivo Diffusion Tractography

Short association fibers (U-fibers) connect proximal cortical areas and constitute the majority of white matter connections in the human brain. U-fibers play an important role in brain development, function, and pathology but are underrepresented in current descriptions of the human brain connectome, primarily due to methodological challenges in diffusion magnetic resonance imaging (dMRI) of these fibers. High spatial resolution and dedicated fiber and tractography models are required to reliably map the U-fibers. Moreover, limited quantitative knowledge of their geometry and distribution makes validation of U-fiber tractography challenging. Submillimeter resolution diffusion MRI—facilitated by a cutting-edge MRI scanner with 300 mT/m maximum gradient amplitude—was used to map U-fiber connectivity between primary and secondary visual cortical areas (V1 and V2, respectively) in vivo. V1 and V2 retinotopic maps were obtained using functional MRI at 7T. The mapped V1–V2 connectivity was retinotopically organized, demonstrating higher connectivity for retinotopically corresponding areas in V1 and V2 as expected. The results were highly reproducible, as demonstrated by repeated measurements in the same participants and by an independent replication group study. This study demonstrates a robust U-fiber connectivity mapping in vivo and is an important step toward construction of a more complete human brain connectome.

The emergence of dyslexia in the developing brain

Developmental dyslexia, a severe deficit in literacy learning, is a neurodevelopmental learning disorder. Yet, it is not clear whether existing neurobiological accounts of dyslexia capture potential predispositions of the deficit or consequences of reduced reading experience. Here, we longitudinally followed 32 children from preliterate to school age using functional and structural magnetic resonance imaging techniques. Based on standardised and age-normed reading and spelling tests administered at school age, children were classified as 16 dyslexic participants and 16 controls. This longitudinal design allowed us to disentangle possible neurobiological predispositions for developing dyslexia from effects of individual differences in literacy experience. In our sample, the disorder can be predicted already before literacy learning from auditory cortex gyrification and aberrant downstream connectivity within the speech processing system. These results provide evidence for the notion that dyslexia may originate from an atypical maturation of the speech network that precedes literacy instruction.

Tractography reproducibility challenge with empirical data (TraCED): The 2017 ISMRM diffusion study group challenge

BACKGROUND

Fiber tracking with diffusion-weighted MRI has become an essential tool for estimating in vivo brain white matter architecture. Fiber tracking results are sensitive to the choice of processing method and tracking criteria.

PURPOSE

To assess the variability for an algorithm in group studies reproducibility is of critical context. However, reproducibility does not assess the validity of the brain connections. Phantom studies provide concrete quantitative comparisons of methods relative to absolute ground truths, yet do no capture variabilities because of in vivo physiological factors. The ISMRM 2017 TraCED challenge was created to fulfill the gap.

STUDY TYPE

A systematic review of algorithms and tract reproducibility studies.

SUBJECTS

Single healthy volunteers.

FIELD STRENGTH/SEQUENCE

3.0T, two different scanners by the same manufacturer. The multishell acquisition included b-values of 1000, 2000, and 3000 s/mm2 with 20, 45, and 64 diffusion gradient directions per shell, respectively.

ASSESSMENT

Nine international groups submitted 46 tractography algorithm entries each consisting 16 tracts per scan. The algorithms were assessed using intraclass correlation (ICC) and the Dice similarity measure.

STATISTICAL TESTS

Containment analysis was performed to assess if the submitted algorithms had containment within tracts of larger volume submissions. This also serves the purpose to detect if spurious submissions had been made.

RESULTS

The top five submissions had high ICC and Dice >0.88. Reproducibility was high within the top five submissions when assessed across sessions or across scanners: 0.87-0.97. Containment analysis shows that the top five submissions are contained within larger volume submissions. From the total of 16 tracts as an outcome relatively the number of tracts with high, moderate, and low reproducibility were 8, 4, and 4.

DATA CONCLUSION

The different methods clearly result in fundamentally different tract structures at the more conservative specificity choices. Data and challenge infrastructure remain available for continued analysis and provide a platform for comparison.

LEVEL OF EVIDENCE

5 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:234-249.

Acute ex vivo changes in brain white matter diffusion tensor metrics

Purpose

Diffusion magnetic resonance imaging and tractography has an important role in the visualization of brain white matter and assessment of tissue microstructure. There is a lack of correspondence between diffusion metrics of live tissue, ex vivo tissue, and histological findings. The objective of this study is to elucidate this connection by determining the specific diffusion alterations between live and ex vivo brain tissue. This may have an important role in the incorporation of diffusion imaging in ex vivo studies as a complement to histological sectioning as well as investigations of novel neurosurgical techniques.

Methods

This study presents a method of high angular resolution diffusion imaging and tractography of intact and non-fixed ex vivo piglet brains. Most studies involving ex vivo brain specimens have been formalin-fixed or excised from their original biological environment, processes both of which are known to affect diffusion parameters. Thus, non-fixed ex vivo tissue is used. A region-of-interest based analysis of diffusion tensor metrics are compared to in vivo subjects in a selection of major white matter bundles in order to assess the translatability of ex vivo diffusion measurements.

Results

Tractography was successfully achieved in both in vivo and ex vivo groups. No significant differences were found in tract connectivity, average streamline length, or apparent fiber density. Significantly decreased diffusivity (mean, axial, and radial; p<0.0005) in the non-fixed ex vivo group and unaltered fractional anisotropy (p>0.059) between groups were observed.

Conclusion

This study validates the extrapolation of non-fixed fractional anisotropy measurements to live tissue and the potential use of ex vivo tissue for methodological development.

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

Evaluation of the structural connectivity (SC) of the brain based on tractography has mainly focused on the choice of diffusion model, tractography algorithm, and their respective parameter settings. Here, we systematically validate SC derived from a post mortem monkey brain, while varying key acquisition parameters such as the b-value, gradient angular resolution and image resolution. As gold standard we use the connectivity matrix obtained invasively with histological tracers by Markov et al. (2014). As performance metric, we use cross entropy as a measure that enables comparison of the relative tracer labeled neuron counts to the streamline counts from tractography. We find that high angular resolution and high signal-to-noise ratio are important to estimate SC, and that SC derived from low image resolution (1.0³ mm³) are in better agreement with the tracer network, than those derived from high image resolution (0.5³ mm³) or at an even lower image resolution (2.0³ mm³). In contradiction, sensitivity and specificity analyses suggest that if the angular resolution is sufficient, the balanced compromise in which sensitivity and specificity are identical remains 60-64% regardless of the other scanning parameters. Interestingly, the tracer graph is assumed to be the gold standard but by thresholding, the balanced compromise increases to 70-75%. Hence, by using performance metrics based on binarized tracer graphs, one risks losing important information, changing the performance of SC graphs derived by tractography and their dependence of different scanning parameters.

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.

A Web-Based Atlas Combining MRI and Histology of the Squirrel Monkey Brain

The squirrel monkey (Saimiri sciureus) is a commonly-used surrogate for humans in biomedical research. In the neuroimaging community, MRI and histological atlases serve as valuable resources for anatomical, physiological, and functional studies of the brain; however, no digital MRI/histology atlas is currently available for the squirrel monkey. This paper describes the construction of a web-based multi-modal atlas of the squirrel monkey brain. The MRI-derived information includes anatomical MRI contrast (i.e., T2-weighted and proton-density-weighted) and diffusion MRI metrics (i.e., fractional anisotropy and mean diffusivity) from data acquired both in vivo and ex vivo on a 9.4 Tesla scanner. The histological images include Nissl and myelin stains, co-registered to the corresponding MRI, allowing identification of cyto- and myelo-architecture. In addition, a bidirectional neuronal tracer, biotinylated dextran amine (BDA) was injected into the primary motor cortex, enabling highly specific identification of regions connected to the injection location. The atlas integrates the results of common image analysis methods including diffusion tensor imaging glyphs, labels of 57 white-matter tracts identified using DTI-tractography, and 18 cortical regions of interest identified from Nissl-revealed cyto-architecture. All data are presented in a common space, and all image types are accessible through a web-based atlas viewer, which allows visualization and interaction of user-selectable contrasts and varying resolutions. By providing an easy to use reference system of anatomical information, our web-accessible multi-contrast atlas forms a rich and convenient resource for comparisons of brain findings across subjects or modalities. The atlas is called the Combined Histology-MRI Integrated Atlas of the Squirrel Monkey (CHIASM). All images are accessible through our web-based viewer ( https://chiasm.vuse.vanderbilt.edu /), and data are available for download at ( https://www.nitrc.org/projects/smatlas/ ).

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.

Fingerprinting Orientation Distribution Functions in diffusion MRI detects smaller crossing angles

Diffusion tractography is routinely used to study white matter architecture and brain connectivity in vivo. A key step for successful tractography of neuronal tracts is the correct identification of tract directions in each voxel. Here we propose a fingerprinting-based methodology to identify these fiber directions in Orientation Distribution Functions, dubbed ODF-Fingerprinting (ODF-FP). In ODF-FP, fiber configurations are selected based on the similarity between measured ODFs and elements in a pre-computed library. In noisy ODFs, the library matching algorithm penalizes the more complex fiber configurations. ODF simulations and analysis of bootstrapped partial and whole-brain in vivo datasets show that the ODF-FP approach improves the detection of fiber pairs with small crossing angles while maintaining fiber direction precision, which leads to better tractography results. Rather than focusing on the ODF maxima, the ODF-FP approach uses the whole ODF shape to infer fiber directions to improve the detection of fiber bundles with small crossing angle. The resulting fiber directions aid tractography algorithms in accurately displaying neuronal tracts and calculating brain connectivity.