Biomedical subjects

Habib Benali

Publications and source records attributed to Habib Benali.

At least 19 recordsLinked to original sources

PubMed2006

CORSICA: correction of structured noise in fMRI by automatic identification of ICA components.

When applied to functional magnetic resonance imaging (fMRI) data, spatial independent component analysis (sICA), a data-driven technique that addresses the blind source separation problem, seems able to extract components specifically related to physiological noise and brain movements. These components should be removed from the data to achieve structured noise reduction and improve any subsequent detection and analysis of signal fluctuations related to neural activity. We propose a new automatic method called CORSICA (CORrection of Structured noise using spatial Independent Component Analysis) to identify the components related to physiological noise, using prior information on the spatial localization of the main physiological fluctuations in fMRI data. As opposed to existing spectral priors, which may be subject to aliasing effects for long-TR data sets (typically acquired with TR >1 s), such spatial priors can be applied to fMRI data, regardless of the TR of the acquisitions. By comparing the proposed automatic selection to a manual selection performed visually by a human operator, we first show that CORSICA is able to identify the noise-related components for long-TR data with a high sensitivity and a specificity of 1. On short-TR data sets, we validate that the proposed method of noise reduction allows a substantial improvement of the signal-to-noise ratio evaluated at the cardiac and respiratory frequencies, even in the gray matter, while preserving the main fluctuations related to neural activity.

Biometry↗

PubMed2006

Preoperative estimation of residual volume for WHO grade II glioma resected with intraoperative functional mapping.

Despite the lack of class I evidence, it is widely agreed that surgery can improve the functional and vital prognosis for WHO grade II gliomas when the resection is at least subtotal radiologically, that is, leaving less than 10 cm(3) of visible residual tumor. Because these tumors frequently invade functional areas, the preoperative estimation of the probable residual volume remains challenging. This article presents a probabilistic map of postoperative residues, with the aim of predicting before the decision for surgical intervention whether the resection could be subtotal. We selected 65 patients who underwent surgery with intraoperative functional mapping between 1999 and 2004 for a WHO grade II glioma located in a sensorimotor and/or language area. For each case, the postoperative image was normalized on a standard atlas, and the residual tumor was segmented. A probabilistic map of residues was then computed. The fusion between the map and a preoperative image allowed a preoperative estimation of the expected extent of resection. The map enhances the regions where grade II glioma cannot be resected. The success rate for the preoperative classification of partial versus subtotal resection is 82%. Although both its reliability and accuracy have to be improved, this probabilistic map gives preoperatively an objective estimation of the expected extent of resection for grade II glioma resected under intraoperative functional mapping. This rationale will assist in decisions regarding surgical resection and may thus contribute to the elaboration of a therapeutic consensus for WHO grade II glioma.

Brain Mapping↗

PubMed2006

Partial correlation for functional brain interactivity investigation in functional MRI.

Examination of functional interactions through effective connectivity requires the determination of three distinct levels of information: (1) the regions involved in the process and forming the spatial support of the network, (2) the presence or absence of interactions between each pair of regions, and (3) the directionality of the existing interactions. While many methods exist to select regions (Step 1), very little is available to complete Step 2. The two main methods developed so far, structural equation modeling (SEM) and dynamical causal modeling (DCM), usually require precise prior information to be used, while such information is sometimes lacking. Assuming that Step 1 was successfully completed, we here propose a data-driven method to deal with Step 2 and extract functional interactions from fMRI datasets through partial correlations. Partial correlation is more closely related to effective connectivity than marginal correlation and provides a convenient graphical representation for functional interactions. As an instance of brain interactivity investigation, we consider how simple hand movements are processed by the bihemispheric cortical motor network. In the proposed framework, Bayesian analysis makes it possible to estimate and test the partial statistical dependencies between regions without any prior model on the underlying functional interactions. We demonstrate the interest of this approach on real data.

Cerebral Cortex↗

PubMed2006

Bayesian spatio-temporal approach for EEG source reconstruction: conciliating ECD and distributed models.

Characterizing the cortical activity sources of electroencephalography (EEG)/magnetoencephalography data is a critical issue since it requires solving an ill-posed inverse problem that does not admit a unique solution. Two main different and complementary source models have emerged: equivalent current dipoles (ECD) and distributed linear (DL) models. While ECD models remain highly popular since they provide an easy way to interpret the solutions, DL models (also referred to as imaging techniques) are known to be more realistic and flexible. In this paper, we show how those two representations of the brain electromagnetic activity can be cast into a common general framework yielding an optimal description and estimation of the EEG sources. From this extended source mixing model, we derive a hybrid approach whose key aspect is the separation between temporal and spatial characteristics of brain activity, which allows to dramatically reduce the number of DL model parameters. Furthermore, the spatial profile of the sources, as a temporal invariant map, is estimated using the entire time window data, allowing to significantly enhance the information available about the spatial aspect of the EEG inverse problem. A Bayesian framework is introduced to incorporate distinct temporal and spatial constraints on the solution and to estimate both parameters and hyperparameters of the model. Using simulated EEG data, the proposed inverse approach is evaluated and compared with standard distributed methods using both classical criteria and ROC curves.

Algorithms↗

PubMed2005

Identification of large-scale networks in the brain using fMRI.

Cognition is thought to result from interactions within large-scale networks of brain regions. Here, we propose a method to identify these large-scale networks using functional magnetic resonance imaging (fMRI). Regions belonging to such networks are defined as sets of strongly interacting regions, each of which showing a homogeneous temporal activity. Our method of large-scale network identification (LSNI) proceeds by first detecting functionally homogeneous regions. The networks of functional interconnections are then found by comparing the correlations among these regions against a model of the correlations in the noise. To test the LSNI method, we first evaluated its specificity and sensitivity on synthetic data sets. Then, the method was applied to four real data sets with a block-designed motor task. The LSNI method correctly recovered the regions whose temporal activity was locked to the stimulus. In addition, it detected two other main networks highly reproducible across subjects, whose activity was dominated by slow fluctuations (0-0.1 Hz). One was located in medial and dorsal regions, and mostly overlapped the "default" network of the brain at rest [Greicius, M.D., Krasnow, B., Reiss, A.L., Menon, V., 2003. Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proceedings of the National Academy of Sciences of the U.S.A. 100, 253-258]; the other was composed of lateral frontal and posterior parietal regions. The LSNI method we propose allows to detect in an exploratory and systematic way all the regions and large-scale networks activated in the working brain.

Adult↗

PubMed2005

Distinct basal ganglia territories are engaged in early and advanced motor sequence learning.

In this study, we used functional MRI (fMRI) at high field (3T) to track the time course of activation in the entire basal ganglia circuitry, as well as other motor-related structures, during the explicit learning of a sequence of finger movements over a month of training. Fourteen right-handed healthy volunteers had to practice 15 min daily a sequence of eight moves using the left hand. MRI sessions were performed on days 1, 14 and 28. In both putamen, activation decreased with practice in rostrodorsal (associative) regions. In contrast, there was a significant signal increase in more caudoventral (sensorimotor) regions of the putamen. Subsequent correlation analyses between signal variations and behavioral variables showed that the error rate (movement accuracy) was positively correlated with signal changes in areas activated during early learning, whereas reaction time (movement speed) was negatively correlated with signal changes in areas activated during advanced learning stages, including the sensorimotor putamen and globus pallidus. These results suggest the possibility that motor representations shift from the associative to the sensorimotor territories of the striato-pallidal complex during the explicit learning of motor sequences, suggesting that motor skills are stored in the sensorimotor territory of the basal ganglia that supports a speedy performance.

Adult↗

PubMed2005

Simulation of anisotropic growth of low-grade gliomas using diffusion tensor imaging.

A recent computational model of brain tumor growth, developed to better describe how gliomas invade through the adjacent brain parenchyma, is based on two major elements: cell proliferation and isotropic cell diffusion. On the basis of this model, glioma growth has been simulated in a virtual brain, provided by a 3D segmented MRI atlas. However, it is commonly accepted that glial cells preferentially migrate along the direction of fiber tracts. Therefore, in this paper, the model has been improved by including anisotropic extension of gliomas. The method is based on a cell diffusion tensor derived from water diffusion tensor (as given by MRI diffusion tensor imaging). Results of simulations have been compared with two clinical examples demonstrating typical growth patterns of low-grade gliomas centered around the insula. The shape and the kinetic evolution are better simulated with anisotropic rather than isotropic diffusion. The best fit is obtained when the anisotropy of the cell diffusion tensor is increased to greater anisotropy than the observed water diffusion tensor. The shape of the tumor is also influenced by the initial location of the tumor. Anisotropic brain tumor growth simulations provide a means to determine the initial location of a low-grade glioma as well as its cell diffusion tensor, both of which might reflect the biological characteristics of invasion.

Anisotropy↗

PubMed2005

From anatomic standardization analysis of perfusion SPECT data to perfusion pattern modeling: evidence of functional networks in healthy subjects and temporal lobe epilepsy patients.

RATIONALE AND OBJECTIVES: In the general context of perfusion pattern modeling from single-photon emission computed tomographic (SPECT) data, the purpose of this study is to characterize interindividual functional variability and functional connectivity between anatomic structures in a set of SPECT data acquired from a homogeneous population of subjects. MATERIALS AND METHODS: From volume of interest (VOI)-perfusion measurements performed on anatomically standardized SPECT data, we proposed to use correspondence analysis (CA) and hierarchical clustering (HC) to explore the structure of statistical dependencies among these measurements. The method was applied to study the perfusion pattern in two populations of subjects; namely, SPECT data from 27 healthy subjects and ictal SPECT data from 10 patients with mesio-temporal lobe epilepsy (MTLE). RESULTS: For healthy subjects, anatomic structures showing statistically dependent perfusion patterns were classified into four groups; namely, temporomesial structures, internal structures, posterior structures, and remaining cortex. For patients with MTLE, they were classified as temporomesial structures, surrounding temporal structures, internal structures, and remaining cortex. Anatomic structures of each group showed similar perfusion behavior so that they may be functionally connected and may belong to the same network. Our main result is that the temporal pole and lenticular nucleus seemed to be highly relevant to characterize ictal perfusion in patients with MTLE. This exploratory analysis suggests that a network involving temporal structures, lenticular nucleus, brainstem, and cerebellum seems to be involved during MTLE seizures. CONCLUSION: CA followed by HC is a promising approach to explore brain perfusion patterns from SPECT VOI measurements.

Adult↗

PubMed2005

Reorganization and plasticity in the adult brain during learning of motor skills.

On the basis of brain imaging studies, Doyon and Ungerleider recently proposed a model describing the cerebral plasticity that occurs in both cortico-striatal and cortico-cerebellar systems of the adult brain during learning of new motor skilled behaviors. This theoretical framework makes several testable predictions with regards to the contribution of these neural systems based on the phase (fast, slow, consolidation, automatization, and retention) and nature of the motor learning processes (motor sequence versus motor adaptation) acquired through repeated practice. There has been recent behavioral, lesion and additional neuroimaging studies that have addressed the assumptions made in this theory that will help in the revision of this model.

Adult↗

PubMed2004

Estimation of the hemodynamic response in event-related functional MRI: Bayesian networks as a framework for efficient Bayesian modeling and inference.

A convenient way to analyze blood-oxygen-level-dependent functional magnetic resonance imaging data consists of modeling the whole brain as a stationary, linear system characterized by its transfer function: the hemodynamic response function (HRF). HRF estimation, though of the greatest interest, is still under investigation, for the problem is ill-conditioned. In this paper, we recall the most general Bayesian model for HRF estimation and show how it can beneficially be translated in terms of Bayesian graphical models, leading to 1) a clear and efficient representation of all structural and functional relationships entailed by the model, and 2) a straightforward numerical scheme to approximate the joint posterior distribution, allowing for estimation of the HRF, as well as all other model parameters. We finally apply this novel technique on both simulations and real data.

Adolescent↗

PubMed2003

Robust Bayesian estimation of the hemodynamic response function in event-related BOLD fMRI using basic physiological information.

In BOLD fMRI data analysis, robust and accurate estimation of the Hemodynamic Response Function (HRF) is still under investigation. Parametric methods assume the shape of the HRF to be known and constant throughout the brain, whereas non-parametric methods mostly rely on artificially increasing the signal-to-noise ratio. We extend and develop a previously proposed method that makes use of basic yet relevant temporal information about the underlying physiological process of the brain BOLD response in order to infer the HRF in a Bayesian framework. A general hypothesis test is also proposed, allowing to take advantage of the knowledge gained regarding the HRF to perform activation detection. The performances of the method are then evaluated by simulation. Great improvement is shown compared to the Maximum-Likelihood estimate in terms of estimation error, variance, and bias. Robustness of the estimators with regard to the actual noise structure or level, as well as the stimulus sequence, is also proven. Lastly, fMRI data with an event-related paradigm are analyzed. As suspected, the regions selected from highly discriminating activation maps resulting from the method exhibit a certain inter-regional homogeneity in term of HRF shape, as well as noticeable inter-regional differences.

Adult↗

PubMed2003

Localization Estimation Algorithm (LEA): a supervised prior-based approach for solving the EEG/MEG inverse problem.

Localizing and quantifying the sources of ElectroEncephaloGraphy (EEG) and MagnetoEncephaloGraphy (MEG) measurements is an ill-posed inverse problem, whose solution requires a spatial regularization involving both anatomical and functional priors. The distributed source model enables the introduction of such constraints. However, the resulting solution is unstable since the equation system one has to solve is badly conditioned and under-determined. We propose an original approach for solving the inverse problem, that allows to deal with a better-determined system and to temper the influence of priors according to their consistency with the measured EEG/MEG data. This Localization Estimation Algorithm (LEA) estimates the amplitude of a selected subset of sources, which are localized based on a prior distribution of activation probability. LEA is evaluated through numerical simulations and compared to a classical Weighted Minimum Norm estimation.

Algorithms↗

PubMed2003

Multivariate group effect analysis in functional Magnetic Resonance Imaging.

In functional MRI (fMRI), analysis of multisubject data typically involves spatially normalizing (i.e. co-registering in a common standard space) all data sets and summarizing results in a single group activation map. This widely used approach does not explicitely account for between-subject anatomo-functional variability. Therefore, we propose a group effect analysis method which makes use of a multivariate model to select the main signal variations that are common to all subjects, while allowing final statistical inference on the individual scale. The normalization step is thus avoided and individual anatomo-functional features are preserved. The approach is evaluated by using simulated data and it is shown that sensitivity is drastically improved compared to more conventional individual analysis.

Algorithms↗

PubMed2003

Estimation of the Hemodynamic Response Function in event-related functional MRI: directed acyclic graphs for a general Bayesian inference framework.

A convenient way to analyze BOLD fMRI data consists of modeling the whole brain as a stationary, linear system characterized by its transfer function: the Hemodynamic Response Function (HRF). HRF estimation, though of the greatest interest, is still under investigation, for the problem is ill-conditioned. In this paper, we recall the most general Bayesian model for HRF estimation and show how it can beneficially be translated in terms of graphical models, leading to (i) a clear and efficient representation of all structural and functional relationships entailed by the model, and (ii) a straightforward numerical scheme to approximate the joint posterior distribution, allowing for estimation of the HRF, as well as all other model parameters. We finally apply this novel technique on both simulations and real data.

Adolescent↗

PubMed2003

Group analysis in functional neuroimaging: selecting subjects using similarity measures.

Standard group analyses of fMRI data rely on spatial and temporal averaging of individuals. This averaging operation is only sensible when the mean is a good representation of the group. This is not the case if subjects are not homogeneous, and it is therefore a major concern in fMRI studies to assess this group homogeneity. We present a method that provides relevant distances or similarity measures between temporal series of brain functional images belonging to different subjects. The method allows a multivariate comparison between data sets of several subjects in the time or in the space domain. These analyses assess the global intersubject variability before averaging subjects and drawing conclusions across subjects, at the population level. We adapt the RV coefficient to measure meaningful spatial or temporal similarities and use multidimensional scaling to give a visual representation of each subject's position with respect to other subjects in the group. We also provide a measure for detecting subjects that may be outliers. Results show that the method is a powerful tool to detect subjects with specific temporal or spatial patterns, and that, despite the apparent loss of information, restricting the analysis to a homogeneous subgroup of subjects does not reduce the statistical sensitivity of standard group fMRI analyses.

Algorithms↗

PubMed2003

Unsupervised robust nonparametric estimation of the hemodynamic response function for any fMRI experiment.

This paper deals with the estimation of the blood oxygen level-dependent response to a stimulus, as measured in functional magnetic resonance imaging (fMRI) data. A precise estimation is essential for a better understanding of cerebral activations. The most recent works have used a nonparametric framework for this estimation, considering each brain region as a system characterized by its impulse response, the so-called hemodynamic response function (HRF). However, the use of these techniques has remained limited since they are not well-adapted to real fMRI data. Here, we develop a threefold extension to previous works. We consider asynchronous event-related paradigms, account for different trial types and integrate several fMRI sessions into the estimation. These generalizations are simultaneously addressed through a badly conditioned observation model. Bayesian formalism is used to model temporal prior information of the underlying physiological process of the brain hemodynamic response. By this way, the HRF estimate results from a tradeoff between information brought by the data and by our prior knowledge. This tradeoff is modeled with hyperparameters that are set to the maximum-likelihood estimate using an expectation conditional maximization algorithm. The proposed unsupervised approach is validated on both synthetic and real fMRI data, the latter originating from a speech perception experiment.

Algorithms↗

PubMed2002

Multivariate model specification for fMRI data.

We present a general method-denoted MoDef-to help specify (or define) the model used to analyze brain imaging data. This method is based on the use of the multivariate linear model on a training data set. We show that when the a priori knowledge about the expected brain response is not too precise, the method allows for the specification of a model that yields a better sensitivity in the statistical results. This obviously relies on the validity of the a priori information, in our case the representativity of the training set, an issue addressed using a cross-validation technique. We propose a fast implementation that allows the use of the method on large data sets as found with functional Magnetic Resonance Images. An example of application is given on an experimental fMRI data set that includes nine subjects who performed a mental computation task. Results show that the method increases the statistical sensitivity of fMRI analyses.

Brain↗

PubMed2002

ICA of fMRI group study data.

This paper proposes to extend independent component analysis (ICA) of functional magnetic resonance imaging (fMRI) data from single subjects to simultaneous analysis of data from a group of subjects. This results in a set of time courses which are common to the whole group, together with an individual spatial response pattern for each of the subjects in the group. The method is illustrated using data from two fMRI experiments. The results show that: (a) ICA is capable of extracting nontrivial task related components without any a priori information about the fMRI experiment; (b) in analysis of group data, ICA identifies components common to the whole group as well as components manifested in single subjects only.

Brain↗