Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Articles, Systems/Circuits

Shaping Intrinsic Neural Oscillations with Periodic Stimulation

Christoph S. Herrmann, Micah M. Murray, Silvio Ionta, Axel Hutt and Jérémie Lefebvre
Journal of Neuroscience 11 May 2016, 36 (19) 5328-5337; DOI: https://doi.org/10.1523/JNEUROSCI.0236-16.2016
Christoph S. Herrmann
1Experimental Psychology Laboratory, Department of Psychology, Cluster of Excellence Hearing4all and European Medical School and
2Research Center Neurosensory Science, Carl-von-Ossietzky University Oldenburg, 26111 Oldenburg, Germany,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Christoph S. Herrmann
Micah M. Murray
3The Laboratory for Investigative Neurophysiology (The LINE), Department of Clinical Neurosciences and Department of Radiology, University Hospital Center and University of Lausanne, Lausanne 1011, Switzerland,
4Electroencephalography Brain Mapping Core, Center for Biomedical Imaging (CIBM), Lausanne 1011, Switzerland,
5Department of Ophthalmology, University of Lausanne, Jules-Gonin Eye Hospital, Lausanne 1000, Switzerland,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Micah M. Murray
Silvio Ionta
3The Laboratory for Investigative Neurophysiology (The LINE), Department of Clinical Neurosciences and Department of Radiology, University Hospital Center and University of Lausanne, Lausanne 1011, Switzerland,
6Rehabilitation Engineering Laboratory, Department of Health Sciences and Technology, ETH Zürich, 8092 Zürich, Switzerland,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Axel Hutt
7Deutscher Wetterdienst, Section FE12 Data Assimilation, 63067 Offenbach am Main, Germany,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jérémie Lefebvre
3The Laboratory for Investigative Neurophysiology (The LINE), Department of Clinical Neurosciences and Department of Radiology, University Hospital Center and University of Lausanne, Lausanne 1011, Switzerland,
8Krembil Research Institute, University Health Network, Toronto, Ontario M5T 2S8, Canada, and
9Department of Mathematics, University of Toronto, Toronto, Ontario M5S 3G3, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Rhythmic brain activity plays an important role in neural processing and behavior. Features of these oscillations, including amplitude, phase, and spectrum, can be influenced by internal states (e.g., shifts in arousal, attention or cognitive ability) or external stimulation. Electromagnetic stimulation techniques such as transcranial magnetic stimulation, transcranial direct current stimulation, and transcranial alternating current stimulation are used increasingly in both research and clinical settings. Currently, the mechanisms whereby time-dependent external stimuli influence population-scale oscillations remain poorly understood. Here, we provide computational insights regarding the mapping between periodic pulsatile stimulation parameters such as amplitude and frequency and the response dynamics of recurrent, nonlinear spiking neural networks. Using a cortical model built of excitatory and inhibitory neurons, we explored a wide range of stimulation intensities and frequencies systematically. Our results suggest that rhythmic stimulation can form the basis of a control paradigm in which one can manipulate the intrinsic oscillatory properties of driven networks via a plurality of input-driven mechanisms. Our results show that, in addition to resonance and entrainment, nonlinear acceleration is involved in shaping the rhythmic response of our modeled network. Such nonlinear acceleration of spontaneous and synchronous oscillatory activity in a neural network occurs in regimes of intense, high-frequency rhythmic stimulation. These results open new perspectives on the manipulation of synchronous neural activity for basic and clinical research.

SIGNIFICANCE STATEMENT Oscillatory activity is widely recognized as a core mechanism for information transmission within and between brain circuits. Noninvasive stimulation methods can shape this activity, something that is increasingly capitalized upon in basic research and clinical practice. Here, we provide computational insights on the mechanistic bases for such effects. Our results show that rhythmic stimulation forms the basis of a control paradigm in which one can manipulate the intrinsic oscillatory properties of driven networks via a plurality of input-driven mechanisms. In addition to resonance and entrainment, nonlinear acceleration is involved in shaping the rhythmic response of our modeled network, particularly in regimes of high-frequency rhythmic stimulation. These results open new perspectives on the manipulation of synchronous neural activity for basic and clinical research.

  • dynamics
  • networks
  • oscillations
  • spectrum
  • stimulation
  • synchrony

Introduction

Rhythmic neural dynamics are a mechanism through which internal states exercise (top-down) influences on stimulus processing and in turn impact perception (Engel et al., 2001; Varela et al., 2001). Synchronization is considered one such mechanism, operating across multiple scales of brain circuitry from single units to interregional coherence (Singer and Gray, 1995; Engel and Singer, 2001). The functional role of oscillations within the gamma band (30 Hz–100 Hz) has received considerable attention (Whittington et al., 1995; Wang and Buzsáki, 1996; Jadi and Sejnowski, 2014). More recently, lower frequencies, which likely rely on different cellular mechanisms and operate over longer timescales than their high-frequency counterparts (Kajikawa and Schroeder, 2015; Podvalny et al., 2015), have been the focus of increasingly intensive research (Lorinscz et al., 2009; Haegens et al., 2014). In particular, the alpha frequency (8–12 Hz) is known to correlate strongly with cognitive performance and serves as a robust predictor of the degree of cognitive decline between healthy and pathological groups (Klimesch, 1999; Başar, 2012; Basar and Guntekin, 2012; Vlahou et al., 2014). The alpha peak frequency has also been found to predict the temporal resolution of visual perception (Samaha and Postle, 2015). Therefore, numerous studies have investigated the consequences of variations in alpha oscillations on cognition, attention, and perception to increase our understanding of the reciprocal interactions between slow oscillatory activity and population neural coding (Lakatos et al., 2008; Siegel et al., 2012). For example, shifts in the peak alpha frequency of different subjects (also called individual alpha frequency) are observed reliably when attentional demands increase during visual tasks (Samaha et al., 2015), overt behavior (Gutmann et al., 2015), and/or sensorimotor inputs (Hülsdünker et al., 2015).

These observations have raised the prospect of manipulating (alpha) rhythmic activity in a controlled manner either through attention-related processes or brain stimulation methods (Thut et al., 2012; Frohlich, 2015). It has indeed been demonstrated that transcranial magnetic stimulation (TMS) and transcranial direct current stimulation and transcranial alternating current stimulation (tACS) can perturb individual alpha oscillations and have a direct effect on visual stimulus perception (Chanes et al., 2013; Cecere et al., 2015) and task performance by reinforcing endogenous slow-wave rhythms (Klimesch, 1999; Henry et al., 2014).

One implication of such findings is that brain stimulation can result in the selective engagement of endogenous (intrinsic) oscillations, effectively making it a potential means of manipulating and controlling cognition and treating neurobiological disorders. This critical question has received much attention in recent studies (Frohlich, 2015; Witkowski et al., 2015) in which mechanisms such as resonance (Thut et al., 2012; Ali et al., 2013) and entrainment (Thut et al., 2011) have been implicated and found to depend strongly on the ongoing dynamics of the stimulated networks (Neuling et al., 2013; Alagapan et al., 2016).

In order for brain stimulation approaches to be optimally effective, it is necessary to have a better mechanistic understanding of their consequences on brain activity, something that can be garnered from computational modeling. However, computational modeling of this issue has rarely considered a wide range of stimulation parameters (Reato et al., 2010). Due to the nonlinearity of neural circuits and the large number of possible stimulation patterns, the mechanisms involved in the combination of exogenous and endogenous oscillatory signals are furthermore poorly understood. To address these limitations, we here consider the response of a simplified cortical network to periodic pulsatile stimuli. Specifically, using a generic computational model of spiking excitatory and inhibitory neurons that exhibits spontaneous alpha oscillatory activity (here set at 10 Hz), we explore different combinations of stimulation frequencies and amplitudes systematically, revealing those in which resonance and/or rhythmic entrainment can be evoked. Then, building on recent advances in neurodynamics (Lefebvre et al., 2015), we report novel regimes of oscillatory acceleration that characterize network responses to input stimuli of high frequencies (50–100 Hz). Our results provide new perspectives on the response of synchronous nonlinear neural systems, in which a plurality of linear and nonlinear mechanisms are combined to shape intrinsic alpha-like oscillatory activity.

Materials and Methods

Cortical network.

To highlight the influence of stimulation on local oscillatory dynamics, we deliberately chose a relatively well understood hybrid Amari model (Amari, 1977) of a cortical network and evaluated its behavior when subjected to dynamic stimulation, that is, during forcing. Our model takes into account the variability of neuronal responses, yet is simplistic enough to remain both computationally and analytically tractable. It consists of a network of pyramidal excitatory neurons (e) and inhibitory interneurons (i), the individual spiking activity of which follows the nonhomogeneous Poisson processes as follows: Embedded Image where Xe,ij(t) = ∑{tl}δe,ij(t − tl) is the spike train of the jth neuron in the population e or i. For simplicity, the activation function is defined by f[u] = 1 whenever u > h and 0 otherwise. It represents how the firing rate of the neurons relates to the cellular potential, u. The excitatory and inhibitory somatic membrane potentials uej(t) and uij(t) evolve according to the set of scalar nonlinear stochastic equations as follows: Embedded Image Here, Ne and Ni refer to the number of e versus i cells, which follow the 4:1 ratio generally observed in cortical tissue (Koch, 1999). The rates αeu, αiu define the timescale of the respective somatic membrane. Adaptation currents vej(t) and vij(t) with gain b and rate a are also included to represent the habituation of neurons to both steady and fluctuating inputs (Prescott and Sejnowski, 2008).

Both excitatory and inhibitory populations are driven by an external stimulation signal I(t), which consists of pulsatile stimuli. In this model, recurrent inputs to the neurons take the form of postsynaptic potentials due to incoming spikes. The cross-population recurrent inputs Gnmj(t) are defined by the following: Embedded Image where EPSPk(t) and IPSPk(t) are afferent postsynaptic excitatory and inhibitory potentials, respectively. They are computed by convoluting spike trains with exponential synapses with time constant τm as follows: Embedded Image Both populations are further subjected to endogenous sources of noise ϕe,ij(t), which is assumed to follow spatially and temporally independent Gaussian white noise profiles with fixed amplitude D. Synaptic weights within (Weejk(c), Wiijk(c)) and between (Weijk(c), Wiejk(c)) excitatory and inhibitory populations exhibit sparse exponential profiles (Hellwig, 2000) with connection probability c as follows: Embedded Image where σn,m2 = σe2, σi2 is the range of the excitatory and/or inhibitory interactions and x(k) refers to the spatial location of neuron k in the network. To represent the sparseness of cortical connections, synaptic weights wnmo(c) are set to 0 with probability of 1 − c. The matrix of synaptic weights is plotted in Figure 1A. The network extends in the one-dimensional spatial domain Ω and the neural lattice is built by distributing cells randomly according to a uniform distribution. Given the spatial organization of the network, synaptic interactions between neurons k and l are thus subjected to propagation delays τjk = |x(k) − x(j)|v−1, with v being the axonal conduction velocity, set here to 0.15 m/s. This value for the conduction velocity is realistic for intracortical connections (Hutt et al., 2003). This gives rise to a distribution of propagation latencies plotted in Figure 1B. For the parameters selected, neurons engage spontaneously in synchronized oscillatory alpha activity at a frequency of ∼10 Hz. Sample spiking responses for randomly selected inhibitory and excitatory cells are plotted in Figure 1C, alongside membrane potential traces exhibiting subthreshold oscillations within the alpha band. These oscillations result from the combined influence of propagation delays and slow cellular adaptation, which is known to lead recurrent nonlinear networks in a state of robust synchrony (Lefebvre et al., 2011). This can be seen from the full network activity plotted in Figure 1D, where the spiking of the neurons is modulated by a global rhythm of 10 Hz. Model parameters are summarized in Table 1 below. We note that, throughout, the network parameters were left unchanged and only stimulation settings such as input frequency and amplitude were varied. Being a single-compartment, conductance-based model, the system above shares many similarities with the well known two-population integrate-and-fire networks that have also been shown to sustain global synchronous activity (Brunel, 2000; Buzsáki and Wang, 2012). Key differences lie in the spike initiation mechanism, which in our case adds an additional source of variability in the firing of the modeled neurons, and in the presence of slow inhibition, resulting in slower oscillations frequencies.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Network connectivity and intrinsic activity without stimulation. A, Matrix of synaptic weights that define the connectivity of the network within (Wee, Wii) and between (Wie, Wei) excitatory and inhibitory neural populations. As defined, the connectivity is spatially profiled, where proximal (resp. distal) neurons share stronger (resp. weaker) synaptic connections. B, Distribution of propagation delays in the network. Given a finite conduction velocity of 0.15 m/s, the timing of synaptic interactions is delayed. The distribution shows that the network dynamics are dominated by many short latencies and few longer ones. C, Sample activity of randomly selected excitatory (black) and inhibitory (blue) neurons. The spiking activities Xej and Xik of these neurons are shown above the corresponding membrane fluctuations, which exhibit clear, albeit noisy, subthreshold oscillations. D, Network intrinsic spiking activity without stimulation. Both populations are set in stable synchronous activity with firing rates modulations of ∼10 Hz.

View this table:
  • View inline
  • View popup
Table 1.

Model parameters

Spectral analysis.

Spectral analysis was performed using a fast Fourier transform routine using freely available C++ scripts (Press et al., 2007).

Stimulation.

The stimulation I(t) in Equation 1 is a train of phasic pulses of duration δt = 300 μs of given forcing frequency ωf and intensity If. In the analysis, both frequency and intensity are changed to explore a variety of independent combinations. This stimulus offers a richer spectral signature compared with sinusoidal signals, for which only one frequency is represented. Later on, it is useful to compute the time average of the stimulation pattern μI = <I(t) > T over a given epoch that is given by the following: Embedded Image which also allows also the computation of the variance as Var(I) = <(I(t) − μI)2 > T = IfμI − μI2. We deliberately chose short stimulation epochs (i.e., 1 s) compared with the timescale usually involved during continuous stimulation paradigms in clinical settings (lasting several minutes) because we were investigating the online effect during stimulation, which has been shown to exist for brief stimulation duration of a few seconds (Fröhlich and McCormick, 2010). For the investigation of plastic effects, such as long-term synaptic potentiation (Huerta and Volpe, 2009), longer stimulation duration would be required.

Neuroelectric activity.

In our model, the network neuroelectric activity is determined by a weighted sum over somatic excitatory and inhibitory potentials as follows: Embedded Image where φe,ik are real positive coefficients. Here, we assume that the network fine scale structure is unknown, and thus consider random weights as follows: φe,ik = [0,1]. This is done to take into account, in our calculations, of different sources of observational variability and constraints that we do not model explicitly while avoiding any filtering issues in the resulting spectral calculations. However, specific choices of the φe,ik distributions can be made to increase the similarity of the network activity to physiological signals such as LFPs and EEG (Lindén et al., 2010; Mazzoni et al., 2015).

Spike coherence.

To compute spike coherence between spike trains (Gabbianni and Koch, 1998), we selected randomly two excitatory neurons, and j, within the cortical populations and binned their respective responses over a given time window ΔT. Therefore, for neuron j, Xjbin(kΔT) = 1 if a spike occurred during the interval [kΔT, (k + 1)ΔT] and zero otherwise, where k = T/ΔT. Coherence, θ(ΔT), was then computed as follows: Embedded Image where ∑T denotes the sum over time intervals. For θ(ΔT) = 0, two spike trains are independent and θ(ΔT) = 1 reflects maximum coherence. A nonlinear (resp. linear) dependence of the coherence on the time window will appear whenever synchronous (resp. asynchronous) spiking exists in the data (Wang and Buzsaki, 1996; Hutt and Buhry, 2014). For the purposes of the current analysis, we have chosen ΔT = 10 ms. This value was selected to optimize the coherence difference between original (synchronous) and shuffled spike trains.

Mean coupling strength.

To measure the relative level of recurrence and nonlinearity in the system, we computed the network-average recurrent inputs received by neurons during stimulation epochs and how they changed with respect to different stimulation parameters. The mean coupling strength, referred here as the κ index, is defined by the following: Embedded Image where ≪.≫N,T stands for the mean taken over the population of size N and during a time interval of duration T. For the purposes of the analysis, we have chosen to compute the mean coupling strength over the full duration of the stimulation epochs to average out time fluctuations conveyed by the oscillatory nature of the network output.

Results

To understand the response evoked during stimulation, we investigated the behavior of a generic cortical network expressing alpha synchrony and subjected to periodic stimulation of various combinations of frequencies and amplitudes and observed the evolution of the spectral properties exposed by the network's neuroelectric output. The goal here was not to derive a precise mathematical formulation of the resulting dynamics, but instead to provide a general overview of driven oscillatory responses unfolded in stimulation parameter space. Limited by the dimensionality of possible stimulation waveforms, we focused our attention on repetitive pulsatile stimuli that closely resemble transcranial magnetic stimulation waveforms and/or direct cortical microstimulation paradigms (see Materials and Methods). The stimulation was assumed to be global: all neurons in the network were subjected to the same driving pattern. The aim of the study was to extract the mapping that links input to output spectral signatures. Therefore, we assumed stimulation epochs to be short compared with synaptic potentiation and/or depression such that plasticity had no significant impact on the dynamics. Therefore, our findings characterize instantaneous responses arising due to the recurrent architecture of the network and the rhythmic nature of its ongoing state.

Figure 2 details different response mechanisms that can be evoked during pulsatile forcing in our model network. Without stimulation, the network is found in a synchronized state, exhibiting strongly correlated activity where firing rate fluctuations oscillate at an intrinsic—also called “natural”—frequency of ωo ≅ 10 Hz. Different stimulation forcing amplitudes (If) and frequencies (ωf) lead the network through different response regimes such as resonance, in which intrinsic oscillations are enhanced (Fig. 2B); entrainment, where the systems' dynamics are locked to the stimulation (Fig. 2C); and nonlinear acceleration (NLA), in which the stimulation causes a shift in the system's intrinsic frequency (Fig. 2D).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Responses of synchronized neural population to periodic pulsatile stimulation. The network neuroelectric output A(t) is shown on the left. The associated power spectral density (PSD) is plotted in the middle next to the excitatory neurons spiking activity, which are shown on the right. A, According to the set of selected parameters, the network engages spontaneously into synchronous alpha activity. In absence of stimulation, the endogenous rhythm stabilizes around an unperturbed intrinsic frequency of ωo ≅ 10 Hz, which modulates the noisy spiking activity of the neurons. Once global stimulation is applied to the network, different mechanisms are recruited and explain the observed responses. B, If the stimulation frequency ωf is close to the intrinsic frequency, its harmonics, and/or subharmonics, resonance occurs. Although the network response frequency ω′ remains stable—that is, ω′ ≈ Embedded Image ωo—the associated spectral power increases significantly. This is also reflected in the spiking activity of driven neurons: firing rates increase into tightly synchronized bursts of spike discharges. C, Whenever the stimulation amplitude increases beyond a certain threshold, intrinsic ongoing activity in entrained by the stimulation drive and the network oscillates at a frequency equal to the stimulations, its harmonics, and/or subharmonics; that is, ω′ ≈ Embedded Image ωf. In this example, a strong rhythmic stimulation of ωf ≅ 17.5 Hz was applied and harmonic entrainment can be observed in the 1:2 Arnold tongue as the network peak frequency shifts to ω′ ≅ 8.8 Hz. D, For high-frequency stimulation, nonlinear interactions provide the means of accelerating ongoing synchronous activity. A spectral shift Δ(ωfIf), the magnitude of which depends on stimulus parameters, arises due to a change of the network intrinsic oscillatory properties. In this example, a stimulation of ωf ≅ 97 Hz can be seen to accelerate the baseline activity from ωo ≅ 10 Hz up to ω′ ≅ 15 Hz.

To disambiguate network responses and see how they are shaped by the system's nonlinearity, we analyzed the network response during fully independent trial epochs of 1 s duration. During each epoch, we computed the power spectral density of the network neuroelectric activity (See Materials and Methods). Varying the stimulation frequency alone (i.e., for a fixed stimulation intensity), we found that the different mechanisms detailed in Figure 2 were evoked sequentially, yet in a nontrivial manner. Figure 3A displays the network‘s neuroelectric power spectral density (see Materials and Methods) evolving as a function of increasing stimulation frequency while keeping the stimulation intensity constant. As detailed in Figure 3B, alongside the fundamental frequency peak seen on the diagonal, one can see both harmonic and subharmonic peaks. For lower frequencies, the stimulation was found to tune the peak frequency expressed by synchronous neurons via a sequence of subharmonic entrainments. For high-frequency stimulation, the network was found in a state of fast oscillations in which ongoing cyclic activity was not entrained to the stimulation, but instead was accelerated with respect to the baseline intrinsic frequency.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Power spectrum of the network neuroelectric output A(t) during stimulation of fixed intensity but variable frequency. A, Increasing the frequency of the stimulation (ωf) shifts the peak frequency (ω′) and alters the spectral profile of the network's activity. The stimulation fundamental frequency is seen aligned on the diagonal. In addition, both harmonics (Hm:nω) and subharmonic (SHn:mω) of the stimulation frequency are visible (above and below the diagonal, respectively). For low-frequency forcing, the peak frequency is close to the system's intrinsic frequency ωo. However, entrainment occurs as ωf increases further, where the network peak frequency successively tracks the stimulation frequency ωf and its subharmonics. Once ωf exceeds 50 Hz, the peak frequency stabilizes to an accelerated frequency. B, Schematic representation of the spectral patterns seen in A. Shifts in the network peak frequency (ω′) are plotted in red. Subharmonic resonance can be observed as indicated by the black lines below the diagonal. Here, If = 1.5. In Figure 4, one can see how these patterns unfold as stimulation intensity is also changed.

How do the peak frequency and power change when stimulation intensity is varied? To better understand how the spectral patterns seen in Figure 3 unfold, we systematically varied both If and ωf and identified in each case the network peak frequency (ω′) and peak power at that same frequency. In Figure 4, multiple spectral patterns can be seen in both the power and frequency plots, indicating a plurality of network responses. For lower stimulation frequencies, the system's responses were expectedly found to be shaped by entrainment and resonance. Indeed, as seen in Figure 4A where the power of the peak frequency is plotted, one can see characteristic resonance peaks (R) aligned at multiples of the network intrinsic frequency (ωo) where the stimulation triggered the amplification of the endogenous alpha rhythm. This effect spanned all stimulus intensities. Harmonic resonances can also be seen where peaks appear at integer multiples of ωo (e.g., nωo) with gradually decreasing amplitudes. Nonlinear effects can also be seen readily. Resonance peaks, which would be vertically aligned at every multiple of the systems' intrinsic frequency in the case of a linear oscillator, are here arched. The effect, which is more salient for higher harmonics, is due to a shift in network intrinsic frequency due to the action of stimulation. In addition, one subharmonic resonance peak can be observed at ∼5 Hz. The different resonances expressed by the network are labeled in Figure 4B.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Spectral mapping between input stimulation parameters and the network neuroelectric output A(t). Stimulation frequency (ωf) and amplitude (If) were varied systematically while the network response was analyzed. A, Spectral peak power expressed by the network as a function of stimulation parameters. Characteristic resonance peaks, vertically aligned, are present at the system's intrinsic frequency (ωo ≅ 10 Hz) and at every integer multiple of it. Subharmonic resonance can also be seen at 5 Hz. Through these resonances, the system's output power is magnified due to the proximity of the stimulation frequency to the intrinsic frequency. The peaks, however, gradually vanish as the stimulation frequency increases >50 Hz, leading to responses of relative weak power for all remaining stimulation conditions. B, Labeled resonance (R) curves seen in A, shown along the subharmonic resonance at 5 Hz (SR). C, Network peak frequency (ω′) as a function of stimulation parameters. Arnold tongues representing entrainment are aligned vertically and indicate regimes where the system's frequency is locked to the stimulation frequency. Harmonic entrainment can also be seen. Entrainment occurs mainly for frequencies close the system's intrinsic frequency and its harmonics and vanishes as the stimulation frequency increases. For higher frequencies and amplitudes, however, nonlinear acceleration of intrinsic oscillatory activity can be observed. Such acceleration occurs for higher stimulation frequencies, where ω′ is seen to transit gradually to frequencies in the beta band by the action of forcing alone. D, Labeled entrainment (E) and NLA regimes as seen in C. Note that Arnold tongues overlap each other. The horizontal dashed lines at If = 1.5 seen in all panels represent the specific example detailed in Figure 3 in which the frequency, but not the intensity, was varied (within the interval 0–50 Hz). One can then see how peak power and frequency are shaped as stimulation intensity is also varied.

Other patterns emerge in Figure 4C, where frequency of the peak frequency is plotted with respect to varying stimulation parameters. For weaker stimulation frequencies, characteristic entrainment (E) patterns, also called Arnold tongues (Jensen et al., 1983), dominate the system's response: the network oscillatory behavior is locked to the stimulation frequency; that is: ω′ ≅ ωf. Subharmonic entrainment also occurs as the stimulation frequency increases: the network is entrained to oscillate at frequencies that are subharmonics of the driving frequency. Although similar to resonance peaks, the mechanism involved in generating these Arnold tongues is fundamentally different. We also note that. whereas the entrainment effect scales with stimulation intensity, it is nonetheless limited to driving frequencies that are close to ωo and/or its harmonics (Fig. 4A,B) and vanishes at higher stimulation frequencies.

For stimulation frequencies >50 Hz, another mechanism emerged and was found to dominate the response of the network, the dynamics of which enter a state of increasing recurrence. In this regime, high-frequency pulses increase the network response frequency via an NLA of the endogenous rhythmic activity. The intrinsic network frequency is shifted, leading to a gradual transition from the alpha (8–12 Hz) to the beta (12–30 Hz) band. In contrast to high-intensity/low-frequency stimulation, where resonance and entrainment dominate, NLA occurs in response to stimuli that have more subtle effects on the individual firing rates of driven neurons. Despite high forcing frequencies and amplitudes, NLA triggers rhythmic responses expressing peak power levels comparable to those seen in the unperturbed state (Fig. 2A) and dominates the peak frequency mapping at high frequencies. Visible entrainment and NLA regimes are labeled in Figure 4D.

Close scrutiny of the network firing rates during stimulation epochs also revealed the influence of stimulation on the spiking activity of individual neurons. We computed the normalized mean firing rates of excitatory neurons by first computing the network mean spiking activity in every stimulation condition and then comparing it with resting-state values. As seen in Figure 5A, for weaker stimulation frequencies (i.e., <50 Hz), neural discharge rates are strongly driven by network resonances: neurons' spiking activities track stimulation-induced changes in network peak power. Therefore, firing rates are strongly tied to the stimulation signal. However, for higher stimulation frequencies (i.e., >50 Hz), notably where NLA occurs, the network's mean firing rates still increase despite the weak neuroelectric power expressed (cf. Fig. 4A).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Firing rates and correlations as a function of varying stimulation parameters. A, The rhythmic stimulation was found to modulate individual firing rates non-monotonically. Changes in stimulation parameters triggered increases in spiking activity that mirrored the resonance peaks pattern found in Figure 3A. For higher frequencies and/or amplitude, a regime of high spiking activity can be found despite the absence of significant spectral power found in the same region. The mean firing rate of an excitatory neuron, averaged over a stimulation epoch of 1 s, is plotted. B, Spike coherence was also found to reflect resonance-induced modulations. Resonance was thus recruited to generate responses with highly correlated firing. Spike coherence was found to be low for high-frequency and/or high-amplitude stimuli despite strong spiking activity.

Therefore, to disambiguate the contributions of firing rates, power, and peak frequency to the emergence of synchronous dynamics in presence of stimulation, we further computed the response spike coherence between neurons as stimulation parameters were varied systematically. As seen in Figure 5B, the shadow of the system's resonance peaks were also found, indicating that rhythmic stimulation triggers, not only high firing rates via resonance, but also strongly correlated states, which is consistent with reported links between expressed neural correlations and the magnitude of spiking activity (De La Rocha et al., 2007). In contrast, weak to no spike coherence was found between the neuron's spike trains for stimuli of high frequency (Fig. 5B).

This stimulus-driven phenomenon echoes previous results that revealed the nonlinear mapping between stochastic input statistics and the spectral features of nonlinear recurrent networks (Lefebvre et al., 2015). Building on the timescale separation that subsists amid slow endogenous activity and high-frequency fluctuations, synergistic interactions between stimulation-induced fluctuations and network nonlinearity support a gain control mechanism in which the network susceptibility is enhanced with respect to baseline. In this framework, the stimulation recruits the network's recurrent nonlinearity and distorts its oscillatory properties, an effect that further scales with stimulation intensity; that is, ω′ ≈ ωo + Δ(If, ωf). To support this perspective, we calculated the stimulation time-average, [μ]I (see Materials and Methods) as stimulation settings were varied. As can be seen from Figure 6A, the stimulation mean intensity expectedly increases with both stimulation amplitude and frequency: maximal recruitment of network recurrent nonlinearity should thus mirror the μI distribution. To verify this, we computed in Figure 6B the mean coupling strength index κ (see Materials and Methods) as stimulation parameters were varied systematically. The κ index gauges the average intensity of recurrent interactions across the network and scales with network-induced intrinsic activity, as opposed to stimulation-driven activity; therefore, it is an indirect measure of the system's effective recurrence. The mean coupling strength was found to be shaped by stimulation distinctively compared with the neuron's firing rates and/or spike coherence: the κ index distribution correlated strongly with the NLA pattern shown in Figure 4C and did not build on from resonance and entrainment patterns despite being modulated by those as well. More importantly, maximal recruitment of network recurrence was achieved alongside maximal stimulation intensity, indicating that nonlinear effects support the oscillatory shifts observed in this region of parameter space. Together, these results suggest that, in addition to resonance and entrainment, rhythmic stimulation can also be used to recruit the intrinsic nonlinearity of cortical networks to shape the oscillatory features displayed by synchronous populations.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

A shift in recurrence mediates the nonlinear acceleration observed for high-frequency stimulation. A, Stimulation time-average μI (see Materials and Methods) as stimulation settings are varied. As expected, μI scales with stimulation frequency and amplitude. B, The mean coupling strength κ (see Materials and Methods) follows the same trend and scales with stimulation parameters. Maximal recruitment of recurrent interactions is achieved with peak stimulation mean intensity and thus underlies the nonlinear acceleration of intrinsic oscillations.

Discussion

There is growing evidence demonstrating that noninvasive neurostimulation engages neural networks by perturbing endogenous oscillatory brain activity. Purported mechanisms include the generation of bursts of time-locked gamma oscillations (Huerta and Volpe, 2009) or the disruption and/or entrainment of endogenous alpha cyclic dynamics (Thut et al., 2011) that lead to quantifiable alterations in cognitive performance (Romei et al., 2011) and perception (Chanes et al., 2013; Cecere et al., 2015). Moreover, the influence of stimulation has further been reported to be state dependent (Neuling et al., 2013). The current quantitative assessment of neurostimulation's impact is oftentimes limited to frequency ranges that are close to the system's natural frequency in which resonance plays a dominant role. Here, we have explored the response of a model network of spiking neurons to various combinations of forcing intensity and frequencies to characterize better the plurality of network responses evoked.

Our results confirm that the engagement of nonlinear spiking networks to low-frequency forcing (i.e., <50 Hz) is shaped by resonance and entrainment. In this regime, an enhancement of the amplitude of the intrinsic oscillation (resonance) can be observed reliably if the stimulation frequency is close to the intrinsic frequency or its harmonics and subharmonics. For increased stimulation amplitudes, even if the network is stimulated at frequencies near the intrinsic frequency, the network dynamics are locked to the stimulation frequency (entrainment): this switch from intrinsic to driving frequency is mediated by an increase in forcing amplitude. The network was also found to exhibit subharmonic entrainment; that is, when the driving frequency is a harmonic of the intrinsic frequency, the network responds with a frequency that is subharmonic to the driving frequency and closest to the intrinsic frequency. This phenomenon has been observed previously in human EEG data when subjects were stimulated with different frequencies of flickering light and so-called steady-state visual evoked potentials (SSVEPs) were recorded (Herrmann, 2001). For example, at a stimulation frequency of 80 Hz, the visual cortex responded with a 10 Hz SSVEP oscillation; that is, a 1/8 subharmonic to the driving frequency. Therefore, and despite the simplicity of our model, our findings are consistent with recent tACS experiments showing that sinusoidal electric current applied to the human scalp resulted in entrainment of ongoing alpha oscillations (Helfrich et al., 2014). Our results also show good accordance with findings in animal models in which multiple driving frequencies were found to recruit the underlying Arnold entrainment tongues (Ali et al., 2013).

In contrast, our results indicate that the response dynamics to high-frequency stimulation (i.e., >50 Hz) builds on a different mechanism: network synchronous oscillations were found to be accelerated with respect to baseline through a nonlinear shift of the intrinsic frequency. Scaled by the stimulation intensity, the effect was characteristic of responses with amplified recurrent interactions, supporting a synergetic stability transition known to affect driven nonlinear networks generically (Lefebvre et al., 2013, 2015). Our result thus extend previous theoretical work (Lefebvre et al., 2015) and demonstrate that oscillatory nonlinear shifts can also be observed in sparse networks of spiking neurons using high-frequency stimuli. These computational predictions are also consistent with experimental findings showing such increases in peak alpha frequency during stimulation to flicker flashes of high frequency (Herrmann, 2001; Roberts and Robinson, 2012; cf. Fig. 1A). In addition, our results support the idea that high-frequency, sustained stimulation of cortical networks can shape the intrinsic alpha peak frequency and power via an alternative mechanism that could constitute a substitute and/or complementary strategy to engage cyclic neural dynamics.

At first glance, one might argue that our results do not compare to numerous experiments investigating direct cortical or transcranial brain stimulation because we used impulse trains rather than sinusoidal currents. It has been shown, however, that stimulating the visual cortex with either sinusoidally or rectangularly modulated light both resulted in sinusoidal SSVEPs (Dreyer and Herrmann, 2015). It is therefore plausible to assume that our findings can be used to better understand the human and animal findings on electric stimulation of brain tissue. In addition, our results replicate and extend the physiological findings of a recent SSVEP study that found a 1:1 Arnold tongue in the vicinity of the intrinsic frequency of ∼10 Hz in the human EEG (Notbohm et al., 2016). However, the number of parameter variations is limited in human and animal experiments, so network simulations that scan a broader range of stimulation frequencies and intensities are necessary.

Many studies that have used tACS have reported after-effects of enhanced amplitudes of ongoing brain oscillations after the end of rhythmical stimulation. It has been argued that such after-effects are due to synaptic plasticity because they do not appear at short stimulation durations (Strüber et al., 2015). Our simulation study did not model synaptic plasticity and thus did not investigate poststimulation after-effects. Nonetheless, our model provided a rich array of potential responses that could recruit and modulate plastic processes differentially. Future work should shed light on the stimulation parameters that are required to achieve such after-effects, which are likely to be critical if transcranial brain stimulation is to be used in treatment of diseases with unbalanced brain oscillations.

Footnotes

  • This work was supported by the Natural Sciences and Engineering Research Council of Canada (J.L.), the Swiss National Science Foundation (Grant PZ00P1_148186 to S.I. and Grants 320030-149982 and the National Centre of Competence in Research project SYNAPSY, The Synaptic Bases of Mental Disease, Project 51AU40-125759 to M.M.M.), the European Research Council under the European Union's Seventh Framework Programme [(FP7/2007–2013)/ERC Grant 257253 to A.H.], and the German Research Foundation (Grant HE3353/8-1).

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Prof. Jérémie Lefebvre, Experimental Psychology Laboratory, Department of Psychology, Krembil Discovery Tower, 60 Leonard Avenue, Toronto, Ontario M5T 2S8, Canada. jeremie.lefebvre{at}uhnresearch.ca

References

  1. ↵
    1. Alagapan S,
    2. Schmidt SL,
    3. Lefebvre J,
    4. Hadar E,
    5. Shin HW,
    6. Fröhlich F
    (2016) Modulation of cortical oscillations by low-frequency direct cortical stimulation is state-dependent. PLoS Biol 14:e1002424, doi:10.1371/journal.pbio.1002424, pmid:27023427.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Ali MM,
    2. Sellers KK,
    3. Fröhlich F
    (2013) Transcranial alternating current stimulation modulates large-scale cortical network activity by network resonance. J Neurosci 33:11262–11275, doi:10.1523/JNEUROSCI.5867-12.2013, pmid:23825429.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Amari S
    (1977) Dynamics of pattern formation in lateral inhibition type neural fields. Biol Cybern 27:77–87, doi:10.1007/bf00337259, pmid:911931.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Başar E
    (2012) A review of alpha activity in integrative brain function: fundamental physiology, sensory coding, cognition and pathology. Int J Psychophysiol 86:1–24, doi:10.1016/j.ijpsycho.2012.07.002, pmid:22820267.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Başar E,
    2. Güntekin B
    (2012) A short review of alpha activity in cognitive processes and in cognitive impairment. Int J Psychophysiol 86:25–38, doi:10.1016/j.ijpsycho.2012.07.001, pmid:22801250.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Brunel N
    (2000) Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons. J Comput Neurosci 8:183–208.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Buzsáki G,
    2. Wang XJ
    (2012) Mechanisms of gamma oscillations. Annu Rev Neurosci 35:203–225, doi:10.1146/annurev-neuro-062111-150444, pmid:22443509.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Cecere R,
    2. Rees G,
    3. Romei V
    (2015) Individual differences in alpha frequency drive crossmodal illusory perception. Curr Biol 25:231–235, doi:10.1016/j.cub.2014.11.034, pmid:25544613.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Chanes L,
    2. Quentin R,
    3. Tallon-Baudry C,
    4. Valero-Cabré A
    (2013) Causal frequency-specific contributions of frontal spatiotemporal patterns induced by non-invasive neurostimulation to human visual performance. J Neurosci 33:5000–5005, doi:10.1523/JNEUROSCI.4401-12.2013, pmid:23486970.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. de la Rocha J,
    2. Doiron B,
    3. Shea-Brown E,
    4. Josic K,
    5. Reyes A
    (2007) Correlation between neural spike trains increases with firing rate. Nature 448:802–806, doi:10.1038/nature06028.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Dreyer AM,
    2. Herrmann CS
    (2015) Frequency-modulated steady-state visual evoked potentials: a new stimulation method for brain-computer interfaces. J Neurosci Methods 241:1–9, doi:10.1016/j.jneumeth.2014.12.004, pmid:25522824.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Engel AK,
    2. Singer W
    (2001) Temporal binding and the neural correlates of sensory awareness. Trends Cogn Sci 5:16–25, doi:10.1016/S1364-6613(00)01568-0, pmid:11164732.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Engel AK,
    2. Fries P,
    3. Singer W
    (2001) Dynamic predictions: oscillations and synchrony in top-down processing. Nat Rev Neurosci 2:704–716, doi:10.1038/35094565, pmid:11584308.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Fröhlich F
    (2015) Experiments and models of cortical oscillations as a target for non-invasive brain stimulation. Prog Brain Res 222:41–73, doi:10.1016/bs.pbr.2015.07.025, pmid:26541376.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Fröhlich F,
    2. McCormick DA
    (2010) Endogenous electric fields may guide neocortical network activity. Neuron 67:129–143, doi:10.1016/j.neuron.2010.06.005, pmid:20624597.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Gabbianni A,
    2. Koch C
    (1998) Methods in neuronal modeling, Principles of spike train analysis (MIT, Cambridge, MA), pp 313–360.
  17. ↵
    1. Gutmann B,
    2. Mierau A,
    3. Hülsdünker T,
    4. Hildebrand C,
    5. Przyklenk A,
    6. Hollmann W,
    7. Strüder HK
    (2015) Effects of physical exercise on individual resting state EEG alpha peak frequency. Neural Plast 2015:717312, doi:10.1155/2015/717312, pmid:25759762.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Haegens S,
    2. Cousijn H,
    3. Wallis G,
    4. Harrison PJ,
    5. Nobre AC
    (2014) Inter- and intra-individual variability in alpha peak frequency. Neuroimage 92:46–55, doi:10.1016/j.neuroimage.2014.01.049, pmid:24508648.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Helfrich RF,
    2. Schneider TR,
    3. Rach S,
    4. Trautmann-Lengsfeld SA,
    5. Engel AK,
    6. Herrmann CS
    (2014) Entrainment of brain oscillations by transcranial alternating current stimulation. Curr Biol 24:333–339, doi:10.1016/j.cub.2013.12.041, pmid:24461998.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Hellwig B
    (2000) A quantitative analysis of the local connectivity between pyramidal neurons in layers 2/3 of the rat visual cortex. Biol Cybern 82:111–121, doi:10.1007/PL00007964, pmid:10664098.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Henry MJ,
    2. Herrmann B,
    3. Obleser J
    (2014) Entrained neural oscillations in multiple frequency bands comodulate behavior. Proc Natl Acad Sci U S A 111:14935–14940, doi:10.1073/pnas.1408741111, pmid:25267634.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Herrmann CS
    (2001) Human EEG responses to 1–100 Hz flicker: resonance phenomena in visual cortex and their potential correlation to cognitive phenomena. Exp Brain Res 137:346–353, doi:10.1007/s002210100682, pmid:11355381.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Huerta PT,
    2. Volpe BT
    (2009) Transcranial magnetic stimulation, synaptic plasticity and network oscillations. J Neuroeng Rehabil 6:7, doi:10.1186/1743-0003-6-7, pmid:19254380.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Hülsdünker T,
    2. Mierau A,
    3. Strüder HK
    (2015) Higher balance task demands are associated with an increase in individual alpha peak frequency. Front Hum Neurosci 9:695, doi:10.3389/fnhum.2015.00695, pmid:26779005.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Hutt A,
    2. Buhry L
    (2014) Study of GABAergic extra-synaptic tonic inhibition in single neurons and neural populations by traversing neural scales: application to propofol-induced anaesthesia. J Comput Neurosci 37:417–437, doi:10.1007/s10827-014-0512-x, pmid:24976146.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Hutt A,
    2. Bestehorn M,
    3. Wennekers T
    (2003) Pattern formation in intracortical neuronal fields. Network 14:351–368, doi:10.1088/0954-898X_14_2_310, pmid:12790189.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Jadi MP,
    2. Sejnowski TJ
    (2014) Cortical oscillations arise from contextual interactions that regulate sparse coding. Proc Natl Acad Sci U S A 111:6780–6785, doi:10.1073/pnas.1405300111, pmid:24742427.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Jensen MH,
    2. Bak P,
    3. Bohr T
    (1983) Complete devil's staircase, fractal dimension and universality of mode-locking structure in the circle map. Phys Rev Lett 50:1637–1639, doi:10.1103/PhysRevLett.50.1637.
    OpenUrlCrossRef
  29. ↵
    1. Kajikawa Y,
    2. Schroeder CE
    (2015) Generation of field potentials and modulation of their dynamics through volume integration of cortical activity. J Neurophysiol 113:339–351, doi:10.1152/jn.00914.2013, pmid:25274348.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Klimesch W
    (1999) EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis. Brain Res Rev 29:169–195, doi:10.1016/S0165-0173(98)00056-3, pmid:10209231.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Koch C
    (1999) Biophysics of computation (OUP, Oxford).
  32. ↵
    1. Lakatos P,
    2. Karmos G,
    3. Mehta AD,
    4. Ulbert I,
    5. Schroeder CE
    (2008) Entrainment of neuronal oscillations as a mechanism of attentional selection. Science 320:110–113, doi:10.1126/science.1154735, pmid:18388295.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Lefebvre J,
    2. Hutt A
    (2013) Additive noise quenches delay-induced oscillations. Europhys Lett 102:60003.
    OpenUrlCrossRef
  34. ↵
    1. Lefebvre J,
    2. Longtin A,
    3. LeBlanc VG
    (2011) Neural adaptation facilitates oscillatory responses to static inputs in a recurrent network of ON and OFF cells. J Comput Neurosci 31:73–86, doi:10.1007/s10827-010-0298-4, pmid:21170577.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Lefebvre J,
    2. Hutt A,
    3. Knebel JF,
    4. Whittingstall K,
    5. Murray MM
    (2015) Stimulus statistics shape oscillations in nonlinear recurrent neural networks. J Neurosci 35:2895–2903, doi:10.1523/JNEUROSCI.3609-14.2015, pmid:25698729.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Lindén H,
    2. Pettersen KH,
    3. Einevoll GT
    (2010) Intrinsic dendritic filtering gives low-pass power spectra of local field potentials. J Comput Neurosci 29:423–444, doi:10.1007/s10827-010-0245-4, pmid:20502952.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Lorincz ML,
    2. Kékesi KA,
    3. Juhász G,
    4. Crunelli V,
    5. Hughes SW
    (2009) Temporal framing of thalamic relay-mode firing by phasic inhibition during the alpha rhythm. Neuron 63:683–696, doi:10.1016/j.neuron.2009.08.012, pmid:19755110.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Mazzoni A,
    2. Lindén H,
    3. Cuntz H,
    4. Lansner A,
    5. Panzeri S,
    6. Einevoll GT
    (2015) Computing the local field potential (LFP) from integrate-and-fire network models. PLoS Comput Biol 11:e1004584, doi:10.1371/journal.pcbi.1004584, pmid:26657024.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Neuling T,
    2. Rach S,
    3. Herrmann CS
    (2013) Orchestrating neuronal networks: sustained after-effects of transcranial alternating current stimulation depend upon brain states. Front Hum Neurosci 7:161, doi:10.3389/fnhum.2013.00161, pmid:23641206.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Notbohm A,
    2. Kurths J,
    3. Herrmann CS
    (2016) Modification of brain oscillations via rhythmic light stimulation provides evidence for entrainment but not for superposition of event-related responses. Front Hum Neurosci 10:10, doi:10.3389/fnhum.2016.00010, pmid:26869898.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Podvalny E,
    2. Noy N,
    3. Harel M,
    4. Bickel S,
    5. Chechik G,
    6. Schroeder CE,
    7. Mehta AD,
    8. Tsodyks M,
    9. Malach R
    (2015) A unifying principle underlying the extracellular field potential spectral responses in the human cortex. J Neurophysiol 114:505–519, doi:10.1152/jn.00943.2014, pmid:25855698.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Prescott SA,
    2. Sejnowski TJ
    (2008) Spike-rate coding and spike-time coding are affected oppositely by different adaptation mechanisms. J Neurosci 28:13649–13661, doi:10.1523/JNEUROSCI.1792-08.2008, pmid:19074038.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Press WH,
    2. Teukolsky SA,
    3. Vetterling WT,
    4. Flannery BP
    (2007) Numerical recipes: the art of scientific computing (Cambridge University, Cambridge), Ed 3.
  44. ↵
    1. Reato D,
    2. Rahman A,
    3. Bikson M,
    4. Parra LC
    (2010) Low-intensity electrical stimulation affects network dynamics by modulating population rate and spike timing. J Neurosci 30:15067–15079, doi:10.1523/JNEUROSCI.2059-10.2010, pmid:21068312.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Roberts JA,
    2. Robinson PA
    (2012) Quantitative theory of driven nonlinear brain dynamics. Neuroimage 62:1947–1955, doi:10.1016/j.neuroimage.2012.05.054, pmid:22652022.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Romei V,
    2. Driver J,
    3. Schyns PG,
    4. Thut G
    (2011) Rhythmic TMS over Parietal cortex links distinct brain frequencies to global versus local visual processing. Curr Biol 21:334–337, doi:10.1016/j.cub.2011.01.035, pmid:21315592.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Samaha J,
    2. Postle BR
    (2015) The speed of alpha-band oscillations predicts the temporal resolution of visual perception. Curr Biol 25:2985–2990, doi:10.1016/j.cub.2015.10.007, pmid:26526370.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Samaha J,
    2. Bauer P,
    3. Cimaroli S,
    4. Postle BR
    (2015) Top-down control of the phase of alpha-band oscillations as a mechanism for temporal prediction. Proc Natl Acad Sci U S A 112:8439–8444, doi:10.1073/pnas.1503686112, pmid:26100913.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Siegel M,
    2. Donner TH,
    3. Engel AK
    (2012) Spectral fingerprints of large-scale neuronal interactions. Nat Rev Neurosci 13:121–134, doi:10.1038/nrn3137.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Singer W,
    2. Gray CM
    (1995) Visual feature integration and the temporal correlation hypothesis. Annu Rev Neurosci 18:555–586, doi:10.1146/annurev.ne.18.030195.003011, pmid:7605074.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Strüber D,
    2. Rach S,
    3. Neuling T,
    4. Herrmann CS
    (2015) On the possible role of stimulation duration for after-effects of transcranial alternating current stimulation. Front Cell Neurosci 9:311, doi:10.3389/fncel.2015.00311, pmid:26321912.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Thut G,
    2. Schyns PG,
    3. Gross J
    (2011) Entrainment of perceptually relevant brain oscillations by non-invasive rhythmic stimulation of the human brain. Front Psychol 2:170, doi:10.3389/fpsyg.2011.00170, pmid:3142861.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Thut G,
    2. Miniussi C,
    3. Gross J
    (2012) The functional importance of rhythmic activity in the brain. Curr Biol 22:R658–R663, doi:10.1016/j.cub.2012.06.061, pmid:22917517.
    OpenUrlCrossRefPubMed
  54. ↵
    1. Varela F,
    2. Lachaux J-P,
    3. Rodriguez E,
    4. Martinerie J
    The brainweb: phase synchronization and large-scale integration. Nat Rev Neurosci 2:229–239.
  55. ↵
    1. Vlahou EL,
    2. Thurm F,
    3. Kolassa IT,
    4. Schlee W
    (2014) Resting-state slow-wave power, healthy aging and cognitive performance. Sci Rep 4:5101, doi:10.1038/srep05101, pmid:24869503.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Wang XJ,
    2. Buzsáki G
    (1996) Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network model. J Neurosci 16:6402–6413, pmid:8815919.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Whittington MA,
    2. Traub RD,
    3. Jefferys JG
    (1995) Synchronized oscillations in interneuron networks driven by metabotropic glutamate receptor activation. Nature 373:612–615, doi:10.1038/373612a0, pmid:7854418.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Witkowski M,
    2. Garcia-Cossio E,
    3. Chander BS,
    4. Braun C,
    5. Birbaumer N,
    6. Robinson SE,
    7. Soekadar SR
    (2015) Mapping entrained brain oscillations during transcranial alternating current stimulation (tACS). NeuroImage doi:10.1016/j.neuroimage.2015.10.024.
    OpenUrlCrossRef
Back to top

In this issue

The Journal of Neuroscience: 36 (19)
Journal of Neuroscience
Vol. 36, Issue 19
11 May 2016
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Shaping Intrinsic Neural Oscillations with Periodic Stimulation
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Shaping Intrinsic Neural Oscillations with Periodic Stimulation
Christoph S. Herrmann, Micah M. Murray, Silvio Ionta, Axel Hutt, Jérémie Lefebvre
Journal of Neuroscience 11 May 2016, 36 (19) 5328-5337; DOI: 10.1523/JNEUROSCI.0236-16.2016

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Shaping Intrinsic Neural Oscillations with Periodic Stimulation
Christoph S. Herrmann, Micah M. Murray, Silvio Ionta, Axel Hutt, Jérémie Lefebvre
Journal of Neuroscience 11 May 2016, 36 (19) 5328-5337; DOI: 10.1523/JNEUROSCI.0236-16.2016
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • dynamics
  • networks
  • oscillations
  • spectrum
  • stimulation
  • synchrony

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Choice Behavior Guided by Learned, But Not Innate, Taste Aversion Recruits the Orbitofrontal Cortex
  • Maturation of Spontaneous Firing Properties after Hearing Onset in Rat Auditory Nerve Fibers: Spontaneous Rates, Refractoriness, and Interfiber Correlations
  • Insulin Treatment Prevents Neuroinflammation and Neuronal Injury with Restored Neurobehavioral Function in Models of HIV/AIDS Neurodegeneration
Show more Articles

Systems/Circuits

  • Face-selective patches in marmosets are involved in dynamic and static facial expression processing
  • Sex differences in the impact of electronic nicotine vapor on corticotropin-releasing factor receptor 1 neurons in the mouse ventral tegmental area
  • The 3D Structural Architecture of the Human Hand Area is Non-Topographic
Show more Systems/Circuits
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.