Characterization of the orientation dependence of magnetization transfer measures in single and crossing-fiber white matter.

PURPOSE: To fully characterize the orientation dependence of magnetization transfer (MT) and inhomogeneous MT (ihMT) measures in the whole white matter (WM), for both single-fiber and crossing-fiber voxels. METHODS: A characterization method was developed using the fiber orientation obtained from diffusion MRI (dMRI) with diffusion tensor imaging (DTI) and constrained spherical deconvolution. This allowed for characterization of the orientation dependence of measures in all of WM, regardless of the number of fiber orientation in a voxel. Furthermore, the orientation dependence inside 31 different WM bundles was characterized to evaluate the homogeneity of the effect. Variation of the results within and between-subject was assessed from a 12-subject dataset. RESULTS: Previous results for single-fiber voxels were reproduced and a novel characterization was produced in voxels of crossing fibers, which seems to follow trends consistent with single-fiber results. Heterogeneity of the orientation dependence across bundles was observed, but homogeneity within similar bundles was also highlighted. Differences in behavior between MT and ihMT measures, as well as the ratio and saturation versions of these, were noted. CONCLUSION: Orientation dependence characterization was proven possible over the entirety of WM. The vast range of effects and subtleties of the orientation dependence on MT measures showed the need for, but also the challenges of, a correction method.

Pore size estimation in axon-mimicking microfibers with diffusion-relaxation MRI.

PURPOSE: This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion-relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power-law approach and a T2-based pore size estimation technique. THEORY AND METHODS: A general diffusion-relaxation theoretical model for the spherical mean signal from water molecules within a distribution of cylinders with varying radii was introduced, encompassing the evaluated models as particular cases. Additionally, a new numerical approach was presented for estimating effective radii (i.e., MRI-visible mean radii) from the ground truth radii distributions, not reliant on previous theoretical approximations and adaptable to various acquisition sequences. The ground truth radii were obtained from scanning electron microscope images. RESULTS: Both methods show a linear relationship between effective radii estimated from MRI data and ground-truth radii distributions, although some discrepancies were observed. The spherical mean power-law method overestimated fiber radii. Conversely, the T2-based method exhibited higher sensitivity to smaller fiber radii, but faced limitations in accurately estimating the radius in one particular phantom, possibly because of material-specific relaxation changes. CONCLUSION: The study demonstrates the feasibility of both techniques to predict pore sizes of hollow microfibers. The T2-based technique, unlike the spherical mean power-law method, does not demand ultra-high diffusion gradients, but requires calibration with known radius distributions. This research contributes to the ongoing development and evaluation of neuroimaging techniques for fiber radius estimation, highlights the advantages and limitations of both methods, and provides datasets for reproducible research.

FIESTA: Autoencoders for accurate fiber segmentation in tractography.

White matter bundle segmentation is a cornerstone of modern tractography to study the brain's structural connectivity in domains such as neurological disorders, neurosurgery, and aging. In this study, we present FIESTA (FIbEr Segmentation in Tractography using Autoencoders), a reliable and robust, fully automated, and easily semi-automatically calibrated pipeline based on deep autoencoders that can dissect and fully populate white matter bundles. This pipeline is built upon previous works that demonstrated how autoencoders can be used successfully for streamline filtering, bundle segmentation, and streamline generation in tractography. Our proposed method improves bundle segmentation coverage by recovering hard-to-track bundles with generative sampling through the latent space seeding of the subject bundle and the atlas bundle. A latent space of streamlines is learned using autoencoder-based modeling combined with contrastive learning. Using an atlas of bundles in standard space (MNI), our proposed method segments new tractograms using the autoencoder latent distance between each tractogram streamline and its closest neighbor bundle in the atlas of bundles. Intra-subject bundle reliability is improved by recovering hard-to-track streamlines, using the autoencoder to generate new streamlines that increase the spatial coverage of each bundle while remaining anatomically correct. Results show that our method is more reliable than state-of-the-art automated virtual dissection methods such as RecoBundles, RecoBundlesX, TractSeg, White Matter Analysis and XTRACT. Our framework allows for the transition from one anatomical bundle definition to another with marginal calibration efforts. Overall, these results show that our framework improves the practicality and usability of current state-of-the-art bundle segmentation framework.

High-frequency longitudinal white matter diffusion- and myelin-based MRI database: Reliability and variability.

Assessing the consistency of quantitative MRI measurements is critical for inclusion in longitudinal studies and clinical trials. Intraclass coefficient correlation and coefficient of variation were used to evaluate the different consistency aspects of diffusion- and myelin-based MRI measures. Multi-shell diffusion and inhomogeneous magnetization transfer data sets were collected from 20 healthy adults at a high-frequency of five MRI sessions. The consistency was evaluated across whole bundles and the track-profile along the bundles. The impact of the fiber populations on the consistency was also evaluated using the number of fiber orientations map. For whole and profile bundles, moderate to high reliability of diffusion and myelin measures were observed. We report higher reliability of measures for multiple fiber populations than single. The overall portrait of the most consistent measurements and bundles drawn from a wide range of MRI techniques presented here will be particularly useful for identifying reliable biomarkers capable of detecting, monitoring and predicting white matter changes in clinical applications and has the potential to inform patient-specific treatment strategies.

A New Advanced MRI Biomarker for Remyelinated Lesions in Multiple Sclerosis.

OBJECTIVES: Neuropathological studies have shown that multiple sclerosis (MS) lesions are heterogeneous in terms of myelin/axon damage and repair as well as iron content. However, it remains a challenge to identify specific chronic lesion types, especially remyelinated lesions, in vivo in patients with MS. METHODS: We performed 3 studies: (1) a cross-sectional study in a prospective cohort of 115 patients with MS and 76 healthy controls, who underwent 3 T magnetic resonance imaging (MRI) for quantitative susceptibility mapping (QSM), myelin water fraction (MWF), and neurite density index (NDI) maps. White matter (WM) lesions in QSM were classified into 5 QSM lesion types (iso-intense, hypo-intense, hyperintense, lesions with hypo-intense rims, and lesions with paramagnetic rim legions [PRLs]); (2) a longitudinal study of 40 patients with MS to study the evolution of lesions over 2 years; (3) a postmortem histopathology-QSM validation study in 3 brains of patients with MS to assess the accuracy of QSM classification to identify neuropathological lesion types in 63 WM lesions. RESULTS: At baseline, hypo- and isointense lesions showed higher mean MWF and NDI values compared to other QSM lesion types (p < 0.0001). Further, at 2-year follow-up, hypo-/iso-intense lesions showed an increase in MWF. Postmortem analyses revealed that QSM highly accurately identifies (1) fully remyelinated areas as hypo-/iso-intense (sensitivity = 88.89% and specificity = 100%), (2) chronic inactive lesions as hyperintense (sensitivity = 71.43% and specificity = 92.00%), and (3) chronic active/smoldering lesions as PRLs (sensitivity = 92.86% and specificity = 86.36%). INTERPRETATION: These results provide the first evidence that it is possible to distinguish chronic MS lesions in a clinical setting, hereby supporting with new biomarkers to develop and assess remyelinating treatments. ANN NEUROL 2022;92:486-502.

A Multicenter Longitudinal MRI Study Assessing LeMan-PV Software Accuracy in the Detection of White Matter Lesions in Multiple Sclerosis Patients.

BACKGROUND: Detecting new and enlarged lesions in multiple sclerosis (MS) patients is needed to determine their disease activity. LeMan-PV is a software embedded in the scanner reconstruction system of one vendor, which automatically assesses new and enlarged white matter lesions (NELs) in the follow-up of MS patients; however, multicenter validation studies are lacking. PURPOSE: To assess the accuracy of LeMan-PV for the longitudinal detection NEL white-matter MS lesions in a multicenter clinical setting. STUDY TYPE: Retrospective, longitudinal. SUBJECTS: A total of 206 patients with a definitive MS diagnosis and at least two follow-up MRI studies from five centers participating in the Swiss Multiple Sclerosis Cohort study. Mean age at first follow-up = 45.2 years (range: 36.9-52.8 years); 70 males. FIELD STRENGTH/SEQUENCE: Fluid attenuated inversion recovery (FLAIR) and T1-weighted magnetization prepared rapid gradient echo (T1-MPRAGE) sequences at 1.5 T and 3 T. ASSESSMENT: The study included 313 MRI pairs of datasets. Data were analyzed with LeMan-PV and compared with a manual "reference standard" provided by a neuroradiologist. A second rater (neurologist) performed the same analysis in a subset of MRI pairs to evaluate the rating-accuracy. The Sensitivity (Se), Specificity (Sp), Accuracy (Acc), F1-score, lesion-wise False-Positive-Rate (aFPR), and other measures were used to assess LeMan-PV performance for the detection of NEL at 1.5 T and 3 T. The performance was also evaluated in the subgroup of 123 MRI pairs at 3 T. STATISTICAL TESTS: Intraclass correlation coefficient (ICC) and Cohen's kappa (CK) were used to evaluate the agreement between readers. RESULTS: The interreader agreement was high for detecting new lesions (ICC = 0.97, Pvalue < 10-20 , CK = 0.82, P value = 0) and good (ICC = 0.75, P value < 10-12 , CK = 0.68, P value = 0) for detecting enlarged lesions. Across all centers, scanner field strengths (1.5 T, 3 T), and for NEL, LeMan-PV achieved: Acc = 61%, Se = 65%, Sp = 60%, F1-score = 0.44, aFPR = 1.31. When both follow-ups were acquired at 3 T, LeMan-PV accuracy was higher (Acc = 66%, Se = 66%, Sp = 66%, F1-score = 0.28, aFPR = 3.03). DATA CONCLUSION: In this multicenter study using clinical data settings acquired at 1.5 T and 3 T, and variations in MRI protocols, LeMan-PV showed similar sensitivity in detecting NEL with respect to other recent 3 T multicentric studies based on neural networks. While LeMan-PV performance is not optimal, its main advantage is that it provides automated clinical decision support integrated into the radiological-routine flow. EVIDENCE LEVEL: 4 TECHNICAL EFFICACY: Stage 2.

Anterior horn atrophy in the cervical spinal cord: A new biomarker in progressive multiple sclerosis.

BACKGROUND: Spinal cord (SC) gray and white matter pathology plays a central role in multiple sclerosis (MS). OBJECTIVE: We aimed to investigate the extent, pattern, and clinical relevance of SC gray and white matter atrophy in vivo. METHODS: 39 relapsing-remitting patients (RRMS), 40 progressive MS patients (PMS), and 24 healthy controls (HC) were imaged at 3T using the averaged magnetization inversion recovery acquisitions sequence. Total and lesional cervical gray and white matter, and posterior (SCPH) and anterior horn (SCAH) areas were automatically quantified. Clinical assessment included the expanded disability status scale, timed 25-foot walk test, nine-hole peg test, and the 12-item MS walking scale. RESULTS: PMS patients had significantly reduced cervical SCAH - but not SCPH - areas compared with HC and RRMS (both p < 0.001). In RRMS and PMS, the cervical SCAH areas increased significantly less in the region of cervical SC enlargement compared with HC (all p < 0.001). This reduction was more pronounced in PMS compared with RRMS (both p < 0.001). In PMS, a lower cervical SCAH area was the most important magnetic resonance imaging (MRI)-variable for higher disability scores. CONCLUSION: MS patients show clinically relevant cervical SCAH atrophy, which is more pronounced in PMS and at the level of cervical SC enlargement.

Optical coherence tomography versus other biomarkers: Associations with physical and cognitive disability in multiple sclerosis.

BACKGROUND: Optical coherence tomography (OCT) is a biomarker of neuroaxonal loss in multiple sclerosis (MS). OBJECTIVE: The objective was to assess the relative role of OCT, next to magnetic resonance imaging (MRI) and serum markers of disability in MS. METHODS: A total of 100 patients and 52 controls underwent OCT to determine peripapillary retinal nerve fiber layer (pRNFL) and ganglion cell-inner plexiform layers (GCIPL). Serum neurofilament light chain (sNfL), total lesion volume (TLV), and brain parenchymal fraction (BPF) were also assessed. The associations of OCT with disability were examined in linear regression models with correction for age, vision, and education. RESULTS: In patients, pRNFL was associated with the Symbol Digit Modalities Test (SDMT; p = 0.030). In the multivariate analysis including sNfL and MRI measures, pRNFL (ß = 0.19, p = 0.044) and TLV (ß = -0.24, p = 0.023) were the only markers associated with the SDMT. pRNFL (p < 0.001) and GCIPL (p < 0.001) showed associations with the Expanded Disability Status Scale (EDSS). In the multivariate analysis, GCIPL showed the strongest association with the EDSS (ß = -0.32, p < 0.001) followed by sNfL (ß = 0.18, p = 0.024). CONCLUSION: The associations of OCT measures with cognitive and physical disability were independent of serum and brain MRI markers of neuroaxonal loss. OCT can be an important tool for stratification in MS, while longitudinal studies using combinations of biomarkers are warranted.

Insights from the IronTract challenge: Optimal methods for mapping brain pathways from multi-shell diffusion MRI.

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Intrathecal Immunoglobulin M Synthesis is an Independent Biomarker for Higher Disease Activity and Severity in Multiple Sclerosis.

OBJECTIVE: We aimed to determine in relapsing multiple sclerosis (MS) whether intrathecal synthesis of immunoglobulin (Ig) M and IgG is associated with outcomes reflecting inflammatory activity and chronic worsening. METHODS: We compared cerebrospinal fluid analysis, clinical and magnetic resonance imaging data, and serum neurofilament light chain (sNfL) levels at baseline and follow-up in 530 patients with relapsing MS. Patients were categorized by the presence of oligoclonal IgG bands (OCGB) and intrathecal synthesis of IgG and IgM (intrathecal fraction [IF]: IgGIF and IgMIF ). Relationships with the time to first relapse, sNfL concentrations, T2-weighted (T2w) lesions, MS Severity Score (MSSS), and time to initiation of high-efficacy therapy were analyzed in covariate-adjusted statistical models. RESULTS: By categorical analysis, in patients with IgMIF the median time to first relapse was 28 months shorter and MSSS on average higher by 1.11 steps compared with patients without intrathecal immunoglobulin synthesis. Moreover, patients with IgMIF had higher sNfL concentrations, more new/enlarging T2w lesions, and higher total T2w lesion counts (all p ≤ 0.01). These associations were absent or equally smaller in patients who were positive for only OCGB or OCGB/IgGIF . Furthermore, quantitative analyses revealed that in patients with IgMIF ≥ median, the time to first relapse and to initiation of high-efficacy therapy was shorter by 32 and by 203 months, respectively (both p < 0.01), in comparison to patients with IgMIF < median. Dose-dependent associations were also found for IgMIF but not for IgGIF with magnetic resonance imaging-defined disease activity and sNfL. INTERPRETATION: This large study supports the value of intrathecal IgM synthesis as an independent biomarker of disease activity and severity in relapsing MS. ANN NEUROL 2021;90:477-489.

Myelin and axon pathology in multiple sclerosis assessed by myelin water and multi-shell diffusion imaging.

Damage to the myelin sheath and the neuroaxonal unit is a cardinal feature of multiple sclerosis; however, a detailed characterization of the interaction between myelin and axon damage in vivo remains challenging. We applied myelin water and multi-shell diffusion imaging to quantify the relative damage to myelin and axons (i) among different lesion types; (ii) in normal-appearing tissue; and (iii) across multiple sclerosis clinical subtypes and healthy controls. We also assessed the relation of focal myelin/axon damage with disability and serum neurofilament light chain as a global biological measure of neuroaxonal damage. Ninety-one multiple sclerosis patients (62 relapsing-remitting, 29 progressive) and 72 healthy controls were enrolled in the study. Differences in myelin water fraction and neurite density index were substantial when lesions were compared to healthy control subjects and normal-appearing multiple sclerosis tissue: both white matter and cortical lesions exhibited a decreased myelin water fraction and neurite density index compared with healthy (P < 0.0001) and peri-plaque white matter (P < 0.0001). Periventricular lesions showed decreased myelin water fraction and neurite density index compared with lesions in the juxtacortical region (P < 0.0001 and P < 0.05). Similarly, lesions with paramagnetic rims showed decreased myelin water fraction and neurite density index relative to lesions without a rim (P < 0.0001). Also, in 75% of white matter lesions, the reduction in neurite density index was higher than the reduction in the myelin water fraction. Besides, normal-appearing white and grey matter revealed diffuse reduction of myelin water fraction and neurite density index in multiple sclerosis compared to healthy controls (P < 0.01). Further, a more extensive reduction in myelin water fraction and neurite density index in normal-appearing cortex was observed in progressive versus relapsing-remitting participants. Neurite density index in white matter lesions correlated with disability in patients with clinical deficits (P < 0.01, beta = -10.00); and neurite density index and myelin water fraction in white matter lesions were associated to serum neurofilament light chain in the entire patient cohort (P < 0.01, beta = -3.60 and P < 0.01, beta = 0.13, respectively). These findings suggest that (i) myelin and axon pathology in multiple sclerosis is extensive in both lesions and normal-appearing tissue; (ii) particular types of lesions exhibit more damage to myelin and axons than others; (iii) progressive patients differ from relapsing-remitting patients because of more extensive axon/myelin damage in the cortex; and (iv) myelin and axon pathology in lesions is related to disability in patients with clinical deficits and global measures of neuroaxonal damage.

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.

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.

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.

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.

ActiveAxADD : Toward non-parametric and orientationally invariant axon diameter distribution mapping using PGSE.

PURPOSE: Non-invasive axon diameter distribution (ADD) mapping using diffusion MRI is an ill-posed problem. Current ADD mapping methods require knowledge of axon orientation before performing the acquisition. Instead, ActiveAx uses a 3D sampling scheme to estimate the orientation from the signal, providing orientationally invariant estimates. The mean diameter is estimated instead of the distribution for the solution to be tractable. Here, we propose an extension (ActiveAxADD ) that provides non-parametric and orientationally invariant estimates of the whole distribution. THEORY: The accelerated microstructure imaging with convex optimization (AMICO) framework accelerates mean diameter estimation using a linear formulation combined with Tikhonov regularization to stabilize the solution. Here, we implement a new formulation (ActiveAxADD ) that uses Laplacian regularization to provide robust estimates of the whole ADD. METHODS: The performance of ActiveAxADD was evaluated using Monte Carlo simulations on synthetic white matter samples mimicking axon distributions reported in histological studies. RESULTS: ActiveAxADD provided robust ADD reconstructions when considering the isolated intra-axonal signal. However, our formulation inherited some common microstructure imaging limitations. When accounting for the extra axonal compartment, estimated ADDs showed spurious peaks and increased variability because of the difficulty of disentangling intra and extra axonal contributions. CONCLUSION: Laplacian regularization solves the ill-posedness regarding the intra axonal compartment. ActiveAxADD can potentially provide non-parametric and orientationally invariant ADDs from isolated intra-axonal signals. However, further work is required before ActiveAxADD can be applied to real data containing extra-axonal contributions, as disentangling the 2 compartment appears to be an overlooked challenge that affects microstructure imaging methods in general.

Tractography reproducibility challenge with empirical data (TraCED): The 2017 ISMRM diffusion study group challenge.

BACKGROUND: Fiber tracking with diffusion-weighted MRI has become an essential tool for estimating in vivo brain white matter architecture. Fiber tracking results are sensitive to the choice of processing method and tracking criteria. PURPOSE: To assess the variability for an algorithm in group studies reproducibility is of critical context. However, reproducibility does not assess the validity of the brain connections. Phantom studies provide concrete quantitative comparisons of methods relative to absolute ground truths, yet do no capture variabilities because of in vivo physiological factors. The ISMRM 2017 TraCED challenge was created to fulfill the gap. STUDY TYPE: A systematic review of algorithms and tract reproducibility studies. SUBJECTS: Single healthy volunteers. FIELD STRENGTH/SEQUENCE: 3.0T, two different scanners by the same manufacturer. The multishell acquisition included b-values of 1000, 2000, and 3000 s/mm2 with 20, 45, and 64 diffusion gradient directions per shell, respectively. ASSESSMENT: Nine international groups submitted 46 tractography algorithm entries each consisting 16 tracts per scan. The algorithms were assessed using intraclass correlation (ICC) and the Dice similarity measure. STATISTICAL TESTS: Containment analysis was performed to assess if the submitted algorithms had containment within tracts of larger volume submissions. This also serves the purpose to detect if spurious submissions had been made. RESULTS: The top five submissions had high ICC and Dice >0.88. Reproducibility was high within the top five submissions when assessed across sessions or across scanners: 0.87-0.97. Containment analysis shows that the top five submissions are contained within larger volume submissions. From the total of 16 tracts as an outcome relatively the number of tracts with high, moderate, and low reproducibility were 8, 4, and 4. DATA CONCLUSION: The different methods clearly result in fundamentally different tract structures at the more conservative specificity choices. Data and challenge infrastructure remain available for continued analysis and provide a platform for comparison. LEVEL OF EVIDENCE: 5 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:234-249.

Sparse wars: A survey and comparative study of spherical deconvolution algorithms for diffusion MRI.

Spherical deconvolution methods are widely used to estimate the brain's white-matter fiber orientations from diffusion MRI data. In this study, eight spherical deconvolution algorithms were implemented and evaluated. These included two model selection techniques based on the extended Bayesian information criterion (i.e., best subset selection and the least absolute shrinkage and selection operator), iteratively reweighted l2- and l1-norm approaches to approximate the l0-norm, sparse Bayesian learning, Cauchy deconvolution, and two accelerated Richardson-Lucy algorithms. Results from our exhaustive evaluation show that there is no single optimal method for all different fiber configurations, suggesting that further studies should be conducted to find the optimal way of combining solutions from different methods. We found l0-norm regularization algorithms to resolve more accurately fiber crossings with small inter-fiber angles. However, in voxels with very dominant fibers, algorithms promoting more sparsity are less accurate in detecting smaller fibers. In most cases, the best algorithm to reconstruct fiber crossings with two fibers did not perform optimally in voxels with one or three fibers. Therefore, simplified validation systems as employed in a number of previous studies, where only two fibers with similar volume fractions were tested, should be avoided as they provide incomplete information. Future studies proposing new reconstruction methods based on high angular resolution diffusion imaging data should validate their results by considering, at least, voxels with one, two, and three fibers, as well as voxels with dominant fibers and different diffusion anisotropies.

Limits to anatomical accuracy of diffusion tractography using modern approaches.

Diffusion MRI fiber tractography is widely used to probe the structural connectivity of the brain, with a range of applications in both clinical and basic neuroscience. Despite widespread use, tractography has well-known pitfalls that limits the anatomical accuracy of this technique. Numerous modern methods have been developed to address these shortcomings through advances in acquisition, modeling, and computation. To test whether these advances improve tractography accuracy, we organized the 3-D Validation of Tractography with Experimental MRI (3D-VoTEM) challenge at the ISBI 2018 conference. We made available three unique independent tractography validation datasets - a physical phantom and two ex vivo brain specimens - resulting in 176 distinct submissions from 9 research groups. By comparing results over a wide range of fiber complexities and algorithmic strategies, this challenge provides a more comprehensive assessment of tractography's inherent limitations than has been reported previously. The central results were consistent across all sub-challenges in that, despite advances in tractography methods, the anatomical accuracy of tractography has not dramatically improved in recent years. Taken together, our results independently confirm findings from decades of tractography validation studies, demonstrate inherent limitations in reconstructing white matter pathways using diffusion MRI data alone, and highlight the need for alternative or combinatorial strategies to accurately map the fiber pathways of the brain.

Down-sampling in diffusion MRI: a bundle-specific DTI and NODDI study.

Introduction: Multi-shell diffusion Magnetic Resonance Imaging (dMRI) data has been widely used to characterise white matter microstructure in several neurodegenerative diseases. The lack of standardised dMRI protocols often implies the acquisition of redundant measurements, resulting in prolonged acquisition times. In this study, we investigate the impact of the number of gradient directions on Diffusion Tensor Imaging (DTI) and on Neurite Orientation Dispersion and Density Imaging (NODDI) metrics. Methods: Data from 124 healthy controls collected in three different longitudinal studies were included. Using an in-house algorithm, we reduced the number of gradient directions in each data shell. We estimated DTI and NODDI measures on six white matter bundles clinically relevant for neurodegenerative diseases. Results: Fractional Anisotropy (FA) measures on bundles where data were sampled at the 30% rate, showed a median L1 distance of up to 3.92% and a 95% CI of (1.74, 8.97)% when compared to those obtained at reference sampling. Mean Diffusivity (MD) reached up to 4.31% and a 95% CI of (1.60, 16.98)% on the same premises. At a sampling rate of 50%, we obtained a median of 3.90% and a 95% CI of (1.99, 16.65)% in FA, and 5.49% with a 95% CI of (2.14, 21.68)% in MD. The Intra-Cellular volume fraction (ICvf) median L1 distance was up to 2.83% with a 95% CI of (1.98, 4.82)% at a 30% sampling rate and 3.95% with a 95% CI of (2.39, 7.81)% at a 50% sampling rate. The volume difference of the reconstructed white matter at reference and 50% sampling reached a maximum of (2.09 ± 0.81)%. Discussion: In conclusion, DTI and NODDI measures reported at reference sampling were comparable to those obtained when the number of dMRI volumes was reduced by up to 30%. Close to reference DTI and NODDI metrics were estimated with a significant reduction in acquisition time using three shells, respectively with: 4 directions at a b value of 700 s/mm2, 14 at 1000 s/mm2, and 32 at 2000 s/mm2. The study revealed aspects that can be important for large-scale clinical studies on bundle-specific diffusion MRI.