Faculty Bibliography

Diffusion Magnetic Resonance Imaging (dMRI) sensitises the MRI signal to spin motion. This includes Brownian diffusion, but also flow across intricate networks of capillaries. This effect, the intra-voxel incoherent motion (IVIM), enables microvasculature characterisation with dMRI, through metrics such as the vascular signal fraction f_V or the vascular Apparent Diffusion Coefficient (ADC) D^∗. The IVIM metrics, while sensitive to perfusion, are protocol-dependent, and their interpretation can change depending on the flow regime spins experience during the dMRI measurements (e.g., diffusive vs ballistic), which is in general not known for a given voxel. These facts hamper their practical clinical utility, and innovative vascular dMRI models are needed to enable the in vivo calculation of biologically meaningful markers of capillary flow. These could have relevant applications in cancer, as in the assessment of the response to anti-angiogenic therapies targeting tumour vessels. This paper tackles this need by introducing SpinFlowSim, an open-source simulator of dMRI signals arising from blood flow within pipe networks. SpinFlowSim, tailored for the laminar flow patterns within capillaries, enables the synthesis of highly-realistic microvascular dMRI signals, given networks reconstructed from histology. We showcase the simulator by generating synthetic signals for 15 networks, reconstructed from liver biopsies, and containing cancerous and non-cancerous tissue. Signals exhibit complex, non-mono-exponential behaviours, consistent with in vivo signal patterns, and pointing towards the co-existence of different flow regimes within the same network, as well as diffusion time dependence. We also demonstrate the potential utility of SpinFlowSim by devising a strategy for microvascular property mapping informed by the synthetic signals, and focussing on the quantification of blood velocity distribution moments and of an apparent network branching index. These were estimated in silico and in vivo, in healthy volunteers scanned at 1.5T and 3T and in 13 cancer patients, scanned at 1.5T. In conclusion, realistic flow simulations, as those enabled by SpinFlowSim, may play a key role in the development of the next-generation of dMRI methods for microvascular mapping, with immediate applications in oncology.

PURPOSE/OBJECTIVE: METHODS: RESULTS: CONCLUSION/CONCLUSIONS:Numerical simulations and theory provide an interpretation of a specific feature of kurtosis time-dependence, offering a potential biomarker for in vivo evaluation of pathology by disentangling the functional (permeability) and structural (cell size) integrity in tissues. This is relevant as the time-dependent diffusion cumulants are sensitive to pathological changes in membrane integrity and cellular structure in diseases, such as ischemic stroke, tumors, and Alzheimer's disease.

Joint modeling of diffusion and relaxation has seen growing interest due to its potential to provide complementary information about tissue microstructure. For brain white matter (WM), we designed an optimal diffusion-relaxometry MRI protocol that samples multiple b-values, B-tensor shapes, and echo times (TE). This variable-TE protocol (27 min) has as subsets a fixed-TE protocol (15 min) and a two-shell dMRI protocol (7 min), both characterizing diffusion only. We assessed the sensitivity, specificity, and reproducibility of these protocols with synthetic experiments and in six healthy volunteers. Compared with the fixed-TE protocol, the variable-TE protocol enables estimation of the free water fraction while also capturing compartmental

The corpus callosum (CC) is the most important interhemispheric white matter (WM) structure composed of several anatomically and functionally distinct WM tracts. Resolving these tracts is a challenge since the callosum appears relatively homogenous in conventional structural imaging. Commonly used callosal parcellation methods such as Hofer and Frahm scheme rely on rigid geometric guidelines to separate the substructures that are limited to consider individual variation. Here we present a novel subject-specific and microstructurally-informed method for callosal parcellation based on axonal water fraction (ƒ) known as a diffusion metric reflective of axon caliber and density. We studied 30 healthy subjects from the Human Connectome Project dataset with multi-shell diffusion MRI. The biophysical parameter ƒ was derived from compartment-specific WM modeling. Inflection points were identified where there were concavity changes in ƒ across the CC to delineate callosal subregions. We observed relatively higher ƒ in anterior and posterior areas known to consist of a greater number of small diameter fibers and lower ƒ in posterior body areas of the CC known to consist of a greater number of large diameter fibers. Based on the degree of change in ƒ along the callosum, seven callosal subregions were consistently delineated for each individual. Therefore, this method provides microstructurally informed callosal parcellation in a subject-specific way, allowing for more accurate analysis in the corpus callosum.

Noninvasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent advances in diffusion magnetic resonance imaging enabled in vivo examination of tissue microstructures well beyond the imaging resolution. Here, we proposed to use diffusion time-dependent diffusion kurtosis imaging (tDKI) to simultaneously assess cellular morphology and transmembrane permeability in hypoxic-ischemic (HI) brain injury. Through numerical simulations and organoid imaging, we demonstrated the feasibility of capturing effective size and permeability changes using tDKI. In vivo MRI of HI-injured mouse brains detected a shift of the tDKI peak to longer diffusion times, suggesting swelling of the cellular processes. Furthermore, we observed a faster decrease of the tDKI tail, reflecting increased transmembrane permeability associated with up-regulated water exchange or necrosis. Such information, unavailable from a single diffusion time, can predict salvageable tissues. Preliminary applications of tDKI in patients with ischemic stroke suggested increased transmembrane permeability in stroke regions, illustrating tDKI's potential for detecting pathological changes in the clinics.

BACKGROUND AND PURPOSE/OBJECTIVE:Because the corpus callosum connects the left and right hemispheres and a variety of WM bundles across the brain in complex ways, damage to the neighboring WM microstructure may specifically disrupt interhemispheric communication through the corpus callosum following mild traumatic brain injury. Here we use a mediation framework to investigate how callosal interhemispheric communication is affected by WM microstructure in mild traumatic brain injury. MATERIALS AND METHODS/METHODS:Multishell diffusion MR imaging was performed on 23 patients with mild traumatic brain injury within 1 month of injury and 17 healthy controls, deriving 11 diffusion metrics, including DTI, diffusional kurtosis imaging, and compartment-specific standard model parameters. Interhemispheric processing speed was assessed using the interhemispheric speed of processing task (IHSPT) by measuring the latency between word presentation to the 2 hemivisual fields and oral word articulation. Mediation analysis was performed to assess the indirect effect of neighboring WM microstructures on the relationship between the corpus callosum and IHSPT performance. In addition, we conducted a univariate correlation analysis to investigate the direct association between callosal microstructures and IHSPT performance as well as a multivariate regression analysis to jointly evaluate both callosal and neighboring WM microstructures in association with IHSPT scores for each group. RESULTS:Several significant mediators in the relationships between callosal microstructure and IHSPT performance were found in healthy controls. However, patients with mild traumatic brain injury appeared to lose such normal associations when microstructural changes occurred compared with healthy controls. CONCLUSIONS:This study investigates the effects of neighboring WM microstructure on callosal interhemispheric communication in healthy controls and patients with mild traumatic brain injury, highlighting that neighboring noncallosal WM microstructures are involved in callosal interhemispheric communication and information transfer. Further longitudinal studies may provide insight into the temporal dynamics of interhemispheric recovery following mild traumatic brain injury.

PURPOSE/OBJECTIVE:To determine whether the spatial scale and magnetic susceptibility of microstructure can be evaluated robustly from the decay of gradient-echo and spin-echo signals. THEORY AND METHODS/METHODS:Gradient-echo and spin-echo images were acquired from suspensions of spherical polystyrene microbeads of 10, 20, and 40 μm nominal diameter. The sizes of the beads and their magnetic susceptibility relative to the medium were estimated from the signal decay curves, using a lookup table generated from Monte Carlo simulations and an analytic model based on the Gaussian phase approximation. RESULTS:Fitting Monte Carlo predictions to spin-echo data yielded acceptable estimates of microstructural parameters for the 20 and 40 μm microbeads. Using gradient-echo data, the Monte Carlo lookup table provided satisfactory parameter estimates for the 20 μm beads but unstable results for the diameter of the largest beads. Neither spin-echo nor gradient-echo data allowed accurate parameter estimation for the smallest beads. The analytic model performed poorly over all bead sizes. CONCLUSIONS:Microstructural sources of magnetic susceptibility produce distinctive non-exponential signatures in the decay of gradient-echo and spin-echo signals. However, inverting the problem to extract microstructural parameters from the signals is nontrivial and, in certain regimes, ill-conditioned. For microstructure with small characteristic length scales, parameter estimation is hampered by the difficulty of acquiring accurate data at very short echo times. For microstructure with large characteristic lengths, the gradient-echo signal approaches the static-dephasing regime, where it becomes insensitive to size. Applicability of the analytic model was further limited by failure of the Gaussian phase approximation for all but the smallest beads.

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings

The corpus callosum (CC) is the most important interhemispheric white matter (WM) structure composed of several anatomically and functionally distinct WM tracts. Resolving these tracts is a challenge since the callosum appears relatively homogenous in conventional structural imaging. Commonly used callosal parcellation methods such as the Hofer/Frahm scheme rely on rigid geometric guidelines to separate the substructures that are limited to consider individual variation. Here we present a novel subject-specific and microstructurally-informed method for callosal parcellation based on axonal water fraction (ƒ) known as a diffusion metric reflective of axon caliber and density. We studied 30 healthy subjects from the Human Connectome Project (HCP) dataset with multi-shell diffusion MRI. The biophysical parameter ƒ was derived from compartment-specific WM modeling. Inflection points were identified where there were concavity changes in ƒ across the CC to delineate callosal subregions. We observed relatively higher ƒ in anterior and posterior areas consisting of a greater number of small diameter fibers and lower ƒ in posterior body areas of the CC consisting of a greater number of large diameter fibers. Based on degree of change in ƒ along the callosum, seven callosal subregions can be consistently delineated for each individual. We observe that ƒ can capture differences in underlying tissue microstructures and seven subregions can be identified across CC. Therefore, this method provides microstructurally informed callosal parcellation in a subject-specific way, allowing for more accurate analysis in the corpus callosum.

Nuclear magnetic resonance (NMR) has been instrumental in deciphering the structure of proteins. Here we show that transverse NMR relaxation, through its time-dependent relaxation rate, is distinctly sensitive to the structure of complex materials or biological tissues at the mesoscopic scale, from micrometers to tens of micrometers. Based on the ideas of universality, we show analytically and numerically that the time-dependent transverse relaxation rate approaches its long-time limit in a power-law fashion, with the dynamical exponent reflecting the universality class of mesoscopic magnetic structure. The spectral line shape acquires the corresponding non-analytic power law singularity at zero frequency. We experimentally detect the change in the dynamical exponent as a result of the transition into maximally random jammed state characterized by hyperuniform correlations. The relation between relaxational dynamics and magnetic structure opens the way for noninvasive characterization of porous media, complex materials and biological tissues.