Advances in diffusion MRI acquisition and processing in the Human Connectome Project

Recovering high-quality fiber orientation distributions from a reduced number of diffusion-weighted images using a model-driven deep learning architecture

PURPOSE

The aim of this study was to develop a model-based deep learning architecture to accurately reconstruct fiber orientation distributions (FODs) from a reduced number of diffusion-weighted images (DWIs), facilitating accurate analysis with reduced acquisition times.

METHODS

Our proposed architecture, Spherical Deconvolution Network (SDNet), performed FOD reconstruction by mapping 30 DWIs to fully sampled FODs, which have been fit to 288 DWIs. SDNet included DWI-consistency blocks within the network architecture, and a fixel-classification penalty within the loss function. SDNet was trained on a subset of the Human Connectome Project, and its performance compared with FOD-Net, and multishell multitissue constrained spherical deconvolution.

RESULTS

SDNet achieved the strongest results with respect to angular correlation coefficient and sum of squared errors. When the impact of the fixel-classification penalty was increased, we observed an improvement in performance metrics reliant on segmenting the FODs into the correct number of fixels.

CONCLUSION

Inclusion of DWI-consistency blocks improved reconstruction performance, and the fixel-classification penalty term offered increased control over the angular separation of fixels in the reconstructed FODs.

Modeling the impact of MRI acquisition bias on structural connectomes: Harmonizing structural connectomes

One way to increase the statistical power and generalizability of neuroimaging studies is to collect data at multiple sites or merge multiple cohorts. However, this usually comes with site-related biases due to the heterogeneity of scanners and acquisition parameters, negatively impacting sensitivity. Brain structural connectomes are not an exception: Being derived from T1-weighted and diffusion-weighted magnetic resonance images, structural connectivity is impacted by differences in imaging protocol. Beyond minimizing acquisition parameter differences, removing bias with postprocessing is essential. In this work we create, from the exhaustive Human Connectome Project Young Adult dataset, a resampled dataset of different b-values and spatial resolutions, modeling a cohort scanned across multiple sites. After demonstrating the statistical impact of acquisition parameters on connectivity, we propose a linear regression with explicit modeling of b-value and spatial resolution, and validate its performance on separate datasets. We show that b-value and spatial resolution affect connectivity in different ways and that acquisition bias can be reduced using a linear regression informed by the acquisition parameters while retaining interindividual differences and hence boosting fingerprinting performance. We also demonstrate the generative potential of our model, and its generalization capability in an independent dataset reflective of typical acquisition practices in clinical settings.

Nonlinear Gradient Field Mapping Using a Spherical Grid Phantom for 3 and 7 Tesla MR Imaging Systems Equipped with High-performance Gradient Coils

Recent high-performance gradient coils are fabricated mainly at the expense of spatial linearity. In this study, we measured the spatial nonlinearity of the magnetic field generated by the gradient coils of two MRI systems with high-performance gradient coils. The nonlinearity of the gradient fields was measured using 3D gradient echo sequences and a spherical phantom with a built-in lattice structure. The spatial variation of the gradient field was approximated to the 3rd order polynomials. The coefficients of the polynomials were calculated using the steepest descent method. The geometric distortion of the acquired 3D MR images was corrected using the polynomials and compared with the 3D images corrected using the harmonic functions provided by the MRI venders. As a result, it was found that the nonlinearity correction formulae provided by the vendors were insufficient and needed to be verified or corrected using a geometric phantom such as used in this study.

An information-theoretic framework for conditional causality analysis of brain networks

Identifying directed network models for multivariate time series is a ubiquitous problem in data science. Granger causality measure (GCM) and conditional GCM (cGCM) are widely used methods for identifying directed connections between time series. Both GCM and cGCM have frequency-domain formulations to characterize the dependence of time series in the spectral domain. However, the original methods were developed using a heuristic approach without rigorous theoretical explanations. To overcome the limitation, the minimum-entropy (ME) estimation approach was introduced in our previous work (Ning & Rathi, 2018) to generalize GCM and cGCM with more rigorous frequency-domain formulations. In this work, this information-theoretic framework is further generalized with three formulations for conditional causality analysis using techniques in control theory, such as state-space representations and spectral factorizations. The three conditional causal measures are developed based on different ME estimation procedures that are motivated by equivalent formulations of the classical minimum mean squared error estimation method. The relationship between the three formulations of conditional causality measures is analyzed theoretically. Their performance is evaluated using simulations and real neuroimaging data to analyze brain networks. The results show that the proposed methods provide more accurate network structures than the original approach.

Evaluation of epilepsy lesion visualisation enhancement in low-field MRI using image quality transfer: a preliminary investigation of clinical potential for applications in developing countries

PURPOSE

Low-field (LF) MRI scanners are common in many Low- and middle-Income countries, but they provide images with worse spatial resolution and contrast than high-field (HF) scanners. Image Quality Transfer (IQT) is a machine learning framework to enhance images based on high-quality references that has recently adapted to LF MRI. In this study we aim to assess if it can improve lesion visualisation compared to LF MRI scans in children with epilepsy.

METHODS

T1-weighted, T2-weighted and FLAIR were acquired from 12 patients (5 to 18 years old, 7 males) with clinical diagnosis of intractable epilepsy on a 0.36T (LF) and a 1.5T scanner (HF). LF images were enhanced with IQT. Seven radiologists blindly evaluated the differentiation between normal grey matter (GM) and white matter (WM) and the extension and definition of epileptogenic lesions in LF, HF and IQT-enhanced images.

RESULTS

When images were evaluated independently, GM-WM differentiation scores of IQT outputs were 26% higher, 17% higher and 12% lower than LF for T1, T2 and FLAIR. Lesion definition scores were 8-34% lower than LF, but became 3% higher than LF for FLAIR and T1 when images were seen side by side. Radiologists with expertise at HF scored IQT images higher than those with expertise at LF.

CONCLUSION

IQT generally improved the image quality assessments. Evaluation of pathology on IQT-enhanced images was affected by familiarity with HF/IQT image appearance. These preliminary results show that IQT could have an important impact on neuroradiology practice where HF MRI is not available.

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.

Comparison of T1-weighted landmark placement and ROI transfer onto diffusion-weighted EPI sequences for targeted tractography tasks in the optic nerve

Diffusion-based tractography in the optic nerve requires sampling strategies assisted by anatomical landmark information (regions of interest [ROIs]). We aimed to investigate the feasibility of expert-placed, high-resolution T1-weighted ROI-data transfer onto lower spatial resolution diffusion-weighted images. Slab volumes from 20 volunteers were acquired and preprocessed including distortion bias correction and artifact reduction. Constrained spherical deconvolution was used to generate a directional diffusion information grid (fibre orientation distribution-model [FOD]). Three neuroradiologists marked landmarks on both diffusion imaging variants and structural datasets. Structural ROI information (volumetric interpolated breath-hold sequence [VIBE]) was respectively registered (linear with 6/12 degrees of freedom [DOF]) onto single-shot EPI (ss-EPI) and readout-segmented EPI (rs-EPI) volumes, respectively. All eight ROI/FOD-combinations were compared in a targeted tractography task of the optic nerve pathway. Inter-rater reliability for placed ROIs among experts was highest in VIBE images (lower confidence interval 0.84 to 0.97, mean 0.91) and lower in both ss-EPI (0.61 to 0.95, mean 0.79) and rs-EPI (0.59 to 0.86, mean 0.70). Tractography success rate based on streamline selection performance was highest in VIBE-drawn ROIs registered (6-DOF) onto rs-EPI FOD (70.0% over 5%-threshold, capped to failed ratio 39/16) followed by both 12-DOF-registered (67.5%; 41/16) and nonregistered VIBE (67.5%; 40/23). On ss-EPI FOD, VIBE-ROI-datasets obtained fewer streamlines overall with each at 55.0% above 5%-threshold and with lower capped to failed ratio (6-DOF: 35/36; 12-DOF: 34/34, nonregistered 33/36). The combination of VIBE-placed ROIs (highest inter-rater reliability) with 6-DOF registration onto rs-EPI targets (best streamline selection performance) is most suitable for white matter template generation required in group studies.

Voxel-Wise Fusion of 3T and 7T Diffusion MRI Data to Extract more Accurate Fiber Orientations

While 7T diffusion magnetic resonance imaging (dMRI) has high spatial resolution, its diffusion imaging quality is usually affected by signal loss due to B1 inhomogeneity, T2 decay, susceptibility, and chemical shift. In contrast, 3T dMRI has relative higher diffusion angular resolution, but lower spatial resolution. Combination of 3T and 7T dMRI, thus, may provide more detailed and accurate information about the voxel-wise fiber orientations to better understand the structural brain connectivity. However, this topic has not yet been thoroughly explored until now. In this study, we explored the feasibility of fusing 3T and 7T dMRI data to extract voxel-wise quantitative parameters at higher spatial resolution. After 3T and 7T dMRI data was preprocessed, respectively, 3T dMRI volumes were coregistered into 7T dMRI space. Then, 7T dMRI data was harmonized to the coregistered 3T dMRI B0 (b = 0) images. Last, harmonized 7T dMRI data was fused with 3T dMRI data according to four fusion rules proposed in this study. We employed high-quality 3T and 7T dMRI datasets (N = 24) from the Human Connectome Project to test our algorithms. The diffusion tensors (DTs) and orientation distribution functions (ODFs) estimated from the 3T-7T fused dMRI volumes were statistically analyzed. More voxels containing multiple fiber populations were found from the fused dMRI data than from 7T dMRI data set. Moreover, extra fiber directions were extracted in temporal brain regions from the fused dMRI data at Otsu's thresholds of quantitative anisotropy, but could not be extracted from 7T dMRI dataset. This study provides novel algorithms to fuse intra-subject 3T and 7T dMRI data for extracting more detailed information of voxel-wise quantitative parameters, and a new perspective to build more accurate structural brain networks.

Language and the Cerebellum: Structural Connectivity to the Eloquent Brain

Neurobiological models of receptive language have focused on the left-hemisphere perisylvian cortex with the assumption that the cerebellum supports peri-linguistic cognitive processes such as verbal working memory. The goal of this study was to identify language-sensitive regions of the cerebellum then map the structural connectivity profile of these regions. Functional imaging data and diffusion-weighted imaging data from the Human Connectome Project (HCP) were analyzed. We found that (a) working memory, motor activity, and language comprehension activated partially overlapping but mostly unique subregions of the cerebellum; (b) the linguistic portion of the cerebello-thalamo-cortical circuit was more extensive than the linguistic portion of the cortico-ponto-cerebellar tract; (c) there was a frontal-lobe bias in the connectivity from the cerebellum to the cerebrum; (d) there was some degree of specificity; and (e) for some cerebellar tracts, individual differences in picture identification ability covaried with fractional anisotropy metrics. These findings yield insights into the structural connectivity of the cerebellum as relates to the uniquely human process of language comprehension.

Discovering the effective connectome of the brain with dynamic Bayesian DAG learning

Understanding the complex mechanisms of the brain can be unraveled by extracting the Dynamic Effective Connectome (DEC). Recently, score-based Directed Acyclic Graph (DAG) discovery methods have shown significant improvements in extracting the causal structure and inferring effective connectivity. However, learning DEC through these methods still faces two main challenges: one with the fundamental impotence of high-dimensional dynamic DAG discovery methods and the other with the low quality of fMRI data. In this paper, we introduce Bayesian Dynamic DAG learning with M-matrices Acyclicity characterization (BDyMA) method to address the challenges in discovering DEC. The presented dynamic DAG enables us to discover direct feedback loop edges as well. Leveraging an unconstrained framework in the BDyMA method leads to more accurate results in detecting high-dimensional networks, achieving sparser outcomes, making it particularly suitable for extracting DEC. Additionally, the score function of the BDyMA method allows the incorporation of prior knowledge into the process of dynamic causal discovery which further enhances the accuracy of results. Comprehensive simulations on synthetic data and experiments on Human Connectome Project (HCP) data demonstrate that our method can handle both of the two main challenges, yielding more accurate and reliable DEC compared to state-of-the-art and traditional methods. Additionally, we investigate the trustworthiness of DTI data as prior knowledge for DEC discovery and show the improvements in DEC discovery when the DTI data is incorporated into the process.

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.

Effect of simultaneous multislice imaging, slice properties, and repetition time on the measured magnetic resonance biexponential intravoxel incoherent motion in the liver

OBJECTIVES

This study aims to investigate the previously reported dependency of intravoxel incoherent motion (IVIM) parameters on simultaneous multislice (SMS) acquisition and repetition time (TR). This includes the influence of slice thickness, slice gaps, and slice order on measured IVIM parameters.

MATERIALS AND METHODS

Diffusion-weighted imaging (DWI) of the liver was performed on 10 healthy volunteers (aged 20-30 years) at 3T with a slice thickness of 5 mm, a slice gap of 5 mm, and a linear slice order. Diffusion-weighted images were acquired with 19 b-values (0-800 s/mm2) using both conventional slice excitation with an acceleration factor of one (AF1) and SMS excitation with an acceleration factor of three (AF3). Each of these measurements were carried out with two repetition times (TRs)- 1,300 ms (prefix s) and 4,500 ms (prefix l)-resulting in four different combinations: sAF1, sAF3, lAF1, and lAF3. Five volunteers underwent additional measurements using a 10 mm slice thickness and with AF1. Median signal values in the liver were used to determine the biexponential IVIM parameters. Statistical significances were assessed using the Kruskal-Wallis test, Wilcoxon signed-rank test, and Student's t-test. In-silico investigations were also used to interpret the data.

RESULTS

There were no significant differences between the biexponential IVIM parameters acquired from sAF1, sAF3, lAF1, and lAF3. Median values of the perfusion fraction f were as follows: 29.9% (sAF1), 26.9% (sAF3), 28.1% (lAF1), and 27.5% (lAF3). In the 10 mm-thick slices, f decreased from 31.3% (lAF1) to 27.4% (sAF1) (p = 0.141).

CONCLUSION

The slice excitation mode did not appear to have any significant influence on the biexponential IVIM parameters. However, our simulations, as well as values reported from previous publications, show that slice thickness, slice gaps, and slice order are relevant and should thus be reported in IVIM studies.

Characterizing white matter connectome abnormalities in patients with temporal lobe epilepsy using threshold-free network-based statistics

INTRODUCTION

Emerging evidence illustrates that temporal lobe epilepsy (TLE) involves network disruptions represented by hyperexcitability and other seizure-related neural plasticity. However, these associations are not well-characterized. Our study characterizes the whole brain white matter connectome abnormalities in TLE patients compared to healthy controls (HCs) from the prospective Epilepsy Connectome Project study. Furthermore, we assessed whether aberrant white matter connections are differentially related to cognitive impairment and a history of focal-to-bilateral tonic-clonic (FBTC) seizures.

METHODS

Multi-shell connectome MRI data were preprocessed using the DESIGNER guidelines. The IIT Destrieux gray matter atlas was used to derive the 162 × 162 structural connectivity matrices (SCMs) using MRTrix3. ComBat data harmonization was applied to harmonize the SCMs from pre- and post-scanner upgrade acquisitions. Threshold-free network-based statistics were used for statistical analysis of the harmonized SCMs. Cognitive impairment status and FBTC seizure status were then correlated with these findings.

RESULTS

We employed connectome measurements from 142 subjects, including 92 patients with TLE (36 males, mean age = 40.1 ± 11.7 years) and 50 HCs (25 males, mean age = 32.6 ± 10.2 years). Our analysis revealed overall significant decreases in cross-sectional area (CSA) of the white matter tract in TLE group compared to controls, indicating decreased white matter tract integrity and connectivity abnormalities in addition to apparent differences in graph theoretic measures of connectivity and network-based statistics. Focal and generalized cognitive impaired TLE patients showcased higher trend-level abnormalities in the white matter connectome via decreased CSA than those with no cognitive impairment. Patients with a positive FBTC seizure history also showed trend-level findings of association via decreased CSA.

CONCLUSIONS

Widespread global aberrant white matter connectome changes were observed in TLE patients and characterized by seizure history and cognitive impairment, laying a foundation for future studies to expand on and validate the novel biomarkers and further elucidate TLE's impact on brain plasticity.

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

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

Diffusion Basis Spectrum Imaging of White Matter in Schizophrenia and Bipolar Disorder

Background

Multiple studies point to the role of neuroinflammation in the pathophysiology of schizophrenia (SCZ), however, there have been few in vivo tools for imaging brain inflammation. Diffusion basis spectrum imaging (DBSI) is an advanced diffusion-based MRI method developed to quantitatively assess microstructural alternations relating to neuroinflammation, axonal fiber, and other white matter (WM) pathologies.

Methods

We acquired one-hour-long high-directional diffusion MRI data from young control (CON, n = 27), schizophrenia (SCZ, n = 21), and bipolar disorder (BPD, n = 21) participants aged 18-30. We applied Tract-based Spatial Statistics (TBSS) to allow whole-brain WM analyses and compare DBSI-derived isotropic and anisotropic diffusion measures between groups. Clinical relationships of DBSI metrics with clinical symptoms were assessed across SCZ and control participants.

Results

In SCZ participants, we found a generalized increase in DBSI-derived cellularity (a putative marker of neuroinflammation), a decrease in restricted fiber fraction (a putative marker of apparent axonal density), and an increase in extra-axonal water (a putative marker of vasogenic edema) across several WM tracts. There were only minimal WM abnormalities noted in BPD, mainly in regions of the corpus callosum (increase in DTI-derived RD and extra-axonal water). DBSI metrics showed significant partial correlations with psychosis and mood symptoms across groups.

Conclusion

Our findings suggest that SCZ involves generalized white matter neuroinflammation, decreased fiber density, and demyelination, which is not seen in bipolar disorder. Larger studies are needed to identify medication-related effects. DBSI metrics could help identify high-risk groups requiring early interventions to prevent the onset of psychosis and improve outcomes.

Mapping short association fibre connectivity up to V3 in the human brain in vivo

Short association fibres (SAF) are the most abundant fibre pathways in the human white matter. Until recently, SAF could not be mapped comprehensively in vivo because diffusion weighted magnetic resonance imaging with sufficiently high spatial resolution needed to map these thin and short pathways was not possible. Recent developments in acquisition hardware and sequences allowed us to create a dedicated in vivo method for mapping the SAF based on sub-millimetre spatial resolution diffusion weighted tractography, which we validated in the human primary (V1) and secondary (V2) visual cortex against the expected SAF retinotopic order. Here, we extended our original study to assess the feasibility of the method to map SAF in higher cortical areas by including SAF up to V3. Our results reproduced the expected retinotopic order of SAF in the V2-V3 and V1-V3 stream, demonstrating greater robustness to the shorter V1-V2 and V2-V3 than the longer V1-V3 connections. The demonstrated ability of the method to map higher-order SAF connectivity patterns in vivo is an important step towards its application across the brain.

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.

Densely sampled stimulus-response map of human cortex with single pulse TMS-EEG and its relation to whole brain neuroimaging measures

Large-scale networks underpin brain functions. How such networks respond to focal stimulation can help decipher complex brain processes and optimize brain stimulation treatments. To map such stimulation-response patterns across the brain non-invasively, we recorded concurrent EEG responses from single-pulse transcranial magnetic stimulation (i.e., TMS-EEG) from over 100 cortical regions with two orthogonal coil orientations from one densely-sampled individual. We also acquired Human Connectome Project (HCP)-styled diffusion imaging scans (six), resting-state functional Magnetic Resonance Imaging (fMRI) scans (120 mins), resting-state EEG scans (108 mins), and structural MR scans (T1- and T2-weighted). Using the TMS-EEG data, we applied network science-based community detection to reveal insights about the brain's causal-functional organization from both a stimulation and recording perspective. We also computed structural and functional maps and the electric field of each TMS stimulation condition. Altogether, we hope the release of this densely sampled (n=1) dataset will be a uniquely valuable resource for both basic and clinical neuroscience research.

Tractometry of the Human Connectome Project: resources and insights

Introduction

The Human Connectome Project (HCP) has become a keystone dataset in human neuroscience, with a plethora of important applications in advancing brain imaging methods and an understanding of the human brain. We focused on tractometry of HCP diffusion-weighted MRI (dMRI) data.

Methods

We used an open-source software library (pyAFQ; https://yeatmanlab.github.io/pyAFQ) to perform probabilistic tractography and delineate the major white matter pathways in the HCP subjects that have a complete dMRI acquisition (n = 1,041). We used diffusion kurtosis imaging (DKI) to model white matter microstructure in each voxel of the white matter, and extracted tract profiles of DKI-derived tissue properties along the length of the tracts. We explored the empirical properties of the data: first, we assessed the heritability of DKI tissue properties using the known genetic linkage of the large number of twin pairs sampled in HCP. Second, we tested the ability of tractometry to serve as the basis for predictive models of individual characteristics (e.g., age, crystallized/fluid intelligence, reading ability, etc.), compared to local connectome features. To facilitate the exploration of the dataset we created a new web-based visualization tool and use this tool to visualize the data in the HCP tractometry dataset. Finally, we used the HCP dataset as a test-bed for a new technological innovation: the TRX file-format for representation of dMRI-based streamlines.

Results

We released the processing outputs and tract profiles as a publicly available data resource through the AWS Open Data program's Open Neurodata repository. We found heritability as high as 0.9 for DKI-based metrics in some brain pathways. We also found that tractometry extracts as much useful information about individual differences as the local connectome method. We released a new web-based visualization tool for tractometry-"Tractoscope" (https://nrdg.github.io/tractoscope). We found that the TRX files require considerably less disk space-a crucial attribute for large datasets like HCP. In addition, TRX incorporates a specification for grouping streamlines, further simplifying tractometry analysis.

Multimodal study of multilevel pulvino-temporal connections: a new piece in the puzzle of lexical retrieval networks

Advanced methods of imaging and mapping the healthy and lesioned brain have allowed for the identification of the cortical nodes and white matter tracts supporting the dual neurofunctional organization of language networks in a dorsal phonological and a ventral semantic stream. Much less understood are the anatomical correlates of the interaction between the two streams; one hypothesis being that of a subcortically mediated interaction, through crossed cortico-striato-thalamo-cortical and cortico-thalamo-cortical loops. In this regard, the pulvinar is the thalamic subdivision that has most regularly appeared as implicated in the processing of lexical retrieval. However, descriptions of its connections with temporal (language) areas remain scarce. Here we assess this pulvino-temporal connectivity using a combination of state-of-the-art techniques: white matter stimulation in awake surgery and postoperative diffusion MRI (n = 4), virtual dissection from the Human Connectome Project 3 and 7 T datasets (n = 172) and operative microscope-assisted post-mortem fibre dissection (n = 12). We demonstrate the presence of four fundamental fibre contingents: (i) the anterior component (Arnold's bundle proper) initially described by Arnold in the 19th century and destined to the anterior temporal lobe; (ii) the optic radiations-like component, which leaves the pulvinar accompanying the optical radiations and reaches the posterior basal temporal cortices; (iii) the lateral component, which crosses the temporal stem orthogonally and reaches the middle temporal gyrus; and (iv) the auditory radiations-like component, which leaves the pulvinar accompanying the auditory radiations to the superomedial aspect of the temporal operculum, just posteriorly to Heschl's gyrus. Each of those components might correspond to a different level of information processing involved in the lexical retrieval process of picture naming.

DIMOND: DIffusion Model OptimizatioN with Deep Learning

Diffusion magnetic resonance imaging is an important tool for mapping tissue microstructure and structural connectivity non-invasively in the in vivo human brain. Numerous diffusion signal models are proposed to quantify microstructural properties. Nonetheless, accurate estimation of model parameters is computationally expensive and impeded by image noise. Supervised deep learning-based estimation approaches exhibit efficiency and superior performance but require additional training data and may be not generalizable. A new DIffusion Model OptimizatioN framework using physics-informed and self-supervised Deep learning entitled "DIMOND" is proposed to address this problem. DIMOND employs a neural network to map input image data to model parameters and optimizes the network by minimizing the difference between the input acquired data and synthetic data generated via the diffusion model parametrized by network outputs. DIMOND produces accurate diffusion tensor imaging results and is generalizable across subjects and datasets. Moreover, DIMOND outperforms conventional methods for fitting sophisticated microstructural models including the kurtosis and NODDI model. Importantly, DIMOND reduces NODDI model fitting time from hours to minutes, or seconds by leveraging transfer learning. In summary, the self-supervised manner, high efficacy, and efficiency of DIMOND increase the practical feasibility and adoption of microstructure and connectivity mapping in clinical and neuroscientific applications.

MR Imaging in Ataxias: Consensus Recommendations by the Ataxia Global Initiative Working Group on MRI Biomarkers

With many viable strategies in the therapeutic pipeline, upcoming clinical trials in hereditary and sporadic degenerative ataxias will benefit from non-invasive MRI biomarkers for patient stratification and the evaluation of therapies. The MRI Biomarkers Working Group of the Ataxia Global Initiative therefore devised guidelines to facilitate harmonized MRI data acquisition in clinical research and trials in ataxias. Recommendations are provided for a basic structural MRI protocol that can be used for clinical care and for an advanced multi-modal MRI protocol relevant for research and trial settings. The advanced protocol consists of modalities with demonstrated utility for tracking brain changes in degenerative ataxias and includes structural MRI, magnetic resonance spectroscopy, diffusion MRI, quantitative susceptibility mapping, and resting-state functional MRI. Acceptable ranges of acquisition parameters are provided to accommodate diverse scanner hardware in research and clinical contexts while maintaining a minimum standard of data quality. Important technical considerations in setting up an advanced multi-modal protocol are outlined, including the order of pulse sequences, and example software packages commonly used for data analysis are provided. Outcome measures most relevant for ataxias are highlighted with use cases from recent ataxia literature. Finally, to facilitate access to the recommendations by the ataxia clinical and research community, examples of datasets collected with the recommended parameters are provided and platform-specific protocols are shared via the Open Science Framework.

FlexDTI: flexible diffusion gradient encoding scheme-based highly efficient diffusion tensor imaging using deep learning

Objective. Most deep neural network-based diffusion tensor imaging methods require the diffusion gradients' number and directions in the data to be reconstructed to match those in the training data. This work aims to develop and evaluate a novel dynamic-convolution-based method called FlexDTI for highly efficient diffusion tensor reconstruction with flexible diffusion encoding gradient scheme.Approach. FlexDTI was developed to achieve high-quality DTI parametric mapping with flexible number and directions of diffusion encoding gradients. The method used dynamic convolution kernels to embed diffusion gradient direction information into feature maps of the corresponding diffusion signal. Furthermore, it realized the generalization of a flexible number of diffusion gradient directions by setting the maximum number of input channels of the network. The network was trained and tested using datasets from the Human Connectome Project and local hospitals. Results from FlexDTI and other advanced tensor parameter estimation methods were compared.Main results. Compared to other methods, FlexDTI successfully achieves high-quality diffusion tensor-derived parameters even if the number and directions of diffusion encoding gradients change. It reduces normalized root mean squared error by about 50% on fractional anisotropy and 15% on mean diffusivity, compared with the state-of-the-art deep learning method with flexible diffusion encoding gradient scheme.Significance. FlexDTI can well learn diffusion gradient direction information to achieve generalized DTI reconstruction with flexible diffusion gradient scheme. Both flexibility and reconstruction quality can be taken into account in this network.

An in vivo Dissection, and Analysis of Socio-Affective Symptoms related to Cerebellum-Midbrain Reward Circuitry in Humans

Emerging research in non-human animals implicates cerebellar projections to the ventral tegmental area (VTA) in appetitive behaviors, but these circuits have not been characterized in humans. Here, we mapped cerebello-VTA white-matter connectivity in humans using probabilistic tractography on diffusion imaging data from the Human Connectome Project. We uncovered the topographical organization of these connections by separately tracking from parcels of cerebellar lobule VI, crus I/II, vermis, paravermis, and cerebrocerebellum. Results revealed that connections from the cerebellum to the VTA predominantly originate in the right hemisphere, interposed nucleus, and paravermal cortex, and terminate mostly ipsilaterally. Paravermal crus I sends the most connections to the VTA compared to other lobules. We discovered a medial-to-lateral gradient of connectivity, such that the medial cerebellum has the highest connectivity with the VTA. Individual differences in microstructure were associated with measures of negative affect and social functioning. By splitting the tracts into quarters, we found that the socio-affective effects were driven by the third quarter of the tract, corresponding to the point at which the fibers leave the deep nuclei. Taken together, we produced detailed maps of cerebello-VTA structural connectivity for the first time in humans and established their relevance for trait differences in socio-affective regulation.

Diffusion-tensor-imaging 1-year-old and 2-year-old infant brain atlases with comprehensive gray and white matter labels

Human infancy is marked by fastest postnatal brain structural changes. It also coincides with the onset of many neurodevelopmental disorders. Atlas-based automated structure labeling has been widely used for analyzing various neuroimaging data. However, the relatively large and nonlinear neuroanatomical differences between infant and adult brains can lead to significant offsets of the labeled structures in infant brains when adult brain atlas is used. Age-specific 1- and 2-year-old brain atlases covering all major gray and white matter (GM and WM) structures with diffusion tensor imaging (DTI) and structural MRI are critical for precision medicine for infant population yet have not been established. In this study, high-quality DTI and structural MRI data were obtained from 50 healthy children to build up three-dimensional age-specific 1- and 2-year-old brain templates and atlases. Age-specific templates include a single-subject template as well as two population-averaged templates from linear and nonlinear transformation, respectively. Each age-specific atlas consists of 124 comprehensively labeled major GM and WM structures, including 52 cerebral cortical, 10 deep GM, 40 WM, and 22 brainstem and cerebellar structures. When combined with appropriate registration methods, the established atlases can be used for highly accurate automatic labeling of any given infant brain MRI. We demonstrated that one can automatically and effectively delineate deep WM microstructural development from 3 to 38 months by using these age-specific atlases. These established 1- and 2-year-old infant brain DTI atlases can advance our understanding of typical brain development and serve as clinical anatomical references for brain disorders during infancy.

Association of Blood Pressure With Brain White Matter Microstructural Integrity Assessed With MRI Diffusion Tensor Imaging in Healthy Young Adults

BACKGROUND

High blood pressure (BP) in middle-aged and older adults is associated with a brain white matter (WM) microstructural abnormality. However, little evidence is available in healthy young adults. We investigated the associations between high BP and WM microstructural integrity in young adults.

METHODS

This study included 1015 healthy young adults (542 women, 22-37 years) from the Human Connectome Project. Brachial systolic and diastolic BP were measured using a semiautomatic or manual sphygmomanometer. Diffusion-weighted magnetic resonance imaging was acquired to obtain diffusion tensor imaging metrics of free water (FW) content, FW-corrected WM fractional anisotropy, axial diffusivity, radial diffusivity, and mean diffusivity. Using whole-brain voxel-wise linear regression models and ANCOVA, we examined associations of BP and hypertension stage with diffusion tensor imaging metrics after adjusting for age, sex, education, body mass index, smoking status, alcohol consumption history, and differences in the b value used for diffusion magnetic resonance imaging.

RESULTS

Systolic and diastolic BP of the sample (mean±SD) were 122.8±13.0 and 76.0±9.9 mm Hg, respectively. Associations of BP with diffusion tensor imaging metrics revealed regional heterogeneity for FW-corrected fractional anisotropy. High BP and high hypertension stage were associated with higher FW and lower FW-corrected axial diffusivity, FW-corrected radial diffusivity, and FW-corrected mean diffusivity. Moreover, associations of high diastolic BP and hypertension stage with high FW were found only in men not in women.

CONCLUSIONS

High BP in young adults is associated with altered brain WM microstructural integrity, suggesting that high BP may have damaging effects on brain WM microstructural integrity in early adulthood, particularly in men.

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Structural Connectivity of the Human Piriform Cortex: an Exploratory Study

BACKGROUND AND OBJECTIVES

The piriform cortex (PC) is part of the primary olfactory network in humans. Recent findings suggest that it plays a role in pathophysiology of epilepsy. Therefore, studying its connectivity can further our understanding of seizure propagation in epilepsy. We aimed to explore the structural connectivity of PC using high-quality human connectome project data coupled with segmentation of PC on anatomic MRI.

METHODS

Twenty subjects were randomly selected from the human connectome project database, and PC was traced on each hemisphere. Probabilistic whole-brain tractography was then used to visualize PC connectivity.

RESULTS

The strongest connectivity was noted between PC and ipsilateral insula in both hemispheres. Specifically, the posterior long gyrus of each insula was predominantly connected to PC. This was followed by connections between PC and basal ganglia as well as orbital frontal cortices.

CONCLUSION

The PC has the strongest connectivity with the insula bilaterally. Specifically, the posterior long gyri of insula have the strongest connectivity. This finding may provide additional insight for localizing and treating temporo-insular epilepsy.

Generalising XTRACT tractography protocols across common macaque brain templates

Non-human primates are extensively used in neuroscience research as models of the human brain, with the rhesus macaque being a prominent example. We have previously introduced a set of tractography protocols (XTRACT) for reconstructing 42 corresponding white matter (WM) bundles in the human and the macaque brain and have shown cross-species comparisons using such bundles as WM landmarks. Our original XTRACT protocols were developed using the F99 macaque brain template. However, additional macaque template brains are becoming increasingly common. Here, we generalise the XTRACT tractography protocol definitions across five macaque brain templates, including the F99, D99, INIA, Yerkes and NMT. We demonstrate equivalence of such protocols in two ways: (a) Firstly by comparing the bodies of the tracts derived using protocols defined across the different templates considered, (b) Secondly by comparing the projection patterns of the reconstructed tracts across the different templates in two cross-species (human-macaque) comparison tasks. The results confirm similarity of all predictions regardless of the macaque brain template used, providing direct evidence for the generalisability of these tractography protocols across the five considered templates.

On the stability of canonical correlation analysis and partial least squares with application to brain-behavior associations

Associations between datasets can be discovered through multivariate methods like Canonical Correlation Analysis (CCA) or Partial Least Squares (PLS). A requisite property for interpretability and generalizability of CCA/PLS associations is stability of their feature patterns. However, stability of CCA/PLS in high-dimensional datasets is questionable, as found in empirical characterizations. To study these issues systematically, we developed a generative modeling framework to simulate synthetic datasets. We found that when sample size is relatively small, but comparable to typical studies, CCA/PLS associations are highly unstable and inaccurate; both in their magnitude and importantly in the feature pattern underlying the association. We confirmed these trends across two neuroimaging modalities and in independent datasets with n ≈ 1000 and n = 20,000, and found that only the latter comprised sufficient observations for stable mappings between imaging-derived and behavioral features. We further developed a power calculator to provide sample sizes required for stability and reliability of multivariate analyses. Collectively, we characterize how to limit detrimental effects of overfitting on CCA/PLS stability, and provide recommendations for future studies.

DiMANI: diffusion MRI for anatomical nuclei imaging-Application for the direct visualization of thalamic subnuclei

The thalamus is a centrally located and heterogeneous brain structure that plays a critical role in various sensory, motor, and cognitive processes. However, visualizing the individual subnuclei of the thalamus using conventional MRI techniques is challenging. This difficulty has posed obstacles in targeting specific subnuclei for clinical interventions such as deep brain stimulation (DBS). In this paper, we present DiMANI, a novel method for directly visualizing the thalamic subnuclei using diffusion MRI (dMRI). The DiMANI contrast is computed by averaging, voxelwise, diffusion-weighted volumes enabling the direct distinction of thalamic subnuclei in individuals. We evaluated the reproducibility of DiMANI through multiple approaches. First, we utilized a unique dataset comprising 8 scans of a single participant collected over a 3-year period. Secondly, we quantitatively assessed manual segmentations of thalamic subnuclei for both intra-rater and inter-rater reliability. Thirdly, we qualitatively correlated DiMANI imaging data from several patients with Essential Tremor with the localization of implanted DBS electrodes and clinical observations. Lastly, we demonstrated that DiMANI can provide similar features at 3T and 7T MRI, using varying numbers of diffusion directions. Our results establish that DiMANI is a reproducible and clinically relevant method to directly visualize thalamic subnuclei. This has significant implications for the development of new DBS targets and the optimization of DBS therapy.

A unified filtering method for estimating asymmetric orientation distribution functions

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

Identification and validation of supervariants reveal novel loci associated with human white matter microstructure

As an essential part of the central nervous system, white matter coordinates communications between different brain regions and is related to a wide range of neurodegenerative and neuropsychiatric disorders. Previous genome-wide association studies (GWASs) have uncovered loci associated with white matter microstructure. However, GWASs suffer from limited reproducibility and difficulties in detecting multi-single-nucleotide polymorphism (multi-SNP) and epistatic effects. In this study, we adopt the concept of supervariants, a combination of alleles in multiple loci, to account for potential multi-SNP effects. We perform supervariant identification and validation to identify loci associated with 22 white matter fractional anisotropy phenotypes derived from diffusion tensor imaging. To increase reproducibility, we use United Kingdom (UK) Biobank White British (n = 30,842) data for discovery and internal validation, and UK Biobank White but non-British (n = 1927) data, Europeans from the Adolescent Brain Cognitive Development study (n = 4399) data, and Europeans from the Human Connectome Project (n = 319) data for external validation. We identify 23 novel loci on the discovery set that have not been reported in the previous GWASs on white matter microstructure. Among them, three supervariants on genomic regions 5q35.1, 8p21.2, and 19q13.32 have P-values lower than 0.05 in the meta-analysis of the three independent validation data sets. These supervariants contain genetic variants located in genes that have been related to brain structures, cognitive functions, and neuropsychiatric diseases. Our findings provide a better understanding of the genetic architecture underlying white matter microstructure.

Prenatal exposure to maternal disadvantage-related inflammatory biomarkers: associations with neonatal white matter microstructure

Prenatal exposure to heightened maternal inflammation has been associated with adverse neurodevelopmental outcomes, including atypical brain maturation and psychiatric illness. In mothers experiencing socioeconomic disadvantage, immune activation can be a product of the chronic stress inherent to such environmental hardship. While growing preclinical and clinical evidence has shown links between altered neonatal brain development and increased inflammatory states in utero, the potential mechanism by which socioeconomic disadvantage differentially impacts neural-immune crosstalk remains unclear. In the current study, we investigated associations between socioeconomic disadvantage, gestational inflammation, and neonatal white matter microstructure in 320 mother-infant dyads over-sampled for poverty. We analyzed maternal serum levels of four cytokines (IL-6, IL-8, IL-10, TNF-α) over the course of pregnancy in relation to offspring white matter microstructure and socioeconomic disadvantage. Higher average maternal IL-6 was associated with very low socioeconomic status (SES; INR < 200% poverty line) and lower neonatal corticospinal fractional anisotropy (FA) and lower uncinate axial diffusivity (AD). No other cytokine was associated with SES. Higher average maternal IL-10 was associated with lower FA and higher radial diffusivity (RD) in corpus callosum and corticospinal tracts, higher optic radiation RD, lower uncinate AD, and lower FA in inferior fronto-occipital fasciculus and anterior limb of internal capsule tracts. SES moderated the relationship between average maternal TNF-α levels during gestation and neonatal white matter diffusivity. When these interactions were decomposed, the patterns indicated that this association was significant and positive among very low SES neonates, whereby TNF-α was inversely and significantly associated with inferior cingulum AD. By contrast, among the more advantaged neonates (lower-to-higher SES [INR ≥ 200% poverty line]), TNF-α was positively and significantly associated with superior cingulum AD. Taken together, these findings suggest that the relationship between prenatal cytokine exposure and white matter microstructure differs as a function of SES. These patterns are consistent with a scenario where gestational inflammation's effects on white matter development diverge depending on the availability of foundational resources in utero.

Impact of brain parcellation on prediction performance in models of cognition and demographics

Brain connectivity analysis begins with the selection of a parcellation scheme that will define brain regions as nodes of a network whose connections will be studied. Brain connectivity has already been used in predictive modelling of cognition, but it remains unclear if the resolution of the parcellation used can systematically impact the predictive model performance. In this work, structural, functional and combined connectivity were each defined with five different parcellation schemes. The resolution and modality of the parcellation schemes were varied. Each connectivity defined with each parcellation was used to predict individual differences in age, education, sex, executive function, self-regulation, language, encoding and sequence processing. It was found that low-resolution functional parcellation consistently performed above chance at producing generalisable models of both demographics and cognition. However, no single parcellation scheme showed a superior predictive performance across all cognitive domains and demographics. In addition, although parcellation schemes impacted the graph theory measures of each connectivity type (structural, functional and combined), these differences did not account for the out-of-sample predictive performance of the models. Taken together, these findings demonstrate that while high-resolution parcellations may be beneficial for modelling specific individual differences, partial voluming of signals produced by the higher resolution of the parcellation likely disrupts model generalisability.

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.

Tract-based analyses of white matter in schizophrenia, bipolar disorder, aging, and dementia using high spatial and directional resolution diffusion imaging: a pilot study

Introduction

Structural brain connectivity abnormalities have been associated with several psychiatric disorders. Schizophrenia (SCZ) is a chronic disabling disorder associated with accelerated aging and increased risk of dementia, though brain findings in the disorder have rarely been directly compared to those that occur with aging.

Methods

We used an automated approach to reconstruct key white matter tracts and assessed tract integrity in five participant groups. We acquired one-hour-long high-directional diffusion MRI data from young control (CON, n =28), bipolar disorder (BPD, n =21), and SCZ (n =22) participants aged 18-30, and healthy elderly (ELD, n =15) and dementia (DEM, n =9) participants. Volume, fractional (FA), radial diffusivity (RD) and axial diffusivity (AD) of seven key white matter tracts (anterior thalamic radiation, ATR; dorsal and ventral cingulum bundle, CBD and CBV; corticospinal tract, CST; and the three superior longitudinal fasciculi: SLF-1, SLF-2 and SLF-3) were analyzed with TRACULA. Group comparisons in tract metrics were performed using multivariate and univariate analyses. Clinical relationships of tract metrics with recent and chronic symptoms were assessed in SCZ and BPD participants.

Results

A MANOVA showed group differences in FA (λ=0.5; p=0.0002) and RD (λ=0.35; p<0.0001) across the seven tracts, but no significant differences in tract AD and volume. Post-hoc analyses indicated lower tract FA and higher RD in ELD and DEM groups compared to CON, BPD and SCZ groups. Lower FA and higher RD in SCZ compared to CON did not meet statistical significance. In SCZ participants, a significant negative correlation was found between chronic psychosis severity and FA in the SLF-1 (r= -0.45; p=0.035), SLF-2 (r= -0.49; p=0.02) and SLF-3 (r= -0.44; p=0.042).

Discussion

Our results indicate impaired white matter tract integrity in elderly populations consistent with myelin damage. Impaired tract integrity in SCZ is most prominent in patients with advanced illness.

Ultrafast diffusion tensor imaging based on deep learning and multi-slice information sharing

Objective. Diffusion tensor imaging (DTI) is excellent for non-invasively quantifying tissue microstructure. Theoretically DTI can be achieved with six different diffusion weighted images and one reference image, but the tensor estimation accuracy is poor in this case. Increasing the number of diffusion directions has benefits for the tensor estimation accuracy, which results in long scan time and makes DTI sensitive to motion. It would be beneficial to decrease the scan time of DTI by using fewer diffusion-weighted images without compromising reconstruction quality.Approach. A novel DTI scan scheme was proposed to achieve fast DTI, where only three diffusion directions per slice was required under a specific direction switching manner, and a deep-learning based reconstruction method was utilized using multi-slice information sharing and correspondingT1-weighted image for high-quality DTI reconstruction. A network with two encoders developed from U-Net was implemented for better utilizing the diffusion data redundancy between neighboring slices. The method performed direct nonlinear mapping from diffusion-weighted images to diffusion tensor.Main results. The performance of the proposed method was verified on the Human Connectome Project public data and clinical patient data. High-quality mean diffusivity, fractional anisotropy, and directionally encoded colormap can be achieved with only three diffusion directions per slice.Significance. High-quality DTI-derived maps can be achieved in less than one minute of scan time. The great reduction of scan time will help push the wider application of DTI in clinical practice.

Structural connectome architecture shapes the maturation of cortical morphology from childhood to adolescence

Cortical thinning is an important hallmark of the maturation of brain morphology during childhood and adolescence. However, the connectome-based wiring mechanism that underlies cortical maturation remains unclear. Here, we show cortical thinning patterns primarily located in the lateral frontal and parietal heteromodal nodes during childhood and adolescence, which are structurally constrained by white matter network architecture and are particularly represented using a network-based diffusion model. Furthermore, connectome-based constraints are regionally heterogeneous, with the largest constraints residing in frontoparietal nodes, and are associated with gene expression signatures of microstructural neurodevelopmental events. These results are highly reproducible in another independent dataset. These findings advance our understanding of network-level mechanisms and the associated genetic basis that underlies the maturational process of cortical morphology during childhood and adolescence.

Cortico-basal ganglia white matter microstructure is linked to restricted repetitive behavior in autism spectrum disorder

BACKGROUND

Restricted repetitive behavior (RRB) is one of two behavioral domains required for the diagnosis of autism spectrum disorder (ASD). Neuroimaging is widely used to study brain alterations associated with ASD and the domain of social and communication deficits, but there has been less work regarding brain alterations linked to RRB.

METHODS

We utilized neuroimaging data from the National Institute of Mental Health Data Archive to assess basal ganglia and cerebellum structure in a cohort of children and adolescents with ASD compared to typically developing (TD) controls. We evaluated regional gray matter volumes from T1-weighted anatomical scans and assessed diffusion-weighted scans to quantify white matter microstructure with free-water imaging. We also investigated the interaction of biological sex and ASD diagnosis on these measures, and their correlation with clinical scales of RRB.

RESULTS

Individuals with ASD had significantly lower free-water corrected fractional anisotropy (FAT) and higher free-water (FW) in cortico-basal ganglia white matter tracts. These microstructural differences did not interact with biological sex. Moreover, both FAT and FW in basal ganglia white matter tracts significantly correlated with measures of RRB. In contrast, we found no significant difference in basal ganglia or cerebellar gray matter volumes.

LIMITATIONS

The basal ganglia and cerebellar regions in this study were selected due to their hypothesized relevance to RRB. Differences between ASD and TD individuals that may occur outside the basal ganglia and cerebellum, and their potential relationship to RRB, were not evaluated.

CONCLUSIONS

These new findings demonstrate that cortico-basal ganglia white matter microstructure is altered in ASD and linked to RRB. FW in cortico-basal ganglia and intra-basal ganglia white matter was more sensitive to group differences in ASD, whereas cortico-basal ganglia FAT was more closely linked to RRB. In contrast, basal ganglia and cerebellar volumes did not differ in ASD. There was no interaction between ASD diagnosis and sex-related differences in brain structure. Future diffusion imaging investigations in ASD may benefit from free-water estimation and correction in order to better understand how white matter is affected in ASD, and how such measures are linked to RRB.

Study protocol - elucidating the neural correlates of functional remediation for older adults with bipolar disorder

Introduction

Beyond mood abnormalities, bipolar disorder (BD) includes cognitive impairments that worsen psychosocial functioning and quality of life. These deficits are especially severe in older adults with BD (OABD), a condition expected to represent most individuals with BD in the upcoming years. Restoring the psychosocial functioning of this population will thus soon represent a public health priority. To help tackle the problem, the Bipolar and Depressive Disorders Unit at the Hospital Clínic of Barcelona has recently adapted its Functional Remediation (FR) program to that population, calling it FROA-BD. However, while scarce previous studies localize the neural mechanisms of cognitive remediation interventions in the dorsal prefrontal cortex, the specific mechanisms are seldom unknown. In the present project, we will investigate the neural correlates of FR-OABD to understand its mechanisms better and inform for potential optimization. The aim is to investigate the brain features and changes associated with FROA-BD efficacy.

Methods

Thirty-two individuals with OABD in full or partial remission will undergo a magnetic resonance imaging (MRI) session before receiving FR-OABD. After completing the FR-OABD intervention, they will undergo another MRI session. The MRI sessions will include structural, diffusion-weighted imaging (DWI), functional MRI (fMRI) with working memory (n-back) and verbal learning tasks, and frontal spectroscopy. We will correlate the pre-post change in dorsolateral and dorsomedial prefrontal cortices activation during the n-back task with the change in psychosocial functioning [measured with the Functioning Assessment Short Test (FAST)]. We will also conduct exploratory whole-brain correlation analyses between baseline or pre-post changes in MRI data and other clinical and cognitive outcomes to provide more insights into the mechanisms and explore potential brain markers that may predict a better treatment response. We will also conduct separate analyses by sex.

Discussion

The results of this study may provide insights into how FROA-BD and other cognitive remediations modulate brain function and thus could optimize these interventions.

Spatiotemporal patterns of cortical microstructural maturation in children and adolescents with diffusion MRI

Neocortical maturation is a dynamic process that proceeds in a hierarchical manner; however, the spatiotemporal organization of cortical microstructure with diffusion MRI has yet to be fully defined. This study characterized cortical microstructural maturation using diffusion MRI (fwe-diffusion tensor imaging [DTI] and neurite orientation dispersion and density imaging [NODDI] multicompartment modeling) in a cohort of 637 children and adolescents between 8 and 21 years of age. We found spatially heterogeneous developmental patterns broadly demarcated into functional domains where NODDI metrics increased, and fwe-DTI metrics decreased with age. By applying nonlinear growth models in a vertex-wise analysis, we observed a general posterior-to-anterior pattern of maturation, where the fwe-DTI measures mean diffusivity and radial diffusivity reached peak maturation earlier than the NODDI metrics neurite density index. Using non-negative matrix factorization, we found occipito-parietal cortical regions that correspond to lower order sensory domains mature earlier than fronto-temporal higher order association domains. Our findings corroborate previous histological and neuroimaging studies that show spatially varying patterns of cortical maturation that may reflect unique developmental processes of cytoarchitectonically determined regional patterns of change.

Age-associated alterations in thalamocortical structural connectivity in youths with a psychosis-spectrum disorder

Psychotic symptoms typically emerge in adolescence. Age-associated thalamocortical connectivity differences in psychosis remain unclear. We analyzed diffusion-weighted imaging data from 1254 participants 8-23 years old (typically developing (TD):N = 626, psychosis-spectrum (PS): N = 329, other psychopathology (OP): N = 299) from the Philadelphia Neurodevelopmental Cohort. We modeled thalamocortical tracts using deterministic fiber tractography, extracted Q-Space Diffeomorphic Reconstruction (QSDR) and diffusion tensor imaging (DTI) measures, and then used generalized additive models to determine group and age-associated thalamocortical connectivity differences. Compared to other groups, PS exhibited thalamocortical reductions in QSDR global fractional anisotropy (GFA, p-values range = 3.0 × 10-6-0.05) and DTI fractional anisotropy (FA, p-values range = 4.2 × 10-4-0.03). Compared to TD, PS exhibited shallower thalamus-prefrontal age-associated increases in GFA and FA during mid-childhood, but steeper age-associated increases during adolescence. TD and OP exhibited decreases in thalamus-frontal mean and radial diffusivities during adolescence; PS did not. Altered developmental trajectories of thalamocortical connectivity may contribute to the disruptions observed in adults with psychosis.

Age-related, multivariate associations between white matter microstructure and behavioral performance in three executive function domains

The executive function (EF) domains of working memory (WM), response inhibition (RI), and set shifting (SS) show maturational gains and are linked to neuroimaging-measured brain changes. This study explored ways in which maturation-linked differences in EF abilities are systematically associated with white matter microstructural differences from adolescence into young adulthood. Diffusion tensor imaging (DTI) and nine neurocognitive tests were collected from 120 healthy subjects ages 12-24. Analyses across the white matter skeleton were performed, focusing on fractional anisotropy (FA). Data were 'fused' using a multivariate technique (CCA+jICA), producing four independent components (ICs) depicting white matter FA values that covaried with test performance. Correlations between age and IC loading coefficients identified three EF-DTI profiles that may change developmentally. In one, SS performance was linked to greater reliance on the FA of ventral brain tracts, and less on dorsal tracts with age. In another, white matter microstructure was related to a pattern of strong WM and weak SS that became more pronounced with age. A final IC revealed that younger individuals with low RI and high WM/SS skills typically matured out of this cognitive imbalance, underscored by white matter changes with age. These novel multivariate results begin to emphasize the complexity of brain structure-cognition relationships in adolescents and young adults.

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.

Connectomes for 40,000 UK Biobank participants: A multi-modal, multi-scale brain network resource

We mapped functional and structural brain networks for more than 40,000 UK Biobank participants. Structural connectivity was estimated with tractography and diffusion MRI. Resting-state functional MRI was used to infer regional functional connectivity. We provide high-quality structural and functional connectomes for multiple parcellation granularities, several alternative measures of interregional connectivity, and a variety of common data pre-processing techniques, yielding more than one million connectomes in total and requiring more than 200,000 h of compute time. For a single subject, we provide 28 out-of-the-box versions of structural and functional brain networks, allowing users to select, e.g., the parcellation and connectivity measure that best suit their research goals. Furthermore, we provide code and intermediate data for the time-efficient reconstruction of more than 1000 different versions of a subject's connectome based on an array of methodological choices. All connectomes are available via the UK Biobank data-sharing platform and our connectome mapping pipelines are openly available. In this report, we describe our connectome resource in detail for users, outline key considerations in developing an efficient pipeline to map an unprecedented number of connectomes, and report on the quality control procedures that were completed to ensure connectome reliability and accuracy. We demonstrate that our structural and functional connectivity matrices meet a number of quality control checks and replicate previously established findings in network neuroscience. We envisage that our resource will enable new studies of the human connectome in health, disease, and aging at an unprecedented scale.

Glymphatic system impairment in sleep disruption: diffusion tensor image analysis along the perivascular space (DTI-ALPS)

PURPOSE

This study aimed to evaluate the relationship between sleep quality as assessed using the Pittsburgh Sleep Quality Index (PSQI) and the index of diffusivity along the perivascular space (ALPS index), a possible indirect indicator of glymphatic system activity.

MATERIALS AND METHODS

This study included the diffusion magnetic resonance imaging (MRI) data of 317 people with sleep disruption and 515 healthy controls (HCs) from the Human Connectome Project (WU-MINN HCP 1200). The ALPS index was calculated automatically based on diffusion tensor image analysis (DTI)-ALPS of diffusion MRI. The ALPS index of the sleep disruption and HC groups was compared using general linear model (GLM) analysis with covariates, such as age, sex, level of education, and intracranial volume. In addition, to confirm the relationship between sleep quality and the ALPS index in the sleep disruption group as well as evaluate the effect of each PSQI component on the ALPS index, correlation analyses between the ALPS indices and PSQI scores of all the components and between the ALPS index and each PSQI component was performed using GLM analysis with the abovementioned covariates, respectively.

RESULTS

The ALPS index was significantly lower in the sleep disruption group than in the HC group (p = 0.001). Moreover, the ALPS indices showed significant negative correlations with the PSQI scores of all the components (false discovery rate [FDR]-corrected p < 0.001). Two significant negative correlations were also found between the ALPS index and PSQI component 2 (sleep latency, FDR-corrected p < 0.001) and 6 (the use of sleep medication, FDR-corrected p < 0.001).

CONCLUSION

Our findings suggest that glymphatic system impairment contributes to sleep disruption in young adults.

Structural connectivity of the multiple demand network in humans and comparison to the macaque brain

Fluid intelligence encompasses a wide range of abilities such as working memory, problem-solving, and relational reasoning. In the human brain, these abilities are associated with the Multiple Demand Network, traditionally thought to involve combined activity of specific regions predominantly in the prefrontal and parietal cortices. However, the structural basis of the interactions between areas in the Multiple Demand Network, as well as their evolutionary basis among primates, remains largely unexplored. Here, we exploit diffusion MRI to elucidate the major white matter pathways connecting areas of the human core and extended Multiple Demand Network. We then investigate whether similar pathways can be identified in the putative homologous areas of the Multiple Demand Network in the macaque monkey. Finally, we contrast human and monkey networks using a recently proposed approach to compare different species' brains within a common organizational space. Our results indicate that the core Multiple Demand Network relies mostly on dorsal longitudinal connections and, although present in the macaque, these connections are more pronounced in the human brain. The extended Multiple Demand Network relies on distinct pathways and communicates with the core Multiple Demand Network through connections that also appear enhanced in the human compared with the macaque.

DENOISING DIFFUSION MRI: CONSIDERATIONS AND IMPLICATIONS FOR ANALYSIS

Development of diffusion MRI (dMRI) denoising approaches has experienced considerable growth over the last years. As noise can inherently reduce accuracy and precision in measurements, its effects have been well characterised both in terms of uncertainty increase in dMRI-derived features and in terms of biases caused by the noise floor, the smallest measurable signal given the noise level. However, gaps in our knowledge still exist in objectively characterising dMRI denoising approaches in terms of both of these effects and assessing their efficacy. In this work, we reconsider what a denoising method should and should not do and we accordingly define criteria to characterise the performance. We propose a comprehensive set of evaluations, including i) benefits in improving signal quality and reducing noise variance, ii) gains in reducing biases and the noise floor and improving, iii) preservation of spatial resolution, iv) agreement of denoised data against a gold standard, v) gains in downstream parameter estimation (precision and accuracy), vi) efficacy in enabling noise-prone applications, such as ultra-high-resolution imaging. We further provide newly acquired complex datasets (magnitude and phase) with multiple repeats that sample different SNR regimes to highlight performance differences under different scenarios. Without loss of generality, we subsequently apply a number of exemplar patch-based denoising algorithms to these datasets, including Non-Local Means, Marchenko-Pastur PCA (MPPCA) in the magnitude and complex domain and NORDIC, and compare them with respect to the above criteria and against a gold standard complex average of multiple repeats. We demonstrate that all tested denoising approaches reduce noise-related variance, but not always biases from the elevated noise floor. They all induce a spatial resolution penalty, but its extent can vary depending on the method and the implementation. Some denoising approaches agree with the gold standard more than others and we demonstrate challenges in even defining such a standard. Overall, we show that dMRI denoising performed in the complex domain is advantageous to magnitude domain denoising with respect to all the above criteria.

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.