System identification of the nonlinear dynamics in the thalamocortical circuit in response to patterned thalamic microstimulation in vivo

• Type of Material: Article
• Authors: Daniel C Millard, Emory University; Qi Wang, Emory University; Clare A Gollnick, Emory University; Garrett Stanley, Emory University
• Language: English
• Date: 2013-12-01
• Publisher: IOP Publishing: Hybrid Open Access
• Publication Version: Author Accepted Manuscript (After Peer Review)
• Copyright Statement: © 2013 IOP Publishing Ltd.
• Final Published Version (URL): http://dx.doi.org/10.1088/1741-2560/10/6/066011
• Title of Journal or Parent Work: Journal of Neural Engineering
• ISSN: 1741-2560
• Volume: 10
• Issue: 6
• Start Page: 066011
• End Page: 066011
• Grant/Funding Information: This work was supported by the National Institutes of Health (NIH Grant R01NS48285); and D.C.M. and C.A.G. were each supported by National Science Foundation Graduate Research Fellowships.
• Keywords: Neurosciences & Neurology; AREA MT; ELECTRICAL-STIMULATION; LAYER 2/3; CORTICAL MICROSTIMULATION; Life Sciences & Biomedicine; Technology; SENSORY RESPONSES; BARREL CORTEX; Engineering, Biomedical; INTRACORTICAL MICROSTIMULATION; Neurosciences; BIOLOGICAL-SYSTEMS; RECEPTIVE-FIELD; Science & Technology; PRIMARY VISUAL-CORTEX; Engineering
• Research Categories: Biology, Neuroscience; Engineering, Biomedical

Abstract

Objective. Nonlinear system identification approaches were used to develop a dynamical model of the network level response to patterns of microstimulation in vivo. Approach. The thalamocortical circuit of the rodent vibrissa pathway was the model system, with voltage sensitive dye imaging capturing the cortical response to patterns of stimulation delivered from a single electrode in the ventral posteromedial thalamus. The results of simple paired stimulus experiments formed the basis for the development of a phenomenological model explicitly containing nonlinear elements observed experimentally. The phenomenological model was fit using datasets obtained with impulse train inputs, Poisson-distributed in time and uniformly varying in amplitude. Main results. The phenomenological model explained 58% of the variance in the cortical response to out of sample patterns of thalamic microstimulation. Furthermore, while fit on trial-averaged data, the phenomenological model reproduced single trial response properties when simulated with noise added into the system during stimulus presentation. The simulations indicate that the single trial response properties were dependent on the relative sensitivity of the static nonlinearities in the two stages of the model, and ultimately suggest that electrical stimulation activates local circuitry through linear recruitment, but that this activity propagates in a highly nonlinear fashion to downstream targets. Significance. The development of nonlinear dynamical models of neural circuitry will guide information delivery for sensory prosthesis applications, and more generally reveal properties of population coding within neural circuits.