Search

Characterizing the brain's dynamic pattern of response to an input in electroencephalography (EEG) is not a trivial task due to the entanglement of the complex spontaneous brain activity. In this context, the brain's response can be defined as (1) the additional neural activity components generated after the input or (2) the changes in the ongoing spontaneous activities induced by the input. Moreover, the response can be manifested in multiple features. Three commonly studied examples of features are (1) transient temporal waveform, (2) time-frequency representation, and (3) phase dynamics. The most extensively used method of average event-related potentials (ERPs) captures the first one, while the latter two and other more complex features are attracting increasing attention. However, there has not been much work providing a systematic illustration and guidance for how to effectively exploit multifaceted features in neural cognitive research. Based on a visual oddball ERPs dataset with 200 participants, this work demonstrates how the information from the above-mentioned features are complementary to each other and how they can be integrated based on stereotypical neural-network-based machine learning approaches to better exploit neural dynamic information in basic and applied cognitive research.

Circadian rhythms are internal processes repeating approximately every 24 hours in living organisms. The dominant circadian pacemaker is synchronized to the environmental light-dark cycle. Other circadian pacemakers, which can have noncanonical circadian mechanisms, are revealed by arousing stimuli, such as scheduled feeding, palatable meals and running wheel access, or methamphetamine administration. Organisms also have ultradian rhythms, which have periods shorter than circadian rhythms. However, the biological mechanism, origin, and functional significance of ultradian rhythms are not well-elucidated. The dominant circadian rhythm often masks ultradian rhythms; therefore, we disabled the canonical circadian clock of mice by knocking out Per1/2/3 genes, where Per1 and Per2 are essential components of the mammalian light-sensitive circadian mechanism. Furthermore, we recorded wheel-running activity every minute under constant darkness for 272 days. We then investigated rhythmic components in the absence of external influences, applying unique multiscale time-resolved methods to analyze the oscillatory dynamics with time-varying frequencies. We found four rhythmic components with periods of ∼17 h, ∼8 h, ∼4 h, and ∼20 min. When the ∼17-h rhythm was prominent, the ∼8-h rhythm was of low amplitude. This phenomenon occurred periodically approximately every 2-3 weeks. We found that the ∼4-h and ∼20-min rhythms were harmonics of the ∼8-h rhythm. Coupling analysis of the ridge-extracted instantaneous frequencies revealed strong and stable phase coupling from the slower oscillations (∼17, ∼8, and ∼4 h) to the faster oscillations (∼20 min), and weak and less stable phase coupling in the reverse direction and between the slower oscillations. Together, this study elucidated the relationship between the oscillators in the absence of the canonical circadian clock, which is critical for understanding their functional significance. These studies are essential as disruption of circadian rhythms contributes to diseases, such as cancer and obesity, as well as mood disorders.

We develop a technique for the multivariate data analysis of perturbed self-sustained oscillators. The approach is based on the reconstruction of the phase dynamics model from observations and on a subsequent exploration of this model. For the system, driven by several inputs, we suggest a dynamical disentanglement procedure, allowing us to reconstruct the variability of the system's output that is due to a particular observed input, or, alternatively, to reconstruct the variability which is caused by all the inputs except for the observed one. We focus on the application of the method to the vagal component of the heart rate variability caused by a respiratory influence. We develop an algorithm that extracts purely respiratory-related variability, using a respiratory trace and times of R-peaks in the electrocardiogram. The algorithm can be applied to other systems where the observed bivariate data can be represented as a point process and a slow continuous signal, e.g. for the analysis of neuronal spiking. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

A foremost challenge in modern network science is the inverse problem of reconstruction (inference) of coupling equations and network topology from the measurements of the network dynamics. Of particular interest are the methods that can operate on real (empirical) data without interfering with the system. One such earlier attempt (Tokuda et al. 2007 Phys. Rev. Lett. 99 , 064101. (doi:10.1103/PhysRevLett.99.064101)) was a method suited for general limit-cycle oscillators, yielding both oscillators' natural frequencies and coupling functions between them (phase equations) from empirically measured time series. The present paper reviews the above method in a way comprehensive to domain-scientists other than physics. It also presents applications of the method to (i) detection of the network connectivity, (ii) inference of the phase sensitivity function, (iii) approximation of the interaction among phase-coherent chaotic oscillators, and (iv) experimental data from a forced Van der Pol electric circuit. This reaffirms the range of applicability of the method for reconstructing coupling functions and makes it accessible to a much wider scientific community. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Modeling and interpreting (partial) synchronous neural activity can be a challenge. We illustrate this by deriving the phase dynamics of two seminal neural mass models: the Wilson-Cowan firing rate model and the voltage-based Freeman model. We established that the phase dynamics of these models differed qualitatively due to an attractive coupling in the first and a repulsive coupling in the latter. Using empirical structural connectivity matrices, we determined that the two dynamics cover the functional connectivity observed in resting state activity. We further searched for two pivotal dynamical features that have been reported in many experimental studies: (1) a partial phase synchrony with a possibility of a transition towards either a desynchronized or a (fully) synchronized state; (2) long-term autocorrelations indicative of a scale-free temporal dynamics of phase synchronization. Only the Freeman phase model exhibited scale-free behavior. Its repulsive coupling, however, let the individual phases disperse and did not allow for a transition into a synchronized state. The Wilson-Cowan phase model, by contrast, could switch into a (partially) synchronized state, but it did not generate long-term correlations although being located close to the onset of synchronization, i.e. in its critical regime. That is, the phase-reduced models can display one of the two dynamical features, but not both., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

The dispersionless Whitham modulation equations in 2+1 (two space dimensions and time) are reviewed and the instabilities identified. The modulation theory is then reformulated, near the Lighthill instability threshold, with a slow phase, moving frame and different scalings. The resulting nonlinear phase modulation equation near the Lighthill surfaces is a geometric form of the 2+1 two-way Boussinesq equation. This equation is universal in the same sense as Whitham theory. Moreover, it is dispersive, and it has a wide range of interesting multi-periodic, quasi-periodic and multi-pulse localized solutions. For illustration the theory is applied to a complex nonlinear 2+1 Klein-Gordon equation which has two Lighthill surfaces in the manifold of periodic travelling waves.This article is part of the theme issue 'Stability of nonlinear waves and patterns and related topics'., (© 2018 The Author(s).)

A dynamic system showing stable rhythmic activity can be represented by the dynamics of phase oscillators. This would provide a useful mathematical framework through which one can understand the system's dynamic properties. A recent study proposed a Bayesian approach capable of extracting the underlying phase dynamics directly from time-series data of a system showing rhythmic activity. Here we extended this method to spike data that otherwise provide only limited phase information. To determine how this method performs with spike data, we applied it to simulated spike data generated by a realistic neuronal network model. We then compared the estimated dynamics obtained based on the spike data with the dynamics theoretically derived from the model. The method successfully extracted the modeled phase dynamics, particularly the interaction function, when the amount of available data was sufficiently large. Furthermore, the method was able to infer synaptic connections based on the estimated interaction function. Thus, the method was found to be applicable to spike data and practical for understanding the dynamic properties of rhythmic neural systems.

Light is the primary signal that calibrates circadian neural circuits and thus coordinates daily physiological and behavioral rhythms with solar entrainment cues. Drosophila and mammalian circadian circuits consist of diverse populations of cellular oscillators that exhibit a wide range of dynamic light responses, periods, phases, and degrees of synchrony. How heterogeneous circadian circuits can generate robust physiological rhythms while remaining flexible enough to respond to synchronizing stimuli has long remained enigmatic. Cryptochrome is a short-wavelength photoreceptor that is endogenously expressed in approximately half of Drosophila circadian neurons. In a previous study, physiological light response was measured using real-time bioluminescence recordings in Drosophila whole-brain explants, which remain intrinsically light-sensitive. Here we apply analysis of real-time bioluminescence experimental data to show detailed dynamic ensemble representations of whole circadian circuit light entrainment at single neuron resolution. Organotypic whole-brain explants were either maintained in constant darkness (DD) for 6 days or exposed to a phase-advancing light pulse on the second day. We find that stronger circadian oscillators support robust overall circuit rhythmicity in DD, whereas weaker oscillators can be pushed toward transient desynchrony and damped amplitude to facilitate a new state of phase-shifted network synchrony. Additionally, we use mathematical modeling to examine how a network composed of distinct oscillator types can give rise to complex dynamic signatures in DD conditions and in response to simulated light pulses. Simulations suggest that complementary coupling mechanisms and a combination of strong and weak oscillators may enable a robust yet flexible circadian network that promotes both synchrony and entrainment. A more complete understanding of how the properties of oscillators and their signaling mechanisms facilitate their distinct roles in light entrainment may allow us to direct and augment the circadian system to speed recovery from jet lag, shift work, and seasonal affective disorder., (© 2016 The Author(s).)

Cerebral palsy (CP) includes a broad range of disorders, which can result in impairment of posture and movement control. Brain-computer interfaces (BCIs) have been proposed as assistive devices for individuals with CP. Better understanding of the neural processing underlying motor control in affected individuals could lead to more targeted BCI rehabilitation and treatment options. We have explored well-known neural correlates of movement, including event-related desynchronization (ERD), phase synchrony, and a recently-introduced measure of phase dynamics, in participants with CP and healthy control participants. Although present, significantly less ERD and phase locking were found in the group with CP. Additionally, inter-group differences in phase dynamics were also significant. Taken together these findings suggest that users with CP exhibit lower levels of motor cortex activation during motor imagery, as reflected in lower levels of ongoing mu suppression and less functional connectivity. These differences indicate that development of BCIs for individuals with CP may pose additional challenges beyond those faced in providing BCIs to healthy individuals.

This paper studies the effects of distributed-delay coupling on the dynamics in a system of non-identical coupled Stuart-Landau oscillators. For uniform and gamma delay distribution kernels, the conditions for amplitude death are obtained in terms of average frequency, frequency detuning and the parameters of the coupling, including coupling strength and phase, as well as the mean time delay and the width of the delay distribution. To gain further insights into the dynamics inside amplitude death regions, the eigenvalues of the corresponding characteristic equations are computed numerically. Oscillatory dynamics of the system is also investigated, using amplitude and phase representation. Various branches of phase-locked solutions are identified, and their stability is analysed for different types of delay distributions.