Modeling brain dynamics in brain tumor patients using the virtual brain

Presurgical planning for brain tumor resection aims at delineating eloquent tissue in the vicinity of the lesion to spare during surgery. To this end, noninvasive neuroimaging techniques such as functional MRI and diffusion-weighted imaging fiber tracking are currently employed. However, taking into account this information is often still insufficient, as the complex nonlinear dynamics of the brain impede straightforward prediction of functional outcome after surgical intervention. Large-scale brain network modeling carries the potential to bridge this gap by integrating neuroimaging data with biophysically based models to predict collective brain dynamics. As a first step in this direction, an appropriate computational model has to be selected, after which suitable model parameter values have to be determined. To this end, we simulated large-scale brain dynamics in 25 human brain tumor patients and 11 human control participants using The Virtual Brain, an open-source neuroinformatics platform. Local and global model parameters of the Reduced Wong-Wang model were individually optimized and compared between brain tumor patients and control subjects. In addition, the relationship between model parameters and structural network topology and cognitive performance was assessed. Results showed (1) significantly improved prediction accuracy of individual functional connectivity when using individually optimized model parameters; (2) local model parameters that can differentiate between regions directly affected by a tumor, regions distant from a tumor, and regions in a healthy brain; and (3) interesting associations between individually optimized model parameters and structural network topology and cognitive performance

Estimation of Fiber Orientations Using Neighborhood Information

Data from diffusion magnetic resonance imaging (dMRI) can be used to reconstruct fiber tracts, for example, in muscle and white matter. Estimation of fiber orientations (FOs) is a crucial step in the reconstruction process and these estimates can be corrupted by noise. In this paper, a new method called Fiber Orientation Reconstruction using Neighborhood Information (FORNI) is described and shown to reduce the effects of noise and improve FO estimation performance by incorporating spatial consistency. FORNI uses a fixed tensor basis to model the diffusion weighted signals, which has the advantage of providing an explicit relationship between the basis vectors and the FOs. FO spatial coherence is encouraged using weighted l1-norm regularization terms, which contain the interaction of directional information between neighbor voxels. Data fidelity is encouraged using a squared error between the observed and reconstructed diffusion weighted signals. After appropriate weighting of these competing objectives, the resulting objective function is minimized using a block coordinate descent algorithm, and a straightforward parallelization strategy is used to speed up processing. Experiments were performed on a digital crossing phantom, ex vivo tongue dMRI data, and in vivo brain dMRI data for both qualitative and quantitative evaluation. The results demonstrate that FORNI improves the quality of FO estimation over other state of the art algorithms.Comment: Journal paper accepted in Medical Image Analysis. 35 pages and 16 figure

Connectivity-enhanced diffusion analysis reveals white matter density disruptions in first episode and chronic schizophrenia.

Reduced fractional anisotropy (FA) is a well-established correlate of schizophrenia, but it remains unclear whether these tensor-based differences are the result of axon damage and/or organizational changes and whether the changes are progressive in the adult course of illness. Diffusion MRI data were collected in 81 schizophrenia patients (54 first episode and 27 chronic) and 64 controls. Analysis of FA was combined with "fixel-based" analysis, the latter of which leverages connectivity and crossing-fiber information to assess both fiber bundle density and organizational complexity (i.e., presence and magnitude of off-axis diffusion signal). Compared with controls, patients with schizophrenia displayed clusters of significantly lower FA in the bilateral frontal lobes, right dorsal centrum semiovale, and the left anterior limb of the internal capsule. All FA-based group differences overlapped substantially with regions containing complex fiber architecture. FA within these clusters was positively correlated with principal axis fiber density, but inversely correlated with both secondary/tertiary axis fiber density and voxel-wise fiber complexity. Crossing fiber complexity had the strongest (inverse) association with FA (r = -0.82). When crossing fiber structure was modeled in the MRtrix fixel-based analysis pipeline, patients exhibited significantly lower fiber density compared to controls in the dorsal and posterior corpus callosum (central, postcentral, and forceps major). Findings of lower FA in patients with schizophrenia likely reflect two inversely related signals: reduced density of principal axis fiber tracts and increased off-axis diffusion sources. Whereas the former confirms at least some regions where myelin and or/axon count are lower in schizophrenia, the latter indicates that the FA signal from principal axis fiber coherence is broadly contaminated by macrostructural complexity, and therefore does not necessarily reflect microstructural group differences. These results underline the need to move beyond tensor-based models in favor of acquisition and analysis techniques that can help disambiguate different sources of white matter disruptions associated with schizophrenia

Improving Fiber Alignment in HARDI by Combining Contextual PDE Flow with Constrained Spherical Deconvolution

We propose two strategies to improve the quality of tractography results computed from diffusion weighted magnetic resonance imaging (DW-MRI) data. Both methods are based on the same PDE framework, defined in the coupled space of positions and orientations, associated with a stochastic process describing the enhancement of elongated structures while preserving crossing structures. In the first method we use the enhancement PDE for contextual regularization of a fiber orientation distribution (FOD) that is obtained on individual voxels from high angular resolution diffusion imaging (HARDI) data via constrained spherical deconvolution (CSD). Thereby we improve the FOD as input for subsequent tractography. Secondly, we introduce the fiber to bundle coherence (FBC), a measure for quantification of fiber alignment. The FBC is computed from a tractography result using the same PDE framework and provides a criterion for removing the spurious fibers. We validate the proposed combination of CSD and enhancement on phantom data and on human data, acquired with different scanning protocols. On the phantom data we find that PDE enhancements improve both local metrics and global metrics of tractography results, compared to CSD without enhancements. On the human data we show that the enhancements allow for a better reconstruction of crossing fiber bundles and they reduce the variability of the tractography output with respect to the acquisition parameters. Finally, we show that both the enhancement of the FODs and the use of the FBC measure on the tractography improve the stability with respect to different stochastic realizations of probabilistic tractography. This is shown in a clinical application: the reconstruction of the optic radiation for epilepsy surgery planning

Gaussian process regression can turn non-uniform and undersampled diffusion MRI data into diffusion spectrum imaging

We propose to use Gaussian process regression to accurately estimate the diffusion MRI signal at arbitrary locations in q-space. By estimating the signal on a grid, we can do synthetic diffusion spectrum imaging: reconstructing the ensemble averaged propagator (EAP) by an inverse Fourier transform. We also propose an alternative reconstruction method guaranteeing a nonnegative EAP that integrates to unity. The reconstruction is validated on data simulated from two Gaussians at various crossing angles. Moreover, we demonstrate on non-uniformly sampled in vivo data that the method is far superior to linear interpolation, and allows a drastic undersampling of the data with only a minor loss of accuracy. We envision the method as a potential replacement for standard diffusion spectrum imaging, in particular when acquistion time is limited.Comment: 5 page