Bundle-o-graphy: improving structural connectivity estimation with adaptive microstructure-informed tractography

Tractography is a powerful tool for the investigation of the complex organization of the brain in vivo, as it allows inferring the macroscopic pathways of the major fiber bundles of the white matter based on non-invasive diffusion-weighted magnetic resonance imaging acquisitions. Despite this unique and compelling ability, some studies have exposed the poor anatomical accuracy of the reconstructions obtained with this technique and challenged its effectiveness for studying brain connectivity. In this work, we describe a novel method to readdress tractography reconstruction problem in a global manner by combining the strengths of so-called generative and discriminative strategies. Starting from an input tractogram, we parameterize the connections between brain regions following a bundle-based representation that allows to drastically reducing the number of parameters needed to model groups of fascicles. The parameters space is explored following an MCMC generative approach, while a discrimininative method is exploited to globally evaluate the set of connections which is updated according to Bayes' rule. Our results on both synthetic and real brain data show that the proposed solution, called bundle-o-graphy, allows improving the anatomical accuracy of the reconstructions while keeping the computational complexity similar to other state-of-the-art methods

Incorporating outlier information into diffusion-weighted MRI modeling for robust microstructural imaging and structural brain connectivity analyses

A B S T R A C T The white matter structures of the human brain can be represented using diffusion-weighted MRI tractography. Unfortunately, tractography is prone to find false-positive streamlines causing a severe decline in its specificity and limiting its feasibility in accurate structural brain connectivity analyses. Filtering algorithms have been pro-posed to reduce the number of invalid streamlines but the currently available filtering algorithms are not suitable to process data that contains motion artefacts which are typical in clinical research. We augmented the Con-vex Optimization Modelling for Microstructure Informed Tractography (COMMIT) algorithm to adjust for these signals drop-out motion artefacts. We demonstrate with comprehensive Monte-Carlo whole brain simulations and in vivo infant data that our robust algorithm is capable of properly filtering tractography reconstructions despite these artefacts. We evaluated the results using parametric and non-parametric statistics and our results demonstrate that if not accounted for, motion artefacts can have severe adverse effects in human brain structural connectivity analyses as well as in microstructural property mappings. In conclusion, the usage of robust filtering methods to mitigate motion related errors in tractogram filtering is highly beneficial, especially in clinical stud-ies with uncooperative patient groups such as infants. With our presented robust augmentation and open-source implementation, robust tractogram filtering is readily available.Peer reviewe

Dual-encoded magnetization transfer and diffusion imaging and its application to tract-specific microstructure mapping

We present a novel dual-encoded magnetization transfer (MT) and diffusion-weighted sequence and demonstrate its potential to resolve distinct properties of white matter fiber tracts at the sub-voxel level. The sequence was designed and optimized for maximal MT contrast efficiency. The resulting whole brain 2.6 mm isotropic protocol to measure tract-specific MT ratio (MTR) has a scan time under 7 minutes. Ten healthy subjects were scanned twice to assess repeatability. Two different analysis methods were contrasted: a technique to extract tract-specific MTR using Convex Optimization Modeling for Microstructure Informed Tractography (COMMIT), a global optimization technique; and conventional MTR tractometry. The results demonstrate that the tract-specific method can reliably resolve the MT ratios of major white matter fiber pathways and is less affected by partial volume effects than conventional multi-modal tractometry. Dual-encoded MT and diffusion is expected to both increase the sensitivity to microstructure alterations of specific tracts due to disease, ageing or learning, as well as lead to weighted structural connectomes with more anatomical specificity.Comment: 26 pages, 7 figure

A Two Alternative Forced Choice Method for Assessing Vibrotactile Discrimination Thresholds in The Lower Limb

The development of an easy to implement, quantitative measure to examine vibration perception would be useful for future application in clinical settings. Vibration sense in the lower limb of younger and older adults was examined using the method of constant stimuli (MCS) and the two-alternative forced choice paradigm. The focus of this experiment was to determine an appropriate stimulation site on the lower limb (tendon versus bone) to assess vibration threshold and to determine if the left and right legs have varying thresholds. Discrimination thresholds obtained at two stimulation sites in the left and right lower limbs showed differences in vibration threshold across the two ages groups, but not across sides of the body nor between stimulation sites within each limb. Overall, the MCS can be implemented simply, reliably, and with minimal time. It can also easily be implemented with low-cost technology. Therefore, it could be a good candidate method to assess the presence of specific deep sensitivity deficits in clinical practice, particularly in populations likely to show the onset of sensory deficits

Predicting MEG resting-state functional connectivity using microstructural information

Understanding how human brain microstructure influences functional connectivity is an important endeavor. In this work, magnetic resonance imaging data from ninety healthy participants were used to calculate structural connectivity matrices using the streamline count, fractional anisotropy, radial diffusivity and a myelin measure (derived from multi-component relaxometry) to assign connection strength. Unweighted binarized structural connectivity matrices were also constructed. Magnetoencephalography resting-state data from those participants were used to calculate functional connectivity matrices, via correlations of the Hilbert envelopes of beamformer timeseries in the delta, theta, alpha and beta frequency bands. Non-negative matrix factorization was performed to identify the components of the functional connectivity. Shortest-path-length and search-information analyses of the structural connectomes were used to predict functional connectivity patterns for each participant. The microstructure-informed algorithms predicted the components of the functional connectivity more accurately than they predicted the total functional connectivity. This provides a methodology to understand functional mechanisms better. The shortest-path-length algorithm exhibited the highest prediction accuracy. Of the weights of the structural connectivity matrices, the streamline count and the myelin measure gave the most accurate predictions, while the fractional anisotropy performed poorly. Overall, different structural metrics paint very different pictures of the structural connectome and its relationship to functional connectivity

Mapping tissue microstructure across the human brain on a clinical scanner with soma and neurite density image metrics

Soma and neurite density image (SANDI) is an advanced diffusion magnetic resonance imaging biophysical signal model devised to probe in vivo microstructural information in the gray matter (GM). This model requires acquisitions that include b values that are at least six times higher than those used in clinical practice. Such high b values are required to disentangle the signal contribution of water diffusing in soma from that diffusing in neurites and extracellular space, while keeping the diffusion time as short as possible to minimize potential bias due to water exchange. These requirements have limited the use of SANDI only to preclinical or cutting‐edge human scanners. Here, we investigate the potential impact of neglecting water exchange in the SANDI model and present a 10‐min acquisition protocol that enables to characterize both GM and white matter (WM) on 3 T scanners. We implemented analytical simulations to (i) evaluate the stability of the fitting of SANDI parameters when diminishing the number of shells; (ii) estimate the bias due to potential exchange between neurites and extracellular space in such reduced acquisition scheme, comparing it with the bias due to experimental noise. Then, we demonstrated the feasibility and assessed the repeatability and reproducibility of our approach by computing microstructural metrics of SANDI with AMICO toolbox and other state‐of‐the‐art models on five healthy subjects. Finally, we applied our protocol to five multiple sclerosis patients. Results suggest that SANDI is a practical method to characterize WM and GM tissues in vivo on performant clinical scanners

Evaluation of tractography-based myelin-weighted connectivity across the lifespan

IntroductionRecent studies showed that the myelin of the brain changes in the life span, and demyelination contributes to the loss of brain plasticity during normal aging. Diffusion-weighted magnetic resonance imaging (dMRI) allows studying brain connectivity in vivo by mapping axons in white matter with tractography algorithms. However, dMRI does not provide insight into myelin; thus, combining tractography with myelin-sensitive maps is necessary to investigate myelin-weighted brain connectivity. Tractometry is designated for this purpose, but it suffers from some serious limitations. Our study assessed the effectiveness of the recently proposed Myelin Streamlines Decomposition (MySD) method in estimating myelin-weighted connectomes and its capacity to detect changes in myelin network architecture during the process of normal aging. This approach opens up new possibilities compared to traditional Tractometry.MethodsIn a group of 85 healthy controls aged between 18 and 68 years, we estimated myelin-weighted connectomes using Tractometry and MySD, and compared their modulation with age by means of three well-known global network metrics.ResultsFollowing the literature, our results show that myelin development continues until brain maturation (40 years old), after which degeneration begins. In particular, mean connectivity strength and efficiency show an increasing trend up to 40 years, after which the process reverses. Both Tractometry and MySD are sensitive to these changes, but MySD turned out to be more accurate.ConclusionAfter regressing the known predictors, MySD results in lower residual error, indicating that MySD provides more accurate estimates of myelin-weighted connectivity than Tractometry

Tractostorm 2 : Optimizing tractography dissection reproducibility with segmentation protocol dissemination

The segmentation of brain structures is a key component of many neuroimaging studies. Consistent anatomical definitions are crucial to ensure consensus on the position and shape of brain structures, but segmentations are prone to variation in their interpretation and execution. White-matter (WM) pathways are global structures of the brain defined by local landmarks, which leads to anatomical definitions being difficult to convey, learn, or teach. Moreover, the complex shape of WM pathways and their representation using tractography (streamlines) make the design and evaluation of dissection protocols difficult and time-consuming. The first iteration of Tractostorm quantified the variability of a pyramidal tract dissection protocol and compared results between experts in neuroanatomy and nonexperts. Despite virtual dissection being used for decades, in-depth investigations of how learning or practicing such protocols impact dissection results are nonexistent. To begin to fill the gap, we evaluate an online educational tractography course and investigate the impact learning and practicing a dissection protocol has on interrater (groupwise) reproducibility. To generate the required data to quantify reproducibility across raters and time, 20 independent raters performed dissections of three bundles of interest on five Human Connectome Project subjects, each with four timepoints. Our investigation shows that the dissection protocol in conjunction with an online course achieves a high level of reproducibility (between 0.85 and 0.90 for the voxel-based Dice score) for the three bundles of interest and remains stable over time (repetition of the protocol). Suggesting that once raters are familiar with the software and tasks at hand, their interpretation and execution at the group level do not drastically vary. When compared to previous work that used a different method of communication for the protocol, our results show that incorporating a virtual educational session increased reproducibility. Insights from this work may be used to improve the future design of WM pathway dissection protocols and to further inform neuroanatomical definitions.Peer reviewe

Variability and reproducibility of multi-echo T2 relaxometry: Insights from multi-site, multi-session and multi-subject MRI acquisitions

Quantitative magnetic resonance imaging (qMRI) can increase the specificity and sensitivity of conventional weighted MRI to underlying pathology by comparing meaningful physical or chemical parameters, measured in physical units, with normative values acquired in a healthy population. This study focuses on multi-echo T2 relaxometry, a qMRI technique that probes the complex tissue microstructure by differentiating compartment-specific T2 relaxation times. However, estimation methods are still limited by their sensitivity to the underlying noise. Moreover, estimating the model's parameters is challenging because the resulting inverse problem is ill-posed, requiring advanced numerical regularization techniques. As a result, the estimates from distinct regularization strategies are different. In this work, we aimed to investigate the variability and reproducibility of different techniques for estimating the transverse relaxation time of the intra- and extra-cellular space (T2IE) in gray (GM) and white matter (WM) tissue in a clinical setting, using a multi-site, multi-session, and multi-run T2 relaxometry dataset. To this end, we evaluated three different techniques for estimating the T2 spectra (two regularized non-negative least squares methods and a machine learning approach). Two independent analyses were performed to study the effect of using raw and denoised data. For both the GM and WM regions, and the raw and denoised data, our results suggest that the principal source of variance is the inter-subject variability, showing a higher coefficient of variation (CoV) than those estimated for the inter-site, inter-session, and inter-run, respectively. For all reconstruction methods studied, the CoV ranged between 0.32 and 1.64%. Interestingly, the inter-session variability was close to the inter-scanner variability with no statistical differences, suggesting that T2IE is a robust parameter that could be employed in multi-site neuroimaging studies. Furthermore, the three tested methods showed consistent results and similar intra-class correlation (ICC), with values superior to 0.7 for most regions. Results from raw data were slightly more reproducible than those from denoised data. The regularized non-negative least squares method based on the L-curve technique produced the best results, with ICC values ranging from 0.72 to 0.92