Improving Fiber Alignment in HARDI by Combining Contextual PDE Flow with Constrained Spherical Deconvolution

Microstructural mapping of neural pathways in Alzheimer's disease using macrostructure-informed normative tractometry

INTRODUCTION

Diffusion-weighted magnetic resonance imaging (dMRI) is sensitive to the microstructural properties of brain tissues and shows great promise in detecting the effects of degenerative diseases. However, many approaches analyze single measures averaged over regions of interest without considering the underlying fiber geometry.

METHODS

We propose a novel macrostructure-informed normative tractometry (MINT) framework to investigate how white matter (WM) microstructure and macrostructure are jointly altered in mild cognitive impairment (MCI) and dementia. We compared MINT-derived metrics with univariate diffusion tensor imaging (DTI) metrics to examine how fiber geometry may impact the interpretation of microstructure.

RESULTS

In two multisite cohorts from North America and India, we find consistent patterns of microstructural and macrostructural anomalies implicated in MCI and dementia; we also rank diffusion metrics' sensitivity to dementia.

DISCUSSION

We show that MINT, by jointly modeling tract shape and microstructure, has the potential to disentangle and better interpret the effects of degenerative disease on the brain's neural pathways.

HIGHLIGHTS

Changes in diffusion tensor imaging metrics may be due to macroscopic changes. Normative models encode normal variability of diffusion metrics in healthy controls. Variational autoencoder applied on tractography can learn patterns of fiber geometry. WM microstructure and macrostructure are modeled with multivariate methods. Transfer learning uses pretraining and fine-tuning for increased efficiency.

BundleAGE: Predicting White Matter Age using Along-Tract Microstructural Profiles from Diffusion MRI

Brain Age Gap Estimation (BrainAGE) is an estimate of the gap between a person's chronological age (CA) and a measure of their brain's 'biological age' (BA). This metric is often used as a marker of accelerated aging, albeit with some caveats. Age prediction models trained on brain structural and functional MRI have been employed to derive BrainAGE biomarkers, for predicting the risk of neurodegeneration. While voxel-based and along-tract microstructural maps from diffusion MRI have been used to study brain aging, no studies have evaluated along-tract microstructure for computing BrainAGE. In this study, we train machine learning models to predict a person's age using along-tract microstructural profiles from diffusion tensor imaging. We were able to demonstrate differential aging patterns across different white matter bundles and microstructural measures. The novel Bundle Age Gap Estimation (BundleAGE) biomarker shows potential in quantifying risk factors for neurodegenerative diseases and aging, while incorporating finer scale information throughout white matter bundles.

Deep Normative Tractometry for Identifying Joint White Matter Macro- and Micro-structural Abnormalities in Alzheimer's Disease

This study introduces the Deep Normative Tractometry (DNT) framework, that encodes the joint distribution of both macrostructural and microstructural profiles of the brain white matter tracts through a variational autoencoder (VAE). By training on data from healthy controls, DNT learns the normative distribution of tract data, and can delineate along-tract micro- and macro-structural abnormalities. Leveraging a large sample size via generative pre-training, we assess DNT's generalizability using transfer learning on data from an independent cohort acquired in India. Our findings demonstrate DNT's capacity to detect widespread diffusivity abnormalities along tracts in mild cognitive impairment and Alzheimer's disease, aligning closely with results from the Bundle Analytics (BUAN) tractometry pipeline. By incorporating tract geometry information, DNT may be able to distinguish disease-related abnormalities in anisotropy from tract macrostructure, and shows promise in enhancing fine-scale mapping and detection of white matter alterations in neurodegenerative conditions.

Microstructural Mapping of Neural Pathways in Alzheimer's Disease using Macrostructure-Informed Normative Tractometry

Introduction

Diffusion MRI is sensitive to the microstructural properties of brain tissues, and shows great promise in detecting the effects of degenerative diseases. However, many approaches analyze single measures averaged over regions of interest, without considering the underlying fiber geometry.

Methods

Here, we propose a novel Macrostructure-Informed Normative Tractometry (MINT) framework, to investigate how white matter microstructure and macrostructure are jointly altered in mild cognitive impairment (MCI) and dementia. We compare MINT-derived metrics with univariate metrics from diffusion tensor imaging (DTI), to examine how fiber geometry may impact interpretation of microstructure.

Results

In two multi-site cohorts from North America and India, we find consistent patterns of microstructural and macrostructural anomalies implicated in MCI and dementia; we also rank diffusion metrics' sensitivity to dementia.

Discussion

We show that MINT, by jointly modeling tract shape and microstructure, has potential to disentangle and better interpret the effects of degenerative disease on the brain's neural pathways.

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).

Deep Normative Tractometry for Identifying Joint White Matter Macro- and Micro-structural Abnormalities in Alzheimer's Disease

This study introduces the Deep Normative Tractometry (DNT) framework, that encodes the joint distribution of both macrostructural and microstructural profiles of the brain white matter tracts through a variational autoencoder (VAE). By training on data from healthy controls, DNT learns the normative distribution of tract data, and can delineate along-tract micro-and macro-structural abnormalities. Leveraging a large sample size via generative pre-training, we assess DNT's generalizability using transfer learning on data from an independent cohort acquired in India. Our findings demonstrate DNT's capacity to detect widespread diffusivity abnormalities along tracts in mild cognitive impairment and Alzheimer's disease, aligning closely with results from the Bundle Analytics (BUAN) tractometry pipeline. By incorporating tract geometry information, DNT may be able to distinguish disease-related abnormalities in anisotropy from tract macrostructure, and shows promise in enhancing fine-scale mapping and detection of white matter alterations in neurodegenerative conditions.

Accelerated Brain Atrophy, Microstructural Decline and Connectopathy in Age-Related Macular Degeneration

Age-related macular degeneration (AMD) has recently been linked to cognitive impairment. We hypothesized that AMD modifies the brain aging trajectory, and we conducted a longitudinal diffusion MRI study on 40 participants (20 with AMD and 20 controls) to reveal the location, extent, and dynamics of AMD-related brain changes. Voxel-based analyses at the first visit identified reduced volume in AMD participants in the cuneate gyrus, associated with vision, and the temporal and bilateral cingulate gyrus, linked to higher cognition and memory. The second visit occurred 2 years after the first and revealed that AMD participants had reduced cingulate and superior frontal gyrus volumes, as well as lower fractional anisotropy (FA) for the bilateral occipital lobe, including the visual and the superior frontal cortex. We detected faster rates of volume and FA reduction in AMD participants in the left temporal cortex. We identified inter-lingual and lingual-cerebellar connections as important differentiators in AMD participants. Bundle analyses revealed that the lingual gyrus had a lower streamline length in the AMD participants at the first visit, indicating a connection between retinal and brain health. FA differences in select inter-lingual and lingual cerebellar bundles at the second visit showed downstream effects of vision loss. Our analyses revealed widespread changes in AMD participants, beyond brain networks directly involved in vision processing.

An atlas of white matter anatomy, its variability, and reproducibility based on constrained spherical deconvolution of diffusion MRI

Virtual dissection of white matter (WM) using diffusion MRI tractography is confounded by its poor reproducibility. Despite the increased adoption of advanced reconstruction models, early region-of-interest driven protocols based on diffusion tensor imaging (DTI) remain the dominant reference for virtual dissection protocols. Here we bridge this gap by providing a comprehensive description of typical WM anatomy reconstructed using a reproducible automated subject-specific parcellation-based approach based on probabilistic constrained-spherical deconvolution (CSD) tractography. We complement this with a WM template in MNI space comprising 68 bundles, including all associated anatomical tract selection labels and associated automated workflows. Additionally, we demonstrate bundle inter- and intra-subject variability using 40 (20 test-retest) datasets from the human connectome project (HCP) and 5 sessions with varying b-values and number of b-shells from the single-subject Multiple Acquisitions for Standardization of Structural Imaging Validation and Evaluation (MASSIVE) dataset. The most reliably reconstructed bundles were the whole pyramidal tracts, primary corticospinal tracts, whole superior longitudinal fasciculi, frontal, parietal and occipital segments of the corpus callosum and middle cerebellar peduncles. More variability was found in less dense bundles, e.g., the fornix, dentato-rubro-thalamic tract (DRTT), and premotor pyramidal tract. Using the DRTT as an example, we show that this variability can be reduced by using a higher number of seeding attempts. Overall inter-session similarity was high for HCP test-retest data (median weighted-dice = 0.963, stdev = 0.201 and IQR = 0.099). Compared to the HCP-template bundles there was a high level of agreement for the HCP test-retest data (median weighted-dice = 0.747, stdev = 0.220 and IQR = 0.277) and for the MASSIVE data (median weighted-dice = 0.767, stdev = 0.255 and IQR = 0.338). In summary, this WM atlas provides an overview of the capabilities and limitations of automated subject-specific probabilistic CSD tractography for mapping white matter fasciculi in healthy adults. It will be most useful in applications requiring a reproducible parcellation-based dissection protocol, and as an educational resource for applied neuroimaging and clinical professionals.

Fourier Transform on the Homogeneous Space of 3D Positions and Orientations for Exact Solutions to Linear PDEs

Fokker–Planck PDEs (including diffusions) for stable Lévy processes (including Wiener processes) on the joint space of positions and orientations play a major role in mechanics, robotics, image analysis, directional statistics and probability theory. Exact analytic designs and solutions are known in the 2D case, where they have been obtained using Fourier transform on SE(2). Here, we extend these approaches to 3D using Fourier transform on the Lie group SE(3) of rigid body motions. More precisely, we define the homogeneous space of 3D positions and orientations R3⋊S2:=SE(3)/({0}×SO(2)) as the quotient in SE(3). In our construction, two group elements are equivalent if they are equal up to a rotation around the reference axis. On this quotient, we design a specific Fourier transform. We apply this Fourier transform to derive new exact solutions to Fokker–Planck PDEs of α-stable Lévy processes on R3⋊S2. This reduces classical analysis computations and provides an explicit algebraic spectral decomposition of the solutions. We compare the exact probability kernel for α=1 (the diffusion kernel) to the kernel for α=12 (the Poisson kernel). We set up stochastic differential equations (SDEs) for the Lévy processes on the quotient and derive corresponding Monte-Carlo methods. We verified that the exact probability kernels arise as the limit of the Monte-Carlo approximations.

Bundle-specific tractography with incorporated anatomical and orientational priors

Anatomical white matter bundles vary in shape, size, length, and complexity, making diffusion MRI tractography reconstruction of some bundles more difficult than others. As a result, bundles reconstruction often suffers from a poor spatial extent recovery. To fill-up the white matter volume as much and as best as possible, millions of streamlines can be generated and filtering techniques applied to address this issue. However, well-known problems and biases are introduced such as the creation of a large number of false positives and over-representation of easy-to-track parts of bundles and under-representation of hard-to-track. To address these challenges, we developed a Bundle-Specific Tractography (BST) algorithm. It incorporates anatomical and orientational prior knowledge during the process of streamline tracing to increase reproducibility, sensitivity, specificity and efficiency when reconstructing certain bundles of interest. BST outperforms classical deterministic, probabilistic, and global tractography methods. The increase in anatomically plausible streamlines, with larger spatial coverage, helps to accurately represent the full shape of bundles, which could greatly enhance and robustify tract-based and connectivity-based neuroimaging studies.

Design and Processing of Invertible Orientation Scores of 3D Images

The enhancement and detection of elongated structures in noisy image data are relevant for many biomedical imaging applications. To handle complex crossing structures in 2D images, 2D orientation scores U : R 2 × S 1 → C were introduced, which already showed their use in a variety of applications. Here we extend this work to 3D orientation scores U : R 3 × S 2 → C . First, we construct the orientation score from a given dataset, which is achieved by an invertible coherent state type of transform. For this transformation we introduce 3D versions of the 2D cake wavelets, which are complex wavelets that can simultaneously detect oriented structures and oriented edges. Here we introduce two types of cake wavelets: the first uses a discrete Fourier transform, and the second is designed in the 3D generalized Zernike basis, allowing us to calculate analytical expressions for the spatial filters. Second, we propose a nonlinear diffusion flow on the 3D roto-translation group: crossing-preserving coherence-enhancing diffusion via orientation scores (CEDOS). Finally, we show two applications of the orientation score transformation. In the first application we apply our CEDOS algorithm to real medical image data. In the second one we develop a new tubularity measure using 3D orientation scores and apply the tubularity measure to both artificial and real medical data.

Stability metrics for optic radiation tractography: Towards damage prediction after resective surgery

BACKGROUND

An accurate delineation of the optic radiation (OR) using diffusion MR tractography may reduce the risk of a visual field deficit after temporal lobe resection. However, tractography is prone to generate spurious streamlines, which deviate strongly from neighboring streamlines and hinder a reliable distance measurement between the temporal pole and the Meyer's loop (ML-TP distance).

NEW METHOD

Stability metrics are introduced for the automated removal of spurious streamlines near the Meyer's loop. Firstly, fiber-to-bundle coherence (FBC) measures can identify spurious streamlines by estimating their alignment with the surrounding streamline bundle. Secondly, robust threshold selection removes spurious streamlines while preventing an underestimation of the extent of the Meyer's loop. Standardized parameter selection is realized through test-retest evaluation of the variability in ML-TP distance.

RESULTS

The variability in ML-TP distance after parameter selection was below 2mm for each of the healthy volunteers studied (N=8). The importance of the stability metrics is illustrated for epilepsy surgery candidates (N=3) for whom the damage to the Meyer's loop was evaluated by comparing the pre- and post-operative OR reconstruction. The difference between predicted and observed damage is in the order of a few millimeters, which is the error in measured ML-TP distance.

COMPARISON WITH EXISTING METHOD(S)

The stability metrics are a novel method for the robust estimate of the ML-TP distance.

CONCLUSIONS

The stability metrics are a promising tool for clinical trial studies, in which the damage to the OR can be related to the visual field deficit that may occur after epilepsy surgery.

The visual white matter: The application of diffusion MRI and fiber tractography to vision science

Visual neuroscience has traditionally focused much of its attention on understanding the response properties of single neurons or neuronal ensembles. The visual white matter and the long-range neuronal connections it supports are fundamental in establishing such neuronal response properties and visual function. This review article provides an introduction to measurements and methods to study the human visual white matter using diffusion MRI. These methods allow us to measure the microstructural and macrostructural properties of the white matter in living human individuals; they allow us to trace long-range connections between neurons in different parts of the visual system and to measure the biophysical properties of these connections. We also review a range of findings from recent studies on connections between different visual field maps, the effects of visual impairment on the white matter, and the properties underlying networks that process visual information supporting visual face recognition. Finally, we discuss a few promising directions for future studies. These include new methods for analysis of MRI data, open datasets that are becoming available to study brain connectivity and white matter properties, and open source software for the analysis of these data.