Traveling waves in visual cortex

A Translaminar Spacetime Code Supports Touch-Evoked Traveling Waves

Linking sensory-evoked traveling waves to underlying circuit patterns is critical to understanding the neural basis of sensory perception. To form this link, we performed simultaneous electrophysiology and two-photon calcium imaging through transparent NeuroGrids and mapped touch-evoked cortical traveling waves and their underlying microcircuit dynamics. In awake mice, both passive and active whisker touch elicited traveling waves within and across barrels, with a fast early component followed by a variable late wave that lasted hundreds of milliseconds post-stimulus. Strikingly, late-wave dynamics were modulated by stimulus value and correlated with task performance. Mechanistically, the late wave component was i) modulated by motor feedback, ii) complemented by a sparse ensemble pattern across layer 2/3, which a balanced-state network model reconciled via inhibitory stabilization, and iii) aligned to regenerative Layer-5 apical dendritic Ca 2+ events. Our results reveal a translaminar spacetime pattern organized by cortical feedback in the sensory cortex that supports touch-evoked traveling waves.

GRAPHICAL ABSTRACT AND HIGHLIGHTS

Whisker touch evokes both early- and late-traveling waves in the barrel cortex over 100's of millisecondsReward reinforcement modulates wave dynamics Late wave emergence coincides with network sparsity in L23 and time-locked L5 dendritic Ca 2+ spikes Experimental and computational results link motor feedback to distinct translaminar spacetime patterns.

A Traveling Waves Perspective on Temporal Binding

Brain oscillations are involved in many cognitive processes, and several studies have investigated their role in cognition. In particular, the phase of certain oscillations has been related to temporal binding and integration processes, with some authors arguing that perception could be an inherently rhythmic process. However, previous research on oscillations mostly overlooked their spatial component: how oscillations propagate through the brain as traveling waves, with systematic phase delays between brain regions. Here, we argue that interpreting oscillations as traveling waves is a useful paradigm shift to understand their role in temporal binding and address controversial results. After a brief definition of traveling waves, we propose an original view on temporal integration that considers this new perspective. We first focus on cortical dynamics, then speculate about the role of thalamic nuclei in modulating the waves, and on the possible consequences for rhythmic temporal binding. In conclusion, we highlight the importance of considering oscillations as traveling waves when investigating their role in cognitive functions.

The direction of theta and alpha travelling waves modulates human memory processing

To support a range of behaviours, the brain must flexibly coordinate neural activity across widespread brain regions. One potential mechanism for this coordination is a travelling wave, in which a neural oscillation propagates across the brain while organizing the order and timing of activity across regions. Although travelling waves are present across the brain in various species, their potential functional relevance has remained unknown. Here, using rare direct human brain recordings, we demonstrate a distinct functional role for travelling waves of theta- and alpha-band (2-13 Hz) oscillations in the cortex. Travelling waves propagate in different directions during separate cognitive processes. In episodic memory, travelling waves tended to propagate in a posterior-to-anterior direction during successful memory encoding and in an anterior-to-posterior direction during recall. Because travelling waves of oscillations correspond to local neuronal spiking, these patterns indicate that rhythmic pulses of activity move across the brain in different directions for separate behaviours. More broadly, our results suggest a fundamental role for travelling waves and oscillations in dynamically coordinating neural connectivity, by flexibly organizing the timing and directionality of network interactions across the cortex to support cognition and behaviour.

Sequentially activated discrete modules appear as traveling waves in neuronal measurements with limited spatiotemporal sampling

Numerous studies have identified traveling waves in the cortex and suggested they play important roles in brain processing. These waves are most often measured using macroscopic methods that are unable to assess the local spiking activity underlying wave dynamics. Here, we investigated the possibility that waves may not be traveling at the single neuron scale. We first show that sequentially activating two discrete brain areas can appear as traveling waves in EEG simulations. We next reproduce these results using an analytical model of two sequentially activated regions. Using this model, we were able to generate wave-like activity with variable directions, velocities, and spatial patterns, and to map the discriminability limits between traveling waves and modular sequential activations. Finally, we investigated the link between field potentials and single neuron excitability using large-scale measurements from turtle cortex ex vivo. We found that while field potentials exhibit wave-like dynamics, the underlying spiking activity was better described by consecutively activated spatially adjacent groups of neurons. Taken together, this study suggests caution when interpreting phase delay measurements as continuously propagating wavefronts in two different spatial scales. A careful distinction between modular and wave excitability profiles across scales will be critical for understanding the nature of cortical computations.

The pulvinar as a hub of visual processing and cortical integration

The pulvinar nucleus of the thalamus is a crucial component of the visual system and plays significant roles in sensory processing and cognitive integration. The pulvinar's extensive connectivity with cortical regions allows for bidirectional communication, contributing to the integration of sensory information across the visual hierarchy. Recent findings underscore the pulvinar's involvement in attentional modulation, feature binding, and predictive coding. In this review, we highlight recent advances in clarifying the pulvinar's circuitry and function. We discuss the contributions of the pulvinar to signal modulation across the global cortical network and place these findings within theoretical frameworks of cortical processing, particularly the global neuronal workspace (GNW) theory and predictive coding.

Neural mechanisms of visual motion extrapolation

Because neural processing takes time, the brain only has delayed access to sensory information. When localising moving objects this is problematic, as an object will have moved on by the time its position has been determined. Here, we consider predictive motion extrapolation as a fundamental delay-compensation strategy. From a population-coding perspective, we outline how extrapolation can be achieved by a forwards shift in the population-level activity distribution. We identify general mechanisms underlying such shifts, involving various asymmetries which facilitate the targeted 'enhancement' and/or 'dampening' of population-level activity. We classify these on the basis of their potential implementation (intra- vs inter-regional processes) and consider specific examples in different visual regions. We consider how motion extrapolation can be achieved during inter-regional signaling, and how asymmetric connectivity patterns which support extrapolation can emerge spontaneously from local synaptic learning rules. Finally, we consider how more abstract 'model-based' predictive strategies might be implemented. Overall, we present an integrative framework for understanding how the brain determines the real-time position of moving objects, despite neural delays.

Perceptual Cycles Travel Across Retinotopic Space

Visual perception waxes and wanes periodically over time at low frequencies (theta: 4-7 Hz; alpha: 8-13 Hz), creating "perceptual cycles." These perceptual cycles can be induced when stimulating the brain with a flickering visual stimulus at the theta or alpha frequency. Here, we took advantage of the well-known organization of the visual system into retinotopic maps (topographic correspondence between visual and cortical spaces) to assess the spatial organization of induced perceptual cycles. Specifically, we tested the hypothesis that they can propagate across the retinotopic space. A disk oscillating in luminance (inducer) at 4, 6, 8, or 10 Hz was presented in the periphery of the visual field to induce perceptual cycles at specific frequencies. EEG recordings verified that the brain responded at the corresponding inducer frequencies and their first harmonics. Perceptual cycles were assessed with a concurrent detection task-target stimuli were displayed at threshold contrast (50% detection) at random times during the inducer. Behavioral results confirmed that perceptual performance was modulated periodically by the inducer at each frequency. We additionally manipulated the distance between the target and the inducer (three possible positions) and showed that the optimal phase, that is, moment of highest target detection, shifted across target distance to the inducer, specifically when its flicker frequency was in the alpha range (8 and 10 Hz). These results demonstrate that induced alpha perceptual cycles travel across the retinotopic space in humans at a propagation speed of 0.3-0.5 m/sec, consistent with the speed of unmyelinated horizontal connections in the visual cortex.

Phantom oscillations in principal component analysis

Principal component analysis (PCA) is a dimensionality reduction method that is known for being simple and easy to interpret. Principal components are often interpreted as low-dimensional patterns in high-dimensional space. However, this simple interpretation fails for timeseries, spatial maps, and other continuous data. In these cases, nonoscillatory data may have oscillatory principal components. Here, we show that two common properties of data cause oscillatory principal components: smoothness and shifts in time or space. These two properties implicate almost all neuroscience data. We show how the oscillations produced by PCA, which we call "phantom oscillations," impact data analysis. We also show that traditional cross-validation does not detect phantom oscillations, so we suggest procedures that do. Our findings are supported by a collection of mathematical proofs. Collectively, our work demonstrates that patterns which emerge from high-dimensional data analysis may not faithfully represent the underlying data.

Beta traveling waves in monkey frontal and parietal areas encode recent reward history

Brain function depends on neural communication, but the mechanisms of this communication are not well understood. Recent studies suggest that one form of neural communication is through traveling waves (TWs)-patterns of neural oscillations that propagate within and between brain areas. We show that TWs are robust in microarray recordings in frontal and parietal cortex and encode recent reward history. Two adult male monkeys made saccades to obtain probabilistic rewards and were sensitive to the (statistically irrelevant) reward on the previous trial. TWs in frontal and parietal areas were stronger in trials that followed a prior reward versus a lack of reward and, in the frontal lobe, correlated with the monkeys' behavioral sensitivity to the prior reward. The findings suggest that neural communication mediated by TWs within the frontal and parietal lobes contribute to maintaining information about recent reward history and mediating the impact of this history on the monkeys' expectations.

Rogue Waves: A Unifying Model of Human Seizure Propagation

Multiple Sources of Fast Traveling Waves During Human Seizures: Resolving a Controversy Schlafly ED, Marshall FA, Merricks EM, Eden UT, Cash SS, Schevon CA, Kramer MA. J Neurosci . 2022;42(36):6966-6982. doi:10.1523/JNEUROSCI.0338-22.2022 During human seizures organized waves of voltage activity rapidly sweep across the cortex. Two contradictory theories describe the source of these fast traveling waves: either a slowly advancing narrow region of multiunit activity (an ictal wavefront) or a fixed cortical location. Limited observations and different analyses prevent resolution of these incompatible theories. Here we address this disagreement by combining the methods and microelectrode array recordings (N = 11 patients, 2 females, N = 31 seizures) from previous human studies to analyze the traveling wave source. We find—inconsistent with both existing theories—a transient relationship between the ictal wavefront and traveling waves, and multiple stable directions of traveling waves in many seizures. Using a computational model that combines elements of both existing theories, we show that interactions between an ictal wavefront and fixed source reproduce the traveling wave dynamics observed in vivo. We conclude that combining both existing theories can generate the diversity of ictal traveling waves. Significance Statement: The source of voltage discharges that propagate across cortex during human seizures remains unknown. Two candidate theories exist, each proposing a different discharge source. Support for each theory consists of observations from a small number of human subject recordings, analyzed with separately developed methods. How the different, limited data and different analysis methods impact the evidence for each theory is unclear. To resolve these differences, we combine the unique, human microelectrode array recordings collected separately for each theory and analyze these combined data with a unified approach. We show that neither existing theory adequately describes the data. We then propose a new theory that unifies existing proposals and successfully reproduces the voltage discharge dynamics observed in vivo.

Waves traveling over a map of visual space can ignite short-term predictions of sensory input

Recent analyses have found waves of neural activity traveling across entire visual cortical areas in awake animals. These traveling waves modulate the excitability of local networks and perceptual sensitivity. The general computational role of these spatiotemporal patterns in the visual system, however, remains unclear. Here, we hypothesize that traveling waves endow the visual system with the capacity to predict complex and naturalistic inputs. We present a network model whose connections can be rapidly and efficiently trained to predict individual natural movies. After training, a few input frames from a movie trigger complex wave patterns that drive accurate predictions many frames into the future solely from the network's connections. When the recurrent connections that drive waves are randomly shuffled, both traveling waves and the ability to predict are eliminated. These results suggest traveling waves may play an essential computational role in the visual system by embedding continuous spatiotemporal structures over spatial maps.

Perceptual echoes as travelling waves may arise from two discrete neuronal sources

Growing evidence suggests that travelling waves are functionally relevant for cognitive operations in the brain. Several electroencephalography (EEG) studies report on a perceptual alpha-echo, representing the brain response to a random visual flicker, propagating as a travelling wave across the cortical surface. In this study, we ask if the propagating activity of the alpha-echo is best explained by a set of discrete sources mixing at the sensor level rather than a cortical travelling wave. To this end, we presented participants with gratings modulated by random noise and simultaneously acquired the ongoing MEG. The perceptual alpha-echo was estimated using the temporal response function linking the visual input to the brain response. At the group level, we observed a spatial decay of the amplitude of the alpha-echo with respect to the sensor where the alpha-echo was the largest. Importantly, the propagation latencies consistently increased with the distance. Interestingly, the propagation of the alpha-echoes was predominantly centro-lateral, while EEG studies reported mainly posterior-frontal propagation. Moreover, the propagation speed of the alpha-echoes derived from the MEG data was around 10 m/s, which is higher compared to the 2 m/s reported in EEG studies. Using source modelling, we found an early component in the primary visual cortex and a phase-lagged late component in the parietal cortex, which may underlie the travelling alpha-echoes at the sensor level. We then simulated the alpha-echoes using realistic EEG and MEG forward models by placing two sources in the parietal and occipital cortices in accordance with our empirical findings. The two-source model could account for both the direction and speed of the observed alpha-echoes in the EEG and MEG data. Our results demonstrate that the propagation of the perceptual echoes observed in EEG and MEG data can be explained by two sources mixing at the scalp level equally well as by a cortical travelling wave. Importantly, these findings should not be directly extrapolated to intracortical recordings, where travelling waves gradually propagate at a sub-millimetre scale.

Near-field potentials index local neural computations more accurately than population spiking

Local field potentials (LFP) are low-frequency extracellular voltage fluctuations thought to primarily arise from synaptic activity. However, unlike highly localized neuronal spiking, LFP is spatially less specific. LFP measured at one location is not entirely generated there due to far-field contributions that are passively conducted across volumes of neural tissue. We sought to quantify how much information within the locally generated, near-field low-frequency activity (nfLFP) is masked by volume-conducted far-field signals. To do so, we measured laminar neural activity in primary visual cortex (V1) of monkeys viewing sequences of multifeatured stimuli. We compared information content of regular LFP and nfLFP that was mathematically stripped of volume-conducted far-field contributions. Information content was estimated by decoding stimulus properties from neural responses via spatiotemporal multivariate pattern analysis. Volume-conducted information differed from locally generated information in two important ways: (1) for stimulus features relevant to V1 processing (orientation and eye-of-origin), nfLFP contained more information. (2) in contrast, the volume-conducted signal was more informative regarding temporal context (relative stimulus position in a sequence), a signal likely to be coming from elsewhere. Moreover, LFP and nfLFP differed both spectrally as well as spatially, urging caution regarding the interpretations of individual frequency bands and/or laminar patterns of LFP. Most importantly, we found that population spiking of local neurons was less informative than either the LFP or nfLFP, with nfLFP containing most of the relevant information regarding local stimulus processing. These findings suggest that the optimal way to read out local computational processing from neural activity is to decode the local contributions to LFP, with significant information loss hampering both regular LFP and local spiking.

Author’s Contributions

Conceptualization, D.A.T., J.A.W, and A.M.; Data Collection, J.A.W., M.A.C., K.D.; Formal Analysis, D.A.T. and J.A.W.; Data Visualization, D.A.T. and J.A.W.; Original Draft, D.A.T., J.A.W., and A.M.; Revisions and Final Draft, D.A.T., J.A.W., M.A.C., K.D., M.T.W., A.M.B., and A.M.

Competing Interests

The authors declare no conflicts of interest.

Human acute microelectrode array recordings with broad cortical access, single-unit resolution, and parallel behavioral monitoring

There are vast gaps in our understanding of the organization and operation of the human nervous system at the level of individual neurons and their networks. Here, we report reliable and robust acute multichannel recordings using planar microelectrode arrays (MEAs) implanted intracortically in awake brain surgery with open craniotomies that grant access to large parts of the cortical hemisphere. We obtained high-quality extracellular neuronal activity at the microcircuit, local field potential level and at the cellular, single-unit level. Recording from the parietal association cortex, a region rarely explored in human single-unit studies, we demonstrate applications on these complementary spatial scales and describe traveling waves of oscillatory activity as well as single-neuron and neuronal population responses during numerical cognition, including operations with uniquely human number symbols. Intraoperative MEA recordings are practicable and can be scaled up to explore cellular and microcircuit mechanisms of a wide range of human brain functions.

Coherent noise enables probabilistic sequence replay in spiking neuronal networks

Animals rely on different decision strategies when faced with ambiguous or uncertain cues. Depending on the context, decisions may be biased towards events that were most frequently experienced in the past, or be more explorative. A particular type of decision making central to cognition is sequential memory recall in response to ambiguous cues. A previously developed spiking neuronal network implementation of sequence prediction and recall learns complex, high-order sequences in an unsupervised manner by local, biologically inspired plasticity rules. In response to an ambiguous cue, the model deterministically recalls the sequence shown most frequently during training. Here, we present an extension of the model enabling a range of different decision strategies. In this model, explorative behavior is generated by supplying neurons with noise. As the model relies on population encoding, uncorrelated noise averages out, and the recall dynamics remain effectively deterministic. In the presence of locally correlated noise, the averaging effect is avoided without impairing the model performance, and without the need for large noise amplitudes. We investigate two forms of correlated noise occurring in nature: shared synaptic background inputs, and random locking of the stimulus to spatiotemporal oscillations in the network activity. Depending on the noise characteristics, the network adopts various recall strategies. This study thereby provides potential mechanisms explaining how the statistics of learned sequences affect decision making, and how decision strategies can be adjusted after learning.

Characterization of spatiotemporal dynamics in EEG data during picture naming with optical flow patterns

In this study, we investigate the spatiotemporal dynamics of the neural oscillations by analyzing the electric potential that arises from neural activity. We identify two types of dynamics based on the frequency and phase of oscillations: standing waves or as out-of-phase and modulated waves, which represent a combination of standing and moving waves. To characterize these dynamics, we use optical flow patterns such as sources, sinks, spirals and saddles. We compare analytical and numerical solutions with real EEG data acquired during a picture-naming task. Analytical approximation of standing waves helps us to establish some properties of pattern location and number. Specifically, sources and sinks are mainly located in the same location, while saddles are positioned between them. The number of saddles correlates with the sum of all the other patterns. These properties are confirmed in both the simulated and real EEG data. In particular, source and sink clusters in the EEG data overlap with each other with median percentages around 60%, and hence have high spatial correlation, while source/sink clusters overlap with saddle clusters in less than 1%, and have different locations. Our statistical analysis showed that saddles account for about 45% of all patterns, while the remaining patterns are present in similar proportions.

Development of top-down cortical propagations in youth

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Dynamic cortical and tractography atlases of proactive and reactive alpha and high-gamma activities

Alpha waves-posterior dominant rhythms at 8-12 Hz reactive to eye opening and closure-are among the most fundamental EEG findings in clinical practice and research since Hans Berger first documented them in the early 20th century. Yet, the exact network dynamics of alpha waves in regard to eye movements remains unknown. High-gamma activity at 70-110 Hz is also reactive to eye movements and a summary measure of local cortical activation supporting sensorimotor or cognitive function. We aimed to build the first-ever brain atlases directly visualizing the network dynamics of eye movement-related alpha and high-gamma modulations, at cortical and white matter levels. We studied 28 patients (age: 5-20 years) who underwent intracranial EEG and electro-oculography recordings. We measured alpha and high-gamma modulations at 2167 electrode sites outside the seizure onset zone, interictal spike-generating areas and MRI-visible structural lesions. Dynamic tractography animated white matter streamlines modulated significantly and simultaneously beyond chance, on a millisecond scale. Before eye-closure onset, significant alpha augmentation occurred at the occipital and frontal cortices. After eye-closure onset, alpha-based functional connectivity was strengthened, while high gamma-based connectivity was weakened extensively in both intra-hemispheric and inter-hemispheric pathways involving the central visual areas. The inferior fronto-occipital fasciculus supported the strengthened alpha co-augmentation-based functional connectivity between occipital and frontal lobe regions, whereas the posterior corpus callosum supported the inter-hemispheric functional connectivity between the occipital lobes. After eye-opening offset, significant high-gamma augmentation and alpha attenuation occurred at occipital, fusiform and inferior parietal cortices. High gamma co-augmentation-based functional connectivity was strengthened, whereas alpha-based connectivity was weakened in the posterior inter-hemispheric and intra-hemispheric white matter pathways involving central and peripheral visual areas. Our results do not support the notion that eye closure-related alpha augmentation uniformly reflects feedforward or feedback rhythms propagating from lower to higher order visual cortex, or vice versa. Rather, proactive and reactive alpha waves involve extensive, distinct white matter networks that include the frontal lobe cortices, along with low- and high-order visual areas. High-gamma co-attenuation coupled to alpha co-augmentation in shared brain circuitry after eye closure supports the notion of an idling role for alpha waves during eye closure. These normative dynamic tractography atlases may improve understanding of the significance of EEG alpha waves in assessing the functional integrity of brain networks in clinical practice; they also may help elucidate the effects of eye movements on task-related brain network measures observed in cognitive neuroscience research.

Subsets of cortico-cortical evoked potentials propagate as traveling waves

Emerging evidence suggests that the temporal dynamics of cortico-cortical evoked potentials (CCEPs) may be used to characterize the patterns of information flow between and within brain networks. At present, however, the spatiotemporal dynamics of CCEP propagation cortically and subcortically are incompletely understood. We hypothesized that CCEPs propagate as an evoked traveling wave emanating from the site of stimulation. To elicit CCEPs, we applied single-pulse stimulation to stereoelectroencephalography (SEEG) electrodes implanted in 21 adult patients with intractable epilepsy. For each robust CCEP, we measured the timing of the maximal descent in evoked local field potentials and broadband high-gamma power (70-150 Hz) envelopes relative to the distance between the recording and stimulation contacts using three different metrics (i.e., Euclidean distance, path length, geodesic distance), representing direct, subcortical, and transcortical propagation, respectively. Many evoked responses to single-pulse electrical stimulation appear to propagate as traveling waves (~17-30%), even in the sparsely sampled, three-dimensional SEEG space. These results provide new insights into the spatiotemporal dynamics of CCEP propagation.

A deep generative adversarial network capturing complex spiral waves in disinhibited circuits of the cerebral cortex

BACKGROUND

In the cerebral cortex, disinhibited activity is characterized by propagating waves that spread across neural tissue. In this pathological state, a widely reported form of activity are spiral waves that travel in a circular pattern around a fixed spatial locus termed the center of mass. Spiral waves exhibit stereotypical activity and involve broad patterns of co-fluctuations, suggesting that they may be of lower complexity than healthy activity.

RESULTS

To evaluate this hypothesis, we performed dense multi-electrode recordings of cortical networks where disinhibition was induced by perfusing a pro-epileptiform solution containing 4-Aminopyridine as well as increased potassium and decreased magnesium. Spiral waves were identified based on a spatially delimited center of mass and a broad distribution of instantaneous phases across electrodes. Individual waves were decomposed into "snapshots" that captured instantaneous neural activation across the entire network. The complexity of these snapshots was examined using a measure termed the participation ratio. Contrary to our expectations, an eigenspectrum analysis of these snapshots revealed a broad distribution of eigenvalues and an increase in complexity compared to baseline networks. A deep generative adversarial network was trained to generate novel exemplars of snapshots that closely captured cortical spiral waves. These synthetic waves replicated key features of experimental data including a tight center of mass, a broad eigenvalue distribution, spatially-dependent correlations, and a high complexity. By adjusting the input to the model, new samples were generated that deviated in systematic ways from the experimental data, thus allowing the exploration of a broad range of states from healthy to pathologically disinhibited neural networks.

CONCLUSIONS

Together, results show that the complexity of population activity serves as a marker along a continuum from healthy to disinhibited brain states. The proposed generative adversarial network opens avenues for replicating the dynamics of cortical seizures and accelerating the design of optimal neurostimulation aimed at suppressing pathological brain activity.

Complexity of cortical wave patterns of the wake mouse cortex

Rich spatiotemporal dynamics of cortical activity, including complex and diverse wave patterns, have been identified during unconscious and conscious brain states. Yet, how these activity patterns emerge across different levels of wakefulness remain unclear. Here we study the evolution of wave patterns utilizing data from high spatiotemporal resolution optical voltage imaging of mice transitioning from barbiturate-induced anesthesia to wakefulness (N = 5) and awake mice (N = 4). We find that, as the brain transitions into wakefulness, there is a reduction in hemisphere-scale voltage waves, and an increase in local wave events and complexity. A neural mass model recapitulates the essential cellular-level features and shows how the dynamical competition between global and local spatiotemporal patterns and long-range connections can explain the experimental observations. These mechanisms possibly endow the awake cortex with enhanced integrative processing capabilities.

Synaptic propagation in neuronal networks with finite-support space-dependent coupling

Traveling waves of electrical activity are ubiquitous in biological neuronal networks. Traveling waves in the brain are associated with sensory processing, phase coding, and sleep. The neuron and network parameters that determine traveling waves' evolution are the synaptic space constant, synaptic conductance, membrane time constant, and synaptic decay time constant. We used an abstract neuron model in a one-dimensional network to investigate the propagation characteristics of traveling wave activity. We formulate a set of evolution equations based on the network connectivity parameters. Using a combination of numerical and analytical approaches, we show that these traveling waves are stable to a series of perturbations with biological relevance.

Cortex-wide spontaneous activity non-linearly steers propagating sensory-evoked activity in awake mice

The brain responds highly variably to identical sensory inputs, but there is no consensus on the nature of this variability. We explore this question using cortex-wide optical voltage imaging and whisker stimulation in awake mice. Clustering analysis reveals that the sensory-evoked activity propagates over the cortex via distinct pathways associated with distinct behavioral states. The pathway taken by each trial is independent of the level of primary sensory-evoked activation but is partially predictable by the spatiotemporal features of the preceding cortical spontaneous activity patterns. The sensory inputs reduce trial-to-trial variability in brain activity and alter temporal autocorrelation in spatial activity pattern evolutions, suggesting non-linear interactions between evoked activities and spontaneous activities. Further, evoked activities and spontaneous activities occupy different positions in the state space, suggesting that sensory inputs can intricately interact with the internal state to generate large-scale evoked activity patterns not frequented by spontaneous brain states.

High resolution, wide field optical imaging of macaque visual cortex with a curved detector

Objective. Cortical activity can be recorded using a variety of tools, ranging in scale from the single neuron (microscopic) to the whole brain (macroscopic). There is usually a trade-off between scale and resolution; optical imaging techniques, with their high spatio-temporal resolution and wide field of view, are best suited to study brain activity at the mesoscale. Optical imaging of cortical areas is however in practice limited by the curvature of the brain, which causes the image quality to deteriorate significantly away from the center of the image.Approach. To address this issue and harness the full potential of optical cortical imaging techniques, we developed a new wide-field optical imaging system adapted to the macaque brain. Our system is composed of a curved detector, an aspherical lens and a ring composed of light emitting diodes providing uniform illumination at wavelengths relevant for the different optical imaging methods, including intrinsic and fluorescence imaging.Main results. The system was characterized and compared with the standard macroscope used for cortical imaging, and a three-fold increase of the area in focus was measured as well as a four-fold increase in the evenness of the optical qualityin vivo.Significance. This new instrument, which is to the best of our knowledge the first use of a curved detector for cortical imaging, should facilitate the observation of wide mesoscale phenomena such as dynamic propagating waves within and between cortical maps, which are otherwise difficult to observe due to technical limitations of the currently available recording tools.

Cortical circuit-based lossless neural integrator for perceptual decision-making: A computational modeling study

The intrinsic uncertainty of sensory information (i.e., evidence) does not necessarily deter an observer from making a reliable decision. Indeed, uncertainty can be reduced by integrating (accumulating) incoming sensory evidence. It is widely thought that this accumulation is instantiated via recurrent rate-code neural networks. Yet, these networks do not fully explain important aspects of perceptual decision-making, such as a subject’s ability to retain accumulated evidence during temporal gaps in the sensory evidence. Here, we utilized computational models to show that cortical circuits can switch flexibly between “retention” and “integration” modes during perceptual decision-making. Further, we found that, depending on how the sensory evidence was readout, we could simulate “stepping” and “ramping” activity patterns, which may be analogous to those seen in different studies of decision-making in the primate parietal cortex. This finding may reconcile these previous empirical studies because it suggests these two activity patterns emerge from the same mechanism.

Learning populations with hubs govern the initiation and propagation of spontaneous bursts in neuronal networks after learning

Spontaneous bursts in neuronal networks with propagation involving a large number of synchronously firing neurons are considered to be a crucial feature of these networks both in vivo and in vitro. Recently, learning has been shown to improve the association and synchronization of spontaneous events in neuronal networks by promoting the firing of spontaneous bursts. However, little is known about the relationship between the learning phase and spontaneous bursts. By combining high-resolution measurement with a 4,096-channel complementary metal-oxide-semiconductor (CMOS) microelectrode array (MEA) and graph theory, we studied how the learning phase influenced the initiation of spontaneous bursts in cultured networks of rat cortical neurons in vitro. We found that a small number of selected populations carried most of the stimulus information and contributed to learning. Moreover, several new burst propagation patterns appeared in spontaneous firing after learning. Importantly, these "learning populations" had more hubs in the functional network that governed the initiation of spontaneous burst activity. These results suggest that changes in the functional structure of learning populations may be the key mechanism underlying increased bursts after learning. Our findings could increase understanding of the important role that synaptic plasticity plays in the regulation of spontaneous activity.

Visual evoked feedforward-feedback traveling waves organize neural activity across the cortical hierarchy in mice

Sensory processing is distributed among many brain regions that interact via feedforward and feedback signaling. Neuronal oscillations have been shown to mediate intercortical feedforward and feedback interactions. Yet, the macroscopic structure of the multitude of such oscillations remains unclear. Here, we show that simple visual stimuli reliably evoke two traveling waves with spatial wavelengths that cover much of the cerebral hemisphere in awake mice. 30-50 Hz feedforward waves arise in primary visual cortex (V1) and propagate rostrally, while 3-6 Hz feedback waves originate in the association cortex and flow caudally. The phase of the feedback wave modulates the amplitude of the feedforward wave and synchronizes firing between V1 and parietal cortex. Altogether, these results provide direct experimental evidence that visual evoked traveling waves percolate through the cerebral cortex and coordinate neuronal activity across broadly distributed networks mediating visual processing.

Propagation Structure of Intrinsic Brain Activity in Migraine without Aura

Previous studies have revealed highly reproducible patterns of temporally lagged brain activity in healthy human adults. However, it is unknown whether temporal organization of intrinsic activity is altered in migraines or if it relates to migraine chronification. In this resting-state functional magnetic resonance imaging study, temporal features of intrinsic activity were investigated using resting-state lag analysis, and 39 episodic migraine patients, 17 chronic migraine patients, and 35 healthy controls were assessed. Temporally earlier intrinsic activity in the hippocampal complex was revealed in the chronic migraine group relative to the other two groups. We also found earlier intrinsic activity in the medial prefrontal cortex in chronic compared with episodic migraines. Both migraine groups showed earlier intrinsic activity in the lateral temporal cortex and sensorimotor cortex compared with the healthy control group. Across all patients, headache frequency negatively correlated with temporal lag of the medial prefrontal cortex and hippocampal complex. Disrupted propagation of intrinsic activity in regions involved in sensory, cognitive and affective processing of pain may contribute to abnormal brain function during migraines. Decreased time latency in the lateral temporal cortex and sensorimotor cortex may be common manifestations in episodic and chronic migraines. The temporal features of the medial prefrontal cortex and hippocampal complex were associated with migraine chronification.

Horizontal connectivity in V1: Prediction of coherence in contour and motion integration

This study demonstrates the functional importance of the Surround context relayed laterally in V1 by the horizontal connectivity, in controlling the latency and the gain of the cortical response to the feedforward visual drive. We report here four main findings: 1) a centripetal apparent motion sequence results in a shortening of the spiking latency of V1 cells, when the orientation of the local inducer and the global motion axis are both co-aligned with the RF orientation preference; 2) this contextual effects grows with visual flow speed, peaking at 150–250°/s when it matches the propagation speed of horizontal connectivity (0.15–0.25 mm/ms); 3) For this speed range, the axial sensitivity of V1 cells is tilted by 90° to become co-aligned with the orientation preference axis; 4) the strength of modulation by the surround context correlates with the spatiotemporal coherence of the apparent motion flow. Our results suggest an internally-generated binding process, linking local (orientation /position) and global (motion/direction) features as early as V1. This long-range diffusion process constitutes a plausible substrate in V1 of the human psychophysical bias in speed estimation for collinear motion. Since it is demonstrated in the anesthetized cat, this novel form of contextual control of the cortical gain and phase is a built-in property in V1, whose expression does not require behavioral attention and top-down control from higher cortical areas. We propose that horizontal connectivity participates in the propagation of an internal “prediction” wave, shaped by visual experience, which links contour co-alignment and global axial motion at an apparent speed in the range of saccade-like eye movements.

A hybrid biological neural network model for solving problems in cognitive planning

A variety of behaviors, like spatial navigation or bodily motion, can be formulated as graph traversal problems through cognitive maps. We present a neural network model which can solve such tasks and is compatible with a broad range of empirical findings about the mammalian neocortex and hippocampus. The neurons and synaptic connections in the model represent structures that can result from self-organization into a cognitive map via Hebbian learning, i.e. into a graph in which each neuron represents a point of some abstract task-relevant manifold and the recurrent connections encode a distance metric on the manifold. Graph traversal problems are solved by wave-like activation patterns which travel through the recurrent network and guide a localized peak of activity onto a path from some starting position to a target state.

Event-Related Phase Synchronization Propagates Rapidly across Human Ventral Visual Cortex

Visual inputs to early visual cortex integrate with semantic, linguistic and memory inputs in higher visual cortex, in a manner that is rapid and accurate, and enables complex computations such as face recognition and word reading. This implies the existence of fundamental organizational principles that enable such efficiency. To elaborate on this, we performed intracranial recordings in 82 individuals while they performed tasks of varying visual and cognitive complexity. We discovered that visual inputs induce highly organized posterior-to-anterior propagating patterns of phase modulation across the ventral occipitotemporal cortex. At individual electrodes there was a stereotyped temporal pattern of phase progression following both stimulus onset and offset, consistent across trials and tasks. The phase of low frequency activity in anterior regions was predicted by the prior phase in posterior cortical regions. This spatiotemporal propagation of phase likely serves as a feed-forward organizational influence enabling the integration of information across the ventral visual stream. This phase modulation manifests as the early components of the event related potential; one of the most commonly used measures in human electrophysiology. These findings illuminate fundamental organizational principles of the higher order visual system that enable the rapid recognition and characterization of a variety of inputs.

Propofol Anesthesia Alters Cortical Traveling Waves

Oscillatory dynamics in cortex seem to organize into traveling waves that serve a variety of functions. Recent studies show that propofol, a widely used anesthetic, dramatically alters cortical oscillations by increasing slow-delta oscillatory power and coherence. It is not known how this affects traveling waves. We compared traveling waves across the cortex of non-human primates before, during, and after propofol-induced loss of consciousness (LOC). After LOC, traveling waves in the slow-delta (∼1 Hz) range increased, grew more organized, and traveled in different directions relative to the awake state. Higher frequency (8-30 Hz) traveling waves, by contrast, decreased, lost structure, and switched to directions where the slow-delta waves were less frequent. The results suggest that LOC may be due, in part, to increases in the strength and direction of slow-delta traveling waves that, in turn, alter and disrupt traveling waves in the higher frequencies associated with cognition.

Individual resting-state alpha peak frequency and within-trial changes in alpha peak frequency both predict visual dual-pulse segregation performance

Although sensory input is continuous, information must be combined over time to guide action and cognition, leading to the proposal of temporal sampling windows. A number of studies have suggested that a 10-Hz sampling window might be involved in the "frame rate" of visual processing. To investigate this, we tested the ability of participants to localize and enumerate 1 or 2 visual flashes presented either at near-threshold or full-contrast intensities, while recording magnetoencephalography. The inter-stimulus interval (ISI) between the 2 flashes was varied across trials. Performance in distinguishing between 1 and 2 flashes was linked to the alpha frequency, both at the individual level and trial-by-trial. Participants with a higher resting-state alpha peak frequency showed the greatest improvement in performance as a function of ISI within a 100-ms time window, while those with slower alpha improved more when ISI exceeded 100 ms. On each trial, correct enumeration (1 vs. 2) performance was paired with faster pre-stimulus instantaneous alpha frequency. Our results suggest that visual sampling/processing speed, linked to peak alpha frequency, is both an individual trait and can vary in a state-dependent manner.

Cortical traveling waves reflect state-dependent hierarchical sequencing of local regions in the human connectome network

Recent human studies using electrocorticography have demonstrated that alpha and theta band oscillations form traveling waves on the cortical surface. According to neural synchronization theories, the cortical traveling waves may group local cortical regions and sequence them by phase synchronization; however these contributions have not yet been assessed. This study aimed to evaluate the functional contributions of traveling waves using connectome-based network modeling. In the simulation, we observed stable traveling waves on the entire cortical surface wherein the topographical pattern of these phases was substantially correlated with the empirically obtained resting-state networks, and local radial waves also appeared within the size of the empirical networks (< 50 mm). Importantly, individual regions in the entire network were instantaneously sequenced by their internal frequencies, and regions with higher intrinsic frequency were seen in the earlier phases of the traveling waves. Based on the communication-through-coherence theory, this phase configuration produced a hierarchical organization of each region by unidirectional communication between the arbitrarily paired regions. In conclusion, cortical traveling waves reflect the intrinsic frequency-dependent hierarchical sequencing of local regions, global traveling waves sequence the set of large-scale cortical networks, and local traveling waves sequence local regions within individual cortical networks.

Magnifying Traveling Waves on the Scalp

Traveling waves appear in various signals that measure neuronal activity. Some signals measured in animals have singles-cell resolution and directly point to neuronal activity. In those cases, activation of distributed neurons forms a wave front, and the front propagates across the cortical surface. Other signals are variants of neuroelectric potentials, i.e. electroencephalography, electrocorticography and field potentials. Instead of having fine spatial resolution, these signals reflect the activity of neuronal populations via volume conduction (VC). Sources of traveling waves in neuroelectric potentials have not been well addressed so far. As animal studies show propagating activation of neurons that spread in measured areas, it is often considered that neuronal activations during scalp waves have similar trajectories of activation, spreading like scalp waves. However, traveling waves on the scalp differ from those found directly on the cortical surface in several dimensions: traveling velocity, traveling distance and areal size occupied by single polarity. We describe that the simplest sources can produce scalp waves with perceived spatial dimensions which are actually a magnification of neuronal activity emanating from local sources due to VC. This viewpoint is not a rigorous proof of our magnification concept. However, we suggest the possibility that the actual dimensions of neuronal activity producing traveling waves is not as large as the dimension of the traveling waves.

Methods for inferring neural circuit interactions and neuromodulation from local field potential and electroencephalogram measures

Electrical recordings of neural mass activity, such as local field potentials (LFPs) and electroencephalograms (EEGs), have been instrumental in studying brain function. However, these aggregate signals lack cellular resolution and thus are not easy to be interpreted directly in terms of parameters of neural microcircuits. Developing tools for a reliable estimation of key neural parameters from these signals, such as the interaction between excitation and inhibition or the level of neuromodulation, is important for both neuroscientific and clinical applications. Over the years, we have developed tools based on neural network modeling and computational analysis of empirical data to estimate neural parameters from aggregate neural signals. This review article gives an overview of the main computational tools that we have developed and employed to invert LFPs and EEGs in terms of circuit-level neural phenomena, and outlines future challenges and directions for future research.

Spatiotemporal structure of sensory-evoked and spontaneous activity revealed by mesoscale imaging in anesthetized and awake mice

Stimuli-evoked and spontaneous brain activity propagates across the cortex in diverse spatiotemporal patterns. Despite extensive studies, the relationship between spontaneous and evoked activity is poorly understood. We investigate this relationship by comparing the amplitude, speed, direction, and complexity of propagation trajectories of spontaneous and evoked activity elicited with visual, auditory, and tactile stimuli using mesoscale wide-field imaging in mice. For both spontaneous and evoked activity, the speed and direction of propagation is modulated by the amplitude. However, spontaneous activity has a higher complexity of the propagation trajectories. For low stimulus strengths, evoked activity amplitude and speed is similar to that of spontaneous activity but becomes dissimilar at higher stimulus strengths. These findings are consistent with observations that primary sensory areas receive widespread inputs from other cortical regions, and during rest, the cortex tends to reactivate traces of complex multisensory experiences that might have occurred in exhibition of different behaviors.

Traveling Waves in Quasi-One-Dimensional Neuronal Minicolumns

Traveling waves of neuronal activity in the cortex have been observed in vivo. These traveling waves have been correlated to various features of observed cortical dynamics, including spike timing variability and correlated fluctuations in neuron membrane potential. Although traveling waves are typically studied as either strictly one-dimensional or two-dimensional excitations, here we investigate the conditions for the existence of quasi-one-dimensional traveling waves that could be sustainable in parts of the brain containing cortical minicolumns. For that, we explore a quasi-one-dimensional network of heterogeneous neurons with a biologically influenced computational model of neuron dynamics and connectivity. We find that background stimulus reliably evokes traveling waves in networks with local connectivity between neurons. We also observe traveling waves in fully connected networks when a model for action potential propagation speed is incorporated. The biological properties of the neurons influence the generation and propagation of the traveling waves. Our quasi-one-dimensional model is not only useful for studying the basic properties of traveling waves in neuronal networks; it also provides a simplified representation of possible wave propagation in columnar or minicolumnar networks found in the cortex.

Spontaneous traveling waves naturally emerge from horizontal fiber time delays and travel through locally asynchronous-irregular states

Studies of sensory-evoked neuronal responses often focus on mean spike rates, with fluctuations treated as internally-generated noise. However, fluctuations of spontaneous activity, often organized as traveling waves, shape stimulus-evoked responses and perceptual sensitivity. The mechanisms underlying these waves are unknown. Further, it is unclear whether waves are consistent with the low rate and weakly correlated “asynchronous-irregular” dynamics observed in cortical recordings. Here, we describe a large-scale computational model with topographically-organized connectivity and conduction delays relevant to biological scales. We find that spontaneous traveling waves are a general property of these networks. The traveling waves that occur in the model are sparse, with only a small fraction of neurons participating in any individual wave. Consequently, they do not induce measurable spike correlations and remain consistent with locally asynchronous irregular states. Further, by modulating local network state, they can shape responses to incoming inputs as observed in vivo.

Spontaneous traveling cortical waves shape neural responses. Using a large-scale computational model, the authors show that transmission delays shape locally asynchronous spiking dynamics into traveling waves without inducing correlations and boost responses to external input, as observed in vivo.

Topological Receptive Field Model for Human Retinotopic Mapping

The mapping between visual inputs on the retina and neuronal activations in the visual cortex, i.e., retinotopic map, is an essential topic in vision science and neuroscience. Human retinotopic maps can be revealed by analyzing the functional magnetic resonance imaging (fMRI) signal responses to designed visual stimuli in vivo. Neurophysiology studies summarized that visual areas are topological (i.e., nearby neurons have receptive fields at nearby locations in the image). However, conventional fMRI-based analyses frequently generate non-topological results because they process fMRI signals on a voxel-wise basis, without considering the neighbor relations on the surface. Here we propose a topological receptive field (tRF) model which imposes the topological condition when decoding retinotopic fMRI signals. More specifically, we parametrized the cortical surface to a unit disk, characterized the topological condition by tRF, and employed an efficient scheme to solve the tRF model. We tested our framework on both synthetic and human fMRI data. Experimental results showed that the tRF model could remove the topological violations, improve model explaining power, and generate biologically plausible retinotopic maps. The proposed framework is general and can be applied to other sensory maps.

Cross Laminar Traveling Components of Field Potentials due to Volume Conduction of Non-Traveling Neuronal Activity in Macaque Sensory Cortices

Field potentials (FPs) reflect neuronal activities in the brain, and often exhibit traveling peaks across recording sites. While traveling FPs are interpreted as propagation of neuronal activity, not all studies directly reveal such propagating patterns of neuronal activation. Neuronal activity is associated with transmembrane currents that form dipoles and produce negative and positive fields. Thereby, FP components reverse polarity between those fields and have minimal amplitudes at the center of dipoles. Although their amplitudes could be smaller, FPs are never flat even around these reversals. What occurs around the reversal has not been addressed explicitly, although those are rationally in the middle of active neurons. We show that sensory FPs around the reversal appeared with peaks traveling across cortical laminae in macaque sensory cortices. Interestingly, analyses of current source density did not depict traveling patterns but lamina-delimited current sinks and sources. We simulated FPs produced by volume conduction of a simplified 2 dipoles' model mimicking sensory cortical laminar current source density components. While FPs generated by single dipoles followed the temporal patterns of the dipole moments without traveling peaks, FPs generated by concurrently active dipole moments appeared with traveling components in the vicinity of dipoles by superimposition of individually non-traveling FPs generated by single dipoles. These results indicate that not all traveling FP are generated by traveling neuronal activity, and that recording positions need to be taken into account to describe FP peak components around active neuronal populations.SIGNIFICANCE STATEMENT Field potentials (FPs) generated by neuronal activity in the brain occur with fields of opposite polarity. Likewise, in the cerebral cortices, they have mirror-imaged waveforms in upper and lower layers. We show that FPs appear like traveling across the cortical layers. Interestingly, the traveling FPs occur without traveling components of current source density, which represents transmembrane currents associated with neuronal activity. These seemingly odd findings are explained using current source density models of multiple dipoles. Concurrently active, non-traveling dipoles produce FPs as mixtures of FPs produced by individual dipoles, and result in traveling FP waveforms as the mixing ratio depends on the distances from those dipoles. The results suggest that not all traveling FP components are associated with propagating neuronal activity.

Go with the FLOW: visualizing spatiotemporal dynamics in optical widefield calcium imaging

Widefield calcium imaging has recently emerged as a powerful experimental technique to record coordinated large-scale brain activity. These measurements present a unique opportunity to characterize spatiotemporally coherent structures that underlie neural activity across many regions of the brain. In this work, we leverage analytic techniques from fluid dynamics to develop a visualization framework that highlights features of flow across the cortex, mapping wavefronts that may be correlated with behavioural events. First, we transform the time series of widefield calcium images into time-varying vector fields using optic flow. Next, we extract concise diagrams summarizing the dynamics, which we refer to as FLOW (flow lines in optical widefield imaging) portraits. These FLOW portraits provide an intuitive map of dynamic calcium activity, including regions of initiation and termination, as well as the direction and extent of activity spread. To extract these structures, we use the finite-time Lyapunov exponent technique developed to analyse time-varying manifolds in unsteady fluids. Importantly, our approach captures coherent structures that are poorly represented by traditional modal decomposition techniques. We demonstrate the application of FLOW portraits on three simple synthetic datasets and two widefield calcium imaging datasets, including cortical waves in the developing mouse and spontaneous cortical activity in an adult mouse.

The impact of a closed-loop thalamocortical model on the spatiotemporal dynamics of cortical and thalamic traveling waves

Propagation of activity in spatially structured neuronal networks has been observed in awake, anesthetized, and sleeping brains. How these wave patterns emerge and organize across brain structures, and how network connectivity affects spatiotemporal neural activity remains unclear. Here, we develop a computational model of a two-dimensional thalamocortical network, which gives rise to emergent traveling waves similar to those observed experimentally. We illustrate how spontaneous and evoked oscillatory activity in space and time emerge using a closed-loop thalamocortical architecture, sustaining smooth waves in the cortex and staggered waves in the thalamus. We further show that intracortical and thalamocortical network connectivity, cortical excitation/inhibition balance, and thalamocortical or corticothalamic delay can independently or jointly change the spatiotemporal patterns (radial, planar and rotating waves) and characteristics (speed, direction, and frequency) of cortical and thalamic traveling waves. Computer simulations predict that increased thalamic inhibition induces slower cortical frequencies and that enhanced cortical excitation increases traveling wave speed and frequency. Overall, our results provide insight into the genesis and sustainability of thalamocortical spatiotemporal patterns, showing how simple synaptic alterations cause varied spontaneous and evoked wave patterns. Our model and simulations highlight the need for spatially spread neural recordings to uncover critical circuit mechanisms for brain functions.

Enhancing visual motion discrimination by desynchronizing bifocal oscillatory activity

Visual motion discrimination involves reciprocal interactions in the alpha band between the primary visual cortex (V1) and mediotemporal areas (V5/MT). We investigated whether modulating alpha phase synchronization using individualized multisite transcranial alternating current stimulation (tACS) over V5 and V1 regions would improve motion discrimination. We tested 3 groups of healthy subjects with the following conditions: (1) individualized In-Phase V1_alpha-V5_alpha tACS (0° lag), (2) individualized Anti-Phase V1_alpha-V5_alpha tACS (180° lag) and (3) sham tACS. Motion discrimination and EEG activity were recorded before, during and after tACS. Performance significantly improved in the Anti-Phase group compared to the In-Phase group 10 and 30 min after stimulation. This result was explained by decreases in bottom-up alpha-V1 gamma-V5 phase-amplitude coupling. One possible explanation of these results is that Anti-Phase V1_alpha-V5_alpha tACS might impose an optimal phase lag between stimulation sites due to the inherent speed of wave propagation, hereby supporting optimized neuronal communication.

Mean-Field Models for EEG/MEG: From Oscillations to Waves

Neural mass models have been used since the 1970s to model the coarse-grained activity of large populations of neurons. They have proven especially fruitful for understanding brain rhythms. However, although motivated by neurobiological considerations they are phenomenological in nature, and cannot hope to recreate some of the rich repertoire of responses seen in real neuronal tissue. Here we consider a simple spiking neuron network model that has recently been shown to admit an exact mean-field description for both synaptic and gap-junction interactions. The mean-field model takes a similar form to a standard neural mass model, with an additional dynamical equation to describe the evolution of within-population synchrony. As well as reviewing the origins of this next generation mass model we discuss its extension to describe an idealised spatially extended planar cortex. To emphasise the usefulness of this model for EEG/MEG modelling we show how it can be used to uncover the role of local gap-junction coupling in shaping large scale synaptic waves.

Learning poly-synaptic paths with traveling waves

Traveling waves are commonly observed across the brain. While previous studies have suggested the role of traveling waves in learning, the mechanism remains unclear. We adopted a computational approach to investigate the effect of traveling waves on synaptic plasticity. Our results indicate that traveling waves facilitate the learning of poly-synaptic network paths when combined with a reward-dependent local synaptic plasticity rule. We also demonstrate that traveling waves expedite finding the shortest paths and learning nonlinear input/output mapping, such as exclusive or (XOR) function.

There are approximately 10¹¹ neurons with 10¹⁴ connections in the human brain. Information transmission among neurons in this large network is considered crucial for our behavior. To achieve this, multiple synaptic connections along a poly-synaptic network path must be adjusted coherently during learning. Because the previously proposed reward-dependent synaptic plasticity rule requires coactivation of presynaptic and postsynaptic neurons, learning can fail if a subset of neurons along a distant network path is inactive at the beginning of learning. We suggest that traveling waves that are initiated at an information source can mitigate this problem. We performed computer simulations of spiking neural networks with reward-dependent local synaptic plasticity rules and traveling waves. Our results show that this combination facilitates the learning and refinement of synaptic network paths. We argue that these features are a general biological strategy for maintaining and optimizing our brain function. Our research provides new insights into how complex neural networks in the brain form during learning and memory consolidation.

Turning the Stimulus On and Off Changes the Direction of α Traveling Waves

Traveling waves have been studied to characterize the complex spatiotemporal dynamics of the brain. Several studies have suggested that the propagation direction of α traveling waves can be task dependent. For example, a recent electroencephalography (EEG) study from our group found that forward waves (i.e., occipital to frontal, FW waves) were observed during visual processing, whereas backward waves (i.e., frontal to occipital, BW waves) mostly occurred in the absence of sensory input. These EEG recordings, however, were obtained from different experimental sessions and different groups of subjects. To further examine how the waves’ direction changes between task conditions, 13 human participants were tested on a target detection task while EEG signals were recorded simultaneously. We alternated visual stimulation (5-s display of visual luminance sequences) and resting state (5 s of black screen) within each single trial, allowing us to monitor the moment-to-moment progression of traveling waves. As expected, the direction of α waves was closely linked with task conditions. First, FW waves from occipital to frontal regions, absent during rest, emerged as a result of visual processing, while BW waves in the opposite direction dominated in the absence of visual inputs, and were reduced (but not eliminated) by external visual inputs. Second, during visual stimulation (but not rest), both waves coexisted on average, but were negatively correlated. In summary, we conclude that the functional role of α traveling waves is closely related with their propagating direction, with stimulus-evoked FW waves supporting visual processing and spontaneous BW waves involved more in top-down control.

DMT alters cortical travelling waves

Psychedelic drugs are potent modulators of conscious states and therefore powerful tools for investigating their neurobiology. N,N, Dimethyltryptamine (DMT) can rapidly induce an extremely immersive state of consciousness characterized by vivid and elaborate visual imagery. Here, we investigated the electrophysiological correlates of the DMT-induced altered state from a pool of participants receiving DMT and (separately) placebo (saline) while instructed to keep their eyes closed. Consistent with our hypotheses, results revealed a spatio-temporal pattern of cortical activation (i.e. travelling waves) similar to that elicited by visual stimulation. Moreover, the typical top-down alpha-band rhythms of closed-eyes rest were significantly decreased, while the bottom-up forward wave was significantly increased. These results support a recent model proposing that psychedelics reduce the 'precision-weighting of priors', thus altering the balance of top-down versus bottom-up information passing. The robust hypothesis-confirming nature of these findings imply the discovery of an important mechanistic principle underpinning psychedelic-induced altered states.

Optimizing Visual Cortex Parameterization with Error-Tolerant Teichmüller Map in Retinotopic Mapping

The mapping between the visual input on the retina to the cortical surface, i.e., retinotopic mapping, is an important topic in vision science and neuroscience. Human retinotopic mapping can be revealed by analyzing cortex functional magnetic resonance imaging (fMRI) signals when the subject is under specific visual stimuli. Conventional methods process, smooth, and analyze the retinotopic mapping based on the parametrization of the (partial) cortical surface. However, the retinotopic maps generated by this approach frequently contradict neuropsychology results. To address this problem, we propose an integrated approach that parameterizes the cortical surface, such that the parametric coordinates linearly relates the visual coordinate. The proposed method helps the smoothing of noisy retinotopic maps and obtains neurophysiological insights in human vision systems. One key element of the approach is the Error-Tolerant Teichmüller Map, which uniforms the angle distortion and maximizes the alignments to self-contradicting landmarks. We validated our overall approach with synthetic and real retinotopic mapping datasets. The experimental results show the proposed approach is superior in accuracy and compatibility. Although we focus on retinotopic mapping, the proposed framework is general and can be applied to process other human sensory maps.

Direction-selective motion discrimination by traveling waves in visual cortex

The majority of neurons in primary visual cortex respond selectively to bars of light that have a specific orientation and move in a specific direction. The spatial and temporal responses of such neurons are non-separable. How neurons accomplish that computational feat without resort to explicit time delays is unknown. We propose a novel neural mechanism whereby visual cortex computes non-separable responses by generating endogenous traveling waves of neural activity that resonate with the space-time signature of the visual stimulus. The spatiotemporal characteristics of the response are defined by the local topology of excitatory and inhibitory lateral connections in the cortex. We simulated the interaction between endogenous traveling waves and the visual stimulus using spatially distributed populations of excitatory and inhibitory neurons with Wilson-Cowan dynamics and inhibitory-surround coupling. Our model reliably detected visual gratings that moved with a given speed and direction provided that we incorporated neural competition to suppress false motion signals in the opposite direction. The findings suggest that endogenous traveling waves in visual cortex can impart direction-selectivity on neural responses without resort to explicit time delays. They also suggest a functional role for motion opponency in eliminating false motion signals.

It is well established that the so-called ‘simple cells’ of the primary visual cortex respond preferentially to oriented bars of light that move across the visual field with a particular speed and direction. The spatiotemporal responses of such neurons are said to be non-separable because they cannot be constructed from independent spatial and temporal neural mechanisms. Contemporary theories of how neurons compute non-separable responses typically rely on finely tuned transmission delays between signals from disparate regions of the visual field. However the existence of such delays is controversial. We propose an alternative neural mechanism for computing non-separable responses that does not require transmission delays. It instead relies on the predisposition of the cortical tissue to spontaneously generate spatiotemporal waves of neural activity that travel with a particular speed and direction. We propose that the endogenous wave activity resonates with the visual stimulus to elicit direction-selective neural responses to visual motion. We demonstrate the principle in computer models and show that competition between opposing neurons robustly enhances their ability to discriminate between visual gratings that move in opposite directions.