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

Type-specific dendritic integration in mouse retinal ganglion cells

Neural computation relies on the integration of synaptic inputs across a neuron’s dendritic arbour. However, it is far from understood how different cell types tune this process to establish cell-type specific computations. Here, using two-photon imaging of dendritic Ca2+ signals, electrical recordings of somatic voltage and biophysical modelling, we demonstrate that four morphologically distinct types of mouse retinal ganglion cells with overlapping excitatory synaptic input (transient Off alpha, transient Off mini, sustained Off, and F-mini Off) exhibit type-specific dendritic integration profiles: in contrast to the other types, dendrites of transient Off alpha cells were spatially independent, with little receptive field overlap. The temporal correlation of dendritic signals varied also extensively, with the highest and lowest correlation in transient Off mini and transient Off alpha cells, respectively. We show that differences between cell types can likely be explained by differences in backpropagation efficiency, arising from the specific combinations of dendritic morphology and ion channel densities

Striatal Proteomic Analysis Suggests that First L-Dopa Dose Equates to Chronic Exposure

L-3,4-dihydroxypheylalanine (L-dopa)-induced dyskinesia represent a debilitating complication of therapy for Parkinson's disease (PD) that result from a progressive sensitization through repeated L-dopa exposures. The MPTP macaque model was used to study the proteome in dopamine-depleted striatum with and without subsequent acute and chronic L-dopa treatment using two-dimensional difference in-gel electrophoresis (2D-DIGE) and mass spectrometry. The present data suggest that the dopamine-depleted striatum is so sensitive to de novo L-dopa treatment that the first ever administration alone would be able (i) to induce rapid post-translational modification-based proteomic changes that are specific to this first exposure and (ii), possibly, lead to irreversible protein level changes that would be not further modified by chronic L-dopa treatment. The apparent equivalence between first and chronic L-dopa administration suggests that priming would be the direct consequence of dopamine loss, the first L-dopa administrations only exacerbating the sensitization process but not inducing it

The Effects of NMDA Subunit Composition on Calcium Influx and Spike Timing-Dependent Plasticity in Striatal Medium Spiny Neurons

Calcium through NMDA receptors (NMDARs) is necessary for the long-term potentiation (LTP) of synaptic strength; however, NMDARs differ in several properties that can influence the amount of calcium influx into the spine. These properties, such as sensitivity to magnesium block and conductance decay kinetics, change the receptor's response to spike timing dependent plasticity (STDP) protocols, and thereby shape synaptic integration and information processing. This study investigates the role of GluN2 subunit differences on spine calcium concentration during several STDP protocols in a model of a striatal medium spiny projection neuron (MSPN). The multi-compartment, multi-channel model exhibits firing frequency, spike width, and latency to first spike similar to current clamp data from mouse dorsal striatum MSPN. We find that NMDAR-mediated calcium is dependent on GluN2 subunit type, action potential timing, duration of somatic depolarization, and number of action potentials. Furthermore, the model demonstrates that in MSPNs, GluN2A and GluN2B control which STDP intervals allow for substantial calcium elevation in spines. The model predicts that blocking GluN2B subunits would modulate the range of intervals that cause long term potentiation. We confirmed this prediction experimentally, demonstrating that blocking GluN2B in the striatum, narrows the range of STDP intervals that cause long term potentiation. This ability of the GluN2 subunit to modulate the shape of the STDP curve could underlie the role that GluN2 subunits play in learning and development

Neural Computation via Neural Geometry: A Place Code for Inter-whisker Timing in the Barrel Cortex?

The place theory proposed by Jeffress (1948) is still the dominant model of how the brain represents the movement of sensory stimuli between sensory receptors. According to the place theory, delays in signalling between neurons, dependent on the distances between them, compensate for time differences in the stimulation of sensory receptors. Hence the location of neurons, activated by the coincident arrival of multiple signals, reports the stimulus movement velocity. Despite its generality, most evidence for the place theory has been provided by studies of the auditory system of auditory specialists like the barn owl, but in the study of mammalian auditory systems the evidence is inconclusive. We ask to what extent the somatosensory systems of tactile specialists like rats and mice use distance dependent delays between neurons to compute the motion of tactile stimuli between the facial whiskers (or ‘vibrissae’). We present a model in which synaptic inputs evoked by whisker deflections arrive at neurons in layer 2/3 (L2/3) somatosensory ‘barrel’ cortex at different times. The timing of synaptic inputs to each neuron depends on its location relative to sources of input in layer 4 (L4) that represent stimulation of each whisker. Constrained by the geometry and timing of projections from L4 to L2/3, the model can account for a range of experimentally measured responses to two-whisker stimuli. Consistent with that data, responses of model neurons located between the barrels to paired stimulation of two whiskers are greater than the sum of the responses to either whisker input alone. The model predicts that for neurons located closer to either barrel these supralinear responses are tuned for longer inter-whisker stimulation intervals, yielding a topographic map for the inter-whisker deflection interval across the surface of L2/3. This map constitutes a neural place code for the relative timing of sensory stimuli

Changing the responses of cortical neurons from sub- to suprathreshold using single spikes in vivo.

Action Potential (APs) patterns of sensory cortex neurons encode a variety of stimulus features, but how can a neuron change the feature to which it responds? Here, we show that in vivo a spike-timing-dependent plasticity (STDP) protocol-consisting of pairing a postsynaptic AP with visually driven presynaptic inputs-modifies a neurons' AP-response in a bidirectional way that depends on the relative AP-timing during pairing. Whereas postsynaptic APs repeatedly following presynaptic activation can convert subthreshold into suprathreshold responses, APs repeatedly preceding presynaptic activation reduce AP responses to visual stimulation. These changes were paralleled by restructuring of the neurons response to surround stimulus locations and membrane-potential time-course. Computational simulations could reproduce the observed subthreshold voltage changes only when presynaptic temporal jitter was included. Together this shows that STDP rules can modify output patterns of sensory neurons and the timing of single-APs plays a crucial role in sensory coding and plasticity.DOI:http://dx.doi.org/10.7554/eLife.00012.001