Publications

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.

Keywords: Ophthalmology

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.

Purpose/Aim of the study: To study changes in retinal blood flow velocity in patients with early and neovascular age-related macular degeneration (AMD). We used the Retinal Function Imager (RFI, Optical Imaging Ltd., Rehovot, Israel), a noninvasive diagnostic approach for measuring blood flow velocity. Materials and Methods: Sixty eyes of 43 AMD patients and 53 eyes of 35 healthy individuals over the age of 50 were recruited for this study. All patients were scanned by the RFI with analysis of blood flow velocity of secondary and tertiary branches of arteries and veins. Differences among groups were assessed by mixed linear models. Results: The average velocity in AMD patients was significantly lower compared to controls in arteries (3.6 +/- 1.4 versus 4.3 +/- 1.0 mm/sec, p=0.009) but not in veins (2.6 +/- 0.9 versus 3.1 +/- 0.6 mm/sec, p = 0.08). When comparing the velocity between low-and high-grade AMD eyes, venous velocity was slower in the high grade AMD eyes only in the "narrow'' group of vessels. Conclusions: Decreased blood flow velocity in retinal arteries in patients with AMD was found. Despite the fact that AMD is essentially a choroidal disease, retinal vessels show a functional abnormality, which may suggest that the vascular abnormality in this disease is more generalized.

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 10 s 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. (C) 2013 Published by Elsevier Inc.

PURPOSE. To study the short-term effects of intravitreal bevacizumab (Avastin) on retinal blood flow velocity and compare them to clinical outcomes assessed by optical coherence tomography (OCT) and tests of visual acuity. METHODS. The Retinal Function Imager (RFI) was used noninvasively and quantitatively to measure retinal blood flow velocity. Eight patients receiving intravitreal injection of Avastin for choroidal neovascularization (CNV) were included in this study. All were imaged by the RFI preinjection and 1 and 7 days postinjection. Visual acuity (VA) and OCT were recorded preinjection and 1 month postinjection. Comparisons were performed using paired Student t test and correlation using Spearman rank test. RESULTS. A good correlation was found between the 1-month change in VA and OCT measurements and the short-term change induced in blood flow velocity. Arterial and venous velocity changes 1 day after the injection correlated with the VA change (p

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.

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.

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.

Little is known about the "inverse" of the receptive field-the region of cortical space whose spatiotemporal pattern of electrical activity is influenced by a given sensory stimulus. We refer to this activated area as the cortical response field, the properties of which remain unexplored. Here, the dynamics of cortical response fields evoked in visual cortex by small, local drifting-oriented gratings were explored using voltage-sensitive dyes. We found that the cortical response field was often characterized by a plateau of activity, beyond the rim of which activity diminished quickly. Plateau rim location was largely independent of stimulus orientation. However, approximately 20 ms following plateau onset, 1-3 peaks emerged on it and were amplified for 25 ms. Spiking was limited to the peak zones, whose location strongly depended on stimulus orientation. Thus, alongside selective amplification of a spatially restricted suprathreshold response, wider activation to just below threshold encompasses all orientation domains within a well-defined retinotopic vicinity of the current stimulus, priming the cortex for processing of subsequent stimuli.

During the last few decades, neuroscientists have benefited from the emergence of many powerful functional imaging techniques that cover broad spatial and temporal scales. We can now image single molecules controlling cell differentiation, growth and death; single cells and their neurites processing electrical inputs and sending outputs; neuronal circuits performing neural computations in vitro; and the intact brain. At present, imaging based on voltage-sensitive dyes (VSDI) offers the highest spatial and temporal resolution for imaging neocortical functions in the living brain, and has paved the way for a new era in the functional imaging of cortical dynamics. It has facilitated the exploration of fundamental mechanisms that underlie neocortical development, function and plasticity at the fundmental level of the cortical column.

The rodent primary somatosensory cortex is spontaneously active in the form of locally synchronous membrane depolarizations (UP states) separated by quiescent hyperpolarized periods (DOWN states) both under anesthesia and during quiet wakefulness. In vivo whole-cell recordings and tetrode unit recordings were combined with voltage-sensitive dye imaging to analyze the relationship of the activity of individual pyramidal neurons in layer 2/3 to the ensemble spatiotemporal dynamics of the spontaneous depolarizations. These were either brief and localized to an area of a barrel column or occurred as propagating waves dependent on local glutamatergic synaptic transmission in layer 2/3. Spontaneous activity inhibited the sensory responses evoked by whisker deflection, accounting almost entirely for the large trial-to-trial variability of sensory-evoked postsynaptic potentials and action potentials. Subthreshold sensory synaptic responses evoked while a cortical area was spontaneously depolarized were smaller, briefer and spatially more confined. Surprisingly, whisker deflections evoked fewer action potentials during the spontaneous depolarizations despite neurons being closer to threshold. The ongoing spontaneous activity thus regulates the amplitude and the time-dependent spread of the sensory response in layer 2/3 barrel cortex.

A novel method of chronic optical imaging based on new voltage-sensitive dyes (VSDs) was developed to facilitate the explorations of the spatial and temporal patterns underlying higher cognitive functions in the neocortex of behaving monkeys. Using this system, we were able to explore cortical dynamics, with high spatial and temporal resolution, over period of less than or equal to1 yr from the same patch of cortex. The visual cortices of trained macaques were stained one to three times a week, and immediately after each staining session, the monkey started to perform the behavioral task, while the primary and secondary visual areas (V1 and V2) were imaged with a fast optical imaging system. Long-term repeated VSD imaging (VSDI) from the same cortical area did not disrupt the normal cortical architecture as confirmed repeatedly by optical imaging based on intrinsic signals. The spatial patterns of functional maps obtained by VSDI were essentially identical to those obtained from the same patch of cortex by imaging based on intrinsic signals. On comparing the relative amplitudes of the evoked signals and differential map obtained using these two different imaging methodologies, we found that VSDI emphasizes subthreshold activity more than imaging based on intrinsic signals, that emphasized more spiking activity. The latency of the VSD-evoked response in V1 ranged from 46 to 68 ms in the different monkeys. The amplitude of the V2 response was only 20-60% of that in V1 As expected from the anatomy, the retinotopic responses to local visual stimuli spread laterally across the cortical surface at a spreading velocity of 0.15-0.19 m/s over a larger area than that expected by the classical magnification factor, reaching its maximal anisotropic spatial extent within similar to40 ms. We correlated the observed dynamics of cortical activation patterns with the monkey's saccadic eye movements and found that due to the slow offset of the cortical response relative to its ons

High-resolution images of the somatotopic hand representation in macaque monkey primary somatosensory cortex (area S-I) were obtained by optical imaging based on intrinsic signals. To visualize somatotopic maps, we imaged optical responses to mild tactile stimulation of each individual fingertip. The activation evoked by stimulation of a single finger was strongest in a narrow transverse band (similar to1 x 4 mm) across the postcentral gyrus. As expected, a sequential organization of these bands was found. However, a significant overlap, especially for the activated areas of fingers 3-5, was found. Surprisingly, in addition to the finger-specific domains, we found that for each of the fingers, weak stimulation activated also a second "common patch" of cortex, located just medially to the representation of the finger. These results were confirmed by targeted multiunit and single-unit recordings guided by the optical maps. The maps remained very stable over many hours of recording. By optimizing the imaging procedures, we were able to obtain the functional maps extremely rapidly (e.g., the map of five fingers in the macaque monkey could be obtained in as little as 5 min). Furthermore, we describe the intraoperative optical imaging of the hand representation in the human brain during neurosurgery and then discuss the implications of the present results for the spatial resolution accomplishable by other neuroimaging techniques, relying on responses of the microcirculation to sensory-evoked electrical activity. This study demonstrates the feasibility of using high-resolution optical imaging to explore reliably short- and long-term plasticity of cortical representations, as well as for applications in the clinical setting.

Explorations of learning and memory, other long-term plastic changes, and additional cognitive functions in the behaving primate brain would greatly benefit from the ability to image the functional architecture within the same patch of cortex, at the columnar level, for a long period of time. We developed methods for long-term optical imaging based on intrinsic signals and repeatedly visualized the same functional domains in behaving macaque cortex for a period extending over 1 year. Using optical imaging and imaging spectroscopy, we first explored the relationship between electrical activity and hemodynamic events in the awake behaving primate and compared it with anesthetized preparations. We found that, whereas the amplitude of the intrinsic signal was much larger in the awake animal, its temporal pattern was similar to that observed in the anesthetized animals. In both groups, deoxyhemoglobin concentration reached a peak 2-3 sec after stimulus onset. Furthermore, the early activity-dependent increase in deoxyhemoglobin concentration (the "initial dip") was far more tightly colocalized with electrical activity than the delayed increase in oxyhemoglobin concentration, known to be associated with an increase in blood flow. The implications of these results for improvement of the spatial resolution of blood oxygenation level-dependent functional magnetic resonance imaging are discussed. After the characterization of the intrinsic signal in the behaving primate, we used this new imaging method to explore the stability of cortical maps in the macaque primary visual cortex. Functional maps of orientation and ocular dominance columns were found to be stable for a period longer than 1 year.

The relation between the activity of a single neocortical neuron and the dynamics of the network in which it is embedded was explored by single-unit recordings and real-time optical imaging. The firing rate of a spontaneously active single neuron strongly depends on the instantaneous spatial pattern of ongoing population activity in a large cortical area. Very similar spatial patterns of population activity were observed both when the neuron fired spontaneously and when it was driven by its optimal stimulus. The evoked patterns could be used to reconstruct the spontaneous activity of single neurons.

Modern functional neuroimaging methods, such as positron-emission tomography (PET), optical imaging of intrinsic signals, and functional MRI (fMRI) utilize activity dependent hemodynamic changes to obtain indirect maps of the evoked electrical activity in the brain. Whereas PET and flow-sensitive MRI map cerebral blood flow (CBF) changes, optical imaging and blood oxygenation level-dependent MRI map areas with changes in the concentration of deoxygenated hemoglobin (HbR), However, the relationship between CBF and HbR during functional activation has never been tested experimentally, Therefore, we investigated this relationship by using imaging spectroscopy and laser-Doppler flowmetry techniques, simultaneously, in the visual cortex of anesthetized cats during sensory stimulation, We found that the earliest microcirculatory change was indeed an increase in HbR, whereas the CBF increase lagged by more than a second after the increase in HbR, The increased HbR was accompanied by a simultaneous increase in total hemoglobin concentration (Hbt), presumably reflecting an early blood volume increase. We found that the CBF changes lagged after Hbt changes by 1 to 2 sec throughout the response, These results support the notion of active neurovascular regulation of blood volume in the capillary bed and the existence of a delayed, passive process of capillary filling.

The goal of this study was to explore the functional organization of direction of motion in cat area 18. Optical imaging was used to record the activity of populations of neurons, We found a patchy distribution of cortical regions exhibiting preference for one direction over the opposite direction of motion, The degree of clustering according to preference of direction was two to four limes smaller than that observed for orientation. In general, direction preference changed smoothly along the cortical surface; however, discontinuities in the direction maps were observed. These discontinuities formed lines that separated pairs of patches with preference for opposite directions. The functional maps for direction and for orientation preference were closely related; typically, an iso-orientation patch was divided into regions that exhibited preference for opposite directions, orthogonal to the orientation. In addition, the lines of discontinuity within the direction map often connected points of singularity in the orientation map, Although the organization of both domains was related, the direction and the orientation selective responses were separable; whereas the selective response according to direction of motion was nearly independent of the length of bars used for visual stimulation, the selective response to orientation decreased significantly with decreasing length of the bars. Extensive single and multiunit electrical recordings, targeted to selected domains of the functional maps, confirmed the features revealed by optical imaging, We conclude that significant processing of direction of motion is performed early in the cat visual pathway.

We have used optical imaging based on intrinsic signals to explore the functional architecture of owl monkey area MT, a cortical region thought to be involved primarily in visual motion processing. As predicted by previous single-unit reports, we found cortical maps specific for the direction of moving visual stimuli. However, these direction maps were not distributed uniformly across all of area MT. Within the direction-specific regions, the activation produced by stimuli moving in opposite directions overlapped significantly. We also found that stimuli of differing shapes, moving in the same direction, activated different cortical regions within area MT, indicating that direction of motion is not the only parameter according to which area MT of owl monkey is organized. Indeed, we found clear evidence for a robust organization for orientation in area MT. Across all of MT, orientation preference changes smoothly, except at isolated line- or point-shaped discontinuities. Generally, paired regions of opposing direction preference were encompassed within a single orientation domain. The degree of segregation in the orientation maps was 3-5 times that found in direction maps. These results suggest that area MT, like V1 and V2, has a rich and multidimensional functional organization, and that orientation, a shape variable, is one of these dimensions.

In primate primary visual cortex, neurons sharing similar response properties are clustered together forming functional domains that appear as a mosaic of patches or bands, often traversing the entire cortical depth from the pia to the white matter. Similarly, each cortical site connects laterally through an extensive network of intrinsic projections that are organized in multiple clusters (patches) and reach distances of up to a few millimeters. The relationship between the functional domains and these laterally connected patches has remained a controversial issue despite intensive research efforts. To investigate this relationship, we obtained high-resolution functional maps of the cortical architecture by in vivo optical imaging. Subsequently, extracellular injections of the sensitive anterograde tracer biocytin were targeted into selected functional domains. Within the ocular dominance system, we found that long-range intrinsic connections tended to link the monocular regions of same-eye ocular dominance columns. Furthermore, we discovered that binocular domains formed a separate set of connections in area V1; binocular regions were selectively connected among themselves but were not connected to strictly monocular regions, suggesting that they constitute a distinct columnar system. In the other subsystem subserving orientation preference, patches of intrinsic connections tended to link domains sharing similar orientation preferences. Analyses of the precision of these connections indicated that in both functional subsystems,

The relationships between cytochrome oxidase blobs, ocular-dominance columns, and iso-orientation domains, subsystems underlying visual perception, were explored in primary visual cortex of macaque monkey. High-resolution maps of these three subsystems were acquired. Optical imaging based on activity-dependent intrinsic signals revealed that the most prominent organizational feature of orientation preference was a radial arrangement, forming a pinwheel-like structure surrounding a singularity point. More than 80% of these pinwheels were centered along the midline of ocular-dominance columns. The iso-orientation contours of adjacent pinwheels crossed borders of ocular-dominance columns at approximately right angles. Pinwheels with the same or opposite directions of orientation-preference change were smoothly connected with each other. On the average, all orientations were equally represented. In exactly the same cortical area, the cytochrome oxidase blobs, thought to be involved in color processing, were also mapped, using cytochrome oxidase histology. Like the centers of pinwheels, the centers of blobs also lie along the midline of ocular-dominance columns. However, the centers of pinwwheels did not coincide with the centers of blobs; these two subsystems are spatially independent. "Hypercolumn" modules, each including two complete pinwheels in two adjacent columns of complementary ocularity, as well as portions of a few blobs, were frequently found but did not seem to be the primary unit of cortical organization. An alternative to hypercolumns is proposed.

Optical imaging of the functional architecture of cortex, based on intrinsic signals, is a useful tool for the study of the development, organization, and function of the living mammalian brain. This relatively noninvasive technique is based on small activity-dependent changes of the optical properties of cortex. Thus far, functional imaging has been performed only on anesthetized animals. Here we establish that this technique is also suitable for exploring the brain of awake behaving primates. We designed a chronic sealed chamber and mounted it on the skull of a cynomolgus monkey (Macaca fascicularis) over the primary visual cortex to permit imaging through a transparent glass window. Restriction of head position alone was sufficient to eliminate movement noise in awake monkey imaging experiments. High-resolution imaging of the ocular dominance columns and the cytochrome oxidase blobs was achieved simply by taking pictures of the exposed cortex when the awake monkey was viewing video movies alternatively with each eye. Furthermore, the functional maps could be obtained without synchronization of the data acquisition to the animal's respiration and the electrocardiogram. The wavelength dependency and time course of the intrinsic signal were similar in anesthetized and awake monkeys, indicating that the signal sources were the same. We therefore conclude that optical imaging is well suited for exploring functional organization related to higher cognitive brain functions of the primate as well as providing a diagnostic tool for delineating functional cortical borders and assessing proper functions of human patients during neurosurgery.

This chapter describes a novel approach to investigate the spatiotemporal distribution of electrical activity in nervous systems. Using voltage-sensitive dyes and an electro-optical measuring system, it has recently become possible to monitor electrical activity simultaneously from multiple sites on the processes of single nerve cells, either in culture or in an intact central nervous system (CNS) in vitro, to detect the activity of many individual neurons controlling a behavioral response in invertebrate ganglia, or to follow the activity of populations of neurons at many neighboring loci in mammalian brain slices or in the intact brain. Employing optical recordings and a display processor, the images of nerve cells light up on a TV monitor when they are electrically active. Thus, the spread of electrical activity can literally be visualized in slow motion. This chapter describes recent progress in the implementation of this new technique.

Publisher Summary Optical signals related to membrane potential have been obtained with several hundred different dyes and thus, one might not expect a single mechanism to be responsible for all of the signals. The dye signals can be divided into two types with basically different mechanisms. It is convenient that dyes that work through the two types of mechanism have in turn been used in separate kinds of applications. In one situation, the optical signals appear to arise from dye that is membrane bound. These optical signals follow changes in membrane potential with time constants of less than 0.01 ms and are, thus, suited for optical monitoring of action potentials and other rapid potential changes that occur in excitable cells. This kind of signal has been termed as fast signal. Optical signals from a different group of dyes originate from the re-equilibration of charged, permeant dye molecules across the cell membrane when the membrane potential changes.