Imaging brain microstructure with diffusion MRI: practicality and applications

SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI

The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called 'b-tensor' encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.

Diffusion tensor imaging reveals neuronal microstructural changes in myalgic encephalomyelitis/chronic fatigue syndrome

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) patients suffer from a variety of physical and neurological complaints indicating the central nervous system plays a role in ME/CFS pathophysiology. Diffusion tensor imaging (DTI) has been used to study microstructural changes in neurodegenerative diseases. In this study, we evaluated DTI parameters to investigate microstructural abnormalities in ME/CFS patients. We estimated DTI parameters in 25 ME/CFS patients who met Fukuda criteria (ME/CFS_Fukuda ), 18 ME/CFS patients who met International Consensus Criteria (ICC) (ME/CFS_ICC ) only and 26 healthy control (HC) subjects. In addition to voxel-based DTI-parameter group comparisons, we performed voxel-based DTI-parameter interaction-with-group regressions with clinical and autonomic measures to test for abnormal regressions. Group comparisons between ME/CFS_ICC and HC detected significant clusters (a) with decreased axial diffusivity (p = .001) and mean diffusivity (p = .01) in the descending cortico-cerebellar tract in the midbrain and pons and (b) with increased transverse diffusivity in the medulla. The mode of anisotropy was significantly decreased (p = .001) in a cluster in the superior longitudinal fasciculus region. Voxel-based group comparisons between ME/CFS_Fukuda and HC did not detect significant clusters. For ME/CFS_ICC and HC, DTI parameter interaction-with-group regressions were abnormal for the clinical measures of information processing score, SF36 physical, sleep disturbance score and respiration rate in both grey and white matter regions. Our study demonstrated that DTI parameters are sensitive to microstructural changes in ME/CFS_ICC and could potentially act as an imaging biomarker of abnormal pathophysiology in ME/CFS. The study also shows that strict case definitions are essential in investigation of the pathophysiology of ME/CFS.

Neurobiological underpinnings of rapid white matter plasticity during intensive reading instruction

Diffusion MRI is a powerful tool for imaging brain structure, but it is challenging to discern the biological underpinnings of plasticity inferred from these and other non-invasive MR measurements. Biophysical modeling of the diffusion signal aims to render a more biologically rich image of tissue microstructure, but the application of these models comes with important caveats. A separate approach for gaining biological specificity has been to seek converging evidence from multi-modal datasets. Here we use metrics derived from diffusion kurtosis imaging (DKI) and the white matter tract integrity (WMTI) model along with quantitative MRI measurements of T1 relaxation to characterize changes throughout the white matter during an 8-week, intensive reading intervention (160 total hours of instruction). Behavioral measures, multi-shell diffusion MRI data, and quantitative T1 data were collected at regular intervals during the intervention in a group of 33 children with reading difficulties (7-12 years old), and over the same period in an age-matched non-intervention control group. Throughout the white matter, mean 'extra-axonal' diffusivity was inversely related to intervention time. In contrast, model estimated axonal water fraction (AWF), overall diffusion kurtosis, and T1 relaxation time showed no significant change over the intervention period. Both diffusion and quantitative T1 based metrics were correlated with pre-intervention reading performance, albeit with distinct anatomical distributions. These results are consistent with the view that rapid changes in diffusion properties reflect phenomena other than widespread changes in myelin density. We discuss this result in light of recent work highlighting non-axonal factors in experience-dependent plasticity and learning.

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.

Design and characterization of a 3D-printed axon-mimetic phantom for diffusion MRI

PURPOSE

To introduce and characterize inexpensive and easily produced 3D-printed axon-mimetic diffusion MRI phantoms in terms of pore geometry and diffusion kurtosis imaging metrics.

METHODS

Phantoms were 3D-printed with a composite printing material that, after the dissolution of the polyvinyl alcohol, exhibits microscopic fibrous pores. Confocal microscopy and synchrotron phase-contrast micro-CT imaging were performed to visualize and assess the pore sizes. Diffusion MRI scans of four identical phantoms and phantoms with varying print parameters in water were performed at 9.4 T. Diffusion kurtosis imaging was fit to both data sets and used to assess the reproducibility between phantoms and effects of print parameters on diffusion kurtosis imaging metrics. Identical scans were performed 25 and 76 days later, to test their stability.

RESULTS

Segmentation of pores in three microscopy images yielded a mean, median, and SD of equivalent pore diameters of 7.57 μm, 3.51 μm, and 12.13 μm, respectively. Phantoms had T₁ /T₂ = 2 seconds/180 ms, and those with identical parameters showed a low coefficient of variation (~10%) in mean diffusivity (1.38 × 10^-3 mm² /s) and kurtosis (0.52) metrics and radial diffusivity (1.01 × 10^-3 mm² /s) and kurtosis (1.13) metrics. Printing temperature and speed had a small effect on diffusion kurtosis imaging metrics (< 16%), whereas infill density had a larger and more variable effect (> 16%). The stability analysis showed small changes over 2.5 months (< 7%).

CONCLUSION

Three-dimension-printed axon-mimetic phantoms can mimic the fibrous structure of axon bundles on a microscopic scale, serving as complex, anisotropic diffusion MRI phantoms.

On the potential for mapping apparent neural soma density via a clinically viable diffusion MRI protocol

•

B-tensor encoding enables estimation of spherical cellular structures in the brain.

•

Spherical compartments may provide markers for apparent neural soma density.

•

Model parameters can be estimated in a fast and robust way using deep learning.

•

Practical acquisition times are achievable on widely available clinical scanners.

Diffusion MRI is a valuable tool for probing tissue microstructure in the brain noninvasively. Today, model-based techniques are widely available and used for white matter characterisation where their development is relatively mature. Conversely, tissue modelling in grey matter is more challenging, and no generally accepted models exist. With advances in measurement technology and modelling efforts, a clinically viable technique that reveals salient features of grey matter microstructure, such as the density of quasi-spherical cell bodies and quasi-cylindrical cell projections, is an exciting prospect. As a step towards capturing the microscopic architecture of grey matter in clinically feasible settings, this work uses a biophysical model that is designed to disentangle the diffusion signatures of spherical and cylindrical structures in the presence of orientation heterogeneity, and takes advantage of B-tensor encoding measurements, which provide additional sensitivity compared to standard single diffusion encoding sequences. For the fast and robust estimation of microstructural parameters, we leverage recent advances in machine learning and replace conventional fitting techniques with an artificial neural network that fits complex biophysical models within seconds. Our results demonstrate apparent markers of spherical and cylindrical geometries in healthy human subjects, and in particular an increased volume fraction of spherical compartments in grey matter compared to white matter. We evaluate the extent to which spherical and cylindrical geometries may be interpreted as correlates of neural soma and neural projections, respectively, and quantify parameter estimation errors in the presence of various departures from the modelling assumptions. While further work is necessary to translate the ideas presented in this work to the clinic, we suggest that biomarkers focussing on quasi-spherical cellular geometries may be valuable for the enhanced assessment of neurodevelopmental disorders and neurodegenerative diseases.

Coaxial electrospun biomimetic copolymer fibres for application in diffusion magnetic resonance imaging

Objective. The use of diffusion magnetic resonance imaging (dMRI) opens the door to characterizing brain microstructure because water diffusion is anisotropic in axonal fibres in brain white matter and is sensitive to tissue microstructural changes. As dMRI becomes more sophisticated and microstructurally informative, it has become increasingly important to use a reference object (usually called an imaging phantom) for validation of dMRI. This study aims to develop axon-mimicking physical phantoms from biocopolymers and assess their feasibility for validating dMRI measurements.Approach. We employed a simple and one-step method-coaxial electrospinning-to prepare axon-mimicking hollow microfibres from polycaprolactone-b-polyethylene glycol (PCL-b-PEG) and poly(D, L-lactide-co-glycolic) acid (PLGA), and used them as building elements to create axon-mimicking phantoms. Electrospinning was firstly conducted using two types of PCL-b-PEG and two types of PLGA with different molecular weights in various solvents, with different polymer concentrations, for determining their spinnability. Polymer/solvent concentration combinations with good fibre spinnability were used as the shell material in the following co-electrospinning process in which the polyethylene oxide polymer was used as the core material. Following the microstructural characterization of both electrospun and co-electrospun fibres using optical and electron microscopy, two prototype phantoms were constructed from co-electrospun anisotropic hollow microfibres after inserting them into water-filled test tubes.Main results. Hollow microfibres that mimic the axon microstructure were successfully prepared from the appropriate core and shell material combinations. dMRI measurements of two phantoms on a 7 tesla (T) pre-clinical scanner revealed that diffusivity and anisotropy measurements are in the range of brain white matter.Significance. This feasibility study showed that co-electrospun PCL-b-PEG and PLGA microfibre-based axon-mimicking phantoms could be used in the validation of dMRI methods which seek to characterize white matter microstructure.

Age-related hearing loss and cognitive decline: MRI and cellular evidence

Extensive evidence supports the association between age-related hearing loss (ARHL) and cognitive decline. It is, however, unknown whether a causal relationship exists between these two, or whether they both result from shared mechanisms. This paper intends to study this relationship through a comprehensive review of MRI findings as well as evidence of cellular alterations. Our review of structural MRI studies demonstrates that ARHL is independently linked to accelerated atrophy of total and regional brain volumes and reduced white matter integrity. Resting-state and task-based fMRI studies on ARHL also show changes in spontaneous neural activity and brain functional connectivity; and alterations in brain areas supporting auditory, language, cognitive, and affective processing independent of age, respectively. Although MRI findings support a causal relationship between ARHL and cognitive decline, the contribution of potential shared mechanisms should also be considered. In this regard, the review of cellular evidence indicates their role as possible common mechanisms underlying both age-related changes in hearing and cognition. Considering existing evidence, no single hypothesis can explain the link between ARHL and cognitive decline, and the contribution of both causal (i.e., the sensory hypothesis) and shared (i.e., the common cause hypothesis) mechanisms is expected.

Can MRI measure myelin? Systematic review, qualitative assessment, and meta-analysis of studies validating microstructural imaging with myelin histology

Recent years have seen an increased understanding of the importance of myelination in healthy brain function and neuropsychiatric diseases. Non-invasive microstructural magnetic resonance imaging (MRI) holds the potential to expand and translate these insights to basic and clinical human research, but the sensitivity and specificity of different MR markers to myelination is a subject of debate. To consolidate current knowledge on the topic, we perform a systematic review and meta-analysis of studies that validate microstructural imaging by combining it with myelin histology. We find meta-analytic evidence for correlations between various myelin histology metrics and markers from different MRI modalities, including fractional anisotropy, radial diffusivity, macromolecular pool, magnetization transfer ratio, susceptibility and longitudinal relaxation rate, but not mean diffusivity. Meta-analytic correlation effect sizes range widely, between R2 = 0.26 and R2 = 0.82. However, formal comparisons between MRI-based myelin markers are limited by methodological variability, inconsistent reporting and potential for publication bias, thus preventing the establishment of a single most sensitive strategy to measure myelin with MRI. To facilitate further progress, we provide a detailed characterisation of the evaluated studies as an online resource. We also share a set of 12 recommendations for future studies validating putative MR-based myelin markers and deploying them in vivo in humans.

In vivo tensor-valued diffusion MRI of focal demyelination in white and deep grey matter of rodents

•

We performed in-vivo tensor-valued diffusion MRI in demyelinating rodents.

•

Lysolecithin was injected in white and deep grey matter to cause focal demyelination.

•

Focal demyelination reduced microscopic fractional anisotropy (µFA).

•

Isotropic kurtosis may be particularly sensitive to deep grey matter lesions.

Background

Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease leading to damage of white matter (WM) and grey matter (GM). Magnetic resonance imaging (MRI) is the modality of choice to assess brain damage in MS, but there is an unmet need in MRI for achieving higher sensitivity and specificity to MS-related microstructural alterations in WM and GM.

Objective

To explore whether tensor-valued diffusion MRI (dMRI) can yield sensitive microstructural read-outs for focal demyelination in cerebral WM and deep GM (DGM).

Methods

Eight rats underwent L-α-Lysophosphatidylcholine (LPC) injections in the WM and striatum to introduce focal demyelination. Multimodal MRI was performed at 7 Tesla after 7 days. Tensor-valued dMRI was complemented by diffusion tensor imaging, quantitative MRI and proton magnetic resonance spectroscopy (MRS).

Results

Quantitative MRI and MRS confirmed that LPC injections caused inflammatory demyelinating lesions in WM and DGM. Tensor-valued dMRI illustrated a significant decline of microscopic fractional anisotropy (µFA) in both LPC-treated WM and DGM (P < 0.005) along with a marked increase of isotropic kurtosis (MK_I) in DGM (P < 0.0001).

Conclusion

Tensor-valued dMRI bears considerable potential for microstructural imaging in MS, suggesting a regional µFA decrease may be a sensitive indicator of MS lesions, while a regional MK_I increase may be particularly sensitive in detecting DGM lesions of MS.

Frontoinsular cortical microstructure is linked to life satisfaction in young adulthood

Life satisfaction is a component of subjective well-being that reflects a global judgement of the quality of life according to an individual's own needs and expectations. As a psychological construct, it has attracted attention due to its relationship to mental health, resilience to stress, and other factors. Neuroimaging studies have identified neurobiological correlates of life satisfaction; however, they are limited to functional connectivity and gray matter morphometry. We explored features of gray matter microstructure obtained through compartmental modeling of multi-shell diffusion MRI data, and we examined cortical microstructure in frontoinsular cortex in a cohort of 807 typical young adults scanned as part of the Human Connectome Project. Our experiments identified the orientation dispersion index (ODI), and analogously fractional anisotropy (FA), of frontoinsular cortex as a robust set of anatomically-specific lateralized diffusion MRI microstructure features that are linked to life satisfaction, independent of other demographic, socioeconomic, and behavioral factors. We further validated our findings in a secondary test-retest dataset and found high reliability of our imaging metrics and reproducibility of outcomes. In our analysis of twin and non-twin siblings, we found basic microstructure in frontoinsular cortex to be strongly genetically determined. We also found a more moderate but still very significant genetic role in determining microstructure as it relates to life satisfaction in frontoinsular cortex. Our findings suggest a potential linkage between well-being and microscopic features of frontoinsular cortex, which may reflect cellular morphology and architecture and may more broadly implicate the integrity of the homeostatic processing performed by frontoinsular cortex as an important component of an individual's judgements of life satisfaction.

Characterizing extracellular diffusion properties using diffusion-weighted MRS of sucrose injected in mouse brain

Brain water and some critically important energy metabolites, such as lactate or glucose, are present in both intracellular and extracellular spaces (ICS/ECS) at significant levels. This ubiquitous nature makes diffusion MRI/MRS data sometimes difficult to interpret and model. While it is possible to glean information on the diffusion properties in ICS by measuring the diffusion of purely intracellular endogenous metabolites (such as NAA), the absence of endogenous markers specific to ECS hampers similar analyses in this compartment. In past experiments, exogenous probes have therefore been injected into the brain to assess their apparent diffusion coefficient (ADC) and thus estimate tortuosity in ECS. Here, we use a similar approach in mice by injecting sucrose, a well-known ECS marker, in either the lateral ventricles or directly in the prefrontal cortex. For the first time, we propose a thorough characterization of ECS diffusion properties encompassing (1) short-range restriction by looking at signal attenuation at high b values, (2) tortuosity and long-range restriction by measuring ADC time-dependence at long diffusion times and (3) microscopic anisotropy by performing double diffusion encoding (DDE) measurements. Overall, sucrose diffusion behavior is strikingly different from that of intracellular metabolites. Acquisitions at high b values not only reveal faster sucrose diffusion but also some sensitivity to restriction, suggesting that the diffusion in ECS is not fully Gaussian at high b. The time evolution of the ADC at long diffusion times shows that the tortuosity regime is not reached yet in the case of sucrose, while DDE experiments suggest that it is not trapped in elongated structures. No major difference in sucrose diffusion properties is reported between the two investigated routes of injection and brain regions. These original experimental insights should be useful to better interpret and model the diffusion signal of molecules that are distributed between ICS and ECS compartments.

Inflammation-driven glial alterations in the cuprizone mouse model probed with diffusion-weighted magnetic resonance spectroscopy at 11.7 T

Inflammation of brain tissue is a complex response of the immune system to the presence of toxic compounds or to cell injury, leading to a cascade of pathological processes that include glial cell activation. Noninvasive MRI markers of glial reactivity would be very useful for in vivo detection and monitoring of inflammation processes in the brain, as well as for evaluating the efficacy of personalized treatments. Due to their specific location in glial cells, myo-inositol (mIns) and choline compounds (tCho) seem to be the best candidates for probing glial-specific intra-cellular compartments. However, their concentrations quantified using conventional proton MRS are not specific for inflammation. In contrast, it has been recently suggested that mIns intra-cellular diffusion, measured using diffusion-weighted MRS (DW-MRS) in a mouse model of reactive astrocytes, could be a specific marker of astrocytic hypertrophy. In order to evaluate the specificity of both mIns and tCho diffusion to inflammation-driven glial alterations, we performed DW-MRS in a volume of interest containing the corpus callosum and surrounding tissue of cuprizone-fed mice after 6 weeks of intoxication, and evaluated the extent of astrocytic and microglial alterations using immunohistochemistry. Both mIns and tCho apparent diffusion coefficients were significantly elevated in cuprizone-fed mice compared with control mice, and histologic evaluation confirmed the presence of severe inflammation. Additionally, mIns and tCho diffusion showed, respectively, strong and moderate correlations with histological measures of astrocytic and microglial area fractions, confirming DW-MRS as a promising tool for specific detection of glial changes under pathological conditions.

Decomposition of high-frequency electrical conductivity into extracellular and intracellular compartments based on two-compartment model using low-to-high multi-b diffusion MRI

Background

As an object’s electrical passive property, the electrical conductivity is proportional to the mobility and concentration of charged carriers that reflect the brain micro-structures. The measured multi-b diffusion-weighted imaging (Mb-DWI) data by controlling the degree of applied diffusion weights can quantify the apparent mobility of water molecules within biological tissues. Without any external electrical stimulation, magnetic resonance electrical properties tomography (MREPT) techniques have successfully recovered the conductivity distribution at a Larmor-frequency.

Methods

This work provides a non-invasive method to decompose the high-frequency conductivity into the extracellular medium conductivity based on a two-compartment model using Mb-DWI. To separate the intra- and extracellular micro-structures from the recovered high-frequency conductivity, we include higher b-values DWI and apply the random decision forests to stably determine the micro-structural diffusion parameters.

Results

To demonstrate the proposed method, we conducted phantom and human experiments by comparing the results of reconstructed conductivity of extracellular medium and the conductivity in the intra-neurite and intra-cell body. The phantom and human experiments verify that the proposed method can recover the extracellular electrical properties from the high-frequency conductivity using a routine protocol sequence of MRI scan.

Conclusion

We have proposed a method to decompose the electrical properties in the extracellular, intra-neurite, and soma compartments from the high-frequency conductivity map, reconstructed by solving the electro-magnetic equation with measured B1 phase signals.

Brain Magnetic Resonance Imaging Phenome-Wide Association Study With Metal Transporter Gene SLC39A8

The SLC39A8 gene encodes a divalent metal transporter, ZIP8. SLC39A8 is associated with pleiotropic effects across multiple tissues, including the brain. We determine the different brain magnetic resonance imaging (MRI) phenotypes associated with SLC39A8. We used a phenome-wide association study approach followed by joint and conditional association analysis. Using the summary statistics datasets from a brain MRI genome-wide association study on adult United Kingdom (UK) Biobank participants, we systematically selected all brain MRI phenotypes associated with single-nucleotide polymorphisms (SNPs) within 500 kb of the SLC39A8 genetic locus. For all significant brain MRI phenotypes, we used GCTA-COJO to determine the number of independent association signals and identify index SNPs for each brain MRI phenotype. Linkage equilibrium for brain phenotypes with multiple independent signals was confirmed by LDpair. We identified 24 brain MRI phenotypes that vary due to MRI type and brain region and contain a SNP associated with the SLC39A8 locus. Missense ZIP8 polymorphism rs13107325 was associated with 22 brain MRI phenotypes. Rare ZIP8 variants present in a published UK Biobank dataset are associated with 6 brain MRI phenotypes also linked to rs13107325. Among the 24 datasets, an additional 4 association signals were identified by GCTA-COJO and confirmed to be in linkage equilibrium with rs13107325 using LDpair. These additional association signals represent new probable causative SNPs in addition to rs13107325. This study provides leads into how genetic variation in SLC39A8, a trace mineral transport gene, is linked to brain structure differences and may affect brain development and nervous system function.

Multiple b-values improve discrimination of cortical gray matter regions using diffusion MRI: an experimental validation with a data-driven approach

Objective

To investigate whether varied or repeated b-values provide better diffusion MRI data for discriminating cortical areas with a data-driven approach.

Methods

Data were acquired from three volunteers at 1.5T with b-values of 800, 1400, 2000 s/mm² along 64 diffusion-encoding directions. The diffusion signal was sampled from gray matter in seven regions of interest (ROIs). Rotational invariants of the local diffusion profile were extracted as features that characterize local tissue properties. Random forest classification experiments assessed whether classification accuracy improved when data with multiple b-values were used over repeated acquisition of the same (1400 s/mm²) b-value to compare all possible pairs of the seven ROIs. Three data sets from the Human Connectome Project were subjected to similar processing and analysis pipelines in eight ROIs.

Results

Three different b-values showed an average improvement in correct classification rates of 5.6% and 4.6%, respectively, in the local and HCP data over repeated measurements of the same b-value. The improvement in correct classification rate reached as high as 16% for individual binary classification experiments between two ROIs. Often using only two of the available three b-values were adequate to make such an improvement in classification rates.

Conclusion

Acquisitions with varying b-values are more suitable for discriminating cortical areas.

Homologous laminar organization of the mouse and human subiculum

The subiculum is the major output component of the hippocampal formation and one of the major brain structures most affected by Alzheimer's disease. Our previous work revealed a hidden laminar architecture within the mouse subiculum. However, the rotation of the hippocampal longitudinal axis across species makes it unclear how the laminar organization is represented in human subiculum. Using in situ hybridization data from the Allen Human Brain Atlas, we demonstrate that the human subiculum also contains complementary laminar gene expression patterns similar to the mouse. In addition, we provide evidence that the molecular domain boundaries in human subiculum correspond to microstructural differences observed in high resolution MRI and fiber density imaging. Finally, we show both similarities and differences in the gene expression profile of subiculum pyramidal cells within homologous lamina. Overall, we present a new 3D model of the anatomical organization of human subiculum and its evolution from the mouse.

Elucidating the complex organization of neural micro-domains in the locust Schistocerca gregaria using dMRI

To understand brain function it is necessary to characterize both the underlying structural connectivity between neurons and the physiological integrity of these connections. Previous research exploring insect brain connectivity has typically used electron microscopy techniques, but this methodology cannot be applied to living animals and so cannot be used to understand dynamic physiological processes. The relatively large brain of the desert locust, Schistercera gregaria (Forksȧl) is ideal for exploring a novel methodology; micro diffusion magnetic resonance imaging (micro-dMRI) for the characterization of neuronal connectivity in an insect brain. The diffusion-weighted imaging (DWI) data were acquired on a preclinical system using a customised multi-shell diffusion MRI scheme optimized to image the locust brain. Endogenous imaging contrasts from the averaged DWIs and Diffusion Kurtosis Imaging (DKI) scheme were applied to classify various anatomical features and diffusion patterns in neuropils, respectively. The application of micro-dMRI modelling to the locust brain provides a novel means of identifying anatomical regions and inferring connectivity of large tracts in an insect brain. Furthermore, quantitative imaging indices derived from the kurtosis model that include fractional anisotropy (FA), mean diffusivity (MD) and kurtosis anisotropy (KA) can be extracted. These metrics could, in future, be used to quantify longitudinal structural changes in the nervous system of the locust brain that occur due to environmental stressors or ageing.

Hierarchical Microstructure Informed Tractography

Background: Tractography uses diffusion magnetic resonance imaging to noninvasively infer the macroscopic pathways of white matter fibers and it is the only available technique to probe in vivo the structural connectivity of the brain. However, despite this unique and compelling ability and its wide range of possible neurological applications, tractography is still limited, lacks anatomical precision, and suffers from a serious sensitivity/specificity trade-off. For this reason, in the past few years, tractography postprocessing techniques have emerged and proved effective for improving the quality of the reconstructions. Among them, the Convex Optimization Modeling for Microstructure Informed Tractography formulation allows incorporating the anatomical prior that fibers are naturally organized in fascicles, and has obtained exceptional results in increasing the accuracy of the estimated tractograms. Methods: We propose an extension to this idea and introduce a multilevel grouping of the streamlines to capture the white matter arrangement in fascicles and subfascicles. We tested our proposed formulation in synthetic and in vivo data. Results: Our experiments show that using multiple levels allows considering information about the white matter organization more adequately and helps to improve further the accuracy of the resulting tractograms. Conclusion: This new formulation represents a further important step toward a more accurate structural connectivity estimation.

Secondary Degeneration of White Matter After Focal Sensorimotor Cortical Ischemic Stroke in Rats

Ischemic lesions could lead to secondary degeneration in remote regions of the brain. However, the spatial distribution of secondary degeneration along with its role in functional deficits is not well understood. In this study, we explored the spatial and connectivity properties of white matter (WM) secondary degeneration in a focal unilateral sensorimotor cortical ischemia rat model, using advanced microstructure imaging on a 14 T MRI system. Significant axonal degeneration was observed in the ipsilateral external capsule and even remote regions including the contralesional external capsule and corpus callosum. Further fiber tractography analysis revealed that only fibers having direct axonal connections with the primary lesion exhibited a significant degeneration. These results suggest that focal ischemic lesions may induce remote WM degeneration, but limited to fibers tied to the primary lesion. These “direct” fibers mainly represent perilesional, interhemispheric, and subcortical axonal connections. At last, we found that primary lesion volume might be the determining factor of motor function deficits.

Longitudinal Reproducibility of Neurite Orientation Dispersion and Density Imaging (NODDI) Derived Metrics in the White Matter

Diffusion-weighted magnetic resonance imaging (DWI) is undergoing constant evolution with the ambitious goal of developing in-vivo histology of the brain. A recent methodological advancement is Neurite Orientation Dispersion and Density Imaging (NODDI), a histologically validated multi-compartment model to yield microstructural features of brain tissue such as geometric complexity and neurite packing density, which are especially useful in imaging the white matter. Since NODDI is increasingly popular in clinical research and fields such as developmental neuroscience and neuroplasticity, it is of vast importance to characterize its reproducibility (or reliability). We acquired multi-shell DWI data in 29 healthy young subjects twice over a rescan interval of 4 weeks to assess the within-subject coefficient of variation (CV_WS), between-subject coefficient of variation (CV_BS) and the intraclass correlation coefficient (ICC), respectively. Using these metrics, we compared regional and voxel-by-voxel reproducibility of the most common image analysis approaches (tract-based spatial statistics [TBSS], voxel-based analysis with different extents of smoothing ["VBM-style"], ROI-based analysis). We observed high test-retest reproducibility for the orientation dispersion index (ODI) and slightly worse results for the neurite density index (NDI). Our findings also suggest that the choice of analysis approach might have significant consequences for the results of a study. Collectively, the voxel-based approach with Gaussian smoothing kernels of ≥4 mm FWHM and ROI-averaging yielded the highest reproducibility across NDI and ODI maps (CV_WS mostly ≤3%, ICC mostly ≥0.8), respectively, whilst smaller kernels and TBSS performed consistently worse. Furthermore, we demonstrate that image quality (signal-to-noise ratio [SNR]) is an important determinant of NODDI metric reproducibility. We discuss the implications of these results for longitudinal and cross-sectional research designs commonly employed in the neuroimaging field.

Towards in vivo g-ratio mapping using MRI: Unifying myelin and diffusion imaging

BACKGROUND

The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease.

NEW METHOD

This is the second review on the topic of g-ratio mapping using MRI.

RESULTS

This review summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches.

COMPARISON WITH EXISTING METHODS

Using simulations based on recently published data, this review reveals caveats to the state-of-the-art calibration methods that have been used for in vivo g-ratio mapping. It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology.

CONCLUSIONS

We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the full potential of many novel techniques yet to be investigated.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

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.

Along-tract statistics of neurite orientation dispersion and density imaging diffusion metrics to enhance MR tractography quantitative analysis in healthy controls and in patients with brain tumors

Along‐tract statistics analysis enables the extraction of quantitative diffusion metrics along specific white matter fiber tracts. Besides quantitative metrics derived from classical diffusion tensor imaging (DTI), such as fractional anisotropy and diffusivities, new parameters reflecting the relative contribution of different diffusion compartments in the tissue can be estimated through advanced diffusion MRI methods as neurite orientation dispersion and density imaging (NODDI), leading to a more specific microstructural characterization. In this study, we extracted both DTI‐ and NODDI‐derived quantitative microstructural diffusion metrics along the most eloquent fiber tracts in 15 healthy subjects and in 22 patients with brain tumors. We obtained a robust intraprotocol reference database of normative along‐tract microstructural metrics, and their corresponding plots, from healthy fiber tracts. Each diffusion metric of individual patient's fiber tract was then plotted and statistically compared to the normative profile of the corresponding metric from the healthy fiber tracts. NODDI‐derived metrics appeared to account for the pathological microstructural changes of the peritumoral tissue more accurately than DTI‐derived ones. This approach may be useful for future studies that may compare healthy subjects to patients diagnosed with other pathological conditions.

Realistic Microstructure Simulator (RMS): Monte Carlo simulations of diffusion in three-dimensional cell segmentations of microscopy images

BACKGROUND

Monte Carlo simulations of diffusion are commonly used as a model validation tool as they are especially suitable for generating the diffusion MRI signal in complicated tissue microgeometries.

NEW METHOD

Here we describe the details of implementing Monte Carlo simulations in three-dimensional (3d) voxelized segmentations of cells in microscopy images. Using the concept of the corner reflector, we largely reduce the computational load of simulating diffusion within and exchange between multiple cells. Precision is further achieved by GPU-based parallel computations.

RESULTS

Our simulation of diffusion in white matter axons segmented from a mouse brain demonstrates its value in validating biophysical models. Furthermore, we provide the theoretical background for implementing a discretized diffusion process, and consider the finite-step effects of the particle-membrane reflection and permeation events, needed for efficient simulation of interactions with irregular boundaries, spatially variable diffusion coefficient, and exchange.

COMPARISON WITH EXISTING METHODS

To our knowledge, this is the first Monte Carlo pipeline for MR signal simulations in a substrate composed of numerous realistic cells, accounting for their permeable and irregularly-shaped membranes.

CONCLUSIONS

The proposed RMS pipeline makes it possible to achieve fast and accurate simulations of diffusion in realistic tissue microgeometry, as well as the interplay with other MR contrasts. Presently, RMS focuses on simulations of diffusion, exchange, and T₁ and T₂ NMR relaxation in static tissues, with a possibility to straightforwardly account for susceptibility-induced T₂^* effects and flow.

CRAFT (Cerclage after full dilatation caesarean section): protocol of a mixed methods study investigating the role of previous in-labour caesarean section in preterm birth risk

BACKGROUND

Full dilatation caesarean sections are associated with recurrent early spontaneous preterm birth and late miscarriage. The risk following first stage caesarean sections, are less well defined, but appears to be increased in late-first stage of labour. The mechanism for this increased risk of late miscarriage and early spontaneous preterm birth in these women is unknown and there are uncertainties with regards to clinical management. Current predictive models of preterm birth (based on transvaginal ultrasound and quantitative fetal fibronectin) have not been validated in these women and it is unknown whether the threshold to define a short cervix (≤25 mm) is reliable in predicting the risk of preterm birth. In addition the efficacy of standard treatments or whether benefit may be derived from prophylactic interventions such as a cervical cerclage is unknown.

METHODS

There are three distinct components to the CRAFT project (CRAFT-OBS, CRAFT-RCT and CRAFT-IMG).

CRAFT-OBS

Observational Study; To evaluate subsequent pregnancy risk of preterm birth in women with a prior caesarean section in established labour. This prospective study of cervical length and quantitative fetal fibronectin data will establish a predictive model of preterm birth.

CRAFT-RCT

Randomised controlled trial arm; To assess treatment for short cervix in women at high risk of preterm birth following a fully dilated caesarean section.

CRAFT-IMG

Imaging sub-study; To evaluate the use of MRI and transvaginal ultrasound imaging of micro and macrostructural cervical features which may predispose to preterm birth in women with a previous fully dilated caesarean section, such as scar position and niche.

DISCUSSION

The CRAFT project will quantify the risk of preterm birth or late miscarriage in women with previous in-labour caesarean section, define the best management and shed light on pathological mechanisms so as to improve the care we offer to women and their babies.

TRIAL REGISTRATION

CRAFT was prospectively registered on 25th November 2019 with the ISRCTN registry ( https://doi.org/10.1186/ISRCTN15068651 ).

Diffusion MRI simulation of realistic neurons with SpinDoctor and the Neuron Module

The diffusion MRI signal arising from neurons can be numerically simulated by solving the Bloch-Torrey partial differential equation. In this paper we present the Neuron Module that we implemented within the Matlab-based diffusion MRI simulation toolbox SpinDoctor. SpinDoctor uses finite element discretization and adaptive time integration to solve the Bloch-Torrey partial differential equation for general diffusion-encoding sequences, at multiple b-values and in multiple diffusion directions. In order to facilitate the diffusion MRI simulation of realistic neurons by the research community, we constructed finite element meshes for a group of 36 pyramidal neurons and a group of 29 spindle neurons whose morphological descriptions were found in the publicly available neuron repository NeuroMorpho.Org. These finite elements meshes range from having 15,163 nodes to 622,553 nodes. We also broke the neurons into the soma and dendrite branches and created finite elements meshes for these cell components. Through the Neuron Module, these neuron and cell components finite element meshes can be seamlessly coupled with the functionalities of SpinDoctor to provide the diffusion MRI signal attributable to spins inside neurons. We make these meshes and the source code of the Neuron Module available to the public as an open-source package. To illustrate some potential uses of the Neuron Module, we show numerical examples of the simulated diffusion MRI signals in multiple diffusion directions from whole neurons as well as from the soma and dendrite branches, and include a comparison of the high b-value behavior between dendrite branches and whole neurons. In addition, we demonstrate that the neuron meshes can be used to perform Monte-Carlo diffusion MRI simulations as well. We show that at equivalent accuracy, if only one gradient direction needs to be simulated, SpinDoctor is faster than a GPU implementation of Monte-Carlo, but if many gradient directions need to be simulated, there is a break-even point when the GPU implementation of Monte-Carlo becomes faster than SpinDoctor. Furthermore, we numerically compute the eigenfunctions and the eigenvalues of the Bloch-Torrey and the Laplace operators on the neuron geometries using a finite elements discretization, in order to give guidance in the choice of the space and time discretization parameters for both finite elements and Monte-Carlo approaches. Finally, we perform a statistical study on the set of 65 neurons to test some candidate biomakers that can potentially indicate the soma size. This preliminary study exemplifies the possible research that can be conducted using the Neuron Module.

NOise reduction with DIstribution Corrected (NORDIC) PCA in dMRI with complex-valued parameter-free locally low-rank processing

Diffusion-weighted magnetic resonance imaging (dMRI) has found great utility for a wide range of neuroscientific and clinical applications. However, high-resolution dMRI, which is required for improved delineation of fine brain structures and connectomics, is hampered by its low signal-to-noise ratio (SNR). Since dMRI relies on the acquisition of multiple different diffusion weighted images of the same anatomy, it is well-suited for denoising methods that utilize correlations across the image series to improve the apparent SNR and the subsequent data analysis. In this work, we introduce and quantitatively evaluate a comprehensive framework, NOise Reduction with DIstribution Corrected (NORDIC) PCA method for processing dMRI. NORDIC uses low-rank modeling of g-factor-corrected complex dMRI reconstruction and non-asymptotic random matrix distributions to remove signal components which cannot be distinguished from thermal noise. The utility of the proposed framework for denoising dMRI is demonstrated on both simulations and experimental data obtained at 3 Tesla with different resolutions using human connectome project style acquisitions. The proposed framework leads to substantially enhanced quantitative performance for estimating diffusion tractography related measures and for resolving crossing fibers as compared to a conventional/state-of-the-art dMRI denoising method.

The connections of the insular VEN area in great apes: A histologically-guided ex vivo diffusion tractography study

We mapped the connections of the insular von Economo neuron (VEN) area in ex vivo brains of a bonobo, an orangutan and two gorillas with high angular resolution diffusion MRI imaging acquired in 36 h imaging sessions for each brain. The apes died of natural causes without neurological disorders. The localization of the insular VEN area was based on cresyl violet-stained histological sections from each brain that were coregistered with structural and diffusion images from the same individuals. Diffusion MRI tractography showed that the insular VEN area is connected with olfactory, gustatory, visual and other sensory systems, as well as systems for the mediation of appetite, reward, aversion and motivation. The insular VEN area in apes is most strongly connected with frontopolar cortex, which could support their capacity to choose voluntarily among alternative courses of action particularly in exploring for food resources. The frontopolar cortex may also support their capacity to take note of potential resources for harvesting in the future (prospective memory). All of these faculties may support insight and volitional choice when contemplating courses of action as opposed to rule-based decision-making.

Microglial Density Alters Measures of Axonal Integrity and Structural Connectivity

Diffusion tensor imaging (DTI) has fundamentally transformed how we interrogate diseases and disorders of the brain in neuropsychiatric illness. DTI and recently developed multicompartment diffusion-weighted imaging (MC-DWI) techniques, such as NODDI (neurite orientation dispersion and density imaging), measure diffusion anisotropy presuming a static neuroglial environment; however, microglial morphology and density are highly dynamic in psychiatric illness, and how alterations in microglial density might influence intracellular measures of diffusion anisotropy in DTI and MC-DWI brain microstructure is unknown. To address this question, DTI and MC-DWI studies of murine brains depleted of microglia were performed, revealing significant alterations in axonal integrity and fiber tractography in DTI and in commonly used MC-DWI models. With accumulating evidence of the role of microglia in neuropsychiatric illness, our findings uncover the unexpected contribution of microglia to measures of axonal integrity and structural connectivity and provide unanticipated insights into the potential influence of microglia in diffusion imaging studies of neuropsychiatric disease.

Double diffusion encoding and applications for biomedical imaging

Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.

Magnetic resonance diffusion tensor tractography of a midbrain auditory circuit in Alligator

Knowledge of brain circuitry is critical for understanding the organization, function, and evolution of central nervous systems. Most commonly, brain connections have been elucidated using histological and experimental methods that require animal sacrifice. On the other hand, magnetic resonance diffusion tensor imaging and associated tractography have emerged as a preferred method to noninvasively visualize brain white matter tracts. However, existing studies have primarily examined large, heavily myelinated fiber tracts. Whether tractography can visualize fiber bundles that contain thin and poorly myelinated axons is uncertain. To address this question, the midbrain auditory pathway to the thalamus was investigated in Alligator. This species was chosen because of its evolutionary importance as it is the reptilian group most closely related to birds and because its brain contains many thin and poorly myelinated tracts. Furthermore, this auditory pathway is well documented in other reptiles, including a related crocodilian. Histological observations and experimental determination of anterograde connections confirmed this path in Alligator. Tractography identified these tracts in Alligator and provided a 3-dimensional picture that accurately identified the neural elements of this circuit. In addition, tractography identified one possible unrecognized pathway. These results demonstrate that tractography can visualize circuits containing thin, poorly myelinated fibers. These findings open the door for future studies to examine these types of pathways in other vertebrates.

In vivo microstructural heterogeneity of white matter lesions in healthy elderly and Alzheimer's disease participants using tissue compositional analysis of diffusion MRI data

White matter hyperintensities (WMH) are regions of high signal intensity typically identified on fluid attenuated inversion recovery (FLAIR). Although commonly observed in elderly individuals, they are more prevalent in Alzheimer's disease (AD) patients. Given that WMH appear relatively homogeneous on FLAIR, they are commonly partitioned into location- or distance-based classes when investigating their relevance to disease. Since pathology indicates that such lesions are often heterogeneous, probing their microstructure in vivo may provide greater insight than relying on such arbitrary classification schemes. In this study, we investigated WMH in vivo using an advanced diffusion MRI method known as single-shell 3-tissue constrained spherical deconvolution (SS3T-CSD), which models white matter microstructure while accounting for grey matter and CSF compartments. Diffusion MRI data and FLAIR images were obtained from AD (n = 48) and healthy elderly control (n = 94) subjects. WMH were automatically segmented, and classified: (1) as either periventricular or deep; or (2) into three distance-based contours from the ventricles. The 3-tissue profile of WMH enabled their characterisation in terms of white matter-, grey matter-, and fluid-like characteristics of the diffusion signal. Our SS3T-CSD findings revealed substantial heterogeneity in the 3-tissue profile of WMH, both within lesions and across the various classes. Moreover, this heterogeneity information indicated that the use of different commonly used WMH classification schemes can result in different disease-based conclusions. We conclude that future studies of WMH in AD would benefit from inclusion of microstructural information when characterising lesions, which we demonstrate can be performed in vivo using SS3T-CSD.

Super-Resolved q-Space deep learning with uncertainty quantification

Diffusion magnetic resonance imaging (dMRI) provides a noninvasive method for measuring brain tissue microstructure. q-Space deep learning(q-DL) methods have been developed to accurately estimate tissue microstructure from dMRI scans acquired with a reduced number of diffusion gradients. In these methods, deep networks are trained to learn the mapping directly from diffusion signals to tissue microstructure. However, the quality of tissue microstructure estimation can be limited not only by the reduced number of diffusion gradients but also by the low spatial resolution of typical dMRI acquisitions. Therefore, in this work we extend q-DL to super-resolved tissue microstructure estimation and propose super-resolvedq-DL (SR-q-DL), where deep networks are designed to map low-resolution diffusion signals undersampled in the q-space to high-resolution tissue microstructure. Specifically, we use a patch-based strategy, where a deep network takes low-resolution patches of diffusion signals as input and outputs high-resolution tissue microstructure patches. The high-resolution patches are then combined to obtain the final high-resolution tissue microstructure map. Motivated by existing q-DL methods, we integrate the sparsity of diffusion signals in the network design, which comprises two functional components. The first component computes sparse representation of diffusion signals for the low-resolution input patch, and the second component maps the low-resolution sparse representation to high-resolution tissue microstructure. The weights in the two components are learned jointly and the trained network performs end-to-end tissue microstructure estimation. In addition to SR-q-DL, we further propose probabilistic SR-q-DL, which can quantify the uncertainty of the network output as well as achieve improved estimation accuracy. In probabilistic SR-q-DL, a deep ensemble strategy is used. Specifically, the deep network for SR-q-DL is revised to produce not only tissue microstructure estimates but also the uncertainty of the estimates. Then, multiple deep networks are trained and their results are fused for the final prediction of high-resolution tissue microstructure and uncertainty quantification. The proposed method was evaluated on two independent datasets of brain dMRI scans. Results indicate that our approach outperforms competing methods in terms of estimation accuracy. In addition, uncertainty measures provided by our method correlate with estimation errors, which indicates potential application of the proposed uncertainty quantification method in brain studies.

The present and the future of microstructure MRI: From a paradigm shift to normal science

The aspiration of imaging tissue microstructure with MRI is to uncover micrometer-scale tissue features within millimeter-scale imaging voxels, in vivo. This kind of super-resolution has fueled a paradigm shift within the biomedical imaging community. However, what feels like an ongoing revolution in MRI, has been conceptually experienced in physics decades ago; from this point of view, our current developments can be seen as Thomas Kuhn's "normal science" stage of progress. While the concept of model-based quantification below the nominal imaging resolution is not new, its possibilities in neuroscience and neuroradiology are only beginning to be widely appreciated. This disconnect calls for communicating the progress of tissue microstructure MR imaging to its potential users. Here, a number of recent research developments are outlined in terms of the overarching concept of coarse-graining the tissue structure over an increasing diffusion length. A variety of diffusion models and phenomena are summarized on the phase diagram of diffusion MRI, with the unresolved problems and future directions corresponding to its unexplored domains.

Low-frequency dominant electrical conductivity imaging of in vivo human brain using high-frequency conductivity at Larmor-frequency and spherical mean diffusivity without external injection current

Diffusion weighted imaging based on random Brownian motion of water molecules within a voxel provides information on the micro-structure of biological tissues through water molecule diffusivity. As the electrical conductivity is primarily determined by the concentration and mobility of ionic charge carriers, the macroscopic electrical conductivity of biological tissues is also related to the diffusion of electrical ions. This paper aims to investigate the low-frequency electrical conductivity by relying on a pre-defined biological model that separates the brain into the intracellular (restricted) and extracellular (hindered) compartments. The proposed method uses B1 mapping technique, which provides a high-frequency conductivity distribution at Larmor frequency, and the spherical mean technique, which directly estimates the microscopic tissue structure based on the water molecule diffusivity and neurite orientation distribution. The total extracellular ion concentration, which is separated from the high-frequency conductivity, is recovered using the estimated diffusivity parameters and volume fraction in each compartment. We propose a method to reconstruct the low-frequency dominant conductivity tensor by taking into consideration the extracted extracellular diffusion tensor map and the reconstructed electrical parameters. To demonstrate the reliability of the proposed method, we conducted two phantom experiments. The first one used a cylindrical acrylic cage filled with an agar in the background region and four anomalies for the effect of ion concentration on the electrical conductivity. The other experiment, in which the effect of ion mobility on the conductivity was verified, used cell-like materials with thin insulating membranes suspended in an electrolyte. Animal and human brain experiments were conducted to visualize the low-frequency dominant conductivity tensor images. The proposed method using a conventional MRI scanner can predict the internal current density map in the brain without directly injected external currents.

Neuroanatomical underpinning of diffusion kurtosis measurements in the cerebral cortex of healthy macaque brains

PURPOSE

To investigate the neuroanatomical underpinning of healthy macaque brain cortical microstructure measured by diffusion kurtosis imaging (DKI), which characterizes non-Gaussian water diffusion.

METHODS

High-resolution DKI was acquired from 6 postmortem macaque brains. Neurofilament density (ND) was quantified based on structure tensor from neurofilament histological images of a different macaque brain sample. After alignment of DKI-derived mean kurtosis (MK) maps to the histological images, MK and histology-based ND were measured at corresponding regions of interests characterized by distinguished cortical MK values in the prefrontal/precentral-postcentral and temporal cortices. Pearson correlation was performed to test significant correlation between these cortical MK and ND measurements.

RESULTS

Heterogeneity of cortical MK across different cortical regions was revealed, with significantly and consistently higher MK measurements in the prefrontal/precentral-postcentral cortex compared to those in the temporal cortex across all six scanned macaque brains. Corresponding higher ND measurements in the prefrontal/precentral-postcentral cortex than in the temporal cortex were also found. The heterogeneity of cortical MK is associated with heterogeneity of histology-based ND measurements, with significant correlation between cortical MK and corresponding ND measurements (P < .005).

CONCLUSION

These findings suggested that DKI-derived MK can potentially be an effective noninvasive biomarker quantifying underlying neuroanatomical complexity inside the cerebral cortical mantle for clinical and neuroscientific research.

White matter microstructure across the adult lifespan: A mixed longitudinal and cross-sectional study using advanced diffusion models and brain-age prediction

The macro- and microstructural architecture of human brain white matter undergoes substantial alterations throughout development and ageing. Most of our understanding of the spatial and temporal characteristics of these lifespan adaptations come from magnetic resonance imaging (MRI), including diffusion MRI (dMRI), which enables visualisation and quantification of brain white matter with unprecedented sensitivity and detail. However, with some notable exceptions, previous studies have relied on cross-sectional designs, limited age ranges, and diffusion tensor imaging (DTI) based on conventional single-shell dMRI. In this mixed cross-sectional and longitudinal study (mean interval: 15.2 months) including 702 multi-shell dMRI datasets, we combined complementary dMRI models to investigate age trajectories in healthy individuals aged 18 to 94 years (57.12% women). Using linear mixed effect models and machine learning based brain age prediction, we assessed the age-dependence of diffusion metrics, and compared the age prediction accuracy of six different diffusion models, including diffusion tensor (DTI) and kurtosis imaging (DKI), neurite orientation dispersion and density imaging (NODDI), restriction spectrum imaging (RSI), spherical mean technique multi-compartment (SMT-mc), and white matter tract integrity (WMTI). The results showed that the age slopes for conventional DTI metrics (fractional anisotropy [FA], mean diffusivity [MD], axial diffusivity [AD], radial diffusivity [RD]) were largely consistent with previous research, and that the highest performing advanced dMRI models showed comparable age prediction accuracy to conventional DTI. Linear mixed effects models and Wilk's theorem analysis showed that the 'FA fine' metric of the RSI model and 'orientation dispersion' (OD) metric of the NODDI model showed the highest sensitivity to age. The results indicate that advanced diffusion models (DKI, NODDI, RSI, SMT mc, WMTI) provide sensitive measures of age-related microstructural changes of white matter in the brain that complement and extend the contribution of conventional DTI.

Lattice Boltzmann method for simulation of diffusion magnetic resonance imaging physics in multiphase tissue models

We report an implementation of the lattice Boltzmann method (LBM) to integrate the Bloch-Torrey equation, which describes the evolution of the transverse magnetization vector and the fate of the signal of diffusion magnetic resonance imaging (dMRI). Motivated by the need to interpret dMRI experiments in biological tissues, and to offset the small time-step limitation of classical LBM, a hybrid LBM scheme is introduced and implemented to solve the Bloch-Torrey equation. A membrane boundary condition is presented which is able to accurately represent the effects of thin curvilinear membranes typically found in biological tissues. As implemented, the hybrid LBM scheme accommodates piece-wise uniform transport, dMRI parameters, periodic and mirroring outer boundary conditions, and finite membrane permeabilities on non-boundary-conforming inner boundaries. By comparing with analytical solutions of limiting cases, we demonstrate that the hybrid LBM scheme is more accurate than the classical LBM scheme. The proposed explicit LBM scheme maintains second-order spatial accuracy, stability, and first-order temporal accuracy for a wide range of parameters. The parallel implementation of the hybrid LBM code in a multi-CPU computer system, as well as on GPUs, is straightforward and efficient. Along with offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a flexible scheme that can by easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.

Characterizing Intra-soma Diffusion with Spherical Mean Spectrum Imaging

Most brain microstructure models are dedicated to the quantification of white matter microstructure, using for example sticks, cylinders, and zeppelins to model intra- and extra-axonal environments. Gray matter presents unique micro-architecture with cell bodies (somas) exhibiting diffusion characteristics that differ from axons in white matter. In this paper, we introduce a method to quantify soma microstructure, giving measures such as volume fraction, diffusivity, and kurtosis. Our method captures a spectrum of diffusion patterns and scales and does not rely on restrictive model assumptions. We show that our method yields unique and meaningful contrasts that are in agreement with histological data. We demonstrate its application in the mapping of the distinct spatial patterns of soma density in the cortex.

Challenges for biophysical modeling of microstructure

The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Diffusion-weighted MRI in neurodegenerative and psychiatric animal models: Experimental strategies and main outcomes

Preclinical MRI approaches constitute a key tool to study a wide variety of neurological and psychiatric illnesses, allowing a more direct investigation of the disorder substrate and, at the same time, the possibility of back-translating such findings to human subjects. However, the lack of consensus on the optimal experimental scheme used to acquire the data has led to relatively high heterogeneity in the choice of protocols, which can potentially impact the comparison between results obtained by different groups, even using the same animal model. This is especially true for diffusion-weighted MRI data, where certain experimental choices can impact not only on the accuracy and precision of the extracted biomarkers, but also on their biological meaning. With this in mind, we extensively examined preclinical imaging studies that used diffusion-weighted MRI to investigate neurodegenerative, neurodevelopmental and psychiatric disorders in rodent models. In this review, we discuss the main findings for each preclinical model, with a special focus on the analysis and comparison of the different acquisition strategies used across studies and their impact on the heterogeneity of the findings.

NODDI in clinical research

Diffusion MRI (dMRI) has proven to be a useful imaging approach for both clinical diagnosis and research investigating the microstructures of nervous tissues, and it has helped us to better understand the neurophysiological mechanisms of many diseases. Though diffusion tensor imaging (DTI) has long been the default tool to analyze dMRI data in clinical research, acquisition with stronger diffusion weightings beyond the DTI regimen is now possible with modern clinical scanners, potentially enabling even more detailed characterization of tissue microstructures. To take advantage of such data, neurite orientation dispersion and density imaging (NODDI) has been proposed as a way to relate the dMRI signal to tissue features via biophysically inspired modeling. The number of reports demonstrating the potential clinical utility of NODDI is rapidly increasing. At the same time, the pitfalls and limitations of NODDI, and general challenges in microstructure modeling, are becoming increasingly recognized by clinicians. dMRI microstructure modeling is a rapidly evolving field with great promise, where people from different scientific backgrounds, such as physics, medicine, biology, neuroscience, and statistics, are collaborating to build novel tools that contribute to improving human healthcare. Here, we review the applications of NODDI in clinical research and discuss future perspectives for investigations toward the implementation of dMRI microstructure imaging in clinical practice.

ConFiG: Contextual Fibre Growth to generate realistic axonal packing for diffusion MRI simulation

This paper presents Contextual Fibre Growth (ConFiG), an approach to generate white matter numerical phantoms by mimicking natural fibre genesis. ConFiG grows fibres one-by-one, following simple rules motivated by real axonal guidance mechanisms. These simple rules enable ConFiG to generate phantoms with tuneable microstructural features by growing fibres while attempting to meet morphological targets such as user-specified density and orientation distribution. We compare ConFiG to the state-of-the-art approach based on packing fibres together by generating phantoms in a range of fibre configurations including crossing fibre bundles and orientation dispersion. Results demonstrate that ConFiG produces phantoms with up to 20% higher densities than the state-of-the-art, particularly in complex configurations with crossing fibres. We additionally show that the microstructural morphology of ConFiG phantoms is comparable to real tissue, producing diameter and orientation distributions close to electron microscopy estimates from real tissue as well as capturing complex fibre cross sections. Signals simulated from ConFiG phantoms match real diffusion MRI data well, showing that ConFiG phantoms can be used to generate realistic diffusion MRI data. This demonstrates the feasibility of ConFiG to generate realistic synthetic diffusion MRI data for developing and validating microstructure modelling approaches.

•

We present ConFiG, a biologically motivated numerical phantom generator for white matter.

•

ConFiG produces phantoms with state-of-the-art density and realistic microstructure.

•

Diffusion MRI simulations in ConFiG phantoms are comparable to real dMRI signals.

A new method for accurate in vivo mapping of human brain connections using microstructural and anatomical information

Diffusion magnetic resonance imaging is a noninvasive imaging modality that has been extensively used in the literature to study the neuronal architecture of the brain in a wide range of neurological conditions using tractography. However, recent studies highlighted that the anatomical accuracy of the reconstructions is inherently limited and challenged its appropriateness. Several solutions have been proposed to tackle this issue, but none of them proved effective to overcome this fundamental limitation. In this work, we present a novel processing framework to inject into the reconstruction problem basic prior knowledge about brain anatomy and its organization and evaluate its effectiveness using both simulated and real human brain data. Our results indicate that our proposed method dramatically increases the accuracy of the estimated brain networks and, thus, represents a major step forward for the study of connectivity.

Diffusion Imaging in the Post HCP Era

Diffusion imaging is a critical component in the pursuit of developing a better understanding of the human brain. Recent technical advances promise enabling the advancement in the quality of data that can be obtained. In this review the context for different approaches relative to the Human Connectome Project are compared. Significant new gains are anticipated from the use of high-performance head gradients. These gains can be particularly large when the high-performance gradients are employed together with ultrahigh magnetic fields. Transmit array designs are critical in realizing high accelerations in diffusion-weighted (d)MRI acquisitions, while maintaining large field of view (FOV) coverage, and several techniques for optimal signal-encoding are now available. Reconstruction and processing pipelines that precisely disentangle the acquired neuroanatomical information are established and provide the foundation for the application of deep learning in the advancement of dMRI for complex tissues. Level of Evidence: 3 Technical Efficacy Stage: Stage 3.

Disentangling white-matter damage from physiological fibre orientation dispersion in multiple sclerosis

Multiple sclerosis leads to diffuse damage of the central nervous system, affecting also the normal-appearing white matter. Demyelination and axonal degeneration reduce regional fractional anisotropy in normal-appearing white matter, which can be routinely mapped with diffusion tensor imaging. However, the standard fractional anisotropy metric is also sensitive to physiological variations in orientation dispersion of white matter fibres. This complicates the detection of disease-related damage in large parts of cerebral white matter where microstructure physiologically displays a high degree of fibre dispersion. To resolve this ambiguity, we employed a novel tensor-valued encoding method for diffusion MRI, which yields a microscopic fractional anisotropy metric that is unaffected by regional variations in orientation dispersion. In 26 patients with relapsing-remitting multiple sclerosis, 14 patients with primary-progressive multiple sclerosis and 27 age-matched healthy controls, we compared standard fractional anisotropy mapping with the novel microscopic fractional anisotropy mapping method, focusing on normal-appearing white matter. Mean microscopic fractional anisotropy and standard fractional anisotropy of normal-appearing white matter were significantly reduced in both patient groups relative to healthy controls, but microscopic fractional anisotropy yielded a better reflection of disease-related white-matter alterations. The reduction in mean microscopic fractional anisotropy showed a significant positive linear relationship with physical disability, as reflected by the expanded disability status scale. Mean reduction of microscopic fractional anisotropy in normal-appearing white matter also scaled positively with individual cognitive dysfunction, as measured with the symbol digit modality test. Mean microscopic fractional anisotropy reduction in normal-appearing white matter also showed a positive relationship with total white-matter lesion load as well as lesion load in specific tract systems. None of these relationships between normal-appearing white-matter microstructure and clinical, cognitive or structural measures emerged when using mean fractional anisotropy. Together, the results provide converging evidence that microscopic fractional anisotropy mapping substantially advances the assessment of cerebral white matter in multiple sclerosis by disentangling microstructure damage from variations in physiological fibre orientation dispersion at the stage of data acquisition. Since tensor-valued encoding can be implemented in routine diffusion MRI, microscopic fractional anisotropy mapping bears considerable potential for the future assessment of disease progression in normal-appearing white matter in both relapsing-remitting and progressive forms of multiple sclerosis as well as other white-matter-related brain diseases.