Parallel Transport Tractography

Log-Cholesky filtering of diffusion tensor fields: Impact on noise reduction

Diffusion tensor imaging (DTI) is a powerful neuroimaging technique that provides valuable insights into the microstructure and connectivity of the brain. By measuring the diffusion of water molecules along neuronal fibers, DTI allows the visualization and study of intricate networks of neural pathways. DTI is a noise-sensitive method, where a low signal-to-noise ratio (SNR) results in significant errors in the estimated tensor field. Tensor field regularization is an effective solution for noise reduction. Diffusion tensors are represented by symmetric positive-definite (SPD) matrices. The space of SPD matrices may be viewed as a Riemannian manifold after defining a suitable metric on its tangent bundle. The Log-Cholesky metric is a recently developed concept with advantages over previously defined Riemannian metrics, such as the affine-invariant and Log-Euclidean metrics. The utility of the Log-Cholesky metric for tensor field regularization and noise reduction has not been investigated in detail. This manuscript provides a quantitative investigation of the impact of Log-Cholesky filtering on noise reduction in DTI. It also provides sufficient details of the linear algebra and abstract differential geometry concepts necessary to implement this technique as a simple and effective solution to filtering diffusion tensor fields.

Preoperative diffusion tensor imaging: Fiber-trajectory-distribution-based tractography to identify facial nerve in vestibular schwannoma

PURPOSE

Tractography of the facial nerve based on diffusion MRI is instrumental before surgery for the resection of vestibular schwannoma, but no excellent methods usable for the suppression of motion and image noise have been proposed. The aim of this study was to effectively suppress noise and provide accurate facial nerve reconstruction by extend a fiber trajectory distribution function based on the fourth-order streamline differential equations.

METHODS

Preoperative MRI from 33 patients with vestibular schwannoma who underwent surgical resection were utilized in this study. First, T1WI and T2WI were used to obtain mask images and regions of interest. Second, probabilistic tractography was employed to obtain the fibers representing the approximate facial nerve pathway, and these fibers were subsequently translated into orientation information for each voxel. Last, the voxel orientation information and the peaks of the fiber orientation distribution were combined to generate a fiber trajectory distribution function, which was used to parameterize the anatomical information. The parameters were determined by minimizing the cost between the trajectory of fibers and the estimated directions.

RESULTS

Qualitative and visual analyses were used to compare facial nerve reconstruction with intraoperative recordings. Compared with other methods (SD_Stream, iFOD1, iFOD2, unscented Kalman filter, parallel transport tractography), the fiber-trajectory-distribution-based tractography provided the most accurate facial nerve reconstructions.

CONCLUSION

The fiber-trajectory-distribution-based tractography can effectively suppress the effect of noise. It is a more valuable aid for surgeons before vestibular schwannoma resection, which may ultimately improve the postsurgical patient's outcome.

Altered interhemispheric connectivity in Huntington's Disease

Pyramidal cells give rise to the corpus callosum, interhemispheric fibers that constitute the associations between the left and the right hemispheres. These interconnections are the substrates of important neurological functions, such as perception, memory, emotion, and movement control, which are all affected in Huntington's disease (HD). In this study we used directional tract density patterns (dTDPs) to evaluate changes in interhemispheric connectivity in gene-expanded individuals, which included presymptomatic and early symptomatic HD subjects. Our results demonstrated regionally selective and progressive differences in dTDPs between distinct regions of the corpus callosum (subdivided by Hofer-Frahm scheme) in the gene-expanded cohorts. In the presymptomatic HD cohort, we found trends, such that the density of fibers was reduced in CC regions IIb, III, and IV (p < 0.05); fibers from these regions project to sensory, premotor, and motor cortical regions, respectively. In the HD cohort, we found reduction in the density of fibers in all CC regions, including in fibers extending to the cortical surface (p < 0.002). Our results support the use of dTDPs to evaluate individual and progressive changes in interhemispheric connectivity in HD.

Bundle-specific tractogram distribution estimation using higher-order streamline differential equation

Streamline tractography locally traces peak directions extracted from fiber orientation distribution (FOD) functions, lacking global information about the trend of the whole fiber bundle. Therefore, it is prone to producing erroneous tracks while missing true positive connections. In this work, we propose a new bundle-specific tractography (BST) method based on a bundle-specific tractogram distribution (BTD) function, which directly reconstructs the fiber trajectory from the start region to the termination region by incorporating the global information in the fiber bundle mask. A unified framework for any higher-order streamline differential equation is presented to describe the fiber bundles with disjoint streamlines defined based on the diffusion vectorial field. At the global level, the tractography process is simplified as the estimation of BTD coefficients by minimizing the energy optimization model, and is used to characterize the relations between BTD and diffusion tensor vector under the prior guidance by introducing the tractogram bundle information to provide anatomic priors. Experiments are performed on simulated Hough, Sine, Circle data, ISMRM 2015 Tractography Challenge data, FiberCup data, and in vivo data from the Human Connectome Project (HCP) for qualitative and quantitative evaluation. Results demonstrate that our approach reconstructs complex fiber geometry more accurately. BTD reduces the error deviation and accumulation at the local level and shows better results in reconstructing long-range, twisting, and large fanning tracts.

Tractometry of Human Visual White Matter Pathways in Health and Disease

Diffusion-weighted MRI (dMRI) provides a unique non-invasive view of human brain tissue properties. The present review article focuses on tractometry analysis methods that use dMRI to assess the properties of brain tissue within the long-range connections comprising brain networks. We focus specifically on the major white matter tracts that convey visual information. These connections are particularly important because vision provides rich information from the environment that supports a large range of daily life activities. Many of the diseases of the visual system are associated with advanced aging, and tractometry of the visual system is particularly important in the modern aging society. We provide an overview of the tractometry analysis pipeline, which includes a primer on dMRI data acquisition, voxelwise model fitting, tractography, recognition of white matter tracts, and calculation of tract tissue property profiles. We then review dMRI-based methods for analyzing visual white matter tracts: the optic nerve, optic tract, optic radiation, forceps major, and vertical occipital fasciculus. For each tract, we review background anatomical knowledge together with recent findings in tractometry studies on these tracts and their properties in relation to visual function and disease. Overall, we find that measurements of the brain's visual white matter are sensitive to a range of disorders and correlate with perceptual abilities. We highlight new and promising analysis methods, as well as some of the current barriers to progress toward integration of these methods into clinical practice. These barriers, such as variability in measurements between protocols and instruments, are targets for future development.

Diffusion Model-based FOD Restoration from High Distortion in dMRI

Fiber orientation distributions (FODs) is a popular model to represent the diffusion MRI (dMRI) data. However, imaging artifacts such as susceptibility-induced distortion in dMRI can cause signal loss and lead to the corrupted reconstruction of FODs, which prohibits successful fiber tracking and connectivity analysis in affected brain regions such as the brain stem. Generative models, such as the diffusion models, have been successfully applied in various image restoration tasks. However, their application on FOD images poses unique challenges since FODs are 4-dimensional data represented by spherical harmonics (SPHARM) with the 4-th dimension exhibiting order-related dependency. In this paper, we propose a novel diffusion model for FOD restoration that can recover the signal loss caused by distortion artifacts. We use volume-order encoding to enhance the ability of the diffusion model to generate individual FOD volumes at all SPHARM orders. Moreover, we add cross-attention features extracted across all SPHARM orders in generating every individual FOD volume to capture the order-related dependency across FOD volumes. We also condition the diffusion model with low-distortion FODs surrounding high-distortion areas to maintain the geometric coherence of the generated FODs. We trained and tested our model using data from the UK Biobank (n = 1315). On a test set with ground truth (n = 43), we demonstrate the high accuracy of the generated FODs in terms of root mean square errors of FOD volumes and angular errors of FOD peaks. We also apply our method to a test set with large distortion in the brain stem area (n = 1172) and demonstrate the efficacy of our method in restoring the FOD integrity and, hence, greatly improving tractography performance in affected brain regions.

Optimal shrinkage denoising breaks the noise floor in high-resolution diffusion MRI

The spatial resolution attainable in diffusion magnetic resonance (MR) imaging is inherently limited by noise. The weaker signal associated with a smaller voxel size, especially at a high level of diffusion sensitization, is often buried under the noise floor owing to the non-Gaussian nature of the MR magnitude signal. Here, we show how the noise floor can be suppressed remarkably via optimal shrinkage of singular values associated with noise in complex-valued k-space data from multiple receiver channels. We explore and compare different low-rank signal matrix recovery strategies to utilize the inherently redundant information from multiple channels. In combination with background phase removal, the optimal strategy reduces the noise floor by 11 times. Our framework enables imaging with substantially improved resolution for precise characterization of tissue microstructure and white matter pathways without relying on expensive hardware upgrades and time-consuming acquisition repetitions, outperforming other related denoising methods.

Blurred streamlines: A novel representation to reduce redundancy in tractography

Tractography is a powerful tool to study brain connectivity in vivo, but it is well known to suffer from an intrinsic trade-off between sensitivity and specificity. A critical - but usually underrated - parameter to choose that can heavily impact the quality of the estimates is the number of streamlines to be reconstructed for a given data set. In fact, sensitivity can be improved by generating more and more streamlines, as all real anatomical connections are likely reconstructed, but lots of false positives are inevitably introduced, too. Consequently, so-called tractography filtering techniques have become increasingly popular to get rid of these false positives and improve specificity. However, increasing number of streamlines introduces redundancy in tractography reconstructions, which may negatively impact the performance of filtering algorithms, especially those based on linear formulations. To address this problem, we introduce a novel streamlines representation, called "blurred streamlines", which drastically reduces the redundancy among streamlines by (i) clustering similar trajectories and (ii) spatially blurring the corresponding signal contributions. We tested the effectiveness of the blurred streamlines both on synthetic and in vivo data. Our results clearly show that this new representation is as accurate as state-of-the-art methods despite using only 5% of the input streamlines, thus significantly decreasing the computational complexity of filtering algorithms as well as storage requirements of the resulting reconstructions.

A unified filtering method for estimating asymmetric orientation distribution functions

Numerous filtering methods have been proposed for estimating asymmetric orientation distribution functions (ODFs) for diffusion magnetic resonance imaging (dMRI). It can be hard to make sense of all these different methods, which share similar features and result in similar outputs. In this work, we disentangle these many filtering methods proposed in the past and combine them into a novel, unified filtering equation. We also propose a self-supervised data-driven approach for calibrating the filtering parameter values. Our equation is implemented in an open-source GPU-accelerated python software to facilitate its integration into any existing dMRI processing pipeline. Our method is applied on multi-shell multi-tissue fiber ODFs from the Human Connectome Project dataset (1.25 mm3 native resolution) and on single-shell single-tissue fiber ODFs from the Bilingualism and the Brain dataset (2.0 mm3 isotropic resolution) to evaluate the occurrence of asymmetric patterns on different spatial resolutions, representing cutting-edge and "clinical" research data. Asymmetry measures such as the asymmetric index (ASI) and our novel number of fiber directions (NuFiD) are then used to explain the behaviour of our method in these images. The contributions of this work are: (i) the disentanglement and unification of filtering methods for estimating asymmetric ODFs; (ii) a calibration method for automatically fixing the parameters governing the filtering; (iii) an open-source, efficient implementation of our unified filtering method for estimating asymmetric ODFs; (iv) a novel number of fiber directions (NuFiD) index for explaining asymmetric fiber configurations; and (v) a novel template of asymmetries, revealing that our filtering method estimates asymmetric configurations in at least 50% of the brain voxels (∼31% of the white matter and ∼63% of the gray matter).

Neurite orientation dispersion and density imaging quantifies microstructural impairment in the thalamus and its connectivity in amyotrophic lateral sclerosis

AIMS

To evaluate microstructural impairment in the thalamus and thalamocortical connectivity using neurite orientation dispersion and density imaging (NODDI) in amyotrophic lateral sclerosis (ALS).

METHODS

This study included 47 healthy controls and 43 ALS patients, whose structural and diffusion-weighted data were collected. We used state-of-the-art parallel transport tractography to identify thalamocortical pathways in individual spaces. Thalamus was then parcellated into six subregions based on its connectivity pattern with the priori defined cortical (i.e., prefrontal/motor/somatosensory/temporal/posterior-parietal/occipital) regions. For each of the thalamic and cortical subregions and thalamo-cortical tracts, we compared the following NODDI metrics between groups: orientation dispersion index (ODI), neurite density index (NDI), and isotropic volume fraction (ISO). We also used these metrics to conduct receiver operating characteristic curve (ROC) analyses and Spearman correlation.

RESULTS

In ALS patients, we found decreased ODI and increased ISO in the thalamic subregion connecting the left motor cortex and other extramotor (e.g., somatosensory and occipital) cortex (Bonferroni-corrected p < 0.05). NDI decreased in the bilateral thalamo-motor and thalamo-somatosensory tracts and in the right thalamo-posterior-parietal and thalamo-occipital tracts (Bonferroni-corrected p < 0.05). NDI reduction in the bilateral thalamo-motor tract (p = 0.017 and 0.009) and left thalamo-somatosensory tract (p = 0.029) was correlated with disease severity. In thalamo-cortical tracts, NDI yielded a higher effect size during between-group comparisons and a greater area under ROC (p < 0.05) compared with conventional diffusion tensor imaging metrics.

CONCLUSIONS

Microstructural impairment in the thalamus and thalamocortical connectivity is the hallmark of ALS. NODDI improved the detection of disrupted thalamo-cortical connectivity in ALS.

Retinotopic degeneration of the retina and optic tracts in autosomal dominant Alzheimer's disease

INTRODUCTION

We investigated the correlation between retinal thickness and optic tract integrity in subjects with autosomal dominant Alzheimer's disease (ADAD) causing mutations.

METHODS

Retinal thicknesses and diffusion tensor images (DTI) were obtained using optical coherence tomography and magnetic resonance imaging, respectively. The association between retinal thickness and DTI measures was adjusted for age, sex, retinotopy, and correlation between eyes.

RESULTS

Optic tract mean diffusivity and axial diffusivity were negatively correlated with retinotopically defined ganglion cell inner plexiform thickness (GCIPL). Fractional anisotropy was negatively correlated with retinotopically defined retinal nerve fiber layer thickness. There was no correlation between outer nuclear layer (ONL) thickness and any DTI measure.

DISCUSSION

In ADAD, GCIPL thickness is significantly associated with retinotopic optic tract DTI measures even in minimally symptomatic subjects. Similar associations were not present with ONL thickness or when ignoring retinotopy. We provide in vivo evidence for optic tract changes resulting from ganglion cell pathology in ADAD.

Evaluation of tractogram filtering methods using human-like connectome phantoms

Tractography algorithms are prone to reconstructing spurious connections. The set of streamlines generated with tractography can be post-processed to retain the streamlines that are most biologically plausible. Several microstructure-informed filtering algorithms are available for this purpose, however, the comparative performance of these methods has not been extensively evaluated. In this study, we aim to evaluate streamline filtering and post-processing algorithms using simulated connectome phantoms. We first establish a framework for generating connectome phantoms featuring brain-like white matter fiber architectures. We then use our phantoms to systematically evaluate the performance of a range of streamline filtering algorithms, including SIFT, COMMIT, and LiFE. We find that all filtering methods successfully improve connectome accuracy, although filter performance depends on the complexity of the underlying white matter fiber architecture. Filtering algorithms can markedly improve tractography accuracy for simple tubular fiber bundles (F-measure deterministic- unfiltered: 0.49 and best filter: 0.72; F-measure probabilistic- unfiltered: 0.37 and best filter: 0.81), but for more complex brain-like fiber architectures, the improvement is modest (F-measure deterministic- unfiltered: 0.53 and best filter: 0.54; F-measure probabilistic- unfiltered: 0.46 and best filter: 0.50). Overall, filtering algorithms have the potential to improve the accuracy of connectome mapping pipelines, particularly for weighted connectomes and pipelines using probabilistic tractography methods. Our results highlight the need for further advances tractography and streamline filtering to improve the accuracy of connectome mapping.

FASSt : Filtering via Symmetric Autoencoder for Spherical Superficial White Matter Tractography

Superficial white matter (SWM) plays an important role in functioning of the human brain, and it contains a large amount of cortico-cortical connections. However, the difficulties of generating complete and reliable U-fibers make SWM-related analysis lag behind relatively matured Deep white matter (DWM) analysis. With the aid of some newly proposed surface-based SWM tractography algorithms, we have developed a specialized SWM filtering method based on a symmetric variational autoencoder (VAE). In this work, we first demonstrate the advantage of the spherical representation and generate these spherical tracts using the triangular mesh and the registered spherical surface. We then introduce the Filtering via symmetric Autoencoder for Spherical Superficial White Matter tractography (FASSt) framework with a novel symmetric weights module to perform the filtering task in a latent space. We evaluate and compare our method with the state-of-the-art clustering-based method on diffusion MRI data from Human Connectome Project (HCP). The results show that our proposed method outperform these clustering methods and achieves excellent performance in groupwise consistency and topographic regularity.

Flow-based Geometric Interpolation of Fiber Orientation Distribution Functions

The fiber orientation distribution function (FOD) is an advanced model for high angular resolution diffusion MRI representing complex fiber geometry. However, the complicated mathematical structures of the FOD function pose challenges for FOD image processing tasks such as interpolation, which plays a critical role in the propagation of fiber tracts in tractography. In FOD-based tractography, linear interpolation is commonly used for numerical efficiency, but it is prone to generate false artificial information, leading to anatomically incorrect fiber tracts. To overcome this difficulty, we propose a flowbased and geometrically consistent interpolation framework that considers peak-wise rotations of FODs within the neighborhood of each location. Our method decomposes a FOD function into multiple components and uses a smooth vector field to model the flows of each peak in its neighborhood. To generate the interpolated result along the flow of each vector field, we develop a closed-form and efficient method to rotate FOD peaks in neighboring voxels and realize geometrically consistent interpolation of FOD components. By combining the interpolation results from each peak, we obtain the final interpolation of FODs. Experimental results on Human Connectome Project (HCP) data demonstrate that our method produces anatomically more meaningful FOD interpolations and significantly enhances tractography performance.

Limb-specific thalamocortical tracts are impaired differently in hemiplegic and diplegic subtypes of cerebral palsy

Thalamocortical pathways are considered crucial in the sensorimotor functioning of children with cerebral palsy (CP). However, previous research has been limited by non-specific tractography seeding and the lack of comparison between different CP subtypes. We compared limb-specific thalamocortical tracts between children with hemiplegic (HP, N = 15) or diplegic (DP, N = 10) CP and typically developed peers (N = 19). The cortical seed-points for the upper and lower extremities were selected (i) manually based on anatomical landmarks or (ii) using functional magnetic resonance imaging (fMRI) activations following proprioceptive-limb stimulation. Correlations were investigated between tract structure (mean diffusivity, MD; fractional anisotropy, FA; apparent fiber density, AFD) and sensorimotor performance (hand skill and postural stability). Compared to controls, our results revealed increased MD in both upper and lower limb thalamocortical tracts in the non-dominant hemisphere in HP and bilaterally in DP subgroup. MD was strongly lateralized in participants with hemiplegia, while AFD seemed lateralized only in controls. fMRI-based tractography results were comparable. The correlation analysis indicated an association between the white matter structure and sensorimotor performance. These findings suggest distinct impairment of functionally relevant thalamocortical pathways in HP and DP subtypes. Thus, the organization of thalamocortical white matter tracts may offer valuable guidance for targeted, life-long rehabilitation in children with CP.

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.

brainlife.io: A decentralized and open source cloud platform to support neuroscience research

Neuroscience research has expanded dramatically over the past 30 years by advancing standardization and tool development to support rigor and transparency. Consequently, the complexity of the data pipeline has also increased, hindering access to FAIR data analysis to portions of the worldwide research community. brainlife.io was developed to reduce these burdens and democratize modern neuroscience research across institutions and career levels. Using community software and hardware infrastructure, the platform provides open-source data standardization, management, visualization, and processing and simplifies the data pipeline. brainlife.io automatically tracks the provenance history of thousands of data objects, supporting simplicity, efficiency, and transparency in neuroscience research. Here brainlife.io's technology and data services are described and evaluated for validity, reliability, reproducibility, replicability, and scientific utility. Using data from 4 modalities and 3,200 participants, we demonstrate that brainlife.io's services produce outputs that adhere to best practices in modern neuroscience research.

Structural and functional connectivity reconstruction with CATO - A Connectivity Analysis TOolbox

We describe a Connectivity Analysis TOolbox (CATO) for the reconstruction of structural and functional brain connectivity based on diffusion weighted imaging and resting-state functional MRI data. CATO is a multimodal software package that enables researchers to run end-to-end reconstructions from MRI data to structural and functional connectome maps, customize their analyses and utilize various software packages to preprocess data. Structural and functional connectome maps can be reconstructed with respect to user-defined (sub)cortical atlases providing aligned connectivity matrices for integrative multimodal analyses. We outline the implementation and usage of the structural and functional processing pipelines in CATO. Performance was calibrated with respect to simulated diffusion weighted imaging from the ITC2015 challenge, test-retest diffusion weighted imaging data and resting-state functional MRI data from the Human Connectome Project. CATO is open-source software distributed under the MIT License and available as a MATLAB toolbox and as a stand-alone application at www.dutchconnectomelab.nl/CATO.

A taxonomy of the brain's white matter: twenty-one major tracts for the 21st century

The functional and computational properties of brain areas are determined, in large part, by their connectivity profiles. Advances in neuroimaging and network neuroscience allow us to characterize the human brain noninvasively, but a comprehensive understanding of the human brain demands an account of the anatomy of brain connections. Long-range anatomical connections are instantiated by white matter, which itself is organized into tracts. These tracts are often disrupted by central nervous system disorders, and they can be targeted by neuromodulatory interventions, such as deep brain stimulation. Here, we characterized the connections, morphology, traversal, and functions of the major white matter tracts in the brain. There are major discrepancies across different accounts of white matter tract anatomy, hindering our attempts to accurately map the connectivity of the human brain. However, we are often able to clarify the source(s) of these discrepancies through careful consideration of both histological tract-tracing and diffusion-weighted tractography studies. In combination, the advantages and disadvantages of each method permit novel insights into brain connectivity. Ultimately, our synthesis provides an essential reference for neuroscientists and clinicians interested in brain connectivity and anatomy, allowing for the study of the association of white matter's properties with behavior, development, and disorders.

DORIS: A diffusion MRI-based 10 tissue class deep learning segmentation algorithm tailored to improve anatomically-constrained tractography

Modern tractography algorithms such as anatomically-constrained tractography (ACT) are based on segmentation maps of white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF). These maps are generally estimated from a T1-weighted (T1w) image and then registered in diffusion weighted images (DWI) space. Registration of T1w to diffusion space and partial volume estimation are challenging and rarely voxel-perfect. Diffusion-based segmentation would, thus, potentially allow not to have higher quality anatomical priors injected in the tractography process. On the other hand, even if FA-based tractography is possible without T1 registration, the literature shows that this technique suffers from multiple issues such as holes in the tracking mask and a high proportion of generated broken and anatomically implausible streamlines. Therefore, there is an important need for a tissue segmentation algorithm that works directly in the native diffusion space. We propose DORIS, a DWI-based deep learning segmentation algorithm. DORIS outputs 10 different tissue classes including WM, GM, CSF, ventricles, and 6 other subcortical structures (putamen, pallidum, hippocampus, caudate, amygdala, and thalamus). DORIS was trained and validated on a wide range of subjects, including 1,000 individuals from 22 to 90 years old from clinical and research DWI acquisitions, from 5 public databases. In the absence of a "true" ground truth in diffusion space, DORIS used a silver standard strategy from Freesurfer output registered onto the DWI. This strategy is extensively evaluated and discussed in the current study. Segmentation maps provided by DORIS are quantitatively compared to Freesurfer and FSL-fast and the impacts on tractography are evaluated. Overall, we show that DORIS is fast, accurate, and reproducible and that DORIS-based tractograms produce bundles with a longer mean length and fewer anatomically implausible streamlines.

Bundle-o-graphy: improving structural connectivity estimation with adaptive microstructure-informed tractography

Tractography is a powerful tool for the investigation of the complex organization of the brain in vivo, as it allows inferring the macroscopic pathways of the major fiber bundles of the white matter based on non-invasive diffusion-weighted magnetic resonance imaging acquisitions. Despite this unique and compelling ability, some studies have exposed the poor anatomical accuracy of the reconstructions obtained with this technique and challenged its effectiveness for studying brain connectivity. In this work, we describe a novel method to readdress tractography reconstruction problem in a global manner by combining the strengths of so-called generative and discriminative strategies. Starting from an input tractogram, we parameterize the connections between brain regions following a bundle-based representation that allows to drastically reducing the number of parameters needed to model groups of fascicles. The parameters space is explored following an MCMC generative approach, while a discrimininative method is exploited to globally evaluate the set of connections which is updated according to Bayes' rule. Our results on both synthetic and real brain data show that the proposed solution, called bundle-o-graphy, allows improving the anatomical accuracy of the reconstructions while keeping the computational complexity similar to other state-of-the-art methods.

The role of neuronavigation in TMS-EEG studies: current applications and future perspectives

Transcranial magnetic stimulation combined with electroencephalography (TMS-EEG) allows measuring non-invasively the electrical response of the human cerebral cortex to a direct perturbation. Complementing TMS-EEG with a structural neuronavigation tool (nTMS-EEG) is key for accurately selecting cortical areas, targeting them, and adjusting the stimulation parameters based on some relevant anatomical priors. This step, together with the employment of visualization tools designed to perform a quality check of TMS-evoked potentials (TEPs) in real-time during acquisition, is key for maximizing the impact of the TMS pulse on the cortex and in ensuring highly reproducible measurements within sessions and across subjects. Moreover, storing stimulation parameters in the neuronavigation system can help in reproducing the stimulation parameters within and across experimental sessions and sharing them across research centers. Finally, the systematic employment of neuronavigation in TMS-EEG studies is also key to standardize measurements in clinical populations in search for reliable diagnostic and prognostic TMS-EEG-based biomarkers for neurological and psychiatric disorders.

Retinal ganglion cell endowment is correlated with optic tract fiber cross section, not density

There is substantial variation between healthy individuals in the number of retinal ganglion cells (RGC) in the eye, with commensurate variation in the number of axons in the optic tracts. Fixel-based analysis of diffusion MR produces estimates of fiber density (FD) and cross section (FC). Using these fixel measurements along with retinal imaging, we asked if individual differences in RGC tissue volume are correlated with individual differences in FD and FC measurements obtained from the optic tracts, and subsequent structures along the cortical visual pathway. We find that RGC endowment is correlated with optic tract FC, but not with FD. RGC volume had a decreasing relationship with measurements from subsequent regions of the visual system (LGN volume, optic radiation FC/FD, and V1 surface area). However, we also found that the variations in each visual area were correlated with the variations in its immediately adjacent visual structure. We only observed these serial correlations when FC is used as the measure of interest for the optic tract and radiations, but no significant relationship was found when FD represented these white matter structures. From these results, we conclude that the variations in RGC endowment, LGN volume, and V1 surface area are better predicted by the overall cross section of the optic tract and optic radiations as compared to the intra-axonal restricted signal component of these white matter pathways. Additionally, the presence of significant correlations between adjacent, but not distant, anatomical structures suggests that there are multiple, local sources of anatomical variation along the visual pathway.

Unsupervised Deep Learning for FOD-Based Susceptibility Distortion Correction in Diffusion MRI

Susceptibility induced distortion is a major artifact that affects the diffusion MRI (dMRI) data analysis. In the Human Connectome Project (HCP), the state-of-the-art method adopted to correct this kind of distortion is to exploit the displacement field from the B0 image in the reversed phase encoding images. However, both the traditional and learning-based approaches have limitations in achieving high correction accuracy in certain brain regions, such as brainstem. By utilizing the fiber orientation distribution (FOD) computed from the dMRI, we propose a novel deep learning framework named DistoRtion Correction Net (DrC-Net), which consists of the U-Net to capture the latent information from the 4D FOD images and the spatial transformer network to propagate the displacement field and back propagate the losses between the deformed FOD images. The experiments are performed on two datasets acquired with different phase encoding (PE) directions including the HCP and the Human Connectome Low Vision (HCLV) dataset. Compared to two traditional methods topup and FODReg and two deep learning methods S-Net and flow-net, the proposed method achieves significant improvements in terms of the mean squared difference (MSD) of fractional anisotropy (FA) images and minimum angular difference between two PEs in white matter and also brainstem regions. In the meantime, the proposed DrC-Net takes only several seconds to predict a displacement field, which is much faster than the FODReg method.

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

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

Identification of Proprioceptive Thalamocortical Tracts in Children: Comparison of fMRI, MEG, and Manual Seeding of Probabilistic Tractography

Studying white matter connections with tractography is a promising approach to understand the development of different brain processes, such as proprioception. An emerging method is to use functional brain imaging to select the cortical seed points for tractography, which is considered to improve the functional relevance and validity of the studied connections. However, it is unknown whether different functional seeding methods affect the spatial and microstructural properties of the given white matter connection. Here, we compared functional magnetic resonance imaging, magnetoencephalography, and manual seeding of thalamocortical proprioceptive tracts for finger and ankle joints separately. We showed that all three seeding approaches resulted in robust thalamocortical tracts, even though there were significant differences in localization of the respective proprioceptive seed areas in the sensorimotor cortex, and in the microstructural properties of the obtained tracts. Our study shows that the selected functional or manual seeding approach might cause systematic biases to the studied thalamocortical tracts. This result may indicate that the obtained tracts represent different portions and features of the somatosensory system. Our findings highlight the challenges of studying proprioception in the developing brain and illustrate the need for using multimodal imaging to obtain a comprehensive view of the studied brain process.