The conduction velocity (CV) of action potentials along axons is a key neurophysiological property central to neural communication. The ability to estimate CV in humans in vivo from non-invasive MRI ...methods would therefore represent a significant advance in neuroscience. However, there are two major challenges that this paper aims to address: (1) Much of the complexity of the neurophysiology of action potentials cannot be captured with currently available MRI techniques. Therefore, we seek to establish the variability in CV that can be captured when predicting CV purely from parameters that have been reported to be estimatable from MRI: inner axon diameter (AD) and g-ratio. (2) errors inherent in existing MRI-based biophysical models of tissue will propagate through to estimates of CV, the extent to which is currently unknown. Issue (1) is investigated by performing a sensitivity analysis on a comprehensive model of axon electrophysiology and determining the relative sensitivity to various morphological and electrical parameters. The investigations suggest that 85% of the variance in CV is accounted for by variation in AD and g-ratio. The observed dependency of CV on AD and g-ratio is well characterised by the previously reported model by Rushton. Issue (2) is investigated through simulation of diffusion and relaxometry MRI data for a range of axon morphologies, applying models of restricted diffusion and relaxation processes to derive estimates of axon volume fraction (AVF), AD and g-ratio and estimating CV from the derived parameters. The results show that errors in the AVF have the biggest detrimental impact on estimates of CV, particularly for sparse fibre populations (AVF<0.3). For our equipment set-up and acquisition protocol, CV estimates are most accurate (below 5% error) where AVF is above 0.3, g-ratio is between 0.6 and 0.85 and AD is high (above 4μm). CV estimates are robust to errors in g-ratio estimation but are highly sensitive to errors in AD estimation, particularly where ADs are small. We additionally show CV estimates in human corpus callosum in a small number of subjects. In conclusion, we demonstrate accurate CV estimates are possible in regions of the brain where AD is sufficiently large. Problems with estimating ADs for smaller axons presents a problem for estimating CV across the whole CNS and should be the focus of further study.
•85% of the variance in CV is accounted for by axon diameter and g-ratio, which are potentially accessible from in vivo MRI.•CV estimates from MRI are robust to errors in myelin and axonal volume estimates, but sensitive to errors in axon diameter.•CV estimates are feasible for large axons but limitations of in vivo imaging of small axons poses a significant challenge.
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 B 0 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.
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.
Display omitted
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.
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.
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.
Display omitted
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.
Abstract Background Recent genome-wide association studies have identified genetic loci that jointly make a considerable contribution to risk of developing Alzheimer’s disease (AD). Because ...neuropathological features of AD can be present several decades before disease onset, we investigated whether effects of polygenic risk are detectable by neuroimaging in young adults. We hypothesised that higher polygenic risk scores (PRS) for AD would be associated with reduced volume of the hippocampus, and other limbic and paralimbic areas. We further hypothesised that AD PRS would affect the microstructure of fibre tracts connecting the hippocampus with other brain areas. Methods We analysed the association between AD PRS and brain imaging parameters using T1-weighted structural (n=272) and diffusion-weighted scans (n=197). Results We found a significant association between AD PRS and left hippocampal volume, with higher risk associated with lower left hippocampal volume (p=0.001). This effect remained when the APOE gene was excluded (p=0.031), suggesting that the relationship between hippocampal volume and AD is the result of multiple genetic factors, not exclusively variability in the APOE gene. The diffusion tensor imaging analysis revealed that fractional anisotropy of the right cingulum was inversely correlated with AD PRS (p=0.009). We thus show that polygenic effects of AD risk variants on brain structure can already be detected in young adults. Conclusions This finding paves the way for further investigation of the effects of AD risk variants and may become useful for efforts to combine genotypic and phenotypic data for risk prediction and to enrich future prevention trials of AD.
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.
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.
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.