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

In vivo cortical microstructure mapping using high-gradient diffusion MRI accounting for intercompartmental water exchange effects

In recent years, mapping tissue microstructure in the cortex using high gradient diffusion MRI has received growing attention. The Soma And Neurite Density Imaging (SANDI) explicitly models the soma compartment in the cortex assuming impermeable membranes. As such, it does not account for diffusion time dependence due to water exchange in the estimated microstructural properties, as neurites in gray matter are much less myelinated than in white matter. In this work, we performed a systematic evaluation of an extended SANDI model for in vivo human cortical microstructural mapping that accounts for water exchange effects between the neurite and extracellular compartments using the anisotropic Kärger model. We refer to this model as in vivo SANDIX, adapting the nomenclature from previous publications. As in the original SANDI model, the soma compartment is modeled as an impermeable sphere due to the much smaller surface-to-volume ratio compared to the neurite compartment. A Monte Carlo simulation study was performed to examine the sensitivity of the in vivo SANDIX model to sphere radii, compartment fractions, and water exchange times. The simulation results indicate that the proposed in vivo SANDIX framework can account for the water exchange effect and provide measures of intra-soma and intra-neurite signal fractions without spurious time-dependence in estimated parameters, whereas the measured water exchange times need to be interpreted with caution. The model was then applied to in vivo diffusion MRI data acquired in 13 healthy adults on the 3-Tesla Connectome MRI scanner equipped with 300 mT/m gradients. The in vivo results exhibited patterns that were consistent with corresponding anatomical characteristics in both cortex and white matter. In particular, the estimated water exchange times in gray and white matter were distinct and differentiated between the two tissue types. Our results show the SANDIX approach applied to high-gradient diffusion MRI data achieves cortical microstructure mapping of the in vivo human brain with the evaluation of water exchange effects. This approach potentially provides a more appropriate description of in vivo cortical microstructure for improving data interpretation in future neurobiological studies.

Considerations and recommendations from the ISMRM diffusion study group for preclinical diffusion MRI: Part 1: In vivo small-animal imaging

Small-animal diffusion MRI (dMRI) has been used for methodological development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. The steps from animal setup and monitoring, to acquisition, analysis, and interpretation are complex, with many decisions that may ultimately affect what questions can be answered using the resultant data. This work aims to present selected considerations and recommendations from the diffusion community on best practices for preclinical dMRI of in vivo animals. We describe the general considerations and foundational knowledge that must be considered when designing experiments. We briefly describe differences in animal species and disease models and discuss why some may be more or less appropriate for different studies. We, then, give recommendations for in vivo acquisition protocols, including decisions on hardware, animal preparation, and imaging sequences, followed by advice for data processing including preprocessing, model-fitting, and tractography. Finally, we provide an online resource that lists publicly available preclinical dMRI datasets and software packages to promote responsible and reproducible research. In each section, we attempt to provide guides and recommendations, but also highlight areas for which no guidelines exist (and why), and where future work should focus. Although we mainly cover the central nervous system (on which most preclinical dMRI studies are focused), we also provide, where possible and applicable, recommendations for other organs of interest. An overarching goal is to enhance the rigor and reproducibility of small animal dMRI acquisitions and analyses, and thereby advance biomedical knowledge.

The Duke Mouse Brain Atlas: MRI and light sheet microscopy stereotaxic atlas of the mouse brain

Atlases of the brain are critical resources that make it possible to share data in a common reference frame. Unexpectedly, there is no three-dimensional (3D) stereotaxic atlas of the mouse brain that provides whole brain coverage at macro to single-cell levels. Diffusion tensor images from five perfusion-fixed (in skull) specimens were acquired at 15 micrometers, the highest resolution ever reported. Diffusion tensor imaging yields multiple 3D volumes, each of which highlights unique cytoarchitecture. The averages were mapped into micro-computed tomography of the mouse skull to create external landmarks (bregma and lambda). Light sheet images of the same brains were coregistered, providing cell maps in the same stereotaxic space. The Allen Reference Atlas was registered to the volume to correct the geometric distortion in that atlas and bring it into the stereotaxic space. The resulting multiscalar (13 terabytes) atlas provides a common spatial framework to anneal data across molecular, structural, and functional studies of mice.

MRI signatures of cortical microstructure in human development align with oligodendrocyte cell-type expression

Neuroanatomical changes to the cortex during adolescence have been well documented using MRI, revealing ongoing cortical thinning and volume loss. Recent advances in MRI hardware and biophysical models of tissue informed by diffusion MRI data hold promise for identifying the cellular changes driving these morphological observations. Using ultra-strong gradient MRI, this study quantifies cortical neurite and soma microstructure in typically developing youth. Across domain-specific networks, cortical neurite signal fraction, attributed to neuronal and glial processes, increases with age. The apparent soma radius, attributed to the apparent radius of glial and neuronal cell bodies, decreases with age. Analyses of two independent post-mortem datasets reveal that genes increasing in expression through adolescence are significantly enriched in cortical oligodendrocytes and Layer 5-6 neurons. In our study, we show spatial and temporal alignment of oligodendrocyte cell-type gene expression with neurite and soma microstructural changes, suggesting that ongoing cortical myelination processes drive adolescent cortical development.

Soma and neurite density imaging detects brain microstructural impairments in amyotrophic lateral sclerosis

OBJECTIVE

To investigate whole-brain microstructural changes in amyotrophic lateral sclerosis (ALS) using soma and neurite density imaging (SANDI), a novel multicompartment model of diffusion-weighted imaging that estimates apparent soma and neurite density.

METHODS

This study consists of 41 healthy controls and 43 patients with ALS, whose diffusion-weighted data were acquired. The SANDI-derived (including signal fractions of soma (fsoma), neurite (fneurite), and extra-cellular space (fextra)) and diffusion tensor imaging (DTI)-derived metrics were obtained. Voxel-based analyses were performed to evaluate intergroup differences and the correlation of SANDI and DTI metrics with clinical parameters.

RESULTS

In ALS patients, fneurite reduction involved both gray matter (primarily the bilateral precentral gyri, supplementary motor area, medial frontal gyrus, anterior cingulate cortex, inferior frontal gyrus, orbital gyrus, paracentral lobule, postcentral gyrus, middle cingulate cortex, hippocampus and parahippocampal gyrus, and insula, and left anterior parts of the temporal lobe) and white matter (primarily the bilateral corticospinal tract, body of corpus callosum, and brainstem) (P <0.05 after false discovery rate correction). The fextra increment showed a similar spatial distribution in ALS patients. Interestingly, the decreased fsoma in ALS primarily located in gray matter; while, the increased fsoma primarily involved white matter. The spatial distribution of fneurite/fextra/fsoma changes was larger than that detected by conventional DTI metrics, and the fneurite/fextra/fsoma were correlated with disease severity.

CONCLUSIONS

SANDI may serve as a clinically relevant model, superior to conventional DTI, for characterizing microstructural impairments such as neurite degeneration and soma alteration in ALS.

Soma and neurite density abnormalities of paramagnetic rim lesions and core-sign lesions in multiple sclerosis

BACKGROUND

In multiple sclerosis (MS), susceptibility-weighted imaging (SWI) may reveal white matter lesions (WML) with a paramagnetic rim ("paramagnetic rim lesions" [PRLs]) or diffuse hypointensity ("core-sign lesions"), reflecting different stages of WML evolution.

OBJECTIVE

Using the soma and neurite density imaging (SANDI) model on diffusion-weighted magnetic resonance imaging (MRI), we characterized microstructural abnormalities of MS PRLs and core-sign lesions and their clinical relevance.

METHODS

Forty MS patients and 20 healthy controls (HC) underwent a 3 T brain MRI. Using SANDI, the fractions of neurite (fneurite) and soma (fsoma) and size of soma (rsoma) were quantified in PRLs (including their core and rim separately), and core-sign lesions identified on SWI-phase.

RESULTS

Among 1811 WMLs, 122 (6.7%) core-sign lesions and 97 (5.4%) PRLs were identified. Compared to HC and MS normal-appearing white matter, all MS WML showed significantly lower fneurite and fsoma and higher rsoma (FDR-p < 0.001). Compared to SWI-isointense WML, core-sign lesions showed a significantly higher fneurite, and lower fsoma and rsoma (FDR-p ≤ 0.005). Compared to SWI-isointense WML and core-sign lesions, PRLs showed a significantly lower fneurite, higher fsoma, and higher rsoma (FDR-p ≤ 0.001). The PRL-core showed significantly lower fneurite, and higher rsoma than PRL-rim (FDR-p < 0.001). Lower PRL fneurite (β ≤ -0.006, FDR-p ≤ 0.015) and higher rsoma (β ≥ 0.032, FDR-p ≤ 0.024) were significantly associated with a longer disease duration and more severe disability.

CONCLUSIONS

In PRLs, the significant and clinically relevant neurite loss and increased soma fraction and size possibly reflect increased astrogliosis and activated microglia. Core-sign lesions exhibit milder axonal loss, microglia density and astrogliosis, supporting their less destructive nature.

In vivo human neurite exchange imaging (NEXI) at 500 mT/m diffusion gradients

Evaluating tissue microstructure and membrane integrity in the living human brain through diffusion-water exchange imaging is challenging due to requirements for a high signal-to-noise ratio and short diffusion times dictated by relatively fast exchange processes. The goal of this work was to demonstrate the feasibility of in vivo imaging of tissue micro-geometries and water exchange within the brain gray matter using the state-of-the-art Connectome 2.0 scanner equipped with an ultra-high-performance gradient system (maximum gradient strength=500 mT/m, maximum slew rate=600 T/m/s). We performed diffusion MRI measurements in 15 healthy volunteers at multiple diffusion times (13-30 ms) and b -values up to 17.5 ms/μm2. The anisotropic Kärger model was applied to estimate the exchange time between intra-neurite and extracellular water in gray matter. The estimated exchange time across the cortical ribbon was around (median±interquartile range) 13±8 ms on Connectome 2.0, substantially faster than that measured using an imaging protocol compatible with Connectome 1.0-alike systems on the same cohort. Our investigation suggested that the NEXI exchange time estimation using a Connectome 1.0 compatible protocol was more prone to residual noise floor biases due to the small time-dependent signal contrasts across diffusion times when the exchange is fast (≤20 ms). Furthermore, spatial variation of exchange time was observed across the cortex, where the motor cortex, somatosensory cortex and visual cortex exhibit longer exchange times compared to other cortical regions. Non-linear fitting for the anisotropic Kärger model was accelerated 100 times using a GPU-based pipeline compared to the conventional CPU-based approach. This study highlighted the importance of the chosen diffusion times and measures to address Rician noise in dMRI data, which can have a substantial impact on the estimated NEXI exchange time and require extra attention when comparing NEXI results between various hardware setups.

Age-related alterations in human cortical microstructure across the lifespan: Insights from high-gradient diffusion MRI

The human brain undergoes age-related microstructural alterations across the lifespan. Soma and Neurite Density Imaging (SANDI), a novel biophysical model of diffusion MRI, provides estimates of cell body (soma) radius and density, and neurite density in gray matter. The goal of this cross-sectional study was to assess the sensitivity of high-gradient diffusion MRI toward age-related alterations in cortical microstructure across the adult lifespan using SANDI. Seventy-two cognitively unimpaired healthy subjects (ages 19-85 years; 40 females) were scanned on the 3T Connectome MRI scanner with a maximum gradient strength of 300mT/m using a multi-shell diffusion MRI protocol incorporating 8 b-values and diffusion time of 19 ms. Intra-soma signal fraction obtained from SANDI model-fitting to the data was strongly correlated with age in all major cortical lobes (r = -0.69 to -0.60, FDR-p < 0.001). Intra-soma signal fraction (r = 0.48-0.63, FDR-p < 0.001) and soma radius (r = 0.28-0.40, FDR-p < 0.04) were significantly correlated with cortical volume in the prefrontal cortex, frontal, parietal, and temporal lobes. The strength of the relationship between SANDI metrics and age was greater than or comparable to the relationship between cortical volume and age across the cortical regions, particularly in the occipital lobe and anterior cingulate gyrus. In contrast to the SANDI metrics, all associations between diffusion tensor imaging (DTI) and diffusion kurtosis imaging metrics and age were low to moderate. These results suggest that high-gradient diffusion MRI may be more sensitive to underlying substrates of neurodegeneration in the aging brain than DTI and traditional macroscopic measures of neurodegeneration such as cortical volume and thickness.

The developing hippocampus: Microstructural evolution through childhood and adolescence

The hippocampus is a structure in the medial temporal lobe which serves multiple cognitive functions. While important, the development of the hippocampus in the formative period of childhood and adolescence has not been extensively investigated, with most contemporary research focusing on macrostructural measures of volume. Thus, there has been little research on the development of the micron-scale structures (i.e., microstructure) of the hippocampus, which engender its cognitive functions. The current study examined age-related changes of hippocampal microstructure using diffusion MRI data acquired with an ultra-strong gradient (300 mT/m) MRI scanner in a sample of children and adolescents (N=88; 8-19 years). Surface-based hippocampal modelling was combined with established microstructural approaches, such as Diffusion Tensor Imaging (DTI) and Neurite Orientation Dispersion Density Imaging (NODDI), and a more advanced gray matter diffusion model Soma And Neurite Density Imaging (SANDI). No significant changes in macrostructural measures (volume, gyrification, and thickness) were found between 8-19 years, while significant changes in microstructure measures related to neurites (from NODDI and SANDI), soma (from SANDI), and mean diffusivity (from DTI) were found. In particular, there was a significant increase across age in neurite MR signal fraction and a significant decrease in extracellular MR signal fraction and mean diffusivity across the hippocampal subfields and long-axis. A significant negative correlation between age and MR apparent soma radius was found in the subiculum and CA1 throughout the anterior and body of the hippocampus. Further surface-based analyses uncovered variability in age-related microstructural changes between the subfields and long-axis, which may reflect ostensible developmental differences along these two axes. Finally, correlation of hippocampal surfaces representing age-related changes of microstructure with maps derived from histology allowed for postulation of the potential underlying microstructure that diffusion changes across age may be capturing. Overall, distinct neurite and soma developmental profiles in the human hippocampus during late childhood and adolescence are reported for the first time.

MRI signatures of cortical microstructure in human development align with oligodendrocyte cell-type expression

Neuroanatomical changes to the cortex during adolescence have been well documented using MRI, revealing ongoing cortical thinning and volume loss with age. However, the underlying cellular mechanisms remain elusive with conventional neuroimaging. Recent advances in MRI hardware and new biophysical models of tissue informed by diffusion MRI data hold promise for identifying the cellular changes driving these morphological observations. This study used ultra-strong gradient MRI to obtain high-resolution, in vivo estimates of cortical neurite and soma microstructure in sample of typically developing children and adolescents. Cortical neurite signal fraction, attributed to neuronal and glial processes, increased with age (mean R2 fneurite=.53, p<3.3e-11, 11.91% increase over age), while apparent soma radius decreased (mean R2 Rsoma=.48, p<4.4e-10, 1% decrease over age) across domain-specific networks. To complement these findings, developmental patterns of cortical gene expression in two independent post-mortem databases were analysed. This revealed increased expression of genes expressed in oligodendrocytes, and excitatory neurons, alongside a relative decrease in expression of genes expressed in astrocyte, microglia and endothelial cell-types. Age-related genes were significantly enriched in cortical oligodendrocytes, oligodendrocyte progenitors and Layer 5-6 neurons (pFDR<.001) and prominently expressed in adolescence and young adulthood. The spatial and temporal alignment of oligodendrocyte cell-type gene expression with neurite and soma microstructural changes suggest that ongoing cortical myelination processes contribute to adolescent cortical development. These findings highlight the role of intra-cortical myelination in cortical maturation during adolescence and into adulthood.

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.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

Human evolution has seen the development of higher-order cognitive and social capabilities in conjunction with the unique laminar cytoarchitecture of the human cortex. Moreover, early-life cortical maldevelopment has been associated with various neurodevelopmental diseases. Despite these connections, there is currently no noninvasive technique available for imaging the detailed cortical laminar structure. This study aims to address this scientific and clinical gap by introducing an approach for imaging human cortical lamina. This method combines diffusion-relaxation multidimensional MRI with a tailored unsupervised machine learning approach that introduces enhanced microstructural sensitivity. This new imaging method simultaneously encodes the microstructure, the local chemical composition and importantly their correlation within complex and heterogenous tissue. To validate our approach, we compared the intra-cortical layers obtained using our ex vivo MRI-based method with those derived from Nissl staining of postmortem human brain specimens. The integration of unsupervised learning with diffusion-relaxation correlation MRI generated maps that demonstrate sensitivity to areal differences in cytoarchitectonic features observed in histology. Significantly, our observations revealed layer-specific diffusion-relaxation signatures, showing reductions in both relaxation times and diffusivities at the deeper cortical levels. These findings suggest a radial decrease in myelin content and changes in cell size and anisotropy, reflecting variations in both cytoarchitecture and myeloarchitecture. Additionally, we demonstrated that 1D relaxation and high-order diffusion MRI scalar indices, even when aggregated and used jointly in a multimodal fashion, cannot disentangle the cortical layers. Looking ahead, our technique holds the potential to open new avenues of research in human neurodevelopment and the vast array of disorders caused by disruptions in neurodevelopment.

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.

Detection of grey matter microstructural substrates of neurodegeneration in multiple sclerosis

Multiple sclerosis features complex pathological changes in grey matter that begin early and eventually lead to diffuse atrophy. Novel approaches to image grey-matter microstructural alterations in vivo are highly sought after and would enable more sensitive monitoring of disease activity and progression. This cross-sectional study aimed to assess the sensitivity of high-gradient diffusion MRI for microstructural tissue damage in cortical and deep grey matter in people with multiple sclerosis and test the hypothesis that reduced cortical cell body density is associated with cortical and deep grey-matter volume loss. Forty-one people with multiple sclerosis (age 24-72, 14 females) and 37 age- and sex-matched healthy controls were scanned on a 3 T Connectom MRI scanner equipped with 300 mT/m gradients using a multi-shell diffusion MRI protocol. The soma and neurite density imaging model was fitted to high-gradient diffusion MRI data to obtain estimates of intra-neurite, intra-cellular and extra-cellular signal fractions and apparent soma radius. Cortical and deep grey-matter microstructural imaging metrics were compared between multiple sclerosis and healthy controls and correlated with grey-matter volume, clinical disability and cognitive outcomes. People with multiple sclerosis showed significant cortical and deep grey-matter volume loss compared with healthy controls. People with multiple sclerosis showed trends towards lower cortical intra-cellular signal fraction and significantly lower intra-cellular and higher extra-cellular signal fractions in deep grey matter, especially the thalamus and caudate, compared with healthy controls. Changes were most pronounced in progressive disease and correlated with the Expanded Disability Status Scale, but not the Symbol Digit Modalities Test. In multiple sclerosis, normalized thalamic volume was associated with thalamic microstructural imaging metrics. Whereas thalamic volume loss did not correlate with cortical volume loss, cortical microstructural imaging metrics were significantly associated with thalamic volume, and not with cortical volume. Compared with the short diffusion time (Δ = 19 ms) achievable on the Connectom scanner, at the longer diffusion time of Δ = 49 ms attainable on clinical scanners, multiple sclerosis-related changes in imaging metrics were generally less apparent with lower effect sizes in cortical and deep grey matter. Soma and neurite density imaging metrics obtained from high-gradient diffusion MRI data provide detailed grey-matter characterization beyond cortical and thalamic volumes and distinguish multiple sclerosis-related microstructural pathology from healthy controls. Cortical cell body density correlates with thalamic volume, appears sensitive to the microstructural substrate of neurodegeneration and reflects disability status in people with multiple sclerosis, becoming more pronounced as disability worsens.

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.

Tensor denoising of multidimensional MRI data

PURPOSE

To develop a denoising strategy leveraging redundancy in high-dimensional data.

THEORY AND METHODS

The SNR fundamentally limits the information accessible by MRI. This limitation has been addressed by a host of denoising techniques, recently including the so-called MPPCA: principal component analysis of the signal followed by automated rank estimation, exploiting the Marchenko-Pastur distribution of noise singular values. Operating on matrices comprised of data patches, this popular approach objectively identifies noise components and, ideally, allows noise to be removed without introducing artifacts such as image blurring, or nonlocal averaging. The MPPCA rank estimation, however, relies on a large number of noise singular values relative to the number of signal components to avoid such ill effects. This condition is unlikely to be met when data patches and therefore matrices are small, for example due to spatially varying noise. Here, we introduce tensor MPPCA (tMPPCA) for the purpose of denoising multidimensional data, such as from multicontrast acquisitions. Rather than combining dimensions in matrices, tMPPCA uses each dimension of the multidimensional data's inherent tensor-structure to better characterize noise, and to recursively estimate signal components.

RESULTS

Relative to matrix-based MPPCA, tMPPCA requires no additional assumptions, and comparing the two in a numerical phantom and a multi-TE diffusion MRI data set, tMPPCA dramatically improves denoising performance. This is particularly true for small data patches, suggesting that tMPPCA can be especially beneficial in such cases.

CONCLUSIONS

The MPPCA denoising technique can be extended to high-dimensional data with improved performance for smaller patch sizes.

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.

Mapping of magnetic resonance imaging's transverse relaxation time at low signal-to-noise ratio using Bloch simulations and principal component analysis image denoising

High-resolution mapping of magnetic resonance imaging (MRI)'s transverse relaxation time (T₂ ) can benefit many clinical applications by offering improved anatomic details, enhancing the ability to probe tissues' microarchitecture, and facilitating the identification of early pathology. Increasing spatial resolutions, however, decreases data's signal-to-noise ratio (SNR), particularly at clinical scan times. This impairs imaging quality, and the accuracy of subsequent radiological interpretation. Recently, principal component analysis (PCA) was employed for denoising diffusion-weighted MR images and was shown to be effective for improving parameter estimation in multiexponential relaxometry. This study combines the Marchenko-Pastur PCA (MP-PCA) signal model with the echo modulation curve (EMC) algorithm for denoising multiecho spin-echo (MESE) MRI data and improving the precision of EMC-generated single T₂ relaxation maps. The denoising technique was validated on simulations, phantom scans, and in vivo brain and knee data. MESE scans were performed on a 3-T Siemens scanner. The acquired images were denoised using the MP-PCA algorithm and were then provided as input for the EMC T₂ -fitting algorithm. Quantitative analysis of the denoising quality included comparing the standard deviation and coefficient of variation of T₂ values, along with gold standard SNR estimation of the phantom scans. The presented denoising technique shows an increase in T₂ maps' precision and SNR, while successfully preserving the morphological features of the tissue. Employing MP-PCA denoising as a preprocessing step decreases the noise-related variability of T₂ maps produced by the EMC algorithm and thus increases their precision. The proposed method can be useful for a wide range of clinical applications by facilitating earlier detection of pathologies and improving the accuracy of patients' follow-up.