Diffusion tensor distribution imaging of an in vivo mouse brain at ultrahigh magnetic field by spatiotemporal encoding

Resiliency to Alzheimer's disease neuropathology can be distinguished from dementia using cortical astrogliosis imaging

Despite the presence of significant Alzheimer's disease (AD) pathology, characterized by amyloid β (Aβ) plaques and phosphorylated tau (pTau) tangles, some cognitively normal elderly individuals do not inevitably develop dementia. These findings give rise to the notion of cognitive 'resilience', suggesting maintained cognitive function despite the presence of AD neuropathology, highlighting the influence of factors beyond classical pathology. Cortical astroglial inflammation, a ubiquitous feature of symptomatic AD, shows a strong correlation with cognitive impairment severity, potentially contributing to the diversity of clinical presentations. However, noninvasively imaging neuroinflammation, particularly astrogliosis, using MRI remains a significant challenge. Here we sought to address this challenge and to leverage multidimensional (MD) MRI, a powerful approach that combines relaxation with diffusion MR contrasts, to map cortical astrogliosis in the human brain by accessing sub-voxel information. Our goal was to test whether MD-MRI can map astroglial pathology in the cerebral cortex, and if so, whether it can distinguish cognitive resiliency from dementia in the presence of hallmark AD neuropathological changes. We adopted a multimodal approach by integrating histological and MRI analyses using human postmortem brain samples. Ex vivo cerebral cortical tissue specimens derived from three groups comprised of non-demented individuals with significant AD pathology postmortem, individuals with both AD pathology and dementia, and non-demented individuals with minimal AD pathology postmortem as controls, underwent MRI at 7 T. We acquired and processed MD-MRI, diffusion tensor, and quantitative T 1 and T 2 MRI data, followed by histopathological processing on slices from the same tissue. By carefully co-registering MRI and microscopy data, we performed quantitative multimodal analyses, leveraging targeted immunostaining to assess MD-MRI sensitivity and specificity towards Aβ, pTau, and glial fibrillary acidic protein (GFAP), a marker for astrogliosis. Our findings reveal a distinct MD-MRI signature of cortical astrogliosis, enabling the creation of predictive maps for cognitive resilience amid AD neuropathological changes. Multiple linear regression linked histological values to MRI changes, revealing that the MD-MRI cortical astrogliosis biomarker was significantly associated with GFAP burden (standardized β=0.658, pFDR<0.0001), but not with Aβ (standardized β=0.009, p FDR =0.913) or pTau (standardized β=-0.196, p FDR =0.051). Conversely, none of the conventional MRI parameters showed significant associations with GFAP burden in the cortex. While the extent to which pathological glial activation contributes to neuronal damage and cognitive impairment in AD is uncertain, developing a noninvasive imaging method to see its affects holds promise from a mechanistic perspective and as a potential predictor of cognitive outcomes.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Nonparametric 5D D-R2 distribution imaging with single-shot EPI at 21.1 T: Initial results for in vivo rat brain

In vivo human diffusion MRI is by default performed using single-shot EPI with greater than 50-ms echo times and associated signal loss from transverse relaxation. The individual benefits of the current trends of increasing B₀ to boost SNR and employing more advanced signal preparation schemes to improve the specificity for selected microstructural properties eventually may be cancelled by increased relaxation rates at high B₀ and echo times with advanced encoding. Here, initial attempts to translate state-of-the-art diffusion-relaxation correlation methods from 3 T to 21.1 T are made to identify hurdles that need to be overcome to fulfill the promises of both high SNR and readily interpretable microstructural information.

Whole-Brain Imaging of Subvoxel T1-Diffusion Correlation Spectra in Human Subjects

T1 relaxation and water mobility generate eloquent MRI tissue contrasts with great diagnostic value in many neuroradiological applications. However, conventional methods do not adequately quantify the microscopic heterogeneity of these important biophysical properties within a voxel, and therefore have limited biological specificity. We describe a new correlation spectroscopic (CS) MRI method for measuring how T1 and mean diffusivity (MD) co-vary in microscopic tissue environments. We develop a clinical pulse sequence that combines inversion recovery (IR) with single-shot isotropic diffusion encoding (IDE) to efficiently acquire whole-brain MRIs with a wide range of joint T1-MD weightings. Unlike conventional diffusion encoding, the IDE preparation ensures that all subvoxel water pools are weighted by their MDs regardless of the sizes, shapes, and orientations of their corresponding microscopic diffusion tensors. Accordingly, IR-IDE measurements are well-suited for model-free, quantitative spectroscopic analysis of microscopic water pools. Using numerical simulations, phantom experiments, and data from healthy volunteers we demonstrate how IR-IDE MRIs can be processed to reconstruct maps of two-dimensional joint probability density functions, i.e., correlation spectra, of subvoxel T1-MD values. In vivo T1-MD spectra show distinct cerebrospinal fluid and parenchymal tissue components specific to white matter, cortical gray matter, basal ganglia, and myelinated fiber pathways, suggesting the potential for improved biological specificity. The one-dimensional marginal distributions derived from the T1-MD correlation spectra agree well with results from other relaxation spectroscopic and quantitative MRI studies, validating the T1-MD contrast encoding and the spectral reconstruction. Mapping subvoxel T1-diffusion correlations in patient populations may provide a more nuanced, comprehensive, sensitive, and specific neuroradiological assessment of the non-specific changes seen on fluid-attenuated inversion recovery (FLAIR) and diffusion-weighted MRIs (DWIs) in cancer, ischemic stroke, or brain injury.

Toward nonparametric diffusion- T1 characterization of crossing fibers in the human brain

PURPOSE

To estimate T 1 for each distinct fiber population within voxels containing multiple brain tissue types.

METHODS

A diffusion- T 1 correlation experiment was carried out in an in vivo human brain using tensor-valued diffusion encoding and multiple repetition times. The acquired data were inverted using a Monte Carlo algorithm that retrieves nonparametric distributions P ( D , R 1 ) of diffusion tensors and longitudinal relaxation rates R 1 = 1 / T 1 . Orientation distribution functions (ODFs) of the highly anisotropic components of P ( D , R 1 ) were defined to visualize orientation-specific diffusion-relaxation properties. Finally, Monte Carlo density-peak clustering (MC-DPC) was performed to quantify fiber-specific features and investigate microstructural differences between white matter fiber bundles.

RESULTS

Parameter maps corresponding to P ( D , R 1 ) 's statistical descriptors were obtained, exhibiting the expected R 1 contrast between brain tissue types. Our ODFs recovered local orientations consistent with the known anatomy and indicated differences in R 1 between major crossing fiber bundles. These differences, confirmed by MC-DPC, were in qualitative agreement with previous model-based works but seem biased by the limitations of our current experimental setup.

CONCLUSIONS

Our Monte Carlo framework enables the nonparametric estimation of fiber-specific diffusion- T 1 features, thereby showing potential for characterizing developmental or pathological changes in T 1 within a given fiber bundle, and for investigating interbundle T 1 differences.

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.

High-Resolution 3D in vivo Brain Diffusion Tensor Imaging at Ultrahigh Fields: Following Maturation on Juvenile and Adult Mice

Diffusion tensor imaging (DTI) is a well-established technique for mapping brain microstructure and white matter tracts in vivo. High resolution DTI, however, is usually associated with low intrinsic sensitivity and therefore long acquisition times. By increasing sensitivity, high magnetic fields can alleviate these demands, yet high fields are also typically associated with significant susceptibility-induced image distortions. This study explores the potential arising from employing new pulse sequences and emerging hardware at ultrahigh fields, to overcome these limitations. To this end, a 15.2 T MRI instrument equipped with a cryocooled surface transceiver coil was employed, and DTI experiments were compared between SPatiotemporal ENcoding (SPEN), a technique that tolerates well susceptibility-induced image distortions, and double-sampled Spin-Echo Echo-Planar Imaging (SE-EPI) methods. Following optimization, SE-EPI afforded whole brain DTI maps at 135 μm isotropic resolution that possessed higher signal-to-noise ratios (SNRs) than SPEN counterparts. SPEN, however, was a better alternative to SE-EPI when focusing on challenging regions of the mouse brain -including the olfactory bulb and the cerebellum. In these instances, the higher robustness of fully refocused SPEN acquisitions coupled to its built-in zooming abilities, provided in vivo DTI maps with 75 μm nominal isotropic spatial resolution. These DTI maps, and in particular the mean diffusion direction (MDD) details, exhibited variations that matched very well the anatomical features known from histological brain Atlases. Using these capabilities, the development of the olfactory bulb (OB) in live mice was followed from week 1 post-partum, until adulthood. The diffusivity of this organ showed a systematic decrease in its overall isotropic value and increase in its fractional anisotropy with age; this maturation was observed for all regions used in the OB's segmentation but was most evident for the lobules' centers, in particular for the granular cell layer. The complexity of the OB neuronal connections also increased during maturation, as evidenced by the growth in directionalities arising in the mean diffusivity direction maps.