Using double pulsed-field gradient MRI to study tissue microstructure in traumatic brain injury (TBI)

Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms

PURPOSE

We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a "powder-averaged" emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues.

METHODS

Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter.

RESULTS

Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation.

CONCLUSION

This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

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.

Double pulsed field gradient diffusion MRI to assess skeletal muscle microstructure

PURPOSE

Preliminary study to determine whether double pulsed field gradient (PFG) diffusion MRI is sensitive to key features of muscle microstructure related to function.

METHODS

The restricted diffusion profile of molecules in models of muscle microstructure derived from histology were systematically simulated using a numerical simulation approach. Diffusion tensor subspace imaging analysis of the diffusion signal was performed, and spherical anisotropy (SA) was calculated for each model. Linear regression was used to determine the predictive capacity of SA on the fiber area, fiber diameter, and surface area to volume ratio of the models. Additionally, a rat model of muscle hypertrophy was scanned using a single PFG and a double PFG pulse sequence, and the restricted diffusion measurements were compared with histological measurements of microstructure.

RESULTS

Excellent agreement between SA and muscle fiber area (r2  = 0.71; p < 0.0001), fiber diameter (r2  = 0.83; p < 0.0001), and surface area to volume ratio (r2  = 0.97; p < 0.0001) in simulated models was found. In a scanned rat leg, the distribution of these microstructural features measured from histology was broad and demonstrated that there is a wide variance in the microstructural features observed, similar to the SA distributions. However, the distribution of fractional anisotropy measurements in the same tissue was narrow.

CONCLUSIONS

This study demonstrates that SA-a scalar value from diffusion tensor subspace imaging analysis-is highly sensitive to muscle microstructural features predictive of function. Furthermore, these techniques and analysis tools can be translated to real experiments in skeletal muscle. The increased dynamic range of SA compared with fractional anisotropy in the same tissue suggests increased sensitivity to detecting changes in tissue microstructure.

Hyperpolarized 13C metabolic imaging detects long-lasting metabolic alterations following mild repetitive traumatic brain injury

Career athletes, active military, and head trauma victims are at increased risk for mild repetitive traumatic brain injury (rTBI), a condition that contributes to the development of epilepsy and neurodegenerative diseases. Standard clinical imaging fails to identify rTBI-induced lesions, and novel non-invasive methods are needed. Here, we evaluated if hyperpolarized 13C magnetic resonance spectroscopic imaging (HP 13C MRSI) could detect long-lasting changes in brain metabolism 3.5 months post-injury in a rTBI mouse model. Our results show that this metabolic imaging approach can detect changes in cortical metabolism at that timepoint, whereas multimodal MR imaging did not detect any structural or contrast alterations. Using Machine Learning, we further show that HP 13C MRSI parameters can help classify rTBI vs. Sham and predict long-term rTBI-induced behavioral outcomes. Altogether, our study demonstrates the potential of metabolic imaging to improve detection, classification and outcome prediction of previously undetected rTBI.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

Meningeal and visual pathway MRI analysis after single and repetitive Closed-Head Impact Model of Engineered Rotational Acceleration (CHIMERA)-induced disruption in male and female mice

The consequences of forceful rotational acceleration on the central nervous system are not fully understood. While traumatic brain injury (TBI) research primarily has focused on effects related to the brain parenchyma, reports of traumatic meningeal enhancement in TBI patients may possess clinical significance. The objective of this study was to evaluate the meninges and brain for changes in dynamic contrast enhancement (DCE) magnetic resonance imaging (MRI) following closed-head impact model of engineered rotational acceleration (CHIMERA)–induced cerebral insult. Adult male and female mice received one (1 × ; n = 19 CHIMERA, n = 19 Sham) or four (4 × one/day; n = 18 CHIMERA, n = 12 Sham) injuries. Each animal underwent three MRI scans: 1 week before injury, immediately after the final injury, and 1 week post-injury. Compared with baseline readings and measures in sham animals, meningeal DCE in males was increased after single impact and repetitive injury. In female mice, DCE was elevated relative to their baseline level after a single impact. One week after CHIMERA, the meningeal enhancement returned to below baseline for single injured male mice, but compared with uninjured mice remained elevated in both sexes in the multiple impact groups. Pre-DCE meningeal T2-weighted relaxation time was increased only after 1 × CHIMERA in injured mice. Since vision is impaired after CHIMERA, visual pathway regions were analyzed through imaging and glial fibrillary acidic protein (GFAP) histology. Initial DCE in the lateral geniculate nucleus (LGN) and superior colliculus (SC) and T2 increases in the optic tract (OPT) and LGN were observed after injury with decreases in DCE and T2 1 week later. Astrogliosis was apparent in the OPT and SC with increased GFAP staining 7 days post-injury. To our knowledge, this is the first study to examine meningeal integrity after CHIMERA in both male and female rodents. DCE-MRI may serve as a useful approach for pre-clinical models of meningeal injury that will enable further evaluation of the underlying mechanisms.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

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.

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.

The closed-head impact model of engineered rotational acceleration (CHIMERA) as an application for traumatic brain injury pre-clinical research: A status report

Closed-head traumatic brain injury (TBI) is a worldwide concern with increasing prevalence and cost to society. Rotational acceleration is a primary mechanism in TBI that results from tissue strains that give rise to diffuse axonal injury. The Closed-Head Impact Model of Engineered Rotational Acceleration (CHIMERA) was recently introduced as a method for the study of impact acceleration effects in pre-clinical TBI research. This review provides a survey of the published literature implementing the CHIMERA device and describes pathological, imaging, neurophysiological, and behavioral findings. Findings show CHIMERA inflicts damage in white matter tracts as a key area of injury. Behaviorally, repeated studies have shown motor deficits and more chronic cognitive effects after CHIMERA injury. Good progress with model application has been accomplished by investigators attending to what is required for model validation. However, the majority of CHIMERA studies only utilize adult male mice. To further establish this model, more work with female animals and various age groups need to be performed, as well as studies to further establish and standardize methodologies for validation of the models for clinical relevance. Common data elements to standardize the reporting methodology for the CHIMERA literature are suggested.

Diffusion MR imaging acquisition and analytics for microstructural delineation in preclinical models of TBI

Significant progress has been made toward improving both the acquisition of clinical diffusion-weighted imaging (DWI) data and its analysis in the uninjured brain, through various techniques including a large number of model-based solutions that have been proposed to fit for multiple tissue compartments, and multiple fibers per voxel. While some of these techniques have been applied to clinical traumatic brain injury (TBI) research, the majority of these technological enhancements have yet to be fully implemented in the preclinical arena of TBI animal model-based research. In this review, we describe the requirement for preclinical, MRI-based efforts to provide systematic confirmation of the applicability of some of these models as indicators of tissue pathology within the injured brain. We review how current DWI techniques are currently being used in animal TBI models, and describe how both acquisition and analytic techniques could be extended to leverage the progress made in clinical work. Finally, we highlight remaining gaps in the preclinical pipeline from data acquisition to final analysis that currently have no real, preclinical-based correlate.