Practical computation of the diffusion MRI signal based on Laplace eigenfunctions: permeable interfaces.

The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium such as brain tissue can be modeled by the Bloch-Torrey partial differential equation. The spatial integral of the solution of this equation in realistic geometry provides a gold-standard reference model for the diffusion MRI signal arising from different tissue micro-structures of interest. A closed form representation of this reference diffusion MRI signal, called matrix formalism, which makes explicit the link between the Laplace eigenvalues and eigenfunctions of the tissue geometry and its diffusion MRI signal, was derived 20 years ago. In addition, once the Laplace eigendecomposition has been computed and saved, the diffusion MRI signal can be calculated for arbitrary diffusion-encoding sequences and b-values at negligible additional cost. In a previous publication, we presented a simulation framework that we implemented inside the MATLAB-based diffusion MRI simulator SpinDoctor that efficiently computes the matrix formalism representation for biological cells subject to impermeable membrane boundary conditions. In this work, we extend our simulation framework to include geometries that contain permeable cell membranes. We describe the new computational techniques that allowed this generalization and we analyze the effects of the magnitude of the permeability coefficient on the eigendecomposition of the diffusion and Bloch-Torrey operators. This work is another step in bringing advanced mathematical tools and numerical method development to the simulation and modeling of diffusion MRI.

Diffusion MRI simulation of realistic neurons with SpinDoctor and the Neuron Module.

The diffusion MRI signal arising from neurons can be numerically simulated by solving the Bloch-Torrey partial differential equation. In this paper we present the Neuron Module that we implemented within the Matlab-based diffusion MRI simulation toolbox SpinDoctor. SpinDoctor uses finite element discretization and adaptive time integration to solve the Bloch-Torrey partial differential equation for general diffusion-encoding sequences, at multiple b-values and in multiple diffusion directions. In order to facilitate the diffusion MRI simulation of realistic neurons by the research community, we constructed finite element meshes for a group of 36 pyramidal neurons and a group of 29 spindle neurons whose morphological descriptions were found in the publicly available neuron repository NeuroMorpho.Org. These finite elements meshes range from having 15,163 nodes to 622,553 nodes. We also broke the neurons into the soma and dendrite branches and created finite elements meshes for these cell components. Through the Neuron Module, these neuron and cell components finite element meshes can be seamlessly coupled with the functionalities of SpinDoctor to provide the diffusion MRI signal attributable to spins inside neurons. We make these meshes and the source code of the Neuron Module available to the public as an open-source package. To illustrate some potential uses of the Neuron Module, we show numerical examples of the simulated diffusion MRI signals in multiple diffusion directions from whole neurons as well as from the soma and dendrite branches, and include a comparison of the high b-value behavior between dendrite branches and whole neurons. In addition, we demonstrate that the neuron meshes can be used to perform Monte-Carlo diffusion MRI simulations as well. We show that at equivalent accuracy, if only one gradient direction needs to be simulated, SpinDoctor is faster than a GPU implementation of Monte-Carlo, but if many gradient directions need to be simulated, there is a break-even point when the GPU implementation of Monte-Carlo becomes faster than SpinDoctor. Furthermore, we numerically compute the eigenfunctions and the eigenvalues of the Bloch-Torrey and the Laplace operators on the neuron geometries using a finite elements discretization, in order to give guidance in the choice of the space and time discretization parameters for both finite elements and Monte-Carlo approaches. Finally, we perform a statistical study on the set of 65 neurons to test some candidate biomakers that can potentially indicate the soma size. This preliminary study exemplifies the possible research that can be conducted using the Neuron Module.

Practical computation of the diffusion MRI signal of realistic neurons based on Laplace eigenfunctions.

The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium such as brain tissue can be modeled by the Bloch-Torrey partial differential equation. The spatial integral of the solution of this equation in realistic geometry provides a gold-standard reference model for the diffusion MRI signal arising from different tissue micro-structures of interest. A closed form representation of this reference diffusion MRI signal called matrix formalism, which makes explicit the link between the Laplace eigenvalues and eigenfunctions of the biological cell and its diffusion MRI signal, was derived 20 years ago. In addition, once the Laplace eigendecomposition has been computed and saved, the diffusion MRI signal can be calculated for arbitrary diffusion-encoding sequences and b-values at negligible additional cost. Up to now, this representation, though mathematically elegant, has not been often used as a practical model of the diffusion MRI signal, due to the difficulties of calculating the Laplace eigendecomposition in complicated geometries. In this paper, we present a simulation framework that we have implemented inside the MATLAB-based diffusion MRI simulator SpinDoctor that efficiently computes the matrix formalism representation for realistic neurons using the finite element method. We show that the matrix formalism representation requires a few hundred eigenmodes to match the reference signal computed by solving the Bloch-Torrey equation when the cell geometry originates from realistic neurons. As expected, the number of eigenmodes required to match the reference signal increases with smaller diffusion time and higher b-values. We also convert the eigenvalues to a length scale and illustrate the link between the length scale and the oscillation frequency of the eigenmode in the cell geometry. We give the transformation that links the Laplace eigenfunctions to the eigenfunctions of the Bloch-Torrey operator and compute the Bloch-Torrey eigenfunctions and eigenvalues. This work is another step in bringing advanced mathematical tools and numerical method development to the simulation and modeling of diffusion MRI.

SpinDoctor: A MATLAB toolbox for diffusion MRI simulation.

The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium can be modeled by the multiple compartment Bloch-Torrey partial differential equation. Under the assumption of negligible water exchange between compartments, the time-dependent apparent diffusion coefficient can be directly computed from the solution of a diffusion equation subject to a time-dependent Neumann boundary condition. This paper describes a publicly available MATLAB toolbox called SpinDoctor that can be used 1) to solve the Bloch-Torrey partial differential equation in order to simulate the diffusion magnetic resonance imaging signal; 2) to solve a diffusion partial differential equation to obtain directly the apparent diffusion coefficient; 3) to compare the simulated apparent diffusion coefficient with a short-time approximation formula. The partial differential equations are solved by P1 finite elements combined with built-in MATLAB routines for solving ordinary differential equations. The finite element mesh generation is performed using an external package called Tetgen. SpinDoctor provides built-in options of including 1) spherical cells with a nucleus; 2) cylindrical cells with a myelin layer; 3) an extra-cellular space enclosed either a) in a box or b) in a tight wrapping around the cells; 4) deformation of canonical cells by bending and twisting; 5) permeable membranes; Built-in diffusion-encoding pulse sequences include the Pulsed Gradient Spin Echo and the Oscillating Gradient Spin Echo. We describe in detail how to use the SpinDoctor toolbox. We validate SpinDoctor simulations using reference signals computed by the Matrix Formalism method. We compare the accuracy and computational time of SpinDoctor simulations with Monte-Carlo simulations and show significant speed-up of SpinDoctor over Monte-Carlo simulations in complex geometries. We also illustrate several extensions of SpinDoctor functionalities, including the incorporation of T2 relaxation, the simulation of non-standard diffusion-encoding sequences, as well as the use of externally generated geometrical meshes.

Incorporating interface permeability into the diffusion MRI signal representation while using impermeable Laplace eigenfunctions.

Objective. The complex-valued transverse magnetization due to diffusion-encoding magnetic field gradients acting on a permeable medium can be modeled by the Bloch-Torrey partial differential equation. The diffusion magnetic resonance imaging (MRI) signal has a representation in the basis of the Laplace eigenfunctions of the medium. However, in order to estimate the permeability coefficient from diffusion MRI data, it is desirable that the forward solution can be calculated efficiently for many values of permeability.Approach. In this paper we propose a new formulation of the permeable diffusion MRI signal representation in the basis of the Laplace eigenfunctions of the same medium where the interfaces are made impermeable.Main results.We proved the theoretical equivalence between our new formulation and the original formulation in the case that the full eigendecomposition is used. We validated our method numerically and showed promising numerical results when a partial eigendecomposition is used. Two diffusion MRI sequences were used to illustrate the numerical validity of our new method.Significance.Our approach means that the same basis (the impermeable set) can be used for all permeability values, which reduces the computational time significantly, enabling the study of the effects of the permeability coefficient on the diffusion MRI signal in the future.

A simulation-driven supervised learning framework to estimate brain microstructure using diffusion MRI.

We propose a framework to train supervised learning models on synthetic data to estimate brain microstructure parameters using diffusion magnetic resonance imaging (dMRI). Although further validation is necessary, the proposed framework aims to seamlessly incorporate realistic simulations into dMRI microstructure estimation. Synthetic data were generated from over 1,000 neuron meshes converted from digital neuronal reconstructions and linked to their neuroanatomical parameters (such as soma volume and neurite length) using an optimized diffusion MRI simulator that produces intracellular dMRI signals from the solution of the Bloch-Torrey partial differential equation. By combining random subsets of simulated neuron signals with a free diffusion compartment signal, we constructed a synthetic dataset containing dMRI signals and 40 tissue microstructure parameters of 1.45 million artificial brain voxels. To implement supervised learning models we chose multilayer perceptrons (MLPs) and trained them on a subset of the synthetic dataset to estimate some microstructure parameters, namely, the volume fractions of soma, neurites, and the free diffusion compartment, as well as the area fractions of soma and neurites. The trained MLPs perform satisfactorily on the synthetic test sets and give promising in-vivo parameter maps on the MGH Connectome Diffusion Microstructure Dataset (CDMD). Most importantly, the estimated volume fractions showed low dependence on the diffusion time, the diffusion time independence of the estimated parameters being a desired property of quantitative microstructure imaging. The synthetic dataset we generated will be valuable for the validation of models that map between the dMRI signals and microstructure parameters. The surface meshes and microstructures parameters of the aforementioned neurons have been made publicly available.

Three-dimensional micro-structurally informed in silico myocardium-Towards virtual imaging trials in cardiac diffusion weighted MRI.

In silico tissue models (viz. numerical phantoms) provide a mechanism for evaluating quantitative models of magnetic resonance imaging. This includes the validation and sensitivity analysis of imaging biomarkers and tissue microstructure parameters. This study proposes a novel method to generate a realistic numerical phantom of myocardial microstructure. The proposed method extends previous studies by accounting for the variability of the cardiomyocyte shape, water exchange between the cardiomyocytes (intercalated discs), disorder class of myocardial microstructure, and four sheetlet orientations. In the first stage of the method, cardiomyocytes and sheetlets are generated by considering the shape variability and intercalated discs in cardiomyocyte-cardiomyocyte connections. Sheetlets are then aggregated and oriented in the directions of interest. The morphometric study demonstrates no significant difference (p>0.01) between the distribution of volume, length, and primary and secondary axes of the numerical and real (literature) cardiomyocyte data. Moreover, structural correlation analysis validates that the in-silico tissue is in the same class of disorderliness as the real tissue. Additionally, the absolute angle differences between the simulated helical angle (HA) and input HA (reference value) of the cardiomyocytes (4.3°±3.1°) demonstrate a good agreement with the absolute angle difference between the measured HA using experimental cardiac diffusion tensor imaging (cDTI) and histology (reference value) reported by (Holmes et al., 2000) (3.7°±6.4°) and (Scollan et al. 1998) (4.9°±14.6°). Furthermore, the angular distance between eigenvectors and sheetlet angles of the input and simulated cDTI is much smaller than those between measured angles using structural tensor imaging (as a gold standard) and experimental cDTI. Combined with the qualitative results, these results confirm that the proposed method can generate richer numerical phantoms for the myocardium than previous studies.

Apparent diffusion coefficient measured by diffusion MRI of moving and deforming domains.

The modeling of the diffusion MRI signal from moving and deforming organs such as the heart is challenging due to significant motion and deformation of the imaged medium during the signal acquisition. Recently, a mathematical formulation of the Bloch-Torrey equation, describing the complex transverse magnetization due to diffusion-encoding magnetic field gradients, was developed to account for the motion and deformation. In that work, the motivation was to cancel the effect of the motion and deformation in the MRI image and the space scale of interest spans multiple voxels. In the present work, we adapt the mathematical equation to study the diffusion MRI signal at the much smaller scale of biological cells. We start with the Bloch-Torrey equation defined on a cell that is moving and deforming and linearize the equation around the magnitude of the diffusion-encoding gradient. The result is a second order signal model in which the linear term gives the imaginary part of the diffusion MRI signal and the quadratic term gives the apparent diffusion coefficient (ADC) attributable to the biological cell. We numerically validate this model for a variety of motions and deformations.

Microstructural organization of human insula is linked to its macrofunctional circuitry and predicts cognitive control.

The human insular cortex is a heterogeneous brain structure which plays an integrative role in guiding behavior. The cytoarchitectonic organization of the human insula has been investigated over the last century using postmortem brains but there has been little progress in noninvasive in vivo mapping of its microstructure and large-scale functional circuitry. Quantitative modeling of multi-shell diffusion MRI data from 413 participants revealed that human insula microstructure differs significantly across subdivisions that serve distinct cognitive and affective functions. Insular microstructural organization was mirrored in its functionally interconnected circuits with the anterior cingulate cortex that anchors the salience network, a system important for adaptive switching of cognitive control systems. Furthermore, insular microstructural features, confirmed in Macaca mulatta, were linked to behavior and predicted individual differences in cognitive control ability. Our findings open new possibilities for probing psychiatric and neurological disorders impacted by insular cortex dysfunction, including autism, schizophrenia, and fronto-temporal dementia.

Diffusion MRI simulation in thin-layer and thin-tube media using a discretization on manifolds.

The Bloch-Torrey partial differential equation can be used to describe the evolution of the transverse magnetization of the imaged sample under the influence of diffusion-encoding magnetic field gradients inside the MRI scanner. The integral of the magnetization inside a voxel gives the simulated diffusion MRI signal. This paper proposes a finite element discretization on manifolds in order to efficiently simulate the diffusion MRI signal in domains that have a thin layer or a thin tube geometrical structure. The variable thickness of the three-dimensional domains is included in the weak formulation established on the manifolds. We conducted a numerical study of the proposed approach by simulating the diffusion MRI signals from the extracellular space (a thin layer medium) and from neurons (a thin tube medium), comparing the results with the reference signals obtained using a standard three-dimensional finite element discretization. We show good agreements between the simulated signals using our proposed method and the reference signals for a wide range of diffusion MRI parameters. The approximation becomes better as the diffusion time increases. The method helps to significantly reduce the required simulation time, computational memory, and difficulties associated with mesh generation, thus opening the possibilities to simulating complicated structures at low cost for a better understanding of diffusion MRI in the brain.

Portable simulation framework for diffusion MRI.

The numerical simulation of the diffusion MRI signal arising from complex tissue micro-structures is helpful for understanding and interpreting imaging data as well as for designing and optimizing MRI sequences. The discretization of the Bloch-Torrey equation by finite elements is a more recently developed approach for this purpose, in contrast to random walk simulations, which has a longer history. While finite element discretization is more difficult to implement than random walk simulations, the approach benefits from a long history of theoretical and numerical developments by the mathematical and engineering communities. In particular, software packages for the automated solutions of partial differential equations using finite element discretization, such as FEniCS, are undergoing active support and development. However, because diffusion MRI simulation is a relatively new application area, there is still a gap between the simulation needs of the MRI community and the available tools provided by finite element software packages. In this paper, we address two potential difficulties in using FEniCS for diffusion MRI simulation. First, we simplified software installation by the use of FEniCS containers that are completely portable across multiple platforms. Second, we provide a portable simulation framework based on Python and whose code is open source. This simulation framework can be seamlessly integrated with cloud computing resources such as Google Colaboratory notebooks working on a web browser or with Google Cloud Platform with MPI parallelization. We show examples illustrating the accuracy, the computational times, and parallel computing capabilities. The framework contributes to reproducible science and open-source software in computational diffusion MRI with the hope that it will help to speed up method developments and stimulate research collaborations.

Quantitative DLA-based compressed sensing for T1-weighted acquisitions.

High resolution Manganese Enhanced Magnetic Resonance Imaging (MEMRI), which uses manganese as a T1 contrast agent, has great potential for functional imaging of live neuronal tissue at single neuron scale. However, reaching high resolutions often requires long acquisition times which can lead to reduced image quality due to sample deterioration and hardware instability. Compressed Sensing (CS) techniques offer the opportunity to significantly reduce the imaging time. The purpose of this work is to test the feasibility of CS acquisitions based on Diffusion Limited Aggregation (DLA) sampling patterns for high resolution quantitative T1-weighted imaging. Fully encoded and DLA-CS T1-weighted images of Aplysia californica neural tissue were acquired on a 17.2T MRI system. The MR signal corresponding to single, identified neurons was quantified for both versions of the T1 weighted images. For a 50% undersampling, DLA-CS can accurately quantify signal intensities in T1-weighted acquisitions leading to only 1.37% differences when compared to the fully encoded data, with minimal impact on image spatial resolution. In addition, we compared the conventional polynomial undersampling scheme with the DLA and showed that, for the data at hand, the latter performs better. Depending on the image signal to noise ratio, higher undersampling ratios can be used to further reduce the acquisition time in MEMRI based functional studies of living tissues.

A two-pool model to describe the IVIM cerebral perfusion.

IntraVoxel Incoherent Motion (IVIM) is a magnetic resonance imaging (MRI) technique capable of measuring perfusion-related parameters. In this manuscript, we show that the mono-exponential model commonly used to process IVIM data might be challenged, especially at short diffusion times. Eleven rat datasets were acquired at 7T using a diffusion-weighted pulsed gradient spin echo sequence with b-values ranging from 7 to 2500 s/mm2 at three diffusion times. The IVIM signals, obtained by removing the diffusion component from the raw MR signal, were fitted to the standard mono-exponential model, a bi-exponential model and the Kennan model. The Akaike information criterion used to find the best model to fit the data demonstrates that, at short diffusion times, the bi-exponential IVIM model is most appropriate. The results obtained by comparing the experimental data to a dictionary of numerical simulations of the IVIM signal in microvascular networks support the hypothesis that such a bi-exponential behavior can be explained by considering the contribution of two vascular pools: capillaries and somewhat larger vessels.

Numerical study of a cylinder model of the diffusion MRI signal for neuronal dendrite trees.

We study numerically how the neuronal dendrite tree structure can affect the diffusion magnetic resonance imaging (dMRI) signal in brain tissue. For a large set of randomly generated dendrite trees, synthetic dMRI signals are computed and fitted to a cylinder model to estimate the effective longitudinal diffusivity D(L) in the direction of neurites. When the dendrite branches are short compared to the diffusion length, D(L) depends significantly on the ratio between the average branch length and the diffusion length. In turn, D(L) has very weak dependence on the distribution of branch lengths and orientations of a dendrite tree, and the number of branches per node. We conclude that the cylinder model which ignores the connectivity of the dendrite tree, can still be adapted to describe the apparent diffusion coefficient in brain tissue.

DLA based compressed sensing for high resolution MR microscopy of neuronal tissue.

In this work we present the implementation of compressed sensing (CS) on a high field preclinical scanner (17.2 T) using an undersampling trajectory based on the diffusion limited aggregation (DLA) random growth model. When applied to a library of images this approach performs better than the traditional undersampling based on the polynomial probability density function. In addition, we show that the method is applicable to imaging live neuronal tissues, allowing significantly shorter acquisition times while maintaining the image quality necessary for identifying the majority of neurons via an automatic cell segmentation algorithm.

Parameter estimation using macroscopic diffusion MRI signal models.

Macroscopic models of the diffusion MRI (dMRI) signal can be helpful to understanding the relationship between the tissue microstructure and the dMRI signal. We study the least squares problem associated with estimating tissue parameters such as the cellular volume fraction, the residence times and the effective diffusion coefficients using a recently developed macroscopic model of the dMRI signal called the Finite Pulse Kärger model that generalizes the original Kärger model to non-narrow gradient pulses. In order to analyze the quality of the estimation in a controlled way, we generated synthetic noisy dMRI signals by including the effect of noise on the exact signal produced by the Finite Pulse Kärger model. The noisy signals were then fitted using the macroscopic model. Minimizing the least squares, we estimated the model parameters. The bias and standard deviations of the estimated model parameters as a function of the signal to noise ratio (SNR) were obtained. We discuss the choice of the b-values, the least square weights, the extension to experimentally obtained dMRI data as well noise correction.

Exploring diffusion across permeable barriers at high gradients. I. Narrow pulse approximation.

The adaptive variation of the gradient intensity with the diffusion time at a constant optimal b-value is proposed to enhance the contribution of the nuclei diffusing across permeable barriers, to the pulsed-gradient spin-echo (PGSE) signal. An exact simple formula the PGSE signal is derived under the narrow pulse approximation in the case of one-dimensional diffusion across a single permeable barrier. The barrier contribution to the signal is shown to be maximal at a particular b-value. The exact formula is then extended to multiple permeable barriers, while the PGSE signal is shown to be sensitive to the permeability and to the inter-barrier distance. Potential applications of the protocol to survey diffusion in three-dimensional domains with permeable membranes are illustrated through numerical simulations.

Numerical simulation of diffusion MRI signals using an adaptive time-stepping method.

The effect on the MRI signal of water diffusion in biological tissues in the presence of applied magnetic field gradient pulses can be modelled by a multiple compartment Bloch-Torrey partial differential equation. We present a method for the numerical solution of this equation by coupling a standard Cartesian spatial discretization with an adaptive time discretization. The time discretization is done using the explicit Runge-Kutta-Chebyshev method, which is more efficient than the forward Euler time discretization for diffusive-type problems. We use this approach to simulate the diffusion MRI signal from the extra-cylindrical compartment in a tissue model of the brain gray matter consisting of cylindrical and spherical cells and illustrate the effect of cell membrane permeability.

Numerical study of a macroscopic finite pulse model of the diffusion MRI signal.

Diffusion magnetic resonance imaging (dMRI) is an imaging modality that probes the diffusion characteristics of a sample via the application of magnetic field gradient pulses. The dMRI signal from a heterogeneous sample includes the contribution of the water proton magnetization from all spatial positions in a voxel. If the voxel can be spatially divided into different Gaussian diffusion compartments with inter-compartment exchange governed by linear kinetics, then the dMRI signal can be approximated using the macroscopic Karger model, which is a system of coupled ordinary differential equations (ODEs), under the assumption that the duration of the diffusion-encoding gradient pulses is short compared to the diffusion time (the narrow pulse assumption). Recently, a new macroscopic model of the dMRI signal, without the narrow pulse restriction, was derived from the Bloch-Torrey partial differential equation (PDE) using periodic homogenization techniques. When restricted to narrow pulses, this new homogenized model has the same form as the Karger model. We conduct a numerical study of the new homogenized model for voxels that are made up of periodic copies of a representative volume that contains spherical and cylindrical cells of various sizes and orientations and show that the signal predicted by the new model approaches the reference signal obtained by solving the full Bloch-Torrey PDE in O(ε(2)), where ε is the ratio between the size of the representative volume and a measure of the diffusion length. When the narrow gradient pulse assumption is not satisfied, the new homogenized model offers a much better approximation of the full PDE signal than the Karger model. Finally, preliminary results of applying the new model to a voxel that is not made up of periodic copies of a representative volume are shown and discussed.

Characterization of glioma microcirculation and tissue features using intravoxel incoherent motion magnetic resonance imaging in a rat brain model.

PURPOSE: Our aim was to investigate the pertinence of diffusion and perfusion magnetic resonance imaging (MRI) parameters obtained at 17.2 T in a 9L glioma rat brain tumor model to evaluate tumor tissue characteristics. MATERIALS AND METHODS: The local animal ethics advisory committee approved this study. 9L glioma cells were injected intracerebrally to 14 Fischer rats. The animals were imaged at 7 or 12 days after implantation on a 17.2-T MRI scanner, using 72 different b values (2-3025 s/mm(2)). The signal attenuation, S/So, was fitted using a kurtosis diffusion model (ADCo and K) and a biexponential diffusion model (fractions ffast and fslow and diffusion coefficients Dfast and Dslow) using b values greater than 300 s/mm(2). To bridge the 2 models, an average diffusion coefficient and a biexponential index were estimated from the biexponential model as ADCo and K equivalents, respectively. Intravoxel incoherent motion perfusion-related parameters were obtained from the residual signal at low b values, after the diffusion component has been removed. Diffusion and perfusion maps were generated for each fitted parameter on a pixel-by-pixel basis, and regions of interest were drawn in the tumor and contralateral side to retrieve diffusion and perfusion parameters. All rats were killed and cellularity and vascularity were quantitatively assessed using histology for comparison with diffusion and perfusion parameters. RESULTS: Intravoxel incoherent motion maps clearly highlighted tumor areas as generally heterogeneous, as confirmed by histology. For diffusion parameters, ADCo and were not significantly different between the tumor and contralateral side, whereas K in the tumor was significantly higher than in contralateral basal ganglia (P < 0.0001), as well as biexponential index (P < 0.001). ADCo and in the tumor at day 7 were significantly higher than at day 12 (P < 0.01 and P < 0.001, respectively). fIVIM in the tumor from the kurtosis diffusion model was significantly higher than in contralateral basal ganglia (P < 0.001). fIVIM in the tumor at day 7 was significantly higher than in the tumor at day 12 (P < 0.0001). There was no significant difference for D* between the tumor and contralateral side (P = 0.06). A significant negative correlation was found between tumor vascularity and fIVIM (P < 0.05) as well as between tumor cell count and (P < 0.01). CONCLUSION: Quantitative non-Gaussian diffusion and perfusion MRI can provide valuable information on microvasculature and tissue structure to improve characterization of brain tumors.