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

Sex-related variability of white matter tracts in the whole HCP cohort

Behavioral differences between men and women have been studied extensively, as have differences in brain anatomy. However, most studies have focused on differences in gray matter, while white matter has been much less studied. We conducted a comprehensive study of 77 deep white matter tracts to analyze their volumetric and microstructural variability between men and women in the full Human Connectome Project (HCP) cohort of 1065 healthy individuals aged 22-35 years. We found a significant difference in total brain volume between men and women (+ 12.6% in men), consistent with the literature. 16 tracts showed significant volumetric differences between men and women, one of which stood out due to a larger effect size: the corpus callosum genu, which was larger in women (+ 7.3% in women, p = 5.76 × 10-19). In addition, we found several differences in microstructural parameters between men and women, both using standard Diffusion Tensor Imaging (DTI) parameters and more complex microstructural parameters from the Neurite Orientation Dispersion and Density Imaging (NODDI) model, with the tracts showing the greatest differences belonging to motor (cortico-spinal tracts, cortico-cerebellar tracts) or limbic (cingulum, fornix, thalamo-temporal radiations) systems. These microstructural differences may be related to known behavioral differences between the sexes in timed motor performance, aggressiveness/impulsivity, and social cognition.

Microstructural Properties of the Cerebellar Peduncles in Children With Developmental Language Disorder

Children with developmental language disorder (DLD) struggle to learn their native language for no apparent reason. While research on the neurobiological underpinnings of the disorder has focused on the role of corticostriatal systems, little is known about the role of the cerebellum in DLD. Corticocerebellar circuits might be involved in the disorder as they contribute to complex sensorimotor skill learning, including the acquisition of spoken language. Here, we used diffusion-weighted imaging data from 77 typically developing and 54 children with DLD and performed probabilistic tractography to identify the cerebellum's white matter tracts: the inferior, middle, and superior cerebellar peduncles. Children with DLD showed lower fractional anisotropy (FA) in the inferior cerebellar peduncles (ICP), fiber tracts that carry motor and sensory input via the inferior olive to the cerebellum. Lower FA in DLD was driven by lower axial diffusivity. Probing this further with more sophisticated modeling of diffusion data, we found higher orientation dispersion but no difference in neurite density in the ICP of children with DLD. Reduced FA is therefore unlikely to be reflecting microstructural differences in myelination, rather the organization of axons in these pathways is disrupted. ICP microstructure was not associated with language or motor coordination performance in our sample. We also found no differences in the middle and superior peduncles, the main pathways connecting the cerebellum with the cortex. To conclude, it is not corticocerebellar but atypical olivocerebellar white matter connections that characterize DLD and suggest the involvement of the olivocerebellar system in speech and language acquisition and development.

Methodological considerations on diffusion MRI tractography in infants aged 0-2 years: a scoping review

Diffusion MRI (dMRI) enables studying the complex architectural organization of the brain's white matter (WM) through virtual reconstruction of WM fiber tracts (tractography). Despite the anticipated clinical importance of applying tractography to study structural connectivity and tract development during the critical period of rapid infant brain maturation, detailed descriptions on how to approach tractography in young infants are limited. Over the past two decades, tractography from infant dMRI has mainly been applied in research settings and focused on diffusion tensor imaging (DTI). Only few studies used techniques superior to DTI in terms of disentangling information on the brain's organizational complexity, including crossing fibers. While more advanced techniques may enhance our understanding of the intricate processes of normal and abnormal brain development and extensive knowledge has been gained from application on adult scans, their applicability in infants has remained underexplored. This may partially be due to the higher technical requirements versus the need to limit scan time in young infants. We review various previously described methodological practices for tractography in the infant brain (0-2 years-of-age) and provide recommendations to optimize advanced tractography approaches to enable more accurate reconstructions of the brain WM's complexity. IMPACT: Diffusion tensor imaging is the technique most frequently used for fiber tracking in the developing infant brain but is limited in capability to disentangle the complex white matter organization. Advanced tractography techniques allow for reconstruction of crossing fiber bundles to better reflect the brain's complex organization. Yet, they pose practical and technical challenges in the fast developing young infant's brain. Methods on how to approach advanced tractography in the young infant's brain have hardly been described. Based on a literature review, recommendations are provided to optimize tractography for the developing infant brain, aiming to advance early diagnosis and neuroprotective strategies.

Reproducibility of arterial spin labeling cerebral blood flow image processing: A report of the ISMRM open science initiative for perfusion imaging (OSIPI) and the ASL MRI challenge

PURPOSE

Arterial spin labeling (ASL) is a widely used contrast-free MRI method for assessing cerebral blood flow (CBF). Despite the generally adopted ASL acquisition guidelines, there is still wide variability in ASL analysis. We explored this variability through the ISMRM-OSIPI ASL-MRI Challenge, aiming to establish best practices for more reproducible ASL analysis.

METHODS

Eight teams analyzed the challenge data, which included a high-resolution T1-weighted anatomical image and 10 pseudo-continuous ASL datasets simulated using a digital reference object to generate ground-truth CBF values in normal and pathological states. We compared the accuracy of CBF quantification from each team's analysis to the ground truth across all voxels and within predefined brain regions. Reproducibility of CBF across analysis pipelines was assessed using the intra-class correlation coefficient (ICC), limits of agreement (LOA), and replicability of generating similar CBF estimates from different processing approaches.

RESULTS

Absolute errors in CBF estimates compared to ground-truth synthetic data ranged from 18.36 to 48.12 mL/100 g/min. Realistic motion incorporated into three datasets produced the largest absolute error and variability between teams, with the least agreement (ICC and LOA) with ground-truth results. Fifty percent of the submissions were replicated, and one produced three times larger CBF errors (46.59 mL/100 g/min) compared to submitted results.

CONCLUSIONS

Variability in CBF measurements, influenced by differences in image processing, especially to compensate for motion, highlights the significance of standardizing ASL analysis workflows. We provide a recommendation for ASL processing based on top-performing approaches as a step toward ASL standardization.

Imaging of developing human brains with ex vivo PSOCT and dMRI

The human brain undergoes substantial developmental changes in the first five years of life. Particularly in the white matter, myelination of axons occurs near birth and continues at a rapid pace during the first 2 to 3 years. Diffusion MRI (dMRI) has revolutionized our understanding of developmental trajectories in white matter. However, the mm-resolution of in vivo techniques bears significant limitation in revealing the microstructure of the developing brain. Polarization sensitive optical coherence tomography (PSOCT) is a three-dimensional (3D) optical imaging technique that uses polarized light interferometry to target myelinated fiber tracts with micrometer resolution. Previous studies have shown that PSOCT contributes significantly to the elucidation of myelin content and quantification of fiber orientation in adult human brains. In this study, we utilized the PSOCT technique to study developing brains during the first 5 years of life in combination with ex vivo dMRI. The results showed that the optical properties of PSOCT quantitatively reveal the myelination process in young children. The imaging contrast of the optic axis orientation is a sensitive measure of fiber orientations in largely unmyelinated brains as young as 3-months-old. The micrometer resolution of PSOCT provides substantially enriched information about complex fiber networks and complements submillimeter dMRI. This new optical tool offers great potential to reveal the white matter structures in normal neurodevelopment and developmental disorders in unprecedented detail.

White matter tract crossing and bottleneck regions in the fetal brain

There is a growing interest in using diffusion MRI to study the white matter tracts and structural connectivity of the fetal brain. Recent progress in data acquisition and processing suggests that this imaging modality has a unique role in elucidating the normal and abnormal patterns of neurodevelopment in utero. However, there have been no efforts to quantify the prevalence of crossing tracts and bottleneck regions, important issues that have been extensively researched for adult brains. In this work, we determined the brain regions with crossing tracts and bottlenecks between 23 and 36 gestational weeks. We performed probabilistic tractography on 59 fetal brain scans and extracted a set of 51 distinct white tracts, which we grouped into 10 major tract bundle groups. We analyzed the results to determine the patterns of tract crossings and bottlenecks. Our results showed that 20-25% of the white matter voxels included two or three crossing tracts. Bottlenecks were more prevalent. Between 75-80% of the voxels were characterized as bottlenecks, with more than 40% of the voxels involving four or more tracts. The results of this study highlight the challenge of fetal brain tractography and structural connectivity assessment and call for innovative image acquisition and analysis methods to mitigate these problems.

Cholinergic modulation supports dynamic switching of resting state networks through selective DMN suppression

Brain activity during the resting state is widely used to examine brain organization, cognition and alterations in disease states. While it is known that neuromodulation and the state of alertness impact resting-state activity, neural mechanisms behind such modulation of resting-state activity are unknown. In this work, we used a computational model to demonstrate that change in excitability and recurrent connections, due to cholinergic modulation, impacts resting-state activity. The results of such modulation in the model match closely with experimental work on direct cholinergic modulation of Default Mode Network (DMN) in rodents. We further extended our study to the human connectome derived from diffusion-weighted MRI. In human resting-state simulations, an increase in cholinergic input resulted in a brain-wide reduction of functional connectivity. Furthermore, selective cholinergic modulation of DMN closely captured experimentally observed transitions between the baseline resting state and states with suppressed DMN fluctuations associated with attention to external tasks. Our study thus provides insight into potential neural mechanisms for the effects of cholinergic neuromodulation on resting-state activity and its dynamics.

Individualized brain mapping for navigated neuromodulation

ABSTRACT

The brain is a complex organ that requires precise mapping to understand its structure and function. Brain atlases provide a powerful tool for studying brain circuits, discovering biological markers for early diagnosis, and developing personalized treatments for neuropsychiatric disorders. Neuromodulation techniques, such as transcranial magnetic stimulation and deep brain stimulation, have revolutionized clinical therapies for neuropsychiatric disorders. However, the lack of fine-scale brain atlases limits the precision and effectiveness of these techniques. Advances in neuroimaging and machine learning techniques have led to the emergence of stereotactic-assisted neurosurgery and navigation systems. Still, the individual variability among patients and the diversity of brain diseases make it necessary to develop personalized solutions. The article provides an overview of recent advances in individualized brain mapping and navigated neuromodulation and discusses the methodological profiles, advantages, disadvantages, and future trends of these techniques. The article concludes by posing open questions about the future development of individualized brain mapping and navigated neuromodulation.

Hemispheric lateralization of language processing: insights from network-based symptom mapping and patient subgroups

The hemispheric laterality of language processing has become a hot topic in modern neuroscience. Although most previous studies have reported left-lateralized language processing, other studies found it to be bilateral. A previous neurocomputational model has proposed a unified framework to explain that the above discrepancy might be from healthy and patient individuals. This model posits an initial symmetry but imbalanced capacity in language processing for healthy individuals, with this imbalance contributing to language recovery disparities following different hemispheric injuries. The present study investigated this model by analyzing the lateralization patterns of language subnetworks across multiple attributes with a group of 99 patients (compared to nonlanguage processing) and examining the lateralization patterns of language subnetworks in subgroups with damage to different hemispheres. Subnetworks were identified using a whole-brain network-based lesion-symptom mapping method, and the lateralization index was quantitatively measured. We found that all the subnetworks in language processing were left-lateralized, while subnetworks in nonlanguage processing had different lateralization patterns. Moreover, diverse hemisphere-injury subgroups exhibited distinct language recovery effects. These findings provide robust support for the proposed neurocomputational model of language processing.

The human brain connectome weighted by the myelin content and total intra-axonal cross-sectional area of white matter tracts

A central goal in neuroscience is the development of a comprehensive mapping between structural and functional brain features, which facilitates mechanistic interpretation of brain function. However, the interpretability of structure-function brain models remains limited by a lack of biological detail. Here, we characterize human structural brain networks weighted by multiple white matter microstructural features including total intra-axonal cross-sectional area and myelin content. We report edge-weight-dependent spatial distributions, variance, small-worldness, rich club, hubs, as well as relationships with function, edge length, and myelin. Contrasting networks weighted by the total intra-axonal cross-sectional area and myelin content of white matter tracts, we find opposite relationships with functional connectivity, an edge-length-independent inverse relationship with each other, and the lack of a canonical rich club in myelin-weighted networks. When controlling for edge length, networks weighted by either fractional anisotropy, radial diffusivity, or neurite density show no relationship with whole-brain functional connectivity. We conclude that the co-utilization of structural networks weighted by total intra-axonal cross-sectional area and myelin content could improve our understanding of the mechanisms mediating the structure-function brain relationship.

Violations of the fluctuation-dissipation theorem reveal distinct nonequilibrium dynamics of brain states

The brain is a nonequilibrium system whose dynamics change in different brain states, such as wakefulness and deep sleep. Thermodynamics provides the tools for revealing these nonequilibrium dynamics. We used violations of the fluctuation-dissipation theorem to describe the hierarchy of nonequilibrium dynamics associated with different brain states. Together with a whole-brain model fitted to empirical human neuroimaging data, and deriving the appropriate analytical expressions, we were able to capture the deviation from equilibrium in different brain states that arises from asymmetric interactions and hierarchical organization.

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.

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.

Functional MRI but not white matter fibre dissection identifies language dominance

OBJECTIVES

Lateralisation of some language pathways has been reported in the literature using diffusion tractography, which is more feasible than functional magnetic resonance imaging (fMRI) in challenging patients. Our retrospective study investigates whether a correlation exists between threshold-independent fMRI language lateralisation and structural lateralisation using tractography in healthy controls and brain tumour patients.

METHODS

Fifteen healthy subjects and 61 patients underwent language fMRI and diffusion-weighted MRI. A regional fMRI laterality index (LI) was calculated. Tracts dissected were the arcuate fasciculus (long direct and short indirect tracts), uncinate fasciculus, inferior longitudinal fasciculus, inferior fronto-occipital fasciculus and frontal aslant tract. An asymmetry index (AI) for each tract was calculated using tract volume analysed with single tensor (ST) and spherical deconvolution (SD) models, as well as hindrance modulated orientational anisotropy (HMOA) for SD tracts. Linear regression assessed the correlation between LI and AI.

RESULTS

In all subjects, there was no significant correlation between LI and AI for any of the dissected tracts. Significant correlations were only found when handedness for controls and tumour volume for patients were included as covariates. In handedness subgroups, the average AI of some tracts showed the same laterality as LI, and some the opposite. Discordant results were observed for ST- and SD-based AIs.

CONCLUSIONS

Our results do not support using tractography in the assessment of language lateralisation. The discordant results between ST and SD indicate that either the structural lateralisation of dissected tracts is less robust than functional lateralisation, or tractography is not sensitive methodology. Other diffusion analysis approaches should be developed.

CLINICAL RELEVANCE STATEMENT

Although diffusion tractography may be more feasible than fMRI in challenging tumour patients and where sedation or anaesthesia is required, our results do not currently recommend replacing fMRI with tractography using volume or HMOA in the assessment of language lateralisation.

KEY POINTS

• No correlation found between fMRI and tractography in language lateralisation. • Discordance between asymmetry indices of different tractography models and metrics. • Tractography not currently recommended in language lateralisation assessment.

Functional Connectivity and Microstructural Network Correlates of Interindividual Variability in Distinct Executive Functions of Healthy Older Adults

Executive functions, essential for daily life, are known to be impaired in older age. Some executive functions, including working memory updating and value-based decision-making, are specifically sensitive to age-related deterioration. While their neural correlates in young adults are well-described, a comprehensive delineation of the underlying brain substrates in older populations, relevant to identify targets for modulation against cognitive decline, is missing. Here, we assessed letter updating and Markov decision-making task performance to operationalize these trainable functions in 48 older adults. Resting-state functional magnetic resonance imaging was acquired to quantify functional connectivity (FC) in task-relevant frontoparietal and default mode networks. Microstructure in white matter pathways mediating executive functions was assessed with diffusion tensor imaging and quantified by tract-based fractional anisotropy (FA). Superior letter updating performance correlated with higher FC between dorsolateral prefrontal cortex and left frontoparietal and hippocampal areas, while superior Markov decision-making performance correlated with decreased FC between basal ganglia and right angular gyrus. Furthermore, better working memory updating performance was related to higher FA in the cingulum bundle and the superior longitudinal fasciculus. Stepwise linear regression showed that cingulum bundle FA added significant incremental contribution to the variance explained by fronto-angular FC alone. Our findings provide a characterization of distinct functional and structural connectivity correlates associated with performance of specific executive functions. Thereby, this study contributes to the understanding of the neural correlates of updating and decision-making functions in older adults, paving the way for targeted modulation of specific networks by modulatory techniques such as behavioral interventions and non-invasive brain stimulation.

Microstructural Properties of the Cerebellar Peduncles in Children with Developmental Language Disorder

Children with developmental language disorder (DLD) struggle to learn their native language for no apparent reason. While research on the neurobiological underpinnings of the disorder has focused on the role of cortico-striatal systems, little is known about the role of the cerebellum in DLD. Cortico-cerebellar circuits might be involved in the disorder as they contribute to complex sensorimotor skill learning, including the acquisition of spoken language. Here, we used diffusion-weighted imaging data from 77 typically developing and 54 children with DLD and performed probabilistic tractography to identify the cerebellum's white matter tracts: the inferior, middle, and superior cerebellar peduncles. Children with DLD showed lower fractional anisotropy (FA) in the inferior cerebellar peduncles (ICP), fiber tracts that carry motor and sensory input via the inferior olive to the cerebellum. Lower FA in DLD was driven by lower axial diffusivity. Probing this further with more sophisticated modeling of diffusion data, we found higher orientation dispersion but no difference in neurite density in the ICP of DLD. Reduced FA is therefore unlikely to be reflecting microstructural differences in myelination in this tract, rather the organization of axons in these pathways is disrupted. ICP microstructure was not associated with language or motor coordination performance in our sample. We also found no differences in the middle and superior peduncles, the main pathways connecting the cerebellum with the cortex. To conclude, it is not cortico-cerebellar but atypical olivocerebellar white matter connections that characterize DLD and suggest the involvement of the olivocerebellar system in speech acquisition and development.

Direct segmentation of brain white matter tracts in diffusion MRI

The brain white matter consists of a set of tracts that connect distinct regions of the brain. Segmentation of these tracts is often needed for clinical and research studies. Diffusion-weighted MRI offers unique contrast to delineate these tracts. However, existing segmentation methods rely on intermediate computations such as tractography or estimation of fiber orientation density. These intermediate computations, in turn, entail complex computations that can result in unnecessary errors. Moreover, these intermediate computations often require dense multi-shell measurements that are unavailable in many clinical and research applications. As a result, current methods suffer from low accuracy and poor generalizability. Here, we propose a new deep learning method that segments these tracts directly from the diffusion MRI data, thereby sidestepping the intermediate computation errors. Our experiments show that this method can achieve segmentation accuracy that is on par with the state of the art methods (mean Dice Similarity Coefficient of 0.826). Compared with the state of the art, our method offers far superior generalizability to undersampled data that are typical of clinical studies and to data obtained with different acquisition protocols. Moreover, we propose a new method for detecting inaccurate segmentations and show that it is more accurate than standard methods that are based on estimation uncertainty quantification. The new methods can serve many critically important clinical and scientific applications that require accurate and reliable non-invasive segmentation of white matter tracts.

One ring to rule them all: The unifying role of prefrontal cortex in steering task-related brain dynamics

Surviving and thriving in a complex world require intricate balancing of higher order brain functions with essential survival-related behaviours. Exactly how this is achieved is not fully understood but a large body of work has shown that different regions in the prefrontal cortex (PFC) play key roles for diverse cognitive and emotional tasks including emotion, control, response inhibition, mental set shifting and working memory. We hypothesised that the key regions are hierarchically organised and we developed a framework for discovering the driving brain regions at the top of the hierarchy, responsible for steering the brain dynamics of higher brain function. We fitted a time-dependent whole-brain model to the neuroimaging data from large-scale Human Connectome Project with over 1,000 participants and computed the entropy production for rest and seven tasks (covering the main domains of cognition). This thermodynamics framework allowed us to identify the main common, unifying drivers steering the orchestration of brain dynamics during difficult tasks; located in key regions of the PFC (inferior frontal gyrus, lateral orbitofrontal cortex, rostral and caudal frontal cortex and rostral anterior cingulate cortex). Selectively lesioning these regions in the whole-brain model demonstrated their causal mechanistic importance. Overall, this shows the existence of a 'ring' of specific PFC regions ruling over the orchestration of higher brain function.

Fiber-specific white matter alterations in Parkinson's disease patients with freezing of gait

Freezing of gait (FOG) is a gait disorder that usually occurs in advanced stages of Parkinson's disease (PD). Understanding the underlying mechanism of FOG is important for treatment and prevention. Previous studies have investigated white matter (WM) structure to explore the pathology of FOG. However, the pathology is still unclear, possibly due to the methodological limitation in identifying specific fiber tracts. This study aimed to investigate tract-specific WM structural changes in FOG patients and their relationships with clinical characteristics. We enrolled 19 PD patients with FOG (PD-FOG), 19 without FOG (PD-woFOG) and 21 controls. Fixel-based analysis is a novel framework to avoid the effect of crossing fibers, which provides the metrics to assess WM morphology. By combining a method for segmenting fibers, we identified abnormalities in the specific fiber tracts. Compared to PD-woFOG, PD-FOG showed significant increased fiber-bundle cross-section (FC) in the corpus callosum (CC), fornix (FX), inferior longitudinal fasciculus (ILF), striato-premotor (ST_PREM), superior thalamic radiation (STR), thalamo-premotor (T_PREM), increased fiber density and cross-section (FDC) in the STR, and decreased fiber density (FD) in the CC and ILF. Additionally, the ILF was correlated with motor, cognition and memory, the CC was correlated with anxiety, and the T_PREM was also correlated with cognition. In conclusion, in addition to impairments of WM found in PD-FOG, we found enhancements in WM, which may imply compensatory mechanisms. Furthermore, multiple fiber tracts were correlated with clinical characteristics, especially the ILF, validating the involvement of transmission circuits of multiple distinct information in mechanisms of FOG.

ModelArray: An R package for statistical analysis of fixel-wise data

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Generative Sampling in Bundle Tractography using Autoencoders (GESTA)

Current tractography methods use the local orientation information to propagate streamlines from seed locations. Many such seeds provide streamlines that stop prematurely or fail to map the true white matter pathways because some bundles are "harder-to-track" than others. This results in tractography reconstructions with poor white and gray matter spatial coverage. In this work, we propose a generative, autoencoder-based method, named GESTA (Generative Sampling in Bundle Tractography using Autoencoders), that produces streamlines achieving better spatial coverage. Compared to other deep learning methods, our autoencoder-based framework uses a single model to generate streamlines in a bundle-wise fashion, and does not require to propagate local orientations. GESTA produces new and complete streamlines for any given white matter bundle, including hard-to-track bundles. Applied on top of a given tractogram, GESTA is shown to be effective in improving the white matter volume coverage in poorly populated bundles, both on synthetic and human brain in vivo data. Our streamline evaluation framework ensures that the streamlines produced by GESTA are anatomically plausible and fit well to the local diffusion signal. The streamline evaluation criteria assess anatomy (white matter coverage), local orientation alignment (direction), and geometry features of streamlines, and optionally, gray matter connectivity. The GESTA framework offers considerable gains in bundle overlap using a reduced set of seeding streamlines with a 1.5x improvement for the "Fiber Cup", and 6x for the ISMRM 2015 Tractography Challenge datasets. Similarly, it provides a 4x white matter volume increase on the BIL&GIN callosal homotopic dataset, and successfully populates bundles on the multi-subject, multi-site, whole-brain in vivo TractoInferno dataset. GESTA is thus a novel deep generative bundle tractography method that can be used to improve the tractography reconstruction of the white matter.

Spatially regularized low-rank tensor approximation for accurate and fast tractography

Tractography based on diffusion Magnetic Resonance Imaging (dMRI) is the prevalent approach to the in vivo delineation of white matter tracts in the human brain. Many tractography methods rely on models of multiple fiber compartments, but the local dMRI information is not always sufficient to reliably estimate the directions of secondary fibers. Therefore, we introduce two novel approaches that use spatial regularization to make multi-fiber tractography more stable. Both represent the fiber Orientation Distribution Function (fODF) as a symmetric fourth-order tensor, and recover multiple fiber orientations via low-rank approximation. Our first approach computes a joint approximation over suitably weighted local neighborhoods with an efficient alternating optimization. The second approach integrates the low-rank approximation into a current state-of-the-art tractography algorithm based on the unscented Kalman filter (UKF). These methods were applied in three different scenarios. First, we demonstrate that they improve tractography even in high-quality data from the Human Connectome Project, and that they maintain useful results with a small fraction of the measurements. Second, on the 2015 ISMRM tractography challenge, they increase overlap, while reducing overreach, compared to low-rank approximation without joint optimization or the traditional UKF, respectively. Finally, our methods permit a more comprehensive reconstruction of tracts surrounding a tumor in a clinical dataset. Overall, both approaches improve reconstruction quality. At the same time, our modified UKF significantly reduces the computational effort compared to its traditional counterpart, and to our joint approximation. However, when used with ROI-based seeding, joint approximation more fully recovers fiber spread.

Tractography-Enhanced Biopsy of Central Core Motor Eloquent Tumours: A Simulation-Based Study

Safe Trajectory planning for navigation guided biopsy (nBx) of motor eloquent tumours (METs) is important to minimise neurological morbidity. Preliminary clinical data suggest that visualisation of the corticospinal tract (CST) and its relation to the tumour may aid in planning a safe trajectory. In this article we assess the impact of tractography in nBx planning in a simulation-based exercise. This single centre cross-sectional study was performed in March 2021 including 10 patients with METs divided into 2 groups: (1) tractography enhanced group (T-nBx; n = 5; CST merged with volumetric MRI); (2) anatomy-based group (A-nBx; n = 5; volumetric MRI only). A biopsy target was chosen on each tumour. Volunteer neurosurgical trainees had to plan a suitable biopsy trajectory on a Stealth S8® workstation for all patients in a single session. A trajectory safety index (TSI) was devised for each trajectory. Data collection and analysis included a comparison of trajectory planning time, trajectory/lobe changes and TSI. A total of 190 trajectories were analysed based on participation from 19 trainees. Mean trajectory planning time for the entire cohort was 225.1 ± 21.97 s. T-nBx required shorter time for planning (p = 0.01). Mean trajectory changes and lobe changes made per biopsy were 3.28 ± 0.29 and 0.45 ± 0.08, respectively. T-nBx required fewer trajectory/lobe changes (p = 0.01). TSI was better in the presence of tractography than A-nBx (p = 0.04). Neurosurgical experience of trainees had no significant impact on the measured parameters despite adjusted analysis. Irrespective of the level of neurosurgical training, surgical planning of navigation guided biopsy for METs may be achieved in less time with a safer trajectory if tractography imaging is available.

Coordinated anatomical and functional variability in the human brain during adolescence

Adolescence represents a time of unparalleled brain development. In particular, developmental changes in morphometric and cytoarchitectural features are accompanied by maturation in the functional connectivity (FC). Here, we examined how three facets of the brain, including myelination, cortical thickness (CT), and resting-state FC, interact in children between the ages of 10 and 15. We investigated the pattern of coordination in these measures by computing correlation matrices for each measure as well as meta-correlations among them both at the regional and network levels. The results revealed consistently higher meta-correlations among myelin, CT, and FC in the sensory-motor cortical areas than in the association cortical areas. We also found that these meta-correlations were stable and little affected by age-related changes in each measure. In addition, regional variations in the meta-correlations were consistent with the previously identified gradient in the FC and therefore reflected the hierarchy of cortical information processing, and this relationship persists in the adult brain. These results demonstrate that heterogeneity in FC among multiple cortical areas are closely coordinated with the development of cortical myelination and thickness during adolescence.

Prediction of the Topography of the Corticospinal Tract on T1-Weighted MR Images Using Deep-Learning-Based Segmentation

INTRODUCTION

Tractography is an invaluable tool in the planning of tumor surgery in the vicinity of functionally eloquent areas of the brain as well as in the research of normal development or of various diseases. The aim of our study was to compare the performance of a deep-learning-based image segmentation for the prediction of the topography of white matter tracts on T1-weighted MR images to the performance of a manual segmentation.

METHODS

T1-weighted MR images of 190 healthy subjects from 6 different datasets were utilized in this study. Using deterministic diffusion tensor imaging, we first reconstructed the corticospinal tract on both sides. After training a segmentation model on 90 subjects of the PIOP2 dataset using the nnU-Net in a cloud-based environment with graphical processing unit (Google Colab), we evaluated its performance using 100 subjects from 6 different datasets.

RESULTS

Our algorithm created a segmentation model that predicted the topography of the corticospinal pathway on T1-weighted images in healthy subjects. The average dice score was 0.5479 (0.3513-0.7184) on the validation dataset.

CONCLUSIONS

Deep-learning-based segmentation could be applicable in the future to predict the location of white matter pathways in T1-weighted scans.

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.

Anatomically informed multi-level fiber tractography for targeted virtual dissection

OBJECTIVES

Diffusion-weighted MRI can assist preoperative planning by reconstructing the trajectory of eloquent fiber pathways, such as the corticospinal tract (CST). However, accurate reconstruction of the full extent of the CST remains challenging with existing tractography methods. We suggest a novel tractography algorithm exploiting unused fiber orientations to produce more complete and reliable results.

METHODS

Our novel approach, referred to as multi-level fiber tractography (MLFT), reconstructs fiber pathways by progressively considering previously unused fiber orientations at multiple levels of tract propagation. Anatomical priors are used to minimize the number of false-positive pathways. The MLFT method was evaluated on synthetic data and in vivo data by reconstructing the CST while compared to conventional tractography approaches.

RESULTS

The radial extent of MLFT reconstructions is comparable to that of probabilistic reconstruction: [Formula: see text] for the left and [Formula: see text] for the right hemisphere according to Wilcoxon test, while achieving significantly higher topography preservation compared to probabilistic tractography: [Formula: see text].

DISCUSSION

MLFT provides a novel way to reconstruct fiber pathways by adding the capability of including branching pathways in fiber tractography. Thanks to its robustness, feasible reconstruction extent and topography preservation, our approach may assist in clinical practice as well as in virtual dissection studies.

Disruptions in Structural and Functional Connectivity Relate to Poststroke Fatigue

Introduction: Poststroke fatigue (PSF) is a disabling condition with unclear etiology. The brain lesion is thought to be an important causal factor in PSF, although focal lesion characteristics such as size and location have not proven to be predictive. Given that the stroke lesion results not only in focal tissue death but also in widespread changes in brain networks that are structurally and functionally connected to damaged tissue, we hypothesized that PSF relates to disruptions in structural and functional connectivity. Materials and Methods: Twelve patients who incurred an ischemic stroke in the middle cerebral artery (MCA) territory 1-3 years prior, and currently experiencing a range of fatigue severity, were enrolled. The patients underwent structural and resting-state functional magnetic resonance imaging (MRI). The structural MRI data were used to measure structural disconnection of gray matter resulting from lesion to white matter pathways. The functional MRI data were used to measure network functional connectivity. Results: The patients showed structural disconnection in varying cortical and subcortical regions. Fatigue severity correlated significantly with structural disconnection of several frontal cortex regions in the ipsilesional (IL) and contralesional hemispheres. Fatigue-related structural disconnection was most severe in the IL rostral middle frontal cortex. Greater structural disconnection of a subset of fatigue-related frontal cortex regions, including the IL rostral middle frontal cortex, trended toward correlating significantly with greater loss in functional connectivity. Among identified fatigue-related frontal cortex regions, only the IL rostral middle frontal cortex showed loss in functional connectivity correlating significantly with fatigue severity. Conclusion: Our results provide evidence that loss in structural and functional connectivity of bihemispheric frontal cortex regions plays a role in PSF after MCA stroke, with connectivity disruptions of the IL rostral middle frontal cortex having a central role. Impact statement Poststroke fatigue (PSF) is a common disabling condition with unclear etiology. We hypothesized that PSF relates to disruptions in structural and functional connectivity secondary to the focal lesion. Using structural and resting-state functional connectivity magnetic resonance imaging (MRI) in patients with chronic middle cerebral artery (MCA) stroke, we found frontal cortex regions in the ipsilesional (IL) and contralesional hemispheres with greater structural disconnection correlating with greater fatigue. Among these fatigue-related cortices, the IL rostral middle frontal cortex showed loss in functional connectivity correlating with fatigue. These findings suggest that disruptions in structural and functional connectivity play a role in PSF after MCA stroke.

Short superficial white matter and aging: a longitudinal multi-site study of 1293 subjects and 2711 sessions

It is estimated that short association fibers running immediately beneath the cortex may make up as much as 60% of the total white matter volume. However, these have been understudied relative to the long-range association, projection, and commissural fibers of the brain. This is largely because of limitations of diffusion MRI fiber tractography, which is the primary methodology used to non-invasively study the white matter connections. Inspired by recent anatomical considerations and methodological improvements in superficial white matter (SWM) tractography, we aim to characterize changes in these fiber systems in cognitively normal aging, which provide insight into the biological foundation of age-related cognitive changes, and a better understanding of how age-related pathology differs from healthy aging. To do this, we used three large, longitudinal and cross-sectional datasets (N = 1293 subjects, 2711 sessions) to quantify microstructural features and length/volume features of several SWM systems. We find that axial, radial, and mean diffusivities show positive associations with age, while fractional anisotropy has negative associations with age in SWM throughout the entire brain. These associations were most pronounced in the frontal, temporal, and temporoparietal regions. Moreover, measures of SWM volume and length decrease with age in a heterogenous manner across the brain, with different rates of change in inter-gyri and intra-gyri SWM, and at slower rates than well-studied long-range white matter pathways. These features, and their variations with age, provide the background for characterizing normal aging, and, in combination with larger association pathways and gray matter microstructural features, may provide insight into fundamental mechanisms associated with aging and cognition.

Karawun: a software package for assisting evaluation of advances in multimodal imaging for neurosurgical planning and intraoperative neuronavigation

PURPOSE

The neuroimaging research community-which includes a broad range of scientific, medical, statistical, and engineering disciplines-has developed many tools to advance our knowledge of brain structure, function, development, aging, and disease. Past research efforts have clearly shaped clinical practice. However, translation of new methodologies into clinical practice is challenging. Anything that can reduce these barriers has the potential to improve the rate at which research outcomes can contribute to clinical practice. In this article, we introduce Karawun, a file format conversion tool, that has become a key part of our work in translating advances in diffusion imaging acquisition and analysis into neurosurgical practice at our institution.

METHODS

Karawun links analysis workflows created using open-source neuroimaging software, to Brainlab (Brainlab AG, Munich, Germany), a commercially available surgical planning and navigation suite. Karawun achieves this using DICOM standards supporting representation of 3D structures, including tractography streamlines, and thus offers far more than traditional screenshot or color overlay approaches.

RESULTS

We show that neurosurgical planning data, created from multimodal imaging data using analysis methods implemented in open-source research software, can be imported into Brainlab. The datasets can be manipulated as if they were created by Brainlab, including 3D visualizations of white matter tracts and other objects.

CONCLUSION

Clinicians can explore and interact with the results of research neuroimaging pipelines using familiar tools within their standard clinical workflow, understand the impact of the new methods on their practice and provide feedback to methods developers. This capability has been important to the translation of advanced analysis techniques into practice at our institution.

Structural and functional connectivity of the inferior temporal numeral area

A growing body of evidence suggests that in adults, there is a spatially consistent "inferior temporal numeral area" (ITNA) in the occipitotemporal cortex that appears to preferentially process Arabic digits relative to non-numerical symbols and objects. However, very little is known about why the ITNA is spatially segregated from regions that process other orthographic stimuli such as letters, and why it is spatially consistent across individuals. In the present study, we used diffusion-weighted imaging and functional magnetic resonance imaging to contrast structural and functional connectivity between left and right hemisphere ITNAs and a left hemisphere letter-preferring region. We found that the left ITNA had stronger structural and functional connectivity than the letter region to inferior parietal regions involved in numerical magnitude representation and arithmetic. Between hemispheres, the left ITNA showed stronger structural connectivity with the left inferior frontal gyrus (Broca's area), while the right ITNA showed stronger structural connectivity to the ipsilateral inferior parietal cortex and stronger functional coupling with the bilateral IPS. Based on their relative connectivity, our results suggest that the left ITNA may be more readily involved in mapping digits to verbal number representations, while the right ITNA may support the mapping of digits to quantity representations.

TractoInferno - A large-scale, open-source, multi-site database for machine learning dMRI tractography

TractoInferno is the world's largest open-source multi-site tractography database, including both research- and clinical-like human acquisitions, aimed specifically at machine learning tractography approaches and related ML algorithms. It provides 284 samples acquired from 3 T scanners across 6 different sites. Available data includes T1-weighted images, single-shell diffusion MRI (dMRI) acquisitions, spherical harmonics fitted to the dMRI signal, fiber ODFs, and reference streamlines for 30 delineated bundles generated using 4 tractography algorithms, as well as masks needed to run tractography algorithms. Manual quality control was additionally performed at multiple steps of the pipeline. We showcase TractoInferno by benchmarking the learn2track algorithm and 5 variations of the same recurrent neural network architecture. Creating the TractoInferno database required approximately 20,000 CPU-hours of processing power, 200 man-hours of manual QC, 3,000 GPU-hours of training baseline models, and 4 Tb of storage, to produce a final database of 350 Gb. By providing a standardized training dataset and evaluation protocol, TractoInferno is an excellent tool to address common issues in machine learning tractography.

Structural Changes in the Arcuate Fasciculus and Recovery of Post-stroke Aphasia: A 6-Month Follow-up Study using Diffusion Tensor Imaging

BACKGROUND

Temporal changes in the structural connectivity of major language tracts after stroke and their contribution to aphasia recovery are unclear.

OBJECTIVE

To investigate longitudinal arcuate fasciculus (AF) integrity changes and their relationship with post-stroke aphasia recovery using diffusion tensor imaging (DTI).

METHODS

Thirty-five patients with aphasia due to first-ever left hemispheric stroke underwent the Korean version of the Western Aphasia Battery and DTI at 1- and 6-month post stroke onset. Fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD) of both AF tracts were analyzed to evaluate the temporal changes in tract integrity and determine the correlation between changes (Δ; follow-up - initial) in DTI parameters and language scores.

RESULTS

At 6 months post-stroke, the mean FA decreased, and mean MD and RD increased in both hemispheres; however, compared with mean AD observed after 1 month, the mean observed at 6 months increased only in the left hemisphere (P < .05). ΔFA of the left AF and proportional change in the aphasia quotient showed a significant positive correlation (r = 0.365, P = .031). No correlation was found between changes in the right AF parameters and language score. The group with increased FA in the left AF showed more significant language improvement than the group with decreased FA.

CONCLUSIONS

During the subacute stage, the integrity of AF decreased in both hemispheres in patients with aphasia, and the change in structural connectivity of the left AF was associated with language improvement.

A multi-scale probabilistic atlas of the human connectome

The human brain is a complex system that can be efficiently represented as a network of structural connectivity. Many imaging studies would benefit from such network information, which is not always available. In this work, we present a whole-brain multi-scale structural connectome atlas. This tool has been derived from a cohort of 66 healthy subjects imaged with optimal technology in the setting of the Human Connectome Project. From these data we created, using extensively validated diffusion-data processing, tractography and gray-matter parcellation tools, a multi-scale probabilistic atlas of the human connectome. In addition, we provide user-friendly and accessible code to match this atlas to individual brain imaging data to extract connection-specific quantitative information. This can be used to associate individual imaging findings, such as focal white-matter lesions or regional alterations, to specific connections and brain circuits. Accordingly, network-level consequences of regional changes can be analyzed even in absence of diffusion and tractography data. This method is expected to broaden the accessibility and lower the yield for connectome research.

Human brain structural connectivity matrices-ready for modelling

The human brain represents a complex computational system, the function and structure of which may be measured using various neuroimaging techniques focusing on separate properties of the brain tissue and activity. We capture the organization of white matter fibers acquired by diffusion-weighted imaging using probabilistic diffusion tractography. By segmenting the results of tractography into larger anatomical units, it is possible to draw inferences about the structural relationships between these parts of the system. This pipeline results in a structural connectivity matrix, which contains an estimate of connection strength among all regions. However, raw data processing is complex, computationally intensive, and requires expert quality control, which may be discouraging for researchers with less experience in the field. We thus provide brain structural connectivity matrices in a form ready for modelling and analysis and thus usable by a wide community of scientists. The presented dataset contains brain structural connectivity matrices together with the underlying raw diffusion and structural data, as well as basic demographic data of 88 healthy subjects.

Measurement(s)	Diffusion Weighted Imaging
Technology Type(s)	Magnetic Resonance Imaging
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Location	Czech Republic

Informative and Reliable Tract Segmentation for Preoperative Planning

Identifying white matter (WM) tracts to locate eloquent areas for preoperative surgical planning is a challenging task. Manual WM tract annotations are often used but they are time-consuming, suffer from inter- and intra-rater variability, and noise intrinsic to diffusion MRI may make manual interpretation difficult. As a result, in clinical practice direct electrical stimulation is necessary to precisely locate WM tracts during surgery. A measure of WM tract segmentation unreliability could be important to guide surgical planning and operations. In this study, we use deep learning to perform reliable tract segmentation in combination with uncertainty quantification to measure segmentation unreliability. We use a 3D U-Net to segment white matter tracts. We then estimate model and data uncertainty using test time dropout and test time augmentation, respectively. We use a volume-based calibration approach to compute representative predicted probabilities from the estimated uncertainties. In our findings, we obtain a Dice of ≈0.82 which is comparable to the state-of-the-art for multi-label segmentation and Hausdorff distance <10mm. We demonstrate a high positive correlation between volume variance and segmentation errors, which indicates a good measure of reliability for tract segmentation ad uncertainty estimation. Finally, we show that calibrated predicted volumes are more likely to encompass the ground truth segmentation volume than uncalibrated predicted volumes. This study is a step toward more informed and reliable WM tract segmentation for clinical decision-making.

Assessment of structural disconnections in gliomas: comparison of indirect and direct approaches

Gliomas are amongst the most common primary brain tumours in adults and are often associated with poor prognosis. Understanding the extent of white matter (WM) which is affected outside the tumoral lesion may be of paramount importance to explain cognitive deficits and the clinical progression of the disease. To this end, we explored both direct (i.e., tractography based) and indirect (i.e., atlas-based) approaches to quantifying WM structural disconnections in a cohort of 44 high- and low-grade glioma patients. While these methodologies have recently gained popularity in the context of stroke and other pathologies, to our knowledge, this is the first time they are applied in patients with brain tumours. More specifically, in this work, we present a quantitative comparison of the disconnection maps provided by the two methodologies by applying well-known metrics of spatial similarity, extension, and correlation. Given the important role the oedematous tissue plays in the physiopathology of tumours, we performed these analyses both by including and excluding it in the definition of the tumoral lesion. This was done to investigate possible differences determined by this choice. We found that direct and indirect approaches offer two distinct pictures of structural disconnections in patients affected by brain gliomas, presenting key differences in several regions of the brain. Following the outcomes of our analysis, we eventually discuss the strengths and pitfalls of these two approaches when applied in this critical field.

The future perspectives of psychiatric neurosurgery

The future of psychiatric neurosurgery can be viewed from two separate perspectives: the immediate future and the distant future. Both show promise, but the treatment strategy for mental diseases and the technology utilized during these separate periods will likely differ dramatically. It can be expected that the initial advancements will be built upon progress of neuroimaging and stereotactic targeting while surgical technology becomes adapted to patient-specific symptomatology and structural/functional imaging parameters. This individualized approach has already begun to show significant promise when applied to deep brain stimulation for treatment-resistant depression and obsessive-compulsive disorder. If effectiveness of these strategies is confirmed by well designed, double-blind, placebo-controlled clinical studies, further technological advances will continue into the distant future, and will likely involve precise neuromodulation at the cellular level, perhaps using wireless technology with or without closed-loop design. This approach, being theoretically less invasive and carrying less risk, may ultimately propel psychiatric neurosurgery to the forefront in the treatment algorithm of mental illness. Despite prominent development of non-invasive therapeutic options, such as stereotactic radiosurgery or transcranial magnetic resonance-guided focused ultrasound, chances are there will still be a need in surgical management of patients with most intractable psychiatric conditions.

Tracing in vivo the dorsal loop of the optic radiation: convergent perspectives from tractography and electrophysiology compared to a neuroanatomical ground truth

The temporo-parietal junction (TPJ) is a cortical area contributing to a multiplicity of visual, language-related, and cognitive functions. In line with this functional richness, also the organization of the underlying white matter is highly complex and includes several bundles. The few studies tackling the outcome and neurological burdens of surgical operations addressing TPJ document the presence of language disturbances and visual field damages, with the latter hardly recovered in time. This observation advocates for identifying and functionally monitoring the optic radiation (OR) bundles that cross the white matter below the TPJ. In the present study, we adopted a multimodal approach to address the anatomo-functional correlates of the OR's dorsal loop. In particular, we combined cadavers' dissection with tractographic and electrophysiological data collected in drug-resistant epileptic patients explored by stereoelectroencephalography (SEEG). Cadaver dissection allowed us to appreciate the course and topography of the dorsal loop. More surprisingly, both tractographic and electrophysiological observations converged on a unitary picture highly coherent with the data obtained by neuroanatomical observation. The combination of diverse and multimodal observations allows overcoming the limitations intrinsic to single methodologies, defining a unitary picture which makes it possible to investigate the dorsal loop both presurgically and at the individual patient level, ultimately contributing to limit the postsurgical damages. Notwithstanding, such a combined approach could serve as a model of investigation for future neuroanatomical inquiries tackling white matter fibers anatomy and function through SEEG-derived neurophysiological data.

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.

Brain Connectometry Changes in Space Travelers After Long-Duration Spaceflight

Humans undergo extreme physiological changes when subjected to long periods of weightlessness, and as we continue to become a space-faring species, it is imperative that we fully understand the physiological changes that occur in the human body, including the brain. In this study, we present findings of brain structural changes associated with long-duration spaceflight based on diffusion magnetic resonance imaging (dMRI) data. Twelve cosmonauts who spent an average of six months aboard the International Space Station (ISS) were scanned in an MRI scanner pre-flight, ten days after flight, and at a follow-up time point seven months after flight. We performed differential tractography, a technique that confines white matter fiber tracking to voxels showing microstructural changes. We found significant microstructural changes in several large white matter tracts, such as the corpus callosum, arcuate fasciculus, corticospinal, corticostriatal, and cerebellar tracts. This is the first paper to use fiber tractography to investigate which specific tracts exhibit structural changes after long-duration spaceflight and may direct future research to investigate brain functional and behavioral changes associated with these white matter pathways.

Assessing the Reliability of Template-Based Clustering for Tractography in Healthy Human Adults

Tractography is a non-invasive technique to investigate the brain’s structural pathways (also referred to as tracts) that connect different brain regions. A commonly used approach for identifying tracts is with template-based clustering, where unsupervised clustering is first performed on a template in order to label corresponding tracts in unseen data. However, the reliability of this approach has not been extensively studied. Here, an investigation into template-based clustering reliability was performed, assessing the output from two datasets: Human Connectome Project (HCP) and MyConnectome project. The effect of intersubject variability on template-based clustering reliability was investigated, as well as the reliability of both deep and superficial white matter tracts. Identified tracts were evaluated by assessing Euclidean distances from a dataset-specific tract average centroid, the volumetric overlap across corresponding tracts, and along-tract agreement of quantitative values. Further, two template-based techniques were employed to evaluate the reliability of different clustering approaches. Reliability assessment can increase the confidence of a tract identifying technique in future applications to study pathways of interest. The two different template-based approaches exhibited similar reliability for identifying both deep white matter tracts and the superficial white matter.

Modelling the perception of music in brain network dynamics

We analyze the influence of music in a network of FitzHugh-Nagumo oscillators with empirical structural connectivity measured in healthy human subjects. We report an increase of coherence between the global dynamics in our network and the input signal induced by a specific music song. We show that the level of coherence depends crucially on the frequency band. We compare our results with experimental data, which also describe global neural synchronization between different brain regions in the gamma-band range in a time-dependent manner correlated with musical large-scale form, showing increased synchronization just before transitions between different parts in a musical piece (musical high-level events). The results also suggest a separation in musical form-related brain synchronization between high brain frequencies, associated with neocortical activity, and low frequencies in the range of dance movements, associated with interactivity between cortical and subcortical regions.

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.

White Matter Microstructural Differences in Youth With Classical Congenital Adrenal Hyperplasia

CONTEXT

Gray matter morphology in the prefrontal cortex and subcortical regions, including the hippocampus and amygdala, are affected in youth with classical congenital adrenal hyperplasia (CAH). It remains unclear if white matter connecting these aforementioned brain regions is compromised in youth with CAH.

OBJECTIVE

To examine brain white matter microstructure in youth with CAH compared to controls.

DESIGN

A cross-sectional sample of 23 youths with CAH due to 21-hydroxylase deficiency (12.9 ± 3.5 year; 61% female) and 33 healthy controls (13.1 ± 2.8 year; 61% female) with 3T multishell diffusion-weighted magnetic resonance brain scans.

MAIN OUTCOME MEASURES

Complementary modeling approaches, including diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI), to examine in vivo white matter microstructure in six white matter tracts that innervate the prefrontal and subcortical regions.

RESULTS

DTI showed CAH youth had lower fractional anisotropy in both the fornix and stria terminalis and higher mean diffusivity in the fornix compared to controls. NODDI modeling revealed that CAH youth have a significantly higher orientation dispersion index in the stria terminalis compared to controls. White matter microstructural integrity was associated with smaller hippocampal and amygdala volumes in CAH youth.

CONCLUSIONS

These patterns of microstructure reflect less restricted water diffusion likely due to less coherency in oriented microstructure. These results suggest that white matter microstructural integrity in the fornix and stria terminalis is compromised and may be an additional related brain phenotype alongside affected hippocampus and amygdala neurocircuitry in individuals with CAH.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Results of the 2020 fastMRI Challenge for Machine Learning MR Image Reconstruction

Accelerating MRI scans is one of the principal outstanding problems in the MRI research community. Towards this goal, we hosted the second fastMRI competition targeted towards reconstructing MR images with subsampled k-space data. We provided participants with data from 7,299 clinical brain scans (de-identified via a HIPAA-compliant procedure by NYU Langone Health), holding back the fully-sampled data from 894 of these scans for challenge evaluation purposes. In contrast to the 2019 challenge, we focused our radiologist evaluations on pathological assessment in brain images. We also debuted a new Transfer track that required participants to submit models evaluated on MRI scanners from outside the training set. We received 19 submissions from eight different groups. Results showed one team scoring best in both SSIM scores and qualitative radiologist evaluations. We also performed analysis on alternative metrics to mitigate the effects of background noise and collected feedback from the participants to inform future challenges. Lastly, we identify common failure modes across the submissions, highlighting areas of need for future research in the MRI reconstruction community.