Retaining information from multidimensional correlation MRI using a spectral regions of interest generator

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.

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 astrogliosis in the individual human brain using multidimensional MRI

There are currently no non-invasive imaging methods available for astrogliosis assessment or mapping in the central nervous system despite its essential role in the response to many disease states, such as infarcts, neurodegenerative conditions, traumatic brain injury and infection. Multidimensional MRI is an increasingly employed imaging modality that maximizes the amount of encoded chemical and microstructural information by probing relaxation (T1 and T2) and diffusion mechanisms simultaneously. Here, we harness the exquisite sensitivity of this imagining modality to derive a signature of astrogliosis and disentangle it from normative brain at the individual level using machine learning. We investigated ex vivo cerebral cortical tissue specimens derived from seven subjects who sustained blast-induced injuries, which resulted in scar-border forming astrogliosis without being accompanied by other types of neuropathological abnormality, and from seven control brain donors. By performing a combined post-mortem radiology and histopathology correlation study we found that astrogliosis induces microstructural and chemical changes that are robustly detected with multidimensional MRI, and which can be attributed to astrogliosis because no axonal damage, demyelination or tauopathy were histologically observed in any of the cases in the study. Importantly, we showed that no one-dimensional T1, T2 or diffusion MRI measurement can disentangle the microscopic alterations caused by this neuropathology. Based on these findings, we developed a within-subject anomaly detection procedure that generates MRI-based astrogliosis biomarker maps ex vivo, which were significantly and strongly correlated with co-registered histological images of increased glial fibrillary acidic protein deposition (r = 0.856, P < 0.0001; r = 0.789, P < 0.0001; r = 0.793, P < 0.0001, for diffusion-T2, diffusion-T1 and T1-T2 multidimensional data sets, respectively). Our findings elucidate the underpinning of MRI signal response from astrogliosis, and the demonstrated high spatial sensitivity and specificity in detecting reactive astrocytes at the individual level, and if reproduced in vivo, will significantly impact neuroimaging studies of injury, disease, repair and aging, in which astrogliosis has so far been an invisible process radiologically.

Abnormal brain gray matter volume in patients with major depressive disorder: Associated with childhood trauma?

BACKGROUND

Patients with major depressive disorders (MDD) have abnormalities in the frontal-limbic structures of the brain. Childhood trauma is a risk factor for both structural brain alterations and MDD. However, the relationships among the three have not been confirmed.

METHODS

Sixty-four patients with MDD and sixty-one healthy controls (HC) were checked by using MRI, the Hamilton Depression Scale (HAMD) and Childhood Trauma Questionnaire (CTQ). Voxel-based morphometry (VBM) was used to compare gray matter volume (GMV) differences between the two groups. Moreover, partial correlation and mediation analyses were conducted to test for potential associations between CTQ scores, different GMV, and clinical variables.

RESULTS

Compared to the HC group, the MDD patients showed decreased GMV in the right middle frontal gyrus (rMFG) and right precentral gyrus (rPreCG). In the patient group, reduced GMV in rMFG was associated with CTQ scores (r = -0.30, P = 0.019) and HAMD scores (r = -0.53, P < 0.001). Finally, in the patient group, mediation analysis revealed that the indirect effect of rMFG GMV in CTQ scores and HAMD scores was 0.115 and the proportion of indirect effect to total effect was 23.86%.

LIMITATIONS

This study used a cross-sectional collection, and it is unclear whether at the longitudinal level the brain GMV mediates the relationship between childhood trauma and depression.

CONCLUSIONS

Abnormalities in the frontal GMV were presented in the MDD patients. It is possible that childhood traumatic experiences cause inefficient GMV and ultimately lead to an increased susceptibility to depression.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

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.

Quantitative diffusion MRI of the abdomen and pelvis

Diffusion MRI has enormous potential and utility in the evaluation of various abdominal and pelvic disease processes including cancer and noncancer imaging of the liver, prostate, and other organs. Quantitative diffusion MRI is based on acquisitions with multiple diffusion encodings followed by quantitative mapping of diffusion parameters that are sensitive to tissue microstructure. Compared to qualitative diffusion-weighted MRI, quantitative diffusion MRI can improve standardization of tissue characterization as needed for disease detection, staging, and treatment monitoring. However, similar to many other quantitative MRI methods, diffusion MRI faces multiple challenges including acquisition artifacts, signal modeling limitations, and biological variability. In abdominal and pelvic diffusion MRI, technical acquisition challenges include physiologic motion (respiratory, peristaltic, and pulsatile), image distortions, and low signal-to-noise ratio. If unaddressed, these challenges lead to poor technical performance (bias and precision) and clinical outcomes of quantitative diffusion MRI. Emerging and novel technical developments seek to address these challenges and may enable reliable quantitative diffusion MRI of the abdomen and pelvis. Through systematic validation in phantoms, volunteers, and patients, including multicenter studies to assess reproducibility, these emerging techniques may finally demonstrate the potential of quantitative diffusion MRI for abdominal and pelvic imaging applications.

Multidimensional MRI for Characterization of Subtle Axonal Injury Accelerated Using an Adaptive Nonlocal Multispectral Filter

Multidimensional MRI is an emerging approach that simultaneously encodes water relaxation (T₁ and T₂) and mobility (diffusion) and replaces voxel-averaged values with subvoxel distributions of those MR properties. While conventional (i.e., voxel-averaged) MRI methods cannot adequately quantify the microscopic heterogeneity of biological tissue, using subvoxel information allows to selectively map a specific T₁-T₂-diffusion spectral range that corresponds to a group of tissue elements. The major obstacle to the adoption of rich, multidimensional MRI protocols for diagnostic or monitoring purposes is the prolonged scan time. Our main goal in the present study is to evaluate the performance of a nonlocal estimation of multispectral magnitudes (NESMA) filter on reduced datasets to limit the total acquisition time required for reliable multidimensional MRI characterization of the brain. Here we focused and reprocessed results from a recent study that identified potential imaging biomarkers of axonal injury pathology from the joint analysis of multidimensional MRI, in particular voxelwise T₁-T₂ and diffusion-T₂ spectra in human Corpus Callosum, and histopathological data. We tested the performance of NESMA and its effect on the accuracy of the injury biomarker maps, relative to the co-registered histological reference. Noise reduction improved the accuracy of the resulting injury biomarker maps, while permitting data reduction of 35.7 and 59.6% from the full dataset for T₁-T₂ and diffusion-T₂ cases, respectively. As successful clinical proof-of-concept applications of multidimensional MRI are continuously being introduced, reliable and robust noise removal and consequent acquisition acceleration would advance the field towards clinically-feasible diagnostic multidimensional MRI protocols.

Diffuse axonal injury has a characteristic multidimensional MRI signature in the human brain

Axonal injury is a major contributor to the clinical symptomatology in patients with traumatic brain injury. Conventional neuroradiological tools, such as CT and MRI, are insensitive to diffuse axonal injury (DAI) caused by trauma. Diffusion tensor MRI parameters may change in DAI lesions; however, the nature of these changes is inconsistent. Multidimensional MRI is an emerging approach that combines T1, T2, and diffusion, and replaces voxel-averaged values with distributions, which allows selective isolation of specific potential abnormal components. By performing a combined post-mortem multidimensional MRI and histopathology study, we aimed to investigate T1-T2-diffusion changes linked to DAI and to define their histopathological correlates. Corpora callosa derived from eight subjects who had sustained traumatic brain injury, and three control brain donors underwent post-mortem ex vivo MRI at 7 T. Multidimensional, diffusion tensor, and quantitative T1 and T2 MRI data were acquired and processed. Following MRI acquisition, slices from the same tissue were tested for amyloid precursor protein (APP) immunoreactivity to define DAI severity. A robust image co-registration method was applied to accurately match MRI-derived parameters and histopathology, after which 12 regions of interest per tissue block were selected based on APP density, but blind to MRI. We identified abnormal multidimensional T1-T2, diffusion-T2, and diffusion-T1 components that are strongly associated with DAI and used them to generate axonal injury images. We found that compared to control white matter, mild and severe DAI lesions contained significantly larger abnormal T1-T2 component (P = 0.005 and P < 0.001, respectively), and significantly larger abnormal diffusion-T2 component (P = 0.005 and P < 0.001, respectively). Furthermore, within patients with traumatic brain injury the multidimensional MRI biomarkers differentiated normal-appearing white matter from mild and severe DAI lesions, with significantly larger abnormal T1-T2 and diffusion-T2 components (P = 0.003 and P < 0.001, respectively, for T1-T2; P = 0.022 and P < 0.001, respectively, for diffusion-T2). Conversely, none of the conventional quantitative MRI parameters were able to differentiate lesions and normal-appearing white matter. Lastly, we found that the abnormal T1-T2, diffusion-T1, and diffusion-T2 components and their axonal damage images were strongly correlated with quantitative APP staining (r = 0.876, P < 0.001; r = 0.727, P < 0.001; and r = 0.743, P < 0.001, respectively), while producing negligible intensities in grey matter and in normal-appearing white matter. These results suggest that multidimensional MRI may provide non-invasive biomarkers for detection of DAI, which is the pathological substrate for neurological disorders ranging from concussion to severe traumatic 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.

Multidimensional Diffusion Magnetic Resonance Imaging for Characterization of Tissue Microstructure in Breast Cancer Patients: A Prospective Pilot Study

Diffusion-weighted imaging is a non-invasive functional imaging modality for breast tumor characterization through apparent diffusion coefficients. Yet, it has so far been unable to intuitively inform on tissue microstructure. In this IRB-approved prospective study, we applied novel multidimensional diffusion (MDD) encoding across 16 patients with suspected breast cancer to evaluate its potential for tissue characterization in the clinical setting. Data acquired via custom MDD sequences was processed using an algorithm estimating non-parametric diffusion tensor distributions. The statistical descriptors of these distributions allow us to quantify tissue composition in terms of metrics informing on cell densities, shapes, and orientations. Additionally, signal fractions from specific cell types, such as elongated cells (bin1), isotropic cells (bin2), and free water (bin3), were teased apart. Histogram analysis in cancers and healthy breast tissue showed that cancers exhibited lower mean values of "size" (1.43 ± 0.54 × 10-3 mm2/s) and higher mean values of "shape" (0.47 ± 0.15) corresponding to bin1, while FGT (fibroglandular breast tissue) presented higher mean values of "size" (2.33 ± 0.22 × 10-3 mm2/s) and lower mean values of "shape" (0.27 ± 0.11) corresponding to bin3 (p < 0.001). Invasive carcinomas showed significant differences in mean signal fractions from bin1 (0.64 ± 0.13 vs. 0.4 ± 0.25) and bin3 (0.18 ± 0.08 vs. 0.42 ± 0.21) compared to ductal carcinomas in situ (DCIS) and invasive carcinomas with associated DCIS (p = 0.03). MDD enabled qualitative and quantitative evaluation of the composition of breast cancers and healthy glands.

Diffusion-PEPTIDE: Distortion- and blurring-free diffusion imaging with self-navigated motion-correction and relaxometry capabilities

PURPOSE

To implement the time-resolved relaxometry PEPTIDE technique into a diffusion acquisition to provide self-navigated, distortion- and blurring-free diffusion imaging that is robust to motion, while simultaneously providing T₂ and T 2 ∗ mapping.

THEORY AND METHODS

The PEPTIDE readout was implemented into a spin-echo diffusion acquisition, enabling reconstruction of a time-series of T₂ - and T 2 ∗ -weighted images, free from conventional echo planar imaging (EPI) distortion and blurring, for each diffusion-encoding. Robustness of PEPTIDE to motion and shot-to-shot phase variation was examined through a deliberate motion-corrupted diffusion experiment. Two diffusion-relaxometry in vivo brain protocols were also examined: (1)1 × 1 × 3 mm³ across 32 diffusion directions in 20 min, (2)1.5 × 1.5 × 3.0 mm³ across 6 diffusion-weighted images in 3.4 min. T₂ , T 2 ∗ , and diffusion parameter maps were calculated from these data. As initial exploration of the rich diffusion-relaxometry data content for use in multi-compartment modeling, PEPTIDE data were acquired of a gadolinium-doped asparagus phantom. These datasets contained two compartments with different relaxation parameters and different diffusion orientation properties, and T₂ relaxation variations across these diffusion directions were explored.

RESULTS

Diffusion-PEPTIDE showed the capability to provide high quality diffusion images and T₂ and T 2 ∗ maps from both protocols. The reconstructions were distortion-free, avoided potential resolution losses exceeding 100% in equivalent EPI acquisitions, and showed tolerance to nearly 30° of rotational motion. Expected variation in T₂ values as a function of diffusion direction was observed in the two-compartment asparagus phantom (P < .01), demonstrating potential to explore diffusion-PEPTIDE data for multi-compartment modeling.

CONCLUSIONS

Diffusion-PEPTIDE provides highly robust diffusion and relaxometry data and offers potential for future applications in diffusion-relaxometry multi-compartment modeling.

Direct and specific assessment of axonal injury and spinal cord microenvironments using diffusion correlation imaging

We describe a practical two-dimensional (2D) diffusion MRI framework to deliver specificity and improve sensitivity to axonal injury in the spinal cord. This approach provides intravoxel distributions of correlations of water mobilities in orthogonal directions, revealing sub-voxel diffusion components. Here we use it to investigate water diffusivities along axial and radial orientations within spinal cord specimens with confirmed, tract-specific axonal injury. First, we show using transmission electron microscopy and immunohistochemistry that tract-specific axonal beading occurs following Wallerian degeneration in the cortico-spinal tract as direct sequelae to closed head injury. We demonstrate that although some voxel-averaged diffusion tensor imaging (DTI) metrics are sensitive to this axonal injury, they are non-specific, i.e., they do not reveal an underlying biophysical mechanism of injury. Then we employ 2D diffusion correlation imaging (DCI) to improve discrimination of different water microenvironments by measuring and mapping the joint water mobility distributions perpendicular and parallel to the spinal cord axis. We determine six distinct diffusion spectral components that differ according to their microscopic anisotropy and mobility. We show that at the injury site a highly anisotropic diffusion component completely disappears and instead becomes more isotropic. Based on these findings, an injury-specific MR image of the spinal cord was generated, and a radiological-pathological correlation with histological silver staining % area was performed. The resulting strong and significant correlation (r=0.70,p < 0.0001) indicates the high specificity with which DCI detects injury-induced tissue alterations. We predict that the ability to selectively image microstructural changes following axonal injury in the spinal cord can be useful in clinical and research applications by enabling specific detection and increased sensitivity to injury-induced microstructural alterations. These results also encourage us to translate DCI to higher spatial dimensions to enable assessment of traumatic axonal injury, and possibly other diseases and disorders in the brain.

MTE-NODDI: Multi-TE NODDI for disentangling non-T2-weighted signal fractions from compartment-specific T2 relaxation times

Neurite orientation dispersion and density imaging (NODDI) has become a popular diffusion MRI technique for investigating microstructural alternations during brain development, maturation and aging in health and disease. However, the NODDI model of diffusion does not explicitly account for compartment-specific T2 relaxation and its model parameters are usually estimated from data acquired with a single echo time (TE). Thus, the NODDI-derived measures, such as the intra-neurite signal fraction, also known as the neurite density index, could be T2-weighted and TE-dependent. This may confound the interpretation of studies as one cannot disentangle differences in diffusion from those in T2 relaxation. To address this challenge, we propose a multi-TE NODDI (MTE-NODDI) technique, inspired by recent studies exploiting the synergy between diffusion and T2 relaxation. MTE-NODDI could give robust estimates of the non-T2-weighted signal fractions and compartment-specific T2 values, as demonstrated by both simulation and in vivo data experiments. Results showed that the estimated non-T2 weighted intra-neurite fraction and compartment-specific T2 values in white matter were consistent with previous studies. The T2-weighted intra-neurite fractions from the original NODDI were found to be overestimated compared to their non-T2-weighted estimates; the overestimation increases with TE, consistent with the reported intra-neurite T2 being larger than extra-neurite T2. Finally, the inclusion of the free water compartment reduces the estimation error in intra-neurite T2 in the presence of cerebrospinal fluid contamination. With the ability to disentangle non-T2-weighted signal fractions from compartment-specific T2 relaxation, MTE-NODDI could help improve the interpretability of future neuroimaging studies, especially those in brain development, maturation and aging.