An ARX model-based approach to trial by trial identification of fMRI-BOLD responses

Identification of Negative BOLD Responses in Epilepsy Using Windkessel Models

Alongside positive blood oxygenation level–dependent (BOLD) responses associated with interictal epileptic discharges, a variety of negative BOLD responses (NBRs) are typically found in epileptic patients. Previous studies suggest that, in general, up to four mechanisms might underlie the genesis of NBRs in the brain: (i) neuronal disruption of network activity, (ii) altered balance of neurometabolic/vascular couplings, (iii) arterial blood stealing, and (iv) enhanced cortical inhibition. Detecting and classifying these mechanisms from BOLD signals are pivotal for the improvement of the specificity of the electroencephalography–functional magnetic resonance imaging (EEG-fMRI) image modality to identify the seizure-onset zones in refractory local epilepsy. This requires models with physiological interpretation that furnish the understanding of how these mechanisms are fingerprinted by their BOLD responses. Here, we used a Windkessel model with viscoelastic compliance/inductance in combination with dynamic models of both neuronal population activity and tissue/blood O₂ to classify the hemodynamic response functions (HRFs) linked to the above mechanisms in the irritative zones of epileptic patients. First, we evaluated the most relevant imprints on the BOLD response caused by variations of key model parameters. Second, we demonstrated that a general linear model is enough to accurately represent the four different types of NBRs. Third, we tested the ability of a machine learning classifier, built from a simulated ensemble of HRFs, to predict the mechanism underlying the BOLD signal from irritative zones. Cross-validation indicates that these four mechanisms can be classified from realistic fMRI BOLD signals. To demonstrate proof of concept, we applied our methodology to EEG-fMRI data from five epileptic patients undergoing neurosurgery, suggesting the presence of some of these mechanisms. We concluded that a proper identification and interpretation of NBR mechanisms in epilepsy can be performed by combining general linear models and biophysically inspired models.

State-space models of impulse hemodynamic responses over motor, somatosensory, and visual cortices

THE PAPER PRESENTS STATE SPACE MODELS OF THE HEMODYNAMIC RESPONSE (HR) OF FNIRS TO AN IMPULSE STIMULUS IN THREE BRAIN REGIONS: motor cortex (MC), somatosensory cortex (SC), and visual cortex (VC). Nineteen healthy subjects were examined. For each cortex, three impulse HRs experimentally obtained were averaged. The averaged signal was converted to a state space equation by using the subspace method. The activation peak and the undershoot peak of the oxy-hemoglobin (HbO) in MC are noticeably higher than those in SC and VC. The time-to-peaks of the HbO in three brain regions are almost the same (about 6.76 76 ± 0.2 s). The time to undershoot peak in VC is the largest among three. The HbO decreases in the early stage (~0.46 s) in MC and VC, but it is not so in SC. These findings were well described with the developed state space equations. Another advantage of the proposed method is its easy applicability in generating the expected HR to arbitrary stimuli in an online (or real-time) imaging. Experimental results are demonstrated.

A Framework for Networked Experiments in Global E-Science

This chapter presents a framework that creates, uses, and communicates information whose organizational dynamics allow individuals to perform a distributed cooperative enterprise in public educational environments. The approach presented here assumes Web services (possibly offered over a grid) are the enacting paradigm used to formalize educational interactions as cooperative services on various computational nodes of a network. By examining a case study involving a well known micro-array experiment in the growing field of bioinformatics, this chapter will detail how specific classes of interactions can be mapped into a service-oriented model that can be implemented in a variety of e-learning contexts. This framework illustrated by this case study allows for a sophisticated degree of e-learning that can be applied to a range of local or international contexts.

A stochastic linear model for fMRI activation analyses

PURPOSE

The debate regarding how best to model variability of the hemodynamic response function in fMRI data has focussed on the linear vs. nonlinear nature of the optimal signal model, with few studies exploring the deterministic vs. stochastic nature of the dynamics. We propose a stochastic linear model (SLM) of the hemodynamic signal and noise dynamics to more robustly infer fMRI activation estimates.

METHODS

The SLM models the hemodynamic signal by an exogenous input autoregressive model driven by Gaussian state noise. Activation weights are inferred by a joint state-parameter iterative coordinate descent algorithm based on the Kalman smoother.

RESULTS

The SLM produced more accurate parameter estimates than the GLM for event-design simulated data. In application to block-design experimental visuo-motor task fMRI data, the SLM resulted in more punctate and well-defined motor cortex activation maps than the GLM, and was able to track variations in the hemodynamics, as expected from a stochastic model.

CONCLUSIONS

We demonstrate in application to both simulated and experimental fMRI data that in comparison to the GLM, the SLM produces more flexible, consistent and enhanced fMRI activation estimates.

Biomedical signal and image processing

Generally, physiological modeling and biomedical signal processing constitute two important paradigms of biomedical engineering (BME): their fundamental concepts are taught starting from undergraduate studies and are more completely dealt with in the last years of graduate curricula, as well as in Ph.D. courses. Traditionally, these two cultural aspects were separated, with the first one more oriented to physiological issues and how to model them and the second one more dedicated to the development of processing tools or algorithms to enhance useful information from clinical data. A practical consequence was that those who did models did not do signal processing and vice versa. However, in recent years,the need for closer integration between signal processing and modeling of the relevant biological systems emerged very clearly [1], [2]. This is not only true for training purposes(i.e., to properly prepare the new professional members of BME) but also for the development of newly conceived research projects in which the integration between biomedical signal and image processing (BSIP) and modeling plays a crucial role. Just to give simple examples, topics such as brain–computer machine or interfaces,neuroengineering, nonlinear dynamical analysis of the cardiovascular (CV) system,integration of sensory-motor characteristics aimed at the building of advanced prostheses and rehabilitation tools, and wearable devices for vital sign monitoring and others do require an intelligent fusion of modeling and signal processing competences that are certainly peculiar of our discipline of BME.

Spatially adaptive mixture modeling for analysis of FMRI time series

Within-subject analysis in fMRI essentially addresses two problems, the detection of brain regions eliciting evoked activity and the estimation of the underlying dynamics. In Makni et aL, 2005 and Makni et aL, 2008, a detection-estimation framework has been proposed to tackle these problems jointly, since they are connected to one another. In the Bayesian formalism, detection is achieved by modeling activating and nonactivating voxels through independent mixture models (IMM) within each region while hemodynamic response estimation is performed at a regional scale in a nonparametric way. Instead of IMMs, in this paper we take advantage of spatial mixture models (SMM) for their nonlinear spatial regularizing properties. The proposed method is unsupervised and spatially adaptive in the sense that the amount of spatial correlation is automatically tuned from the data and this setting automatically varies across brain regions. In addition, the level of regularization is specific to each experimental condition since both the signal-to-noise ratio and the activation pattern may vary across stimulus types in a given brain region. These aspects require the precise estimation of multiple partition functions of underlying Ising fields. This is addressed efficiently using first path sampling for a small subset of fields and then using a recently developed fast extrapolation technique for the large remaining set. Simulation results emphasize that detection relying on supervised SMM outperforms its IMM counterpart and that unsupervised spatial mixture models achieve similar results without any hand-tuning of the correlation parameter. On real datasets, the gain is illustrated in a localizer fMRI experiment: brain activations appear more spatially resolved using SMM in comparison with classical general linear model (GLM)-based approaches, while estimating a specific parcel-based HRF shape. Our approach therefore validates the treatment of unsmoothed fMRI data without fixed GLM definition at the subject level and makes also the classical strategy of spatial Gaussian filtering deprecated.

Biomedical applications of Piece-Wise Affine identification for hybrid systems

Modeling switching processes for control purpose takes advantage from the Piece-Wise Affine identification of hybrid dynamical systems briefly recalled in this paper. A couple of applications are addressed, namely to discriminate hormone pulses from background noise, in a physiologically switching process, and to identify sleep apneas, as pathological switching among healthy and potentially risky states. Other potential applications are proposed.

Model estimation of cerebral hemodynamics between blood flow and volume changes: a data-based modeling approach

It is well known that there is a dynamic relationship between cerebral blood flow (CBF) and cerebral blood volume (CBV). With increasing applications of functional MRI, where the blood oxygen-level-dependent signals are recorded, the understanding and accurate modeling of the hemodynamic relationship between CBF and CBV becomes increasingly important. This study presents an empirical and data-based modeling framework for model identification from CBF and CBV experimental data. It is shown that the relationship between the changes in CBF and CBV can be described using a parsimonious autoregressive with exogenous input model structure. It is observed that neither the ordinary least-squares (LS) method nor the classical total least-squares (TLS) method can produce accurate estimates from the original noisy CBF and CBV data. A regularized total least-squares (RTLS) method is thus introduced and extended to solve such an error-in-the-variables problem. Quantitative results show that the RTLS method works very well on the noisy CBF and CBV data. Finally, a combination of RTLS with a filtering method can lead to a parsimonious but very effective model that can characterize the relationship between the changes in CBF and CBV.