Surface-based tracking for short association fibre tractography.

It is estimated that in the human brain, short association fibres (SAF) represent more than half of the total white matter volume and their involvement has been implicated in a range of neurological and psychiatric conditions. This population of fibres, however, remains relatively understudied in the neuroimaging literature. Some of the challenges pertinent to the mapping of SAF include their variable anatomical course and proximity to the cortical mantle, leading to partial volume effects and potentially affecting streamline trajectory estimation. This work considers the impact of seeding and filtering strategies and choice of scanner, acquisition, data resampling to propose a whole-brain, surface-based short (≤30-40 mm) SAF tractography approach. The framework is shown to produce longer streamlines with a predilection for connecting gyri as well as high cortical coverage. We further demonstrate that certain areas of subcortical white matter become disproportionally underrepresented in diffusion-weighted MRI data with lower angular and spatial resolution and weaker diffusion weighting; however, collecting data with stronger gradients than are usually available clinically has minimal impact, making our framework translatable to data collected on commonly available hardware. Finally, the tractograms are examined using voxel- and surface-based measures of consistency, demonstrating moderate reliability, low repeatability and high between-subject variability, urging caution when streamline count-based analyses of SAF are performed.

Mutation-related magnetization-transfer, not axon density, drives white matter differences in premanifest Huntington disease: Evidence from in vivo ultra-strong gradient MRI.

White matter (WM) alterations have been observed in Huntington disease (HD) but their role in the disease-pathophysiology remains unknown. We assessed WM changes in premanifest HD by exploiting ultra-strong-gradient magnetic resonance imaging (MRI). This allowed to separately quantify magnetization transfer ratio (MTR) and hindered and restricted diffusion-weighted signal fractions, and assess how they drove WM microstructure differences between patients and controls. We used tractometry to investigate region-specific alterations across callosal segments with well-characterized early- and late-myelinating axon populations, while brain-wise differences were explored with tract-based cluster analysis (TBCA). Behavioral measures were included to explore disease-associated brain-function relationships. We detected lower MTR in patients' callosal rostrum (tractometry: p = .03; TBCA: p = .03), but higher MTR in their splenium (tractometry: p = .02). Importantly, patients' mutation-size and MTR were positively correlated (all p-values < .01), indicating that MTR alterations may directly result from the mutation. Further, MTR was higher in younger, but lower in older patients relative to controls (p = .003), suggesting that MTR increases are detrimental later in the disease. Finally, patients showed higher restricted diffusion signal fraction (FR) from the composite hindered and restricted model of diffusion (CHARMED) in the cortico-spinal tract (p = .03), which correlated positively with MTR in the posterior callosum (p = .033), potentially reflecting compensatory mechanisms. In summary, this first comprehensive, ultra-strong gradient MRI study in HD provides novel evidence of mutation-driven MTR alterations at the premanifest disease stage which may reflect neurodevelopmental changes in iron, myelin, or a combination of these.

MICRA: Microstructural image compilation with repeated acquisitions.

We provide a rich multi-contrast microstructural MRI dataset acquired on an ultra-strong gradient 3T Connectom MRI scanner comprising 5 repeated sets of MRI microstructural contrasts in 6 healthy human participants. The availability of data sets that support comprehensive simultaneous assessment of test-retest reliability of multiple microstructural contrasts (i.e., those derived from advanced diffusion, multi-component relaxometry and quantitative magnetisation transfer MRI) in the same population is extremely limited. This unique dataset is offered to the imaging community as a test-bed resource for conducting specialised analyses that may assist and inform their current and future research. The Microstructural Image Compilation with Repeated Acquisitions (MICRA) dataset includes raw data and computed microstructure maps derived from multi-shell and multi-direction encoded diffusion, multi-component relaxometry and quantitative magnetisation transfer acquisition protocols. Our data demonstrate high reproducibility of several microstructural MRI measures across scan sessions as shown by intra-class correlation coefficients and coefficients of variation. To illustrate a potential use of the MICRA dataset, we computed sample sizes required to provide sufficient statistical power a priori across different white matter pathways and microstructure measures for different statistical comparisons. We also demonstrate whole brain white matter voxel-wise repeatability in several microstructural maps. The MICRA dataset will be of benefit to researchers wishing to conduct similar reliability tests, power estimations or to evaluate the robustness of their own analysis pipelines.

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data.

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson-Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b-matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b-values in contrast to the perhaps common assumption that only high b-value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300 mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

A comparative study of gradient nonlinearity correction strategies for processing diffusion data obtained with ultra-strong gradient MRI scanners.

PURPOSE: The analysis of diffusion data obtained under large gradient nonlinearities necessitates corrections during data reconstruction and analysis. While two such preprocessing pipelines have been proposed, no comparative studies assessing their performance exist. Furthermore, both pipelines neglect the impact of subject motion during acquisition, which, in the presence of gradient nonlinearities, induces spatio-temporal B-matrix variations. Here, spatio-temporal B-matrix tracking (STB) is proposed and its performance compared to established pipelines. METHODS: Diffusion tensor MRI (DT-MRI) was performed using a 300 mT/m gradient system. Data were acquired with volunteers positioned in regions with pronounced gradient nonlinearities, and used to compare the performance of six different processing pipelines, including STB. RESULTS: Up to 30% errors were observed in DT-MRI parameter estimates when neglecting gradient nonlinearities. Moreover, the order in which B0 inhomogeneity, eddy current and gradient nonlinearity corrections were performed was found to impact the consistency of parameter estimates significantly. Although, no pipeline emerged as a clear winner, the STB approach seemed to yield the most consistent parameter estimates under large gradient nonlinearities. CONCLUSIONS: Under large gradient nonlinearities, the choice of preprocessing pipeline significantly impacts the estimated diffusion parameters. Motion-induced spatio-temporal B-matrix variations can lead to systematic bias in the parameter estimates, that can be ameliorated using the proposed STB framework.

Mammillothalamic Disconnection Alters Hippocampocortical Oscillatory Activity and Microstructure: Implications for Diencephalic Amnesia.

Diencephalic amnesia can be as debilitating as the more commonly known temporal lobe amnesia, yet the precise contribution of diencephalic structures to memory processes remains elusive. Across four cohorts of male rats, we used discrete lesions of the mammillothalamic tract to model aspects of diencephalic amnesia and assessed the impact of these lesions on multiple measures of activity and plasticity within the hippocampus and retrosplenial cortex. Lesions of the mammillothalamic tract had widespread indirect effects on hippocampocortical oscillatory activity within both theta and gamma bands. Both within-region oscillatory activity and cross-regional synchrony were altered. The network changes were state-dependent, displaying different profiles during locomotion and paradoxical sleep. Consistent with the associations between oscillatory activity and plasticity, complementary analyses using several convergent approaches revealed microstructural changes, which appeared to reflect a suppression of learning-induced plasticity in lesioned animals. Together, these combined findings suggest a mechanism by which damage to the medial diencephalon can impact upon learning and memory processes, highlighting an important role for the mammillary bodies in the coordination of hippocampocortical activity.SIGNIFICANCE STATEMENT Information flow within the Papez circuit is critical to memory. Damage to ascending mammillothalamic projections has consistently been linked to amnesia in humans and spatial memory deficits in animal models. Here we report on the changes in hippocampocortical oscillatory dynamics that result from chronic lesions of the mammillothalamic tract and demonstrate, for the first time, that the mammillary bodies, independently of the supramammillary region, contribute to frequency modulation of hippocampocortical theta oscillations. Consistent with the associations between oscillatory activity and plasticity, the lesions also result in a suppression of learning-induced plasticity. Together, these data support new functional models whereby mammillary bodies are important for coordinating hippocampocortical activity rather than simply being a relay of hippocampal information as previously assumed.

Tractography in the presence of multiple sclerosis lesions.

Accurate anatomical localisation of specific white matter tracts and the quantification of their tract-specific microstructural damage in conditions such as multiple sclerosis (MS) can contribute to a better understanding of symptomatology, disease evolution and intervention effects. Diffusion MRI-based tractography is being used increasingly to segment white matter tracts as regions-of-interest for subsequent quantitative analysis. Since MS lesions can interrupt the tractography algorithm's tract reconstruction, clinical studies frequently resort to atlas-based approaches, which are convenient but ignorant to individual variability in tract size and shape. Here, we revisit the problem of individual tractography in MS, comparing tractography algorithms using: (i) The diffusion tensor framework; (ii) constrained spherical deconvolution (CSD); and (iii) damped Richardson-Lucy (dRL) deconvolution. Firstly, using simulated and in vivo data from 29 MS patients and 19 healthy controls, we show that the three tracking algorithms respond differentially to MS pathology. While the tensor-based approach is unable to deal with crossing fibres, CSD produces spurious streamlines, in particular in tissue with high fibre loss and low diffusion anisotropy. With dRL, streamlines are increasingly interrupted in pathological tissue. Secondly, we demonstrate that despite the effects of lesions on the fibre orientation reconstruction algorithms, fibre tracking algorithms are still able to segment tracts that pass through areas with a high prevalence of lesions. Combining dRL-based tractography with an automated tract segmentation tool on data from 131 MS patients, the cortico-spinal tracts and arcuate fasciculi could be reconstructed in more than 90% of individuals. Comparing tract-specific microstructural parameters (fractional anisotropy, radial diffusivity and magnetisation transfer ratio) in individually segmented tracts to those from a tract probability map, we show that there is no systematic disease-related bias in the individually reconstructed tracts, suggesting that lesions and otherwise damaged parts are not systematically omitted during tractography. Thirdly, we demonstrate modest anatomical correspondence between the individual and tract probability-based approach, with a spatial overlap between 35 and 55%. Correlations between tract-averaged microstructural parameters in individually segmented tracts and the probability-map approach ranged between r=.53 (p<.001) for radial diffusivity in the right cortico-spinal tract and r=.97 (p<.001) for magnetisation transfer ratio in the arcuate fasciculi. Our results show that MS white matter lesions impact fibre orientation reconstructions but this does not appear to hinder the ability to anatomically reconstruct white matter tracts in MS. Individual tract segmentation in MS is feasible on a large scale and could prove a powerful tool for investigating diagnostic and prognostic markers.

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 â€‹mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies.

Impact of b-value on estimates of apparent fibre density.

Recent advances in diffusion magnetic resonance imaging (dMRI) analysis techniques have improved our understanding of fibre-specific variations in white matter microstructure. Increasingly, studies are adopting multi-shell dMRI acquisitions to improve the robustness of dMRI-based inferences. However, the impact of b-value choice on the estimation of dMRI measures such as apparent fibre density (AFD) derived from spherical deconvolution is not known. Here, we investigate the impact of b-value sampling scheme on estimates of AFD. First, we performed simulations to assess the correspondence between AFD and simulated intra-axonal signal fraction across multiple b-value sampling schemes. We then studied the impact of sampling scheme on the relationship between AFD and age in a developmental population (n = 78) aged 8-18 (mean = 12.4, SD = 2.9 years) using hierarchical clustering and whole brain fixel-based analyses. Multi-shell dMRI data were collected at 3.0T using ultra-strong gradients (300 mT/m), using 6 diffusion-weighted shells ranging from b = 0 to 6,000 s/mm2 . Simulations revealed that the correspondence between estimated AFD and simulated intra-axonal signal fraction was improved with high b-value shells due to increased suppression of the extra-axonal signal. These results were supported by in vivo data, as sensitivity to developmental age-relationships was improved with increasing b-value (b = 6,000 s/mm2 , median R2 = .34; b = 4,000 s/mm2 , median R2 = .29; b = 2,400 s/mm2 , median R2 = .21; b = 1,200 s/mm2 , median R2 = .17) in a tract-specific fashion. Overall, estimates of AFD and age-related microstructural development were better characterised at high diffusion-weightings due to improved correspondence with intra-axonal properties.

Dimensionality reduction of diffusion MRI measures for improved tractometry of the human brain.

Various diffusion MRI (dMRI) measures have been proposed for characterising tissue microstructure over the last 15 years. Despite the growing number of experiments using different dMRI measures in assessments of white matter, there has been limited work on: 1) examining their covariance along specific pathways; and on 2) combining these different measures to study tissue microstructure. Indeed, it quickly becomes intractable for existing analysis pipelines to process multiple measurements at each voxel and at each vertex forming a streamline, highlighting the need for new ways to visualise or analyse such high-dimensional data. In a sample of 36 typically developing children aged 8-18 years, we profiled various commonly used dMRI measures across 22 brain pathways. Using a data-reduction approach, we identified two biologically-interpretable components that capture 80% of the variance in these dMRI measures. The first derived component captures properties related to hindrance and restriction in tissue microstructure, while the second component reflects characteristics related to tissue complexity and orientational dispersion. We then demonstrate that the components generated by this approach preserve the biological relevance of the original measurements by showing age-related effects across developmentally sensitive pathways. In summary, our findings demonstrate that dMRI analyses can benefit from dimensionality reduction techniques, to help disentangling the neurobiological underpinnings of white matter organisation.