A novel imaging marker of cortical "cellularity" in multiple sclerosis patients

Pathological data showed focal inflammation and regions of diffuse neuronal loss in the cortex of people with multiple sclerosis (MS). In this work, we applied a novel model ("soma and neurite density imaging (SANDI)") to multishell diffusion-weighted MRI data acquired in healthy subjects and people with multiple sclerosis (pwMS), in order to investigate inflammation and degeneration-related changes in the cortical tissue of pwMS. We aimed to (i) establish whether SANDI is applicable in vivo clinical data; (ii) investigate inflammatory and degenerative changes using SANDI soma fraction (fsoma)-a marker of cellularity-in both cortical lesions and in the normal-appearing-cortex and (iii) correlate SANDI fsoma with clinical and biological measures in pwMS. We applied a simplified version of SANDI to a clinical scanners. We then provided evidence that pwMS exhibited an overall decrease in cortical SANDI fsoma compared to healthy subjects, suggesting global degenerative processes compatible with neuronal loss. On the other hand, we have found that progressive pwMS showed a higher SANDI fsoma in the outer part of the cortex compared to relapsing-remitting pwMS, possibly supporting current pathological knowledge of increased innate inflammatory cells in these regions. A similar finding was obtained in subpial lesions in relapsing-remitting patients, reflecting existing pathological data in these lesion types. A significant correlation was found between SANDI fsoma and serum neurofilament light chain-a biomarker of inflammatory axonal damage-suggesting a relationship between SANDI soma fraction and inflammatory processes in pwMS again. Overall, our data show that SANDI fsoma is a promising biomarker to monitor changes in cellularity compatible with neurodegeneration and neuroinflammation in the cortex of MS patients.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

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

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

Microglial morphometric analysis: so many options, so little consistency

Quantification of microglial activation through morphometric analysis has long been a staple of the neuroimmunologist's toolkit. Microglial morphological phenomics can be conducted through either manual classification or constructing a digital skeleton and extracting morphometric data from it. Multiple open-access and paid software packages are available to generate these skeletons via semi-automated and/or fully automated methods with varying degrees of accuracy. Despite advancements in methods to generate morphometrics (quantitative measures of cellular morphology), there has been limited development of tools to analyze the datasets they generate, in particular those containing parameters from tens of thousands of cells analyzed by fully automated pipelines. In this review, we compare and critique the approaches using cluster analysis and machine learning driven predictive algorithms that have been developed to tackle these large datasets, and propose improvements for these methods. In particular, we highlight the need for a commitment to open science from groups developing these classifiers. Furthermore, we call attention to a need for communication between those with a strong software engineering/computer science background and neuroimmunologists to produce effective analytical tools with simplified operability if we are to see their wide-spread adoption by the glia biology community.

In Vivo Myelin Water Quantification Using Diffusion-Relaxation Correlation MRI: A Comparison of 1D and 2D Methods

Multidimensional Magnetic Resonance Imaging (MRI) is a versatile tool for microstructure mapping. We use a diffusion weighted inversion recovery spin echo (DW-IR-SE) sequence with spiral readouts at ultra-strong gradients to acquire a rich diffusion-relaxation data set with sensitivity to myelin water. We reconstruct 1D and 2D spectra with a two-step convex optimization approach and investigate a variety of multidimensional MRI methods, including 1D multi-component relaxometry, 1D multi-component diffusometry, 2D relaxation correlation imaging, and 2D diffusion-relaxation correlation spectroscopic imaging (DR-CSI), in terms of their potential to quantify tissue microstructure, including the myelin water fraction (MWF). We observe a distinct spectral peak that we attribute to myelin water in multi-component T1 relaxometry, T1-T2 correlation, T1-D correlation, and T2-D correlation imaging. Due to lower achievable echo times compared to diffusometry, MWF maps from relaxometry have higher quality. Whilst 1D multi-component T1 data allows much faster myelin mapping, 2D approaches could offer unique insights into tissue microstructure and especially myelin diffusion.

Higher order diffusion imaging as a putative index of human sleep-related microstructural changes and glymphatic clearance

The brain has a unique macroscopic waste clearance system, termed the glymphatic system which utilises perivascular tunnels surrounded by astroglia to promote cerebrospinal-interstitial fluid exchange. Rodent studies have demonstrated a marked increase in glymphatic clearance during sleep which has been linked to a sleep-induced expansion of the extracellular space and concomitant reduction in intracellular volume. However, despite being implicated in the pathophysiology of multiple human neurodegenerative disorders, non-invasive techniques for imaging glymphatic clearance in humans are currently limited. Here we acquired multi-shell diffusion weighted MRI (dwMRI) in twenty-one healthy young participants (6 female, 22.3 ± 3.2 years) each scanned twice, once during wakefulness and once during sleep induced by a combination of one night of sleep deprivation and 10 mg of the hypnotic zolpidem 30 min before scanning. To capture hypothesised sleep-associated changes in intra/extracellular space, dwMRI were analysed using higher order diffusion modelling with the prediction that sleep-associated increases in interstitial (extracellular) fluid volume would result in a decrease in diffusion kurtosis, particularly in areas associated with slow wave generation at the onset of sleep. In line with our hypothesis, we observed a global reduction in diffusion kurtosis (t15=2.82, p = 0.006) during sleep as well as regional reductions in brain areas associated with slow wave generation during early sleep and default mode network areas that are highly metabolically active during wakefulness. Analysis with a higher-order representation of diffusion (MAP-MRI) further indicated that changes within the intra/extracellular domain rather than membrane permeability likely underpin the observed sleep-associated decrease in kurtosis. These findings identify higher-order modelling of dwMRI as a potential new non-invasive method for imaging glymphatic clearance and extend rodent findings to suggest that sleep is also associated with an increase in interstitial fluid volume in humans.

Comprehensive Brain Tumour Characterisation with VERDICT-MRI: Evaluation of Cellular and Vascular Measures Validated by Histology

The aim of this work was to extend the VERDICT-MRI framework for modelling brain tumours, enabling comprehensive characterisation of both intra- and peritumoural areas with a particular focus on cellular and vascular features. Diffusion MRI data were acquired with multiple b-values (ranging from 50 to 3500 s/mm2), diffusion times, and echo times in 21 patients with brain tumours of different types and with a wide range of cellular and vascular features. We fitted a selection of diffusion models that resulted from the combination of different types of intracellular, extracellular, and vascular compartments to the signal. We compared the models using criteria for parsimony while aiming at good characterisation of all of the key histological brain tumour components. Finally, we evaluated the parameters of the best-performing model in the differentiation of tumour histotypes, using ADC (Apparent Diffusion Coefficient) as a clinical standard reference, and compared them to histopathology and relevant perfusion MRI metrics. The best-performing model for VERDICT in brain tumours was a three-compartment model accounting for anisotropically hindered and isotropically restricted diffusion and isotropic pseudo-diffusion. VERDICT metrics were compatible with the histological appearance of low-grade gliomas and metastases and reflected differences found by histopathology between multiple biopsy samples within tumours. The comparison between histotypes showed that both the intracellular and vascular fractions tended to be higher in tumours with high cellularity (glioblastoma and metastasis), and quantitative analysis showed a trend toward higher values of the intracellular fraction (fic) within the tumour core with increasing glioma grade. We also observed a trend towards a higher free water fraction in vasogenic oedemas around metastases compared to infiltrative oedemas around glioblastomas and WHO 3 gliomas as well as the periphery of low-grade gliomas. In conclusion, we developed and evaluated a multi-compartment diffusion MRI model for brain tumours based on the VERDICT framework, which showed agreement between non-invasive microstructural estimates and histology and encouraging trends for the differentiation of tumour types and sub-regions.

Dysfunction of the hypothalamic-pituitary adrenal axis and its influence on aging: the role of the hypothalamus

As part of the hypothalamic-pituitary adrenal (HPA) axis, the hypothalamus exerts pivotal influence on metabolic and endocrine homeostasis. With age, these processes are subject to considerable change, resulting in increased prevalence of physical disability and cardiac disorders. Yet, research on the aging human hypothalamus is lacking. To assess detailed hypothalamic microstructure in middle adulthood, 39 healthy participants (35-65 years) underwent comprehensive structural magnetic resonance imaging. In addition, we studied HPA axis dysfunction proxied by hair cortisol and waist circumference as potential risk factors for hypothalamic alterations. We provide first evidence of regionally different hypothalamic microstructure, with age effects in its anterior-superior subunit, a critical area for HPA axis regulation. Further, we report that waist circumference was related to increased free water and decreased iron content in this region. In age, hair cortisol was additionally associated with free water content, such that older participants with higher cortisol levels were more vulnerable to free water content increase than younger participants. Overall, our results suggest no general age-related decline in hypothalamic microstructure. Instead, older individuals could be more susceptible to risk factors of hypothalamic decline especially in the anterior-superior subregion, including HPA axis dysfunction, indicating the importance of endocrine and stress management in age.

Temporal Diffusion Ratio (TDR) for imaging restricted diffusion: Optimisation and pre-clinical demonstration

Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.

Joint estimation of relaxation and diffusion tissue parameters for prostate cancer with relaxation-VERDICT MRI

This work presents a biophysical model of diffusion and relaxation MRI for prostate called relaxation vascular, extracellular and restricted diffusion for cytometry in tumours (rVERDICT). The model includes compartment-specific relaxation effects providing T1/T2 estimates and microstructural parameters unbiased by relaxation properties of the tissue. 44 men with suspected prostate cancer (PCa) underwent multiparametric MRI (mp-MRI) and VERDICT-MRI followed by targeted biopsy. We estimate joint diffusion and relaxation prostate tissue parameters with rVERDICT using deep neural networks for fast fitting. We tested the feasibility of rVERDICT estimates for Gleason grade discrimination and compared with classic VERDICT and the apparent diffusion coefficient (ADC) from mp-MRI. The rVERDICT intracellular volume fraction f_ic discriminated between Gleason 3 + 3 and 3 + 4 (p = 0.003) and Gleason 3 + 4 and ≥ 4 + 3 (p = 0.040), outperforming classic VERDICT and the ADC from mp-MRI. To evaluate the relaxation estimates we compare against independent multi-TE acquisitions, showing that the rVERDICT T2 values are not significantly different from those estimated with the independent multi-TE acquisition (p > 0.05). Also, rVERDICT parameters exhibited high repeatability when rescanning five patients (R² = 0.79-0.98; CV = 1-7%; ICC = 92-98%). The rVERDICT model allows for accurate, fast and repeatable estimation of diffusion and relaxation properties of PCa sensitive enough to discriminate Gleason grades 3 + 3, 3 + 4 and ≥ 4 + 3.

Microstructural parameters from DW-MRI for tumour characterization and local recurrence prediction in particle therapy of skull-base chordoma

BACKGROUND

Quantitative imaging such as Diffusion-Weighted MRI (DW-MRI) can be exploited to non-invasively derive patient-specific tumor microstructure information for tumor characterization and local recurrence risk prediction in radiotherapy.

PURPOSE

To characterize tumor microstructure according to proliferative capacity and predict local recurrence through microstructural markers derived from pre-treatment conventional DW-MRI, in skull-base chordoma (SBC) patients treated with proton (PT) and carbon ion (CIRT) radiotherapy.

METHODS

Forty-eight patients affected by SBC, who underwent conventional DW-MRI before treatment and were enrolled for CIRT (n = 25) or PT (n = 23), were retrospectively selected. Clinically verified local recurrence information (LR) and histological information (Ki-67, proliferation index) were collected. Apparent diffusion coefficient (ADC) maps were calculated from pre-treatment DW-MRI and, from these, a set of microstructural parameters (cellular radius R, volume fraction vf, diffusion D) were derived by applying a fine-tuning procedure to a framework employing Monte Carlo simulations on synthetic cell substrates. In addition, apparent cellularity (ρ_app ) was estimated from vf and R for an easier clinical interpretation. Histogram-based metrics (mean, median, variance, entropy) from estimated parameters were considered to investigate differences (Mann-Whitney U-test, α = 0.05) in estimated tumor microstructure in SBCs characterized by low or high cell proliferation (Ki-67). Recurrence-free survival analyses were also performed to assess the ability of the microstructural parameters to stratify patients according to the risk of local recurrence (Kaplan-Meier curves, log-rank test α = 0.05).

RESULTS

Refined microstructural markers revealed optimal capabilities in discriminating patients according to cell proliferation, achieving best results with mean values (p-values were 0.0383, 0.0284, 0.0284, 0.0468, and 0.0088 for ADC, R, vf, D, and ρ_app, respectively). Recurrence-free survival analyses showed significant differences between populations at high and low risk of local recurrence as stratified by entropy values of estimated microstructural parameters (p = 0.0110).

CONCLUSION

Patient-specific microstructural information was non-invasively derived providing potentially useful tools for SBC treatment personalization and optimization in particle therapy.

In vivo quantification of brain soma and neurite density abnormalities in multiple sclerosis

BACKGROUND

Soma and neurite density imaging (SANDI) is a new biophysical model that incorporates soma in addition to neurite density, thus possibly providing more specific information about the complex pathological processes of multiple sclerosis (MS).

PURPOSE

To discriminate the pathological abnormalities of MS white matter (WM) lesions, normal-appearing (NA) WM and cortex and to evaluate the associations among SANDI-derived measures, clinical disability, and conventional MRI variables.

METHODS

Twenty healthy controls (HC) and 23 MS underwent a 3 T brain MRI. Using SANDI on diffusion-weighted sequence, the fractions of neurite (fneurite) and soma (fsoma) were assessed in WM lesions, NAWM, and cortex.

RESULTS

Compared to HC WM, MS NAWM showed lower fneurite (false discovery rate [FDR]-p = 0.011). In MS patients, WM lesions showed lower fneurite and fsoma compared to both HC and MS NAWM (FDR-p < 0.001 for all). In the cortex, MS patients had lower fneurite and fsoma compared to HC (FDR-p ≤ 0.009). Compared to both HC and RRMS, PMS patients had lower fneurite in NAWM (vs HC: FDR-p < 0.001; vs RRMS: FDR-p = 0.003) and cortex (vs HC: FDR-p < 0.001; vs RRMS: p = 0.031, not surviving FDR correction), and lower cortical fsoma (vs HC: FDR-p < 0.001; vs RRMS: FDR-p = 0.009). Compared to HC, PMS also showed a higher fsoma in NAWM (FDR-p = 0.015). Fneurite and fsoma in the different brain compartments were correlated with age, phenotype, disease duration, disability, WM lesion volumes, normalized brain, cortical, and WM volumes (r from - 0.761 to 0.821, FDR-p ≤ 0.4).

CONCLUSIONS

SANDI may represent a clinically relevant model to discriminate different neurodegenerative phenomena that gradually accumulate through MS disease course.

Minimal number of sampling directions for robust measures of the spherical mean diffusion weighted signal: Effects of sampling directions, b-value, signal-to-noise ratio, hardware, and fitting strategy

Several recent multi-compartment diffusion MRI investigations and modeling strategies have utilized the orientationally-averaged, or spherical mean, diffusion-weighted signal to study tissue microstructure of the central nervous system. Most experimental designs sample a large number of diffusion weighted directions in order to calculate the spherical mean signal, however, sampling a subset of these directions may increase scanning efficiency and enable either a decrease in scan time or the ability to sample more diffusion weightings. Here, we aim to determine the minimum number of gradient directions needed for a robust measurement of the spherical mean signal. We used computer simulations to characterize the variation of the measured spherical mean signal as a function of the number of gradient directions, while also investigating the effects of diffusion weighting (b-value), signal-to-noise ratio (SNR), available hardware, and spherical mean fitting strategy. We then utilize empirically acquired data in the brain and spinal cord to validate simulations, showing experimental results are in good agreement with simulations. We summarize these results by providing an intuitive lookup table to facilitate the determination of the minimal number of sampling directions needed for robust spherical mean measurements, and give recommendations based on SNR and experimental conditions.

Longitudinal multimodal MRI characterization of a knock-in mouse model of Huntington's disease reveals early gray and white matter alterations

Pathogenesis of the inherited neurodegenerative disorder Huntington's disease (HD) is progressive with a long presymptomatic phase in which subtle changes occur up to 15 years before the onset of symptoms. Thus, there is a need for early, functional biomarker to better understand disease progression and to evaluate treatment efficacy far from onset. Recent studies have shown that white matter may be affected early in mutant HTT gene carriers. A previous study performed on 12 months old Ki140CAG mice showed reduced glutamate level measured by Chemical Exchange Saturation Transfer of glutamate (gluCEST), especially in the corpus callosum. In this study, we scanned longitudinally Ki140CAG mice with structural MRI, diffusion tensor imaging, gluCEST and magnetization transfer imaging, in order to assess white matter integrity over the life of this mouse model characterized by slow progression of symptoms. Our results show early defects of diffusion properties in the anterior part of the corpus callosum at 5 months of age, preceding gluCEST defects in the same region at 8 and 12 months that spread to adjacent regions. At 12 months, frontal and piriform cortices showed reduced gluCEST, as well as the pallidum. MT imaging showed reduced signal in the septum at 12 months. Cortical and striatal atrophy then appear at 18 months. Vulnerability of the striatum and motor cortex, combined with alterations of anterior corpus callosum, seems to point out the potential role of white matter in the brain dysfunction that characterizes HD and the pertinence of gluCEST and DTI as biomarkers in HD.

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange

Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.

Soma and Neurite Density MRI (SANDI) of the in-vivo mouse brain and comparison with the Allen Brain Atlas

Diffusion MRI (dMRI) provides unique insights into the neural tissue milieu by probing interactions between diffusing molecules and tissue microstructure. Most dMRI techniques focus on white matter (WM) tissues, nevertheless, interest in gray matter characterizations is growing. The Soma and Neurite Density MRI (SANDI) methodology harnesses a model incorporating water diffusion in spherical objects (assumed to be associated with cell bodies) and in impermeable "sticks" (assumed to represent neurites), which potentially enables the characterization of cellular and neurite densities. Recognising the importance of rodents in animal models of development, aging, plasticity, and disease, we here employ SANDI for in-vivo preclinical imaging and provide a first validation of the methodology by comparing SANDI metrics with cellular density reflected by the Allen mouse brain atlas. SANDI was implemented on a 9.4T scanner equipped with a cryogenic coil, and in-vivo experiments were carried out on N = 6 mice. Pixelwise, ROI-based, and atlas comparisons were performed, magnitude vs. real-valued analyses were compared, and shorter acquisitions with reduced the number of b-value shells were investigated. Our findings reveal good reproducibility of the SANDI parameters, including the sphere and stick fractions, as well as sphere size (CoV < 7%, 12% and 3%, respectively). Additionally, we find a very good rank correlation between SANDI-driven sphere fraction and Allen mouse brain atlas contrast that represents cellular density. We conclude that SANDI is a viable preclinical MRI technique that can greatly contribute to research on brain tissue microstructure.

Transient Anomalous Diffusion MRI Measurement Discriminates Porous Polymeric Matrices Characterized by Different Sub-Microstructures and Fractal Dimension

Considering the current development of new nanostructured and complex materials and gels, it is critical to develop a sub-micro-scale sensitivity tool to quantify experimentally new parameters describing sub-microstructured porous systems. Diffusion NMR, based on the measurement of endogenous water's diffusion displacement, offers unique information on the structural features of materials and tissues. In this paper, we applied anomalous diffusion NMR protocols to quantify the subdiffusion of water and to measure, in an alternative, non-destructive and non-invasive modality, the fractal dimension d_w of systems characterized by micro and sub-micro geometrical structures. To this end, three highly heterogeneous porous-polymeric matrices were studied. All the three matrices composed of glycidylmethacrylate-divynilbenzene porous monoliths obtained through the High Internal Phase Emulsion technique were characterized by pores of approximately spherical symmetry, with diameters in the range of 2-10 μm. Pores were interconnected by a plurality of window holes present on pore walls, which were characterized by size coverings in the range of 0.5-2 μm. The walls were characterized by a different degree of surface roughness. Moreover, complementary techniques, namely Field Emission Scanning Electron Microscopy (FE-SEM) and dielectric spectroscopy, were used to corroborate the NMR results. The experimental results showed that the anomalous diffusion α parameter that quantifies subdiffusion and d_w = 2/α changed in parallel to the specific surface area S (or the surface roughness) of the porous matrices, showing a submicroscopic sensitivity. The results reported here suggest that the anomalous diffusion NMR method tested may be a valid experimental tool to corroborate theoretical and simulation results developed and performed for describing highly heterogeneous and complex systems. On the other hand, non-invasive and non-destructive anomalous subdiffusion NMR may be a useful tool to study the characteristic features of new highly heterogeneous nanostructured and complex functional materials and gels useful in cultural heritage applications, as well as scaffolds useful in tissue engineering.

Neural networks for parameter estimation in microstructural MRI: Application to a diffusion-relaxation model of white matter

Specific features of white matter microstructure can be investigated by using biophysical models to interpret relaxation-diffusion MRI brain data. Although more intricate models have the potential to reveal more details of the tissue, they also incur time-consuming parameter estimation that may converge to inaccurate solutions due to a prevalence of local minima in a degenerate fitting landscape. Machine-learning fitting algorithms have been proposed to accelerate the parameter estimation and increase the robustness of the attained estimates. So far, learning-based fitting approaches have been restricted to microstructural models with a reduced number of independent model parameters where dense sets of training data are easy to generate. Moreover, the degree to which machine learning can alleviate the degeneracy problem is poorly understood. For conventional least-squares solvers, it has been shown that degeneracy can be avoided by acquisition with optimized relaxation-diffusion-correlation protocols that include tensor-valued diffusion encoding. Whether machine-learning techniques can offset these acquisition requirements remains to be tested. In this work, we employ artificial neural networks to vastly accelerate the parameter estimation for a recently introduced relaxation-diffusion model of white matter microstructure. We also develop strategies for assessing the accuracy and sensitivity of function fitting networks and use those strategies to explore the impact of the acquisition protocol. The developed learning-based fitting pipelines were tested on relaxation-diffusion data acquired with optimal and sub-optimal acquisition protocols. Networks trained with an optimized protocol were observed to provide accurate parameter estimates within short computational times. Comparing neural networks and least-squares solvers, we found the performance of the former to be less affected by sub-optimal protocols; however, model fitting networks were still susceptible to degeneracy issues and their use could not fully replace a careful design of the acquisition protocol.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Analysis of tumour microstructure estimation from conventional diffusion MRI and application to skull-base chordoma. Annu Int Conf IEEE Eng Med Biol Soc 2021;2021:3761-3764. [PMID: 34892054 DOI: 10.1109/embc46164.2021.9630129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Skull-base chordoma (SBC) is a rare tumour whose molecular and radiological characteristics are still being investigated. In neuro-oncology microstructural imaging techniques, like diffusion-weighted MRI (DW-MRI), have been widely investigated, with the apparent diffusion coefficient (ADC) being one of the most used DW-MRI parameters due to its ease of acquisition and computation. ADC is a potential biomarker without a clear link to microstructure. The aim of this work was to derive microstructural information from conventional ADC, showing its potential for the characterisation of skull-base chordomas. Sixteen patients affected by SBC, who underwent conventional DW-MRI were retrospectively selected. From mono-exponential fits of DW-MRI, ADC maps were estimated using different sets of b-values. DW-MRI signals were simulated from synthetic substrates , which mimic the cellular packing of a tumour tissue with well-defined microstructural features. Starting from a published method, an error-driven procedure was evaluated to improve the estimates of microstructural parameters obtained through the simulated signals. A quantitative description of the tumour microstructure was then obtained from the DW-MRI images. This allowed successfully differentiating patients according to histologically-verified cell proliferation information.Clinical Relevance - The impact on cancer management derives from the expected improvement of radiation treatment quality tailored to a patient-specific non-invasive description of tumour microstructure.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Training data distribution significantly impacts the estimation of tissue microstructure with machine learning

PURPOSE

Supervised machine learning (ML) provides a compelling alternative to traditional model fitting for parameter mapping in quantitative MRI. The aim of this work is to demonstrate and quantify the effect of different training data distributions on the accuracy and precision of parameter estimates when supervised ML is used for fitting.

METHODS

We fit a two- and three-compartment biophysical model to diffusion measurements from in-vivo human brain, as well as simulated diffusion data, using both traditional model fitting and supervised ML. For supervised ML, we train several artificial neural networks, as well as random forest regressors, on different distributions of ground truth parameters. We compare the accuracy and precision of parameter estimates obtained from the different estimation approaches using synthetic test data.

RESULTS

When the distribution of parameter combinations in the training set matches those observed in healthy human data sets, we observe high precision, but inaccurate estimates for atypical parameter combinations. In contrast, when training data is sampled uniformly from the entire plausible parameter space, estimates tend to be more accurate for atypical parameter combinations but may have lower precision for typical parameter combinations.

CONCLUSION

This work highlights that estimation of model parameters using supervised ML depends strongly on the training-set distribution. We show that high precision obtained using ML may mask strong bias, and visual assessment of the parameter maps is not sufficient for evaluating the quality of the estimates.

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.

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.

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.

Data-Driven multi-Contrast spectral microstructure imaging with InSpect: INtegrated SPECTral component estimation and mapping

•

Unsupervised learning technique for spectroscopic analysis of quantitative MRI.

•

Shares information across voxels to improve estimation of multi-dimensional or single-dimensional spectra.

•

Spectral maps are dramatically improved compared to existing approaches.

•

Can potentially identify and map tissue environments; in placental diffusion-relaxometry MRI we demonstrate that it identifies components that correspond to distinct tissue types.

We introduce and demonstrate an unsupervised machine learning technique for spectroscopic analysis of quantitative MRI experiments. Our algorithm supports estimation of one-dimensional spectra from single-contrast data, and multidimensional correlation spectra from simultaneous multi-contrast data. These spectrum-based approaches allow model-free investigation of tissue properties, but require regularised inversion of a Laplace transform or Fredholm integral, which is an ill-posed calculation. Here we present a method that addresses this limitation in a data-driven way. The algorithm simultaneously estimates a canonical basis of spectral components and voxelwise maps of their weightings, thereby pooling information across whole images to regularise the ill-posed problem. We show in simulations that our algorithm substantially outperforms current voxelwise spectral approaches. We demonstrate the method on multi-contrast diffusion-relaxometry placental MRI scans, revealing anatomically-relevant sub-structures, and identifying dysfunctional placentas. Our algorithm vastly reduces the data required to reliably estimate spectra, opening up the possibility of quantitative MRI spectroscopy in a wide range of new applications. Our InSpect code is available at github.com/paddyslator/inspect.

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.

Improving the characterization of meningioma microstructure in proton therapy from conventional apparent diffusion coefficient measurements using Monte Carlo simulations of diffusion MRI

PURPOSE

Proton therapy could benefit from noninvasively gaining tumor microstructure information, at both planning and monitoring stages. The anatomical location of brain tumors, such as meningiomas, often hinders the recovery of such information from histopathology, and conventional noninvasive imaging biomarkers, like the apparent diffusion coefficient (ADC) from diffusion-weighted MRI (DW-MRI), are nonspecific. The aim of this study was to retrieve discriminative microstructural markers from conventional ADC for meningiomas treated with proton therapy. These markers were employed for tumor grading and tumor response assessment.

METHODS

DW-MRIs from patients affected by meningioma and enrolled in proton therapy were collected before (n = 35) and 3 months after (n = 25) treatment. For the latter group, the risk of an adverse outcome was inferred by their clinical history. Using Monte Carlo methods, DW-MRI signals were simulated from packings of synthetic cells built with well-defined geometrical and diffusion properties. Patients' ADC was modeled as a weighted sum of selected simulated signals. The weights that best described a patient's ADC were determined through an optimization procedure and used to estimate a set of markers of tumor microstructure: diffusion coefficient (D), volume fraction (vf), and radius (R). Apparent cellularity (ρ_app ) was estimated from vf and R for an easier clinical interpretability. Differences between meningothelial and atypical subtypes, and low- and high-grade meningiomas were assessed with nonparametric statistical tests, whereas sensitivity and specificity with ROC analyses. Similar analyses were performed for patients showing low or high risk of an adverse outcome to preliminary evaluate response to treatment.

RESULTS

Significant (P < 0.05) differences in median ADC, D, vf, R, and ρ_app values were found when comparing meningiomas' subtypes and grades. ROC analyses showed that estimated microstructural parameters reached higher specificity than ADC for subtyping (0.93 for D and vf vs 0.80 for ADC) and grading (0.75 for R vs 0.67 for ADC). High- and low-risk patients showed significant differences in ADC and microstructural parameters. The skewness of ρ_app was the parameter with highest AUC (0.90) and sensitivity (0.75).

CONCLUSIONS

Matching measured with simulated ADC yielded a set of potential imaging markers for meningiomas grading and response monitoring in proton therapy, showing higher specificity than conventional ADC. These markers can provide discriminative information about spatial patterns of tumor microstructure implying important advantages for patient-specific proton therapy workflows.

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.

Machine learning based white matter models with permeability: An experimental study in cuprizone treated in-vivo mouse model of axonal demyelination

The intra-axonal water exchange time (τ_i), a parameter associated with axonal permeability, could be an important biomarker for understanding and treating demyelinating pathologies such as Multiple Sclerosis. Diffusion-Weighted MRI (DW-MRI) is sensitive to changes in permeability; however, the parameter has so far remained elusive due to the lack of general biophysical models that incorporate it. Machine learning based computational models can potentially be used to estimate such parameters. Recently, for the first time, a theoretical framework using a random forest (RF) regressor suggests that this is a promising new approach for permeability estimation. In this study, we adopt such an approach and for the first time experimentally investigate it for demyelinating pathologies through direct comparison with histology. We construct a computational model using Monte Carlo simulations and an RF regressor in order to learn a mapping between features derived from DW-MRI signals and ground truth microstructure parameters. We test our model in simulations, and find strong correlations between the predicted and ground truth parameters (intra-axonal volume fraction f: R² =0.99, τ_i: R² =0.84, intrinsic diffusivity d: R² =0.99). We then apply the model in-vivo, on a controlled cuprizone (CPZ) mouse model of demyelination, comparing the results from two cohorts of mice, CPZ (N=8) and healthy age-matched wild-type (WT, N=8). We find that the RF model estimates sensible microstructure parameters for both groups, matching values found in literature. Furthermore, we perform histology for both groups using electron microscopy (EM), measuring the thickness of the myelin sheath as a surrogate for exchange time. Histology results show that our RF model estimates are very strongly correlated with the EM measurements (ρ = 0.98 for f, ρ = 0.82 for τ_i). Finally, we find a statistically significant decrease in τ_i in all three regions of the corpus callosum (splenium/genu/body) of the CPZ cohort (<τ_i>=310ms/330ms/350ms) compared to the WT group (<τ_i>=370ms/370ms/380ms). This is in line with our expectations that τ_i is lower in regions where the myelin sheath is damaged, as axonal membranes become more permeable. Overall, these results demonstrate, for the first time experimentally and in vivo, that a computational model learned from simulations can reliably estimate microstructure parameters, including the axonal permeability .

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.

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 ​mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies.

SANDI: A compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI

This work introduces a compartment-based model for apparent cell body (namely soma) and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging through DW-MRI presents water diffusion in white (WM) and gray (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b ​≤ ​3,000 ​s/mm² (or 3 ​ms/μm²), it has been also shown to fail in GM at high b values (b≫3,000 ​s/mm² or 3 ​ms/μm²). Here we hypothesise that the unmodelled soma compartment (i.e. cell body of any brain cell type: from neuroglia to neurons) may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma of any brain cell type is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (b_max ​= ​40,000 ​s/mm², or 40 ​ms/μm²) and human (b_max ​= ​10,000 ​s/mm², or 10 ​ms/μm²) brain. We show that SANDI defines new contrasts representing complementary information on the brain cyto- and myelo-architecture. Indeed, we show maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

Revisiting double diffusion encoding MRS in the mouse brain at 11.7T: Which microstructural features are we sensitive to?

Brain metabolites, such as N-acetylaspartate or myo-inositol, are constantly probing their local cellular environment under the effect of diffusion. Diffusion-weighted NMR spectroscopy therefore presents unparalleled potential to yield cell-type specific microstructural information. Double diffusion encoding (DDE) consists in applying two diffusion blocks, where gradient's direction in the second block is varied during the course of the experiment. Unlike single diffusion encoding, DDE measurements at long mixing time display some angular modulation of the signal amplitude which reflects microscopic anisotropy (μA), while requiring relatively low gradient strength. This angular dependence has been formerly used to quantify cell fiber diameter using a model of isotropically oriented infinite cylinders. However, how additional features of the cell microstructure (such as cell body diameter, fiber length and branching) may also influence the DDE signal has been little explored. Here, we used a cryoprobe as well as state-of-the-art post-processing to perform DDE acquisitions with high accuracy and precision in the mouse brain at 11.7 ​T. We then compared our results to simulated DDE datasets obtained in various 3D cell models in order to pinpoint which features of cell morphology may influence the most the angular dependence of the DDE signal. While the infinite cylinder model poorly fits our experimental data, we show that incorporating branched fiber structure in our model allows more realistic interpretation of the DDE signal. Lastly, data acquired in the short mixing time regime suggest that some sensitivity to cell body diameter might be retrieved, although additional experiments would be required to further support this statement.

Combined diffusion-relaxometry MRI to identify dysfunction in the human placenta

Purpose

A combined diffusion‐relaxometry MR acquisition and analysis pipeline for in vivo human placenta, which allows for exploration of coupling between T2* and apparent diffusion coefficient (ADC) measurements in a sub 10‐minute scan time.

Methods

We present a novel acquisition combining a diffusion prepared spin echo with subsequent gradient echoes. The placentas of 17 pregnant women were scanned in vivo, including both healthy controls and participants with various pregnancy complications. We estimate the joint T2*‐ADC spectra using an inverse Laplace transform.

Results

T2*‐ADC spectra demonstrate clear quantitative separation between normal and dysfunctional placentas.

Conclusions

Combined T2*‐diffusivity MRI is promising for assessing fetal and maternal health during pregnancy. The T2*‐ADC spectrum potentially provides additional information on tissue microstructure, compared to measuring these two contrasts separately. The presented method is immediately applicable to the study of other organs.

Age-related microstructural and physiological changes in normal brain measured by MRI γ-metrics derived from anomalous diffusion signal representation

Nowadays, increasing longevity associated with declining cerebral nervous system functions, suggests the need for continued development of new imaging contrast mechanisms to support the differential diagnosis of age-related decline. In our previous papers, we developed a new imaging contrast metrics derived from anomalous diffusion signal representation and obtained from diffusion-weighted (DW) data collected by varying diffusion gradient strengths. Recently, we highlighted that the new metrics, named γ-metrics, depended on the local inhomogeneity due to differences in magnetic susceptibility between tissues and diffusion compartments in young healthy subjects, thus providing information about myelin orientation and iron content within cerebral regions. The major structural modifications occurring in brain aging are myelinated fibers damage in nerve fibers and iron accumulation in gray matter nuclei. Therefore, we investigated the potential of γ-metrics in relation to other conventional diffusion metrics such as DTI, DKI and NODDI in detecting age-related structural changes in white matter (WM) and subcortical gray matter (scGM). DW-images were acquired in 32 healthy subjects, adults and elderly (age range 20-77 years) using 3.0T and 12 b-values up to 5000 s/mm². Association between diffusion metrics and subjects' age was assessed using linear regression. A decline in mean γ (Mγ) in the scGM and a complementary increase in radial γ (γ_⊥) in frontal WM, genu of corpus callosum and anterior corona radiata with advancing age were found. We suggested that the increase in γ⊥ might reflect declined myelin density, and Mγ decrease might mirror iron accumulation. An increase in D_// and a decrease in the orientation dispersion index (ODI) were associated with axonal loss in the pyramidal tracts, while their inverted trends within the thalamus were thought to be linked to reduced architectural complexity of nerve fibers. γ-metrics together with conventional diffusion-metrics can more comprehensively characterize the complex mechanisms underlining age-related changes than conventional diffusion techniques alone.

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

Brain Metabolite Diffusion from Ultra-Short to Ultra-Long Time Scales: What Do We Learn, Where Should We Go?

In vivo diffusion-weighted MR spectroscopy (DW-MRS) allows measuring diffusion properties of brain metabolites. Unlike water, most metabolites are confined within cells. Hence, their diffusion is expected to purely reflect intracellular properties, opening unique possibilities to use metabolites as specific probes to explore cellular organization and structure. However, interpretation and modeling of DW-MRS, and more generally of intracellular diffusion, remains difficult. In this perspective paper, we will focus on the study of the time-dependency of brain metabolite apparent diffusion coefficient (ADC). We will see how measuring ADC over several orders of magnitude of diffusion times, from less than 1 ms to more than 1 s, allows clarifying our understanding of brain metabolite diffusion, by firmly establishing that metabolites are neither massively transported by active mechanisms nor massively confined in subcellular compartments or cell bodies. Metabolites appear to be instead diffusing in long fibers typical of neurons and glial cells such as astrocytes. Furthermore, we will evoke modeling of ADC time-dependency to evaluate the effect of, and possibly quantify, some structural parameters at various spatial scales, departing from a simple model of hollow cylinders and introducing additional complexity, either short-ranged (such as dendritic spines) or long-ranged (such as cellular fibers ramification). Finally, we will discuss the experimental feasibility and expected benefits of extending the range of diffusion times toward even shorter and longer values.

Multimodal-3D imaging based on μMRI and μCT techniques bridges the gap with histology in visualization of the bone regeneration process

Bone repair/regeneration is usually investigated through X-ray computed microtomography (μCT) supported by histology of extracted samples, to analyse biomaterial structure and new bone formation processes. Magnetic resonance imaging (μMRI) shows a richer tissue contrast than μCT, despite at lower resolution, and could be combined with μCT in the perspective of conducting non-destructive 3D investigations of bone. A pipeline designed to combine μMRI and μCT images of bone samples is here described and applied on samples of extracted human jawbone core following bone graft. We optimized the coregistration procedure between μCT and μMRI images to avoid bias due to the different resolutions and contrasts. Furthermore, we used an Adaptive Multivariate Clustering, grouping homologous voxels in the coregistered images, to visualize different tissue types within a fused 3D metastructure. The tissue grouping matched the 2D histology applied only on 1 slice, thus extending the histology labelling in 3D. Specifically, in all samples, we could separate and map 2 types of regenerated bone, calcified tissue, soft tissues, and/or fat and marrow space. Remarkably, μMRI and μCT alone were not able to separate the 2 types of regenerated bone. Finally, we computed volumes of each tissue in the 3D metastructures, which might be exploited by quantitative simulation. The 3D metastructure obtained through our pipeline represents a first step to bridge the gap between the quality of information obtained from 2D optical microscopy and the 3D mapping of the bone tissue heterogeneity and could allow researchers and clinicians to non-destructively characterize and follow-up bone regeneration.

Anti-Inflammatory, Immunomodulatory, and Tissue Repair Activity on Human Keratinocytes by Green Innovative Nanocomposites

The use of raw materials obtained by waste and processed through innovative industrial methodologies has generated an industry of about a trillion dollars in a short time, and in the near future will provide resources and services for the conservation and sustainable use of natural resources in order to ensure a better and fairer welfare for the human race. The production of nano-fiber chitin non-woven tissue is in accordance with the Organization for Economic Co-operation and Development (OECD) and European Union (EU) bio-economic programs: 100% biodegradable, ecological, and therefore useful in decreasing dependence on fossil fuel resources. The aim of our study is the evaluation of different formulations of a non-woven tissue obtained from electrospinning of a mixture of nanochitin fibrils, lignin, and poly (ethylene) oxide (PEO) on the restoration of damaged tissues. Wound repair is a complex process that involves epithelial and immune cells and includes the induction of metalloproteinases, inflammatory mediators, and angiogenic factors. Our in vitro results have shown that all of the realized chitin nanofibrils-bio-lignin non-woven tissues tested as nontoxic for human keratinocytes (HaCat) cells. Furthermore, the bio-composites that included bio-lignin at 0.1% have been able to modulate the expression of pro-inflammatory cytokines (Tumor Necrosis Factor-α, IL-1α, and IL8), lipopolysaccharide (LPS)-induced, and matrix metalloproteinases (MMPs) and human beta-defensin 2 (HBD-2) expression in HaCat cells, suggesting an anti-inflammatory and immunomodulatory role. Taken together, our results suggest that our chitin nanofibrils-bio-lignin non-woven tissue represents a skin-friendly tool that is able to favor a correct and fast wound repair.

Modeling diffusion of intracellular metabolites in the mouse brain up to very high diffusion-weighting: Diffusion in long fibers (almost) accounts for non-monoexponential attenuation

PURPOSE

To investigate how intracellular metabolites diffusion measured in vivo up to very high q/b in the mouse brain can be explained in terms of simple geometries.

METHODS

10 mice were scanned using our new STE-LASER sequence, at 11.7 Tesla (T), up to qmax = 1 μm-1 at diffusion time td = 63.2 ms, corresponding to bmax = 60 ms/µm². We model cell fibers as randomly oriented cylinders, with radius a and intracellular diffusivity Dintracyl, and fit experimental data as a function of q to estimate Dintracyl and a.

RESULTS

Randomly oriented cylinders account well for measured attenuation, giving fiber radii and Dintracyl in the expected ranges (0.5-1.5 µm and 0.30-0.45 µm2/ms, respectively). The only exception is N-acetyl-aspartate (NAA) (extracted a∼0), which we show to be compatible with a small fraction of the NAA pool being confined in highly restricted compartments (with short T2).

CONCLUSION

The non-monoexponential signal attenuation of intracellular metabolites in the mouse brain can be described by diffusion in long and thin cylinders, yielding realistic Dintra and fiber diameters. However, this simple model may require small "corrections" for NAA, in the form of a small fraction of the NAA signal originating from a highly restricted compartment. Magn Reson Med, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

The γ-parameter of anomalous diffusion quantified in human brain by MRI depends on local magnetic susceptibility differences

Motivated by previous results obtained in vitro, we investigated the dependence of the anomalous diffusion (AD) MRI technique on local magnetic susceptibility differences (Δχ) driven by magnetic field inhomogeneity in human brains. The AD-imaging contrast investigated here is quantified by the stretched-exponential parameter γ, extracted from diffusion weighted (DW) data collected by varying diffusion gradient strengths. We performed T₂* and DW experiments in eight healthy subjects at 3.0T. T₂*-weighted images at different TEs=(10,20,35,55)ms and DW-EPI images with fourteen b-values from 0 to 5000s/mm² were acquired. AD-metrics and Diffusion Tensor Imaging (DTI) parameters were compared and correlated to R₂* and to Δχ values taken from literature for the gray (GM) and the white (WM) matter. Pearson's correlation test and Analysis of Variance with Bonferroni post-hoc test were used. Significant strong linear correlations were found between AD γ-metrics and R₂* in both GM and WM of the human brain, but not between DTI-metrics and R₂*. Depending on Δχ driven magnetic field inhomogeneity, the new contrast provided by AD-γ imaging reflects Δχ due to differences in myelin orientation and iron content within selected regions in the WM and GM, respectively. This feature of the AD-γ imaging due to the fact that γ is quantified by using MRI, may be an alternative strategy to investigate, at high magnetic fields, microstructural changes in myelin, and alterations due to iron accumulation. Possible clinical applications might be in the field of neurodegenerative diseases.

Metabolite diffusion up to very high b in the mouse brain in vivo: Revisiting the potential correlation between relaxation and diffusion properties

PURPOSE

To assess the potential correlation between metabolites diffusion and relaxation in the mouse brain, which is of importance for interpreting and modeling metabolite diffusion based on pure geometry, irrespective of relaxation properties (multicompartmental relaxation or surface relaxivity).

METHODS

A new diffusion-weighted magnetic resonance spectroscopy sequence is introduced, dubbed "STE-LASER," which presents several nice properties, in particular the absence of cross-terms with selection gradients and a very clean localization. Metabolite diffusion is then measured in a large voxel in the mouse brain at 11.7 Tesla using a cryoprobe, resulting in excellent signal-to-noise ratio, up to very high b-values under different echo time, mixing time, and diffusion time combinations.

RESULTS

Our results suggest that the correlation between relaxation and diffusion properties is extremely small or even nonexistent for metabolites in the mouse brain.

CONCLUSION

The present work strongly supports the interpretation and modeling of metabolite diffusion primarily based on geometry, irrespective of relaxation properties, at least under current experimental conditions. Magn Reson Med 77:1390-1398, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.