Interneurons contribute to the hemodynamic/metabolic response to epileptiform discharges.

Interpretation of hemodynamic responses in epilepsy is hampered by an incomplete understanding of the underlying neurovascular coupling, especially the contributions of excitation and inhibition. We made simultaneous multimodal recordings of local field potentials (LFPs), firing of individual neurons, blood flow, and oxygen level in the somatosensory cortex of anesthetized rats. Epileptiform discharges induced by bicuculline injections were used to trigger large local events. LFP and blood flow were robustly coupled, as were LFP and tissue oxygen. In a parametric linear model, LFP and the baseline activities of cerebral blood flow and tissue partial oxygen tension contributed significantly to blood flow and oxygen responses. In an analysis of recordings from 402 neurons, blood flow/tissue oxygen correlated with the discharge of putative interneurons but not of principal cells. Our results show that interneuron activity is important in the vascular and metabolic responses during epileptiform discharges.

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.

Motion clouds: model-based stimulus synthesis of natural-like random textures for the study of motion perception.

Choosing an appropriate set of stimuli is essential to characterize the response of a sensory system to a particular functional dimension, such as the eye movement following the motion of a visual scene. Here, we describe a framework to generate random texture movies with controlled information content, i.e., Motion Clouds. These stimuli are defined using a generative model that is based on controlled experimental parametrization. We show that Motion Clouds correspond to dense mixing of localized moving gratings with random positions. Their global envelope is similar to natural-like stimulation with an approximate full-field translation corresponding to a retinal slip. We describe the construction of these stimuli mathematically and propose an open-source Python-based implementation. Examples of the use of this framework are shown. We also propose extensions to other modalities such as color vision, touch, and audition.

Two-photon GCaMP6f imaging of infrared neural stimulation evoked calcium signals in mouse cortical neurons in vivo.

Infrared neural stimulation is a promising tool for stimulating the brain because it can be used to excite with high spatial precision without the need of delivering or inserting any exogenous agent into the tissue. Very few studies have explored its use in the brain, as most investigations have focused on sensory or motor nerve stimulation. Using intravital calcium imaging with the genetically encoded calcium indicator GCaMP6f, here we show that the application of infrared neural stimulation induces intracellular calcium signals in Layer 2/3 neurons in mouse cortex in vivo. The number of neurons exhibiting infrared-induced calcium response as well as the amplitude of those signals are shown to be both increasing with the energy density applied. By studying as well the spatial extent of the stimulation, we show that reproducibility of the stimulation is achieved mainly in the central part of the infrared beam path. Stimulating in vivo at such a degree of precision and without any exogenous chromophores enables multiple applications, from mapping the brain's connectome to applications in systems neuroscience and the development of new therapeutic tools for investigating the pathological brain.

Non-thermal Electroporation Ablation of Epileptogenic Zones Stops Seizures in Mice While Providing Reduced Vascular Damage and Accelerated Tissue Recovery.

In epilepsy, the most frequent surgical procedure is the resection of brain tissue in the temporal lobe, with seizure-free outcomes in approximately two-thirds of cases. However, consequences of surgery can vary strongly depending on the brain region targeted for removal, as surgical morbidity and collateral damage can lead to significant complications, particularly when bleeding and swelling are located near delicate functional cortical regions. Although focal thermal ablations are well-explored in epilepsy as a minimally invasive approach, hemorrhage and edema can be a consequence as the blood-brain barrier is still disrupted. Non-thermal irreversible electroporation (NTIRE), common in many other medical tissue ablations outside the brain, is a relatively unexplored method for the ablation of neural tissue, and has never been reported as a means for ablation of brain tissue in the context of epilepsy. Here, we present a detailed visualization of non-thermal ablation of neural tissue in mice and report that NTIRE successfully ablates epileptic foci in mice, resulting in seizure-freedom, while causing significantly less hemorrhage and edema compared to conventional thermal ablation. The NTIRE approach to ablation preserves the blood-brain barrier while pathological circuits in the same region are destroyed. Additionally, we see the reinnervation of fibers into ablated brain regions from neighboring areas as early as day 3 after ablation. Our evidence demonstrates that NTIRE could be utilized as a precise tool for the ablation of surgically challenging epileptogenic zones in patients where the risk of complications and hemorrhage is high, allowing not only reduced tissue damage but potentially accelerated recovery as vessels and extracellular matrix remain intact at the point of ablation.

Investigation of linear coupling between single-event blood flow responses and interictal discharges in a model of experimental epilepsy.

A successful outcome of epilepsy neurosurgery relies on an accurate delineation of the epileptogenic region to be resected. Functional magnetic resonance imaging (fMRI) would allow doing this noninvasively at high spatial resolution. However, a clear, quantitative description of the relationship between hemodynamic changes and the underlying epileptiform neuronal activity is still missing, thereby preventing the systematic use of fMRI for routine epilepsy surgery planning. To this aim, we used a local epilepsy model to record simultaneously cerebral blood flow (CBF) with laser Doppler (LD) and local field potentials (LFP) in rat frontal cortex. CBF responses to individual interictal-like spikes were large and robust. Their amplitude correlated linearly with spike amplitude. Moreover, the CBF response added linearly in time over a large range of spiking rates. CBF responses could thus be predicted by a linear model of the kind currently used for the interpretation of fMRI data, but including also the spikes' amplitudes as additional information. Predicted and measured CBF responses matched accurately. For high spiking frequencies (above approximately 0.2 Hz), the responses saturated but could eventually recover, indicating the presence of multiple neurovascular coupling mechanisms, which might act at different spatiotemporal scales. Spatially, CBF responses peaked at the center of epileptic activity and displayed a spatial specificity at least as good as the millimeter. These results suggest that simultaneous electroencephalographic and blood flow-based fMRI recordings should be suitable for the noninvasive precise localization of hyperexcitable regions in epileptic patients candidate for neurosurgery.

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.

Multimodal Characterization of Neural Networks Using Highly Transparent Electrode Arrays.

Transparent and flexible materials are attractive for a wide range of emerging bioelectronic applications. These include neural interfacing devices for both recording and stimulation, where low electrochemical electrode impedance is valuable. Here the conducting polymer poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) is used to fabricate electrodes that are small enough to allow unencumbered optical access for imaging a large cell population with two-photon (2P) microscopy, yet provide low impedance for simultaneous high quality recordings of neural activity in vivo. To demonstrate this, pathophysiological activity was induced in the mouse cortex using 4-aminopyridine (4AP), and the resulting electrical activity was detected with the PEDOT:PSS-based probe while imaging calcium activity directly below the probe area. The induced calcium activity of the neuronal network as measured by the fluorescence change in the cells correlated well with the electrophysiological recordings from the cortical grid of PEDOT:PSS microelectrodes. Our approach provides a valuable vehicle for complementing classical high temporal resolution electrophysiological analysis with optical imaging.

Visual stimulation quenches global alpha range activity in awake primate V4: a case study.

Increasing evidence suggests that sensory stimulation not only changes the level of cortical activity with respect to baseline but also its structure. Despite having been reported in a multitude of conditions and preparations (for instance, as a quenching of intertrial variability, Churchland et al., 2010), such changes remain relatively poorly characterized. Here, we used optical imaging of voltage-sensitive dyes to explore, in V4 of an awake macaque, the spatiotemporal characteristics of both visually evoked and spontaneously ongoing neuronal activity and their difference. With respect to the spontaneous case, we detected a reduction in large-scale activity ([Formula: see text]) in the alpha range (5 to 12.5 Hz) during sensory inflow accompanied by a decrease in pairwise correlations. Moreover, the spatial patterns of correlation obtained during the different visual stimuli were on the average more similar one to another than they were to that obtained in the absence of stimulation. Finally, these observed changes in activity dynamics approached saturation already at very low stimulus contrasts, unlike the progressive, near-linear increase of the mean raw evoked responses over a wide range of contrast values, which could indicate a specific switching in the presence of a sensory inflow.

Compartment-resolved imaging of activity-dependent dynamics of cortical blood volume and oximetry.

Optical imaging, positron emission tomography, and functional magnetic resonance imaging (fMRI) all rely on vascular responses to image neuronal activity. Although these imaging techniques are used successfully for functional brain mapping, the detailed spatiotemporal dynamics of hemodynamic events in the various microvascular compartments have remained unknown. Here we used high-resolution optical imaging in area 18 of anesthetized cats to selectively explore sensory-evoked cerebral blood-volume (CBV) changes in the various cortical microvascular compartments. To avoid the confounding effects of hematocrit and oximetry changes, we developed and used a new fluorescent blood plasma tracer and combined these measurements with optical imaging of intrinsic signals at a near-isosbestic wavelength for hemoglobin (565 nm). The vascular response began at the arteriolar level, rapidly spreading toward capillaries and venules. Larger veins lagged behind. Capillaries exhibited clear blood-volume changes. Arterioles and arteries had the largest response, whereas the venous response was smallest. Information about compartment-specific oxygen tension dynamics was obtained in imaging sessions using 605 nm illumination, a wavelength known to reflect primarily oximetric changes, thus being more directly related to electrical activity than CBV changes. Those images were radically different: the response began at the parenchyma level, followed only later by the other microvascular compartments. These results have implications for the modeling of fMRI responses (e.g., the balloon model). Furthermore, functional maps obtained by imaging the capillary CBV response were similar but not identical to those obtained using the early oximetric signal, suggesting the presence of different regulatory mechanisms underlying these two hemodynamic processes.

Hemodynamic responses in cortex investigated with optical imaging methods. Implications for functional brain mapping.

During the last 20 years, optical imaging methods - either alone or in combination with other recording techniques - has proven a fruitful approach to explore both the physiological and the functional aspects of activity-evoked hemodynamic responses in cortex. One of the main advantages of optical imaging consists in its high spatio-temporal resolution (in the order of few microns and milliseconds), allowing not only to unambiguously distinguish between activity patterns relating to the underlying functional architecture and those originating from the activation of medium/large blood vessels, but also to investigate the various activity-evoked hemodynamic processes at very fine detail. Here, we briefly review the principal findings obtained by optical imaging about the spatio-temporal properties of the various hemodynamic responses in cortex, i.e., changes in blood-oxygenation, blood-volume, and, to some extent, blood-flow. We also discuss the implications of those findings for non-invasive high-resolution functional brain imaging, in particular for fMRI. Finally, we underscore the importance of novel approaches for high-resolution blood-flow imaging, in the context of the need to gather information at fine spatial detail about the blood-flow response, necessary to constrain the multiple free parameters of hemodynamic response models.

Imaging the Neocortex Functional Architecture Using Multiple Intrinsic Signals: Implications for Hemodynamic-Based Functional Imaging.

Optical imaging based on intrinsic signals has provided a new level of understanding of the principles underlying cortical development, organization, and function, providing a spatial resolution of up to 20 µm for mapping cortical columns in vivo. This introduction briefly reviews the development of this technique, the types of applications that have been pursued, and the general implications of some findings for other neuroimaging techniques based on hemodynamic responses (e.g., functional magnetic resonance imaging).

Voltage-Sensitive Dye Imaging of Neocortical Activity.

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 the dynamics of local and large range populations. Voltage-sensitive dye imaging (VSDI), based on organic voltage probes, reveals neural population activity in areas ranging from a few tens of micrometers to a couple of centimeters, or two areas up to ~10 cm apart. VSDI provides a submillisecond temporal resolution and a spatial resolution of ~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 to explore the lateral spread of retinotopic or somatotopic activation; the dynamic spatiotemporal pattern resulting from sensory activation, including the somatosensory, olfactory, auditory, and visual modalities; and motor preparation and the properties of spontaneously occurring population activity. In this introduction, we focus on VSDI in vivo and review results obtained mostly in the visual system in our laboratory.

Accurate spike estimation from noisy calcium signals for ultrafast three-dimensional imaging of large neuronal populations in vivo.

Extracting neuronal spiking activity from large-scale two-photon recordings remains challenging, especially in mammals in vivo, where large noises often contaminate the signals. We propose a method, MLspike, which returns the most likely spike train underlying the measured calcium fluorescence. It relies on a physiological model including baseline fluctuations and distinct nonlinearities for synthetic and genetically encoded indicators. Model parameters can be either provided by the user or estimated from the data themselves. MLspike is computationally efficient thanks to its original discretization of probability representations; moreover, it can also return spike probabilities or samples. Benchmarked on extensive simulations and real data from seven different preparations, it outperformed state-of-the-art algorithms. Combined with the finding obtained from systematic data investigation (noise level, spiking rate and so on) that photonic noise is not necessarily the main limiting factor, our method allows spike extraction from large-scale recordings, as demonstrated on acousto-optical three-dimensional recordings of over 1,000 neurons in vivo.

Special report: Noninvasive multi-parameter functional optical imaging of the eye.

Advancement in the treatment of blindness depends on the development of new technologies that enable early detection, follow-up, and treatment of disease. The authors describe direct, noninvasive imaging of four parameters: blood flow, blood oximetry, metabolic state, and hidden vasculature, particularly capillaries. These are functional parameters of the retina known to be degraded by retinal disease. The new Retinal Function Imager (Optical Imaging, Ltd., Rehovot, Israel) can image all four parameters as intrinsic reflectance intensity differences over the retina's surface. During the past 2 decades, imaging of small optical signals has been a powerful tool for high-resolution functional mapping in the neocortex. In this article, this technology is applied to the retina and demonstrates a general tool for noninvasively probing retinal function in many modalities. Imaging functional changes before anatomic consequences arise holds promise as a powerful tool for early diagnosis and treatment of retinal disease.

High-resolution functional optical imaging: from the neocortex to the eye.

Much work remains to be done to establish the clinical usefulness of the RFI for early diagnosis and treatment guidance. The discoveries obtained by functional optical imaging of the neocortex in the last 15 years and the recent RFI studies of the eyes of normal subjects and patients with diabetic retinopathy, glaucoma, and age-related macular degeneration suggest that functional optical imaging of the retina is likely to become a multi-modality powerful clinical tool.

Elliptically polarized light for depth resolved optical imaging.

It is shown that using elliptically polarized light permits selecting well-defined subsurface volumes in a turbid medium. This suggests the possibility of probing biological tissues at specific depths. First, we present the method and preliminary results obtained on an Intralipid phantom. We next report on the method's performance on a biological phantom (chicken breast) and, finally, on the exposed cortex of an anesthetized rat.

A new variational method for erythrocyte velocity estimation in wide-field imaging in vivo.

Measuring erythrocyte velocity in individual microvessels has important applications for biomedical and functional imaging. Recent multiphoton fluorescence microscopy approaches require injecting fluorescent tracers; moreover, only one or few vessels can be imaged at a time. To overcome these shortcomings, we used CCD-based optical imaging of intrinsic absorption changes in macroscopic vascular networks to record erythrocytes' trajectories over several mm (2) of cortical surface. We then demonstrate the feasibility of erythrocyte velocity estimation from such wide-field data, using two robust, independent, algorithms. The first one is a recently published Radon transform-based algorithm that estimates erythrocyte velocity locally. We adapt it to data obtained in wide-field imaging and show, for the first time, its performance on such datasets. The second ("fasttrack") algorithm is novel. It is based on global energy minimization techniques to estimate the full spatiotemporal erythrocytes' trajectories inside vessels. We test the two algorithms on both simulated and biological data, obtained in rat cerebral cortex in a spreading depression experiment. On vessels with medium-slow erythrocyte velocities both algorithms performed well, allowing their usage as benchmark one for another. However, our novel fasttrack algorithm outperformed the other one for higher velocities, as encountered in the arterial network.

Coupling between neuronal activity and microcirculation: implications for functional brain imaging.

In the neocortex, neurons with similar response properties are often clustered together in column-like structures, giving rise to what has become known as functional architecture-the mapping of various stimulus feature dimensions onto the cortical sheet. At least partially, we owe this finding to the availability of several functional brain imaging techniques, both post-mortem and in-vivo, which have become available over the last two generations, revolutionizing neuroscience by yielding information about the spatial organization of active neurons in the brain. Here, we focus on how our understanding of such functional architecture is linked to the development of those functional imaging methodologies, especially to those that image neuronal activity indirectly, through metabolic or haemodynamic signals, rather than directly through measurement of electrical activity. Some of those approaches allow exploring functional architecture at higher spatial resolution than others. In particular, optical imaging of intrinsic signals reaches the striking detail of approximately 50 mum, and, together with other methodologies, it has allowed characterizing the metabolic and haemodynamic responses induced by sensory-evoked neuronal activity. Here, we review those findings about the spatio-temporal characteristics of neurovascular coupling and discuss their implications for functional brain imaging, including position emission tomography, and non-invasive neuroimaging techniques, such as funtional magnetic resonance imaging, applicable also to the human brain.

Behavioral receptive field for ocular following in humans: dynamics of spatial summation and center-surround interactions.

Visual neurons integrate information over a finite part of the visual field with high selectivity. This classical receptive field is modulated by peripheral inputs that play a role in both neuronal response normalization and contextual modulations. However, the consequences of these properties for visuomotor transformations are yet incompletely understood. To explore those, we recorded short-latency ocular following responses in humans to large center-only and center-surround stimuli. We found that eye movements are triggered by a mechanism that integrates motion over a restricted portion of the visual field, the size of which depends on stimulus contrast and increases as a function of time after response onset. We also found evidence for a strong nonisodirectional center-surround organization, responsible for normalizing the central, driving input so that motor responses are set to their most linear contrast dynamics. Such response normalization is delayed about 20 ms relative to tracking onset, gradually builds up over time, and is partly tuned for surround orientation/direction. These results outline the spatiotemporal organization of a behavioral receptive field, which might reflect a linear integration among subpopulations of cortical visual motion detectors.