Histological validation of per-bundle water diffusion metrics within a region of fiber crossing following axonal degeneration

Multi-tensor fixel-based metrics in tractometry: application to multiple sclerosis

Traditional Diffusion Tensor Imaging (DTI) metrics are affected by crossing fibers and lesions. Most of the previous tractometry works use the single diffusion tensor, which leads to limited sensitivity and challenging interpretation of the results in crossing fiber regions. In this work, we propose a tractometry pipeline that combines white matter tractography with multi-tensor fixel-based metrics. These multi-tensors are estimated using the stable, accurate and robust to noise Multi-Resolution Discrete Search method (MRDS). The spatial coherence of the multi-tensor field estimated with MRDS, which includes up to three anisotropic and one isotropic tensors, is tractography-regularized using the Track Orientation Density Imaging method. Our end-to-end tractometry pipeline goes from raw data to track-specific multi-tensor-metrics tract profiles that are robust to noise and crossing fibers. A comprehensive evaluation conducted in a phantom simulating healthy and damaged tissue with the standard model, as well as in a healthy cohort of 20 individuals scanned along 5 time points, demonstrates the advantages of using multi-tensor metrics over traditional single-tensor metrics in tractometry. Qualitative assessment in a cohort of patients with relapsing-remitting multiple sclerosis reveals that the pipeline effectively detects white matter anomalies in the presence of crossing fibers and lesions.

White matter structural bases for phase accuracy during tapping synchronization

We determined the intersubject association between the rhythmic entrainment abilities of human subjects during a synchronization-continuation tapping task (SCT) and the macro- and microstructural properties of their superficial (SWM) and deep (DWM) white matter. Diffusion-weighted images were obtained from 32 subjects who performed the SCT with auditory or visual metronomes and five tempos ranging from 550 to 950 ms. We developed a method to determine the density of short-range fibers that run underneath the cortical mantle, interconnecting nearby cortical regions (U-fibers). Notably, individual differences in the density of U-fibers in the right audiomotor system were correlated with the degree of phase accuracy between the stimuli and taps across subjects. These correlations were specific to the synchronization epoch with auditory metronomes and tempos around 1.5 Hz. In addition, a significant association was found between phase accuracy and the density and bundle diameter of the corpus callosum (CC), forming an interval-selective map where short and long intervals were behaviorally correlated with the anterior and posterior portions of the CC. These findings suggest that the structural properties of the SWM and DWM in the audiomotor system support the tapping synchronization abilities of subjects, as cortical U-fiber density is linked to the preferred tapping tempo and the bundle properties of the CC define an interval-selective topography.

Structural and functional connectivity in relation to executive functions in antipsychotic-naïve patients with first episode schizophrenia and levels of glutamatergic metabolites

Patients with schizophrenia exhibit structural and functional dysconnectivity but the relationship to the well-documented cognitive impairments is less clear. This study investigates associations between structural and functional connectivity and executive functions in antipsychotic-naïve patients experiencing schizophrenia. Sixty-four patients with schizophrenia and 95 matched controls underwent cognitive testing, diffusion weighted imaging and resting state functional magnetic resonance imaging. In the primary analyses, groupwise interactions between structural connectivity as measured by fixel-based analyses and executive functions were investigated using multivariate linear regression analyses. For significant structural connections, secondary analyses examined whether functional connectivity and associations with executive functions also differed for the two groups. In group comparisons, patients exhibited cognitive impairments across all executive functions compared to controls (p < 0.001), but no group difference were observed in the fixel-based measures. Primary analyses revealed a groupwise interaction between planning abilities and fixel-based measures in the left anterior thalamic radiation (p = 0.004), as well as interactions between cognitive flexibility and fixel-based measures in the isthmus of corpus callosum and cingulum (p = 0.049). Secondary analyses revealed increased functional connectivity between grey matter regions connected by the left anterior thalamic radiation (left thalamus with pars opercularis p = 0.018, and pars orbitalis p = 0.003) in patients compared to controls. Moreover, a groupwise interaction was observed between cognitive flexibility and functional connectivity between contralateral regions connected by the isthmus (precuneus p = 0.028, postcentral p = 0.012), all p-values corrected for multiple comparisons. We conclude that structural and functional connectivity appear to associate with executive functions differently in antipsychotic-naïve patients with schizophrenia compared to controls.

Fibre density and fibre-bundle cross-section of the corticospinal tract are distinctly linked to psychosis-specific symptoms in antipsychotic-naïve patients with first-episode schizophrenia

Multiple lines of research support the dysconnectivity hypothesis of schizophrenia. However, findings on white matter (WM) alterations in patients with schizophrenia are widespread and non-specific. Confounding factors from magnetic resonance image (MRI) processing, clinical diversity, antipsychotic exposure, and substance use may underlie some of the variability. By application of refined methodology and careful sampling, we rectified common confounders investigating WM and symptom correlates in a sample of strictly antipsychotic-naïve first-episode patients with schizophrenia. Eighty-six patients and 112 matched controls underwent diffusion MRI. Using fixel-based analysis (FBA), we extracted fibre-specific measures such as fibre density and fibre-bundle cross-section. Group differences on fixel-wise measures were examined with multivariate general linear modelling. Psychopathology was assessed with the Positive and Negative Syndrome Scale. We separately tested multivariate correlations between fixel-wise measures and predefined psychosis-specific versus anxio-depressive symptoms. Results were corrected for multiple comparisons. Patients displayed reduced fibre density in the body of corpus callosum and in the middle cerebellar peduncle. Fibre density and fibre-bundle cross-section of the corticospinal tract were positively correlated with suspiciousness/persecution, and negatively correlated with delusions. Fibre-bundle cross-section of isthmus of corpus callosum and hallucinatory behaviour were negatively correlated. Fibre density and fibre-bundle cross-section of genu and splenium of corpus callosum were negative correlated with anxio-depressive symptoms. FBA revealed fibre-specific properties of WM abnormalities in patients and differentiated associations between WM and psychosis-specific versus anxio-depressive symptoms. Our findings encourage an itemised approach to investigate the relationship between WM microstructure and clinical symptoms in patients with schizophrenia.

Differentiation of white matter histopathology using b-tensor encoding and machine learning

Diffusion-weighted magnetic resonance imaging (DW-MRI) is a non-invasive technique that is sensitive to microstructural geometry in neural tissue and is useful for the detection of neuropathology in research and clinical settings. Tensor-valued diffusion encoding schemes (b-tensor) have been developed to enrich the microstructural data that can be obtained through DW-MRI. These advanced methods have proven to be more specific to microstructural properties than conventional DW-MRI acquisitions. Additionally, machine learning methods are particularly useful for the study of multidimensional data sets. In this work, we have tested the reach of b-tensor encoding data analyses with machine learning in different histopathological scenarios. We achieved this in three steps: 1) We induced different levels of white matter damage in rodent optic nerves. 2) We obtained ex vivo DW-MRI data with b-tensor encoding schemes and calculated quantitative metrics using Q-space trajectory imaging. 3) We used a machine learning model to identify the main contributing features and built a voxel-wise probabilistic classification map of histological damage. Our results show that this model is sensitive to characteristics of microstructural damage. In conclusion, b-tensor encoded DW-MRI data analyzed with machine learning methods, have the potential to be further developed for the detection of histopathology and neurodegeneration.

Investigating Tissue-Specific Abnormalities in Alzheimer's Disease with Multi-Shell Diffusion MRI

BACKGROUND

Most studies using diffusion-weighted MRI (DW-MRI) in Alzheimer's disease (AD) have focused their analyses on white matter (WM) microstructural changes using the diffusion (kurtosis) tensor model. Although recent works have addressed some limitations of the tensor model, such as the representation of crossing fibers and partial volume effects with cerebrospinal fluid (CSF), the focus remains in modeling and analyzing the WM.

OBJECTIVE

In this work, we present a brain analysis approach for DW-MRI that disentangles multiple tissue compartments as well as micro- and macroscopic effects to investigate differences between groups of subjects in the AD continuum and controls.

METHODS

By means of the multi-tissue constrained spherical deconvolution of multi-shell DW-MRI, underlying brain tissue is modeled with a WM fiber orientation distribution function along with the contributions of gray matter (GM) and CSF to the diffusion signal. From this multi-tissue model, a set of measures capturing tissue diffusivity properties and morphology are extracted. Group differences were interrogated following fixel-, voxel-, and tensor-based morphometry approaches while including strong FWE control across multiple comparisons.

RESULTS

Abnormalities related to AD stages were detected in WM tracts including the splenium, cingulum, longitudinal fasciculi, and corticospinal tract. Changes in tissue composition were identified, particularly in the medial temporal lobe and superior longitudinal fasciculus.

CONCLUSION

This analysis framework constitutes a comprehensive approach allowing simultaneous macro and microscopic assessment of WM, GM, and CSF, from a single DW-MRI dataset.

Using Advanced Diffusion-Weighted Imaging to Predict Cell Counts in Gray Matter: Potential and Pitfalls

Recent advances in diffusion imaging have given it the potential to non-invasively detect explicit neurobiological properties, beyond what was previously possible with conventional structural imaging. However, there is very little known about what cytoarchitectural properties these metrics, especially those derived from newer multi-shell models like Neurite Orientation Dispersion and Density Imaging (NODDI) correspond to. While these diffusion metrics do not promise any inherent cell type specificity, different brain cells have varying morphologies, which could influence the diffusion signal in distinct ways. This relationship is currently not well-characterized. Understanding the possible cytoarchitectural signatures of diffusion measures could allow them to estimate important neurobiological properties like cell counts, potentially resulting in a powerful clinical diagnostic tool. Here, using advanced diffusion imaging (NODDI) in the mouse brain, we demonstrate that different regions have unique relationships between cell counts and diffusion metrics. We take advantage of this exclusivity to introduce a framework to predict cell counts of different types of cells from the diffusion metrics alone, in a region-specific manner. We also outline the challenges of reliably developing such a model and discuss the precautions the field must take when trying to tie together medical imaging modalities and histology.

Longitudinal changes of ADHD symptoms in association with white matter microstructure: A tract-specific fixel-based analysis

•

HI symptom remission is associated with more follow-up lCST FD.

•

Combined symptom remission is associated with more follow-up lCST FC.

•

Altered white matter development may be moderated by preceding symptom trajectory.

Background

Variation in the longitudinal course of childhood attention deficit/hyperactivity disorder (ADHD) coincides with neurodevelopmental maturation of brain structure and function. Prior work has attempted to determine how alterations in white matter (WM) relate to changes in symptom severity, but much of that work has been done in smaller cross-sectional samples using voxel-based analyses. Using standard diffusion-weighted imaging (DWI) methods, we previously showed WM alterations were associated with ADHD symptom remission over time in a longitudinal sample of probands, siblings, and unaffected individuals. Here, we extend this work by further assessing the nature of these changes in WM microstructure by including an additional follow-up measurement (aged 18 – 34 years), and using the more physiologically informative fixel-based analysis (FBA).

Methods

Data were obtained from 139 participants over 3 clinical and 2 follow-up DWI waves, and analyzed using FBA in regions-of-interest based on prior findings. We replicated previously reported significant models and extended them by adding another time-point, testing whether changes in combined ADHD and hyperactivity-impulsivity (HI) continuous symptom scores are associated with fixel metrics at follow-up.

Results

Clinical improvement in HI symptoms over time was associated with more fiber density at follow-up in the left corticospinal tract (lCST) (t_max = 1.092, standardized effect[SE] = 0.044, p_FWE = 0.016). Improvement in combined ADHD symptoms over time was associated with more fiber cross-section at follow-up in the lCST (t_max = 3.775, SE = 0.051, p_FWE = 0.019).

Conclusions

Aberrant white matter development involves both lCST micro- and macrostructural alterations, and its path may be moderated by preceding symptom trajectory.

Brain microstructural antecedents of visual difficulties in infants born very preterm

Infants born very preterm (VPT) are at risk of later visual problems. Although neonatal screening can identify ophthalmologic abnormalities, subtle perinatal brain injury and/or delayed brain maturation may be significant contributors to complex visual-behavioral problems. Our aim was to assess the micro and macrostructural antecedents of early visual-behavioral difficulties in VPT infants by using diffusion MRI (dMRI) at term-equivalent age. We prospectively recruited a cohort of 262 VPT infants (≤32 weeks gestational age [GA]) from five neonatal intensive care units. We obtained structural and diffusion MRI at term-equivalent age and administered the Preverbal Visual Assessment (PreViAs) questionnaire to parents at 3-4 months corrected age. We used constrained spherical deconvolution to reconstruct nine white matter tracts of the visual pathways with high reliability and performed fixel-based analysis to derive fiber density (FD), fiber-bundle cross-section (FC), and combined fiber density and cross-section (FDC). In multiple logistic regression analyses, we related these tract metrics to visual-behavioral function. Of 262 infants, 191 had both high-quality dMRI and completed PreViAs, constituting the final cohort: mean (SD) GA was 29.3 (2.4) weeks, 90 (47.1%) were males, and postmenstrual age (PMA) at MRI was 42.8 (1.3) weeks. FD and FC of several tracts were altered in infants with (N = 59) versus those without retinopathy of prematurity (N = 132). FDC of the left posterior thalamic radiations (PTR), left inferior longitudinal fasciculus (ILF), right superior longitudinal fasciculus (SLF), and left inferior fronto-occipital fasciculus (IFOF) were significantly associated with visual attention scores, prior to adjusting for confounders. After adjustment for PMA at MRI, GA, severe retinopathy of prematurity, and total brain volume, FDC of the left PTR, left ILF, and left IFOF remained significantly associated with visual attention. Early visual-behavioral difficulties in VPT infants are preceded by micro and macrostructural abnormalities in several major visual pathways at term-equivalent age.

White matter degeneration revealed by fiber-specific analysis relates to recovery of hand function after stroke

Recent developments of higher‐order diffusion‐weighted imaging models have enabled the estimation of specific white matter fiber populations within a voxel, addressing limitations of traditional imaging markers of white matter integrity. We applied fixel based analysis (FBA) to investigate the evolution of fiber‐specific white matter changes in a prospective study of stroke patients and upper limb motor deficit over 1 year after stroke. We studied differences in fiber density and macrostructural changes in fiber cross‐section. Motor function was assessed by grip strength. We conducted a whole‐brain analysis of fixel metrics and predefined corticospinal tract (CST) region of interest in relation to changes in motor functions. In 30 stroke patients (mean age 62.3 years, SD ±16.9; median NIHSS 4, IQR 2–5), whole‐brain FBA revealed progressing loss of fiber density and cross‐section in the ipsilesional corticospinal tract and long‐range fiber tracts such as the superior longitudinal fascicle and trans‐callosal tracts extending towards contralesional white matter tracts. Lower FBA metrics measured at the brainstem section of the CST 1 month after stroke were significantly associated with lower grip strength 3 months (p = .009, adjusted R² = 0.259) and 1 year (T4: p < .001, adj. R² = 0.515) after stroke. Compared to FA, FBA metrics showed a comparably strong association with grip strength at later time points. Using FBA, we demonstrate progressive fiber‐specific white matter loss after stroke and association with functional motor outcome. Our results promote the application of fiber‐specific analysis to detect secondary neurodegeneration after stroke in relation to clinical recovery.

Early micro- and macrostructure of sensorimotor tracts and development of cerebral palsy in high risk infants

Infants born very preterm (VPT) are at high risk of motor impairments such as cerebral palsy (CP), and diagnosis can take 2 years. Identifying in vivo determinants of CP could facilitate presymptomatic detection and targeted intervention. Our objectives were to derive micro‐ and macrostructural measures of sensorimotor white matter tract integrity from diffusion MRI at term‐equivalent age, and determine their association with early diagnosis of CP. We enrolled 263 VPT infants (≤32 weeks gestational age) as part of a large prospective cohort study. Diffusion and structural MRI were acquired at term. Following consensus guidelines, we defined early diagnosis of CP based on abnormal structural MRI at term and abnormal neuromotor exam at 3–4 months corrected age. Using Constrained Spherical Deconvolution, we derived a white matter fiber orientation distribution (fOD) for subjects, performed probabilistic whole‐brain tractography, and segmented nine sensorimotor tracts of interest. We used the recently developed fixel‐based (FB) analysis to compute fiber density (FD), fiber‐bundle cross‐section (FC), and combined fiber density and cross‐section (FDC) for each tract. Of 223 VPT infants with high‐quality diffusion MRI data, 14 (6.3%) received an early diagnosis of CP. The cohort's mean (SD) gestational age was 29.4 (2.4) weeks and postmenstrual age at MRI scan was 42.8 (1.3) weeks. FD, FC, and FDC for each sensorimotor tract were significantly associated with early CP diagnosis, with and without adjustment for confounders. Measures of sensorimotor tract integrity enhance our understanding of white matter changes that antecede and potentially contribute to the development of CP in VPT infants.

Use of computational fluid dynamics for 3D fiber tract visualization on human high-thickness histological slices: histological mesh tractography

Understanding the intricate three-dimensional relationship between fiber bundles and subcortical nuclei is not a simple task. It is of paramount importance in neurosciences, especially in the field of functional neurosurgery. The current methods for in vivo and post mortem fiber tract visualization have shortcomings and contributions to the field are welcome. Several tracts were chosen to implement a new technique to help visualization of white matter tracts, using high-thickness histology and dark field images. Our study describes the use of computational fluid dynamic simulations for visualization of 3D fiber tracts segmented from dark field microscopy in high-thickness histological slices (histological mesh tractography). A post mortem human brain was MRI scanned prior to skull extraction, histologically processed and serially cut at 430 µm thickness as previously described by our group. High-resolution dark field images were used to segment the outlines of the structures. These outlines served as basis for the construction of a 3D structured mesh, were a Finite Volume Method (FVM) simulation of water flow was performed to generate streamlines representing the geometry. The simulations were accomplished by an open source computer fluid dynamics software. The resulting simulation rendered a realistic 3D impression of the segmented anterior commissure, the left anterior limb of the internal capsule, the left uncinate fascicle, and the dentato-rubral tracts. The results are in line with clinical findings, diffusion MR imaging and anatomical dissection methods.

Assessment of the structural complexity of diffusion MRI voxels using 3D electron microscopy in the rat brain

Validation and interpretation of diffusion magnetic resonance imaging (dMRI) requires detailed understanding of the actual microstructure restricting the diffusion of water molecules. In this study, we used serial block-face scanning electron microscopy (SBEM), a three-dimensional electron microscopy (3D-EM) technique, to image seven white and grey matter volumes in the rat brain. SBEM shows excellent contrast of cellular membranes, which are the major components restricting the diffusion of water in tissue. Additionally, we performed 3D structure tensor (3D-ST) analysis on the SBEM volumes and parameterised the resulting orientation distributions using Watson and angular central Gaussian (ACG) probability distributions as well as spherical harmonic (SH) decomposition. We analysed how these parameterisations described the underlying orientation distributions and compared their orientation and dispersion with corresponding parameters from two dMRI methods, neurite orientation dispersion and density imaging (NODDI) and constrained spherical deconvolution (CSD). Watson and ACG parameterisations and SH decomposition captured well the 3D-ST orientation distributions, but ACG and SH better represented the distributions due to its ability to model asymmetric dispersion. The dMRI parameters corresponded well with the 3D-ST parameters in the white matter volumes, but the correspondence was less evident in the more complex grey matter. SBEM imaging and 3D-ST analysis also revealed that the orientation distributions were often not axially symmetric, a property neatly captured by the ACG distribution. Overall, the ability of SBEM to image diffusion barriers in intricate detail, combined with 3D-ST analysis and parameterisation, provides a step forward toward interpreting and validating the dMRI signals in complex brain tissue microstructure.

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/mm². 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/mm², median R² = .34; b = 4,000 s/mm², median R² = .29; b = 2,400 s/mm², median R² = .21; b = 1,200 s/mm², median R² = .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.

Histological and diffusion-weighted magnetic resonance imaging data from normal and degenerated optic nerve and chiasm of the rat

Diffusion-weighted magnetic resonance imaging (dMRI) is widely used to infer microstructural characteristics of tissue, particularly in cerebral white matter. Histological validation of the metrics derived from dMRI methods are needed to fully characterize their ability to capture biologically-relevant histological features non-invasively. The data described here were used to correlate metrics derived from dMRI and quantitative histology in an animal model of axonal degeneration (“Histological validation of per-bundle water diffusion metrics within a region of fiber crossing following axonal degeneration” [1]). Unilateral retinal ischemia/reperfusion was induced in 10 rats, by the elevation of pressure of the anterior chamber of the eye for 90 min. Five rats were used as controls. After five weeks, injured animals were intracardially perfused to analyze the optic nerves and chiasm with dMRI and histology. This resulted in 15 brain scans, each with 80 diffusion-sensitizing gradient directions with b = 2000 and 2500 s/mm² and 20 non-diffusion-weighted images (b = 0 s/mm²), with isometric voxel resolution of 125 μm³. Histological sections were obtained after dMRI. Optical microscopy photomicrographs of the optic nerves (stained with toluidine blue) are available, as well as their corresponding automatic segmentations of axons and myelin.