Along-axon diameter variation and axonal orientation dispersion revealed with 3D electron microscopy: implications for quantifying brain white matter microstructure with histology and diffusion MRI

Quantifying axonal features of human superficial white matter from three-dimensional multibeam serial electron microscopy data assisted by deep learning

Short-range association fibers located in the superficial white matter play an important role in mediating higher-order cognitive function in humans. Detailed morphological characterization of short-range association fibers at the microscopic level promises to yield important insights into the axonal features driving cortico-cortical connectivity in the human brain yet has been difficult to achieve to date due to the challenges of imaging at nanometer-scale resolution over large tissue volumes. This work presents results from multi-beam scanning electron microscopy (EM) data acquired at 4 × 4 × 33 nm3 resolution in a volume of human superficial white matter measuring 200 × 200 × 112 μm3, leveraging automated analysis methods. Myelin and myelinated axons were automatically segmented using deep convolutional neural networks (CNNs), assisted by transfer learning and dropout regularization techniques. A total of 128,285 myelinated axons were segmented, of which 70,321 and 2102 were longer than 10 and 100 μm, respectively. Marked local variations in diameter (i.e., beading) and direction (i.e., undulation) were observed along the length of individual axons. Myelinated axons longer than 10 μm had inner diameters around 0.5 µm, outer diameters around 1 µm, and g-ratios around 0.5. This work fills a gap in knowledge of axonal morphometry in the superficial white matter and provides a large 3D human EM dataset and accurate segmentation results for a variety of future studies in different fields.

Engineering clinical translation of OGSE diffusion MRI

Oscillating gradient spin echo (OGSE) diffusion MRI (dMRI) can probe the diffusive dynamics on short time scales ≲10 ms, which translates into the sensitivity to tissue microstructure at the short length scales≲ 10 μ $$ \lesssim 10\kern0.3em \upmu $$ m. OGSE-based tissue microstructure imaging techniques able to characterize the cell diameter and cellular density have been established in pre-clinical studies. The unique image contrast of OGSE dMRI has been shown to differentiate tumor types and malignancies, enable early diagnosis of treatment effectiveness, and reveal different pathophysiology of lesions in stroke and neurological diseases. Recent innovations in high-performance gradient human MRI systems provide an opportunity to translate OGSE research findings in pre-clinical studies to human research and the clinic. The implementation of OGSE dMRI in human studies has the promise to advance our understanding of human brain microstructure and improve patient care. Compared to the clinical standard (pulsed gradient spin echo), engineering OGSE diffusion encoding for human imaging is more challenging. This review summarizes the impact of hardware and human biophysical safety considerations on the waveform design, imaging parameter space, and image quality of OGSE dMRI. Here we discuss the effects of the gradient amplitude, slew rate, peripheral nerve stimulation, cardiac stimulation, gradient driver, acoustic noise and mechanical vibration, eddy currents, gradient nonlinearity, concomitant gradient, motion and flow, and signal-to-noise ratio. We believe that targeted engineering for safe, high-quality, and reproducible imaging will enable the translation of OGSE dMRI techniques into the clinic.

Comparison of different neurite density metrics with brain asymmetry evaluation

The standard diffusion MRI model with intra- and extra-axonal water pools offers a set of microstructural parameters describing brain white matter architecture. However, non-linearities in the standard model and diffusion data contamination by noise and imaging artefacts make estimation of diffusion metrics challenging. In order to develop reliable diffusion approaches and to avoid computational model degeneracy, additional theoretical assumptions allowing stable numerical implementations are required. Advanced diffusion approaches allow for estimation of intra-axonal water fraction (AWF), describing a key structural characteristic of brain tissue. AWF can be interpreted as an indirect measure or proxy of neurite density and has a potential as useful clinical biomarker. Established diffusion approaches such as white matter tract integrity, neurite orientation dispersion and density imaging (NODDI), and spherical mean technique provide estimates of AWF within their respective theoretical frameworks. In the present study, we estimated AWF metrics using different diffusion approaches and compared measures of brain asymmetry between the different metrics in a sub-sample of 182 subjects from the UK Biobank. Multivariate decomposition by mean of linked independent component analysis revealed that the various AWF proxies derived from the different diffusion approaches reflect partly non-overlapping variance of independent components, with distinct anatomical distributions and sensitivity to age. Further, voxel-wise analysis revealed age-related differences in AWF-based brain asymmetry, indicating less apparent left-right hemisphere difference with higher age. Finally, we demonstrated that NODDI metrics suffer from a quite strong dependence on used numerical algorithms and post-processing pipeline. The analysis based on AWF metrics strongly depends on the used diffusion approach and leads to poorly reproducible results.

Large-scale EM data reveals myelinated axonal changes and altered connectivity in the corpus callosum of an autism mouse model

Introduction

Autism spectrum disorder (ASD) encompasses a diverse range of neurodevelopmental disorders with complex etiologies, including genetic, environmental, and neuroanatomical factors. While the exact mechanisms underlying ASD remain unclear, structural abnormalities in the brain offer valuable insights into its pathophysiology. The corpus callosum, the largest white matter tract in the brain, plays a crucial role in interhemispheric communication, and its structural abnormalities may contribute to ASD-related phenotypes.

Methods

To investigate the ultrastructural alterations in the corpus callosum associated with ASD, we utilized serial scanning electron microscopy (sSEM) in mice. A dataset of the entire sagittal sections of the corpus callosum from wild-type and Shank3B mutant mice was acquired at 4 nm resolution, enabling precise comparisons of myelinated axon properties. Leveraging a fine-tuned EM-SAM model for automated segmentation, we quantitatively analyzed key metrics, including G-ratio, myelin thickness, and axonal density.

Results

In the corpus callosum of Shank3B autism model mouse, we observed a significant increase in myelinated axon density, accompanied by thinner myelin sheaths compared to wild-type. Additionally, we identified abnormalities in the diameter distribution of myelinated axons and deviations in G-ratio. Notably, these ultrastructural alterations were widespread across the corpus callosum, suggesting a global disruption of myelinated axon integrity.

Discussion

This study provides novel insights into the microstructural abnormalities of the corpus callosum in ASD mouse, supporting the hypothesis that myelination deficits contribute to ASD-related communication impairments between brain hemispheres. However, given the structural focus of this study, further research integrating functional assessments is necessary to establish a direct link between these morphological changes and ASD-related neural dysfunction.

Bridging the 3D geometrical organisation of white matter pathways across anatomical length scales and species

We used diffusion MRI and x-ray synchrotron imaging on monkey and mice brains to examine the organisation of fibre pathways in white matter across anatomical scales. We compared the structure in the corpus callosum and crossing fibre regions and investigated the differences in cuprizone-induced demyelination in mouse brains versus healthy controls. Our findings revealed common principles of fibre organisation that apply despite the varying patterns observed across species; small axonal fasciculi and major bundles formed laminar structures with varying angles, according to the characteristics of major pathways. Fasciculi exhibited non-straight paths around obstacles like blood vessels, comparable across the samples of varying fibre complexity and demyelination. Quantifications of fibre orientation distributions were consistent across anatomical length scales and modalities, whereas tissue anisotropy had a more complex relationship, both dependent on the field-of-view. Our study emphasises the need to balance field-of-view and voxel size when characterising white matter features across length scales.

Scattering approach to diffusion quantifies axonal damage in brain injury

Early diagnosis and noninvasive monitoring of neurological disorders require sensitivity to elusive cellular-level alterations that occur much earlier than volumetric changes observable with the millimeter-resolution of medical imaging modalities. Morphological changes in axons, such as axonal varicosities or beadings, are observed in neurological disorders, as well as in development and aging. Here, we reveal the sensitivity of time-dependent diffusion MRI (dMRI) to axonal morphology at the micrometer scale. Scattering theory uncovers the two parameters that determine the diffusive dynamics of water in axons: the average reciprocal cross-section and the variance of long-range cross-sectional fluctuations. This theoretical development allowed us to predict dMRI metrics sensitive to axonal alterations across tens of thousands of axons in seconds rather than months of simulations in a rat model of traumatic brain injury. Our approach bridges the gap between micrometers and millimeters in resolution, offering quantitative, objective biomarkers applicable to a broad spectrum of neurological disorders.

Susceptibility-induced internal gradients reveal axon morphology and cause anisotropic effects in the diffusion-weighted MRI signal

Diffusion-weighted MRI is our most promising method for estimating microscopic tissue morphology in vivo. The signal acquisition is based on scanner-generated external magnetic gradients. However, it will also be affected by susceptibility-induced internal magnetic gradients caused by interactions between the tissue and the static magnetic field of the scanner. With 3D in silico experiments, we show how internal gradients cause morphology-, compartment-, and orientation-dependence of spin-echo and pulsed-gradient spin-echo experiments in myelinated axons. These effects surpass those observed with previous 2D modelling corresponding to straight cylinders. For an ex vivo monkey brain, we observe the orientation-dependence generated only when including non-circular cross-sections in the in silico morphological configurations, and find orientation-dependent deviation of up to 17% for diffusion tensor metrics. Interestingly, we find that the orientation-dependence not only biases the signal across different brain regions, but also carries a sensitivity to the morphology of axonal cross-sections which is not attainable by the idealised theoretical diffusion-weighted MRI signal.

Time to bury the chisel: a continuous dorsal association tract system

The arcuate fasciculus may be subdivided into a tract directly connecting frontal and temporal lobes and a pair of indirect subtracts in which the fronto-temporal connection is mediated by connections to the inferior parietal lobe. This tripartition has been advanced as an improvement over the centuries-old consensus that the lateral dorsal association fibers form a continuous system with no discernible discrete parts. Moreover, it has been used as the anatomical basis for functional hypotheses regarding linguistic abilities. Ex hypothesi, damage to the indirect subtracts leads to deficits in the repetition of multi-word sequences, whereas damage to the direct subtract leads to deficits in the immediate reproduction of single multisyllabic words. We argue that this partitioning of the dorsal association tract system enjoys no special anatomical status, and the search for the anatomical substrates of linguistic abilities should not be constrained by it. Instead, the merit of any postulated partitioning should primarily be judged on the basis of whether it enlightens or obfuscates our understanding of the behavior of patients in which individual subtracts are damaged.

A survey on integral equations for bioelectric modeling

Bioelectric modeling problems, such as electroencephalography, magnetoencephalography, transcranial electrical stimulation, deep brain stimulation, and transcranial magnetic stimulation, among others, can be approached through the formulation and resolution of integral equations of theboundary element method(BEM). Recently, it has been realized that thecharge-based formulationof the BEM is naturally well-suited for the application of thefast multipole method(FMM). The FMM is a powerful algorithm for the computation of many-body interactions and is widely applied in electromagnetic modeling problems. With the introduction of BEM-FMM in the context of bioelectromagnetism, the BEM can now compete with thefinite element method(FEM) in a number of application cases. This survey has two goals: first, to give a modern account of the main BEM formulations in the literature and their integration with FMM, directed to general researchers involved in development of BEM software for bioelectromagnetic applications. Second, to survey different techniques and available software, and to contrast different BEM and FEM approaches. As a new contribution, we showcase that the charge-based formulation is dual to the more common surface potential formulation.

Exploring white matter dynamics and morphology through interactive numerical phantoms: the White Matter Generator

Brain white matter is a dynamic environment that continuously adapts and reorganizes in response to stimuli and pathological changes. Glial cells, especially, play a key role in tissue repair, inflammation modulation, and neural recovery. The movements of glial cells and changes in their concentrations can influence the surrounding axon morphology. We introduce the White Matter Generator (WMG) tool to enable the study of how axon morphology is influenced through such dynamical processes, and how this, in turn, influences the diffusion-weighted MRI signal. This is made possible by allowing interactive changes to the configuration of the phantom generation throughout the optimization process. The phantoms can consist of myelinated axons, unmyelinated axons, and cell clusters, separated by extra-cellular space. Due to morphological flexibility and computational advantages during the optimization, the tool uses ellipsoids as building blocks for all structures; chains of ellipsoids for axons, and individual ellipsoids for cell clusters. After optimization, the ellipsoid representation can be converted to a mesh representation which can be employed in Monte-Carlo diffusion simulations. This offers an effective method for evaluating tissue microstructure models for diffusion-weighted MRI in controlled bio-mimicking white matter environments. Hence, the WMG offers valuable insights into white matter's adaptive nature and implications for diffusion-weighted MRI microstructure models, and thereby holds the potential to advance clinical diagnosis, treatment, and rehabilitation strategies for various neurological disorders and injuries.

In vivo mapping of cellular resolution neuropathology in brain ischemia with diffusion MRI

Noninvasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent advances in diffusion magnetic resonance imaging enabled in vivo examination of tissue microstructures well beyond the imaging resolution. Here, we proposed to use diffusion time-dependent diffusion kurtosis imaging (tDKI) to simultaneously assess cellular morphology and transmembrane permeability in hypoxic-ischemic (HI) brain injury. Through numerical simulations and organoid imaging, we demonstrated the feasibility of capturing effective size and permeability changes using tDKI. In vivo MRI of HI-injured mouse brains detected a shift of the tDKI peak to longer diffusion times, suggesting swelling of the cellular processes. Furthermore, we observed a faster decrease of the tDKI tail, reflecting increased transmembrane permeability associated with up-regulated water exchange or necrosis. Such information, unavailable from a single diffusion time, can predict salvageable tissues. Preliminary applications of tDKI in patients with ischemic stroke suggested increased transmembrane permeability in stroke regions, illustrating tDKI's potential for detecting pathological changes in the clinics.

The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

A translational MRI approach to validate acute axonal damage detection as an early event in multiple sclerosis

Axonal degeneration is a central pathological feature of multiple sclerosis and is closely associated with irreversible clinical disability. Current noninvasive methods to detect axonal damage in vivo are limited in their specificity and clinical applicability, and by the lack of proper validation. We aimed to validate an MRI framework based on multicompartment modeling of the diffusion signal (AxCaliber) in rats in the presence of axonal pathology, achieved through injection of a neurotoxin damaging the neuronal terminal of axons. We then applied the same MRI protocol to map axonal integrity in the brain of multiple sclerosis relapsing-remitting patients and age-matched healthy controls. AxCaliber is sensitive to acute axonal damage in rats, as demonstrated by a significant increase in the mean axonal caliber along the targeted tract, which correlated with neurofilament staining. Electron microscopy confirmed that increased mean axonal diameter is associated with acute axonal pathology. In humans with multiple sclerosis, we uncovered a diffuse increase in mean axonal caliber in most areas of the normal-appearing white matter, preferentially affecting patients with short disease duration. Our results demonstrate that MRI-based axonal diameter mapping is a sensitive and specific imaging biomarker that links noninvasive imaging contrasts with the underlying biological substrate, uncovering generalized axonal damage in multiple sclerosis as an early event.

Estimations of Charge Deposition Onto Convoluted Axon Surfaces Within Extracellular Electric Fields

OBJECTIVE

Biophysical models of neural stimulation are a valuable approach to explaining the mechanisms of neuronal recruitment via applied extracellular electric fields. Typically, the applied electric field is estimated via a macroscopic finite element method solution and then applied to cable models as an extracellular voltage source. However, the field resolution is limited by the finite element size (typically 10's-100's of times greater than average neuronal cross-section). As a result, induced charges deposited onto anatomically realistic curved membrane interfaces are not taken into consideration. However, these details may alter estimates of the applied electric field and predictions of neural tissue activation.

METHODS

To estimate microscopic variations of the electric field, data for intra-axonal space segmented from 3D scanning electron microscopy of the mouse brain genu of corpus callosum were used. The boundary element fast multipole method was applied to accurately compute the extracellular solution. Neuronal recruitment was then estimated via an activating function.

RESULTS

Taking the physical structure of the arbor into account generally predicts higher values of the activating function. The relative integral 2-norm difference is 90% on average when the entire axonal arbor is present. A large fraction of this difference might be due to the axonal body itself. When an isolated physical axon is considered with all other axons removed, the relative integral 2-norm difference between the single-axon solution and the complete solution is 25% on average.

CONCLUSION

Our result may provide an explanation as to why Deep Brain Stimulation experiments typically predict lower activation thresholds than commonly used FEM/Cable model approaches to predicting neuronal responses to extracellular electrical stimulation.

SIGNIFICANCE

These results may change methods for bi-domain neural modeling and neural excitation.

Multiscale label-free imaging of myelin in human brain tissue with polarization-sensitive optical coherence tomography and birefringence microscopy

The combination of polarization-sensitive optical coherence tomography (PS-OCT) and birefringence microscopy (BRM) enables multiscale assessment of myelinated axons in postmortem brain tissue, and these tools are promising for the study of brain connectivity and organization. We demonstrate label-free imaging of myelin structure across the mesoscopic and microscopic spatial scales by performing serial-sectioning PS-OCT of a block of human brain tissue and periodically sampling thin sections for high-resolution imaging with BRM. In co-registered birefringence parameter maps, we observe good correspondence and demonstrate that BRM enables detailed validation of myelin (hence, axonal) organization, thus complementing the volumetric information content of PS-OCT.

Origins of brain tissue elasticity under multiple loading modes by analyzing the microstructure-based models

Constitutive behaviors and material properties of brain tissue play an essential role in accurately modeling its mechanical responses. However, the measured mechanical behaviors of brain tissue exhibit a large variability, and the reported elastic modulus can differ by orders of magnitude. Here we develop the micromechanical models based on the actual microstructure of the longitudinally anisotropic plane of brain tissue to investigate the microstructural origins of the large variability. Specifically, axonal fiber bundles with the specified configurations are distributed in an equivalent matrix. All micromechanical models are subjected to multiple loading modes, such as tensile, compressive, and shear loading, under periodic boundary conditions. The predicted results agree well with the experimental results. Furthermore, we investigate how brain tissue elasticity varies with its microstructural features. It is revealed that the large variability in brain tissue elasticity stems from the volume fraction of axonal fiber, the aspect ratio of axonal fiber, and the distribution of axonal fiber orientation. The volume fraction has the greatest impact on the mechanical behaviors of brain tissue, followed by the distribution of axonal fiber orientation, then the aspect ratio. This study provides critical insights for understanding the microstructural origins of the large variability in brain tissue elasticity.

A Method for In-Vivo Mapping of Axonal Diameter Distributions in the Human Brain Using Diffusion-Based Axonal Spectrum Imaging (AxSI)

In this paper we demonstrate a generalized and simplified pipeline called axonal spectrum imaging (AxSI) for in-vivo estimation of axonal characteristics in the human brain. Whole-brain estimation of the axon diameter, in-vivo and non-invasively, across all fiber systems will allow exploring uncharted aspects of brain structure and function relations with emphasis on connectivity and connectome analysis. While axon diameter mapping is important in and of itself, its correlation with conduction velocity will allow, for the first time, the explorations of information transfer mechanisms within the brain. We demonstrate various well-known aspects of axonal morphometry (e.g., the corpus callosum axon diameter variation) as well as other aspects that are less explored (e.g., axon diameter-based separation of the superior longitudinal fasciculus into segments). Moreover, we have created an MNI based mean axon diameter map over the entire brain for a large cohort of subjects providing the reference basis for future studies exploring relation between axon properties, its connectome representation, and other functional and behavioral aspects of the brain.

Quantitative MRI in leukodystrophies

Leukodystrophies constitute a large and heterogeneous group of genetic diseases primarily affecting the white matter of the central nervous system. Different disorders target different white matter structural components. Leukodystrophies are most often progressive and fatal. In recent years, novel therapies are emerging and for an increasing number of leukodystrophies trials are being developed. Objective and quantitative metrics are needed to serve as outcome measures in trials. Quantitative MRI yields information on microstructural properties, such as myelin or axonal content and condition, and on the chemical composition of white matter, in a noninvasive fashion. By providing information on white matter microstructural involvement, quantitative MRI may contribute to the evaluation and monitoring of leukodystrophies. Many distinct MR techniques are available at different stages of development. While some are already clinically applicable, others are less far developed and have only or mainly been applied in healthy subjects. In this review, we explore the background, current status, potential and challenges of available quantitative MR techniques in the context of leukodystrophies.

Merged magnetic resonance and light sheet microscopy of the whole mouse brain

We have developed workflows to align 3D magnetic resonance histology (MRH) of the mouse brain with light sheet microscopy (LSM) and 3D delineations of the same specimen. We start with MRH of the brain in the skull with gradient echo and diffusion tensor imaging (DTI) at 15 μm isotropic resolution which is ~ 1,000 times higher than that of most preclinical MRI. Connectomes are generated with superresolution tract density images of ~5 μm. Brains are cleared, stained for selected proteins, and imaged by LSM at 1.8 μm/pixel. LSM data are registered into the reference MRH space with labels derived from the ABA common coordinate framework. The result is a high-dimensional integrated volume with registration (HiDiver) with alignment precision better than 50 µm. Throughput is sufficiently high that HiDiver is being used in quantitative studies of the impact of gene variants and aging on mouse brain cytoarchitecture and connectomics.

Artificial intelligence for diffusion MRI-based tissue microstructure estimation in the human brain: an overview

Artificial intelligence (AI) has made significant advances in the field of diffusion magnetic resonance imaging (dMRI) and other neuroimaging modalities. These techniques have been applied to various areas such as image reconstruction, denoising, detecting and removing artifacts, segmentation, tissue microstructure modeling, brain connectivity analysis, and diagnosis support. State-of-the-art AI algorithms have the potential to leverage optimization techniques in dMRI to advance sensitivity and inference through biophysical models. While the use of AI in brain microstructures has the potential to revolutionize the way we study the brain and understand brain disorders, we need to be aware of the pitfalls and emerging best practices that can further advance this field. Additionally, since dMRI scans rely on sampling of the q-space geometry, it leaves room for creativity in data engineering in such a way that it maximizes the prior inference. Utilization of the inherent geometry has been shown to improve general inference quality and might be more reliable in identifying pathological differences. We acknowledge and classify AI-based approaches for dMRI using these unifying characteristics. This article also highlighted and reviewed general practices and pitfalls involving tissue microstructure estimation through data-driven techniques and provided directions for building on them.

The influence of axonal beading and undulation on axonal diameter mapping

We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Larmor frequency shift from magnetized cylinders with arbitrary orientation distribution

The magnetic susceptibility of tissue can provide valuable information about its chemical composition and microstructural organization. However, the relation between the magnetic microstructure and the measurable Larmor frequency shift is understood only for a few idealized cases. Here we analyze the microstructure formed by magnetized, NMR-invisible infinite cylinders suspended in an NMR-reporting fluid. Through simulations, we scrutinize various geometries of mesoscopic Lorentz cavities and inclusions, and show that the cavity size should be approximately one order of magnitude larger than the width of the inclusions. We also analytically derive the Larmor frequency shift for a population of cylinders with arbitrary orientation dispersion and show that it is determined by the l = 2 Laplace expansion coefficients p 2 m of the cylinders' orientation distribution function. Our work underscores the need to account for microstructural organization when estimating magnetic tissue properties.

Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI

We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.

Harmonic viscoelastic response of 3D histology-informed white matter model

White matter (WM) consists of bundles of long axons embedded in a glial matrix, which lead to anisotropic mechanical properties of brain tissue, and this complicates direct numerical simulations of WM viscoelastic response. The detailed axonal geometry contains scales that range from axonal diameter (microscale) to many diameters (mesoscale) imposing an additional challenge to numerical simulations. Here we describe the development of a 3D homogenization model for the central nervous system (CNS) that accounts for the anisotropy introduced by the axon/neuroglia composite, the axonal trace curvature, and the tissue dynamic response in the frequency domain. Homogenized models that allow the incorporation of all the above factors are important for accurately simulating the tissue's mechanical behavior, and this in turn is essential in interpreting non-invasive elastography measurements. Geometric and material parameters affect the material properties and thus the response of the brain tissue. More complex, orthotropic, or anisotropic material properties are to be considered as necessitated by the 3D tissue structure. An assembly of micro-scale 3D representative elemental volumes (REVs) is constructed, leading to an integrated mesoscale WM finite element model. Assemblies of microscopic REVs, with orientations based on geometrical reconstructions driven by confocal microscopy data are employed to form the elements of the WM model. Each REV carries local material properties based on a finite element model of biphasic (axon-glial matrix) unidirectional composite. The viscoelastic response of the microscopic REVs is extracted based on geometric information and fiber volume fractions calculated from the relative distance between the local elements and global axonal trace. The response of the WM tissue model is homogenized by averaging the shear moduli over the total volume (thus deriving effective properties) under realistic external loading conditions. Under harmonic shear loading, it is proven that that the effective transverse shear moduli are higher than the axial moduli when the axon moduli are higher than the glial. Methodologically, the process of using micro-scale 3D REVs to describe more complex axon geometries avoids the partition process in traditional composite finite element methods (based on partition of finite element grids) and constitutes a robust algorithm to automatically build a WM model based on available axonal trace information. Analytically, the model provides unmatched simulation flexibility and computational power as the position, orientation, and the magnitude of each tissue building block is calculated using real tissue data, as are the training and testing processes at each level of the multiscale WM tissue.

Viscoelastic damage evaluation of the axon

In this manuscript, we have studied the microstructure of the axonal cytoskeleton and adopted a bottom-up approach to evaluate the mechanical responses of axons. The cytoskeleton of the axon includes the microtubules (MT), Tau proteins (Tau), neurofilaments (NF), and microfilaments (MF). Although most of the rigidity of the axons is due to the MT, the viscoelastic response of axons comes from the Tau. Early studies have shown that NF and MF do not provide significant elasticity to the overall response of axons. Therefore, the most critical aspect of the mechanical response of axons is the microstructural topology of how MT and Tau are connected and construct the cross-linked network. Using a scanning electron microscope (SEM), the cross-sectional view of the axons revealed that the MTs are organized in a hexagonal array and cross-linked by Tau. Therefore, we have developed a hexagonal Representative Volume Element (RVE) of the axonal microstructure with MT and Tau as fibers. The matrix of the RVE is modeled by considering a combined effect of NF and MF. A parametric study is done by varying fiber geometric and mechanical properties. The Young’s modulus and spacing of MT are varied between 1.5 and 1.9 GPa and 20–38 nm, respectively. Tau is modeled as a 3-parameter General Maxwell viscoelastic material. The failure strains for MT and Tau are taken to be 50 and 40%, respectively. A total of 4 RVEs are prepared for finite element analysis, and six loading cases are inspected to quantify the three-dimensional (3D) viscoelastic relaxation response. The volume-averaged stress and strain are then used to fit the relaxation Prony series. Next, we imposed varying strain rates (between 10/sec to 50/sec) on the RVE and analyzed the axonal failure process. We have observed that the 40% failure strain of Tau is achieved in all strain rates before the MT reaches its failure strain of 50%. The corresponding axonal failure strain and stress vary between 6 and 11% and 5–19.8 MPa, respectively. This study can be used to model macroscale axonal aggregate typical of the white matter region of the brain tissue.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

Orientation dependence of inhomogeneous magnetization transfer and dipolar order relaxation rate in phospholipid bilayers

Inhomogeneous magnetization transfer (ihMT) is a novel MRI technique used to measure white matter myelination in the brain and spinal cord. In the brain, ihMT has a strong orientation dependence which is likely to arise from the anisotropy of dipolar couplings between protons on oriented lipids in the myelin bilayers. We measured the orientation dependence of the second moment (M₂) of the lineshape, dipolar order relaxation rate (R_1D), and ihMT ratio (ihMTR) in an oriented phospholipid bilayer at 9.4 T. We found a strong orientation dependence in all three parameters. ihMTR and R_1D were maximized when the bilayers were aligned perpendicular to B₀ and minimized near the magic angle (∼54.7°). M₂ followed an orientation dependence given by the second Legendre polynomial squared as predicted by the form of the secular dipolar Hamiltonian. These results were used to calculate the orientation dependence of R_1D and ihMTR in a diffusionless myelin sheath model, which showed ihMTR was maximised for fibers perpendicular to B₀ and minimised at 45°, similar to ex-vivo spinal cord with a larger prepulse frequency offset, but in contrast to in vivo brain findings. Adding fiber dispersion to this model smoothed the orientation dependence curve as expected. Our results suggest the importance of the effects of lipid diffusion and prepulse offset frequency on ihMTR.

Resolution and b value dependent Structural Connectome in ex vivo Mouse Brain

Diffusion magnetic resonance imaging has been widely used in both clinical and preclinical studies to characterize tissue microstructure and structural connectivity. The diffusion MRI protocol for the Human Connectome Project (HCP) has been developed and optimized to obtain high-quality, high-resolution diffusion MRI (dMRI) datasets. However, such efforts have not been fully explored in preclinical studies, especially for rodents. In this study, high quality dMRI datasets of mouse brains were acquired at 9.4T system from two vendors. In particular, we acquired a high-spatial resolution dMRI dataset (25 μm isotropic with 126 diffusion encoding directions), which we believe to be the highest spatial resolution yet obtained; and a high-angular resolution dMRI dataset (50 μm isotropic with 384 diffusion encoding directions), which we believe to be the highest angular resolution compared to the dMRI datasets at the microscopic resolution. We systematically investigated the effects of three important parameters that affect the final outcome of the connectome: b value (1000s/mm² to 8000 s/mm²), angular resolution (10 to 126), and spatial resolution (25 μm to 200 μm). The stability of tractography and connectome increase with the angular resolution, where more than 50 angles is necessary to achieve consistent results. The connectome and quantitative parameters derived from graph theory exhibit a linear relationship to the b value (R² > 0.99); a single-shell acquisition with b value of 3000 s/mm² shows comparable results to the multi-shell high angular resolution dataset. The dice coefficient decreases and both false positive rate and false negative rate gradually increase with coarser spatial resolution. Our study provides guidelines and foundations for exploration of tradeoffs among acquisition parameters for the structural connectome in ex vivo mouse brain.

Bundle myelin fraction (BMF) mapping of different white matter connections using microstructure informed tractography

To date, we have scarce information about the relative myelination level of different fiber bundles in the human brain. Indirect evidence comes from postmortem histology data but histological stainings are unable to follow a specific bundle and determine its intrinsic myelination. In this context, quantitative MRI, and diffusion MRI tractography may offer a viable solution by providing, respectively, voxel-wise myelin sensitive maps and the pathways of the major tracts of the brain. Then, "tractometry" can be used to combine these two pieces of information by averaging tissue features (obtained from any voxel-wise map) along the streamlines recovered with diffusion tractography. Although this method has been widely used in the literature, in cases of voxels containing multiple fiber populations (each with different levels of myelination), tractometry provides biased results because the same value will be attributed to any bundle passing through the voxel. To overcome this bias, we propose a new method - named "myelin streamline decomposition" (MySD) - which extends convex optimization modeling for microstructure informed tractography (COMMIT) allowing the actual value measured by a microstructural map to be deconvolved on each individual streamline, thereby recovering unique bundle-specific myelin fractions (BMFs). We demonstrate the advantage of our method with respect to tractometry in well-studied bundles and compare the cortical projection of the obtained bundle-wise myelin values of both methods. We also prove the stability of our approach across different subjects and different MRI sensitive myelin mapping approaches. This work provides a proof-of-concept of in vivo investigations of entire neuronal pathways that, to date, are not possible.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (K_iso), variance due to diffusion anisotropy (K_aniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible μK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of μK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that μK significantly contributes to total diffusional kurtosis both in gray and white matter tissue but, as expected, not in the ventricles. The first μK maps of the human brain are presented, revealing that the spatial distribution of μK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring μK and assuming the multiple gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

Effect of Particle Size and Surface Charge on Nanoparticles Diffusion in the Brain White Matter

Purpose

Brain disorders have become a serious problem for healthcare worldwide. Nanoparticle-based drugs are one of the emerging therapies and have shown great promise to treat brain diseases. Modifications on particle size and surface charge are two efficient ways to increase the transport efficiency of nanoparticles through brain-blood barrier; however, partly due to the high complexity of brain microstructure and limited visibility of Nanoparticles (NPs), our understanding of how these two modifications can affect the transport of NPs in the brain is insufficient.

Methods

In this study, a framework, which contains a stochastic geometric model of brain white matter (WM) and a mathematical particle tracing model, was developed to investigate the relationship between particle size/surface charge of the NPs and their effective diffusion coefficients (D) in WM.

Results

The predictive capabilities of this method have been validated using published experimental tests. For negatively charged NPs, both particle size and surface charge are positively correlated with D before reaching a size threshold. When Zeta potential (Zp) is less negative than -10 mV, the difference between NPs’ D in WM and pure interstitial fluid (IF) is limited.

Conclusion

A deeper understanding on the relationships between particle size/surface charge of NPs and their D in WM has been obtained. The results from this study and the developed modelling framework provide important tools for the development of nano-drugs and nano-carriers to cure brain diseases.

What kind of network is the brain?

The different areas of the cerebral cortex are linked by a network of white matter, comprising the myelinated axons of pyramidal cells. Is this network a neural net, in the sense that representations of the world are embodied in the structure of the net, its pattern of nodes, and connections? Or is it a communications network, where the same physical substrate carries different information from moment to moment? This question is part of the larger question of whether the brain is better modeled by connectionism or by symbolic artificial intelligence (AI), but we review it in the specific context of the psychophysics of stimulus comparison and the format and protocol of information transmission over the long-range tracts of the brain.

Advanced Diffusion MR Imaging for Multiple Sclerosis in the Brain and Spinal Cord

Diffusion tensor imaging (DTI) has been established its usefulness in evaluating normal-appearing white matter (NAWM) and other lesions that are difficult to evaluate with routine clinical MRI in the evaluation of the brain and spinal cord lesions in multiple sclerosis (MS), a demyelinating disease. With the recent advances in the software and hardware of MRI systems, increasingly complex and sophisticated MRI and analysis methods, such as q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, white matter tract integrity, and multiple diffusion encoding, referred to as advanced diffusion MRI, have been proposed. These are capable of capturing in vivo microstructural changes in the brain and spinal cord in normal and pathological states in greater detail than DTI.

This paper reviews the current status of recent advanced diffusion MRI for assessing MS in vivo as part of an issue celebrating two decades of magnetic resonance in medical sciences (MRMS), an official journal of the Japanese Society of Magnetic Resonance in Medicine.

Towards a representative reference for MRI-based human axon radius assessment using light microscopy

Non-invasive assessment of axon radii via MRI bears great potential for clinical and neuroscience research as it is a main determinant of the neuronal conduction velocity. However, there is a lack of representative histological reference data at the scale of the cross-section of MRI voxels for validating the MRI-visible, effective radius (r_eff). Because the current gold standard stems from neuroanatomical studies designed to estimate the bulk-determined arithmetic mean radius (r_arith) on small ensembles of axons, it is unsuited to estimate the tail-weighted r_eff. We propose CNN-based segmentation on high-resolution, large-scale light microscopy (lsLM) data to generate a representative reference for r_eff. In a human corpus callosum, we assessed estimation accuracy and bias of r_arith and r_eff. Furthermore, we investigated whether mapping anatomy-related variation of r_arith and r_eff is confounded by low-frequency variation of the image intensity, e.g., due to staining heterogeneity. Finally, we analyzed the error due to outstandingly large axons in r_eff. Compared to r_arith, r_eff was estimated with higher accuracy (maximum normalized-root-mean-square-error of r_eff: 8.5 %; r_arith: 19.5 %) and lower bias (maximum absolute normalized-mean-bias-error of r_eff: 4.8 %; r_arith: 13.4 %). While r_arith was confounded by variation of the image intensity, variation of r_eff seemed anatomy-related. The largest axons contributed between 0.8 % and 2.9 % to r_eff. In conclusion, the proposed method is a step towards representatively estimating r_eff at MRI voxel resolution. Further investigations are required to assess generalization to other brains and brain areas with different axon radii distributions.

Dysregulated miR-29a-3p/PMP22 Modulates Schwann Cell Proliferation and Migration During Peripheral Nerve Regeneration

Schwann cells switch to a repair phenotype following peripheral nerve injury and create a favorable microenvironment to drive nerve repair. Many microRNAs (miRNAs) are differentially expressed in the injured peripheral nerves and play essential roles in regulating Schwann cell behaviors. Here, we examine the temporal expression patterns of miR-29a-3p after peripheral nerve injury and demonstrate significant up-regulation of miR-29a-3p in injured sciatic nerves. Elevated miR-29a-3p inhibits Schwann cell proliferation and migration, while suppressed miR-29a-3p executes reverse effects. In vivo injection of miR-29a-3p agomir to rat sciatic nerves hinders the proliferation and migration of Schwann cells, delays the elongation and myelination of axons, and retards the functional recovery of injured nerves. Mechanistically, miR-29a-3p modulates Schwann cell activities via negatively regulating peripheral myelin protein 22 (PMP22), and PMP22 extensively affects Schwann cell metabolism. Our results disclose the vital role of miR-29a-3p/PMP22 in regulating Schwann cell phenotype following sciatic nerve injury and shed light on the mechanistic basis of peripheral nerve regeneration.

Does powder averaging remove dispersion bias in diffusion MRI diameter estimates within real 3D axonal architectures?

Noninvasive estimation of axon diameter with diffusion MRI holds the potential to investigate the dynamic properties of the brain network and pathology of neurodegenerative diseases. Recent studies use powder averaging to account for complex white matter architectures, but these have not been validated for real axonal geometries from regions that contain fibre crossings. Here, we present 120-304μm long segmented axons from X-ray nano-holotomography volumes of a splenium and crossing fibre region of a vervet monkey brain. We show that the axons in the complex crossing fibre region, which contains callosal, association, and corticospinal connections, are larger and exhibit a wider distribution than those of the splenium region. To accurately estimate the axon diameter in these regions, therefore, sensitivity to a wide range of diameters is required. We demonstrate how the q-value, b-value, signal-to-noise ratio and the assumed intra-axonal parallel diffusivity influence the range of measurable diameters with powder average approaches. Furthermore, we show how Gaussian distributed noise results in a wider range of measurable diameter at high b-values than Rician distributed noise, even at high signal-to-noise ratios of 100. The number of gradient directions is also shown to impose a lower bound on measurable diameter. Our results indicate that axon diameter estimation can be performed with only few b-shells, and that additional shells do not improve the accuracy of the estimate. For strong gradients available on human Connectom and preclinical scanners, Monte Carlo simulations of diffusion confirm that powder averaging techniques succeed in providing accurate estimates of axon diameter across a range of sequence parameters and diffusion times, even in complex white matter architectures. At relatively low b-values, the diameter estimate becomes sensitive to axonal microdispersion and the intra-axonal parallel diffusivity shows time dependency at both in vivo and ex vivo intrinsic diffusivities.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Comparative analysis of signal models for microscopic fractional anisotropy estimation using q-space trajectory encoding

Microscopic diffusion anisotropy imaging using diffusion-weighted MRI and multidimensional diffusion encoding is a promising method for quantifying clinically and scientifically relevant microstructural properties of neural tissue. Several methods for estimating microscopic fractional anisotropy (µFA), a normalized measure of microscopic diffusion anisotropy, have been introduced but the differences between the methods have received little attention thus far. In this study, the accuracy and precision of µFA estimation using q-space trajectory encoding and different signal models were assessed using imaging experiments and simulations. Three healthy volunteers and a microfibre phantom were imaged with five non-zero b-values and gradient waveforms encoding linear and spherical b-tensors. Since the ground-truth µFA was unknown in the imaging experiments, Monte Carlo random walk simulations were performed using axon-mimicking fibres for which the ground truth was known. Furthermore, parameter bias due to time-dependent diffusion was quantified by repeating the simulations with tuned waveforms, which have similar power spectra, and with triple diffusion encoding, which, unlike q-space trajectory encoding, is not based on the assumption of time-independent diffusion. The truncated cumulant expansion of the powder-averaged signal, gamma-distributed diffusivities assumption, and q-space trajectory imaging, a generalization of the truncated cumulant expansion to individual signals, were used to estimate µFA. The gamma-distributed diffusivities assumption consistently resulted in greater µFA values than the second order cumulant expansion, 0.1 greater when averaged over the whole brain. In the simulations, the generalized cumulant expansion provided the most accurate estimates. Importantly, although time-dependent diffusion caused significant overestimation of µFA using all the studied methods, the simulations suggest that the resulting bias in µFA is less than 0.1 in human white matter.

Scan-rescan repeatability of axonal imaging metrics using high-gradient diffusion MRI and statistical implications for study design

Axon diameter mapping using diffusion MRI in the living human brain has attracted growing interests with the increasing availability of high gradient strength MRI systems. A systematic assessment of the consistency of axon diameter estimates within and between individuals is needed to gain a comprehensive understanding of how such methods extend to quantifying differences in axon diameter index between groups and facilitate the design of neurobiological studies using such measures. We examined the scan-rescan repeatability of axon diameter index estimation based on the spherical mean technique (SMT) approach using diffusion MRI data acquired with gradient strengths up to 300 mT/m on a 3T Connectom system in 7 healthy volunteers. We performed statistical power analyses using data acquired with the same protocol in a larger cohort consisting of 15 healthy adults to investigate the implications for study design. Results revealed a high degree of repeatability in voxel-wise restricted volume fraction estimates and tract-wise estimates of axon diameter index derived from high-gradient diffusion MRI data. On the region of interest (ROI) level, across white matter tracts in the whole brain, the Pearson’s correlation coefficient of the axon diameter index estimated between scan and rescan experiments was r = 0.72 with an absolute deviation of 0.18 μm. For an anticipated 10% effect size in studies of axon diameter index, most white matter regions required a sample size of less than 15 people to observe a measurable difference between groups using an ROI-based approach. To facilitate the use of high-gradient strength diffusion MRI data for neuroscientific studies of axonal microstructure, the comprehensive multi-gradient strength, multi-diffusion time data used in this work will be made publicly available, in support of open science and increasing the accessibility of such data to the greater scientific community.

The link between gliomas infiltration and white matter architecture investigated with electron microscopy and diffusion tensor imaging

Diffuse low-grade gliomas (DLGG) display different preferential locations in eloquent and secondary associative brain areas. The reason for this tendency is still unknown. We hypothesized that the intrinsic architecture and water diffusion properties of the white matter bundles in these regions may facilitate gliomas infiltration. Magnetic resonance imaging of sixty-seven diffuse low-grade gliomas patients were normalized to/and segmented in MNI space to create three probabilistic infiltration weighted gradient maps according to the molecular status of each tumor group (IDH mutated, IDH wild-type and IDH mutated/1p19q co-deleted). Diffusion tensor imaging (DTI)- based parameters were derived for five major white matter bundles, displaying regional differences in the grade of infiltration, averaged over 20 healthy individuals acquired from the Human connectome project (HCP) database. Transmission electron microscopy (TEM) was used to analyze fiber density, fiber diameter and g-ratio in 100 human white matter regions, sampled from cadaver specimens, reflecting areas with different gliomas infiltration in each white matter bundle. Histological results and DTI-based parameters were compared in anatomical regions of high- and low grade of infiltration (HIF and LIF) respectively. We detected differences in the white matter infiltration of five major white matter bundles in three groups. Astrocytomas IDHm infiltrated left fronto-temporal subcortical areas. Astrocytomas IDHwt were detected in the posterior-temporal and temporo-parietal regions bilaterally. Oligodendrogliomas IDHm/1p19q infiltrated anterior subcortical regions of the frontal lobes bilaterally. Regional differences within the same white matter bundles were detected by both TEM- and DTI analysis linked to different topographical variables. Our multimodal analysis showed that HIF regions, common to all the groups, displayed a smaller fiber diameter, lower FA and higher RD compared with LIF regions. Our results suggest that the both morphological features and diffusion parameters of the white matter may be different in regions linked to the preferential location of DLGG.

Large-scale electron microscopic volume imaging of interfascicular oligodendrocytes in the mouse corpus callosum

Single oligodendrocytes produce myelin sheaths around multiple axons in the central nervous system. Interfascicular oligodendrocytes (IOs) facilitate nerve conduction, but their detailed morphologies remain largely unknown. In the present study, we three-dimensionally reconstructed IOs in the corpus callosum of adult mouse using serial block face scanning electron microscopy. The cell bodies of IOs were morphologically polarized and extended thick processes from the cytoplasm-rich part of the cell. Processes originating from the cell body of each IO can be classified into two types: one myelinates an axon without branching, while the other type branches and each branch myelinates a distinct axon. Myelin sheaths originating from a particular IO have biased thicknesses, wrapping axons of a limited range of diameters. Consistent with this finding, IOs transduced and visualized with a rabies viral vector expressing GFP showed statistically significant variation in their myelination patterns. We further reconstructed the sheath immediately adjacent to that derived from each of the analyzed IOs; the thicknesses of the pair of sheaths were significantly correlated despite emanating from different IOs. These results suggest that a single axon could regulate myelin sheath thicknesses, even if the sheaths are derived from distinct IOs. Collectively, our results indicate that the IOs have their own myelin profiles defined by myelin thickness and axonal diameter although axons may regulate thickness of myelin sheath.

Bundle-Specific Axon Diameter Index as a New Contrast to Differentiate White Matter Tracts

In the central nervous system of primates, several pathways are characterized by different spectra of axon diameters. In vivo methods, based on diffusion-weighted magnetic resonance imaging, can provide axon diameter index estimates non-invasively. However, such methods report voxel-wise estimates, which vary from voxel-to-voxel for the same white matter bundle due to partial volume contributions from other pathways having different microstructure properties. Here, we propose a novel microstructure-informed tractography approach, COMMITAxSize, to resolve axon diameter index estimates at the streamline level, thus making the estimates invariant along trajectories. Compared to previously proposed voxel-wise methods, our formulation allows the estimation of a distinct axon diameter index value for each streamline, directly, furnishing a complementary measure to the existing calculation of the mean value along the bundle. We demonstrate the favourable performance of our approach comparing our estimates with existing histologically-derived measurements performed in the corpus callosum and the posterior limb of the internal capsule. Overall, our method provides a more robust estimation of the axon diameter index of pathways by jointly estimating the microstructure properties of the tissue and the macroscopic organisation of the white matter connectivity.

Activated microglia drive demyelination via CSF1R signaling

Microgliosis is a prominent pathological feature in many neurological diseases including multiple sclerosis (MS), a progressive auto-immune demyelinating disorder. The precise role of microglia, parenchymal central nervous system (CNS) macrophages, during demyelination, and the relative contributions of peripheral macrophages are incompletely understood. Classical markers used to identify microglia do not reliably discriminate between microglia and peripheral macrophages, confounding analyses. Here, we use a genetic fate mapping strategy to identify microglia as predominant responders and key effectors of demyelination in the cuprizone (CUP) model. Colony-stimulating factor 1 (CSF1), also known as macrophage colony-stimulating factor (M-CSF) - a secreted cytokine that regulates microglia development and survival-is upregulated in demyelinated white matter lesions. Depletion of microglia with the CSF1R inhibitor PLX3397 greatly abrogates the demyelination, loss of oligodendrocytes, and reactive astrocytosis that results from CUP treatment. Electron microscopy (EM) and serial block face imaging show myelin sheaths remain intact in CUP treated mice depleted of microglia. However, these CUP-damaged myelin sheaths are lost and robustly phagocytosed upon-repopulation of microglia. Direct injection of CSF1 into CNS white matter induces focal microgliosis and demyelination indicating active CSF1 signaling can promote demyelination. Finally, mice defective in adopting a toxic astrocyte phenotype that is driven by microglia nevertheless demyelinate normally upon CUP treatment implicating microglia rather than astrocytes as the primary drivers of CUP-mediated demyelination. Together, these studies indicate activated microglia are required for and can drive demyelination directly and implicate CSF1 signaling in these events.

Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding

Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.

DeepACSON automated segmentation of white matter in 3D electron microscopy

Tracing the entirety of ultrastructures in large three-dimensional electron microscopy (3D-EM) images of the brain tissue requires automated segmentation techniques. Current segmentation techniques use deep convolutional neural networks (DCNNs) and rely on high-contrast cellular membranes and high-resolution EM volumes. On the other hand, segmenting low-resolution, large EM volumes requires methods to account for severe membrane discontinuities inescapable. Therefore, we developed DeepACSON, which performs DCNN-based semantic segmentation and shape-decomposition-based instance segmentation. DeepACSON instance segmentation uses the tubularity of myelinated axons and decomposes under-segmented myelinated axons into their constituent axons. We applied DeepACSON to ten EM volumes of rats after sham-operation or traumatic brain injury, segmenting hundreds of thousands of long-span myelinated axons, thousands of cell nuclei, and millions of mitochondria with excellent evaluation scores. DeepACSON quantified the morphology and spatial aspects of white matter ultrastructures, capturing nanoscopic morphological alterations five months after the injury.

A shape-adjusted ellipse approach corrects for varied axonal dispersion angles and myelination in primate nerve roots

Segmentation of axons in light and electron micrographs allows for quantitative high-resolution analysis of nervous tissues, but varied axonal dispersion angles result in over-estimates of fiber sizes. To overcome this technical challenge, we developed a novel shape-adjusted ellipse (SAE) determination of axonal size and myelination as an all-inclusive and non-biased tool to correct for oblique nerve fiber presentations. Our new resource was validated by light and electron microscopy against traditional methods of determining nerve fiber size and myelination in rhesus macaques as a model system. We performed detailed segmental mapping and characterized the morphological signatures of autonomic and motor fibers in primate lumbosacral ventral roots (VRs). An en bloc inter-subject variability for the preganglionic parasympathetic fibers within the L7-S2 VRs was determined. The SAE approach allows for morphological ground truth data collection and assignment of individual axons to functional phenotypes with direct implications for fiber mapping and neuromodulation studies.

Sensitivity analysis of effective transverse shear viscoelastic and diffusional properties of myelinated white matter

Motivated by the need to interpret the results from a combined use of in vivo brain Magnetic Resonance Elastography (MRE) and Diffusion Tensor Imaging (DTI), we developed a computational framework to study the sensitivity of single-frequency MRE and DTI metrics to white matter microstructure and cell-level mechanical and diffusional properties. White matter was modeled as a triphasic unidirectional composite, consisting of parallel cylindrical inclusions (axons) surrounded by sheaths (myelin), and embedded in a matrix (glial cells plus extracellular matrix). Only 2D mechanics and diffusion in the transverse plane (perpendicular to the axon direction) was considered, and homogenized (effective) properties were derived for a periodic domain containing a single axon. The numerical solutions of the MRE problem were performed with ABAQUS and by employing a sophisticated boundary-conforming grid generation scheme. Based on the linear viscoelastic response to harmonic shear excitation and steady-state diffusion in the transverse plane, a systematic sensitivity analysis of MRE metrics (effective transverse shear storage and loss moduli) and DTI metric (effective radial diffusivity) was performed for a wide range of microstructural and intrinsic (phase-based) physical properties. The microstructural properties considered were fiber volume fraction, and the myelin sheath/axon diameter ratio. The MRE and DTI metrics are very sensitive to the fiber volume fraction, and the intrinsic viscoelastic moduli of the glial phase. The MRE metrics are nonlinear functions of the fiber volume fraction, but the effective diffusion coefficient varies linearly with it. Finally, the transverse metrics of both MRE and DTI are insensitive to the axon diameter in steady state. Our results are consistent with the limited anisotropic MRE and co-registered DTI measurements, mainly in the corpus callosum, available in the literature. We conclude that isotropic MRE and DTI constitutive models are good approximations for myelinated white matter in the transverse plane. The unidirectional composite model presented here is used for the first time to model harmonic shear stress under MRE-relevant frequency on the cell level. This model can be extended to 3D in order to inform the solution of the inverse problem in MRE, establish the biological basis of MRE metrics, and integrate MRE/DTI with other modalities towards increasing the specificity of neuroimaging.

Axon morphology is modulated by the local environment and impacts the noninvasive investigation of its structure-function relationship

Axonal conduction velocity, which ensures efficient function of the brain network, is related to axon diameter. Noninvasive, in vivo axon diameter estimates can be made with diffusion magnetic resonance imaging, but the technique requires three-dimensional (3D) validation. Here, high-resolution, 3D synchrotron X-ray nano-holotomography images of white matter samples from the corpus callosum of a monkey brain reveal that blood vessels, cells, and vacuoles affect axonal diameter and trajectory. Within single axons, we find that the variation in diameter and conduction velocity correlates with the mean diameter, contesting the value of precise diameter determination in larger axons. These complex 3D axon morphologies drive previously reported 2D trends in axon diameter and g-ratio. Furthermore, we find that these morphologies bias the estimates of axon diameter with diffusion magnetic resonance imaging and, ultimately, impact the investigation and formulation of the axon structure-function relationship.

A simple estimate of axon size with diffusion MRI

Noninvasive estimation of mean axon diameter presents a new opportunity to explore white matter plasticity, development, and pathology. Several diffusion-weighted MRI (DW-MRI) methods have been proposed to measure the average axon diameter in white matter, but they typically require many diffusion encoding measurements and complicated mathematical models to fit the signal to multiple tissue compartments, including intra- and extra-axonal spaces. Here, Monte Carlo simulations uncovered a straightforward DW-MRI metric of axon diameter: the change in radial apparent diffusion coefficient estimated at different effective diffusion times, ΔD_⊥. Simulations indicated that this metric increases monotonically within a relevant range of effective mean axon diameter while being insensitive to changes in extra-axonal volume fraction, axon diameter distribution, g-ratio, and influence of myelin water. Also, a monotonic relationship was found to exist for signals coming from both intra- and extra-axonal compartments. The slope in ΔD_⊥ with effective axon diameter increased with the difference in diffusion time of both oscillating and pulsed gradient diffusion sequences.