Variational, geometric, and statistical methods for modeling brain anatomy and function

Developments in the human machine interface technologies and their applications: a review

Human-machine interface (HMI) techniques use bioelectrical signals to gain real-time synchronised communication between the human body and machine functioning. HMI technology not only provides a real-time control access but also has the ability to control multiple functions at a single instance of time with modest human inputs and increased efficiency. The HMI technologies yield advanced control access on numerous applications such as health monitoring, medical diagnostics, development of prosthetic and assistive devices, automotive and aerospace industry, robotic controls and many more fields. In this paper, various physiological signals, their acquisition and processing techniques along with their respective applications in different HMI technologies have been discussed.

Bayesian Tractography Using Geometric Shape Priors

The problem of estimating neuronal fiber tracts connecting different brain regions is important for various types of brain studies, including understanding brain functionality and diagnosing cognitive impairments. The popular techniques for tractography are mostly sequential—tracts are grown sequentially following principal directions of local water diffusion profiles. Despite several advancements on this basic idea, the solutions easily get stuck in local solutions, and can't incorporate global shape information. We present a global approach where fiber tracts between regions of interest are initialized and updated via deformations based on gradients of a posterior energy. This energy has contributions from diffusion data, global shape models, and roughness penalty. The resulting tracts are relatively immune to issues such as tensor noise and fiber crossings, and achieve more interpretable tractography results. We demonstrate this framework using both simulated and real dMRI and HARDI data.

Awake surgery between art and science. Part I: clinical and operative settings

Awake surgery requires coordinated teamwork and communication between the surgeon and the anesthesiologist, as he monitors the patient, the neuroradiologist as he interprets the images for intraoperative confirmation, and the neuropsychologist and neurophysiologist as they evaluate in real-time the patient's responses to commands and questions. To improve comparison across published studies on clinical assessment and operative settings in awake surgery, we reviewed the literature, focusing on methodological differences and aims. In complex, interdisciplinary medical care, such differences can affect the outcome and the cost-benefit ratio of the treatment. Standardization of intraoperative mapping and related controversies will be discussed in Part II.

A robust variational approach for simultaneous smoothing and estimation of DTI

Estimating diffusion tensors is an essential step in many applications - such as diffusion tensor image (DTI) registration, segmentation and fiber tractography. Most of the methods proposed in the literature for this task are not simultaneously statistically robust and feature preserving techniques. In this paper, we propose a novel and robust variational framework for simultaneous smoothing and estimation of diffusion tensors from diffusion MRI. Our variational principle makes use of a recently introduced total Kullback-Leibler (tKL) divergence for DTI regularization. tKL is a statistically robust dissimilarity measure for diffusion tensors, and regularization by using tKL ensures the symmetric positive definiteness of tensors automatically. Further, the regularization is weighted by a non-local factor adapted from the conventional non-local means filters. Finally, for the data fidelity, we use the nonlinear least-squares term derived from the Stejskal-Tanner model. We present experimental results depicting the positive performance of our method in comparison to competing methods on synthetic and real data examples.

Multiple cortical surface correspondence using pairwise shape similarity

Accurately corresponding a population of human cortical surfaces provides important shape information for the diagnosis of many brain diseases. This problem is very challenging due to the highly convoluted nature of cortical surfaces. Pairwise methods using a fixed template may not handle well the case when a target cortical surface is substantially different from the template. In this paper, we develop a new method to organize the population of cortical surfaces into pairs with high shape similarity and only correspond such similar pairs to achieve a higher accuracy. In particular, we use the geometric information to identify co-located gyri and sulci for defining a new measure of shape similarity. We conduct experiments on 40 instances of the cortical surface, resulting in an improved performance over several existing shape-correspondence methods.

A future for systems and computational neuroscience in France?

This special issue of the Journal of Physiology, Paris, is an outcome of NeuroComp'06, the first French conference in Computational Neuroscience. The preparation for this conference, held at Pont-à-Mousson in October 2006, was accompanied by a survey which has resulted in an up-to-date inventory of human resources and labs in France concerned with this emerging new field of research (see team directory in http://neurocomp.risc.cnrs.fr/). This thematic JPP issue gathers some of the key scientific presentations made on the occasion of this first interdisciplinary meeting, which should soon become recognized as a yearly national conference representative of a new scientific community. The present introductory paper presents the general scientific context of the conference and reviews some of the historical and conceptual foundations of Systems and Computational Neuroscience in France.

Diffusion-based spatial priors for imaging

We describe a Bayesian scheme to analyze images, which uses spatial priors encoded by a diffusion kernel, based on a weighted graph Laplacian. This provides a general framework to formulate a spatial model, whose parameters can be optimized. The application we have in mind is a spatiotemporal model for imaging data. We illustrate the method on a random effects analysis of fMRI contrast images from multiple subjects; this simplifies exposition of the model and enables a clear description of its salient features. Typically, imaging data are smoothed using a fixed Gaussian kernel as a pre-processing step before applying a mass-univariate statistical model (e.g., a general linear model) to provide images of parameter estimates. An alternative is to include smoothness in a multivariate statistical model (Penny, W.D., Trujillo-Barreto, N.J., Friston, K.J., 2005. Bayesian fMRI time series analysis with spatial priors. Neuroimage 24, 350–362). The advantage of the latter is that each parameter field is smoothed automatically, according to a measure of uncertainty, given the data. In this work, we investigate the use of diffusion kernels to encode spatial correlations among parameter estimates. Nonlinear diffusion has a long history in image processing; in particular, flows that depend on local image geometry (Romeny, B.M.T., 1994. Geometry-driven Diffusion in Computer Vision. Kluwer Academic Publishers) can be used as adaptive filters. This can furnish a non-stationary smoothing process that preserves features, which would otherwise be lost with a fixed Gaussian kernel. We describe a Bayesian framework that incorporates non-stationary, adaptive smoothing into a generative model to extract spatial features in parameter estimates. Critically, this means adaptive smoothing becomes an integral part of estimation and inference. We illustrate the method using synthetic and real fMRI data.

Parameter estimate of signal transduction pathways

BACKGROUND

The "inverse" problem is related to the determination of unknown causes on the bases of the observation of their effects. This is the opposite of the corresponding "direct" problem, which relates to the prediction of the effects generated by a complete description of some agencies. The solution of an inverse problem entails the construction of a mathematical model and takes the moves from a number of experimental data. In this respect, inverse problems are often ill-conditioned as the amount of experimental conditions available are often insufficient to unambiguously solve the mathematical model. Several approaches to solving inverse problems are possible, both computational and experimental, some of which are mentioned in this article. In this work, we will describe in details the attempt to solve an inverse problem which arose in the study of an intracellular signaling pathway.

RESULTS

Using the Genetic Algorithm to find the sub-optimal solution to the optimization problem, we have estimated a set of unknown parameters describing a kinetic model of a signaling pathway in the neuronal cell. The model is composed of mass action ordinary differential equations, where the kinetic parameters describe protein-protein interactions, protein synthesis and degradation. The algorithm has been implemented on a parallel platform. Several potential solutions of the problem have been computed, each solution being a set of model parameters. A sub-set of parameters has been selected on the basis on their small coefficient of variation across the ensemble of solutions.

CONCLUSION

Despite the lack of sufficiently reliable and homogeneous experimental data, the genetic algorithm approach has allowed to estimate the approximate value of a number of model parameters in a kinetic model of a signaling pathway: these parameters have been assessed to be relevant for the reproduction of the available experimental data.

Evaluation of the GTRACT diffusion tensor tractography algorithm: A validation and reliability study

Fiber tracking, based on diffusion tensor imaging (DTI), is the only approach available to non-invasively study the three-dimensional structure of white matter tracts. Two major obstacles to this technique are partial volume artifacts and tracking errors caused by image noise. In this paper, a novel fiber tracking algorithm called Guided Tensor Restore Anatomical Connectivity Tractography (GTRACT) is presented. This algorithm utilizes a multi-pass approach to fiber tracking. In the first pass, a 3D graph search algorithm is utilized. The second pass incorporates anatomical connectivity information generated in the first pass to guide the tracking in this stage. This approach improves the ability to reconstruct complex fiber paths as well as the tracking accuracy. Validation and reliability studies using this algorithm were performed on both synthetic phantom data and clinical human brain data. A method is also proposed for the evaluating reliability of fiber tract generation based both on the position of the fiber tracts, as well the anisotropy values along the path. The results demonstrate that the GTRACT algorithm is less sensitive to image noise and more capable of handling areas of complex fiber crossing, compared to conventional streamline methods.

Validation of semiautomated methods for quantifying cingulate cortical metrics in schizophrenia

This paper validates semiautomated methods for reconstructing cortical surfaces of the cingulate gyrus from high-resolution magnetic resonance (MR) images. Bayesian segmentation was used to delineate the image voxels into five tissue types: cerebrospinal fluid (CSF), gray matter (GM), white matter (WM), and partial volumes of CSF/GM and GM/WM; the tissues were then recalibrated as CSF, GM, and WM via the Neyman-Pearson Likelihood Ratio Test. To generate cortical surfaces at the interface of GM and WM, the thresholds between the tissue types were first used to reassign partial volume voxels to CSF, GM, and WM with minimum error (that varied from 0.06 to 0.15 for the 10 subjects). Next, topology-correct cortical surfaces were generated and validated with almost all surface vertices lying within one voxel (0.5 mm) of hand contours. Dynamic programming was used to delineate and extract the cingulate gyrus from the cortical surfaces based on its gyral and sulcal boundaries. The intraclass correlation coefficient for surface area obtained by two raters for all 10 surfaces was 0.82. In addition, by repeating the entire procedure three times in one subject, we obtained a coefficient of variation of 0.0438 for surface area.