Dynamic Patterns of Spontaneous Ongoing Activity in the Visual Cortex of Anesthetized and Awake Monkeys are Different.

Ongoing internal cortical activity plays a major role in perception and behavior both in animals and humans. Previously we have shown that spontaneous patterns resembling orientation-maps appear over large cortical areas in the primary visual-cortex of anesthetized cats. However, it remains unknown 1) whether spontaneous-activity in the primate also displays similar patterns and 2) whether a significant difference exists between cortical ongoing-activity in the anesthetized and awake primate. We explored these questions by combining voltage-sensitive-dye imaging with multiunit and local-field-potential recordings. Spontaneously emerging orientation and ocular-dominance maps, spanning up to 6 × 6 mm2, were readily observed in anesthetized but not in awake monkeys. Nevertheless, spontaneous correlated-activity involving orientation-domains was observed in awake monkeys. Under both anesthetized and awake conditions, spontaneous correlated-activity coincided with traveling waves. We found that spontaneous activity resembling orientation-maps in awake animals spans smaller cortical areas in each instance, but over time it appears across all of V1. Furthermore, in the awake monkey, our results suggest that the synaptic strength had been completely reorganized including connections between dissimilar elements of the functional architecture. These findings lend support to the notion that ongoing-activity has many more fast switching representations playing an important role in cortical function and behavior.

Interhemispheric Synchrony of Spontaneous Cortical States at the Cortical Column Level.

In cat early visual cortex, neural activity patterns resembling evoked orientation maps emerge spontaneously under anesthesia. To test if such patterns are synchronized between hemispheres, we performed bilateral imaging in anesthetized cats using a new improved voltage-sensitive dye. We observed map-like activity patterns spanning early visual cortex in both hemispheres simultaneously. Patterns virtually identical to maps associated with the cardinal and oblique orientations emerged as leading principal components of the spontaneous fluctuations, and the strength of transient orientation states was correlated with their duration, providing evidence that these maps are transiently attracting states. A neural mass model we developed reproduced the dynamics of both smooth and abrupt orientation state transitions observed experimentally. The model suggests that map-like activity arises from slow modulations in spontaneous firing in conjunction with interplay between excitation and inhibition. Our results highlight the efficiency and functional precision of interhemispheric connectivity.

Critical dynamics, anesthesia and information integration: Lessons from multi-scale criticality analysis of voltage imaging data.

Critical dynamics are thought to play an important role in neuronal information-processing: near critical networks exhibit neuronal avalanches, cascades of spatiotemporal activity that are scale-free, and are considered to enhance information capacity and transfer. However, the exact relationship between criticality, awareness, and information integration remains unclear. To characterize this relationship, we applied multi-scale avalanche analysis to voltage-sensitive dye imaging data collected from animals of various species under different anesthetics. We found that anesthesia systematically varied the scaling behavior of neural dynamics, a change that was mirrored in reduced neural complexity. These findings were corroborated by applying the same analyses to a biophysically realistic cortical network model, in which multi-scale criticality measures were associated with network properties and the capacity for information integration. Our results imply that multi-scale criticality measures are potential biomarkers for assessing the level of consciousness.

Milliseconds of Sensory Input Abruptly Modulate the Dynamics of Cortical States for Seconds.

Spontaneous internal activity plays a major role in higher brain functions. The question of how it modulates sensory evoked activity and behavior has been explored in anesthetized rodents, cats, monkeys and in behaving human subjects. However, the complementary question of how a brief sensory input modulates the internally generated activity in vivo remains unresolved, and high-resolution mapping of these bidirectional interactions was never performed. Integrating complementary methodologies, at population and single cells levels, we explored this question. Voltage-sensitive dye imaging of population activity in anesthetized rats' somatosensory cortex revealed that spontaneous up-states were largely diminished for ~2 s, even after a single weak whisker deflection. This effect was maximal at the stimulated barrel but spread across several cortical areas. A higher velocity whisker deflection evoked activity at ~15Hz. Two-photon calcium imaging activity and cell-attached recordings confirmed the VSD results and revealed that for several seconds most single cells decreased their firing, but a small number increased firing. Comparing single deflection with long train stimulation, we found a dominant effect of the first population spike. We suggest that, at the onset of a sensory input, some internal messages are silenced to prevent overloading of the processing of relevant incoming sensory information.

Imaging the Dynamics of Neocortical Population Activity in Behaving and Freely Moving Mammals.

The development of functional imaging techniques applicable to neuroscience and covering a wide range of spatial and temporal scales has greatly facilitated the exploration of the relationships between cognition, behaviour and electrical brain activity. For mammals, the neocortex plays a particularly profound role in generating sensory perception, controlling voluntary movement, higher cognitive functions and planning goal-directed behaviours. Since these remarkable functions of the neocortex cannot be explored in simple model preparations or in anesthetised animals, the neural basis of behaviour must be explored in awake behaving subjects. Because neocortical function is highly distributed across many rapidly interacting regions, it is essential to measure spatiotemporal dynamics of cortical activity in real-time. Extensive work in anesthetised mammals has shown that in vivo Voltage-Sensitive Dye Imaging (VSDI) reveals the neocortical population membrane potential dynamics at millisecond temporal resolution and subcolumnar spatial resolution. Here, we describe recent advances indicating that VSDI is also already well-developed for exploring cortical function in behaving monkeys and mice. The first animal model, the non-human primate, is well-suited for fundamental exploration of higher-level cognitive function and behavior. The second animal model, the mouse, benefits from a rich arsenal of molecular and genetic technologies. In the monkey, imaging from the same patch of cortex, repeatedly, is feasible for a long period of time, up to a year. In the rodent, VSDI is applicable to freely moving and awake head-restrained mice. Interactions between different cortical areas and different cortical columns can therefore now be dynamically mapped through VSDI and related to the corresponding behaviour. Thus by applying VSDI to mice and monkeys one can begin to explore how behaviour emerges from neuronal activity in neuronal networks residing in different cortical areas.

Imaging the Dynamics of Mammalian Neocortical Population Activity In-Vivo.

Neural computations underlying sensory perception, cognition and motor control are performed by populations of neurons at different anatomical and temporal scales. Few techniques are currently available for exploring dynamics of local and large range populations. Voltage-sensitive dye imaging (VSDI) reveals neural population activity in areas ranging from a few tens of microns to a couple of centimeters, or two areas up to ~10 cm apart. VSDI provides a sub-millisecond temporal resolution, and a spatial resolution of about 50 µm. The dye signal emphasizes subthreshold synaptic potentials. VSDI has been applied in the mouse, rat, gerbil, ferret, tree shrew, cat and monkey cortices, in order to explore lateral spread of retinotopic or somatotopic activation, the dynamic spatiotemporal pattern resulting from sensory activation, including the somatosensory, olfactory, auditory, and visual modalities, as well as motor preparation and the properties of spontaneously-occurring population activity. In this chapter we focus on VSDI in-vivo and review results obtained mostly in the visual system in our laboratory.

Imaging Submillisecond Membrane Potential Changes from Individual Regions of Single Axons, Dendrites and Spines.

A central question in neuronal network analysis is how the interaction between individual neurons produces behavior and behavioral modifications. This task depends critically on how exactly signals are integrated by individual nerve cells functioning as complex operational units. Regional electrical properties of branching neuronal processes which determine the input-output function of any neuron are extraordinarily complex, dynamic, and, in the general case, impossible to predict in the absence of detailed measurements. To obtain such a measurement one would, ideally, like to be able to monitor, at multiple sites, subthreshold events as they travel from the sites of origin (synaptic contacts on distal dendrites) and summate at particular locations to influence action potential initiation. It became possible recently to carry out this type of measurement using high-resolution multisite recording of membrane potential changes with intracellular voltage-sensitive dyes. This chapter reviews the development and foundation of the method of voltage-sensitive dye recording from individual neurons. Presently, this approach allows monitoring membrane potential transients from all parts of the dendritic tree as well as from axon collaterals and individual dendritic spines.

Temporally-structured acquisition of multidimensional optical imaging data facilitates visualization of elusive cortical representations in the behaving monkey.

Fundamental understanding of higher cognitive functions can greatly benefit from imaging of cortical activity with high spatiotemporal resolution in the behaving non-human primate. To achieve rapid imaging of high-resolution dynamics of cortical representations of spontaneous and evoked activity, we designed a novel data acquisition protocol for sensory stimulation by rapidly interleaving multiple stimuli in continuous sessions of optical imaging with voltage-sensitive dyes. We also tested a new algorithm for the "temporally structured component analysis" (TSCA) of a multidimensional time series that was developed for our new data acquisition protocol, but was tested only on simulated data (Blumenfeld, 2010). In addition to the raw data, the algorithm incorporates prior knowledge about the temporal structure of the data as well as input from other information. Here we showed that TSCA can successfully separate functional signal components from other signals referred to as noise. Imaging of responses to multiple visual stimuli, utilizing voltage-sensitive dyes, was performed on the visual cortex of awake monkeys. Multiple cortical representations, including orientation and ocular dominance maps as well as the hitherto elusive retinotopic representation of orientation stimuli, were extracted in only 10s of imaging, approximately two orders of magnitude faster than accomplished by conventional methods. Since the approach is rather general, other imaging techniques may also benefit from the same stimulation protocol. This methodology can thus facilitate rapid optical imaging explorations in monkeys, rodents and other species with a versatility and speed that were not feasible before.

A processing work-flow for measuring erythrocytes velocity in extended vascular networks from wide field high-resolution optical imaging data.

Comprehensive information on the spatio-temporal dynamics of the vascular response is needed to underpin the signals used in hemodynamics-based functional imaging. It has recently been shown that red blood cells (RBCs) velocity and its changes can be extracted from wide-field optical imaging recordings of intrinsic absorption changes in cortex. Here, we describe a complete processing work-flow for reliable RBC velocity estimation in cortical networks. Several pre-processing steps are implemented: image co-registration, necessary to correct for small movements of the vasculature, semi-automatic image segmentation for fast and reproducible vessel selection, reconstruction of RBC trajectories patterns for each micro-vessel, and spatio-temporal filtering to enhance the desired data characteristics. The main analysis step is composed of two robust algorithms for estimating the RBCs' velocity field. Vessel diameter and its changes are also estimated, as well as local changes in backscattered light intensity. This full processing chain is implemented with a software suite that is freely distributed. The software uses efficient data management for handling the very large data sets obtained with in vivo optical imaging. It offers a complete and user-friendly graphical user interface with visualization tools for displaying and exploring data and results. A full data simulation framework is also provided in order to optimize the performances of the algorithm with respect to several characteristics of the data. We illustrate the performance of our method in three different cases of in vivo data. We first document the massive RBC speed response evoked by a spreading depression in anesthetized rat somato-sensory cortex. Second, we show the velocity response elicited by a visual stimulation in anesthetized cat visual cortex. Finally, we report, for the first time, visually-evoked RBC speed responses in an extended vascular network in awake monkey extrastriate cortex.

Embedding of cortical representations by the superficial patch system.

Pyramidal cells in layers 2 and 3 of the neocortex of many species collectively form a clustered system of lateral axonal projections (the superficial patch system--Lund JS, Angelucci A, Bressloff PC. 2003. Anatomical substrates for functional columns in macaque monkey primary visual cortex. Cereb Cortex. 13:15-24. or daisy architecture--Douglas RJ, Martin KAC. 2004. Neuronal circuits of the neocortex. Annu Rev Neurosci. 27:419-451.), but the function performed by this general feature of the cortical architecture remains obscure. By comparing the spatial configuration of labeled patches with the configuration of responses to drifting grating stimuli, we found the spatial organizations both of the patch system and of the cortical response to be highly conserved between cat and monkey primary visual cortex. More importantly, the configuration of the superficial patch system is directly reflected in the arrangement of function across monkey primary visual cortex. Our results indicate a close relationship between the structure of the superficial patch system and cortical responses encoding a single value across the surface of visual cortex (self-consistent states). This relationship is consistent with the spontaneous emergence of orientation response-like activity patterns during ongoing cortical activity (Kenet T, Bibitchkov D, Tsodyks M, Grinvald A, Arieli A. 2003. Spontaneously emerging cortical representations of visual attributes. Nature. 425:954-956.). We conclude that the superficial patch system is the physical encoding of self-consistent cortical states, and that a set of concurrently labeled patches participate in a network of mutually consistent representations of cortical input.

Increased retinal blood flow velocity in patients with early diabetes mellitus.

PURPOSE: To compare retinal blood flow velocity in small vessels of patients with early diabetes mellitus (DM), without any morphologic changes related to diabetic retinopathy, with that in a control group. METHODS: The authors used the retinal function imager to measure blood flow velocities, from many small vessels, simultaneously. Twenty-three eyes of 14 patients with early DM and 51 eyes of 31 healthy subjects were enrolled. Differences between the patients and the control group were assessed by mixed linear models. RESULTS: Venous average velocity significantly increased in the DM group (3.8 ± 1.2 vs. 2.9 ± 0.5 mm/second, P < 0.0001) than in the healthy subjects. Arterial velocity of DM patients was also significantly higher (4.7 ± 1.7 vs. 4.1 ± 0.9 mm/second, P = 0.03). There was no statistically significant difference between groups in age, gender, heart rate, and systolic blood pressure. The diastolic blood pressure in the DM patients was lower than that in the healthy group (P = 0.03). CONCLUSION: There was an increase in arterial and venous retinal blood flow velocities of patients with early DM with no diabetic retinopathy. These findings support the notion that abnormalities in vessel function exist in diabetic eyes before the development of structural changes. This noninvasive approach facilitated the assessment of early hemodynamic abnormalities and may assist in screening and monitoring.

A dynamic neural field model of mesoscopic cortical activity captured with voltage-sensitive dye imaging.

A neural field model is presented that captures the essential non-linear characteristics of activity dynamics across several millimeters of visual cortex in response to local flashed and moving stimuli. We account for physiological data obtained by voltage-sensitive dye (VSD) imaging which reports mesoscopic population activity at high spatio-temporal resolution. Stimulation included a single flashed square, a single flashed bar, the line-motion paradigm--for which psychophysical studies showed that flashing a square briefly before a bar produces sensation of illusory motion within the bar--and moving squares controls. We consider a two-layer neural field (NF) model describing an excitatory and an inhibitory layer of neurons as a coupled system of non-linear integro-differential equations. Under the assumption that the aggregated activity of both layers is reflected by VSD imaging, our phenomenological model quantitatively accounts for the observed spatio-temporal activity patterns. Moreover, the model generalizes to novel similar stimuli as it matches activity evoked by moving squares of different speeds. Our results indicate that feedback from higher brain areas is not required to produce motion patterns in the case of the illusory line-motion paradigm. Physiological interpretation of the model suggests that a considerable fraction of the VSD signal may be due to inhibitory activity, supporting the notion that balanced intra-layer cortical interactions between inhibitory and excitatory populations play a major role in shaping dynamic stimulus representations in the early visual cortex.

Retinal blood flow velocity measured by retinal function imaging in retinitis pigmentosa.

BACKGROUND: To measure the retinal blood flow velocity in patients with retinitis pigmentosa using the retinal function imaging technique. METHODS: The clinical observational investigation included a study group of five eyes of five patients (age: 55.7 ± 8.6 years) with retinitis pigmentosa (RP) and a control group of five eyes of five healthy subjects. We used a randomly chosen eye of the RP patients, and compared its results to the normal subjects using a mixed linear model, correcting for heart rate, age, and gender. RESULTS: The mean blood velocity in the narrow retinal veins (1.7 ± 0.35 cm/s versus 3.0 ± 0.35 cm/s; P < 0.001) and wide retinal veins (1.5 ± 0.35 cm/s versus 3.1 ± 0.30 cm/s; P < 0.001) was significantly lower in the study group than in the control group not correcting for heart rate, age or gender. Correspondingly, the arterial blood flow velocity was significantly lower in the study group than in the control group for the narrow arterial vessels (2.3 ± 0.55 versus 4.2 ± 0.5; P = 0.006) and for the wide retinal arteries (2.5 ± 1.05 cm/s versus 4.8 ± 1.0 cm/s; P < 0.001). CONCLUSIONS: Using the retinal function imaging technology revealed significantly lower retinal blood flow velocities in the small and large retinal vessels in patients with retinitis pigmentosa than in healthy subjects. This corresponds with the known decrease in the retinal vessel diameters as observed upon ophthalmoscopy in patients with retinitis pigmentosa. Retinal function imaging technology may hold promise for measurements of retinal blood flow parameters.

Imaging cortical correlates of illusion in early visual cortex.

Exploring visual illusions reveals fundamental principles of cortical processing. Illusory motion perception of non-moving stimuli was described almost a century ago by Gestalt psychologists. However, the underlying neuronal mechanisms remain unknown. To explore cortical mechanisms underlying the 'line-motion' illusion, we used real-time optical imaging, which is highly sensitive to subthreshold activity. We examined, in the visual cortex of the anaesthetized cat, responses to five stimuli: a stationary small square and a long bar; a moving square; a drawn-out bar; and the well-known line-motion illusion, a stationary square briefly preceding a long stationary bar presentation. Whereas flashing the bar alone evoked the expected localized, short latency and high amplitude activity patterns, presenting a square 60-100 ms before a bar induced the dynamic activity patterns resembling that of fast movement. The preceding square, even though physically non-moving, created gradually propagating subthreshold cortical activity that must contribute to illusory motion, because it was indistinguishable from cortical representations of real motion in this area. These findings demonstrate the effect of spatio-temporal patterns of subthreshold synaptic potentials on cortical processing and the shaping of perception.

Spontaneously emerging cortical representations of visual attributes.

Spontaneous cortical activity--ongoing activity in the absence of intentional sensory input--has been studied extensively, using methods ranging from EEG (electroencephalography), through voltage sensitive dye imaging, down to recordings from single neurons. Ongoing cortical activity has been shown to play a critical role in development, and must also be essential for processing sensory perception, because it modulates stimulus-evoked activity, and is correlated with behaviour. Yet its role in the processing of external information and its relationship to internal representations of sensory attributes remains unknown. Using voltage sensitive dye imaging, we previously established a close link between ongoing activity in the visual cortex of anaesthetized cats and the spontaneous firing of a single neuron. Here we report that such activity encompasses a set of dynamically switching cortical states, many of which correspond closely to orientation maps. When such an orientation state emerged spontaneously, it spanned several hypercolumns and was often followed by a state corresponding to a proximal orientation. We suggest that dynamically switching cortical states could represent the brain's internal context, and therefore reflect or influence memory, perception and behaviour.

Reduced retinal blood flow velocity in diabetic retinopathy.

PURPOSE: The purpose of this study was to compare the retinal blood flow velocities of patients with diabetes and healthy control subjects. We used a novel device offering a noninvasive diagnostic of retinal function. METHODS: Flow velocities in retinal arterioles and venules were quantitatively analyzed by retinal function imager scanning in 58 eyes of 42 patients with nonproliferative diabetic retinopathy and 51 eyes of 32 normal subjects. Group differences were assessed by the mixed-model effect. RESULTS: Average velocity in arterial compartments (in mm/s) was 3.74 +/- 1.09 for the diabetic group and 4.19 +/- 0.99 for the control subjects. The average velocity of all segments, taking associated heart rate and individual segment widths into account, was 17% slower in the diabetic group (P < 0.0001). In both groups, average venous compartment velocity was lower than the arterial velocity (2.61 +/- 0.65 for the diabetic group; 3.03 +/- 0.59 for the control subjects). Individual vein velocities, taking heart rate and segment widths into account, was 17% slower, on average, in the diabetic group (P < 0.0001). CONCLUSION: Our measurement showed significantly decreased flow velocities in the retinal arterioles and venules of patients with diabetes compared with healthy control subjects, supporting the view of abnormal vessel function in eyes with nonproliferative diabetic retinopathy.

Arousal increases the representational capacity of cortical tissue.

Arousal patently transforms the faculties of complex organisms. Although typical changes in cortical activity such as seen in EEG and LFP measurements are associated with change in state of arousal, it remains unclear what in the constitution of such state dependent activity enables this profound enhancement of ability. We put forward the hypothesis that arousal modulates cortical activity by rendering it more fit to represent information. We argue that representational capacity is of a dual nature-it requires not only that cortical tissue generate complex activity (i.e. spatiotemporal neuronal events), but also a complex cortical activity space (which is comprised of such spatiotemporal events). We explain that the topological notion of complexity-homology-is the pertinent measure of the complexity of neuronal activity spaces, as homological structure indicates not only the degree to which underlying activity is inherently clustered but also registers the effective dimensionality of the configurations formed by such clusters. Changes of this sort in the structure of cortical activity spaces can serve as the basis of the enhanced capacity to make perceptual/behavioral distinctions brought about by arousal. To show the feasibility of these ideas, we analyzed voltage sensitive dye imaging (VSDI) data acquired from primate visual cortex in disparate states of arousal. Our results lend some support to the theory: first as arousal increased so did the complexity of activity (that is the complexity of VSDI movies). Moreover, the complexity of structure of activity space (that is VSDI movie space) as measured by persistent homology-a multi scale topological measure of complexity-increased with arousal as well.

Removal of spatial biological artifacts in functional maps by local similarity minimization.

Functional maps obtained by various technologies, including optical imaging techniques, f-MRI, PET, and others, may be contaminated with biological artifacts such as vascular patterns or large patches of parenchyma. These artifacts originate mostly from changes in the microcirculation that result from either activity-dependent changes in volume or from oximetric changes that do not co-localize with neuronal activity per se. Standard methods do not always suffice to reduce such artifacts, in which case conspicuous spatial artifacts mask details of the underlying activity patterns. Here we propose a simple algorithm that efficiently removes spatial biological artifacts contaminating high-resolution functional maps. We validated this procedure by applying it to cortical maps resulting from optical imaging, based either on voltage-sensitive dye signals or on intrinsic signals. To remove vascular spatial patterns we first constructed a template of typical artifacts (vascular/cardiac pulsation/vasomotion), using principle components derived from baseline information obtained in the absence of stimulation. Next, we modified this template by means of local similarity minimization (LSM), achieved by measuring neighborhood similarity between contaminated data and the artifact template and then abolishing the similarity. LSM thus removed spatial patterns originating from the cortical vasculature components, including large fields of capillary parenchyma, helping to unveil details of neuronal activity patterns that were otherwise masked by these vascular artifacts. Examples obtained from our imaging experiments with anaesthetized cats and behaving monkeys showed that the LSM method is both general and reproducible, and is often superior to other available procedures.

Spatio-temporal dynamics of odor representations in the mammalian olfactory bulb.

We explored the spatio-temporal dynamics of odor-evoked activity in the rat and mouse main olfactory bulb (MOB) using voltage-sensitive dye imaging (VSDI) with a new probe. The high temporal resolution of VSDI revealed odor-specific sequences of glomerular activation. Increasing odor concentrations reduced response latencies, increased response amplitudes, and recruited new glomerular units. However, the sequence of glomerular activation was maintained. Furthermore, we found distributed MOB activity locked to the nasal respiration cycle. The spatial distribution of its amplitude and phase was heterogeneous and changed by sensory input in an odor-specific manner. Our data show that in the mammalian olfactory bulb, odor identity and concentration are represented by spatio-temporal patterns, rather than spatial patterns alone.

Columnar resolution of blood volume and oximetry functional maps in the behaving monkey; implications for FMRI.

The ultimate goal of high-resolution functional brain mapping is single-condition (stimulus versus no-stimulus maps) rather than differential imaging (comparing two "stimulus maps"), because the appropriate ("orthogonal") stimuli are rarely available. This requires some component(s) of activity-dependent hemodynamic signals to closely colocalize with electrical activity, like the early increase in deoxyhemoglobin, shown previously to yield high-quality functional single-condition maps. Conversely, nonlocal vascular responses dominate in cerebral blood volume (CBV)-based single-condition maps. Differential CBV maps are largely restricted to the parenchyma, implying that part of the CBV response does colocalize with electrical activity at fine spatial scale. By removing surface vascular activation from optical imaging data, we document the existence of a capillary CBV response component, regulated at fine spatial scale and yielding single-condition maps exhibiting approximately 100 microm resolution. Blood volume and -flow based single-condition functional mapping at columnar level should thus be feasible, provided that the capillary response component is selectively imaged.