Bayesian estimation of orientation preference maps

Imaging techniques such as optical imaging of intrinsic signals, 2-photon calcium imaging and voltage sensitive dye imaging can be used to measure the functional organization of visual cortex across different spatial and temporal scales. Here, we present Bayesian methods based on Gaussian processes for extracting topographic maps from functional imaging data. In particular, we focus on the estimation of orientation preference maps (OPMs) from intrinsic signal imaging data. We model the underlying map as a bivariate Gaussian process, with a prior covariance function that reflects known properties of OPMs, and a noise covariance adjusted to the data. The posterior mean can be interpreted as an optimally smoothed estimate of the map, and can be used for model based interpolations of the map from sparse measurements. By sampling from the posterior distribution, we can get error bars on statistical properties such as preferred orientations, pinwheel locations or pinwheel counts. Finally, the use of an explicit probabilistic model facilitates interpretation of parameters and quantitative model comparisons. We demonstrate our model both on simulated data and on intrinsic signaling data from ferret visual cortex

Denoising Two-Photon Calcium Imaging Data

Two-photon calcium imaging is now an important tool for in vivo imaging of biological systems. By enabling neuronal population imaging with subcellular resolution, this modality offers an approach for gaining a fundamental understanding of brain anatomy and physiology. Proper analysis of calcium imaging data requires denoising, that is separating the signal from complex physiological noise. To analyze two-photon brain imaging data, we present a signal plus colored noise model in which the signal is represented as harmonic regression and the correlated noise is represented as an order autoregressive process. We provide an efficient cyclic descent algorithm to compute approximate maximum likelihood parameter estimates by combing a weighted least-squares procedure with the Burg algorithm. We use Akaike information criterion to guide selection of the harmonic regression and the autoregressive model orders. Our flexible yet parsimonious modeling approach reliably separates stimulus-evoked fluorescence response from background activity and noise, assesses goodness of fit, and estimates confidence intervals and signal-to-noise ratio. This refined separation leads to appreciably enhanced image contrast for individual cells including clear delineation of subcellular details and network activity. The application of our approach to in vivo imaging data recorded in the ferret primary visual cortex demonstrates that our method yields substantially denoised signal estimates. We also provide a general Volterra series framework for deriving this and other signal plus correlated noise models for imaging. This approach to analyzing two-photon calcium imaging data may be readily adapted to other computational biology problems which apply correlated noise models.National Institutes of Health (U.S.) (DP1 OD003646-01)National Institutes of Health (U.S.) (R01EB006385-01)National Institutes of Health (U.S.) (EY07023)National Institutes of Health (U.S.) (EY017098

Les cartes fonctionnelles dans le cortex visuel du chat : nouvelles stratégies d’évaluation en imagerie optique et mise en évidence de l’organisation anatomo-fonctionnelle

Le regroupement des neurones de propriétés similaires est à l’origine de modules permettant d’optimiser l’analyse de l’information. La conséquence est la présence de cartes fonctionnelles dans le cortex visuel primaire de certains mammifères pour de nombreux paramètres tels que l’orientation, la direction du mouvement ou la position des stimuli (visuotopie). Le premier volet de cette thèse est consacré à caractériser l’organisation modulaire dans le cortex visuel primaire pour un paramètre fondamental, la suppression centre / pourtour et au delà du cortex visuel primaire (dans l’aire 21a), pour l’orientation et la direction. Toutes les études ont été effectuées à l’aide de l’imagerie optique des signaux intrinsèques sur le cortex visuel du chat anesthésié. La quantification de la modulation par la taille des stimuli à permis de révéler la présence de modules de forte et de faible suppression par le pourtour dans le cortex visuel primaire (aires 17 et 18). Ce type d’organisation n’avait été observé jusqu’ici que dans une aire de plus haut niveau hiérarchique chez le primate. Une organisation modulaire pour l’orientation, similaire à celle observée dans le cortex visuel primaire a été révélée dans l’aire 21a. Par contre, contrairement à l’aire 18, l’aire 21a ne semblait pas être organisée en domaine de direction. L’ensemble de ces résultats pourront permettre d’alimenter les connaissances sur l’organisation anatomo-fonctionnelle du cortex visuel du chat mais également de mieux comprendre les facteurs qui déterminent la présence d’une organisation modulaire. Le deuxième volet abordé dans cette thèse s’est intéressé à l’amélioration de l’aspect quantitatif apporté par l’analyse temporelle en imagerie optique des signaux intrinsèques. Cette nouvelle approche, basée sur l’analyse de Fourier a permis d’augmenter considérablement le rapport signal / bruit des enregistrements. Toutefois, cette analyse ne s’est basée jusqu’ici que sur la quantification d’une seule harmonique ce qui a limité son emploi à la cartographie de l’orientation et de rétinotopie uniquement. En exploitant les plus hautes harmoniques, un modèle a été proposé afin d’estimer la taille des champs récepteurs et la sélectivité à la direction. Ce modèle a par la suite été validé par des approches conventionnelles dans le cortex visuel primaire.The clustering of neurons of similar properties is at the basis of the brain modular architecture and is considered as a strategy to optimized processing. One consequence of this clustering is the presence of functional maps in the primary visual cortex of several mammals based on features such as orientation, direction of motion and stimulus position (retinotopy). The first section of this thesis was aimed at characterizing the modular organization of functions in primary and higher-order areas. First, we investigated the possibility that a fundamental cell property, the receptive field center / surround suppression, could be orderly represented in the primary visual cortex. Second, we determined the level of modular organization in area 21a for two key properties, orientation and direction of motion. All studies were based on the optical imaging of intrinsic signals in anesthetized cats. Results indicate the presence of high and low surround suppression modules in the primary visual cortex (areas 17 and 18). To date, such organization has been discovered only in a higher-order area in primate. A modular organization for orientation, similar to the one observed in areas 17 and 18 was observed in area 21a. On the other hand, in contrast to area 18, no direction modules were discovered in area 21a. Overall, the first part of this thesis increased our knowledge about the anatomo-fonctional organization of cat visual cortex. They will also be instrumental to better understand the factors leading to the presence of a modular organization in the cortex. The second section of this thesis was directed to the development of a novel quantitative tool for the temporal analysis of optical imaging intrinsic signals. This new approach, based on Fourier decomposition, allowed to greatly increase the signal to noise ratio of the recordings. Until now, this analysis was only been based on single harmonic quantification, limiting its application for orientation and rétinotopy mapping only. A model exploiting higher harmonics was then developed to estimate additional parameters such as the receptive field size and direction selectivity. Thereafter, this model was validated with success by conventional approaches on the primary visual cortex

Abstract Spatiotemporal analysis of optical imaging data

Previous methods for analyzing optical imaging data have relied heavily on temporal averaging. However, response dynamics are rich sources of information. Here, we develop and present a method that combines principal component analysis and multitaper harmonic analysis to extract the statistically significant spatial and temporal response from optical imaging data. We apply the method to both simulated data and experimental optical imaging data from the cat primary visual cortex. © 2003 Elsevier Science (USA). All rights reserved