Periodic brain stimulation, such as repetitive transcranial magnetic stimulation (rTMS) and transcranial alternating current stimulation (tACS), can be used to engage cortical rhythms (Fröhlich, 2015). Such findings have raised the fascinating prospect of manipulating rhythmic brain activity in a controlled manner, thereby engaging neural circuits at a functional level to manipulate cognition and to treat disorders of the central nervous system (Cecere et al., 2015; Fröhlich, 2015; Romei et al., 2016; Ezzyat et al., 2017). Yet, the brain is not a passive receiver. Invasive and non-invasive brain stimulation differ in their effect on neuronal dynamics as a function of the state of the targeted network, a robust observation across a variety of stimulation modalities (Neuling et al., 2013; Ruhnau et al., 2016; Alagapan et al., 2016). Indeed, previous studies have reported the seemingly counterintuitive fact that the susceptibility of neural tissue to exogenous control (during stimulation for instance) is increased during task-engaged states (Neuling et al., 2013; Ruhnau et al., 2016; Alagapan et al., 2016). Task-engaged states are generally characterized by higher neuronal firing rates and by population responses that are predominantly asynchronous, and by a suppression of low-frequency electrical brain activity compared to what is seen in the rest condition (Churchland et al., 2010; Rosenbaum et al., 2017). By contrast, the rest state is characterized by strong endogenous oscillations, such as alpha oscillations (8 Hz–12 Hz), which reflect internally driven brain processes (Pfurtscheller et al., 1996; Klimesch, 2012). To compensate for these state-dependent differences in population activity, it has been suggested that stimulation parameters should be calibrated in a closed loop (Boyle and Frohlich, 2013; Lustenberger et al., 2016) in order for stimulation to be optimally effective. But the lack of understanding about the cause of these state-dependent changes and how they interfere with brain stimulation remains one of the main limitations for the development of new paradigms that are based on the selective control of brain rhythms.

We formulated the hypothesis that the dominance of asynchronous activity in the task state — in the form of irregular fluctuations in neural activity — may enhance the susceptibility of cortical neurons to periodic stimulation. Specifically, we asked whether stochastic resonance (SR) was involved. SR is a phenomenon by which the presence of uncorrelated random fluctuations, or so-called ‘noise’, amplifies the salience of a weak periodic signal driving a non-linear system (McDonnell and Abbott, 2009). In absence of noise, a subthreshold signal is too weak to generate an output that reflects the input. For intermediate levels of noise, however, the superposition of the signal and the random fluctuations results in super-threshold — and thus detectable — dynamics. Further increase in noise typically causes the saliency of the signal to decrease. Such a link between SR and state-dependent effects of brain stimulation have been reported experimentally across various modalities (Schwarzkopf et al., 2011; Silvanto and Cattaneo, 2017), in which irregular waveforms have further been shown to improve the effect of rhythmic stimuli (Wuehr et al., 2017). However, it remains unclear how the task-related decrease in alpha oscillations is involved in enhancing the susceptibility of cortical neurons to entrainment by exogenous electric fields.

Computational approaches and modeling are poised to answer many of these fundamental questions, but have received little attention, especially with respect to the clinical applications of brain stimulation. By identifying the network mechanisms involved, modeling can be used to catalyze the development of more efficient stimulation protocols and customized clinical interventions. Using simulations, we have previously used a cortical network model to better understand how periodic stimulation interacts with resting-state alpha (8 Hz–12 Hz) activity, and found that multiple mechanisms where involved (Herrmann et al., 2016). Echoing the theory of non-linear oscillators in physics, stimulations with increasing amplitudes and/or frequencies were found to shape networks responses through different mechanisms: resonance, in which intrinsic oscillations are enhanced by stimuli with frequencies close to the network natural frequency; and entrainment, where the systems’ dynamics become phase locked to the stimulation and thus adopt the same frequency. Entrainment occurs only for specific sets of stimulation amplitudes and frequencies (regions in stimulation parameters space called Arnold tongues). But what happens in the presence of state-induced changes in baseline activity?

To address these questions and to provide predictions about optimal entrainment conditions, we harnessed computational and mathematical techniques to analyze state-dependent effects on a network model of the thalamo-cortical system in the presence of periodic stimulation. We found that increased thalamic drive suppressed endogenous oscillations throughout the thalamo-cortical loop, leading to an asynchronous state characterized by intense spiking activity that was weakly synchronized. We defined regimes of low (with respect to high) thalamic drive and resting (with respect to task-engaged) brain state to represent two externally imposed experimental conditions. In presence of periodic stimulation, the transition from the rest to the task-engaged state was accompanied by an increased susceptibility to entrainment, as demonstrated by amplified power at stimulation frequency and phase alignment between cortical responses and stimulation. Specifically, we found that the oscillatory response of the system switched from alpha oscillations to stimulus-induced activity as thalamic drive was increased, i.e. during a transition from the resting state to the task state. Taken together, these results show that the thalamo-cortical loop implements a gain control mechanism that regulates the robustness of alpha oscillatory activity and, by doing so, modulates cortical susceptibility to rhythmic entrainment by stimulation.

To understand how the effect of periodic stimulation depends on brain state (Alagapan et al., 2016), we examined the hypothesis that cortical susceptibility to entrainment is controlled by sub-cortical (i.e. thalamic) populations that are reciprocally connected with the cortical networks targeted by stimulation. Given the important role of the thalamus in controlling cortical state (Poulet et al., 2012) and resting-state alpha activity (Hughes et al., 2004; Hughes and Crunelli, 2005; Lorincz et al., 2009), we investigated how increased thalamic inputs regulated cortical excitability and the power of endogenous alpha oscillations. By brain state, we here refer to the state of the thalamo-cortical network, which is determined by the afferent sensory drive to thalamus in the model and parallels two different experimental conditions (i.e. rest and task-engaged). We developed a model of the thalamo-cortical network composed of cortical, thalamic relay and reticular populations that exhibit synchronous activity within the alpha band (see Materials and methods). In this model, cortical populations of excitatory and inhibitory cells share both feedforward and feedback connections with the thalamus. Sensory inputs (e.g. visual) are sent through the lateral geniculate nucleus (LGN) to cortical areas for further processing (Poulet et al., 2012), impacting the pattern of activity of cortical neurons. The model structure is illustrated in Figure 1, along with the steady-state responses generated by the different groups of neurons present.

Figure 1

Download asset Open asset

To quantify the impact of thalamic drive on cortical activity and alpha oscillations — first without stimulation — we calculated both cortical and sub-cortical firing rates, and the power of endogenous oscillations as the input to the thalamus was increased (Figure 2). We did this to define quantitatively different brain states, i.e. resting state vs task-engaged state. For low input to the thalamus, all populations across the network were maintained in a highly synchronous state, where recurrent interactions generate strong firing-rate oscillations with an envelope of about 8 Hz. As can be seen in Figure 1B, the activity of the different populations is synchronized, and the phase difference is due to the presence of propagation delays from the thalamus to and back from the cortex.

Figure 2

Download asset Open asset

Increases in thalamic inputs have a destabilizing effect on endogenous oscillations. While increasing thalamic drive also increased firing rates throughout the system, a gradual suppression of synchronous alpha activity could also be observed: the power of endogenous alpha oscillations decreased substantially. This increased thalamic drive also changed the spiking patterns of cortical neurons: spiking activity became irregular and asynchronous, in line with studies showing a transition towards decorrelated cortical dynamics at task-onset (Pfurtscheller et al., 1996; Poulet et al., 2012; Churchland et al., 2010). This negative correlation between alpha power and spike rate is further in line with other experimental findings (Haegens et al., 2011; van Kerkoerle et al., 2014). Taken together, these results demonstrate that changes in thalamic state are sufficient to destabilize alpha oscillations throughout the thalamo-cortical loop, as well as to shape spike correlations. We thus defined rest and task states as limit cases of low (i.e. rest) and high (i.e. task-engaged) thalamic drive, respectively (see Materials and methods), as indicated in Figure 2C. We note that aside from the thalamic input, the network parameters were kept constant. Such a smooth and gradual destabilization of alpha oscillations is suggestive of a noise-induced supercritical Hopf bifurcation (See Materials and methods). According to this well-known mechanism (Horsthemke and Lefever, 1982), noise acts as a gain control parameter by which non-linear oscillations sustained by delayed and recurrent inhibition (as in our model) can be altered and/or suppressed by a change in noise variance (Lefebvre et al., 2015; Hutt et al., 2016). In the process, noise also linearizes the dynamics by smoothing the non-linear response function of threshold systems (Gammaitoni, 1995; McDonnell and Abbott, 2009; Lefebvre et al., 2015). In the context of our model, rest and task-engaged states would then correspond to opposite points above and below the transition point, respectively, where the ‘noise’ is here mediated by increasingly irregular input from thalamus to cortex.

Having defined rest and task-engaged states in our network model and having characterized the influence of thalamic activity on alpha oscillations, we next asked how state-dependent changes in resting state oscillations (the destabilization of alpha oscillations) were impacted by cortical periodic stimulation. To disambiguate the contribution of resonance (the amplification of endogenous oscillations) and entrainment (the phase locking to an externally applied frequency) and to test network responses, we first considered a cortical stimulation with moderate amplitude and frequency of 11 Hz (${\omega }_{stim}$), which does not share any low order harmonic relationship with the endogenous frequency (i.e. ${\omega }_o=$8 Hz). This was done to align our simulation to a previous experimental setting (Alagapan et al., 2016), in order to determine whether we could reproduce the previously observed state-dependent effects. Then, we investigated how spiking activity and cortical oscillations properties changed as thalamic drive was gradually increased in the presence of periodic stimulation (Figure 3). We first computed the correlations in the spiking response amongst excitatory cells in the cortical network, either with 11 Hz stimulation or with no stimulation (sham), see Figure 3A. Without stimulation, thalamic drive suppresses correlated spiking and synchronous activity decreases. By contrast, when stimulation is applied in regimes of weak thalamic drive, correlations are lower than in the sham condition, indicating that stimulation itself interferes with endogenous activity. By increasing thalamic drive, spike-rate correlations gradually increase, indicating that cells are becoming entrained by the stimulation. Taken together, these observations elucidate a gradual transition from recurrent to externally driven dynamics as thalamic input is increased. To explore this further, we plotted, in Figure 3B, the peak power found both at alpha frequency (8 Hz) and at stimulation frequency (here 11 Hz). The power at stimulation frequency did indeed increase during the transition from the rest to the task state and that a shift in oscillatory regime can be observed. Indeed, one can clearly observe a jump in the dominant frequency from 8 Hz (endogenous alpha frequency) to 11 Hz (stimulation frequency) as resting state alpha oscillations are gradually suppressed.

Figure 3

Download asset Open asset

We hypothesized that this enhanced spectral power at the stimulation frequency might occur through SR. SR has been implicated in the enhancement of sub-threshold signals in which random fluctuations — of either external or internal origin — act like a pedestal increasing the sensitivity of stimulated neurons to a given set of low-intensity inputs (Miniussi et al., 2013). Initially formulated in bistable systems with an implicit time-scale (a ‘resonance’), SR is also present in systems with threshold nonlinearities, which we note are present in our model.

To verify whether SR was indeed involved —and to disambiguate the contribution of SR and dithering (McDonnell and Abbott, 2009) — we chose a weak amplitude stimulation and computed the mutual information between stimulation signal and firing-rate response at the stimulation frequency across successive and independent trials, while noise variance was gradually increased across all populations in the network. As shown in Figure 3D, the mutual information was found to be optimized at an intermediate value, indicating that SR — in its most strict definition (McDonnell and Abbott, 2009) — was indeed involved. In contrast to SR, dithering — a similar phenomenon by which the saliency of weak signals is improved by noise — would cause the mutual information to increase monotonically, that is, it would not exhibit a peak at some intermediate value of noise. To confirm the presence of entrainment (to the stimulation and not to another internal oscillation), we performed a phase analysis comparing the phase of the cortical firing rates with the phase of the stimulation across independent trials. The results are shown in Figure 3E and F. There, one can see that weak phase alignment was found in the rest state (Figure 3E), whereas clear phase locking can be seen in the task state (Figure 3F). Our simulations show that cortical spiking responses are phase locked to the stimulation waveform only for intermediate values of thalamic drive. Below this drive, endogenous network activity suppresses cortical entrainment, and beyond this value, irregular noisy fluctuations originating from the thalamic population decreases the saliency of the temporal structure of the stimulation waveform. Taken together, the results above indicate that state-dependent inputs simultaneously amplify stimulation-induced activity while suppressing endogenous rhythmic activity through SR, enabling the stimulation to entrain cortical neurons.

Figure 4 summarizes these findings in the rest and task-engaged states. In the rest state (Figure 4A), the dynamics of the cortical neurons were predominantly governed by internal, recurrent fluctuations. When stimulation is applied in this state (Figure 4B), ongoing alpha oscillations are slightly perturbed but no cortical entrainment occurs. Stimulation has a minimal impact on the activity of cortical neurons: the phasic alternation of network-wide coincident firing and inhibitory recovery greatly reduces susceptibility to entrainment. This absence of phase locking of neuronal activity to the stimulation signal can be most clearly seen from the firing activity of cortical excitatory and inhibitory neurons. The firing activity of these cells occurs by successive bursts at 8 Hz, with little effect from stimulation-induced membrane fluctuations, which remain sub-threshold. The poor susceptibility is also apparent from the relatively small contribution of stimulation power to the overall spectral density distribution of the simulated cortical EEG (Figure 3C). These indicate that the rest state is characterized by highly stable attractor dynamics. By contrast, robust entrainment was observed in the task-engaged state (Figure 4C), as seen from a dominant spectral peak located precisely at the stimulation frequency (Figure 3C). Despite the fact that stimulation amplitude and frequency remained the same, cortical neurons were here fully entrained by the 11 Hz stimulation. The irregular activity that characterizes the task state was here replaced by stimulus-induced fluctuations while neural firing was tightly phased locked to the stimulus dynamics. Spectral power at the stimulation frequency is significantly amplified when compared to that in the rest state, whereas negligible power can be found at the endogenous alpha frequency. This forms yet another indication that the resting state attractor has been destabilized under the action of thalamic drive and the resulting increase in neural noise (Lefebvre et al., 2015; Hutt et al., 2016), also in line with requirements of SR (McDonnell and Abbott, 2009).

Figure 4

Download asset Open asset

But how does this state-dependent susceptibility depend on stimulation parameters? To answer this question, we first fixed stimulation amplitude and varied its frequency between 0 and 50 Hz, while computing the power spectral density in each case (Figure 5). A peak at the stimulation frequency can be seen along the diagonal, but the rest state power spectral density is characterized by a dominant horizontal peak at 8 Hz. According to this analysis, the spectral contribution of stimulation is negligible unless the stimulation frequency is close to the endogenous frequency and/or its harmonics (horizontal lines in Figure 5A), as expected via resonance (Ali et al., 2013; Herrmann et al., 2016). In the task-engaged state, the opposite occurs: reduced power at the endogenous frequency promotes entrainment, and stimulation power increases significantly across the range of stimulation frequency considered.

Figure 5

Download asset Open asset

As a next step, we computed the power-spectral density of the simulated EEG activity of cortical excitatory and inhibitory cells, and systematically measured the peak frequency and power for all combinations of stimulation frequency and amplitudes within a given range. We then identified regions for which the dominant (i.e. peak) frequency was defined either by endogenous oscillations or by the stimulation — thus identifying regimes of entrainment. In Figure 6A, one can see the peak frequency of the simulated cortical EEG activity in the rest state. For most combinations of stimulation amplitudes and frequencies, the peak frequency remains stable and equal to the endogenous alpha frequency (i.e. 8 Hz). However, for higher amplitudes and for stimulation frequencies near 8 Hz, cortical neurons were gradually entrained by the stimulation (1:1 entrainment). A narrow triangular entrainment region (the so-called Arnold tongue, delimited by white dashed lines) (Glass and Belair, 1980; Hunter and Milton, 2003) emerges and gets larger as stimulation amplitude increases. Note that the asymmetrical shape of the Arnold tongue is here due to a stimulation-induced shift of the endogenous frequency (Hutt et al., 2016). The equivalent calculations were made in Figure 6B but for the task-engaged state. Multiple differences compared to the rest state can be observed. The Arnold tongue spans a much larger portion of the stimulus-parameter space, indicating robust entrainment for much weaker stimulation amplitudes. In Figure 6C and D, the power at peak response frequency is plotted. The power found at the endogenous frequency in the rest state outside the Arnold tongue (here also delimited by a white dashed line) is high, and is gradually suppressed as stimulation parameters are changed towards the entrainment region. There, peak power, associated with the stimulation frequency is much smaller (Figure 6C). Under the action of thalamic input in the task state, endogenous alpha oscillations are suppressed compared to those in the rest state, leading to weak power outside the Arnold tongue (Figure 6D). In contrast to the rest state, stimulation at the endogenous frequency enhances the peak power. In sum, the effect of stimulation on peak power when the frequency is close to the resonant frequency (vertical black dashed lines in Figure 6C and D) also depends on the state, that is, the thalamic drive.

Figure 6

Download asset Open asset

What mechanism is responsible for this state-dependence? Mathematical insights can play a key role here in helping us to understand how periodic stimulation interacts with intrinsic limit cycles (i.e. alpha oscillations) and how changes in stability can promote entrainment. In Alagapan et al., 2016), the authors modelled state-dependent entrainment using a population-scale network of cortical neurons interacting through both local and feedback projections in the presence of robust resting state oscillations. To learn more about the mechanism involved, we considered a simplified neural oscillator with delayed feedback (Lefebvre and Hutt, 2013; Lefebvre et al., 2015; Hutt et al., 2016) as a simplified model for the thalamo-cortical network. Networks of neurons that are exposed to delayed feedback and recurrent inhibition commonly display rhythmic activity whose features are tightly linked to input statistics (Dhamala et al. 2001; Doiron et al., 2003; Lefebvre and Hutt, 2013). To better understand the mechanism underlying the state-dependent changes in the amplitude of resting state oscillations, we investigated the dynamics of this simplified oscillator (see Materials and methods) in the presence of periodic forcing (i.e. stimulation) and various levels of noise. In can be shown that as input noise increases during the task state, limit cycle solutions are destabilized through a supercritical Hopf bifurcation. Analyzing the response of this simplified model both above and below the bifurcation threshold, that is, in the vicinity of the point where self-sustained (i.e. intrinsic) oscillations become unstable, we found that limit cycle solutions of this system are entrained by periodic forcing, although their amplitudes are highly sensitive to noise intensity. To see this, we computed resonance curves (Figure 7B) for this oscillator and compared the results for both high and low values of input noise variance. This simple qualitative analysis shows that the linearization induced by noise decreases the amplitude of the resonant solutions and increases the amplitude of all other forcing frequencies by transitioning through a critical state. Spectral clustering, in which the power and/or amplitude of forced solutions are concentrated near the intrinsic resonance, is a state of weak susceptibility to entrainment. The consequences of input-induced linearization in the simple conceptual model above are two-fold: (1) it causes a suppression of resonant oscillations, and (2) it augments the amplitude of non-resonant solutions. This was found to be in agreement with the state-dependency observed in our full spiking model as depicted in Figure 7C.

Figure 7

Download asset Open asset

Alpha activity has been implicated in a wide variety of physiological and cognitive functions (Başar and Basar, 2012; Mierau et al., 2017). Numerous recent studies have demonstrated that periodically modulated electric fields can interfere with alpha oscillations to trigger a measurable effect on stimulus perception (Chanes et al., 2013; Cecere et al., 2015) and task performance by entraining these rhythms (Henry et al., 2014; Samaha and Postle, 2015). However, the efficacy of these approaches has been shown to depend heavily on brain states and to fluctuate according to internally governed dynamics (Zaehle et al., 2010; Miniussi et al., 2013; Ruhnau et al., 2016; Alagapan et al., 2016).

To better understand this phenomenon and to elucidate previous experimental findings (Alagapan et al., 2016), we here combined computational and mathematical approaches and analyzed a computational model of the thalamo-cortical system exhibiting resting state alpha oscillations. The dynamics of our network model were analyzed in two conditions: the rest state, which we have defined as a regime of weak afferent input to the thalamus, and the task-engaged state, where thalamic drive is high to reflect experimentally imposed conditions. We found that a destabilization of alpha oscillations occurred during the transition between rest and task-engaged states, in which both cortical and sub-cortical spiking patterns switched from correlated and synchronous to irregular and asynchronous. Then, applying periodic stimulation to cortical neurons, we found that the susceptibility of the system to external perturbations was strongly controlled by the presence (or absence) of an oscillatory attractor (i.e. the alpha oscillation) and as such depended heavily on state. By comparing the stimulation waveform with the cortical firing-rate responses for weak stimulation amplitudes, we found that SR was involved in improving cortical susceptibility to entrainment. While in the rest state, applied stimulation had a weak impact on the activity of cortical neurons whose dynamics remained locked to internal recurrent oscillations. Arnold tongues (i.e. regions in the stimulation parameter space where network activity is phase locked to the stimulation) were found to be very narrow. By contrast, in the task-engaged state, excitatory and inhibitory cortical neurons were freed from the intrinsic attractor and found to be more susceptible to entrainment, as shown by significantly larger Arnold tongues. Using a simplified neural oscillator model, our mathematical analysis revealed that thalamic input to the cortex, playing a similar role as noise, could control the stability of intrinsic oscillations by provoking a bifurcation, and by doing so, could amplify the response to non-resonant stimuli. Aside from the state-dependent input and/or noise, only stimulation parameters, such as amplitude and frequency, were changed in our study. Taken together, our results show that internal control over attractor states shapes the susceptibility of cortical populations to stimulation and entrainment. Our modelling results further suggest that the thalamo-cortical system behaves like a non-linear oscillator whose susceptibility to periodic forcing is gated by noise through SR and a passage through a bifurcation.

SR is widely appreciated as a mechanism that regulates cortical excitability (Miniussi et al., 2013) and has measurable effects on cognitive performance and perception (van der Groen and Wenderoth, 2016). Our computational results demonstrate the important role of SR, and thus of asynchronous and irregular neural dynamics, in regulating the responses of cortical neurons to stimulation. Using a simplified neural oscillator model, we have shown that noise — playing the role of uncorrelated inputs from the thalamus — can suppress oscillation through a Hopf bifurcation. Practically speaking, this is accomplished by a noise-induced linearization of the response function of the system, a generic effect observed in systems with threshold non-linearities (see McDonnell and Abbott, 2009) and references therein), which can further be shown to impact the frequency and amplitude of oscillatory solutions in non-linear systems (Lefebvre et al., 2015; Hutt et al., 2016). Such linearization also plays a fundamental role in dithering, a mechanism by which noise is added to a signal, prior to its digitalization, to increase its integrity (McDonnell and Abbott, 2009). Formulated in the context of neural systems, irregular and/or uncorrelated fluctuations can thus be used to increase the saliency or fidelity of neural responses. Interestingly, this ‘correlation-induced blockade’ has been previously discussed with respect to alpha oscillations (Klimesch et al., 2007), where large-scale synchrony has been hypothesized to implement a functional inhibition gate through which increased synchronized fluctuations shut down the mutual coding capacity of neural networks (Zohary et al., 1994; Greenberg et al., 2008). In this perspective, we hypothesize that alpha oscillations — as well as other intrinsic attractors — block perturbations of neural activity, triggered either by endogenous (e.g. inter-area signaling) or exogenous (e.g. stimulation) manipulations. Future work will show whether endogenous correlated neural activity, which goes beyond the purely rhythmic activity observed here, could also play a role in regulating susceptibility to entrainment, and more generally in control.

We here note that, in line with previous theoretical (Hutt et al., 2016) and experimental (Mierau et al., 2017) findings, the peak alpha frequency was found to change as a function of thalamic drive between 8 and 8.5 Hz. This can be seen for instance in Figure 3C (top), Figure 5A,B and Figure 6A,B, where the peak frequency is slightly higher in the task-engaged state than in the rest state. This is because inputs from the thalamus to the cortex, whose variance changes as a function of thalamic drive, correspond to a stream of irregular fluctuations (i.e. Poisson noise) that are known to alter the oscillatory properties of non-linear networks (Lefebvre et al., 2015). As here these fluctuations possess stationary statistics (i.e. mean and variance that did not fluctuate in time), no temporal variability was observed in the peak frequency of endogenous oscillations. As such, this non-linear phenomenon, which is well understood mathematically, does not impact our conclusions.

Our results are conceptually connected to experimental findings in studies of tACS. In the task-engaged state, successful entrainment was reflected in a shift of the oscillation frequency to the stimulation frequency. This agrees with the findings of multiple studies of tACS that point towards a successful shift in the alpha peak frequency to the stimulation frequency during stimulation (e.g. Helfrich et al., 2014). Several studies using tACS (see Thut et al., 2017 and references therein) and a recent study (Minami and Amano, 2017) that uses amplitude-modulated tACS (AM-tACS) have employed some type of visual task, ranging from a simple vigilance task to prevent participants from falling asleep all the way to sophisticated visual psychophysics paradigms.

As for any modeling work, a key aspect to discuss is the level of abstraction. The electrical field applied by tACS is the strongest in the cortex due to its closest proximity to the stimulation electrodes. However, any modulation of cortical network dynamics by stimulation must, almost by definition, also alter activity in the regions interconnected with the targeted cortical circuits, including the thalamus. However, dissecting the effects of stimulation on an interconnected network remains a challenge because of the feedback nature of the network dynamics. Here, we have chosen to study a reduced, relatively abstract model of the thalamo-cortical system. By contrast, in a recent study in which we focused on the cellular and synaptic biophysics and their role in shaping the response to stimulation, we examined the response to stimulation of thalamus in isolation (Li et al., 2017). Qualitatively, the findings from the two studies agree as the detailed thalamic model suggested that faster oscillations, typically associated with the task-engaged state, are more responsive to stimulation than the task-negative alpha oscillations.

Our results further raise interesting questions regarding sensory entrainment paradigms (Keitel et al., 2017). In the presence of periodic sensory stimulation, or sensory inputs possessing a rhythmic structure, thalamic drive would also possess a rhythmic structure. It would prove interesting to see how these fluctuations in thalamo-cortical inputs would interact with endogenous oscillations, and whether brain state fluctuations could gate their propagation towards cortical networks. Given the central role of thalamus in sensory gating (e.g. Sherman, 2007; Poulet et al., 2012), sensory rhythmic inputs that reach the thalamus first may have a different effect on thalamo-cortical oscillations than does cortical stimulation. Yet another fundamental difference is the magnitude of the input. Sensory input to the thalamus (e.g., retinal input to the LGN) is strong enough to cause spiking in the thalamus, which is time-locked to the temporal structure of the visual input (Hirsch et al., 2015). Thus such input is super-threshold in comparison to the subthreshold stimulation examined here.

In summary, our results demonstrate that target engagement by stimulation is state-dependent. The fact that the presence of a pronounced endogenous oscillation limits the impact of stimulation in terms of entrainment points towards important questions for the use of such stimulation paradigms in research, and ultimately in the clinic. First, the overall state of the network needs to be measured and documented in studies that investigate brain stimulation for the modulation of brain rhythms. Second, it remains mostly unclear how adjusting the frequency of the alpha oscillation affects cognition and behavior (Romei et al., 2016; Mierau et al., 2017). Yet, a reduced alpha peak frequency is a hallmark of a range of neurological and psychiatric illnesses (e.g. Rossini et al., 2007) and probably an overall marker of brain health (Klimesch, 1999, 2012). In such work, our findings suggest that states of low alpha power represent the ideal state for the application of stimulation, and that modulation of brain oscillations is best achieved in states of low endogenous alpha rhythmic activity. This finding may be important for the next generation of brain stimulation paradigms developed using rational design.

1. 1. Alagapan S
   2. Schmidt SL
   3. Lefebvre J
   4. Hadar E
   5. Shin HW
   6. Frohlich F

   (2016) Modulation of cortical oscillations by low-frequency direct cortical stimulation is state-dependent

   PLOS Biology 14:e1002424.

   https://doi.org/10.1371/journal.pbio.1002424

   • PubMed
   • Google Scholar
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

   Journal of Neuroscience 33:11262–11275.

   https://doi.org/10.1523/JNEUROSCI.5867-12.2013

   • PubMed
   • Google Scholar
3. 1. Başar E
   2. Basar E

   (2012) A review of alpha activity in integrative brain function: fundamental physiology, sensory coding, cognition and pathology

   International Journal of Psychophysiology 86:1–24.

   https://doi.org/10.1016/j.ijpsycho.2012.07.002

   • PubMed
   • Google Scholar
4. 1. Boyle MR
   2. Frohlich F

   (2013) EEG feedback-controlled transcranial alternating current stimulation

   International IEEE EMBS Conference on Neural Engineering pp. 140–143.

   https://doi.org/10.1109/NER.2013.6695891

   • Google Scholar
5. 1. Cecere R
   2. Rees G
   3. Romei V

   (2015) Individual differences in alpha frequency drive crossmodal illusory perception

   Current Biology 25:231–235.

   https://doi.org/10.1016/j.cub.2014.11.034

   • PubMed
   • Google Scholar
6. 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

   Journal of Neuroscience 33:5000–5005.

   https://doi.org/10.1523/JNEUROSCI.4401-12.2013

   • PubMed
   • Google Scholar
7. 1. Churchland MM
   2. Yu BM
   3. Cunningham JP
   4. Sugrue LP
   5. Cohen MR
   6. Corrado GS
   7. Newsome WT
   8. Clark AM
   9. Hosseini P
   10. Scott BB
   11. Bradley DC
   12. Smith MA
   13. Kohn A
   14. Movshon JA
   15. Armstrong KM
   16. Moore T
   17. Chang SW
   18. Snyder LH
   19. Lisberger SG
   20. Priebe NJ
   21. Finn IM
   22. Ferster D
   23. Ryu SI
   24. Santhanam G
   25. Sahani M
   26. Shenoy KV

   (2010) Stimulus onset quenches neural variability: a widespread cortical phenomenon

   Nature Neuroscience 13:369–378.

   https://doi.org/10.1038/nn.2501

   • PubMed
   • Google Scholar
8. 1. Deco G
   2. Jirsa V
   3. McIntosh AR
   4. Sporns O
   5. Kötter R

   (2009) Key role of coupling, delay, and noise in resting brain fluctuations

   Proceedings of the National Academy of Sciences 106:10302–10307.

   https://doi.org/10.1073/pnas.0901831106

   • PubMed
   • Google Scholar
9. 1. Dhamala M
   2. Jirsa VK
   3. Ding M

   (2004) Enhancement of neural synchrony by time delay

   Physical Review Letters 92:074104.

   https://doi.org/10.1103/PhysRevLett.92.074104

   • PubMed
   • Google Scholar
10. 1. Doiron B
    2. Chacron MJ
    3. Maler L
    4. Longtin A
    5. Bastian J

    (2003) Inhibitory feedback required for network oscillatory responses to communication but not prey stimuli

    Nature 421:539–543.

    https://doi.org/10.1038/nature01360

    • PubMed
    • Google Scholar
11. 1. Ezzyat Y
    2. Kragel JE
    3. Burke JF
    4. Levy DF
    5. Lyalenko A
    6. Wanda P
    7. O'Sullivan L
    8. Hurley KB
    9. Busygin S
    10. Pedisich I
    11. Sperling MR
    12. Worrell GA
    13. Kucewicz MT
    14. Davis KA
    15. Lucas TH
    16. Inman CS
    17. Lega BC
    18. Jobst BC
    19. Sheth SA
    20. Zaghloul K
    21. Jutras MJ
    22. Stein JM
    23. Das SR
    24. Gorniak R
    25. Rizzuto DS
    26. Kahana MJ

    (2017) Direct brain stimulation modulates encoding states and memory performance in humans

    Current Biology 27:1251–1258.

    https://doi.org/10.1016/j.cub.2017.03.028

    • PubMed
    • Google Scholar
12. 1. Fröhlich F

    (2015) Experiments and models of cortical oscillations as a target for noninvasive brain stimulation

    Progress in Brain Research 222:41–73.

    https://doi.org/10.1016/bs.pbr.2015.07.025

    • PubMed
    • Google Scholar
13. 1. Gammaitoni L

    (1995) Stochastic resonance and the dithering effect in threshold physical systems

    Physical Review E 52:4691–4698.

    https://doi.org/10.1103/PhysRevE.52.4691

    • PubMed
    • Google Scholar
14. Book

    1. Glass L
    2. Belair J

    (1980)

    Continuation of Arnold tongues in mathematical models of periodically forced biological oscillators

    In: Othmer H. G, editors. Nonlinear Oscillations in Biology and Chemistry. Springer-Verlag. pp. 232–243.

    • Google Scholar
15. 1. Greenberg DS
    2. Houweling AR
    3. Kerr JN

    (2008) Population imaging of ongoing neuronal activity in the visual cortex of awake rats

    Nature Neuroscience 11:749–751.

    https://doi.org/10.1038/nn.2140

    • PubMed
    • Google Scholar
16. 1. Haegens S
    2. Nácher V
    3. Luna R
    4. Romo R
    5. Jensen O

    (2011) α-Oscillations in the monkey sensorimotor network influence discrimination performance by rhythmical inhibition of neuronal spiking

    PNAS 108:19377–19382.

    https://doi.org/10.1073/pnas.1117190108

    • PubMed
    • Google Scholar
17. 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

    Current Biology 24:333–339.

    https://doi.org/10.1016/j.cub.2013.12.041

    • PubMed
    • Google Scholar
18. 1. Hellwig B

    (2000) A quantitative analysis of the local connectivity between pyramidal neurons in layers 2/3 of the rat visual cortex

    Biological Cybernetics 82:111–121.

    https://doi.org/10.1007/PL00007964

    • PubMed
    • Google Scholar
19. 1. Henry MJ
    2. Herrmann B
    3. Obleser J

    (2014) Entrained neural oscillations in multiple frequency bands comodulate behavior

    PNAS 111:14935–14940.

    https://doi.org/10.1073/pnas.1408741111

    • PubMed
    • Google Scholar
20. 1. Herrmann C
    2. Ionta S
    3. Hutt A
    4. Murray MM
    5. Lefebvre J

    (2016) Shaping intrinsic oscillatory neural activity with periodic stimulation

    Journal of Neuroscience 36:5328–5337.

    https://doi.org/10.1523/JNEUROSCI.0236-16.2016

    • Google Scholar
21. 1. Hirsch JA
    2. Wang X
    3. Sommer FT
    4. Martinez LM

    (2015) How inhibitory circuits in the thalamus serve vision

    Annual Review of Neuroscience 38:309–329.

    https://doi.org/10.1146/annurev-neuro-071013-014229

    • PubMed
    • Google Scholar
22. Book

    1. Horsthemke W
    2. Lefever R

    (1982)

    Noise-Induced Tran-sitions: Theory and Applications in Physics, Chemistry and Biology

    Berlin: Springer.

    • Google Scholar
23. 1. Hughes SW
    2. Crunelli V

    (2005) Thalamic mechanisms of EEG alpha rhythms and their pathological implications

    The Neuroscientist 11:357–372.

    https://doi.org/10.1177/1073858405277450

    • PubMed
    • Google Scholar
24. 1. Hughes SW
    2. Lörincz M
    3. Cope DW
    4. Blethyn KL
    5. Kékesi KA
    6. Parri HR
    7. Juhász G
    8. Crunelli V

    (2004) Synchronized oscillations at alpha and theta frequencies in the lateral geniculate nucleus

    Neuron 42:253–268.

    https://doi.org/10.1016/S0896-6273(04)00191-6

    • PubMed
    • Google Scholar
25. 1. Hunter JD
    2. Milton J

    (2003)

    Epilepsy as a Dynamic Disease. Biological and Medical Physics Series

    Controlling neural synchrony with periodic and aperiodic stimuli, Epilepsy as a Dynamic Disease. Biological and Medical Physics Series, Heidelberg, Springer.

    • Google Scholar
26. 1. Hutt A
    2. Bestehorn M
    3. Wennekers T

    (2003) Pattern formation in intracortical neuronal fields

    Network: Computation in Neural Systems 14:351–368.

    https://doi.org/10.1088/0954-898X_14_2_310

    • Google Scholar
27. 1. Hutt A
    2. Mierau A
    3. Lefebvre J

    (2016) Dynamic Control of Synchronous Activity in Networks of Spiking Neurons

    Plos One 11:e0161488.

    https://doi.org/10.1371/journal.pone.0161488

    • PubMed
    • Google Scholar
28. Preprint

    1. Keitel C
    2. Benwell CSY
    3. Thut G
    4. Gross J

    (2017) No changes in parieto-occipital alpha during neural phase locking to visual quasi-periodic theta-, alpha-, and beta-band stimulation

    bioRxiv.

    https://doi.org/10.1101/219766

    • Google Scholar
29. 1. Klimesch W
    2. Sauseng P
    3. Hanslmayr S

    (2007) EEG alpha oscillations: the inhibition-timing hypothesis

    Brain Research Reviews 53:63–88.

    https://doi.org/10.1016/j.brainresrev.2006.06.003

    • PubMed
    • Google Scholar
30. 1. Klimesch W

    (1999) EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis

    Brain Research Reviews 29:169–195.

    https://doi.org/10.1016/S0165-0173(98)00056-3

    • PubMed
    • Google Scholar
31. 1. Klimesch W

    (2012) α-band oscillations, attention, and controlled access to stored information

    Trends in Cognitive Sciences 16:606–617.

    https://doi.org/10.1016/j.tics.2012.10.007

    • PubMed
    • Google Scholar
32. 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

    Journal of Neuroscience 35:2895–2903.

    https://doi.org/10.1523/JNEUROSCI.3609-14.2015

    • PubMed
    • Google Scholar
33. 1. Lefebvre J
    2. Hutt A

    (2013) Additive noise quenches delay-induced oscillations

    EPL 102:60003.

    https://doi.org/10.1209/0295-5075/102/60003

    • Google Scholar
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

    Journal of Computational Neuroscience 31:73–86.

    https://doi.org/10.1007/s10827-010-0298-4

    • PubMed
    • Google Scholar
35. 1. Li G
    2. Henriquez CS
    3. Fröhlich F

    (2017) Unified thalamic model generates multiple distinct oscillations with state-dependent entrainment by stimulation

    PLOS Computational Biology 13:e1005797.

    https://doi.org/10.1371/journal.pcbi.1005797

    • PubMed
    • Google Scholar
36. 1. Linden H
    2. Pettersen KH
    3. Einevoll GT

    (2010) Intrinsic dendritic filtering gives low-pass power spectra of local field potentials

    Journal of Computational Neuroscience 29:423–444.

    https://doi.org/10.1007/s10827-010-0245-4

    • PubMed
    • Google Scholar
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.

    https://doi.org/10.1016/j.neuron.2009.08.012

    • PubMed
    • Google Scholar
38. 1. Lustenberger C
    2. Boyle MR
    3. Alagapan S
    4. Mellin JM
    5. Vaughn BV
    6. Fröhlich F

    (2016) Feedback-controlled transcranial alternating current stimulation reveals a functional role of sleep spindles in motor memory consolidation

    Current Biology 26:2127–2136.

    https://doi.org/10.1016/j.cub.2016.06.044

    • PubMed
    • Google Scholar
39. 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 Computational Biology 11:e1004584.

    https://doi.org/10.1371/journal.pcbi.1004584

    • PubMed
    • Google Scholar
40. 1. McDonnell MD
    2. Abbott D

    (2009) What is stochastic resonance? Definitions, misconceptions, debates, and its relevance to biology

    PLoS Computational Biology 5:e1000348.

    https://doi.org/10.1371/journal.pcbi.1000348

    • PubMed
    • Google Scholar
41. 1. Mierau A
    2. Klimesch W
    3. Lefebvre J

    (2017) State-dependent alpha peak frequency shifts: Experimental evidence, potential mechanisms and functional implications

    Neuroscience 360:146–154.

    https://doi.org/10.1016/j.neuroscience.2017.07.037

    • PubMed
    • Google Scholar
42. 1. Minami S
    2. Amano K

    (2017) Illusory jitter perceived at the frequency of alpha oscillations

    Current Biology 27:2344–2351.

    https://doi.org/10.1016/j.cub.2017.06.033

    • PubMed
    • Google Scholar
43. 1. Miniussi C
    2. Harris JA
    3. Ruzzoli M

    (2013) Modelling non-invasive brain stimulation in cognitive neuroscience

    Neuroscience & Biobehavioral Reviews 37:1702–1712.

    https://doi.org/10.1016/j.neubiorev.2013.06.014

    • PubMed
    • Google Scholar
44. 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

    Frontiers in Human Neuroscience 7:161.

    https://doi.org/10.3389/fnhum.2013.00161

    • PubMed
    • Google Scholar
45. 1. Pfurtscheller G
    2. Stancák A
    3. Neuper C

    (1996) Event-related synchronization (ERS) in the alpha band--an electrophysiological correlate of cortical idling: a review

    International Journal of Psychophysiology 24:39–46.

    https://doi.org/10.1016/S0167-8760(96)00066-9

    • PubMed
    • Google Scholar
46. 1. Poulet JF
    2. Fernandez LM
    3. Crochet S
    4. Petersen CC

    (2012) Thalamic control of cortical states

    Nature Neuroscience 15:370–372.

    https://doi.org/10.1038/nn.3035

    • PubMed
    • Google Scholar
47. Book

    1. Press WH
    2. Teukolsky SA
    3. Vetterling WT
    4. Flannery BP

    (2007)

    Numericalrecipes: The Art of Scientific Computing (3)

    Cambridge: Cambridge University.

    • Google Scholar
48. 1. Roberts JA
    2. Robinson PA

    (2008) Modeling absence seizure dynamics: implications for basic mechanisms and measurement of thalamocortical and corticothalamic latencies

    Journal of Theoretical Biology 253:189–201.

    https://doi.org/10.1016/j.jtbi.2008.03.005

    • PubMed
    • Google Scholar
49. 1. Romei V
    2. Thut G
    3. Silvanto J

    (2016) Information-based approaches of noninvasive transcranial brain stimulation

    Trends in Neurosciences 39:782–795.

    https://doi.org/10.1016/j.tins.2016.09.001

    • PubMed
    • Google Scholar
50. 1. Rosenbaum R
    2. Smith MA
    3. Kohn A
    4. Rubin JE
    5. Doiron B

    (2017) The spatial structure of correlated neuronal variability

    Nature Neuroscience 20:107–114.

    https://doi.org/10.1038/nn.4433

    • PubMed
    • Google Scholar
51. 1. Rossini PM
    2. Rossi S
    3. Babiloni C
    4. Polich J

    (2007) Clinical neurophysiology of aging brain: from normal aging to neurodegeneration

    Progress in Neurobiology 83:375–400.

    https://doi.org/10.1016/j.pneurobio.2007.07.010

    • PubMed
    • Google Scholar
52. 1. Ruhnau P
    2. Neuling T
    3. Fuscá M
    4. Herrmann CS
    5. Demarchi G
    6. Weisz N

    (2016) Eyes wide shut: Transcranial alternating current stimulation drives alpha rhythm in a state dependent manner

    Scientific Reports 6:27138.

    https://doi.org/10.1038/srep27138

    • PubMed
    • Google Scholar
53. 1. Samaha J
    2. Postle BR

    (2015) The speed of alpha-band oscillations predicts the temporal resolution of visual perception

    Current Biology 25:2985–2990.

    https://doi.org/10.1016/j.cub.2015.10.007

    • PubMed
    • Google Scholar
54. 1. Schwarzkopf DS
    2. Silvanto J
    3. Rees G

    (2011) Stochastic resonance effects reveal the neural mechanisms of transcranial magnetic stimulation

    Journal of Neuroscience 31:3143–3147.

    https://doi.org/10.1523/JNEUROSCI.4863-10.2011

    • PubMed
    • Google Scholar
55. 1. Sherman SM

    (2007) The thalamus is more than just a relay

    Current Opinion in Neurobiology 17:417–422.

    https://doi.org/10.1016/j.conb.2007.07.003

    • PubMed
    • Google Scholar
56. 1. Silvanto J
    2. Cattaneo Z

    (2017) Common framework for "virtual lesion" and state-dependent TMS: The facilitatory/suppressive range model of online TMS effects on behavior

    Brain and Cognition 119:32–38.

    https://doi.org/10.1016/j.bandc.2017.09.007

    • PubMed
    • Google Scholar
57. 1. Slézia A
    2. Hangya B
    3. Ulbert I
    4. Acsády L

    (2011) Phase advancement and nucleus-specific timing of thalamocortical activity during slow cortical oscillation

    Journal of Neuroscience 31:607–617.

    https://doi.org/10.1523/JNEUROSCI.3375-10.2011

    • PubMed
    • Google Scholar
58. 1. Thut G
    2. Bergmann TO
    3. Fröhlich F
    4. Soekadar SR
    5. Brittain JS
    6. Valero-Cabré A
    7. Sack AT
    8. Miniussi C
    9. Antal A
    10. Siebner HR
    11. Ziemann U
    12. Herrmann CS

    (2017) Guiding transcranial brain stimulation by EEG/MEG to interact with ongoing brain activity and associated functions: A position paper

    Clinical Neurophysiology 128:843–857.

    https://doi.org/10.1016/j.clinph.2017.01.003

    • PubMed
    • Google Scholar
59. 1. Thut G
    2. Miniussi C
    3. Gross J

    (2012) The functional importance of rhythmic activity in the brain

    Current Biology 22:R658–R663.

    https://doi.org/10.1016/j.cub.2012.06.061

    • PubMed
    • Google Scholar
60. 1. van der Groen O
    2. Wenderoth N

    (2016) Transcranial random noise stimulation of visual cortex: Stochastic resonance enhances central mechanisms of perception

    The Journal of Neuroscience 36:5289–5298.

    https://doi.org/10.1523/JNEUROSCI.4519-15.2016

    • PubMed
    • Google Scholar
61. 1. van Kerkoerle T
    2. Self MW
    3. Dagnino B
    4. Gariel-Mathis MA
    5. Poort J
    6. van der Togt C
    7. Roelfsema PR

    (2014) Alpha and gamma oscillations characterize feedback and feedforward processing in monkey visual cortex

    PNAS 111:14332–14341.

    https://doi.org/10.1073/pnas.1402773111

    • PubMed
    • Google Scholar
62. 1. Wuehr M
    2. Boerner JC
    3. Pradhan C
    4. Decker J
    5. Jahn K
    6. Brandt T
    7. Schniepp R

    (2017) Stochastic resonance in the human vestibular system - Noise-induced facilitation of vestibulospinal reflexes

    Brain Stimulation.

    https://doi.org/10.1016/j.brs.2017.10.016

    • PubMed
    • Google Scholar
63. 1. Zaehle T
    2. Rach S
    3. Herrmann CS

    (2010) Transcranial alternating current stimulation enhances individual alpha activity in human EEG

    PLoS One 5:e13766.

    https://doi.org/10.1371/journal.pone.0013766

    • PubMed
    • Google Scholar
64. 1. Zohary E
    2. Shadlen MN
    3. Newsome WT

    (1994) Correlated neuronal discharge rate and its implications for psychophysical performance

    Nature 370:140–143.

    https://doi.org/10.1038/370140a0

    • PubMed
    • Google Scholar

Author details

1. Jérémie Lefebvre

   Krembil Research Institute, Toronto, Canada

   Contribution

   Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing

   For correspondence

   jeremie.lefebvre@uhnresearch.ca

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-0369-4565

2. Axel Hutt

   FE12 - Data Assimilation, Deutscher Wetterdienst, Offenbach am Main, Germany

   Contribution

   Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-0041-7431

3. Flavio Frohlich

   Department of Psychiatry, University of North Carolina at Chapel Hill, Chapel Hill, United States

   Contribution

   Conceptualization, Resources, Supervision, Writing—original draft, Project administration, Writing—review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-3724-5621

• 1,950

  views

• 374

  downloads

• 37

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.