Phase reorganization leads to transient β-LFP spatial wave patterns in motor cortex during steady-state movement preparation

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

Intrinsic cortical activity forms traveling waves that modulate sensory-evoked responses and perceptual sensitivity. These intrinsic traveling waves (iTWs) may arise from the coordination of synaptic activity through long-range feature-dependent horizontal connectivity within cortical areas. In a spiking network model that incorporates feature-selective patchy connections, we observe iTW motifs that result from shifts in excitatory/inhibitory balance as action potentials traverse these patchy connections. To test whether feature-selective motifs occur in vivo, we examined data recorded in the middle temporal visual area (Area MT) of marmosets performing a visual detection task. We find that some iTWs form motifs that are feature selective, exhibiting direction-selective modulations in spiking activity. Further, motifs modulate the gain of target-evoked responses and perceptual sensitivity if the target matches the preference of the motif. These results suggest that iTWs are shaped by the patchy horizontal fiber projections in the cortex and can regulate neural and perceptual sensitivity in a feature-selective manner.

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.

Evaluation of abstract rule-based associations in the human premotor cortex during passive observation

Decision-making often manifests in behavior, typically yielding overt motor actions. This complex process requires the registration of sensory information with one's internal representation of the current context, before a categorical judgment of the most appropriate motor behavior can be issued. The construct concept of embodied decision-making encapsulates this sequence of complex processes, whereby behaviorally salient information from the environment is represented in an abstracted space of potential motor actions rather than only in an abstract cognitive "decision" space. Theoretical foundations and some empirical evidence account for support the involvement of premotor cortical circuits in embodied cognitive functions. Animal models show that premotor circuits participate in the registration and evaluation of actions performed by peers in social situations, that is, prior to controlling one's voluntary movements guided by arbitrary stimulus-response rules. However, such evidence from human data is currently limited. Here we used time-resolved magnetoencephalography imaging to characterize activations of the premotor cortex as human participants observed arbitrary, non-biological visual stimuli that either respected or violated a simple stimulus-response association rule. The participants had learned this rule previously, either actively, by performing a motor task (active learning), or passively, by observing a computer perform the same task (passive learning). We discovered that the human premotor cortex is activated during the passive observation of the correct execution of a sequence of events according to a rule learned previously. Premotor activation also differs when the subjects observe incorrect stimulus sequences. These premotor effects are present even when the observed events are of a non-motor, abstract nature, and even when the stimulus-response association rule was learned via passive observations of a computer agent performing the task, without requiring overt motor actions from the human participant. We found evidence of these phenomena by tracking cortical beta-band signaling in temporal alignment with the observation of task events and behavior. We conclude that premotor cortical circuits that are typically engaged during voluntary motor behavior are also involved in the interpretation of events of a non-ecological, unfamiliar nature but related to a learned abstract rule. As such, the present study provides the first evidence of neurophysiological processes of embodied decision-making in human premotor circuits when the observed events do not involve motor actions of a third party.

Spatiotemporal organisation of human sensorimotor beta burst activity

Beta oscillations in human sensorimotor cortex are hallmark signatures of healthy and pathological movement. In single trials, beta oscillations include bursts of intermittent, transient periods of high-power activity. These burst events have been linked to a range of sensory and motor processes, but their precise spatial, spectral, and temporal structure remains unclear. Specifically, a role for beta burst activity in information coding and communication suggests spatiotemporal patterns, or travelling wave activity, along specific anatomical gradients. We here show in human magnetoencephalography recordings that burst activity in sensorimotor cortex occurs in planar spatiotemporal wave-like patterns that dominate along two axes either parallel or perpendicular to the central sulcus. Moreover, we find that the two propagation directions are characterised by distinct anatomical and physiological features. Finally, our results suggest that sensorimotor beta bursts occurring before and after a movement can be distinguished by their anatomical, spectral, and spatiotemporal characteristics, indicating distinct functional roles.

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.

Beta oscillations and waves in motor cortex can be accounted for by the interplay of spatially structured connectivity and fluctuating inputs

The beta rhythm (13-30 Hz) is a prominent brain rhythm. Recordings in primates during instructed-delay reaching tasks have shown that different types of traveling waves of oscillatory activity are associated with episodes of beta oscillations in motor cortex during movement preparation. We propose here a simple model of motor cortex based on local excitatory-inhibitory neuronal populations coupled by long-range excitation, where additionally inputs to the motor cortex from other neural structures are represented by stochastic inputs on the different model populations. We show that the model accurately reproduces the statistics of recording data when these external inputs are correlated on a short time scale (25 ms) and have two different components, one that targets the motor cortex locally and another one that targets it in a global and synchronized way. The model reproduces the distribution of beta burst durations, the proportion of the different observed wave types, and wave speeds, which we show not to be linked to axonal propagation speed. When the long-range connectivity or the local input targets are anisotropic, traveling waves are found to preferentially propagate along the axis where connectivity decays the fastest. Different from previously proposed mechanistic explanations, the model suggests that traveling waves in motor cortex are the reflection of the dephasing by external inputs, putatively of thalamic origin, of an oscillatory activity that would otherwise be spatially synchronized by recurrent connectivity.

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.

Grasp-squeeze adaptation to changes in object compliance leads to dynamic beta-band communication between primary somatosensory and motor cortices

In asking the question of how the brain adapts to changes in the softness of manipulated objects, we studied dynamic communication between the primary sensory and motor cortical areas when nonhuman primates grasp and squeeze an elastically deformable manipulandum to attain an instructed force level. We focused on local field potentials recorded from S1 and M1 via intracortical microelectrode arrays. We computed nonparametric spectral Granger Causality to assess directed cortico-cortical interactions between these two areas. We demonstrate that the time-causal relationship between M1 and S1 is bidirectional in the beta-band (15-30 Hz) and that this interareal communication develops dynamically as the subjects adjust the force of hand squeeze to reach the target level. In particular, the directed interaction is strongest when subjects are focused on maintaining the instructed force of hand squeeze in a steady state for several seconds. When the manipulandum's compliance is abruptly changed, beta-band interareal communication is interrupted for a short period (~ 1 s) and then is re-established once the subject has reached a new steady state. These results suggest that transient beta oscillations can provide a communication subspace for dynamic cortico-cortical S1-M1 interactions during maintenance of steady sensorimotor states.

Synchronization, Stochasticity, and Phase Waves in Neuronal Networks With Spatially-Structured Connectivity

Oscillations in the beta/low gamma range (10–45 Hz) are recorded in diverse neural structures. They have successfully been modeled as sparsely synchronized oscillations arising from reciprocal interactions between randomly connected excitatory (E) pyramidal cells and local interneurons (I). The synchronization of spatially distant oscillatory spiking E–I modules has been well-studied in the rate model framework but less so for modules of spiking neurons. Here, we first show that previously proposed modifications of rate models provide a quantitative description of spiking E–I modules of Exponential Integrate-and-Fire (EIF) neurons. This allows us to analyze the dynamical regimes of sparsely synchronized oscillatory E–I modules connected by long-range excitatory interactions, for two modules, as well as for a chain of such modules. For modules with a large number of neurons (> 10⁵), we obtain results similar to previously obtained ones based on the classic deterministic Wilson-Cowan rate model, with the added bonus that the results quantitatively describe simulations of spiking EIF neurons. However, for modules with a moderate (~ 10⁴) number of neurons, stochastic variations in the spike emission of neurons are important and need to be taken into account. On the one hand, they modify the oscillations in a way that tends to promote synchronization between different modules. On the other hand, independent fluctuations on different modules tend to disrupt synchronization. The correlations between distant oscillatory modules can be described by stochastic equations for the oscillator phases that have been intensely studied in other contexts. On shorter distances, we develop a description that also takes into account amplitude modes and that quantitatively accounts for our simulation data. Stochastic dephasing of neighboring modules produces transient phase gradients and the transient appearance of phase waves. We propose that these stochastically-induced phase waves provide an explanative framework for the observations of traveling waves in the cortex during beta oscillations.

Interactions Between Motor Thalamic Field Potentials and Single-Unit Spiking Are Correlated With Behavior in Rats

Field potential (FP) oscillations are believed to coordinate brain activity over large spatiotemporal scales, with specific features (e.g., phase and power) in discrete frequency bands correlated with motor output. Furthermore, complex correlations between oscillations in distinct frequency bands (phase-amplitude, amplitude-amplitude, and phase-phase coupling) are commonly observed. However, the mechanisms underlying FP-behavior correlations and cross-frequency coupling remain unknown. The thalamus plays a central role in generating many circuit-level neural oscillations, and single-unit activity in motor thalamus (Mthal) is correlated with behavioral output. We, therefore, hypothesized that motor thalamic spiking coordinates motor system FPs and underlies FP-behavior correlations. To investigate this possibility, we recorded wideband motor thalamic (Mthal) electrophysiology as healthy rats performed a two-alternative forced-choice task. Delta (1–4 Hz), beta (13–30 Hz), low gamma (30–70 Hz), and high gamma (70–200 Hz) power were strongly modulated by task performance. As in the cortex, the delta phase was correlated with beta/low gamma power and reaction time. Most interestingly, subpopulations of Mthal neurons defined by their relationship to the behavior exhibited distinct relationships with FP features. Specifically, neurons whose activity was correlated with action selection and movement speed were entrained to delta oscillations. Furthermore, changes in their activity anticipated power fluctuations in beta/low gamma bands. These complex relationships suggest mechanisms for commonly observed FP-FP and spike-FP correlations, as well as subcortical influences on motor output.

Is there an Intrinsic Relationship between LFP Beta Oscillation Amplitude and Firing Rate of Individual Neurons in Macaque Motor Cortex?

The properties of motor cortical local field potential (LFP) beta oscillations have been extensively studied. Their relationship to the local neuronal spiking activity was also addressed. Yet, whether there is an intrinsic relationship between the amplitude of beta oscillations and the firing rate of individual neurons remains controversial. Some studies suggest a mapping of spike rate onto beta amplitude, while others find no systematic relationship. To help resolve this controversy, we quantified in macaque motor cortex the correlation between beta amplitude and neuronal spike count during visuomotor behavior. First, in an analysis termed “task-related correlation”, single-trial data obtained across all trial epochs were included. These correlations were significant in up to 32% of cases and often strong. However, a trial-shuffling control analysis recombining beta amplitudes and spike counts from different trials revealed these task-related correlations to reflect systematic, yet independent, modulations of the 2 signals with the task. Second, in an analysis termed “trial-by-trial correlation”, only data from fixed trial epochs were included, and correlations were calculated across trials. Trial-by-trial correlations were weak and rarely significant. We conclude that there is no intrinsic relationship between the firing rate of individual neurons and LFP beta oscillation amplitude in macaque motor cortex.

Propagating Motor Cortical Dynamics Facilitate Movement Initiation

Voluntary movement initiation involves the modulations of large groups of neurons in the primary motor cortex (M1). Yet similar modulations occur during movement planning when no movement occurs. Here, we show that a sequential spatiotemporal pattern of excitability propagates across M1 prior to the movement initiation in one of two oppositely oriented directions along the rostro-caudal axis. Using spatiotemporal patterns of intracortical microstimulation, we find that reaction time increases significantly when stimulation is delivered against, but not with, the natural propagation direction. Functional connections among M1 units emerge at movement that are oriented along the same rostro-caudal axis but not during movement planning. Finally, we show that beta amplitude profiles can more accurately decode muscle activity when they conform to the natural propagating patterns. These findings provide the first causal evidence that large-scale, propagating patterns of cortical excitability are behaviorally relevant and may be a necessary component of movement initiation.