Representing diffusion MRI in 5D for segmentation of white matter tracts with a level set method.

We present a method for segmenting white matter tracts from high angular resolution diffusion MR. images by representing the data in a 5 dimensional space of position and orientation. Whereas crossing fiber tracts cannot be separated in 3D position space, they clearly disentangle in 5D position-orientation space. The segmentation is done using a 5D level set method applied to hyper-surfaces evolving in 5D position-orientation space. In this paper we present a methodology for constructing the position-orientation space. We then show how to implement the standard level set method in such a non-Euclidean high dimensional space. The level set theory is basically defined for N-dimensions but there are several practical implementation details to consider, such as mean curvature. Finally, we will show results from a synthetic model and a few preliminary results on real data of a human brain acquired by high angular resolution diffusion MRI.

Diffusion spectrum magnetic resonance imaging (DSI) tractography of crossing fibers.

MRI tractography is the mapping of neural fiber pathways based on diffusion MRI of tissue diffusion anisotropy. Tractography based on diffusion tensor imaging (DTI) cannot directly image multiple fiber orientations within a single voxel. To address this limitation, diffusion spectrum MRI (DSI) and related methods were developed to image complex distributions of intravoxel fiber orientation. Here we demonstrate that tractography based on DSI has the capacity to image crossing fibers in neural tissue. DSI was performed in formalin-fixed brains of adult macaque and in the brains of healthy human subjects. Fiber tract solutions were constructed by a streamline procedure, following directions of maximum diffusion at every point, and analyzed in an interactive visualization environment (TrackVis). We report that DSI tractography accurately shows the known anatomic fiber crossings in optic chiasm, centrum semiovale, and brainstem; fiber intersections in gray matter, including cerebellar folia and the caudate nucleus; and radial fiber architecture in cerebral cortex. In contrast, none of these examples of fiber crossing and complex structure was identified by DTI analysis of the same data sets. These findings indicate that DSI tractography is able to image crossing fibers in neural tissue, an essential step toward non-invasive imaging of connectional neuroanatomy.

Diffusion-weighted MRI in metastatic gastrointestinal tumors (GIST): A pilot study on the assessment of treatment response in comparison with 18F-FDG PET/CT

Purpose: To evaluate the clinical potential of diffusion-weighted MR imaging with apparent diffusion coefficient (ADC) mapping for the assessment of gastrointestinal stromal tumor (GIST) response to targeted therapy in comparison with 18F-FDG PET/CT. Methods and materials: Five patients (3W/2M, aged 56 ± 13 y) with metastatic GIST underwent both a 18F-FDG PET/CT (Discovery LS, GE Healthcare) and a MRI (VIBE T1 Gd, DWI [b = 50,300,600] and ADC mapping) before and after change in therapy. Exams were first analyzed blindly, then PET/CT images were coregistered to T1 Gd MR images for lesion detection. SUVmax and ADC were measured for the six largest lesions on MRI. The relationship between SUVmax and ADC was analyzed using Spearman's correlation. Results: Altogether, 24 lesions (15 hepatic and 9 non-hepatic) were analyzed on both modalities. Three PET/CT lesions (12.5%) were initially not considered on ADC and 4 lesions on the second PET/CT were excluded because of hepatic vascular activity spillover. SUVmax decreased from 7.2 ± 7.7 g/mL to 5.9 ± 5.9 g/mL (P = 0.53) and ADC increased from 1.2x10-3 mm2/s ± 0.4 to 1.4x10-3 mm2/s ± 0.4 (P = 0.07). There was a significant association between SUVmax decrease and ADC increase (rho= -0.64, P = 0.004). Conclusion: Changes in ADC from diffusion-weighted MRI reflect response of 18F-FDG-avid GIST to therapy. The exact diagnostic value of DWI needs to be investigated further, as well as the effect of lesion size and time under therapy before imaging. Furthermore, the proven association between SUVmax and ADC may be useful for the assessment of treatment response in 18F-FDG non-avid GIST.

Mutual information of IID complex Gaussian signals on block Rayleigh-faded channels

We present a method to compute, quickly and efficiently, the mutual information achieved by an IID (independent identically distributed) complex Gaussian signal on a block Rayleigh-faded channel without side information at the receiver. The method accommodates both scalar and MIMO (multiple-input multiple-output) settings. Operationally, this mutual information represents the highest spectral efficiency that can be attained using Gaussiancodebooks. Examples are provided that illustrate the loss in spectral efficiency caused by fast fading and how that loss is amplified when multiple transmit antennas are used. These examples are further enriched by comparisons with the channel capacity under perfect channel-state information at the receiver, and with the spectral efficiency attained by pilot-based transmission.

Mutual Information of IID complex Gaussian signals on block Rayleigh-faded channels

We present a method to compute, quickly and efficiently, the mutual information achieved by an IID (independent identically distributed) complex Gaussian signal on a block Rayleigh-faded channel without side information at the receiver. The method accommodates both scalar and MIMO (multiple-input multiple-output) settings. Operationally, this mutual information represents the highest spectral efficiency that can be attained using Gaussiancodebooks. Examples are provided that illustrate the loss in spectral efficiency caused by fast fading and how that loss is amplified when multiple transmit antennas are used. These examples are further enriched by comparisons with the channel capacity under perfect channel-state information at the receiver, and with the spectral efficiency attained by pilot-based transmission.

Diffusion-weighted MRI in metastatic gastrointestinal tumours (GIST): a pilot study on the assessment of treatment response in comparison with 18 F-FDG PET/CT

Purpose: To evaluate the clinical potential of diffusion-weighted MR imaging with apparent diffusion coefficient (ADC) mapping for the assessment of gastrointestinal stromal tumour (GIST) response to targeted therapy in comparison with 18F-FDG PET/CT Methods and Materials: Five patients (3 W/2M, aged 56±13 y) with metastatic GIST underwent both a 18F-FDG PET/CT (Discovery LS, GE Healthcare) and a MRI (VIBE T1 Gd, DWI [b = 50,300,600] and ADC mapping) before and after change in therapy. Exams were first analysed blindly and then PET/CT images were coregistered to T1 Gd MR images for lesion detection. SUVmax and ADC were measured for the six largest lesions on MRI. The relationship between SUVmax and ADC was analysed using Spearman's correlation. Results: Altogether, 24 lesions (15 hepatic and 9 non-hepatic) were analysed on both modalities. Three PET/CT lesions (12.5%) were initially not considered on ADC and 4 lesions on the second PET/CT were excluded because of hepatic vascular activity spillover. SUVmax decreased from 7.2±7.7 g/mL to 5.9±5.9 g/mL (P = 0.53) and ADC increased from 1.2x10-3 mm2/s ± 0.4 to 1.4x10-3 mm2/s ± 0.4 (P = 0.07). There was a significant association between SUVmax decrease and ADC increase (rho= -0.64, P = 0.004). Conclusion: Changes in ADC from diffusion-weighted MRI reflect response of 18F-FDG-avid GIST to therapy. The exact diagnostic value of DWI needs to be investigated further, as well as the effect of lesion size and time under therapy before imaging. Furthermore, the proven association between SUVmax and ADC may be useful for the assessment of treatment response in 18F-FDG non-avid GIST.

Functional MRI evaluation of liver tumour response after radiofrequency: short- and mid-term evolution of diffusion parameters

Purpose: To evaluate the short- and mid-term evolutions of the apparent diffusion coefficient of lesions treated with RF, in order to determine if the ADC can be used as a marker of tumour response. Methods and Materials: Twenty patients were treated for a liver malignancy with RF and were examined on a 1.5 T/3.0 T machine with T2, gadolinium-enhanced T1 and diffusion sequences: before treatment (< 1 month), just after treatment (< 1 month) and midterm (3-6 months). The ADC was measured in the whole lesion and in the area with the most restricted diffusion (MRDA). The ROI size was also measured on the diffusion map. The Pearson/ANOVA tests were used. Results: All patients were successfully treated with complete disappearance of CE. The lesional size on T2 showed a negative evolution in time (p < 0.002). The ADC in the whole lesion showed a bell-shaped evolution (increasing just after RF, then decreasing, p = 0.02). The ROI size on the diffusion map followed a similar course (p = 0.01). For the MRDA, such evolutions were also found, but they were not significant. There was a negative correlation between CE and the ADC (p < 0.02) and between the lesional size on T2 and ADC (p = 0.03) in the whole lesion. There were also positive correlations between the ROI size and ADC (p = 0.0008) and between CE and the size on T2 (p = 0.0002). The ADC in MRDA showed some non-significant correlations with other variables. Conclusion: The lesions successfully treated with RF have a clear and predictable evolution in terms of T2 size, CE and ADC.

Accelerated tumor invasion under non-isotropic cell dispersal in glioblastomas

Glioblastomas are highly diffuse, malignant tumors that have so far evaded clinical treatment. The strongly invasive behavior of cells in these tumors makes them very resistant to treatment, and for this reason both experimental and theoretical efforts have been directed toward understanding the spatiotemporal pattern of tumor spreading. Although usual models assume a standard diffusion behavior, recent experiments with cell cultures indicate that cells tend to move in directions close to that of glioblastoma invasion, thus indicating that a biasedrandom walk model may be much more appropriate. Here we show analytically that, for realistic parameter values, the speeds predicted by biased dispersal are consistent with experimentally measured data. We also find that models beyond reaction–diffusion–advection equations are necessary to capture this substantial effect of biased dispersal on glioblastoma spread

A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion

An epidemic model is formulated by a reactionâeuro"diffusion system where the spatial pattern formation is driven by cross-diffusion. The reaction terms describe the local dynamics of susceptible and infected species, whereas the diffusion terms account for the spatial distribution dynamics. For both self-diffusion and cross-diffusion, nonlinear constitutive assumptions are suggested. To simulate the pattern formation two finite volume formulations are proposed, which employ a conservative and a non-conservative discretization, respectively. An efficient simulation is obtained by a fully adaptive multiresolution strategy. Numerical examples illustrate the impact of the cross-diffusion on the pattern formation.

l0-deconvolution for compressive diffusion MRI

Diffusion MRI is a well established imaging modality providing a powerful way to non-invasively probe the structure of the white matter. Despite the potential of the technique, the intrinsic long scan times of these sequences have hampered their use in clinical practice. For this reason, a wide variety of methods have been proposed to shorten acquisition times. [...] We here review a recent work where we propose to further exploit the versatility of compressed sensing and convex optimization with the aim to characterize the fiber orientation distribution sparsity more optimally. We re-formulate the spherical deconvolution problem as a constrained l0 minimization.

Diffusion tensor imaging shows structural remodeling of stroke mirror region: results from a pilot study.

Background: The role of the non-injured hemisphere in stroke recovery is poorly understood. In this pilot study, we sought to explore the presence of structural changes detectable by diffusion tensor imaging (DTI) in the contralesional hemispheres of patients who recovered well from ischemic stroke. Methods: We analyzed serial DTI data from 16 stroke patients who had moderate initial neurological deficits (NIHSS scores 3-12) and good functional outcome at 3-6 months (NIHSS score 0 or modified Rankin Score ≤1). We segmented the brain tissue in gray and white matter (GM and WM) and measured the apparent diffusion coefficient (ADC) and fractional anisotropy in the infarct, in the contralesional infarct mirror region as well as in concentrically expanding regions around them. Results: We found that GM and WM ADC significantly increased in the infarct region (p < 0.01) from acute to chronic time points, whereas in the infarct mirror region, GM and WM ADC increased (p < 0.01) and WM fractional anisotropy decreased (p < 0.05). No significant changes were detected in other regions. Conclusion: DTI-based metrics are sensitive to regional structural changes in the contralesional hemisphere during stroke recovery. Prospective studies in larger cohorts with varying levels of recovery are needed to confirm our findings.

Kinematic reduction of reaction-diffusion fronts with multiplicative noise. Derivation of stochastic sharp-interface equations

We study the dynamics of generic reaction-diffusion fronts, including pulses and chemical waves, in the presence of multiplicative noise. We discuss the connection between the reaction-diffusion Langevin-like field equations and the kinematic (eikonal) description in terms of a stochastic moving-boundary or sharp-interface approximation. We find that the effective noise is additive and we relate its strength to the noise parameters in the original field equations, to first order in noise strength, but including a partial resummation to all orders which captures the singular dependence on the microscopic cutoff associated with the spatial correlation of the noise. This dependence is essential for a quantitative and qualitative understanding of fluctuating fronts, affecting both scaling properties and nonuniversal quantities. Our results predict phenomena such as the shift of the transition point between the pushed and pulled regimes of front propagation, in terms of the noise parameters, and the corresponding transition to a non-Kardar-Parisi-Zhang universality class. We assess the quantitative validity of the results in several examples including equilibrium fluctuations and kinetic roughening. We also predict and observe a noise-induced pushed-pulled transition. The analytical predictions are successfully tested against rigorous results and show excellent agreement with numerical simulations of reaction-diffusion field equations with multiplicative noise.

Mean first-passage time of continuous non-Markovian processes driven by colored noise

An equation for mean first-passage times of non-Markovian processes driven by colored noise is derived through an appropriate backward integro-differential equation. The equation is solved in a Bourret-like approximation. In a weak-noise bistable situation, non-Markovian effects are taken into account by an effective diffusion coefficient. In this situation, our results compare satisfactorily with other approaches and experimental data.

RNA interference against aquaporin-4 decreases the apparent diffusion coefficient in the normal brain

Aim: Diffusion weighted magnetic resonance imaging (MRI) is now widely used in human brain diagnosis.1 To date molecular mechanisms underlying changes in Apparent Diffusion Coefficient (ADC) signals remain poorly understood. AQP4, localized to astrocytes, is one of the most highly expressed cerebral AQPs.2 AQP4 is involved in water movement within the cell membrane of cultured astrocytes.3 We hypothesize that AQP4 contributes to water diffusion and underlying ADC values in normal brain. Methods: We used an RNA interference (RNAi) protocol in vivo, to acutely knockdown expression of AQP4 in rat brain and to determine whether this was associated with changes in brain ADC values using MRI protocols as previously described.4 RNAi was performed using specific small interference RNA (siRNA) against AQP4 (siAQP4) and a non-targeted-siRNA (siGLO) as a control. The specificity and efficiency of the siAQP4 were first tested in vitro in astrocyte and hippocampal slice cultures. In vivo, siRNAs were injected into the rat cortex 3d prior to MRI acquisition and AQP4 was assessed by western blot (n=4) and immunohistochemistry (n=6). Histology was performed on adjacent slices. Results: siAQP4 application on primary astrocyte cultures induced a 76% decrease in AQP4 expression after 4 days. In hippocampal slice cultures; we also found a significant decrease in AQP4 expression in astrocytes after siAQP4. In vivo, injection of non-targeted siRNA (siGLO) tagged with CY3 allowed us to show that GFAP positive cells (astrocytes) were positively stained with CY3-siGLO, showing efficient transfection. Western blot and immunohistochemical analysis showed that siAQP4 induced a ~30% decrease in AQP4 expression without modification of tissue properties or cell death. After siAQP4 treatment, a significant decrease in ADC values (~50%) were observed without altered of T2 values. Conclusions: Together these results suggest that AQP4 reduces water diffusion through the astrocytic plasma membrane and decreases ADC values. Our findings demonstrate for the first time that astrocytic AQP4 contributes significantly to brain water diffusion and ADC values in normal brain. These results open new avenues to interpretation of ADC values under normal physiological conditions and in acute and chronic brain injuries.

Accelerated Microstructure Imaging via Convex Optimization (AMICO) from diffusion MRI data.

Microstructure imaging from diffusion magnetic resonance (MR) data represents an invaluable tool to study non-invasively the morphology of tissues and to provide a biological insight into their microstructural organization. In recent years, a variety of biophysical models have been proposed to associate particular patterns observed in the measured signal with specific microstructural properties of the neuronal tissue, such as axon diameter and fiber density. Despite very appealing results showing that the estimated microstructure indices agree very well with histological examinations, existing techniques require computationally very expensive non-linear procedures to fit the models to the data which, in practice, demand the use of powerful computer clusters for large-scale applications. In this work, we present a general framework for Accelerated Microstructure Imaging via Convex Optimization (AMICO) and show how to re-formulate this class of techniques as convenient linear systems which, then, can be efficiently solved using very fast algorithms. We demonstrate this linearization of the fitting problem for two specific models, i.e. ActiveAx and NODDI, providing a very attractive alternative for parameter estimation in those techniques; however, the AMICO framework is general and flexible enough to work also for the wider space of microstructure imaging methods. Results demonstrate that AMICO represents an effective means to accelerate the fit of existing techniques drastically (up to four orders of magnitude faster) while preserving accuracy and precision in the estimated model parameters (correlation above 0.9). We believe that the availability of such ultrafast algorithms will help to accelerate the spread of microstructure imaging to larger cohorts of patients and to study a wider spectrum of neurological disorders.