Mean-Field Approximations With Adaptive Coupling for Networks With Spike-Timing-Dependent Plasticity.

Understanding the effect of spike-timing-dependent plasticity (STDP) is key to elucidating how neural networks change over long timescales and to design interventions aimed at modulating such networks in neurological disorders. However, progress is restricted by the significant computational cost associated with simulating neural network models with STDP and by the lack of low-dimensional description that could provide analytical insights. Phase-difference-dependent plasticity (PDDP) rules approximate STDP in phase oscillator networks, which prescribe synaptic changes based on phase differences of neuron pairs rather than differences in spike timing. Here we construct mean-field approximations for phase oscillator networks with STDP to describe part of the phase space for this very high-dimensional system. We first show that single-harmonic PDDP rules can approximate a simple form of symmetric STDP, while multiharmonic rules are required to accurately approximate causal STDP. We then derive exact expressions for the evolution of the average PDDP coupling weight in terms of network synchrony. For adaptive networks of Kuramoto oscillators that form clusters, we formulate a family of low-dimensional descriptions based on the mean-field dynamics of each cluster and average coupling weights between and within clusters. Finally, we show that such a two-cluster mean-field model can be fitted to synthetic data to provide a low-dimensional approximation of a full adaptive network with symmetric STDP. Our framework represents a step toward a low-dimensional description of adaptive networks with STDP, and could for example inform the development of new therapies aimed at maximizing the long-lasting effects of brain stimulation.

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.

Predictive coding models for pain perception.

Pain is a complex, multidimensional experience that involves dynamic interactions between sensory-discriminative and affective-emotional processes. Pain experiences have a high degree of variability depending on their context and prior anticipation. Viewing pain perception as a perceptual inference problem, we propose a predictive coding paradigm to characterize evoked and non-evoked pain. We record the local field potentials (LFPs) from the primary somatosensory cortex (S1) and the anterior cingulate cortex (ACC) of freely behaving rats-two regions known to encode the sensory-discriminative and affective-emotional aspects of pain, respectively. We further use predictive coding to investigate the temporal coordination of oscillatory activity between the S1 and ACC. Specifically, we develop a phenomenological predictive coding model to describe the macroscopic dynamics of bottom-up and top-down activity. Supported by recent experimental data, we also develop a biophysical neural mass model to describe the mesoscopic neural dynamics in the S1 and ACC populations, in both naive and chronic pain-treated animals. Our proposed predictive coding models not only replicate important experimental findings, but also provide new prediction about the impact of the model parameters on the physiological or behavioral read-out-thereby yielding mechanistic insight into the uncertainty of expectation, placebo or nocebo effect, and chronic pain.

Next-generation neural mass and field modeling.

The Wilson-Cowan population model of neural activity has greatly influenced our understanding of the mechanisms for the generation of brain rhythms and the emergence of structured brain activity. As well as the many insights that have been obtained from its mathematical analysis, it is now widely used in the computational neuroscience community for building large-scale in silico brain networks that can incorporate the increasing amount of knowledge from the Human Connectome Project. Here, we consider a neural population model in the spirit of that originally developed by Wilson and Cowan, albeit with the added advantage that it can account for the phenomena of event-related synchronization and desynchronization. This derived mean-field model provides a dynamic description for the evolution of synchrony, as measured by the Kuramoto order parameter, in a large population of quadratic integrate-and-fire model neurons. As in the original Wilson-Cowan framework, the population firing rate is at the heart of our new model; however, in a significant departure from the sigmoidal firing rate function approach, the population firing rate is now obtained as a real-valued function of the complex-valued population synchrony measure. To highlight the usefulness of this next-generation Wilson-Cowan style model, we deploy it in a number of neurobiological contexts, providing understanding of the changes in power spectra observed in electro- and magnetoencephalography neuroimaging studies of motor cortex during movement, insights into patterns of functional connectivity observed during rest and their disruption by transcranial magnetic stimulation, and to describe wave propagation across cortex.

A neuromechanistic model for rhythmic beat generation.

When listening to music, humans can easily identify and move to the beat. Numerous experimental studies have identified brain regions that may be involved with beat perception and representation. Several theoretical and algorithmic approaches have been proposed to account for this ability. Related to, but different from the issue of how we perceive a beat, is the question of how we learn to generate and hold a beat. In this paper, we introduce a neuronal framework for a beat generator that is capable of learning isochronous rhythms over a range of frequencies that are relevant to music and speech. Our approach combines ideas from error-correction and entrainment models to investigate the dynamics of how a biophysically-based neuronal network model synchronizes its period and phase to match that of an external stimulus. The model makes novel use of on-going faster gamma rhythms to form a set of discrete clocks that provide estimates, but not exact information, of how well the beat generator spike times match those of a stimulus sequence. The beat generator is endowed with plasticity allowing it to quickly learn and thereby adjust its spike times to achieve synchronization. Our model makes generalizable predictions about the existence of asymmetries in the synchronization process, as well as specific predictions about resynchronization times after changes in stimulus tempo or phase. Analysis of the model demonstrates that accurate rhythmic time keeping can be achieved over a range of frequencies relevant to music, in a manner that is robust to changes in parameters and to the presence of noise.

Relationships Between Neuronal Oscillatory Amplitude and Dynamic Functional Connectivity.

Event-related fluctuations of neural oscillatory amplitude are reported widely in the context of cognitive processing and are typically interpreted as a marker of brain "activity". However, the precise nature of these effects remains unclear; in particular, whether such fluctuations reflect local dynamics, integration between regions, or both, is unknown. Here, using magnetoencephalography, we show that movement induced oscillatory modulation is associated with transient connectivity between sensorimotor regions. Further, in resting-state data, we demonstrate a significant association between oscillatory modulation and dynamic connectivity. A confound with such empirical measurements is that increased amplitude necessarily means increased signal-to-noise ratio (SNR): this means that the question of whether amplitude and connectivity are genuinely coupled, or whether increased connectivity is observed purely due to increased SNR is unanswered. Here, we counter this problem by analogy with computational models which show that, in the presence of global network coupling and local multistability, the link between oscillatory modulation and long-range connectivity is a natural consequence of neural networks. Our results provide evidence for the notion that connectivity is mediated by neural oscillations, and suggest that time-frequency spectrograms are not merely a description of local synchrony but also reflect fluctuations in long-range connectivity.

Order-indeterminant event-based maps for learning a beat.

The process by which humans synchronize to a musical beat is believed to occur through error-correction where an individual's estimates of the period and phase of the beat time are iteratively adjusted to align with an external stimuli. Mathematically, error-correction can be described using a two-dimensional map where convergence to a fixed point corresponds to synchronizing to the beat. In this paper, we show how a neural system, called a beat generator, learns to adapt its oscillatory behavior through error-correction to synchronize to an external periodic signal. We construct a two-dimensional event-based map, which iteratively adjusts an internal parameter of the beat generator to speed up or slow down its oscillatory behavior to bring it into synchrony with the periodic stimulus. The map is novel in that the order of events defining the map are not a priori known. Instead, the type of error-correction adjustment made at each iterate of the map is determined by a sequence of expected events. The map possesses a rich repertoire of dynamics, including periodic solutions and chaotic orbits.

How do spatially distinct frequency specific MEG networks emerge from one underlying structural connectome? The role of the structural eigenmodes.

Functional networks obtained from magnetoencephalography (MEG) from different frequency bands show distinct spatial patterns. It remains to be elucidated how distinct spatial patterns in MEG networks emerge given a single underlying structural network. Recent work has suggested that the eigenmodes of the structural network might serve as a basis set for functional network patterns in the case of functional MRI. Here, we take this notion further in the context of frequency band specific MEG networks. We show that a selected set of eigenmodes of the structural network can predict different frequency band specific networks in the resting state, ranging from delta (1-4 Hz) to the high gamma band (40-70 Hz). These predictions outperform predictions based from surrogate data, suggesting a genuine relationship between eigenmodes of the structural network and frequency specific MEG networks. We then show that the relevant set of eigenmodes can be excited in a network of neural mass models using linear stability analysis only by including delays. Excitation of an eigenmode in this context refers to a dynamic instability of a network steady state to a spatial pattern with a corresponding coherent temporal oscillation. Simulations verify the results from linear stability analysis and suggest that theta, alpha and beta band networks emerge very near to the bifurcation. The delta and gamma bands in the resting state emerges further away from the bifurcation. These results show for the first time how delayed interactions can excite the relevant set of eigenmodes that give rise to frequency specific functional connectivity patterns.

Firing rate models for gamma oscillations.

Gamma oscillations are readily observed in a variety of brain regions during both waking and sleeping states. Computational models of gamma oscillations typically involve simulations of large networks of synaptically coupled spiking units. These networks can exhibit strongly synchronized gamma behavior, whereby neurons fire in near synchrony on every cycle, or weakly modulated gamma behavior, corresponding to stochastic, sparse firing of the individual units on each cycle of the population gamma rhythm. These spiking models offer valuable biophysical descriptions of gamma oscillations; however, because they involve many individual neuronal units they are limited in their ability to communicate general network-level dynamics. Here we demonstrate that few-variable firing rate models with established synaptic timescales can account for both strongly synchronized and weakly modulated gamma oscillations. These models go beyond the classical formulations of rate models by including at least two dynamic variables per population: firing rate and synaptic activation. The models' flexibility to capture the broad range of gamma behavior depends directly on the timescales that represent recruitment of the excitatory and inhibitory firing rates. In particular, we find that weakly modulated gamma oscillations occur robustly when the recruitment timescale of inhibition is faster than that of excitation. We present our findings by using an extended Wilson-Cowan model and a rate model derived from a network of quadratic integrate-and-fire neurons. These biophysical rate models capture the range of weakly modulated and coherent gamma oscillations observed in spiking network models, while additionally allowing for greater tractability and systems analysis. NEW & NOTEWORTHY Here we develop simple and tractable models of gamma oscillations, a dynamic feature observed throughout much of the brain with significant correlates to behavior and cognitive performance in a variety of experimental contexts. Our models depend on only a few dynamic variables per population, but despite this they qualitatively capture features observed in previous biophysical models of gamma oscillations that involve many individual spiking units.

A mean field model for movement induced changes in the beta rhythm.

In electrophysiological recordings of the brain, the transition from high amplitude to low amplitude signals are most likely caused by a change in the synchrony of underlying neuronal population firing patterns. Classic examples of such modulations are the strong stimulus-related oscillatory phenomena known as the movement related beta decrease (MRBD) and post-movement beta rebound (PMBR). A sharp decrease in neural oscillatory power is observed during movement (MRBD) followed by an increase above baseline on movement cessation (PMBR). MRBD and PMBR represent important neuroscientific phenomena which have been shown to have clinical relevance. Here, we present a parsimonious model for the dynamics of synchrony within a synaptically coupled spiking network that is able to replicate a human MEG power spectrogram showing the evolution from MRBD to PMBR. Importantly, the high-dimensional spiking model has an exact mean field description in terms of four ordinary differential equations that allows considerable insight to be obtained into the cause of the experimentally observed time-lag from movement termination to the onset of PMBR (∼ 0.5 s), as well as the subsequent long duration of PMBR (∼ 1 - 10 s). Our model represents the first to predict these commonly observed and robust phenomena and represents a key step in their understanding, in health and disease.

Knowledge and attitudes of allied health professional students regarding the stroke rehabilitation team and the role of the Speech and Language Therapist.

BACKGROUND: One of the major barriers to effective team working among healthcare professionals is a lack of knowledge of each other's roles. The importance of understanding Irish healthcare students' attitudes towards team working and each other's roles led to the development of this study. AIMS: The aims were to investigate allied health professional students' perceptions and experiences of the stroke rehabilitation team and the role of the Speech and Language Therapist (SLT). METHODS & PROCEDURES: A survey first developed by Felsher and Ross (1994) and further developed by Insalaco et al. (2007) was adapted to the Irish healthcare setting. The survey was administered to final-year Occupational Therapy (n = 23), Speech and Language Therapy (21) students and Physiotherapy (20) students (64 in total) (a 98.5% response rate). OUTCOMES & RESULTS: Results indicate that students had a good understanding of teamwork in the healthcare setting and the possible benefits and challenges it presents. Students had a strong appreciation for interprofessional collaboration, with the majority (79%) choosing shared leadership as their preferred option for the stroke rehabilitation team. Further to this, the team approaches that students felt were most appropriate for the stroke rehabilitation setting were the more collaborative approaches of interdisciplinary (43.5%) and transdisciplinary (37.1%). The students had clear perceptions of the SLT's role in aphasia, dysphagia, dysarthria, apraxia and auditory agnosia, but were less knowledgeable of the SLT's role in the acquired disorders of alexia and agraphia (p < 0.05). More than half of all students perceived that the SLT is involved in the treatment of hemispatial neglect (55.5%), depression (71.5%) and visual agnosia (59.4%). CONCLUSIONS & IMPLICATIONS: The results provide valuable information for further developments in interprofessional education at an undergraduate level. Further opportunities should be provided to students to collaborate with each other, particularly in their final year of training as, by then, students have a well-established knowledge of their own roles and would be more capable of sharing this role with other professions. Through this collaboration students would also gain valuable insight into the importance of teamwork, which they could take with them into their professional careers.

Piecewise-constant optimal control strategies for controlling the outbreak of COVID-19 in the Irish population.

We introduce a deterministic SEIR model and fit it to epidemiological data for the COVID-19 outbreak in Ireland. We couple the model to economic considerations - we formulate an optimal control problem in which the cost to the economy of the various non-pharmaceutical interventions is minimized, subject to hospital admissions never exceeding a threshold value corresponding to health-service capacity. Within the framework of the model, the optimal strategy of disease control is revealed to be one of disease suppression, rather than disease mitigation.

Auditory streaming and bistability paradigm extended to a dynamic environment.

We explore stream segregation with temporally modulated acoustic features using behavioral experiments and modelling. The auditory streaming paradigm in which alternating high- A and low-frequency tones B appear in a repeating ABA-pattern, has been shown to be perceptually bistable for extended presentations (order of minutes). For a fixed, repeating stimulus, perception spontaneously changes (switches) at random times, every 2-15 s, between an integrated interpretation with a galloping rhythm and segregated streams. Streaming in a natural auditory environment requires segregation of auditory objects with features that evolve over time. With the relatively idealized ABA-triplet paradigm, we explore perceptual switching in a non-static environment by considering slowly and periodically varying stimulus features. Our previously published model captures the dynamics of auditory bistability and predicts here how perceptual switches are entrained, tightly locked to the rising and falling phase of modulation. In psychoacoustic experiments we find that entrainment depends on both the period of modulation and the intrinsic switch characteristics of individual listeners. The extended auditory streaming paradigm with slowly modulated stimulus features presented here will be of significant interest for future imaging and neurophysiology experiments by reducing the need for subjective perceptual reports of ongoing perception.

Next-generation neural field model: The evolution of synchrony within patterns and waves.

Neural field models are commonly used to describe wave propagation and bump attractors at a tissue level in the brain. Although motivated by biology, these models are phenomenological in nature. They are built on the assumption that the neural tissue operates in a near synchronous regime, and hence, cannot account for changes in the underlying synchrony of patterns. It is customary to use spiking neural network models when examining within population synchronization. Unfortunately, these high-dimensional models are notoriously hard to obtain insight from. In this paper, we consider a network of θ-neurons, which has recently been shown to admit an exact mean-field description in the absence of a spatial component. We show that the inclusion of space and a realistic synapse model leads to a reduced model that has many of the features of a standard neural field model coupled to a further dynamical equation that describes the evolution of network synchrony. Both Turing instability analysis and numerical continuation software are used to explore the existence and stability of spatiotemporal patterns in the system. In particular, we show that this new model can support states above and beyond those seen in a standard neural field model. These states are typified by structures within bumps and waves showing the dynamic evolution of population synchrony.