JEDI: Joint Estimation Diffusion Imaging of macroscopic and microscopic tissue properties

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.

Active Cortex Tractography

Most existing diffusion tractography algorithms are affected by gyral bias, causing the termination of streamlines at gyral crowns instead of sulcal banks. In this paper, we propose a tractography technique, called active cortex tractography (ACT), to overcome gyral bias by enabling fiber streamlines to curve naturally into the cortex. We show that the cortex can play an active role in cortical tractography by providing anatomical information to overcome orientation ambiguities as the streamlines enter the superficial white matter in gyral blades and approach the cortex. This is achieved by devising a direction scouting mechanism that takes into account the white matter surface normal vectors. The scouting mechanism allows probing of directions further in space to prepare the streamlines to turn at appropriate angles. The surface normal vectors guide the streamlines to turn into the cortex, perpendicular to the white-gray matter interface. Evaluation using synthetic, macaque and human data with different streamline seeding schemes demonstrates that ACT improves cortical tractography.

Brain Waves: Emergence of Localized, Persistent, Weakly Evanescent Cortical Loops

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of nonevanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full-brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the commonly agreed paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation because of the cortex wave modes in the highly folded areas will exhibit no apparent correlation with the fiber directions. Nonlinear coupling of those linear weakly evanescent wave modes then provides a universal mechanism for the emergence of synchronized brain wave field activity. The resonant and nonresonant terms of nonlinear coupling between multiple modes produce both synchronous spiking-like high-frequency wave activity as well as low-frequency wave rhythms. Numerical simulation of forced multiple-mode dynamics shows that, as forcing increases, there is a transition from damped to oscillatory regime that can then transition quickly to a nonoscillatory state when a critical excitation threshold is reached. The resonant nonlinear coupling results in the emergence of low-frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these weakly evanescent cortical wave modes have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and magnetoencephalography, and for the understanding of brain activity in general, including mechanisms of memory.