Flow detection of propagating waves with temporospatial correlation of activity

Cortex-Wide Dynamics of Intrinsic Electrical Activities: Propagating Waves and Their Interactions

Cortical circuits generate patterned activities that reflect intrinsic brain dynamics that lay the foundation for any, including stimuli-evoked, cognition and behavior. However, the spatiotemporal organization properties and principles of this intrinsic activity have only been partially elucidated because of previous poor resolution of experimental data and limited analysis methods. Here we investigated continuous wave patterns in the 0.5-4 Hz (delta band) frequency range on data from high-spatiotemporal resolution optical voltage imaging of the upper cortical layers in anesthetized mice. Waves of population activities propagate in heterogeneous directions to coordinate neuronal activities between different brain regions. The complex wave patterns show characteristics of both stereotypy and variety. The location and type of wave patterns determine the dynamical evolution when different waves interact with each other. Local wave patterns of source, sink, or saddle emerge at preferred spatial locations. Specifically, "source" patterns are predominantly found in cortical regions with low multimodal hierarchy such as the primary somatosensory cortex. Our findings reveal principles that govern the spatiotemporal dynamics of spontaneous cortical activities and associate them with the structural architecture across the cortex.SIGNIFICANCE STATEMENT Intrinsic brain activities, as opposed to external stimulus-evoked responses, have increasingly gained attention, but it remains unclear how these intrinsic activities are spatiotemporally organized at the cortex-wide scale. By taking advantage of the high spatiotemporal resolution of optical voltage imaging, we identified five wave pattern types, and revealed the organization properties of different wave patterns and the dynamical mechanisms when they interact with each other. Moreover, we found a relationship between the emergence probability of local wave patterns and the multimodal structure hierarchy across cortical areas. Our findings reveal the principles of spatiotemporal wave dynamics of spontaneous activities and associate them with the underlying hierarchical architecture across the cortex.

Detection and analysis of spatiotemporal patterns in brain activity

There is growing evidence that population-level brain activity is often organized into propagating waves that are structured in both space and time. Such spatiotemporal patterns have been linked to brain function and observed across multiple recording methodologies and scales. The ability to detect and analyze these patterns is thus essential for understanding the working mechanisms of neural circuits. Here we present a mathematical and computational framework for the identification and analysis of multiple classes of wave patterns in neural population-level recordings. By drawing a conceptual link between spatiotemporal patterns found in the brain and coherent structures such as vortices found in turbulent flows, we introduce velocity vector fields to characterize neural population activity. These vector fields are calculated for both phase and amplitude of oscillatory neural signals by adapting optical flow estimation methods from the field of computer vision. Based on these velocity vector fields, we then introduce order parameters and critical point analysis to detect and characterize a diverse range of propagating wave patterns, including planar waves, sources, sinks, spiral waves, and saddle patterns. We also introduce a novel vector field decomposition method that extracts the dominant spatiotemporal structures in a recording. This enables neural data to be represented by the activity of a small number of independent spatiotemporal modes, providing an alternative to existing dimensionality reduction techniques which separate space and time components. We demonstrate the capabilities of the framework and toolbox with simulated data, local field potentials from marmoset visual cortex and optical voltage recordings from whole mouse cortex, and we show that pattern dynamics are non-random and are modulated by the presence of visual stimuli. These methods are implemented in a MATLAB toolbox, which is freely available under an open-source licensing agreement.

Structured activity such as propagating wave patterns at the level of neural circuits can arise from highly variable firing activity of individual neurons. This property makes the brain, a quintessential example of a complex system, analogous to other complex physical systems such as turbulent fluids, in which structured patterns like vortices similarly emerge from molecules that behave irregularly. In this study, by uniquely adapting techniques for the identification of coherent structures in fluid turbulence, we develop new analytical and computational methods for the reliable detection of a diverse range of propagating wave patterns in large-scale neural recordings, for comprehensive analysis and visualization of these patterns, and for analysis of their dominant spatiotemporal modes. We demonstrate that these methods can be used to uncover the essential spatiotemporal properties of neural population activity recorded by different modalities, thus offering new insights into understanding the working mechanisms of neural systems.

High-performance, inexpensive setup for simultaneous multisite recording of electrophysiological signals and mesoscale voltage imaging in the mouse cortex

Simultaneous recording of optical and electrophysiological signals from multiple cortical areas may provide crucial information to expand our understanding of cortical function. However, the insertion of multiple electrodes into the brain may compromise optical imaging by both restricting the field of view and interfering with the approaches used to stabilize the specimen. Existing methods that combine electrophysiological recording and optical imaging in vivo implement either multiple surface electrodes, silicon probes, or a single electrode for deeper recordings. To address such limitation, we built a microelectrode array (hyperdrive, patent US5928143 A) compatible with wide-field imaging that allows insertion of up to 12 probes into a large brain area (8 mm diameter). The hyperdrive is comprised of a circle of individual microdrives where probes are positioned at an angle leaving a large brain area unobstructed for wide-field imaging. Multiple tetrodes and voltage-sensitive dye imaging were used for acute simultaneous registration of spontaneous and evoked cortical activity in anesthetized mice. The electrophysiological signals were used to extract local field potential (LFP) traces, multiunit, and single-unit spiking activity. To demonstrate our approach, we compared LFP and VSD signals over multiple regions of the cortex and analyzed the relationship between single-unit and global cortical population activities. The study of the interactions between cortical activity at local and global scales, such as the one presented in this work, can help to expand our knowledge of brain function.

Catching the voltage gradient-asymmetric boost of cortical spread generates motion signals across visual cortex: a brief review with special thanks to Amiram Grinvald

Wide-field voltage imaging is unique in its capability to capture snapshots of activity-across the full gradient of average changes in membrane potentials from subthreshold to suprathreshold levels-of hundreds of thousands of superficial cortical neurons that are simultaneously active. Here, I highlight two examples where voltage-sensitive dye imaging (VSDI) was exploited to track gradual space-time changes of activity within milliseconds across several millimeters of cortex at submillimeter resolution: the line-motion condition, measured in Amiram Grinvald's Laboratory more than 10 years ago and-coming full circle running VSDI in my laboratory-another motion-inducing condition, in which two neighboring stimuli counterchange luminance simultaneously. In both examples, cortical spread is asymmetrically boosted, creating suprathreshold activity drawn out over primary visual cortex. These rapidly propagating waves may integrate brain signals that encode motion independent of direction-selective circuits.

Improving voltage-sensitive dye imaging: with a little help from computational approaches

Voltage-sensitive dye imaging (VSDI) is a key neurophysiological recording tool because it reaches brain scales that remain inaccessible to other techniques. The development of this technique from in vitro to the behaving nonhuman primate has only been made possible thanks to the long-lasting, visionary work of Amiram Grinvald. This work has opened new scientific perspectives to the great benefit to the neuroscience community. However, this unprecedented technique remains largely under-utilized, and many future possibilities await for VSDI to reveal new functional operations. One reason why this tool has not been used extensively is the inherent complexity of the signal. For instance, the signal reflects mainly the subthreshold neuronal population response and is not linked to spiking activity in a straightforward manner. Second, VSDI gives access to intracortical recurrent dynamics that are intrinsically complex and therefore nontrivial to process. Computational approaches are thus necessary to promote our understanding and optimal use of this powerful technique. Here, we review such approaches, from computational models to dissect the mechanisms and origin of the recorded signal, to advanced signal processing methods to unravel new neuronal interactions at mesoscopic scale. Only a stronger development of interdisciplinary approaches can bridge micro- to macroscales.

Optical-flow analysis toolbox for characterization of spatiotemporal dynamics in mesoscale optical imaging of brain activity

Wide-field optical imaging techniques constitute powerful tools to investigate mesoscale neuronal activity. The sampled data constitutes a sequence of image frames in which one can investigate the flow of brain activity starting and terminating at source and sink locations respectively. Approaches to the analyses of information flow include qualitative assessment to identify sources and sinks of activity as well as their trajectories, and quantitative measurements based on computing the temporal variation of the intensity of pixels. Furthermore, in a few studies estimates of wave motion have been reported using optical-flow techniques from computer vision. However, a comprehensive toolbox for the quantitative analyses of mesoscale brain activity data is still lacking. We present a graphical-user-interface toolbox based in Matlab® for investigating the spatiotemporal dynamics of mesoscale brain activity using optical-flow analyses. The toolbox includes the implementation of three optical-flow methods namely Horn-Schunck, Combined Local-Global, and Temporospatial algorithms for estimating velocity vector fields of flow of mesoscale brain activity. From the velocity vector fields we determined the locations of sources and sinks as well as the trajectories and temporal velocities of flow of activity. Using simulated data as well as experimentally derived sensory-evoked voltage and calcium imaging data from mice, we compared the efficacy of the three optical-flow methods for determining spatiotemporal dynamics. Our results indicate that the combined local-global method we employed, yields the best results for estimating wave motion. The automated approach permits rapid and effective quantification of mesoscale brain dynamics and may facilitate the study of brain function in response to new experiences or pathology.

The stimulus-evoked population response in visual cortex of awake monkey is a propagating wave

Propagating waves occur in many excitable media and were recently found in neural systems from retina to neocortex. While propagating waves are clearly present under anaesthesia, whether they also appear during awake and conscious states remains unclear. One possibility is that these waves are systematically missed in trial-averaged data, due to variability. Here we present a method for detecting propagating waves in noisy multichannel recordings. Applying this method to single-trial voltage-sensitive dye imaging data, we show that the stimulus-evoked population response in primary visual cortex of the awake monkey propagates as a travelling wave, with consistent dynamics across trials. A network model suggests that this reliability is the hallmark of the horizontal fibre network of superficial cortical layers. Propagating waves with similar properties occur independently in secondary visual cortex, but maintain precise phase relations with the waves in primary visual cortex. These results show that, in response to a visual stimulus, propagating waves are systematically evoked in several visual areas, generating a consistent spatiotemporal frame for further neuronal interactions.

Propagating waves of cortical neuronal activity are implicated in various cognitive processes and have been observed in anaesthetised animals. Here, the authors demonstrate the existence of propagating waves in awake monkeys during visual stimulation, and show that they are mediated by horizontal fibres in the cortex.

Spontaneous cortical activity alternates between motifs defined by regional axonal projections

Using millisecond-timescale voltage-sensitive dye imaging in lightly anesthetized or awake adult mice, we show that a palette of sensory-evoked and hemisphere-wide activity motifs are represented in spontaneous activity. These motifs can reflect multiple modes of sensory processing, including vision, audition and touch. We found similar cortical networks with direct cortical activation using channelrhodopsin-2. Regional analysis of activity spread indicated modality-specific sources, such as primary sensory areas, a common posterior-medial cortical sink where sensory activity was extinguished within the parietal association area and a secondary anterior medial sink within the cingulate and secondary motor cortices for visual stimuli. Correlation analysis between functional circuits and intracortical axonal projections indicated a common framework corresponding to long-range monosynaptic connections between cortical regions. Maps of intracortical monosynaptic structural connections predicted hemisphere-wide patterns of spontaneous and sensory-evoked depolarization. We suggest that an intracortical monosynaptic connectome shapes the ebb and flow of spontaneous cortical activity.

Wave propagation of cortical population activity under urethane anesthesia is state dependent

BACKGROUND

Propagating waves of excitation have been observed extensively in the neocortex, during both spontaneous and sensory-evoked activity, and they play a critical role in spatially organizing information processing. However, the state-dependence of these spatiotemporal propagation patterns is largely unexplored. In this report, we use voltage-sensitive dye imaging in the rat visual cortex to study the propagation of spontaneous population activity in two discrete cortical states induced by urethane anesthesia.

RESULTS

While laminar current source density patterns of spontaneous population events in these two states indicate a considerable degree of similarity in laminar networks, lateral propagation in the more active desynchronized state is approximately 20% faster than in the slower synchronized state. Furthermore, trajectories of wave propagation exhibit a strong anisotropy, but the preferred direction is different depending on cortical state.

CONCLUSIONS

Our results show that horizontal wave propagation of spontaneous neural activity is largely dependent on the global activity states of local cortical circuits.

Stimulus-entrained oscillatory activity propagates as waves from area 18 to 17 in cat visual cortex

Previous studies in cat visual cortex reported that area 18 can actively drive neurons in area 17 through cortico-cortical projections. However, the dynamics of such cortico-cortical interaction remains unclear. Here we used multielectrode arrays to examine the spatiotemporal pattern of neuronal activity in cat visual cortex across the 17/18 border. We found that full-field contrast reversal gratings evoked oscillatory wave activity propagating from area 18 to 17. The wave direction was independent of the grating orientation, and could not be accounted for by the spatial distribution of receptive field latencies, suggesting that the waves are largely mediated by intrinsic connections in the cortex. Different from the evoked waves, spontaneous waves propagated along both directions across the 17/18 border. Together, our results suggest that visual stimulation may enhance the flow of information from area 18 to 17.

In vivo imaging of epileptic activity using 2-NBDG, a fluorescent deoxyglucose analog

Accurately locating epileptic foci has great importance in advancing the treatment of epilepsy. In this study, epileptic seizures were first induced by intracortical injection of 4-aminopyridine in rats. A fluorescent deoxyglucose substitute, 2-(N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino)-2-deoxyglucose (2-NBDG), was then continuously injected via the tail vein. Brain glucose metabolism was subsequently monitored by fluorescence imaging of 2-NBDG. The initial uptake rate of 2-NBDG at the injection site of 4-aminopyridine significantly exceeded that of the control injection site, which indicated local hypermetabolism induced by seizures. Our results show that 2-NBDG can be used for localizing epileptic foci.