On the Functional Role of Gamma Synchronization in the Retinogeniculate System of the Cat.

Fast gamma oscillations, generated within the retina, and transmitted to the cortex via the lateral geniculate nucleus (LGN), are thought to carry information about stimulus size and continuity. This hypothesis relies mainly on studies conducted under anesthesia and the extent to which it holds under more naturalistic conditions remains unclear. Using multielectrode recordings of spiking activity in the retina and the LGN of both male and female cats, we show that visually driven gamma oscillations are absent for awake states and are highly dependent on halothane (or isoflurane). Under ketamine, responses were nonoscillatory, as in the awake condition. Response entrainment to the monitor refresh was commonly observed up to 120 Hz and was superseded by the gamma oscillatory responses induced by halothane. Given that retinal gamma oscillations are contingent on halothane anesthesia and absent in the awake cat, such oscillations should be considered artifactual, thus playing no functional role in vision.SIGNIFICANCE STATEMENT Gamma rhythms have been proposed to be a robust encoding mechanism critical for visual processing. In the retinogeniculate system of the cat, many studies have shown gamma oscillations associated with responses to static stimuli. Here, we extend these observations to dynamic stimuli. An unexpected finding was that retinal gamma responses strongly depend on halothane concentration levels and are absent in the awake cat. These results weaken the notion that gamma in the retina is relevant for vision. Notably, retinal gamma shares many of the properties of cortical gamma. In this respect, oscillations induced by halothane in the retina may serve as a valuable preparation, although artificial, for studying oscillatory dynamics.

Cortical maps as a fundamental neural substrate for visual representation.

Visual perception is the product of serial hierarchical processing, parallel processing, and remapping on a dynamic network involving several topographically organized cortical visual areas. Here, we will focus on the topographical organization of cortical areas and the different kinds of visual maps found in the primate brain. We will interpret our findings in light of a broader representational framework for perception. Based on neurophysiological data, our results do not support the notion that vision can be explained by a strict representational model, where the objective visual world is faithfully represented in our brain. On the contrary, we find strong evidence that vision is an active and constructive process from the very initial stages taking place in the eye and from the very initial stages of our development. A constructive interplay between perceptual and motor systems (e.g., during saccadic eye movements) is actively learnt from early infancy and ultimately provides our fluid stable visual perception of the world.

A new robust multivariate mode estimator for eye-tracking calibration.

We propose in this work a new method for estimating the main mode of multivariate distributions, with application to eye-tracking calibration. When performing eye-tracking experiments with poorly cooperative subjects, such as infants or monkeys, the calibration data generally suffer from high contamination. Outliers are typically organized in clusters, corresponding to fixations in the time intervals when subjects were not looking at the calibration points. In this type of multimodal distributions, most central tendency measures fail at estimating the principal fixation coordinates (the first mode), resulting in errors and inaccuracies when mapping the gaze to the screen coordinates. Here, we developed a new algorithm to identify the first mode of multivariate distributions, named BRIL, which relies on recursive depth-based filtering. This novel approach was tested on artificial mixtures of Gaussian and Uniform distributions, and compared to existing methods (conventional depth medians, robust estimators of location and scatter, and clustering-based approaches). We obtained outstanding performances, even for distributions containing very high proportions of outliers, both grouped in clusters and randomly distributed. Finally, we demonstrate the strength of our method in a real-world scenario using experimental data from eye-tracking calibrations with Capuchin monkeys, especially for highly contaminated distributions where other algorithms typically lack accuracy.

Spontaneous variability in gamma dynamics described by a damped harmonic oscillator driven by noise.

Circuits of excitatory and inhibitory neurons generate gamma-rhythmic activity (30-80 Hz). Gamma-cycles show spontaneous variability in amplitude and duration. To investigate the mechanisms underlying this variability, we recorded local-field-potentials (LFPs) and spikes from awake macaque V1. We developed a noise-robust method to detect gamma-cycle amplitudes and durations, which showed a weak but positive correlation. This correlation, and the joint amplitude-duration distribution, is well reproduced by a noise-driven damped harmonic oscillator. This model accurately fits LFP power-spectra, is equivalent to a linear, noise-driven E-I circuit, and recapitulates two additional features of gamma: (1) Amplitude-duration correlations decrease with oscillation strength; (2) amplitudes and durations exhibit strong and weak autocorrelations, respectively, depending on oscillation strength. Finally, longer gamma-cycles are associated with stronger spike-synchrony, but lower spike-rates in both (putative) excitatory and inhibitory neurons. In sum, V1 gamma-dynamics are well described by the simplest possible model of gamma: A damped harmonic oscillator driven by noise.

A Distinct Class of Bursting Neurons with Strong Gamma Synchronization and Stimulus Selectivity in Monkey V1.

Cortical computation depends on interactions between excitatory and inhibitory neurons. The contributions of distinct neuron types to sensory processing and network synchronization in primate visual cortex remain largely undetermined. We show that in awake monkey V1, there exists a distinct cell type (>>30% of neurons) that has narrow-waveform (NW) action potentials and high spontaneous discharge rates and fires in high-frequency bursts. These neurons are more stimulus selective and phase locked to 30- to 80-Hz gamma oscillations than other neuron types. Unlike other neuron types, their gamma-phase locking is highly predictive of orientation tuning. We find evidence for strong rhythmic inhibition in these neurons, suggesting that they interact with interneurons to act as excitatory pacemakers for the V1 gamma rhythm. We did not find a similar class of NW bursting neurons in L2-L4 of mouse V1. Given its properties, this class of NW bursting neurons should be pivotal for the encoding and transmission of stimulus information.

Callosal Influence on Visual Receptive Fields Has an Ocular, an Orientation-and Direction Bias.

One leading hypothesis on the nature of visual callosal connections (CC) is that they replicate features of intrahemispheric lateral connections. However, CC act also in the central part of the binocular visual field. In agreement, early experiments in cats indicated that they provide the ipsilateral eye part of binocular receptive fields (RFs) at the vertical midline (Berlucchi and Rizzolatti, 1968), and play a key role in stereoscopic function. But until today callosal inputs to receptive fields activated by one or both eyes were never compared simultaneously, because callosal function has been often studied by cutting or lesioning either corpus callosum or optic chiasm not allowing such a comparison. To investigate the functional contribution of CC in the intact cat visual system we recorded both monocular and binocular neuronal spiking responses and receptive fields in the 17/18 transition zone during reversible deactivation of the contralateral hemisphere. Unexpectedly from many of the previous reports, we observe no change in ocular dominance during CC deactivation. Throughout the transition zone, a majority of RFs shrink, but several also increase in size. RFs are significantly more affected for ipsi- as opposed to contralateral stimulation, but changes are also observed with binocular stimulation. Noteworthy, RF shrinkages are tiny and not correlated to the profound decreases of monocular and binocular firing rates. They depend more on orientation and direction preference than on eccentricity or ocular dominance of the receiving neuron's RF. Our findings confirm that in binocularly viewing mammals, binocular RFs near the midline are constructed via the direct geniculo-cortical pathway. They also support the idea that input from the two eyes complement each other through CC: Rather than linking parts of RFs separated by the vertical meridian, CC convey a modulatory influence, reflecting the feature selectivity of lateral circuits, with a strong cardinal bias.

Context matters: the illusive simplicity of macaque V1 receptive fields.

Even in V1, where neurons have well characterized classical receptive fields (CRFs), it has been difficult to deduce which features of natural scenes stimuli they actually respond to. Forward models based upon CRF stimuli have had limited success in predicting the response of V1 neurons to natural scenes. As natural scenes exhibit complex spatial and temporal correlations, this could be due to surround effects that modulate the sensitivity of the CRF. Here, instead of attempting a forward model, we quantify the importance of the natural scenes surround for awake macaque monkeys by modeling it non-parametrically. We also quantify the influence of two forms of trial to trial variability. The first is related to the neuron's own spike history. The second is related to ongoing mean field population activity reflected by the local field potential (LFP). We find that the surround produces strong temporal modulations in the firing rate that can be both suppressive and facilitative. Further, the LFP is found to induce a precise timing in spikes, which tend to be temporally localized on sharp LFP transients in the gamma frequency range. Using the pseudo R(2) as a measure of model fit, we find that during natural scene viewing the CRF dominates, accounting for 60% of the fit, but that taken collectively the surround, spike history and LFP are almost as important, accounting for 40%. However, overall only a small proportion of V1 spiking statistics could be explained (R(2)∼5%), even when the full stimulus, spike history and LFP were taken into account. This suggests that under natural scene conditions, the dominant influence on V1 neurons is not the stimulus, nor the mean field dynamics of the LFP, but the complex, incoherent dynamics of the network in which neurons are embedded.

Mapping of visual receptive fields by tomographic reconstruction.

The moving bar experiment is a classic paradigm for characterizing the receptive field (RF) properties of neurons in primary visual cortex (V1). Current approaches for analyzing neural spiking activity recorded from these experiments do not take into account the point-process nature of these data and the circular geometry of the stimulus presentation. We present a novel analysis approach to mapping V1 receptive fields that combines point-process generalized linear models (PPGLM) with tomographic reconstruction computed by filtered-back projection. We use the method to map the RF sizes and orientations of 251 V1 neurons recorded from two macaque monkeys during a moving bar experiment. Our cross-validated goodness-of-fit analyses show that the PPGLM provides a more accurate characterization of spike train data than analyses based on rate functions computed by the methods of spike-triggered averages or first-order Wiener-Volterra kernel. Our analysis leads to a new definition of RF size as the spatial area over which the spiking activity is significantly greater than baseline activity. Our approach yields larger RF sizes and sharper orientation tuning estimates. The tomographic reconstruction paradigm further suggests an efficient approach to choosing the number of directions and the number of trials per direction in designing moving bar experiments. Our results demonstrate that standard tomographic principles for image reconstruction can be adapted to characterize V1 RFs and that two fundamental properties, size and orientation, may be substantially different from what is currently reported.

Orientation selectivity and noise correlation in awake monkey area V1 are modulated by the gamma cycle.

Gamma-band synchronization adjusts the timing of excitatory and inhibitory inputs to a neuron. Neurons in the visual cortex are selective for stimulus orientation because of dynamic interactions between excitatory and inhibitory inputs. We hypothesized that these interactions and hence also orientation selectivity vary during the gamma cycle. We determined for each spike its phase relative to the gamma cycle. As a function of gamma phase, we then determined spike rates and their orientation selectivity. Orientation selectivity was modulated by gamma phase. The firing rate of spiking activity that occurred close to a neuron's mean gamma phase of firing was most orientation selective. This stimulus-selective signal could best be conveyed to postsynaptic neurons if it were not corrupted by noise correlations. Noise correlations between firing rates were modulated by gamma phase such that they were not statistically detectable for the spiking activity occurring close to a neuron's mean gamma phase of firing. Thus, gamma-band synchronization produces spiking activity that carries maximal stimulus selectivity and minimal noise correlation in its firing rate, and at the same time synchronizes this spiking activity for maximal impact on postsynaptic targets.

Gamma responses correlate with temporal expectation in monkey primary visual cortex.

Gamma oscillations have been linked to selective attention. Here, we investigate the effects of expecting a behaviorally relevant event (a change in the fixation point) on the oscillatory patterning of the local field potential and spiking responses in V1. Three protocols were used. In the first protocol, fixation point change occurred at a fixed time point, enabling predictions on task timing. In the second, fixation point change occurred in trial blocks either early or late in the trial, allowing us to compare responses during epochs of low and high expectation. Finally, we used a cue to indicate the upcoming fixation point change. All protocols led to an increase in gamma oscillations associated with alpha suppression when the monkeys expected an event in time. These effects were spatially widespread, since comparable results were observed for both central and peripheral visual representations in V1. Our findings indicate that expectations associated with perceptual decisions, motor responses, or upcoming reward may have a strong effect on the primary visual cortex, causing global, spatially nonselective modulation of gamma activity.

Extraction of Network Topology From Multi-Electrode Recordings: Is there a Small-World Effect?

The simultaneous recording of the activity of many neurons poses challenges for multivariate data analysis. Here, we propose a general scheme of reconstruction of the functional network from spike train recordings. Effective, causal interactions are estimated by fitting generalized linear models on the neural responses, incorporating effects of the neurons' self-history, of input from other neurons in the recorded network and of modulation by an external stimulus. The coupling terms arising from synaptic input can be transformed by thresholding into a binary connectivity matrix which is directed. Each link between two neurons represents a causal influence from one neuron to the other, given the observation of all other neurons from the population. The resulting graph is analyzed with respect to small-world and scale-free properties using quantitative measures for directed networks. Such graph-theoretic analyses have been performed on many complex dynamic networks, including the connectivity structure between different brain areas. Only few studies have attempted to look at the structure of cortical neural networks on the level of individual neurons. Here, using multi-electrode recordings from the visual system of the awake monkey, we find that cortical networks lack scale-free behavior, but show a small, but significant small-world structure. Assuming a simple distance-dependent probabilistic wiring between neurons, we find that this connectivity structure can account for all of the networks' observed small-world ness. Moreover, for multi-electrode recordings the sampling of neurons is not uniform across the population. We show that the small-world-ness obtained by such a localized sub-sampling overestimates the strength of the true small-world structure of the network. This bias is likely to be present in all previous experiments based on multi-electrode recordings.

A new look at gamma? High- (>60 Hz) γ-band activity in cortical networks: function, mechanisms and impairment.

γ-band oscillations are thought to play a crucial role in information processing in cortical networks. In addition to oscillatory activity between 30 and 60 Hz, current evidence from electro- and magnetoencephalography (EEG/MEG) and local-field potentials (LFPs) has consistently shown oscillations >60 Hz (high γ-band) whose function and generating mechanisms are unclear. In the present paper, we summarize data that highlights the importance of high γ-band activity for cortical computations through establishing correlations between the modulation of oscillations in the 60-200 Hz frequency and specific cognitive functions. Moreover, we will suggest that high γ-band activity is impaired in neuropsychiatric disorders, such as schizophrenia and epilepsy. In the final part of the paper, we will review physiological mechanisms underlying the generation of high γ-band oscillations and discuss the functional implications of low vs. high γ-band activity patterns in cortical networks.

Gamma-phase shifting in awake monkey visual cortex.

Gamma-band synchronization is abundant in nervous systems. Typically, the strength or precision of gamma-band synchronization is studied. However, the precise phase with which individual neurons are synchronized to the gamma-band rhythm might have interesting consequences for their impact on further processing and for spike timing-dependent plasticity. Therefore, we investigated whether the spike times of individual neurons shift systematically in the gamma cycle as a function of the neuronal activation strength. We found that stronger neuronal activation leads to spikes earlier in the gamma cycle, i.e., we observed gamma-phase shifting. Gamma-phase shifting occurred on very rapid timescales. It was particularly pronounced for periods in which gamma-band synchronization was relatively weak and for neurons that were only weakly coupled to the gamma rhythm. We suggest that gamma-phase shifting is brought about by an interplay between overall excitation and gamma-rhythmic synaptic input and has interesting consequences for neuronal coding, competition, and plasticity.

Synchronization dynamics in response to plaid stimuli in monkey V1.

Gamma synchronization has generally been associated with grouping processes in the visual system. Here, we examine in monkey V1 whether gamma oscillations play a functional role in segmenting surfaces of plaid stimuli. Local field potentials (LFPs) and spiking activity were recorded simultaneously from multiple sites in the opercular and calcarine regions while the monkeys were presented with sequences of single and superimposed components of plaid stimuli. In accord with the previous studies, responses to the single components (gratings) exhibited strong and sustained gamma-band oscillations (30-65 Hz). The superposition of the second component, however, led to profound changes in the temporal structure of the responses, characterized by a drastic reduction of gamma oscillations in the spiking activity and systematic shifts to higher frequencies in the LFP ( approximately 10% increase). Comparisons between cerebral hemispheres and across monkeys revealed robust subject-specific spectral signatures. A possible interpretation of our results may be that single gratings induce strong cooperative interactions among populations of cells that share similar response properties, whereas plaids lead to competition. Overall, our results suggest that the functional architecture of the cortex is a major determinant of the neuronal synchronization dynamics in V1.

Neural synchrony in cortical networks: history, concept and current status.

Following the discovery of context-dependent synchronization of oscillatory neuronal responses in the visual system, the role of neural synchrony in cortical networks has been expanded to provide a general mechanism for the coordination of distributed neural activity patterns. In the current paper, we present an update of the status of this hypothesis through summarizing recent results from our laboratory that suggest important new insights regarding the mechanisms, function and relevance of this phenomenon. In the first part, we present recent results derived from animal experiments and mathematical simulations that provide novel explanations and mechanisms for zero and nero-zero phase lag synchronization. In the second part, we shall discuss the role of neural synchrony for expectancy during perceptual organization and its role in conscious experience. This will be followed by evidence that indicates that in addition to supporting conscious cognition, neural synchrony is abnormal in major brain disorders, such as schizophrenia and autism spectrum disorders. We conclude this paper with suggestions for further research as well as with critical issues that need to be addressed in future studies.

Directional responses of visual wulst neurones to grating and plaid patterns in the awake owl.

The avian retinothalamofugal pathway reaches the telencephalon in an area known as visual wulst. A close functional analogy between this area and the early visual cortex of mammals has been established in owls. The goal of the present study was to assess quantitatively the directional selectivity and motion integration capability of visual wulst neurones, aspects that have not been previously investigated. We recorded extracellularly from a total of 101 cells in awake burrowing owls. From this sample, 88% of the units exhibited modulated directional responses to sinusoidal gratings, with a mean direction index of 0.74 +/- 0.03 and tuning bandwidth of 28 +/- 1.16 degrees . A direction index higher than 0.5 was observed in 66% of the cells, thereby qualifying them as direction selective. Motion integration was tested with moving plaids, made by adding two sinusoidal gratings of different orientations. We found that 80% of direction-selective cells responded optimally to the motion direction of the component gratings, whereas none responded to the global motion of plaids, whose direction was intermediate to that of the gratings. The remaining 20% were unclassifiable. The strength of component motion selectivity rapidly increased over a 200 ms period following stimulus onset, maintaining a relatively sustained profile thereafter. Overall, our data suggest that, as in the mammalian primary visual cortex, the visual wulst neurones of owls signal the local orientated features of a moving object. How and where these potentially ambiguous signals are integrated in the owl brain might be important for understanding the mechanisms underlying global motion perception.

Brightness induction: rate enhancement and neuronal synchronization as complementary codes.

In cat visual cortex, we investigated with parallel recordings from multiple units the neuronal correlates of perceived brightness. The perceived brightness of a center grating was changed by varying the orientation or the relative spatial phase of a surrounding grating. Brightness enhancement by orientation contrast is associated with an increase of discharge rates of responses to the center grating but not with changes in spike synchronization. In contrast, if brightness enhancement is induced by phase offset, discharge rates are unchanged but synchronization increases between neurons responding to the center grating. The changes in synchronization correlate well with changes in perceived brightness that were assessed in parallel in human subjects using the same stimuli. These results indicate that in cerebral cortex the modulation of synchronicity of responses is used as a mechanism complementary to rate changes to enhance the saliency of neuronal responses.

See globally, spike locally: oscillations in a retinal model encode large visual features.

We show that coherent oscillations among neighboring ganglion cells in a retinal model encode global topological properties, such as size, that cannot be deduced unambiguously from their local, time-averaged firing rates. Whereas ganglion cells may fire similar numbers of spikes in response to both small and large spots, only large spots evoke coherent high frequency oscillations, potentially allowing downstream neurons to infer global stimulus properties from their local afferents. To determine whether such information might be extracted over physiologically realistic spatial and temporal scales, we analyzed artificial spike trains whose oscillatory correlations were similar to those measured experimentally. Oscillatory power in the upper gamma band, extracted on single-trials from multi-unit spike trains, supported good to excellent size discrimination between small and large spots, with performance improving as the number of cells and/or duration of the analysis window was increased. By using Poisson distributed spikes to normalize the firing rate across stimulus conditions, we further found that coincidence detection, or synchrony, yielded substantially poorer performance on identical size discrimination tasks. To determine whether size encoding depended on contiguity independent of object shape, we examined the total oscillatory activity across the entire model retina in response to random binary images. As the ON-pixel probability crossed the percolation threshold, which marks the sudden emergence of large connected clusters, the total gamma-band activity exhibited a sharp transition, a phenomena that may be experimentally observable. Finally, a reanalysis of previously published oscillatory responses from cat ganglion cells revealed size encoding consistent with that predicted by the retinal model.

Activity patterns in human motion-sensitive areas depend on the interpretation of global motion.

Numerous imaging studies have contributed to the localization of motion-sensitive areas in the human brain. It is, however, still unclear how these areas contribute to global motion perception. Here, we investigate with functional MRI whether the motion-sensitive area hMT+/V5 is involved in perceptual segmentation and integration of motion signals. Stimuli were overlapping moving gratings that can be perceived either as two independently moving, transparent surfaces or as a single surface moving in an intermediate direction. We examined whether motion-sensitive area hMT+/V5 is involved in mediating the switches between the two percepts. The data show differential activation of hMT+/V5 with perceptual switches, suggesting that these are associated with a reconfiguration of cell assemblies in this area.

Feed-forward synchronization: propagation of temporal patterns along the retinothalamocortical pathway.

Visual responses in the cortex and lateral geniculate nucleus (LGN) are often associated with synchronous oscillatory patterning. In this short review, we examine the possible relationships between subcortical and cortical synchronization mechanisms. Our results obtained from simultaneous multi-unit recordings show strong synchronization of oscillatory responses between retina, LGN and cortex, indicating that cortical neurons can be synchronized by oscillatory activity relayed through the LGN. This feed-forward synchronization mechanism operating in the 60 to 120 Hz frequency range was observed mostly for static stimuli. In response to moving stimuli, by contrast, cortical synchronization was independent of oscillatory inputs from the LGN, with oscillation frequency in the range of 30 to 60 Hz. The functional implications of synchronization of activity from parallel channels are discussed, in particular its significance for signal transmission and cortical integration processes.