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

Mixed selectivity: Cellular computations for complexity

The property of mixed selectivity has been discussed at a computational level and offers a strategy to maximize computational power by adding versatility to the functional role of each neuron. Here, we offer a biologically grounded implementational-level mechanistic explanation for mixed selectivity in neural circuits. We define pure, linear, and nonlinear mixed selectivity and discuss how these response properties can be obtained in simple neural circuits. Neurons that respond to multiple, statistically independent variables display mixed selectivity. If their activity can be expressed as a weighted sum, then they exhibit linear mixed selectivity; otherwise, they exhibit nonlinear mixed selectivity. Neural representations based on diverse nonlinear mixed selectivity are high dimensional; hence, they confer enormous flexibility to a simple downstream readout neural circuit. However, a simple neural circuit cannot possibly encode all possible mixtures of variables simultaneously, as this would require a combinatorially large number of mixed selectivity neurons. Gating mechanisms like oscillations and neuromodulation can solve this problem by dynamically selecting which variables are mixed and transmitted to the readout.

Stereotyped spatiotemporal dynamics of spontaneous activity in visual cortex prior to eye-opening

Over the course of development, functional sensory representations emerge in the visual cortex. Prior to eye-opening, modular patterns of spontaneous activity form long-range networks that may serve as a precursor for mature network organization. Although the spatial structure of these networks has been well studied, their temporal features, which may contribute to their continued plasticity and development, remain largely uncharacterized. To address this, we imaged hours of spontaneous network activity in the visual cortex of developing ferrets of both sexes utilizing a fast calcium indicator (GCaMP8m) and widefield imaging at high temporal resolution (50Hz), then segmented out spatiotemporal events. The spatial structure of this activity was highly modular, exhibiting spatially segregated active domains consistent with prior work. We found that the vast majority of events showed a clear dynamic component in which modules activated sequentially across the field of view, but only a minority of events were well-fit with a linear traveling wave. We found that spatiotemporal events occur in repeated and stereotyped motifs, reoccurring across hours of imaging. Finally, we found that the most frequently occurring single-frame spatial activity patterns were predictive of future activity patterns over hundreds of milliseconds. Together, our results demonstrate that spontaneous activity in the early developing cortex exhibits a rich spatiotemporal structure, suggesting a potential role in the maturation and refinement of future functional representations.

Brain-state mediated modulation of inter-laminar dependencies in visual cortex

Spatial attention is critical for recognizing behaviorally relevant objects in a cluttered environment. How the deployment of spatial attention aids the hierarchical computations of object recognition remains unclear. We investigated this in the laminar cortical network of visual area V4, an area strongly modulated by attention. We found that deployment of attention strengthened unique dependencies in neural activity across cortical layers. On the other hand, shared dependencies were reduced within the excitatory population of a layer. Surprisingly, attention strengthened unique dependencies within a laminar population. Crucially, these modulation patterns were also observed during successful behavioral outcomes that are thought to be mediated by internal brain state fluctuations. Successful behavioral outcomes were also associated with phases of reduced neural excitability, suggesting a mechanism for enhanced information transfer during optimal states. Our results suggest common computation goals of optimal sensory states that are attained by either task demands or internal fluctuations.

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.

Simulated dynamical transitions in a heterogeneous marmoset pFC cluster

Network analysis of the marmoset cortical connectivity data indicates a significant 3D cluster in and around the pre-frontal cortex. A multi-node, heterogeneous neural mass model of this six-node cluster was constructed. Its parameters were informed by available experimental and simulation data so that each neural mass oscillated in a characteristic frequency band. Nodes were connected with directed, weighted links derived from the marmoset structural connectivity data. Heterogeneity arose from the different link weights and model parameters for each node. Stimulation of the cluster with an incident pulse train modulated in the standard frequency bands induced a variety of dynamical state transitions that lasted in the range of 5-10 s, suggestive of timescales relevant to short-term memory. A short gamma burst rapidly reset the beta-induced transition. The theta-induced transition state showed a spontaneous, delayed reset to the resting state. An additional, continuous gamma wave stimulus induced a new beating oscillatory state. Longer or repeated gamma bursts were phase-aligned with the beta oscillation, delivering increasing energy input and causing shorter transition times. The relevance of these results to working memory is yet to be established, but they suggest interesting opportunities.

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.

Human connectome topology directs cortical traveling waves and shapes frequency gradients

Traveling waves and neural oscillation frequency gradients are pervasive in the human cortex. While the direction of traveling waves has been linked to brain function and dysfunction, the factors that determine this direction remain elusive. We hypothesized that structural connectivity instrength gradients - defined as the gradually varying sum of incoming connection strengths across the cortex - could shape both traveling wave direction and frequency gradients. We confirm the presence of instrength gradients in the human connectome across diverse cohorts and parcellations. Using a cortical network model, we demonstrate how these instrength gradients direct traveling waves and shape frequency gradients. Our model fits resting-state MEG functional connectivity best in a regime where instrength-directed traveling waves and frequency gradients emerge. We further show how structural subnetworks of the human connectome generate opposing wave directions and frequency gradients observed in the alpha and beta bands. Our findings suggest that structural connectivity instrength gradients affect both traveling wave direction and frequency gradients.

Neural assemblies coordinated by cortical waves are associated with waking and hallucinatory brain states

The relationship between sensory stimuli and perceptions is brain-state dependent: in wakefulness, suprathreshold stimuli evoke perceptions; under anesthesia, perceptions are abolished; and during dreaming and in dissociated states, percepts are internally generated. Here, we exploit this state dependence to identify brain activity associated with internally generated or stimulus-evoked perceptions. In awake mice, visual stimuli phase reset spontaneous cortical waves to elicit 3-6 Hz feedback traveling waves. These stimulus-evoked waves traverse the cortex and entrain visual and parietal neurons. Under anesthesia as well as during ketamine-induced dissociation, visual stimuli do not disrupt spontaneous waves. Uniquely, in the dissociated state, spontaneous waves traverse the cortex caudally and entrain visual and parietal neurons, akin to stimulus-evoked waves in wakefulness. Thus, coordinated neuronal assemblies orchestrated by traveling cortical waves emerge in states in which perception can manifest. The awake state is privileged in that this coordination is reliably elicited by external visual stimuli.

Induction of excitatory brain state governs plastic functional changes in visual cortical topology

Adult visual plasticity underlying local remodeling of the cortical circuitry in vivo appears to be associated with a spatiotemporal pattern of strongly increased spontaneous and evoked activity of populations of cells. Here we review and discuss pioneering work by us and others about principles of plasticity in the adult visual cortex, starting with our study which showed that a confined lesion in the cat retina causes increased excitability in the affected region in the primary visual cortex accompanied by fine-tuned restructuring of neuronal function. The underlying remodeling processes was further visualized with voltage-sensitive dye (VSD) imaging that allowed a direct tracking of retinal lesion-induced reorganization across horizontal cortical circuitries. Nowadays, application of noninvasive stimulation methods pursues the idea further of increased cortical excitability along with decreased inhibition as key factors for the induction of adult cortical plasticity. We used high-frequency transcranial magnetic stimulation (TMS), for the first time in combination with VSD optical imaging, and provided evidence that TMS-amplified excitability across large pools of neurons forms the basis for noninvasively targeting reorganization of orientation maps in the visual cortex. Our review has been compiled on the basis of these four own studies, which we discuss in the context of historical developments in the field of visual cortical plasticity and the current state of the literature. Overall, we suggest markers of LTP-like cortical changes at mesoscopic population level as a main driving force for the induction of visual plasticity in the adult. Elevations in excitability that predispose towards cortical plasticity are most likely a common property of all cortical modalities. Thus, interventions that increase cortical excitability are a promising starting point to drive perceptual and potentially motor learning in therapeutic applications.

Common modular architecture across diverse cortical areas in early development

In order to deal with a complex environment, animals form a diverse range of neural representations that vary across cortical areas, ranging from largely unimodal sensory input to higher-order representations of goals, outcomes, and motivation. The developmental origin of this diversity is currently unclear, as representations could arise through processes that are already area-specific from the earliest developmental stages or alternatively, they could emerge from an initially common functional organization shared across areas. Here, we use spontaneous activity recorded with two-photon and widefield calcium imaging to reveal the functional organization across the early developing cortex in ferrets, a species with a well-characterized columnar organization and modular structure of spontaneous activity in the visual cortex. We find that in animals 7 to 14 d prior to eye-opening and ear canal opening, spontaneous activity in both sensory areas (auditory and somatosensory cortex, A1 and S1, respectively), and association areas (posterior parietal and prefrontal cortex, PPC and PFC, respectively) showed an organized and modular structure that is highly similar to the organization in V1. In all cortical areas, this modular activity was distributed across the cortical surface, forming functional networks that exhibit millimeter-scale correlations. Moreover, this modular structure was evident in highly coherent spontaneous activity at the cellular level, with strong correlations among local populations of neurons apparent in all cortical areas examined. Together, our results demonstrate a common distributed and modular organization across the cortex during early development, suggesting that diverse cortical representations develop initially according to similar design principles.

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.

Corrective mechanisms of motion extrapolation

Transmission and processing of sensory information in the visual system takes time. For motion perception, our brain can overcome this intrinsic neural delay through extrapolation mechanisms and accurately predict the current position of a continuously moving object. But how does the system behave when the motion abruptly changes and the prediction becomes wrong? Here we address this question by studying the perceived position of a moving object with various abrupt motion changes by human observers. We developed a task in which a bar is monotonously moving horizontally, and then motion suddenly stops, reverses, or disappears-then-reverses around two vertical stationary reference lines. Our results showed that participants overestimated the position of the stopping bar but did not perceive an overshoot in the motion reversal condition. When a temporal gap was added at the reverse point, the perceptual overshoot of the end point scaled with the gap durations. Our model suggests that the overestimation of the object position when it disappears is not linear as a function of its speeds but gradually fades out. These results can thus be reconciled in a single process where there is an interplay of the cortical motion prediction mechanisms and the late sensory transient visual inputs.

Horizontal cortical connections shape intrinsic traveling waves into feature-selective motifs that regulate perceptual sensitivity

Intrinsic, ongoing fluctuations of cortical activity form traveling waves that modulate the gain of sensory-evoked responses and perceptual sensitivity. Several lines of evidence suggest that intrinsic traveling waves (iTWs) may arise, in part, from the coordination of synaptic activity through the recurrent horizontal connectivity within cortical areas, which include long range patchy connections that link neurons with shared feature preferences. In a spiking network model with anatomical topology that incorporates feature-selective patchy connections, which we call the Balanced Patchy Network (BPN), we observe repeated iTWs, which we refer to as motifs. In the model, motifs stem from fluctuations in the excitability of like-tuned neurons that result from shifts in E/I balance as action potentials traverse these patchy connections. To test if feature-selective motifs occur in vivo, we examined data previously recorded using multielectrode arrays in Area MT of marmosets trained to perform a threshold visual detection task. Using a newly developed method for comparing the similarity of wave patterns we found that some iTWs can be grouped into motifs. As predicted by the BPN, many of these motifs are feature selective, exhibiting direction-selective modulations in ongoing spiking activity. Further, motifs modulate the gain of the response evoked by a target and perceptual sensitivity to the target if the target matches the preference of the motif. These results provide evidence that iTWs are shaped by the patterns of horizontal fiber projections in the cortex and that patchy connections enable iTWs to regulate neural and perceptual sensitivity in a feature selective manner.

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.

Response retention and apparent motion effect in visual cortex models

Apparent motion is a visual illusion in which stationary stimuli, flashing in distinct spatial locations at certain time intervals, are perceived as one stimulus moving between these locations. In the primary visual cortex, apparent-motion stimuli produce smooth spatio-temporal patterns of activity similar to those produced by continuously moving stimuli. An important prerequisite for producing such activity patterns is prolongation of responses to brief stimuli. Indeed, a brief stimulus can evoke in the visual cortex a long response, outlasting the stimulus by hundreds of milliseconds. Here we use firing-rate based models with simple ring structure, and biologically-detailed conductance-based refractory density (CBRD) model with retinotopic space representation to analyze the response retention and the origin of smooth profiles of activity in response to apparent-motion stimuli. We show that the strength of recurrent connectivity is the major factor that endorses neuronal networks with the ability for response retention. The same strengths of recurrent connections mediate the appearance of bump attractor in the ring models. Factors such as synaptic depression, NMDA receptor mediated currents, and conductances regulating spike adaptation influence response retention, but cannot substitute for the weakness of recurrent connections to reproduce response retention in models with weak connectivity. However, the weakness of lateral recurrent connections can be compensated by layering: in multi-layer models even with weaker connections the activity retains due to its feedforward propagation from layer to layer. Using CBRD model with retinotopic space representation we further show that smooth spatio-temporal profiles of activity in response to apparent-motion stimuli are produced in the models expressing response retention, but not in the models that fail to produce response retention. Together, these results demonstrate a link between response retention and the ability of neuronal networks to generate spatio-temporal patterns of activity, which are compatible with perception of apparent motion.

The inhibitory control of traveling waves in cortical networks

Propagating waves of activity can be evoked and can occur spontaneously in vivo and in vitro in cerebral cortex. These waves are thought to be instrumental in the propagation of information across cortical regions and as a means to modulate the sensitivity of neurons to subsequent stimuli. In normal tissue, the waves are sparse and tightly controlled by inhibition and other negative feedback processes. However, alterations of this balance between excitation and inhibition can lead to pathological behavior such as seizure-type dynamics (with low inhibition) or failure to propagate (with high inhibition). We develop a spiking one-dimensional network of neurons to explore the reliability and control of evoked waves and compare this to a cortical slice preparation where the excitability can be pharmacologically manipulated. We show that the waves enhance sensitivity of the cortical network to stimuli in specific spatial and temporal ways. To gain further insight into the mechanisms of propagation and transitions to pathological behavior, we derive a mean-field model for the synaptic activity. We analyze the mean-field model and a piece-wise constant approximation of it and study the stability of the propagating waves as spatial and temporal properties of the inhibition are altered. We show that that the transition to seizure-like activity is gradual but that the loss of propagation is abrupt and can occur via either the loss of existence of the wave or through a loss of stability leading to complex patterns of propagation.

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.

Neural assemblies coordinated by cortical waves are associated with waking and hallucinatory brain states

The relationship between sensory stimuli and perceptions is brain-state dependent: in wakefulness stimuli evoke perceptions; under anesthesia perceptions are abolished; during dreaming and in dissociated states, percepts are internally generated. Here, we exploit this state dependence to identify brain activity associated with internally generated or stimulus-evoked perception. In awake mice, visual stimuli phase reset spontaneous cortical waves to elicit 3-6 Hz feedback traveling waves. These stimulus-evoked waves traverse the cortex and entrain visual and parietal neurons. Under anesthesia and during ketamine-induced dissociation, visual stimuli do not disrupt spontaneous waves. Uniquely in the dissociated state, spontaneous waves traverse the cortex caudally and entrain visual and parietal neurons, akin to stimulus-evoked waves in wakefulness. Thus, coordinated neuronal assemblies orchestrated by traveling cortical waves emerge in states in which perception can manifest. The awake state is privileged in that this coordination is elicited by specifically by external visual stimuli.

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.

Rhythmic attentional scanning

Sensory processing, short-term memory, and decision-making often deal with multiple items, or options, simultaneously. I review evidence suggesting that the brain handles such multiple items by "rhythmic attentional scanning (RAS)": each item is processed in a separate cycle of the theta rhythm, involving several gamma cycles, to reach an internally consistent representation in the form of a gamma-synchronized neuronal group. Within each theta cycle, items that are extended in representational space are scanned by traveling waves. Such scanning might go across small numbers of simple items linked into a chunk.

Studying memory processes at different levels with simultaneous depth and surface EEG recordings

Investigating cognitive brain functions using non-invasive electrophysiology can be challenging due to the particularities of the task-related EEG activity, the depth of the activated brain areas, and the extent of the networks involved. Stereoelectroencephalographic (SEEG) investigations in patients with drug-resistant epilepsy offer an extraordinary opportunity to validate information derived from non-invasive recordings at macro-scales. The SEEG approach can provide brain activity with high spatial specificity during tasks that target specific cognitive processes (e.g., memory). Full validation is possible only when performing simultaneous scalp SEEG recordings, which allows recording signals in the exact same brain state. This is the approach we have taken in 12 subjects performing a visual memory task that requires the recognition of previously viewed objects. The intracranial signals on 965 contact pairs have been compared to 391 simultaneously recorded scalp signals at a regional and whole-brain level, using multivariate pattern analysis. The results show that the task conditions are best captured by intracranial sensors, despite the limited spatial coverage of SEEG electrodes, compared to the whole-brain non-invasive recordings. Applying beamformer source reconstruction or independent component analysis does not result in an improvement of the multivariate task decoding performance using surface sensor data. By analyzing a joint scalp and SEEG dataset, we investigated whether the two types of signals carry complementary information that might improve the machine-learning classifier performance. This joint analysis revealed that the results are driven by the modality exhibiting best individual performance, namely SEEG.

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.

Simulations approaching data: cortical slow waves in inferred models of the whole hemisphere of mouse

The development of novel techniques to record wide-field brain activity enables estimation of data-driven models from thousands of recording channels and hence across large regions of cortex. These in turn improve our understanding of the modulation of brain states and the richness of traveling waves dynamics. Here, we infer data-driven models from high-resolution in-vivo recordings of mouse brain obtained from wide-field calcium imaging. We then assimilate experimental and simulated data through the characterization of the spatio-temporal features of cortical waves in experimental recordings. Inference is built in two steps: an inner loop that optimizes a mean-field model by likelihood maximization, and an outer loop that optimizes a periodic neuro-modulation via direct comparison of observables that characterize cortical slow waves. The model reproduces most of the features of the non-stationary and non-linear dynamics present in the high-resolution in-vivo recordings of the mouse brain. The proposed approach offers new methods of characterizing and understanding cortical waves for experimental and computational neuroscientists.

Chaotic dynamics in spatially distributed neuronal networks generate population-wide shared variability

Neural activity in the cortex is highly variable in response to repeated stimuli. Population recordings across the cortex demonstrate that the variability of neuronal responses is shared among large groups of neurons and concentrates in a low dimensional space. However, the source of the population-wide shared variability is unknown. In this work, we analyzed the dynamical regimes of spatially distributed networks of excitatory and inhibitory neurons. We found chaotic spatiotemporal dynamics in networks with similar excitatory and inhibitory projection widths, an anatomical feature of the cortex. The chaotic solutions contain broadband frequency power in rate variability and have distance-dependent and low-dimensional correlations, in agreement with experimental findings. In addition, rate chaos can be induced by globally correlated noisy inputs. These results suggest that spatiotemporal chaos in cortical networks can explain the shared variability observed in neuronal population responses.

Precise Spiking Motifs in Neurobiological and Neuromorphic Data

Why do neurons communicate through spikes? By definition, spikes are all-or-none neural events which occur at continuous times. In other words, spikes are on one side binary, existing or not without further details, and on the other, can occur at any asynchronous time, without the need for a centralized clock. This stands in stark contrast to the analog representation of values and the discretized timing classically used in digital processing and at the base of modern-day neural networks. As neural systems almost systematically use this so-called event-based representation in the living world, a better understanding of this phenomenon remains a fundamental challenge in neurobiology in order to better interpret the profusion of recorded data. With the growing need for intelligent embedded systems, it also emerges as a new computing paradigm to enable the efficient operation of a new class of sensors and event-based computers, called neuromorphic, which could enable significant gains in computation time and energy consumption-a major societal issue in the era of the digital economy and global warming. In this review paper, we provide evidence from biology, theory and engineering that the precise timing of spikes plays a crucial role in our understanding of the efficiency of neural networks.

Deterministic and Stochastic Components of Cortical Down States: Dynamics and Modulation

Slow oscillations are an emergent activity of the cerebral cortex network consisting of alternating periods of activity (Up states) and silence (Down states). Up states are periods of persistent cortical activity that share properties with that of underlying wakefulness. However, the occurrence of Down states is almost invariably associated with unconsciousness, both in animal models and clinical studies. Down states have been attributed relevant functions, such as being a resetting mechanism or breaking causal interactions between cortical areas. But what do Down states consist of? Here, we explored in detail the network dynamics (e.g., synchronization and phase) during these silent periods in vivo (male mice), in vitro (ferrets, either sex), and in silico, investigating various experimental conditions that modulate them: anesthesia levels, excitability (electric fields), and excitation/inhibition balance. We identified metastability as two complementary phases composing such quiescence states: a highly synchronized "deterministic" period followed by a low-synchronization "stochastic" period. The balance between these two phases determines the dynamical properties of the resulting rhythm, as well as the responsiveness to incoming inputs or refractoriness. We propose detailed Up and Down state cycle dynamics that bridge cortical properties emerging at the mesoscale with their underlying mechanisms at the microscale, providing a key to understanding unconscious states.SIGNIFICANCE STATEMENT The cerebral cortex expresses slow oscillations consisting of Up (active) and Down (silent) states. Such activity emerges not only in slow wave sleep, but also under anesthesia and in brain lesions. Down states functionally disconnect the network, and are associated with unconsciousness. Based on a large collection of data, novel data analysis approaches and computational modeling, we thoroughly investigate the nature of Down states. We identify two phases: a highly synchronized "deterministic" period, followed by a low-synchronization "stochastic" period. The balance between these two phases determines the dynamic properties of the resulting rhythm and responsiveness to incoming inputs. This finding reconciles different theories of slow rhythm generation and provides clues about how the brain switches from conscious to unconscious brain states.

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.

Large-scale, closed-loop interrogation of neural circuits underlying cognition

Cognitive functions are increasingly understood to involve coordinated activity patterns between multiple brain regions, and their disruption by neuropsychiatric disorders is similarly complex. Closed-loop neurostimulation can directly modulate neural signals with temporal and spatial precision. How to leverage such an approach to effectively identify and target distributed neural networks implicated in mediating cognition remains unclear. We review current conceptual and technical advances in this area, proposing that devices that enable large-scale acquisition, integrated processing, and multiregion, arbitrary waveform stimulation will be critical for mechanistically driven manipulation of cognitive processes in physiological and pathological brain networks.

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.

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.

Waveform detection by deep learning reveals multi-area spindles that are selectively modulated by memory load

Sleep is generally considered to be a state of large-scale synchrony across thalamus and neocortex; however, recent work has challenged this idea by reporting isolated sleep rhythms such as slow oscillations and spindles. What is the spatial scale of sleep rhythms? To answer this question, we adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves in high-noise settings for analysis of neural recordings in sleep. We then studied sleep spindles in non-human primate electrocorticography (ECoG), human electroencephalogram (EEG), and clinical intracranial electroencephalogram (iEEG) recordings in the human. Within each recording type, we find widespread spindles occur much more frequently than previously reported. We then analyzed the spatiotemporal patterns of these large-scale, multi-area spindles and, in the EEG recordings, how spindle patterns change following a visual memory task. Our results reveal a potential role for widespread, multi-area spindles in consolidation of memories in networks widely distributed across primate cortex.

The brain processes memories as we sleep, generating rhythms of electrical activity called ‘sleep spindles’. Sleep spindles were long thought to be a state where the entire brain was fully synchronized by this rhythm. This was based on EEG recordings, short for electroencephalogram, a technique that uses electrodes on the scalp to measure electrical activity in the outermost layer of the brain, the cortex. But more recent intracranial recordings of people undergoing brain surgery have challenged this idea and suggested that sleep spindles may not be a state of global brain synchronization, but rather localised to specific areas.

Mofrad et al. sought to clarify the extent to which spindles co-occur at multiple sites in the brain, which could shed light on how networks of neurons coordinate memory storage during sleep. To analyse highly variable brain wave recordings, Mofrad et al. adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves. The resulting algorithm, designed to more sensitively detect spindles amongst other brain activity, was then applied to a range of sleep recordings from humans and macaque monkeys.

The analyses revealed that widespread and complex patterns of spindle rhythms, spanning multiple areas in the cortex of the brain, actually appear much more frequently than previously thought. This finding was consistent across all the recordings analysed, even recordings under the skull, which provide the clearest window into brain circuits. Further analyses found that these multi-area spindles occurred more often in sleep after people had completed tasks that required holding many visual scenes in memory, as opposed to control conditions with fewer visual scenes.

In summary, Mofrad et al. show that neuroscientists had previously not appreciated the complex and dynamic patterns in this sleep rhythm. These patterns in sleep spindles may be able to adapt based on the demands needed for memory storage, and this will be the subject of future work. Moreover, the findings support the idea that sleep spindles help coordinate the consolidation of memories in brain circuits that stretch across the cortex. Understanding this mechanism may provide insights into how memory falters in aging and sleep-related diseases, such as Alzheimer’s disease. Lastly, the algorithm developed by Mofrad et al. stands to be a useful tool for analysing other rhythmic waveforms in noisy recordings.

Spontaneous Spiking Is Governed by Broadband Fluctuations

Populations of cortical neurons generate rhythmic fluctuations in their ongoing spontaneous activity. These fluctuations can be seen in the local field potential (LFP), which reflects summed return currents from synaptic activity in the local population near a recording electrode. The LFP is spectrally broad, and many researchers view this breadth as containing many narrowband oscillatory components that may have distinct functional roles. This view is supported by the observation that the phase of narrowband oscillations is often correlated with cortical excitability and can relate to the timing of spiking activity and the fidelity of sensory evoked responses. Accordingly, researchers commonly tune in to these channels by narrowband filtering the LFP. Alternatively, neural activity may be fundamentally broadband and composed of transient, nonstationary rhythms that are difficult to approximate as oscillations. In this view, the instantaneous state of the broad ensemble relates directly to the excitability of the local population with no particular allegiance to any frequency band. To test between these alternatives, we asked whether the spiking activity of neocortical neurons in marmoset of either sex is better aligned with the phase of the LFP within narrow frequency bands or with a broadband measure. We find that the phase of broadband LFP fluctuations provides a better predictor of spike timing than the phase after filtering in narrow bands. These results challenge the view of the neocortex as a system composed of narrowband oscillators and supports a view in which neural activity fluctuations are intrinsically broadband.SIGNIFICANCE STATEMENT Research into the dynamical state of neural populations often attributes unique significance to the state of narrowband oscillatory components. However, rhythmic fluctuations in cortical activity are nonstationary and broad spectrum. We find that the timing of spontaneous spiking activity is better captured by the state of broadband fluctuations over any latent oscillatory component. These results suggest narrowband interpretations of rhythmic population activity may be limited, and broader representations may provide higher fidelity in describing moment-to-moment fluctuations in cortical activity.

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.

Spontaneous neuronal oscillations in the human insula are hierarchically organized traveling waves

The insula plays a fundamental role in a wide range of adaptive human behaviors, but its electrophysiological dynamics are poorly understood. Here, we used human intracranial electroencephalographic recordings to investigate the electrophysiological properties and hierarchical organization of spontaneous neuronal oscillations within the insula. We analyzed the neuronal oscillations of the insula directly and found that rhythms in the theta and beta frequency oscillations are widespread and spontaneously present. These oscillations are largely organized along the anterior–posterior (AP) axis of the insula. Both the left and right insula showed anterior-­to-posterior decreasing gradients for the power of oscillations in the beta frequency band. The left insula also showed a posterior-to-anterior decreasing frequency gradient and an anterior-to-posterior decreasing power gradient in the theta frequency band. In addition to measuring the power of these oscillations, we also examined the phase of these signals across simultaneous recording channels and found that the insula oscillations in the theta and beta bands are traveling waves. The strength of the traveling waves in each frequency was positively correlated with the amplitude of each oscillation. However, the theta and beta traveling waves were uncoupled to each other in terms of phase and amplitude, which suggested that insular traveling waves in the theta and beta bands operate independently. Our findings provide new insights into the spatiotemporal dynamics and hierarchical organization of neuronal oscillations within the insula, which, given its rich connectivity with widespread cortical regions, indicates that oscillations and traveling waves have an important role in intrainsular and interinsular communications.

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.

Spatially distributed computation in cortical circuits

The traditional view of neural computation in the cerebral cortex holds that sensory neurons are specialized, i.e., selective for certain dimensions of sensory stimuli. This view was challenged by evidence of contextual interactions between stimulus dimensions in which a neuron's response to one dimension strongly depends on other dimensions. Here, we use methods of mathematical modeling, psychophysics, and electrophysiology to address shortcomings of the traditional view. Using a model of a generic cortical circuit, we begin with the simple demonstration that cortical responses are always distributed among neurons, forming characteristic waveforms, which we call neural waves. When stimulated by patterned stimuli, circuit responses arise by interference of neural waves. Results of this process depend on interaction between stimulus dimensions. Comparison of modeled responses with responses of biological vision makes it clear that the framework of neural wave interference provides a useful alternative to the standard concept of neural computation.

Spontaneous variability in gamma dynamics described by a damped harmonic oscillator driven by noise

Circuits of excitatory and inhibitory neurons generate gamma-rhythmic activity (30-80 Hz). Gamma-cycles show spontaneous variability in amplitude and duration. To investigate the mechanisms underlying this variability, we recorded local-field-potentials (LFPs) and spikes from awake macaque V1. We developed a noise-robust method to detect gamma-cycle amplitudes and durations, which showed a weak but positive correlation. This correlation, and the joint amplitude-duration distribution, is well reproduced by a noise-driven damped harmonic oscillator. This model accurately fits LFP power-spectra, is equivalent to a linear, noise-driven E-I circuit, and recapitulates two additional features of gamma: (1) Amplitude-duration correlations decrease with oscillation strength; (2) amplitudes and durations exhibit strong and weak autocorrelations, respectively, depending on oscillation strength. Finally, longer gamma-cycles are associated with stronger spike-synchrony, but lower spike-rates in both (putative) excitatory and inhibitory neurons. In sum, V1 gamma-dynamics are well described by the simplest possible model of gamma: A damped harmonic oscillator driven by noise.

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.

Revisiting horizontal connectivity rules in V1: from like-to-like towards like-to-all

Horizontal connections in the primary visual cortex of carnivores, ungulates and primates organize on a near-regular lattice. Given the similar length scale for the regularity found in cortical orientation maps, the currently accepted theoretical standpoint is that these maps are underpinned by a like-to-like connectivity rule: horizontal axons connect preferentially to neurons with similar preferred orientation. However, there is reason to doubt the rule's explanatory power, since a growing number of quantitative studies show that the like-to-like connectivity preference and bias mostly observed at short-range scale, are highly variable on a neuron-to-neuron level and depend on the origin of the presynaptic neuron. Despite the wide availability of published data, the accepted model of visual processing has never been revised. Here, we review three lines of independent evidence supporting a much-needed revision of the like-to-like connectivity rule, ranging from anatomy to population functional measures, computational models and to theoretical approaches. We advocate an alternative, distance-dependent connectivity rule that is consistent with new structural and functional evidence: from like-to-like bias at short horizontal distance to like-to-all at long horizontal distance. This generic rule accounts for the observed high heterogeneity in interactions between the orientation and retinotopic domains, that we argue is necessary to process non-trivial stimuli in a task-dependent manner.

Traveling waves in the prefrontal cortex during working memory

Neural oscillations are evident across cortex but their spatial structure is not well- explored. Are oscillations stationary or do they form "traveling waves", i.e., spatially organized patterns whose peaks and troughs move sequentially across cortex? Here, we show that oscillations in the prefrontal cortex (PFC) organized as traveling waves in the theta (4-8Hz), alpha (8-12Hz) and beta (12-30Hz) bands. Some traveling waves were planar but most rotated. The waves were modulated during performance of a working memory task. During baseline conditions, waves flowed bidirectionally along a specific axis of orientation. Waves in different frequency bands could travel in different directions. During task performance, there was an increase in waves in one direction over the other, especially in the beta band.

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.

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.

Modulation of the dynamical state in cortical network models

Cortical neural responses can be modulated by various factors, such as stimulus inputs and the behavior state of the animal. Understanding the circuit mechanisms underlying modulations of network dynamics is important to understand the flexibility of circuit computations. Identifying the dynamical state of a network is an important first step to predict network responses to external stimulus and top-down modulatory inputs. Models in stable or unstable dynamical regimes require different analytic tools to estimate the network responses to inputs and the structure of neural variability. In this article, I review recent cortical models of state-dependent responses and their predictions about the underlying modulatory mechanisms.

Brain Activity Fluctuations Propagate as Waves Traversing the Cortical Hierarchy

The brain exhibits highly organized patterns of spontaneous activity as measured by resting-state functional magnetic resonance imaging (fMRI) fluctuations that are being widely used to assess the brain's functional connectivity. Some evidence suggests that spatiotemporally coherent waves are a core feature of spontaneous activity that shapes functional connectivity, although this has been difficult to establish using fMRI given the temporal constraints of the hemodynamic signal. Here, we investigated the structure of spontaneous waves in human fMRI and monkey electrocorticography. In both species, we found clear, repeatable, and directionally constrained activity waves coursed along a spatial axis approximately representing cortical hierarchical organization. These cortical propagations were closely associated with activity changes in distinct subcortical structures, particularly those related to arousal regulation, and modulated across different states of vigilance. The findings demonstrate a neural origin of spatiotemporal fMRI wave propagation at rest and link it to the principal gradient of resting-state fMRI connectivity.

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.

Control of Brain State Transitions with a Photoswitchable Muscarinic Agonist

The ability to control neural activity is essential for research not only in basic neuroscience, as spatiotemporal control of activity is a fundamental experimental tool, but also in clinical neurology for therapeutic brain interventions. Transcranial-magnetic, ultrasound, and alternating/direct current (AC/DC) stimulation are some available means of spatiotemporal controlled neuromodulation. There is also light-mediated control, such as optogenetics, which has revolutionized neuroscience research, yet its clinical translation is hampered by the need for gene manipulation. As a drug-based light-mediated control, the effect of a photoswitchable muscarinic agonist (Phthalimide-Azo-Iper (PAI)) on a brain network is evaluated in this study. First, the conditions to manipulate M2 muscarinic receptors with light in the experimental setup are determined. Next, physiological synchronous emergent cortical activity consisting of slow oscillations-as in slow wave sleep-is transformed into a higher frequency pattern in the cerebral cortex, both in vitro and in vivo, as a consequence of PAI activation with light. These results open the way to study cholinergic neuromodulation and to control spatiotemporal patterns of activity in different brain states, their transitions, and their links to cognition and behavior. The approach can be applied to different organisms and does not require genetic manipulation, which would make it translational to humans.

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.