Limits to anatomical accuracy of diffusion tractography using modern approaches

A reproducible pipeline for parcellation of the anterior limb of the internal capsule

BACKGROUND

The anterior limb of the internal capsule (ALIC) is a white matter structure connecting the prefrontal cortex (PFC) to the brainstem, thalamus, and subthalamic nucleus. It is a target for deep brain stimulation (DBS) for obsessive-compulsive disorder. There is strong interest in improving DBS targeting by using diffusion tractography to reconstruct and target specific ALIC fiber pathways, but this methodology is susceptible to errors and lacks validation. To address these limitations, we developed a novel diffusion tractography pipeline that generates reliable and biologically validated ALIC white matter reconstructions.

METHODS

Following algorithm development and refinement, we analyzed 43 control subjects each with 2 sets of 3T MRI data and a subset of 5 controls with 7T data from the Human Connectome Project. We generated 22 segmented ALIC fiber bundles (11 per hemisphere) based on prefrontal PFC regions of interest, and we analyzed the relationships among bundles.

RESULTS

We successfully reproduced the topographies established by prior anatomical work using images acquired at both 3T and 7T. Quantitative assessment demonstrated significantly smaller intra-subject variability relative to inter-subject variability for both test and retest groups across all but one PFC region. We examined the overlap between fibers from different PFC regions and a response tract for obsessive-compulsive disorder deep brain stimulation, and we reconstructed the PFC hyperdirect pathway using a modified version of our pipeline.

DISCUSSION

Our dMRI algorithm reliably generates biologically validated ALIC white matter reconstructions, allowing for more precise modelling of fibers for neuromodulation therapies.

A multi-modal, asymmetric, weighted, and signed description of anatomical connectivity

The macroscale connectome is the network of physical, white-matter tracts between brain areas. The connections are generally weighted and their values interpreted as measures of communication efficacy. In most applications, weights are either assigned based on imaging features-e.g. diffusion parameters-or inferred using statistical models. In reality, the ground-truth weights are unknown, motivating the exploration of alternative edge weighting schemes. Here, we explore a multi-modal, regression-based model that endows reconstructed fiber tracts with directed and signed weights. We find that the model fits observed data well, outperforming a suite of null models. The estimated weights are subject-specific and highly reliable, even when fit using relatively few training samples, and the networks maintain a number of desirable features. In summary, we offer a simple framework for weighting connectome data, demonstrating both its ease of implementation while benchmarking its utility for typical connectome analyses, including graph theoretic modeling and brain-behavior associations.

A review of diffusion MRI in mood disorders: mechanisms and predictors of treatment response

By measuring the molecular diffusion of water molecules in brain tissue, diffusion MRI (dMRI) provides unique insight into the microstructure and structural connections of the brain in living subjects. Since its inception, the application of dMRI in clinical research has expanded our understanding of the possible biological bases of psychiatric disorders and successful responses to different therapeutic interventions. Here, we review the past decade of diffusion imaging-based investigations with a specific focus on studies examining the mechanisms and predictors of therapeutic response in people with mood disorders. We present a brief overview of the general application of dMRI and key methodological developments in the field that afford increasingly detailed information concerning the macro- and micro-structural properties and connectivity patterns of white matter (WM) pathways and their perturbation over time in patients followed prospectively while undergoing treatment. This is followed by a more in-depth summary of particular studies using dMRI approaches to examine mechanisms and predictors of clinical outcomes in patients with unipolar or bipolar depression receiving pharmacological, neurostimulation, or behavioral treatments. Limitations associated with dMRI research in general and with treatment studies in mood disorders specifically are discussed, as are directions for future research. Despite limitations and the associated discrepancies in findings across individual studies, evolving research strongly indicates that the field is on the precipice of identifying and validating dMRI biomarkers that could lead to more successful personalized treatment approaches and could serve as targets for evaluating the neural effects of novel treatments.

Multimodal cross-examination of progressive apraxia of speech by diffusion tensor imaging-based tractography and Tau-PET scans

Progressive apraxia of speech (PAOS) is a 4R tauopathy characterized by difficulties with motor speech planning. Neurodegeneration in PAOS targets the premotor cortex, particularly the supplementary motor area (SMA), with degeneration of white matter (WM) tracts connecting premotor and motor cortices and Broca's area observed on diffusion tensor imaging (DTI). We aimed to assess flortaucipir uptake across speech-language-related WM tracts identified using DTI tractography in PAOS. Twenty-two patients with PAOS and 26 matched healthy controls were recruited by the Neurodegenerative Research Group (NRG) and underwent MRI and flortaucipir-PET. The patient population included patients with primary progressive apraxia of speech (PPAOS) and non-fluent variant/agrammatic primary progressive aphasia (agPPA). Flortaucipir PET scans and DTI were coregistered using rigid registration with a mutual information cost function in subject space. Alignments between DTI and flortaucipir PET were inspected in all cases. Whole-brain tractography was calculated using deterministic algorithms by a tractography reconstruction tool (DSI-studio) and specific tracts were identified using an automatic fiber tracking atlas-based method. Fractional anisotropy (FA) and flortaucipir standardized uptake value ratios (SUVRs) were averaged across the frontal aslant tract, arcuate fasciculi, inferior frontal-occipital fasciculus, inferior and middle longitudinal fasciculi, as well as the SMA commissural fibers. Reduced FA (p < .0001) and elevated flortaucipir SUVR (p = .0012) were observed in PAOS cases compared to controls across all combined WM tracts. For flortaucipir SUVR, the greatest differentiation of PAOS from controls was achieved with the SMA commissural fibers (area under the receiver operator characteristic curve [AUROC] = 0.83), followed by the left arcuate fasciculus (AUROC = 0.75) and left frontal aslant tract (AUROC = 0.71). Our findings demonstrate that flortaucipir uptake is increased across WM tracts related to speech/language difficulties in PAOS.

Structural underpinnings and long-term effects of resilience in Parkinson's disease

Resilience in neuroscience generally refers to an individual's capacity to counteract the adverse effects of a neuropathological condition. While resilience mechanisms in Alzheimer's disease are well-investigated, knowledge regarding its quantification, neurobiological underpinnings, network adaptations, and long-term effects in Parkinson's disease is limited. Our study involved 151 Parkinson's patients from the Parkinson's Progression Marker Initiative Database with available Magnetic Resonance Imaging, Dopamine Transporter Single-Photon Emission Computed Tomography scans, and clinical information. We used an improved prediction model linking neuropathology to symptom severity to estimate individual resilience levels. Higher resilience levels were associated with a more active lifestyle, increased grey matter volume in motor-associated regions, a distinct structural connectivity network and maintenance of relative motor functioning for up to a decade. Overall, the results indicate that relative maintenance of motor function in Parkinson's patients may be associated with greater neuronal substrate, allowing higher tolerance against neurodegenerative processes through dynamic network restructuring.

Structural-functional brain network coupling predicts human cognitive ability

Individual differences in general cognitive ability (GCA) have a biological basis within the structure and function of the human brain. Network neuroscience investigations revealed neural correlates of GCA in structural as well as in functional brain networks. However, whether the relationship between structural and functional networks, the structural-functional brain network coupling (SC-FC coupling), is related to individual differences in GCA remains an open question. We used data from 1030 adults of the Human Connectome Project, derived structural connectivity from diffusion weighted imaging, functional connectivity from resting-state fMRI, and assessed GCA as a latent g-factor from 12 cognitive tasks. Two similarity measures and six communication measures were used to model possible functional interactions arising from structural brain networks. SC-FC coupling was estimated as the degree to which these measures align with the actual functional connectivity, providing insights into different neural communication strategies. At the whole-brain level, higher GCA was associated with higher SC-FC coupling, but only when considering path transitivity as neural communication strategy. Taking region-specific variations in the SC-FC coupling strategy into account and differentiating between positive and negative associations with GCA, allows for prediction of individual cognitive ability scores in a cross-validated prediction framework (correlation between predicted and observed scores: r = 0.25, p < .001). The same model also predicts GCA scores in a completely independent sample (N = 567, r = 0.19, p < .001). Our results propose structural-functional brain network coupling as a neurobiological correlate of GCA and suggest brain region-specific coupling strategies as neural basis of efficient information processing predictive of cognitive ability.

Fibre orientation atlas guided rapid segmentation of white matter tracts

Fibre tract delineation from diffusion magnetic resonance imaging (MRI) is a valuable clinical tool for neurosurgical planning and navigation, as well as in research neuroimaging pipelines. Several popular methods are used for this task, each with different strengths and weaknesses making them more or less suited to different contexts. For neurosurgical imaging, priorities include ease of use, computational efficiency, robustness to pathology and ability to generalise to new tracts of interest. Many existing methods use streamline tractography, which may require expert neuroimaging operators for setting parameters and delineating anatomical regions of interest, or suffer from as a lack of generalisability to clinical scans involving deforming tumours and other pathologies. More recently, data-driven approaches including deep-learning segmentation models and streamline clustering methods have improved reproducibility and automation, although they can require large amounts of training data and/or computationally intensive image processing at the point of application. We describe an atlas-based direct tract mapping technique called 'tractfinder', utilising tract-specific location and orientation priors. Our aim was to develop a clinically practical method avoiding streamline tractography at the point of application while utilising prior anatomical knowledge derived from only 10-20 training samples. Requiring few training samples allows emphasis to be placed on producing high quality, neuro-anatomically accurate training data, and enables rapid adaptation to new tracts of interest. Avoiding streamline tractography at the point of application reduces computational time, false positives and vulnerabilities to pathology such as tumour deformations or oedema. Carefully filtered training streamlines and track orientation distribution mapping are used to construct tract specific orientation and spatial probability atlases in standard space. Atlases are then transformed to target subject space using affine registration and compared with the subject's voxel-wise fibre orientation distribution data using a mathematical measure of distribution overlap, resulting in a map of the tract's likely spatial distribution. This work includes extensive performance evaluation and comparison with benchmark techniques, including streamline tractography and the deep-learning method TractSeg, in two publicly available healthy diffusion MRI datasets (from TractoInferno and the Human Connectome Project) in addition to a clinical dataset comprising paediatric and adult brain tumour scans. Tract segmentation results display high agreement with established techniques while requiring less than 3 min on average when applied to a new subject. Results also display higher robustness than compared methods when faced with clinical scans featuring brain tumours and resections. As well as describing and evaluating a novel proposed tract delineation technique, this work continues the discussion on the challenges surrounding the white matter segmentation task, including issues of anatomical definitions and the use of quantitative segmentation comparison metrics.

Tractography with T1-weighted MRI and associated anatomical constraints on clinical quality diffusion MRI

Diffusion MRI (dMRI) streamline tractography, the gold-standard for in vivo estimation of white matter (WM) pathways in the brain, has long been considered as a product of WM microstructure. However, recent advances in tractography demonstrated that convolutional recurrent neural networks (CoRNN) trained with a teacher-student framework have the ability to learn to propagate streamlines directly from T1 and anatomical context. Training for this network has previously relied on high resolution dMRI. In this paper, we generalize the training mechanism to traditional clinical resolution data, which allows generalizability across sensitive and susceptible study populations. We train CoRNN on a small subset of the Baltimore Longitudinal Study of Aging (BLSA), which better resembles clinical scans. We define a metric, termed the epsilon ball seeding method, to compare T1 tractography and traditional diffusion tractography at the streamline level. We show that under this metric T1 tractography generated by CoRNN reproduces diffusion tractography with approximately three millimeters of error.

Diffusion tensor imaging-based multi-fiber tracking reconstructions can regionally differentiate phonetic versus prosodic subtypes of progressive apraxia of speech

Two subtypes of progressive apraxia of speech (PAOS) have been recognized: phonetic PAOS (PAOS_ph) where speech output is dominated by distorted sound substitutions and prosodic PAOS (PAOS_pr) which is dominated by segmented speech. We investigate whether these PAOS subtypes have different white matter microstructural abnormalities measured by diffusion tensor tractography. Thirty-three patients with PAOS (21 PAOS_ph and 12 PAOS_pr) and 19 healthy controls were recruited by the Neurodegenerative Research Group (NRG) and underwent diffusion MRI. Using a whole-brain tractography approach, fractional anisotropy (FA) and mean diffusivity (MD) were extracted for cortico-cortical, cortico-subcortical, cortical-projection, and cerebello-cortical white matter tracts. A hierarchical linear model was applied to assess tract-level FA and MD across groups. Both PAOS_ph and PAOS_pr showed degeneration of cortico-cortical, cortico-subcortical, cortical-projection, and cerebello-cortical white matter tracts compared to controls. However, degeneration of the body of corpus callosum, superior thalamic radiation, and superior cerebellar peduncle was greater in PAOS_pr compared to PAOS_ph, and degeneration of the inferior segment of the superior longitudinal fasciculus (SLF) was greater in PAOS_ph compared to PAOS_pr. Worse parkinsonism correlated with greater degeneration of cortico-cortical and cortico-subcortical tracts in PAOS_ph. Apraxia of speech articulatory error score correlated with degeneration of the superior cerebellar peduncle tracts in PAOS_pr. Phonetic and prosodic PAOS involve the compromise of a similar network of tracts, although there are connectivity differences between types. Whereas clinical parameters are the current gold standard to distinguish PAOS subtypes, our results allege the use of DTI-based tractography as a supplementary method to investigate such variants.

Resting-state functional magnetic resonance imaging indices are related to electrophysiological dysfunction in degenerative cervical myelopathy

The age-related degenerative pathologies of the cervical spinal column that comprise degenerative cervical myelopathy (DCM) cause myelopathy due spinal cord compression. Functional neurological assessment of DCM can potentially reveal the severity and pathological mechanism of DCM. However, functional assessment by conventional MRI remains difficult. This study used resting-state functional MRI (rs-fMRI) to investigate the relationship between functional connectivity (FC) strength and neurophysiological indices and examined the feasibility of functional assessment by FC for DCM. Preoperatively, 34 patients with DCM underwent rs-fMRI scans. Preoperative central motor conduction time (CMCT) reflecting motor functional disability and intraoperative somatosensory evoked potentials (SEP) reflecting sensory functional disability were recorded as electrophysiological indices of severity of the cervical spinal cord impairment. We performed seed-to-voxel FC analysis and correlation analyses between FC strength and the two electrophysiological indices. We found that FC strength between the primary motor cortex and the precuneus correlated significantly positively with CMCT, and that between the lateral part of the sensorimotor cortex and the lateral occipital cortex also showed a significantly positive correlation with SEP amplitudes. These results suggest that we can evaluate neurological and electrophysiological severity in patients with DCM by analyzing FC strengths between certain brain regions.

Advancements in Diffusion MRI Tractography for Neurosurgery

ABSTRACT

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

Bridging functional and anatomical neural connectivity through cluster synchronization

The dynamics of the brain results from the complex interplay of several neural populations and is affected by both the individual dynamics of these areas and their connection structure. Hence, a fundamental challenge is to derive models of the brain that reproduce both structural and functional features measured experimentally. Our work combines neuroimaging data, such as dMRI, which provides information on the structure of the anatomical connectomes, and fMRI, which detects patterns of approximate synchronous activity between brain areas. We employ cluster synchronization as a tool to integrate the imaging data of a subject into a coherent model, which reconciles structural and dynamic information. By using data-driven and model-based approaches, we refine the structural connectivity matrix in agreement with experimentally observed clusters of brain areas that display coherent activity. The proposed approach leverages the assumption of homogeneous brain areas; we show the robustness of this approach when heterogeneity between the brain areas is introduced in the form of noise, parameter mismatches, and connection delays. As a proof of concept, we apply this approach to MRI data of a healthy adult at resting state.

Supervised tractogram filtering using Geometric Deep Learning

A tractogram is a virtual representation of the brain white matter. It is composed of millions of virtual fibers, encoded as 3D polylines, which approximate the white matter axonal pathways. To date, tractograms are the most accurate white matter representation and thus are used for tasks like presurgical planning and investigations of neuroplasticity, brain disorders, or brain networks. However, it is a well-known issue that a large portion of tractogram fibers is not anatomically plausible and can be considered artifacts of the tracking procedure. With Verifyber, we tackle the problem of filtering out such non-plausible fibers using a novel fully-supervised learning approach. Differently from other approaches based on signal reconstruction and/or brain topology regularization, we guide our method with the existing anatomical knowledge of the white matter. Using tractograms annotated according to anatomical principles, we train our model, Verifyber, to classify fibers as either anatomically plausible or non-plausible. The proposed Verifyber model is an original Geometric Deep Learning method that can deal with variable size fibers, while being invariant to fiber orientation. Our model considers each fiber as a graph of points, and by learning features of the edges between consecutive points via the proposed sequence Edge Convolution, it can capture the underlying anatomical properties. The output filtering results highly accurate and robust across an extensive set of experiments, and fast; with a 12GB GPU, filtering a tractogram of 1M fibers requires less than a minute.

Reconstructing the somatotopic organization of the corticospinal tract remains a challenge for modern tractography methods

The corticospinal tract (CST) is a critically important white matter fiber tract in the human brain that enables control of voluntary movements of the body. The CST exhibits a somatotopic organization, which means that the motor neurons that control specific body parts are arranged in order within the CST. Diffusion magnetic resonance imaging (MRI) tractography is increasingly used to study the anatomy of the CST. However, despite many advances in tractography algorithms over the past decade, modern, state-of-the-art methods still face challenges. In this study, we compare the performance of six widely used tractography methods for reconstructing the CST and its somatotopic organization. These methods include constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, unscented Kalman filter (UKF) tractography methods including multi-fiber (UKF2T) and single-fiber (UKF1T) models, the generalized q-sampling imaging (GQI) based deterministic tractography method, and the TractSeg method. We investigate CST somatotopy by dividing the CST into four subdivisions per hemisphere that originate in the leg, trunk, hand, and face areas of the primary motor cortex. A quantitative and visual comparison is performed using diffusion MRI data (N = 100 subjects) from the Human Connectome Project. Quantitative evaluations include the reconstruction rate of the eight anatomical subdivisions, the percentage of streamlines in each subdivision, and the coverage of the white matter-gray matter (WM-GM) interface. CST somatotopy is further evaluated by comparing the percentage of streamlines in each subdivision to the cortical volumes for the leg, trunk, hand, and face areas. Overall, UKF2T has the highest reconstruction rate and cortical coverage. It is the only method with a significant positive correlation between the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex. However, our experimental results show that all compared tractography methods are biased toward generating many trunk streamlines (ranging from 35.10% to 71.66% of total streamlines across methods). Furthermore, the coverage of the WM-GM interface in the largest motor area (face) is generally low (under 40%) for all compared tractography methods. Different tractography methods give conflicting results regarding the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex, indicating that there is generally no clear relationship, and that reconstruction of CST somatotopy is still a large challenge. Overall, we conclude that while current tractography methods have made progress toward the well-known challenge of improving the reconstruction of the lateral projections of the CST, the overall problem of performing a comprehensive CST reconstruction, including clinically important projections in the lateral (hand and face areas) and medial portions (leg area), remains an important challenge for diffusion MRI tractography.

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.

Connectome-based modelling of neurodegenerative diseases: towards precision medicine and mechanistic insight

Neurodegenerative diseases are the most common cause of dementia. Although their underlying molecular pathologies have been identified, there is substantial heterogeneity in the patterns of progressive brain alterations across and within these diseases. Recent advances in neuroimaging methods have revealed that pathological proteins accumulate along specific macroscale brain networks, implicating the network architecture of the brain in the system-level pathophysiology of neurodegenerative diseases. However, the extent to which 'network-based neurodegeneration' applies across the wide range of neurodegenerative disorders remains unclear. Here, we discuss the state-of-the-art of neuroimaging-based connectomics for the mapping and prediction of neurodegenerative processes. We review findings supporting brain networks as passive conduits through which pathological proteins spread. As an alternative view, we also discuss complementary work suggesting that network alterations actively modulate the spreading of pathological proteins between connected brain regions. We conclude this Perspective by proposing an integrative framework in which connectome-based models can be advanced along three dimensions of innovation: incorporating parameters that modulate propagation behaviour on the basis of measurable biological features; building patient-tailored models that use individual-level information and allowing model parameters to interact dynamically over time. We discuss promises and pitfalls of these strategies for improving disease insights and moving towards precision medicine.

Structural disconnection relates to functional changes after temporal lobe epilepsy surgery

Epilepsy surgery consists of surgical resection of the epileptic focus and is recommended for patients with drug-resistant focal epilepsy. However, focal brain lesions can lead to effects in distant brain regions. Similarly, the focal resection in temporal lobe epilepsy surgery has been shown to lead to functional changes distant from the resection. Here we hypothesize that there are changes in brain function caused by temporal lobe epilepsy surgery in regions distant from the resection that are due to their structural disconnection from the resected epileptic focus. Therefore, the goal of this study was to localize changes in brain function caused by temporal lobe epilepsy surgery and relate them to the disconnection from the resected epileptic focus. This study takes advantage of the unique opportunity that epilepsy surgery provides to investigate the effects of focal disconnections on brain function in humans, which has implications in epilepsy and broader neuroscience. Changes in brain function from pre- to post-epilepsy surgery were quantified in a group of temporal lobe epilepsy patients (n = 36) using a measure of resting state functional MRI activity fluctuations. We identified regions with significant functional MRI changes that had high structural connectivity to the resected region in healthy controls (n = 96) and patients based on diffusion MRI. The structural disconnection from the resected epileptic focus was then estimated using presurgical diffusion MRI and related to the functional MRI changes from pre- to post-surgery in these regions. Functional MRI activity fluctuations increased from pre- to post-surgery in temporal lobe epilepsy in the two regions most highly structurally connected to the resected epileptic focus in healthy controls and patients-the thalamus and the fusiform gyrus ipsilateral to the side of surgery (PFWE < 0.05). Broader surgeries led to larger functional MRI changes in the thalamus than more selective surgeries (P < 0.05), but no other clinical variables were related to functional MRI changes in either the thalamus or fusiform. The magnitude of the functional MRI changes in both the thalamus and fusiform increased with a higher estimated structural disconnection from the resected epileptic focus when controlling for the type of surgery (P < 0.05). These results suggest that the structural disconnection from the resected epileptic focus may contribute to the functional changes seen after epilepsy surgery. Broadly, this study provides a novel link between focal disconnections in the structural brain network and downstream effects on function in distant brain regions.

Development of an educational method to rethink and learn oncological brain surgery in an "a la carte" connectome-based perspective

BACKGROUND

Understanding the structural connectivity of white matter tracts (WMT) and their related functions is a prerequisite to implementing an "a la carte" "connectomic approach" to glioma surgery. However, accessible resources facilitating such an approach are lacking. Here we present an educational method that is readily accessible, simple, and reproducible that enables the visualization of WMTs on individual patient images via an atlas-based approach.

METHODS

Our method uses the patient's own magnetic resonance imaging (MRI) images and consists of three main steps: data conversion, normalization, and visualization; these are accomplished using accessible software packages and WMT atlases. We implement our method on three common cases encountered in glioma surgery: a right supplementary motor area tumor, a left insular tumor, and a left temporal tumor.

RESULTS

Using patient-specific perioperative MRIs with open-sourced and co-registered atlas-derived WMTs, we highlight the critical subnetworks requiring specific surgical monitoring identified intraoperatively using direct electrostimulation mapping with cognitive monitoring. The aim of this didactic method is to provide the neurosurgical oncology community with an accessible and ready-to-use educational tool, enabling neurosurgeons to improve their knowledge of WMTs and to better learn their oncologic cases, especially in glioma surgery using awake mapping.

CONCLUSIONS

Taking no more than 3-5 min per patient and irrespective of their resource settings, we believe that this method will enable junior surgeons to develop an intuition, and a robust 3-dimensional imagery of WMT by regularly applying it to their cases both before and after surgery to develop an "a la carte" connectome-based perspective to glioma surgery.

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.

Statistical Learning Methods for Neuroimaging Data Analysis with Applications

The aim of this review is to provide a comprehensive survey of statistical challenges in neuroimaging data analysis, from neuroimaging techniques to large-scale neuroimaging studies and statistical learning methods. We briefly review eight popular neuroimaging techniques and their potential applications in neuroscience research and clinical translation. We delineate four themes of neuroimaging data and review major image processing analysis methods for processing neuroimaging data at the individual level. We briefly review four large-scale neuroimaging-related studies and a consortium on imaging genomics and discuss four themes of neuroimaging data analysis at the population level. We review nine major population-based statistical analysis methods and their associated statistical challenges and present recent progress in statistical methodology to address these challenges.

CACTUS: a computational framework for generating realistic white matter microstructure substrates

Monte-Carlo diffusion simulations are a powerful tool for validating tissue microstructure models by generating synthetic diffusion-weighted magnetic resonance images (DW-MRI) in controlled environments. This is fundamental for understanding the link between micrometre-scale tissue properties and DW-MRI signals measured at the millimetre-scale, optimizing acquisition protocols to target microstructure properties of interest, and exploring the robustness and accuracy of estimation methods. However, accurate simulations require substrates that reflect the main microstructural features of the studied tissue. To address this challenge, we introduce a novel computational workflow, CACTUS (Computational Axonal Configurator for Tailored and Ultradense Substrates), for generating synthetic white matter substrates. Our approach allows constructing substrates with higher packing density than existing methods, up to 95% intra-axonal volume fraction, and larger voxel sizes of up to 500μm3 with rich fibre complexity. CACTUS generates bundles with angular dispersion, bundle crossings, and variations along the fibres of their inner and outer radii and g-ratio. We achieve this by introducing a novel global cost function and a fibre radial growth approach that allows substrates to match predefined targeted characteristics and mirror those reported in histological studies. CACTUS improves the development of complex synthetic substrates, paving the way for future applications in microstructure imaging.

Role of interhemispheric connectivity in recovery from postoperative supplementary motor area syndrome in glioma patients

OBJECTIVE

Surgical resection of gliomas involving the supplementary motor area (SMA) frequently results in SMA syndrome, a symptom complex characterized by transient akinesia and mutism. Because the factors influencing patient functional outcomes after surgery remain elusive, the authors investigated network-based predictors in a multicentric cohort of glioma patients.

METHODS

The participants were 50 patients treated for glioma located in the SMA at one of the three centers participating in the study. Postoperative functional outcomes (motor deficits, mutism) and duration of symptoms were assessed during hospitalization. Long-term outcome was assessed 3 months after surgery. MRI-based lesion-symptom mapping was performed to estimate the severity of gray matter damage and white matter disconnection.

RESULTS

The median duration of acute symptoms was 3 days (range 1-42 days). Long-term deficits involving fine motor movements and speech were found at follow-up in 27 patients (54%). Disconnection of the central callosal fibers was associated with prolonged acute symptoms (p < 0.05). Postoperative mutism was significantly related to disconnection severity of the left frontopontine tract, frontal aslant tract, cingulum, and corticostriatal tract (p < 0.05). Disconnection of midposterior callosal fibers and lesion loads within the left medial Brodmann area 4 were associated with long-term motor deficits (p < 0.05).

CONCLUSIONS

This study provides evidence for the pathophysiology and predictive factors of postoperative SMA syndrome by demonstrating the relation of the disconnection of callosal fibers with prolonged symptom duration (central segment) and long-term motor deficits (midposterior segment). These data may be useful for presurgical risk assessment and adequate consultation for patients prior to undergoing resection of glioma located within the SMA region.

Assessing informative tract segmentation and nTMS for pre-operative planning

BACKGROUND

Deep learning-based (DL) methods are the best-performing methods for white matter tract segmentation in anatomically healthy subjects. However, tract annotations are variable or absent in clinical data and manual annotations are especially difficult in patients with tumors where normal anatomy may be distorted. Direct cortical and subcortical stimulation is the gold standard ground truth to determine the cortical and sub-cortical lo- cation of motor-eloquent areas intra-operatively. Nonetheless, this technique is invasive, prolongs the surgical procedure, and may cause patient fatigue. Navigated Transcranial Magnetic Stimulation (nTMS) has a well-established correlation to direct cortical stimulation for motor mapping and the added advantage of being able to be acquired pre-operatively.

NEW METHOD

In this work, we evaluate the feasibility of using nTMS motor responses as a method to assess corticospinal tract (CST) binary masks and estimated uncertainty generated by a DL-based tract segmentation in patients with diffuse gliomas.

RESULTS

Our results show CST binary masks have a high overlap coefficient (OC) with nTMS response masks. A strong negative correlation is found between estimated uncertainty and nTMS response mask distance to the CST binary mask.

COMPARISON WITH EXISTING METHODS

We compare our approach (UncSeg) with the state-of-the-art TractSeg in terms of OC between the CST binary masks and nTMS response masks.

CONCLUSIONS

In this study, we demonstrate that estimated uncertainty from UncSeg is a good measure of the agreement between the CST binary masks and nTMS response masks distance to the CST binary mask boundary.

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.

Exploring interhemispheric connectivity using the directional tract density patterns of the corpus callosum

The corpus callosum (CC) is one of the most important interhemispheric white matter tracts that connects interrelated regions of the cerebral cortex. Its disruption has been investigated in previous studies and has been found to play an important role in several neurodegenerative disorders. Currently available methods to assess the interhemispheric connectivity of the CC have several limitations: i) they require the a priori identification of specific cortical regions as targets or seeds, ii) they are limited by the characterization of only small components of the structure, primarily voxels that constitute the mid-sagittal slice, and iii) they use global measures of microstructural integrity, which provide only limited characterization. In order to address some of these limitations, we developed a novel method that enables the characterization of white matter tracts covering the structure of CC, from the mid-sagittal plane to corresponding regions of cortex, using directional tract density patterns (dTDPs). We demonstrate that different regions of CC have distinctive dTDPs that reflect a unique regional topology. We conducted a pilot study using this approach to evaluate two different datasets collected from healthy subjects, and we demonstrate that this method is reliable, reproducible, and independent of diffusion acquisition parameters, suggesting its potential applicability to clinical applications.

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

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

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

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

Comparison of Anatomical Pathway Models with Tractography Estimates of the Pallidothalamic, Cerebellothalamic, and Corticospinal Tracts

Introduction: Models of structural connectivity in the human brain are typically simulated using tractographic approaches. However, the nonlinear fitting of anatomical pathway atlases to de novo subject brains represents a simpler alternative that is hypothesized to provide more anatomically realistic results. Therefore, the goal of this study was to perform a side-by-side comparison of the streamline estimates generated by either pathway atlas fits or tractographic reconstructions in the same subjects. Methods: Our analyses focused on reconstruction of the corticospinal tract (CST), cerebellothalamic (CBT), and pallidothalamic (PT) pathways using example datasets from the Human Connectome Project (HCP). We used MRtrix3 to explore whole brain, as well as manual seed-to-target, tractography approaches. In parallel, we performed nonlinear fits of an axonal pathway atlas to each HCP dataset using Advanced Normalization Tools (ANTs). Results: The different methods produced notably different estimates for each pathway in each subject. The fitted atlas pathways were highly stereotyped and exhibited low variability in their streamline trajectories. Manual tractography resulted in pathway estimates that generally corresponded with the fitted atlas pathways, but with a higher degree of variability in the individual streamlines. Pathway reconstructions derived from whole-brain tractography exhibited the highest degree of variability and struggled to create anatomically realistic representations for either the CBT or PT pathways. Conclusion: The speed, simplicity, reproducibility, and realism of anatomical pathway model fits makes them an appealing option for some forms of structural connectivity modeling in the human brain. Impact statement Axonal pathway modeling is an important component of deep brain stimulation (DBS) research studies that seek to identify the brain connections that are directly activated by stimulation. The corticospinal tract, cerebellothalamic (CBT), and pallidothalamic (PT) pathways are specifically relevant to the study of subthalamic DBS for the treatment of Parkinson's disease. Our results suggest that anatomical pathway model fits of the CBT and PT pathways to de novo subject brains represent a more anatomically realistic option than tractographic approaches when studying subthalamic DBS.

Using mesoscopic tract-tracing data to guide the estimation of fiber orientation distributions in the mouse brain from diffusion MRI

Diffusion MRI (dMRI) tractography is the only tool for non-invasive mapping of macroscopic structural connectivity over the entire brain. Although it has been successfully used to reconstruct large white matter tracts in the human and animal brains, the sensitivity and specificity of dMRI tractography remained limited. In particular, the fiber orientation distributions (FODs) estimated from dMRI signals, key to tractography, may deviate from histologically measured fiber orientation in crossing fibers and gray matter regions. In this study, we demonstrated that a deep learning network, trained using mesoscopic tract-tracing data from the Allen Mouse Brain Connectivity Atlas, was able to improve the estimation of FODs from mouse brain dMRI data. Tractography results based on the network generated FODs showed improved specificity while maintaining sensitivity comparable to results based on FOD estimated using a conventional spherical deconvolution method. Our result is a proof-of-concept of how mesoscale tract-tracing data can guide dMRI tractography and enhance our ability to characterize brain connectivity.

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.

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.

Presurgical Executive Functioning in Low-Grade Glioma Patients Cannot Be Topographically Mapped

Executive dysfunctions have a high prevalence in low-grade glioma patients and may be the result of structural disconnections of particular subcortical tracts and/or networks. However, little research has focused on preoperative low-grade glioma patients. The frontotemporoparietal network has been closely linked to executive functions and is substantiated by the superior longitudinal fasciculus. The aim of this study was to investigate their role in executive functions in low-grade glioma patients. Patients from two neurological centers were included with IDH-mutated low-grade gliomas. The sets of preoperative predictors were (i) distance between the tumor and superior longitudinal fasciculus, (ii) structural integrity of the superior longitudinal fasciculus, (iii) overlap between tumor and cortical networks, and (iv) white matter disconnection of the same networks. Linear regression and random forest analyses were performed. The group of 156 patients demonstrated significantly lower performance than normative samples and had a higher prevalence of executive impairments. However, both regression and random forest analyses did not demonstrate significant results, meaning that neither structural, cortical network overlap, nor network disconnection predictors explained executive performance. Overall, our null results indicate that there is no straightforward topographical explanation of executive performance in low-grade glioma patients. We extensively discuss possible explanations, including plasticity-induced network-level equipotentiality. Finally, we stress the need for the development of novel methods to unveil the complex and interacting mechanisms that cause executive deficits in low-grade glioma patients.

Connectivity-based parcellation of the amygdala and identification of its main white matter connections

The amygdala plays a role in emotion, learning, and memory and has been implicated in behavioral disorders. Better understanding of the amygdala circuitry is crucial to develop new therapies for these disorders. We used data from 200 healthy-subjects from the human connectome project. Using probabilistic tractography, we created population statistical maps of amygdala connectivity to brain regions involved in limbic, associative, memory, and reward circuits. Based on the amygdala connectivity with these regions, we applied k-means clustering to parcellate the amygdala into three clusters. The resultant clusters were averaged across all subjects and the main white-matter pathways of the amygdala from each averaged cluster were generated. Amygdala parcellation into three clusters showed a medial-to-lateral pattern. The medial cluster corresponded with the centromedial and cortical nuclei, the basal cluster with the basal nuclei and the lateral cluster with the lateral nuclei. The connectivity analysis revealed different white-matter pathways consistent with the anatomy of the amygdala circuit. This in vivo connectivity-based parcellation of the amygdala delineates three clusters of the amygdala in a mediolateral pattern based on its connectivity with brain areas involved in cognition, memory, emotion, and reward. The human amygdala circuit presented in this work provides the first step for personalized amygdala circuit mapping for patients with behavioral disorders.

Diffusion tensor imaging of superficial prefrontal white matter in healthy aging

The prefrontal cortex (PFC) is a heterogenous structure that is highly susceptible to the effects of aging. Few studies have investigated age effects on the superficial white matter (WM) contained within the PFC using in-vivo magnetic resonance imaging (MRI). This study used diffusion tensor imaging (DTI) tractography to examine the effects of age, sex, and intracranial volume (ICV) on superficial WM within specific PFC subregions, and to model the relationships with age using higher order polynomial regression modelling. PFC WM of 140 healthy individuals, aged 18-85, was segmented into medial and lateral orbitofrontal, medial prefrontal, and dorsolateral prefrontal subregions. Differences due to age in microstructural parameters such as fractional anisotropy (FA), axial and radial diffusivities, and macrostructural measures of tract volumes, fiber counts, average fiber lengths, and average number of fibers per voxel were examined. We found that most prefrontal subregions demonstrated age effects, with decreases in FA, tract volume, and fiber counts, and increases in all diffusivity measures. Age relationships were mostly non-linear, with higher order regressions chosen in most cases. Declines in PFC FA began at the onset of adulthood while the greatest changes in diffusivity and volume did not occur until middle age. The effects of age were most prominent in medial tracts while the lateral orbitofrontal tracts were less affected. Significant effects of sex and ICV were also observed in certain parameters. The patterns mostly followed myelination order, with late-myelinating prefrontal subregions experiencing earlier and more pronounced age effects, further supporting the frontal theory of aging.

A structural connectivity atlas of limbic brainstem nuclei

Background

Understanding the structural connectivity of key brainstem nuclei with limbic cortical regions is essential to the development of therapeutic neuromodulation for depression, chronic pain, addiction, anxiety and movement disorders. Several brainstem nuclei have been identified as the primary central nervous system (CNS) source of important monoaminergic ascending fibers including the noradrenergic locus coeruleus, serotonergic dorsal raphe nucleus, and dopaminergic ventral tegmental area. However, due to practical challenges to their study, there is limited data regarding their in vivo anatomic connectivity in humans.

Objective

To evaluate the structural connectivity of the following brainstem nuclei with limbic cortical areas: locus coeruleus, ventral tegmental area, periaqueductal grey, dorsal raphe nucleus, and nucleus tractus solitarius. Additionally, to develop a group average atlas of these limbic brainstem structures to facilitate future analyses.

Methods

Each nucleus was manually masked from 197 Human Connectome Project (HCP) structural MRI images using FSL software. Probabilistic tractography was performed using FSL's FMRIB Diffusion Toolbox. Connectivity with limbic cortical regions was calculated and compared between brainstem nuclei. Results were aggregated to produce a freely available MNI structural atlas of limbic brainstem structures.

Results

A general trend was observed for a high probability of connectivity to the amygdala, hippocampus and DLPFC with relatively lower connectivity to the orbitofrontal cortex, NAc, hippocampus and insula. The locus coeruleus and nucleus tractus solitarius demonstrated significantly greater connectivity to the DLPFC than amygdala while the periaqueductal grey, dorsal raphe nucleus, and ventral tegmental area did not demonstrate a significant difference between these two structures.

Conclusion

Monoaminergic and other modulatory nuclei in the brainstem project widely to cortical limbic regions. We describe the structural connectivity across the several key brainstem nuclei theorized to influence emotion, reward, and cognitive functions. An increased understanding of the anatomic basis of the brainstem's role in emotion and other reward-related processing will support targeted neuromodulatary therapies aimed at alleviating the symptoms of neuropsychiatric disorders.

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.

Maximal safe resection of diffuse lower grade gliomas primarily within central lobe using cortical/subcortical direct electrical stimulation under awake craniotomy

Background

Diffuse lower-grade glioma (DLGG) in the central lobe is a challenge for safe resection procedures. To improve the extent of resection and reduce the risk of postoperative neurological deficits, we performed an awake craniotomy with cortical-subcortical direct electrical stimulation (DES) mapping for patients with DLGG located primarily within the central lobe. We investigated the outcomes of cortical-subcortical brain mapping using DES in an awake craniotomy for central lobe DLGG resection.

Methods

We performed a retrospective analysis of clinical data of a cohort of consecutively treated patients from February 2017 to August 2021 with diffuse lower-grade gliomas located primarily within the central lobe. All patients underwent awake craniotomy with DES for cortical and subcortical mapping of eloquent brain areas, neuronavigation, and/or ultrasound to identify tumor location. Tumors were removed according to functional boundaries. Maximum safe tumor resection was the surgical objective for all patients.

Results

Thirteen patients underwent 15 awake craniotomies with intraoperative mapping of eloquent cortices and subcortical fibers using DES. Maximum safe tumor resection was achieved according to functional boundaries in all patients. The pre-operative tumor volumes ranged from 4.3 cm³ to 137.3 cm³ (median 19.2 cm³). The mean extent of tumor resection was 94.6%, with eight cases (53.3%) achieving total resection, four (26.7%) subtotal and three (20.0%) partial. The mean tumor residue was 1.2 cm³. All patients experienced early postoperative neurological deficits or worsening conditions. Three patients (20.0%) experienced late postoperative neurological deficits at the 3-month follow-up, including one moderate and two mild neurological deficits. None of the patients experienced late onset severe neurological impairments post-operatively. Ten patients with 12 tumor resections (80.0%) had resumed activities of daily living at the 3-month follow-up. Among 14 patients with pre-operative epilepsy, 12 (85.7%) were seizure-free after treatment with antiepileptic drugs 7 days after surgery up to the last follow-up.

Conclusions

DLGG located primarily in the central lobe deemed inoperable can be safely resected using awake craniotomy with intraoperative DES without severe permanent neurological sequelae. Patients experienced an improved quality of life in terms of seizure control.

Tractography-based navigated TMS language mapping protocol

Introduction

This study explores the feasibility of implementing a tractography-based navigated transcranial magnetic stimulation (nTMS) language mapping protocol targeting cortical terminations of the arcuate fasciculus (AF). We compared the results and distribution of errors from the new protocol to an established perisylvian nTMS protocol that stimulated without any specific targeting over the entire perisylvian cortex.

Methods

Sixty right-handed patients with language-eloquent brain tumors were examined in this study with one half of the cohort receiving the tractographybased protocol and the other half receiving the perisylvian protocol. Probabilistic tractography using MRtrix3 was performed for patients in the tractography-based group to identify the AF's cortical endpoints. nTMS mappings were performed and resulting language errors were classified into five psycholinguistic groups.

Results

Tractography and nTMS were successfully performed in all patients. The tractogram-based group showed a significantly higher median overall ER than the perisylvian group (3.8% vs. 2.9% p <.05). The median ER without hesitation errors in the tractogram-based group was also significantly higher than the perisylvian group (2.0% vs. 1.4%, p <.05). The ERs by error type showed no significant differences between protocols except in the no response ER, with a higher median ER in the tractogram-based group (0.4% vs. 0%, p <.05). Analysis of ERs based on the Corina cortical parcellation system showed especially high nTMS ERs over the posterior middle temporal gyrus (pMTG) in the perisylvian protocol and high ERs over the middle and ventral postcentral gyrus (vPoG), the opercular inferior frontal gyrus (opIFG) and the ventral precentral gyrus (vPrG) in the tractography-based protocol.

Discussion

By considering the white matter anatomy and performing nTMS on the cortical endpoints of the AF, the efficacy of nTMS in disrupting patients' object naming abilities was increased. The newly introduced method showed proof of concept and resulted in AF-specific ERs and noninvasive cortical language maps, which could be applied to additional fiber bundles related to the language network in future nTMS studies.

Automation and standardization of subject-specific region-of-interest segmentation for investigation of diffusion imaging in clinical populations

Diffusion weighted imaging (DWI) is key in clinical neuroimaging studies. In recent years, DWI has undergone rapid evolution and increasing applications. Diffusion magnetic resonance imaging (dMRI) is widely used to analyse group-level differences in white matter (WM), but suffers from limitations that can be particularly impactful in clinical groups where 1) structural abnormalities may increase erroneous inter-subject registration and 2) subtle differences in WM microstructure between individuals can be missed. It also lacks standardization protocols for analyses at the subject level. Region of Interest (ROI) analyses in native diffusion space can help overcome these challenges, with manual segmentation still used as the gold standard. However, robust automated approaches for the analysis of ROI-extracted native diffusion characteristics are limited. Subject-Specific Diffusion Segmentation (SSDS) is an automated pipeline that uses pre-existing imaging analysis methods to carry out WM investigations in native diffusion space, while overcoming the need to interpolate diffusion images and using an intermediate T1 image to limit registration errors and guide segmentation. SSDS is validated in a cohort of healthy subjects scanned three times to derive test-retest reliability measures and compared to other methods, namely manual segmentation and tract-based spatial statistics as an example of group-level method. The performance of the pipeline is further tested in a clinical population of patients with traumatic brain injury and structural abnormalities. Mean FA values obtained from SSDS showed high test-retest and were similar to FA values estimated from the manual segmentation of the same ROIs (p-value > 0.1). The average dice similarity coefficients (DSCs) comparing results from SSDS and manual segmentations was 0.8 ± 0.1. Case studies of TBI patients showed robustness to the presence of significant structural abnormalities, indicating its potential clinical application in the identification and diagnosis of WM abnormalities. Further recommendation is given regarding the tracts used with SSDS.

A connectomics-based taxonomy of mammals

Mammalian taxonomies are conventionally defined by morphological traits and genetics. How species differ in terms of neural circuits and whether inter-species differences in neural circuit organization conform to these taxonomies is unknown. The main obstacle to the comparison of neural architectures has been differences in network reconstruction techniques, yielding species-specific connectomes that are not directly comparable to one another. Here, we comprehensively chart connectome organization across the mammalian phylogenetic spectrum using a common reconstruction protocol. We analyse the mammalian MRI (MaMI) data set, a database that encompasses high-resolution ex vivo structural and diffusion MRI scans of 124 species across 12 taxonomic orders and 5 superorders, collected using a unified MRI protocol. We assess similarity between species connectomes using two methods: similarity of Laplacian eigenspectra and similarity of multiscale topological features. We find greater inter-species similarities among species within the same taxonomic order, suggesting that connectome organization reflects established taxonomic relationships defined by morphology and genetics. While all connectomes retain hallmark global features and relative proportions of connection classes, inter-species variation is driven by local regional connectivity profiles. By encoding connectomes into a common frame of reference, these findings establish a foundation for investigating how neural circuits change over phylogeny, forging a link from genes to circuits to behaviour.

Multimodal analysis demonstrating the shaping of functional gradients in the marmoset brain

The discovery of functional gradients introduce a new perspective in understanding the cortical spectrum of intrinsic dynamics, as it captures major axes of functional connectivity in low-dimensional space. However, how functional gradients arise and dynamically vary remains poorly understood. In this study, we investigated the biological basis of functional gradients using awake resting-state fMRI, retrograde tracing and gene expression datasets in marmosets. We found functional gradients in marmosets showed a sensorimotor-to-visual principal gradient followed by a unimodal-to-multimodal gradient, resembling functional gradients in human children. Although strongly constrained by structural wirings, functional gradients were dynamically modulated by arousal levels. Utilizing a reduced model, we uncovered opposing effects on gradient dynamics by structural connectivity (inverted U-shape) and neuromodulatory input (U-shape) with arousal fluctuations, and dissected the contribution of individual neuromodulatory receptors. This study provides insights into biological basis of functional gradients by revealing the interaction between structural connectivity and ascending neuromodulatory system.

Tractography indicates lateralized differences between trigeminal and olfactory pathways

Odorous sensations are based on trigeminal and olfactory perceptions. Both trigeminal and olfactory stimuli generate overlapping as well as distinctive activations in the olfactory cortex including the piriform cortex. Orbitofrontal cortex (OFC), an integrative center for all senses, is directly activated in the presence of olfactory stimulations. In contrast, the thalamus, a very important midbrain structure, is not directly activated in the presence of odors, but rather acts as a relay for portions of olfactory information between primary olfactory cortex and higher-order processing centers. The aims of the study were (1) to examine the number of streamlines between the piriform cortex and the OFC and also between the piriform cortex and the thalamus and (2) to explore potential correlations between these streamlines and trigeminal and olfactory chemosensory perceptions. Thirty-eight healthy subjects were recruited for the study and underwent diffusion MRI using a 3T MRI scanner with 67 diffusion directions. ROIs were adapted from two studies looking into olfaction in terms of functional and structural properties of the olfactory system. The "waytotal number" was used which corresponds to number of streamlines between two regions of interests. We found the number of streamlines between the piriform cortex and the thalamus to be higher in the left hemisphere, whereas the number of streamlines between the piriform cortex and the OFC were higher in the right hemisphere. We also found streamlines between the piriform cortex and the thalamus to be positively correlated with the intensity of irritating (trigeminal) odors. On the other hand, streamlines between the piriform cortex and the OFC were correlated with the threshold scores for these trigeminal odors. This is the first studying the correlations between streamlines and olfactory scores using tractography. Results suggest that different chemosensory stimuli are processed through different networks in the chemosensory system involving the thalamus.

Concurrent mapping of brain ontogeny and phylogeny within a common space: Standardized tractography and applications

Developmental and evolutionary effects on brain organization are complex, yet linked, as evidenced by the correspondence in cortical area expansion across these vastly different time scales. However, it is still not possible to study concurrently the ontogeny and phylogeny of cortical areal connections, which is arguably more relevant to brain function than allometric measurements. Here, we propose a novel framework that allows the integration of structural connectivity maps from humans (adults and neonates) and nonhuman primates (macaques) onto a common space. We use white matter bundles to anchor the common space and use the uniqueness of cortical connection patterns to these bundles to probe area specialization. This enabled us to quantitatively study divergences and similarities in connectivity over evolutionary and developmental scales, to reveal brain maturation trajectories, including the effect of premature birth, and to translate cortical atlases between diverse brains. Our findings open new avenues for an integrative approach to imaging neuroanatomy.

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

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

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

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

A groupwise registration and tractography framework for cardiac myofiber architecture description by diffusion MRI: An application to the ventricular junctions

Background

Knowledge of the normal myocardial–myocyte orientation could theoretically allow the definition of relevant quantitative biomarkers in clinical routine to diagnose heart pathologies. A whole heart diffusion tensor template representative of the global myofiber organization over species is therefore crucial for comparisons across populations. In this study, we developed a groupwise registration and tractography framework to resolve the global myofiber arrangement of large mammalian sheep hearts. To demonstrate the potential application of the proposed method, a novel description of sub-regions in the intraventricular septum is presented.

Methods

Three explanted sheep (ovine) hearts (size ~12×8×6 cm3, heart weight ~ 150 g) were perfused with contrast agent and fixative and imaged in a 9.4T magnet. A group-wise registration of high-resolution anatomical and diffusion-weighted images were performed to generate anatomical and diffusion tensor templates. Diffusion tensor metrics (eigenvalues, eigenvectors, fractional anisotropy …) were computed to provide a quantitative and spatially-resolved analysis of cardiac microstructure. Then tractography was performed using deterministic and probabilistic algorithms and used for different purposes: i) Visualization of myofiber architecture, ii) Segmentation of sub-area depicting the same fiber organization, iii) Seeding and Tract Editing. Finally, dissection was performed to confirm the existence of macroscopic structures identified in the diffusion tensor template.

Results

The template creation takes advantage of high-resolution anatomical and diffusion-weighted images obtained at an isotropic resolution of 150 μm and 600 μm respectively, covering ventricles and atria and providing information on the normal myocardial architecture. The diffusion metric distributions from the template were found close to the one of the individual samples validating the registration procedure. Small new sub-regions exhibiting spatially sharp variations in fiber orientation close to the junctions of the septum and ventricles were identified. Each substructure was defined and represented using streamlines. The existence of a fiber-bundles in the posterior junction was validated by anatomical dissection. A complex structural organization of the anterior junction in comparison to the posterior junction was evidenced by the high-resolution acquisition.

Conclusions

A new framework combining cardiac template generation and tractography was applied on the whole sheep heart. The framework can be used for anatomical investigation, characterization of microstructure and visualization of myofiber orientation across samples. Finally, a novel description of the ventricular junction in large mammalian sheep hearts was proposed.

Ultra-high field (10.5T) diffusion-weighted MRI of the macaque brain

Diffu0sion-weighted magnetic resonance imaging (dMRI) is a non-invasive imaging technique that provides information about the barriers to the diffusion of water molecules in tissue. In the brain, this information can be used in several important ways, including to examine tissue abnormalities associated with brain disorders and to infer anatomical connectivity and the organization of white matter bundles through the use of tractography algorithms. However, dMRI also presents certain challenges. For example, historically, the biological validation of tractography models has shown only moderate correlations with anatomical connectivity as determined through invasive tract-tracing studies. Some of the factors contributing to such issues are low spatial resolution, low signal-to-noise ratios, and long scan times required for high-quality data, along with modeling challenges like complex fiber crossing patterns. Leveraging the capabilities provided by an ultra-high field scanner combined with denoising, we have acquired whole-brain, 0.58 mm isotropic resolution dMRI with a 2D-single shot echo planar imaging sequence on a 10.5 Tesla scanner in anesthetized macaques. These data produced high-quality tractograms and maps of scalar diffusion metrics in white matter. This work demonstrates the feasibility and motivation for in-vivo dMRI studies seeking to benefit from ultra-high fields.

In-vivo data-driven parcellation of Heschl's gyrus using structural connectivity

The human auditory cortex around Heschl’s gyrus (HG) exhibits diverging patterns across individuals owing to the heterogeneity of its substructures. In this study, we investigated the subregions of the human auditory cortex using data-driven machine-learning techniques at the individual level and assessed their structural and functional profiles. We studied an openly accessible large dataset of the Human Connectome Project and identified the subregions of the HG in humans using data-driven clustering techniques with individually calculated imaging features of cortical folding and structural connectivity information obtained via diffusion magnetic resonance imaging tractography. We characterized the structural and functional profiles of each HG subregion according to the cortical morphology, microstructure, and functional connectivity at rest. We found three subregions. The first subregion (HG₁) occupied the central portion of HG, the second subregion (HG₂) occupied the medial-posterior-superior part of HG, and the third subregion (HG₃) occupied the lateral-anterior-inferior part of HG. The HG₃ exhibited strong structural and functional connectivity to the association and paralimbic areas, and the HG₁ exhibited a higher myelin density and larger cortical thickness than other subregions. A functional gradient analysis revealed a gradual axis expanding from the HG₂ to the HG₃. Our findings clarify the individually varying structural and functional organization of human HG subregions and provide insights into the substructures of the human auditory cortex.