White matter tract microstructure, macrostructure, and associated cortical gray matter morphology across the lifespan.

Characterizing how, when and where the human brain changes across the lifespan is fundamental to our understanding of developmental processes of childhood and adolescence, degenerative processes of aging, and divergence from normal patterns in disease and disorders. We aimed to provide detailed descriptions of white matter pathways across the lifespan by thoroughly characterizing white matter microstructure, white matter macrostructure, and morphology of the cortex associated with white matter pathways. We analyzed 4 large, high-quality, publicly-available datasets comprising 2789 total imaging sessions, and participants ranging from 0 to 100 years old, using advanced tractography and diffusion modeling. We first find that all microstructural, macrostructural, and cortical features of white matter bundles show unique lifespan trajectories, with rates and timing of development and degradation that vary across pathways - describing differences between types of pathways and locations in the brain, and developmental milestones of maturation of each feature. Second, we show cross-sectional relationships between different features that may help elucidate biological changes occurring during different stages of the lifespan. Third, we show unique trajectories of age-associations across features. Finally, we find that age associations during development are strongly related to those during aging. Overall, this study reports normative data for several features of white matter pathways of the human brain that will be useful for studying normal and abnormal white matter development and degeneration.

Regional cortical thinning, demyelination and iron loss in cerebral small vessel disease.

The link between white matter hyperintensities (WMH) and cortical thinning is thought to be an important pathway by which WMH contributes to cognitive deficits in cerebral small vessel disease (SVD). However, the mechanism behind this association and the underlying tissue composition abnormalities are unclear. The objective of this study is to determine the association between WMH and cortical thickness, and the in vivo tissue composition abnormalities in the WMH-connected cortical regions. In this cross-sectional study, we included 213 participants with SVD who underwent standardized protocol including multimodal neuroimaging scans and cognitive assessment (i.e. processing speed, executive function and memory). We identified the cortex connected to WMH using probabilistic tractography starting from the WMH and defined the WMH-connected regions at three connectivity levels (low, medium and high connectivity level). We calculated the cortical thickness, myelin and iron of the cortex based on T1-weighted, quantitative R1, R2* and susceptibility maps. We used diffusion-weighted imaging to estimate the mean diffusivity of the connecting white matter tracts. We found that cortical thickness, R1, R2* and susceptibility values in the WMH-connected regions were significantly lower than in the WMH-unconnected regions (all Pcorrected < 0.001). Linear regression analyses showed that higher mean diffusivity of the connecting white matter tracts were related to lower thickness (ß = -0.30, Pcorrected < 0.001), lower R1 (ß = -0.26, Pcorrected = 0.001), lower R2* (ß = -0.32, Pcorrected < 0.001) and lower susceptibility values (ß = -0.39, Pcorrected < 0.001) of WMH-connected cortical regions at high connectivity level. In addition, lower scores on processing speed were significantly related to lower cortical thickness (ß = 0.20, Pcorrected = 0.030), lower R1 values (ß = 0.20, Pcorrected = 0.006), lower R2* values (ß = 0.29, Pcorrected = 0.006) and lower susceptibility values (ß = 0.19, Pcorrected = 0.024) of the WMH-connected regions at high connectivity level, independent of WMH volumes and the cortical measures of WMH-unconnected regions. Together, our study demonstrated that the microstructural integrity of white matter tracts passing through WMH is related to the regional cortical abnormalities as measured by thickness, R1, R2* and susceptibility values in the connected cortical regions. These findings are indicative of cortical thinning, demyelination and iron loss in the cortex, which is most likely through the disruption of the connecting white matter tracts and may contribute to processing speed impairment in SVD, a key clinical feature of SVD. These findings may have implications for finding intervention targets for the treatment of cognitive impairment in SVD by preventing secondary degeneration.

Dissociable Contributions of Thalamic-Subregions to Cognitive Impairment in Small Vessel Disease.

BACKGROUND: Structural network damage is a potentially important mechanism by which cerebral small vessel disease (SVD) can cause cognitive impairment. As a central hub of the structural network, the role of thalamus in SVD-related cognitive impairments remains unclear. We aimed to determine the associations between the structural alterations of thalamic subregions and cognitive impairments in SVD. METHODS: In this cross-sectional study, 205 SVD participants without thalamic lacunes from the third follow-up (2020) of the prospective RUN DMC study (Radboud University Nijmegen Diffusion Tensor and Magnetic Resonance Cohort), which was initiated in 2006, Nijmegen, were included. Cognitive functions included processing speed, executive function, and memory. Probabilistic tractography was performed from thalamus to 6 cortical regions, followed by connectivity-based thalamic segmentation to assess each thalamic subregion volume and connectivity (measured by mean diffusivity [MD] of the connecting white matter tracts) with the cortex. Least absolute shrinkage and selection operator regression analysis was conducted to identify the volumes or connectivity of the total thalamus and 6 thalamic subregions that have the strongest association with cognitive performance. Linear regression and mediation analyses were performed to test the association of least absolute shrinkage and selection operator-selected thalamic subregion volume or MD with cognitive performance, while adjusting for age and education. RESULTS: We found that higher MD of the thalamic-motor tract was associated with worse processing speed (ß=-0.27; P<0.001), higher MD of the thalamic-frontal tract was associated with worse executive function (ß=-0.24; P=0.001), and memory (ß=-0.28; P<0.001), respectively. The mediation analysis showed that MD of thalamocortical tracts mediated the association between corresponding thalamic subregion volumes and the cognitive performances in 3 domains. CONCLUSIONS: Our results suggest that the structural alterations of thalamus are linked to cognitive impairment in SVD, largely depending on the damage pattern of the white matter tracts connecting specific thalamic subregions and cortical regions.

A role for the fornix in temporal sequence memory.

Converging evidence from studies of human and nonhuman animals suggests that the hippocampus contributes to sequence learning by using temporal context to bind sequentially occurring items. The fornix is a white matter pathway containing the major input and output pathways of the hippocampus, including projections from medial septum and to diencephalon, striatum, lateral septum and prefrontal cortex. If the fornix meaningfully contributes to hippocampal function, then individual differences in fornix microstructure might predict sequence memory. Here, we tested this prediction by performing tractography in 51 healthy adults who had undertaken a sequence memory task. Microstructure properties of the fornix were compared with those of tracts connecting medial temporal lobe regions but not predominantly the hippocampus: the Parahippocampal Cingulum bundle (PHC) (conveying retrosplenial projections to parahippocampal cortex) and the Inferior Longitudinal Fasciculus (ILF) (conveying occipital projections to perirhinal cortex). Using principal components analysis, we combined Free-Water Elimination Diffusion Tensor Imaging and Neurite Orientation Dispersion and Density Imaging measures obtained from multi-shell diffusion MRI into two informative indices: the first (PC1) capturing axonal packing/myelin and the second (PC2) capturing microstructural complexity. We found a significant correlation between fornix PC2 and implicit reaction-time indices of sequence memory, indicating that greater fornix microstructural complexity is associated with better sequence memory. No such relationship was found with measures from the PHC and ILF. This study highlights the importance of the fornix in aiding memory for objects within a temporal context, potentially reflecting a role in mediating inter-regional communication within an extended hippocampal system.

On the performance of multi-compartment relaxometry for myelin water imaging (MCR-MWI) - test-retest repeatability and inter-protocol reproducibility.

In this study, we optimized the variable flip angle (VFA) acquisition scheme using numerical simulations to shorten the acquisition time of multicompartment relaxometry for myelin water imaging (MCR-MWI) to a clinically practical range in the absence of advanced image reconstruction methods. As the primary objective of this study, the test-retest repeatability of myelin water fraction (MWF) measurements of MCR-MWI is evaluated on three gradient echo (GRE) sequence settings using the optimized VFA schemes with different echo times and repetition times, emulating various scanner setups. The cross-protocol reproducibility of MCR-MWI and MCR with diffusion-informed myelin water imaging (MCR-DIMWI) is also examined. As a secondary objective, we explore the bundle-specific profiles of various microstructural parameters from MCR-(DI)MWI and their cross-correlations to determine if these parameters possess supplementary microstructure information beyond myelin concentration. Numerical simulations indicate that MCR-MWI can be performed with a minimum of three flip angles covering a wide range of T1 weightings without adding significant bias. This is supported by the results of an in vivo experiment, allowing whole-brain 1.5 mm isotropic MWF maps to be acquired in 9 min, reducing the total scan time to 40% of the original implementation without significant quality degradation. Good test-retest repeatability is observed for MCR-MWI for all three GRE protocols. While good correlations can also be found in MWF across protocols, systematic differences are observed. Bundle-specific MWF analysis reveals that certain white matter bundles are similar in all participants. We also found that microstructure relaxation parameters have low linear correlations with MWF. MCR-MWI is a reproducible measure of myelin. However, attention should be paid to the protocol related MWF differences when comparing different studies, as the MWF bias up to 0.5% can be observed across the protocols examined in this work.

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.

Tractography dissection variability: What happens when 42 groups dissect 14 white matter bundles on the same dataset?

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.

Tract-specific MRI measures explain learning and recall differences in multiple sclerosis.

Cognitive difficulties are common and a key concern for people with multiple sclerosis. Advancing knowledge of the role of white matter pathology in multiple sclerosis-related cognitive impairment is essential as both occur early in the disease with implications for early intervention. Consequently, this cross-sectional study asked whether quantifying the relationships between lesions and specific white matter structures could better explain co-existing cognitive differences than whole brain imaging measures. Forty participants with relapse-onset multiple sclerosis underwent cognitive testing and MRI at 3 Tesla. They were classified as cognitively impaired (n = 24) or unimpaired (n = 16) and differed across verbal fluency, learning and recall tasks corrected for intelligence and education (corrected P-values = 0.007-0.04). The relationships between lesions and white matter were characterized across six measures: conventional voxel-based T2 lesion load, whole brain tractogram load (lesioned volume/whole tractogram volume), whole bundle volume, bundle load (lesioned volume/whole bundle volume), Tractometry (diffusion-tensor and high angular resolution diffusion measures sampled from all bundle streamlines) and lesionometry (diffusion measures sampled from streamlines traversing lesions only). The tract-specific measures were extracted from corpus callosum segments (genu and isthmus), striato-prefrontal and -parietal pathways, and the superior longitudinal fasciculi (sections I, II and III). White matter measure-task associations demonstrating at least moderate evidence against the null hypothesis (Bayes Factor threshold < 0.2) were examined using independent t-tests and covariate analyses (significance level P < 0.05). Tract-specific measures were significant predictors (all P-values < 0.05) of task-specific clinical scores and diminished the significant effect of group as a categorical predictor in Story Recall (isthmus bundle load), Figure Recall (right striato-parietal lesionometry) and Design Learning (left superior longitudinal fasciculus III volume). Lesion load explained the difference in List Learning, whereas Letter Fluency was not associated with any of the imaging measures. Overall, tract-specific measures outperformed the global lesion and tractogram load measures. Variation in regional lesion burden translated to group differences in tract-specific measures, which in turn, attenuated differences in individual cognitive tasks. The structural differences converged in temporo-parietal regions with particular influence on tasks requiring visuospatial-constructional processing. We highlight that measures quantifying the relationships between tract-specific structure and multiple sclerosis lesions uncovered associations with cognition masked by overall tract volumes and global lesion and tractogram loads. These tract-specific white matter quantifications show promise for elucidating the relationships between neuropathology and cognition in multiple sclerosis.

Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI.

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

Detecting microstructural deviations in individuals with deep diffusion MRI tractometry.

Most diffusion magnetic resonance imaging studies of disease rely on statistical comparisons between large groups of patients and healthy participants to infer altered tissue states in the brain; however, clinical heterogeneity can greatly challenge their discriminative power. There is currently an unmet need to move away from the current approach of group-wise comparisons to methods with the sensitivity to detect altered tissue states at the individual level. This would ultimately enable the early detection and interpretation of microstructural abnormalities in individual patients, an important step towards personalized medicine in translational imaging. To this end, Detect was developed to advance diffusion magnetic resonance imaging tractometry towards single-patient analysis. By operating on the manifold of white-matter pathways and learning normative microstructural features, our framework captures idiosyncrasies in patterns along white-matter pathways. Our approach paves the way from traditional group-based comparisons to true personalized radiology, taking microstructural imaging from the bench to the bedside.

Resolving bundle-specific intra-axonal T2 values within a voxel using diffusion-relaxation tract-based estimation.

At the typical spatial resolution of MRI in the human brain, approximately 60-90% of voxels contain multiple fiber populations. Quantifying microstructural properties of distinct fiber populations within a voxel is therefore challenging but necessary. While progress has been made for diffusion and T1-relaxation properties, how to resolve intra-voxel T2 heterogeneity remains an open question. Here a novel framework, named COMMIT-T2, is proposed that uses tractography-based spatial regularization with diffusion-relaxometry data to estimate multiple intra-axonal T2 values within a voxel. Unlike previously-proposed voxel-based T2 estimation methods, which (when applied in white matter) implicitly assume just one fiber bundle in the voxel or the same T2 for all bundles in the voxel, COMMIT-T2 can recover specific T2 values for each unique fiber population passing through the voxel. In this approach, the number of recovered unique T2 values is not determined by a number of model parameters set a priori, but rather by the number of tractography-reconstructed streamlines passing through the voxel. Proof-of-concept is provided in silico and in vivo, including a demonstration that distinct tract-specific T2 profiles can be recovered even in the three-way crossing of the corpus callosum, arcuate fasciculus, and corticospinal tract. We demonstrate the favourable performance of COMMIT-T2 compared to that of voxelwise approaches for mapping intra-axonal T2 exploiting diffusion, including a direction-averaged method and AMICO-T2, a new extension to the previously-proposed Accelerated Microstructure Imaging via Convex Optimization (AMICO) framework.

Computing and visualising intra-voxel orientation-specific relaxation-diffusion features in the human brain.

Diffusion MRI techniques are used widely to study the characteristics of the human brain connectome in vivo. However, to resolve and characterise white matter (WM) fibres in heterogeneous MRI voxels remains a challenging problem typically approached with signal models that rely on prior information and constraints. We have recently introduced a 5D relaxation-diffusion correlation framework wherein multidimensional diffusion encoding strategies are used to acquire data at multiple echo-times to increase the amount of information encoded into the signal and ease the constraints needed for signal inversion. Nonparametric Monte Carlo inversion of the resulting datasets yields 5D relaxation-diffusion distributions where contributions from different sub-voxel tissue environments are separated with minimal assumptions on their microscopic properties. Here, we build on the 5D correlation approach to derive fibre-specific metrics that can be mapped throughout the imaged brain volume. Distribution components ascribed to fibrous tissues are resolved, and subsequently mapped to a dense mesh of overlapping orientation bins to define a smooth orientation distribution function (ODF). Moreover, relaxation and diffusion measures are correlated to each independent ODF coordinate, thereby allowing the estimation of orientation-specific relaxation rates and diffusivities. The proposed method is tested on a healthy volunteer, where the estimated ODFs were observed to capture major WM tracts, resolve fibre crossings, and, more importantly, inform on the relaxation and diffusion features along with distinct fibre bundles. If combined with fibre-tracking algorithms, the methodology presented in this work has potential for increasing the depth of characterisation of microstructural properties along individual WM pathways.

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.

Multimodal principal component analysis to identify major features of white matter structure and links to reading.

The role of white matter in reading has been established by diffusion tensor imaging (DTI), but DTI cannot identify specific microstructural features driving these relationships. Neurite orientation dispersion and density imaging (NODDI), inhomogeneous magnetization transfer (ihMT) and multicomponent driven equilibrium single-pulse observation of T1/T2 (mcDESPOT) can be used to link more specific aspects of white matter microstructure and reading due to their sensitivity to axonal packing and fiber coherence (NODDI) and myelin (ihMT and mcDESPOT). We applied principal component analysis (PCA) to combine DTI, NODDI, ihMT and mcDESPOT measures (10 in total), identify major features of white matter structure, and link these features to both reading and age. Analysis was performed for nine reading-related tracts in 46 neurotypical 6-16 year olds. We identified three principal components (PCs) which explained 79.5% of variance in our dataset. PC1 probed tissue complexity, PC2 described myelin and axonal packing, while PC3 was related to axonal diameter. Mixed effects regression models did not identify any significant relationships between principal components and reading skill. Bayes factor analysis revealed that the absence of relationships was not due to low power. Increasing PC1 in the left arcuate fasciculus with age suggest increases in tissue complexity, while increases of PC2 in the bilateral arcuate, inferior longitudinal, inferior fronto-occipital fasciculi, and splenium suggest increases in myelin and axonal packing with age. Multimodal white matter imaging and PCA provide microstructurally informative, powerful principal components which can be used by future studies of development and cognition. Our findings suggest major features of white matter undergo development during childhood and adolescence, but changes are not linked to reading during this period in our typically-developing sample.

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.

Tractostorm: The what, why, and how of tractography dissection reproducibility.

Investigative studies of white matter (WM) brain structures using diffusion MRI (dMRI) tractography frequently require manual WM bundle segmentation, often called "virtual dissection." Human errors and personal decisions make these manual segmentations hard to reproduce, which have not yet been quantified by the dMRI community. It is our opinion that if the field of dMRI tractography wants to be taken seriously as a widespread clinical tool, it is imperative to harmonize WM bundle segmentations and develop protocols aimed to be used in clinical settings. The EADC-ADNI Harmonized Hippocampal Protocol achieved such standardization through a series of steps that must be reproduced for every WM bundle. This article is an observation of the problematic. A specific bundle segmentation protocol was used in order to provide a real-life example, but the contribution of this article is to discuss the need for reproducibility and standardized protocol, as for any measurement tool. This study required the participation of 11 experts and 13 nonexperts in neuroanatomy and "virtual dissection" across various laboratories and hospitals. Intra-rater agreement (Dice score) was approximately 0.77, while inter-rater was approximately 0.65. The protocol provided to participants was not necessarily optimal, but its design mimics, in essence, what will be required in future protocols. Reporting tractometry results such as average fractional anisotropy, volume or streamline count of a particular bundle without a sufficient reproducibility score could make the analysis and interpretations more difficult. Coordinated efforts by the diffusion MRI tractography community are needed to quantify and account for reproducibility of WM bundle extraction protocols in this era of open and collaborative science.

Author Correction: The challenge of mapping the human connectome based on diffusion tractography.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Reducing variability in along-tract analysis with diffusion profile realignment.

Diffusion weighted magnetic resonance imaging (dMRI) provides a non invasive virtual reconstruction of the brain's white matter structures through tractography. Analyzing dMRI measures along the trajectory of white matter bundles can provide a more specific investigation than considering a region of interest or tract-averaged measurements. However, performing group analyses with this along-tract strategy requires correspondence between points of tract pathways across subjects. This is usually achieved by creating a new common space where the representative streamlines from every subject are resampled to the same number of points. If the underlying anatomy of some subjects was altered due to, e.g., disease or developmental changes, such information might be lost by resampling to a fixed number of points. In this work, we propose to address the issue of possible misalignment, which might be present even after resampling, by realigning the representative streamline of each subject in this 1D space with a new method, coined diffusion profile realignment (DPR). Experiments on synthetic datasets show that DPR reduces the coefficient of variation for the mean diffusivity, fractional anisotropy and apparent fiber density when compared to the unaligned case. Using 100 in vivo datasets from the human connectome project, we simulated changes in mean diffusivity, fractional anisotropy and apparent fiber density. Independent Student's t-tests between these altered subjects and the original subjects indicate that regional changes are identified after realignment with the DPR algorithm, while preserving differences previously detected in the unaligned case. This new correction strategy contributes to revealing effects of interest which might be hidden by misalignment and has the potential to improve the specificity in longitudinal population studies beyond the traditional region of interest based analysis and along-tract analysis workflows.

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.