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 in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • 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
Research Articles, Systems/Circuits

No Evidence for Entrainment: Endogenous Gamma Oscillations and Rhythmic Flicker Responses Coexist in Visual Cortex

Katharina Duecker, Tjerk P. Gutteling, Christoph S. Herrmann and Ole Jensen
Journal of Neuroscience 4 August 2021, 41 (31) 6684-6698; DOI: https://doi.org/10.1523/JNEUROSCI.3134-20.2021
Katharina Duecker
1Centre for Human Brain Health, School of Psychology, University of Birmingham, Birmingham B15 2SA, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Katharina Duecker
Tjerk P. Gutteling
1Centre for Human Brain Health, School of Psychology, University of Birmingham, Birmingham B15 2SA, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tjerk P. Gutteling
Christoph S. Herrmann
2Department of Psychology, Faculty VI-Medicine and Health Sciences, Carl-von-Ossietzky University of Oldenburg, Oldenburg 26129, 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
Ole Jensen
1Centre for Human Brain Health, School of Psychology, University of Birmingham, Birmingham B15 2SA, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ole Jensen
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Over the past decades, numerous studies have linked cortical gamma oscillations (∼30–100 Hz) to neurocomputational mechanisms. Their functional relevance, however, is still passionately debated. Here, we asked whether endogenous gamma oscillations in the human brain can be entrained by a rhythmic photic drive >50 Hz. Such a noninvasive modulation of endogenous brain rhythms would allow conclusions about their causal involvement in neurocognition. To this end, we systematically investigated oscillatory responses to a rapid sinusoidal flicker in the absence and presence of endogenous gamma oscillations using magnetoencephalography (MEG) in combination with a high-frequency projector. The photic drive produced a robust response over visual cortex to stimulation frequencies of up to 80 Hz. Strong, endogenous gamma oscillations were induced using moving grating stimuli as repeatedly done in previous research. When superimposing the flicker and the gratings, there was no evidence for phase or frequency entrainment of the endogenous gamma oscillations by the photic drive. Unexpectedly, we did not observe an amplification of the flicker response around participants' individual gamma frequencies (IGFs); rather, the magnitude of the response decreased monotonically with increasing frequency. Source reconstruction suggests that the flicker response and the gamma oscillations were produced by separate, coexistent generators in visual cortex. The presented findings challenge the notion that cortical gamma oscillations can be entrained by rhythmic visual stimulation. Instead, the mechanism generating endogenous gamma oscillations seems to be resilient to external perturbation.

SIGNIFICANCE STATEMENT We aimed to investigate to what extent ongoing, high-frequency oscillations in the gamma-band (30–100 Hz) in the human brain can be entrained by a visual flicker. Gamma oscillations have long been suggested to coordinate neuronal firing and enable interregional communication. Our results demonstrate that rhythmic visual stimulation cannot hijack the dynamics of ongoing gamma oscillations; rather, the flicker response and the endogenous gamma oscillations coexist in different visual areas. Therefore, while a visual flicker evokes a strong neuronal response even at high frequencies in the gamma-band, it does not entrain endogenous gamma oscillations in visual cortex. This has important implications for interpreting studies investigating the causal and neuroprotective effects of rhythmic sensory stimulation in the gamma-band.

  • entrainment
  • flicker
  • frequency tagging
  • gamma oscillations
  • magnetoencephalography
  • neuronal oscillations

Introduction

Cortical gamma oscillations have been repeatedly linked to the formation of neuronal ensembles through synchronization of spiking activity in rodents and primates (Eckhorn et al., 1988; Gray and Singer, 1989; Gray et al., 1992; Wehr and Laurent, 1996; Brosch et al., 2002), including humans (Tallon et al., 1995; Müller et al., 1997; Rodriguez et al., 1999; Hoogenboom et al., 2006). Accordingly, they have been ascribed a supporting role for neuronal computations within populations (Singer and Gray, 1995; Singer, 1999, 2009; Von der Malsburg, 1999; Engel et al., 2001; Nikolić et al., 2013) as well as interregional functional connectivity (Bressler, 1990; Varela et al., 2001; Fries et al., 2007). Indeed, numerous studies have been able to link gamma oscillations in the human brain to cognitive processes and perception (for review, see Başar-Eroglu et al., 1996; Herrmann and Mecklinger, 2001; Jensen et al., 2007; Tallon-Baudry, 2009; Uhlhaas et al., 2009), whereas anomalous gamma-band activity has been associated with impaired cognition and awareness, as in, e.g., autism spectrum disorder, schizophrenia and Alzheimer's dementia (for review, see Herrmann and Demiralp, 2005; Uhlhaas and Singer, 2006; Uhlhaas et al., 2009; Traub and Whittington, 2010; Grützner et al., 2013).

In this study, we aimed to entrain, i.e., synchronize, gamma oscillations in the human visual cortex to a rhythmic photic drive at frequencies above 50 Hz. Stimulation at such high frequencies has recently been applied in rapid frequency tagging (RFT) protocols, to investigate spatial attention (Zhigalov et al., 2019) and audiovisual integration in speech (Drijvers et al., 2021), with minimal visibility of the flicker. The ability to non-invasively modulate gamma rhythms would allow to study their causal role in neuronal processing and cognition, as well as their therapeutic potential, as recently proposed by (Iaccarino et al., 2016; Adaikkan et al., 2019).

It is widely accepted that rhythmic inhibition imposed by inhibitory interneurons forms the backbone of neuronal gamma oscillations (Traub et al., 1996; Lozano-Soldevilla et al., 2014; for review, see Bartos et al., 2007; Buzsáki and Wang, 2012). Indeed, Cardin et al. (2009) demonstrate evidence for resonance, i.e., a targeted amplification, in the gamma-band, in response to optogenetic stimulation of GABAergic interneurons, but not when driving excitatory pyramidal cells (also see Tiesinga, 2012). Here, we ask whether a rapid photic flicker can hijack human visual gamma oscillations; a positive outcome would suggest that visual stimulation can modulate pyramidal-inhibitory-network-gamma (PING) activity. To this end, we designed a paradigm that embraces the definition of resonance and entrainment as stated in dynamical systems theory. While neuroscientific studies widely rely on this terminology (Hutcheon and Yarom, 2000; Schwab et al., 2006; Notbohm et al., 2016; Lakatos et al., 2019), the prerequisites of entrainment are often not sufficiently accounted for, as pointed out by Helfrich et al. (2019). Entrainment requires the presence of a self-sustained oscillator that synchronizes to an external drive (Pikovsky et al., 2003; Thut et al., 2011). This synchronization is reflected by a convergence of the frequency and phase of the endogenous oscillator to the driving force (Pikovsky et al., 2003). Similarly, resonance is reflected by periodic responses to a rhythmic drive and an amplification of individually preferred rhythms, but does not require the presence of self-sustained oscillations per se (Pikovsky et al., 2003; Helfrich et al., 2019). Indeed, studies on photic stimulation at a broad range of frequencies (Herrmann, 2001; Gulbinaite et al., 2019), including the alpha-band (Notbohm et al., 2016), have provided evidence for both resonance and entrainment in the visual system (also see Rager and Singer, 1998, for resonance phenomena in cat visual cortex).

In this study, oscillatory responses to photic stimulation from 52 to 90 Hz were investigated in the presence and absence of visually induced gamma oscillations. In the flicker condition, a rhythmic flicker was applied to a circular, invisible patch. In the flicker&gratings condition, the flicker was superimposed on moving grating stimuli that have been shown to reliably induce strong, narrow-band gamma oscillations (Hoogenboom et al., 2006, 2010; Van Pelt and Fries, 2013). These oscillations reflect individual neuronal dynamics (Hoogenboom et al., 2006; Van Pelt and Fries, 2013) and have been shown to propagate to downstream areas in the visual hierarchy (Buffalo et al., 2011; Bosman et al., 2012; Bastos et al., 2015; Michalareas et al., 2016). Therefore, we will use the terms induced and endogenous gamma oscillations interchangeably in the following. We chose moving grating stimuli to elicit narrow-band endogenous gamma oscillations since more complex stimuli induce a broad-band gamma response which might not reflect oscillations (Hermes et al., 2015a,b).

We expected the visual system to resonate to frequencies close to the endogenous gamma rhythm elicited by the gratings, as well as a synchronization of the gamma oscillations and the rhythmic flicker. As we will demonstrate, the moving gratings did generate strong endogenous gamma oscillations, and the photic drive did produce robust responses at frequencies up to 80 Hz. However, to our great surprise, there was no evidence that the rhythmic stimulation entrains endogenous gamma oscillations.

Materials and Methods

Experimental procedure and apparatus

The MEG data were recorded using a MEGIN Triux system housed in a magnetically shielded room (MSR; Vacuumschmelze GmbH & Co.). Neuromagnetic signals were acquired from 204 orthogonal planar gradiometers and 102 magnetometers at 102 sensor positions. Horizontal and vertical EOG, the cardiac ECG signals, stimulus markers as well as luminance changes recorded by a photodiode were acquired together with the neuromagnetic signal. The data were sampled at 1000 Hz and lowpass filtered online at 330 Hz. Structural magnetic resonance images (MRIs), for later co-registration with the MEG data, were acquired using a 3 Tesla Siemens MAGNETOM Prisma whole-body scanner (Siemens AG; TE = 2 ms, and TR = 2 s). For two subjects, the T1-weighted images obtained in previous experiments, using a 3 Tesla Philips Achieva Scanner (Philips North America Corporation), were used (scanned at the former Birmingham University Imaging Center). Participants were invited to two separate sessions during which the MEG data and the anatomic images were acquired, respectively. Whenever possible, the MEG recording preceded the MRI scan; otherwise, the MEG session was scheduled at least 48 h after the MRI session to avoid any residual magnetization from the MRI system. Volunteers were requested to remove all metal items (e.g., jewelry) before entering the MSR. To enable later co-registration between MRI and MEG data, four to five head-position-indicator (HPI) coils were attached to the participants' foreheads. Along with the position of the coils, three fiducial landmarks (nasion, left and right tragus) and over 200 head-shape samples were digitized using a Polhemus Fastrak (Polhemus, Colchester, USA). Following the preparation, the participants were seated in upright position under the dewar, with orientation set to 60°. The MEG experiment consisted of fifteen blocks lasting 4 min 30 s each. Participants were offered breaks every ∼20 min but remained seated. At the beginning of each of these recording blocks, subjects were instructed to sit with the top and backside of their head touching the sensor helmet. The positions of the HPI coils relative to the sensors was gathered at the beginning of each recording block, but not continuously. The MEG experiment lasted ∼75 min in total.

Rapid photic stimulation

Stimuli were presented using a Propixx lite projector (VPixx Technologies Inc) which allows refresh rates of up to 1440 Hz. To achieve this high-frequency mode, the projector separates the screen (initial resolution: 1920 × 1080 pixels) into quadrants and treats them as separate frames, resulting in a display resolution of 960 × 540 pixels. The RGB color codes for each quadrant, red, green and blue, are converted to a gray scale, separately for each frame and color, and presented consecutively within one refresh interval. The twelve frames are presented at a refresh rate of 120 Hz, resulting in 12 × 120 Hz = 1440 Hz. This approach allows to drive the luminance of each pixel with high temporal precision, allowing for smooth sinusoidal modulations, reducing unwanted harmonics (see Fig. 1C,D). In this study, we applied rapid rhythmic stimulation at frequencies ranging from 52 to 90 Hz in 2-Hz increments.

Experimental paradigm

Stimuli were created in MATLAB 2017a (The MathWorks) and presented using the Psychophysics Toolbox version 3 (Brainard, 1997).

Conditions

The experiment consisted of two conditions that will be referred to as the flicker and the flicker&gratings condition, respectively. Each trial began with a 1-s interval, in which a central white fixation cross was presented on a dark gray background. In the flicker trials, a photic drive in the shape of a circular patch of diameter 2.62° was presented for 2 s. Therefor, the patch's luminance was modulated sinusoidally at frequencies between 52 and 90 Hz (Fig. 1A). To minimize the visibility of the flicker, the mean luminance of the patch was matched with the background (33% luminance, 213.5 cd/m2, RGB [84 84 84]). Frequencies were randomized and balanced across trials. The patch was centered on the fixation cross, such that it was presented both foveally and parafoveally. Each trial ended with a 2-s interval in which only the fixation cross was presented. In the flicker&gratings condition, the baseline interval was followed by a 2-s presentation of a moving grating stimulus that has been shown to reliably elicit gamma oscillations in visual cortex (Hoogenboom et al., 2006, 2010; Muthukumaraswamy and Singh, 2013; Tan et al., 2016). The stimulus was the same size as the patch (2.62°) and had a spatial frequency of 9.1 rings/° (see Fig. 1B), the individual rings' width was 0.11°. The rings contracted toward the center of the screen with a velocity of 0.56°/s, i.e., ∼4 cycles/s. In the subsequent 2-s interval, the gratings were flickered at the respective frequencies, by sinusoidally modulating the luminance of the entire stimulus with each screen refresh. The trial concluded with a 2-s interval in which the concentric moving circles remained on the screen without photic stimulation. To keep the overall brightness of the stimulation similar between conditions, the luminance of the circular patch in the flicker condition ranged from 0% to 66% (of the projector's maximum), while the brightness of the gratings in the flicker&gratings ranged from 33% to 99%, with an average luminance of 214.5 cd/m2 during the presentation of the flicker. The resulting contrast between the gray and black rings, of 66%, has been previously demonstrated to induce clearly identifiable gamma oscillations (Self et al., 2016). The flicker was replicated in the lower right corner of the screen, to acquire the stimulation signal with a photodiode. The rationale of this design was to investigate whether and how the resonance properties of the visual system change when an endogenous gamma oscillator in visual cortex is activated; and whether the flicker response modulates the ongoing oscillatory activity. Studying these two phenomena in the flicker&gratings condition required a characterization of both the gamma oscillations and flicker response in isolation. The former was achieved by presenting the gratings without the flicker. To extract the flicker response, we aimed to avoid any gamma-band activity in visual cortex. This was implemented by applying the flicker to a texture-free, invisible patch. Given the filter properties of the visual system (for review, see Cormack, 2005), we were further interested in identifying an upper limit of the frequencies inducing reliable responses. As we expect these results to guide future studies employing the rapid flicker for frequency tagging, we chose an invisible patch to avoid any confounds by response enhancement, e.g., by object-based attention or figure-ground segregation (Self et al., 2016).

Task and time course

Participants were kept vigilant by performing a simple visual detection task that required them to respond to a 45° rotation of the fixation cross at the center of the screen, which occurred once every minute (Zaehle et al., 2010). Data including the target and/or the responses were discarded and not considered in the analysis. The rotation took place after a trial in the majority, i.e., 60%, of the cases. The remaining 40% of rotations took place at any point during a trial. The experiment was divided into 15 blocks of 4.5 min, resulting in a recording time of 75 min in total. The 40 frequency × condition combinations were presented once in each block, in randomized order, resulting in a total of 15 trials per flicker frequency and condition. To minimize the amount of trials rejected by eye-blink artifacts, 3-s breaks, indicated by a motivating catchphrase or happy face on the screen, were incorporated every five trials, i.e., every 25–35 s. Participants were instructed to use these breaks to rest their eyes.

Participants

This project was reviewed and approved by the local Ethics Committee at the University of Birmingham. Thirty-one students of the University of Birmingham participated in the experiment. One experimental session was terminated prematurely due to the participant not being cooperative, resulting in a sample of 30 participants (15 female), aged 25.7 ± 3.4 years. This sample size was decided on based on a conceptually similar study investigating entrainment of neuronal α oscillations by Notbohm et al. (2016). All volunteers declared not to have had a history of neuropsychiatric or psychological disorders, reported to be medication-free and had normal or corrected-to-normal vision. For safety reasons, volunteers with metal items inside their bodies were excluded at the selection state. Before taking part in the study, participants gave informed consent, in accordance with the declaration of Helsinki, to both the MEG recording and the MRI scan and were explicitly apprised of their right to abort the experiment at any point. The reimbursement amounted to £15/h. To allow analysis of flicker responses at frequencies with a sufficient distance to the individual gamma frequency (IGF; see Results, Identifying IGFs) i.e., ±6 Hz, eight participants were excluded because of their IGF being below 58 Hz. Thus, the data of 22 participants were included in the following analyses (11 female; mean age 25.7 years).

MEG data analysis

Analyses were performed in MATLAB 2017a and 2019b (The MathWorks) using the FieldTrip toolbox (Oostenveld et al., 2011).

Sensor analysis

At the sensor level, the analysis was confined to the planar gradiometer signals, as these provided the best signal-to-noise ratio. The sensor positions relative to the HPI coils were loaded in from the data files and averaged for each participant.

MEG preprocessing

Trials containing the target or button presses were excluded. The data were read into MATLAB as 5- and 7-s trials for the flicker and flicker&gratings conditions, respectively. Artefactual sensors were identified visually during and after the recordings for each participant, and interpolated with the data of their neighboring sensors (0–2 sensors per participant). The individual trials were linearly detrended. Trials containing head movements and/or multiple eye blinks were discarded using a semi-automatic approach. An ICA approach (“runica” implemented in FieldTrip) was used to project out cardiac signals, eye blinks and eye movement.

Time-frequency representation (TFR) of power

TFRs of power were calculated using a sliding time-window approach (ΔT = 0.5 s; 0.05-s steps). A Hanning taper (0.5 s) was applied before the Fourier-transform. This approach induced spectral smoothing of ±3 Hz. Relative power change in response to the stimulation, i.e., the moving grating and/or the photic drive, was calculated as: Pnormalized=PstimPbase−1,(1) with Pstim being the power during stimulation and Pbase being the power in the baseline interval. The baseline interval was 0.75–0.25 s before the onset of the flicker (flicker condition) or the moving grating stimulus (flicker&gratings condition).

IGF

The frequency band of the oscillatory activity elicited in response to the moving grating stimulus was identified individually per participant. TFRs of power were calculated over the first 3 seconds of each trial, that is, the baseline interval and the presentation of the moving grating, in the flicker&gratings condition and averaged over trials. The normalized power was averaged over the 0.25- to 1.75-s interval, and the frequency bin with the maximum relative power was considered the IGF. For each participant, the four to six gradiometers with the strongest gamma response to the moving gratings were selected as the sensors of interest (SOIs).

Phase-locking

The average phase-synchrony between the photodiode (recording the visual flicker) and the neuromagnetic signal at the SOI was quantified by the phase-locking value (PLV; Lachaux et al., 1999; Bastos and Schoffelen, 2016) calculated using a 0.5-s sliding window multiplied with a Hanning taper. The phases of both signals were calculated from Fourier transformations, applied to the tapered segments. The PLV was computed separately for each frequency × condition combination: PLV=1n|∑n=1Nexp(jθ(t,n))|,(2) where θ(t,n)=ϕm(t,n)−ϕp(t,n) is the phase difference between the MEG (m) and the photodiode (p) signal at time bin t in trial n (see Lachaux et al., 1999; p. 195).

Phase difference as a measure of entrainment

Additionally, we investigated changes in phase difference between the photodiode and neuromagnetic signal over time for flicker frequencies of IGF ± 6 Hz, to identify intervals of strong synchrony, so-called phase plateaus. MEG and photodiode signals (ΔT = 3 cycles = 3/fflickers) were convolved with a complex Hanning taper using the sliding time window approach. Phase angles were derived from the Fourier transformed time series, unwrapped and subtracted to estimate the phase difference over time for each trial. Plateaus were defined as a constant phase angle (maximum average gradient <0.01 radians/ms) over the duration of one cycle of the stimulation frequency: ∑i=1ΔT|∇θi|n ⩽ 0.01rad/ms,(3) with ∇θi being the gradient, i.e., slope, of the phase angle between MEG and photodiode signal at a given sample i; n being the length of the cycle in ms, rounded up to the next integer, e.g., 17 ms for a flicker frequency of 60 Hz. This approach allowed to identify intermittent phase plateaus in each trial. In comparison, the PLV analysis described above quantifies the phase-similarity of the two signals over trials, and is therefore not feasible to capture brief episodes of synchrony between the MEG signal and the stimulation.

Statistical analysis

Statistical analysis was performed in RStudio version 1.2.1355 (RStudio Inc.; R version 3.6.1., The R Foundation for Statistical Computing).

Source analysis

MRI preprocessing

The raw T1 weighted images were converted from DICOM to NIFTI. The coordinate system of the participants' individual MRI was aligned to the anatomic landmarks using the head-surface obtained from the MRI and the scalp shapes digitized before the recordings. Realignment was done automatically using the iterative closest point (ICP) algorithm (Besl and McKay, 1992) implemented in the FieldTrip toolbox and corrected manually as necessary. The digitized head shape of one participant, for whom there was no anatomic image available, was aligned to a standardized template brain.

Linearly constrained minimum variance (LCMV) beamforming

The neuroanatomical origins of the visually induced gamma oscillations and the response induced by the photic drive condition were estimated using LCMV spatial filters (Veen et al., 1992), implemented in the FieldTrip toolbox (Oostenveld et al., 2011). The MEG forward model was calculated using single-shell head-models, estimated based on the aligned anatomic images, and an equally spaced 4 mm grid, warped into Montreal Neurologic Institute (MNI) space (Nolte, 2003; also see Oostenveld et al., 2011; Stenroos et al., 2012), yielding 37,163 dipoles inside the brain. The preprocessed data, epoched in 7- and 5-s trials for the respective conditions, were bandpass filtered at 50–92 Hz, by applying second order Butterworth two-pass high-pass and low-pass filters. To identify the peak locations of the endogenous gamma oscillations and flicker response, respectively, segments of 0.5 s of the baseline interval (0.75–0.25 s before stimulation) and the stimulation interval (0.75–1.25 s after flicker/grating onset) were extracted from the data in both conditions. The peak source of the flicker response to the flickering gratings was isolated based on the 2.75–3.25 interval, when the photic drive was superimposed on the gratings, contrasted with the 0.75–1.25 interval during which the gratings were presented. For each participant, a common covariance matrix for the 204 planar gradiometers was computed based on the extracted time series and used to estimate the spatial filter coefficients for each dipole location, whereby only the direction with the highest dipole moment was considered. Data in the baseline and stimulation intervals were projected to source space by multiplying each filter coefficient with the sensor time series. Fast Fourier transforms of the resulting time series, multiplied with a Hanning taper, were computed for each of the 37,163 virtual channels, separately for the baseline and stimulation intervals, and averaged over trials. Relative power change at the IGF and flicker frequencies was computed by applying Equation 1 to the Fourier-transformed baseline and stimulation intervals. The source-localized power change values at flicker frequencies up to 78 Hz were averaged to identify a common source for the oscillatory response to the photic drive.

Experimental design and statistical analyses

Using the experimental set up outlined above, this study aimed to explore resonance properties of the visual cortex, reflecting oscillatory dynamics in each participant. Furthermore, we asked whether responses to a visual flicker close to and at the IGF are enhanced when the flicker is superimposed on the moving grating stimuli. This would reflect a change in the oscillatory dynamics in presence of the endogenous gamma oscillations. In this context, we hypothesized that these oscillations would synchronize to the flicker. The 40 frequency × condition combinations were tested in a within-subject design. Resonance at individually preferred rhythms would be revealed by a relatively high response magnitude to stimulation at the preferred frequency in comparison to the surrounding frequencies (Herrmann, 2001; Schwab et al., 2006; Notbohm et al., 2016; H1). A general decrease in response to the flicker as a function of frequency would suggest an absence of such an amplification (H0). Entrainment of the ongoing gamma rhythm by the flicker response would result in the peak frequency of the gamma oscillator being synchronized to the stimulation frequency. This is reflected by a reduction in power at the IGF during the application of the flicker to the gratings, at frequencies different from the IGF, compared with the presentation of the gratings alone (H1). Statistical analyses were performed in R (R Core Team, 2020; version 3.6.3., using RStudio version 1.2.5033, RStudio Inc.). The statistical power of the individual tests was evaluated using Bayes factors, computed using the BayesFactor package in R (Morey and Rouder, 2018). As the identified IGF was found to be higher than the frequency inducing the strongest flicker response in the majority of participants, we quantified their relationship using a simple binomial test with an a priori defined α level of 0.01. The linearity of the flicker response power as a function of flicker frequency, i.e., evidence for the H0 as observed in the results reported below, was corroborated using linear regression models implemented in the R base package. Changes in the power at the IGF, with the onset of the flicker in the flicker&gratings condition, were examined using a repeated measures ANOVA on the factors time (pre and during flicker) and flicker frequency (above and below IGF), as implemented in package ez in R (Lawrence, 2016). Lastly, we compared the peak sources of the gamma oscillations and flicker responses, identified using LCMV beamforming, in both conditions using dependent sample t tests. As the direction of the distances was not known a priori, the α level was set to 0.025. To reduce the dimensionality of the comparisons, the obtained 3D coordinates were first projected along their first Principal Component (Herrmann et al., 2011). The p values of the three comparisons were corrected using the Benjamini–Hochberg procedure.

Results

The aim of the current study was to characterize entrainment and resonance properties in the visual cortex in absence and presence of gamma-band oscillations induced by visual gratings. To this end, we drove the visual cortex with a rapid flicker at frequencies ranging from 52 to 90 Hz, in steps of 2 Hz. The photic drive was applied either to a circular patch (the flicker condition; Fig. 1A,C) or to the light gray rings of a moving grating stimulus (the flicker&gratings condition; Fig. 1B,D). We hypothesized that a photic drive in the flicker&gratings condition would entrain the grating-induced oscillations. This would be observed as the endogenous gamma oscillation synchronizing with the flicker. Synchronization would be reflected by a constant phase angle between the neuromagnetic signal and the stimulation (“phase entrainment”), as well as a reduction in power at the IGF, indicating a change in the peak frequency of the gamma oscillator toward the flicker frequency (“frequency entrainment”; Pikovsky et al., 2003). Moreover, we expected the presence of the induced gamma oscillator to change the resonance properties (compared with the flicker condition), reflected by an amplification of responses to stimulation frequencies equal to the endogenous gamma rhythm. Response magnitudes in the flicker condition were expected to reveal resonance properties of the visual system in absence of gamma oscillations, demonstrating favorable stimulation frequencies to be used in future experiments applying RFT (Zhigalov et al., 2019; Drijvers et al., 2021).

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

The experimental paradigm. A, Trials in the flicker condition. A 1-s baseline interval with a central fixation cross was followed by a 2-s interval of the rapid flicker applied to a circular patch of size 2.62°. The average luminance in the flickering patch was equal to the surrounding gray color, making the photic drive almost unperceivable. The trials ended with 2 s of the fixation cross only. B, The trials in the flicker&gratings condition. The 1-s baseline interval was followed by 2 s of grating stimuli presented centrally on the screen, contracting inwards. Subsequently, the flicker was imposed onto the stimuli for 2 s. The trial ended with a 2-s presentation of the moving gratings without photic stimulation. C, Sinusoidal luminance change in one pixel induced by the photic drive at 52 Hz in the flicker condition. D, Luminance change in one pixel as a result of the flicker and the gratings moving concentrically with a velocity of four cycles per second. To maintain a similar mean luminance between conditions, photic modulation of the invisible patch in A ranged from 0% to 66% (mean RGB [84 84 84]), while the light gray rings of the grating, that is 50% of the stimulus' surface, were flickered between 33% and 99% (mean RGB [168 168 168] per ring).

Identifying IGFs

The frequency of the endogenous gamma rhythm is known to vary between participants (Hoogenboom et al., 2006, 2010; Muthukumaraswamy et al., 2010; Van Pelt et al., 2012). Therefore, each subject's IGF was identified first, based on the 0- to 2-s interval in the flicker&gratings condition during which the moving grating stimuli were presented without the visual flicker (Fig. 1C). The TFRs of power are depicted in Figure 2A,B for two representative participants. The center column shows the power averaged over time (0.25–1.75 s after the stimulus onset to avoid any event-related field confounds) demonstrating distinct peaks at 58 and 74 Hz for these participants. The topographies in the right column depict relative power change at the identified frequencies, focally in sensors over the occipital cortex. For each subject, the two to three combined planar gradiometers showing maximum relative power change in the gamma-band were selected for further analysis (SOIs) per visual inspection. These sensors strongly overlapped between participants. The data of participants with an IGF closer than 6 Hz to the lowest (52 Hz) drive, i.e., IGF <58 Hz, were not considered for further analyses.

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

Identification of the Individual Gamma Frequencies (IGFs) and Sensors-of-Interest (SOIs). A, B, The Time-Frequency representations (TFRs) of power, power spectra (averaged over 0.25–1.75 s) and topographic representations (combined planar gradiometers) of the IGF for two representative participants. The TFRs of power were calculated from the Fourier transforms using a 500-ms sliding window, resulting in spectral smoothing of ±3 Hz. The IGFs were identified from the spectral peak in 0.25- to 1.75-s interval of the TFRs. Identified IGFs are indicated by dashed lines. C, The grand average of the power analysis after aligning the individual TFRs and spectra to the IGF (N = 22).

Figure 2C depicts the averaged TFRs of power as well as the power spectrum for the remaining subjects (N = 22), aligned to each participant's IGF before averaging. The moving grating stimulus induced sustained oscillatory activity constrained to the IGF ±8 Hz, with an average relative power change of 80% in the 0.25- to 1.75-s interval compared with baseline. In short, the moving gratings produced robust gamma oscillations observable in the individual participants which reliably allowed us to identify the IGFs.

Photic drive induces responses up to 80 Hz

We next set out to quantify the rhythmic response to the flicker as a function of frequency in the flicker condition, in which stimulation was applied to an invisible patch. Figure 3 A,B, left panel, depicts the overlaid power spectra for the different stimulation frequencies in two representative participants (the same as in Fig. 2). The spectra were estimated by averaging the TFRs of power in the 0.25- to 1.75-s interval after flicker onset. Because of the overlap of the sensors detecting the gamma oscillations and photic drive response (compare Figs. 2, 3, right columns) the same SOI were used as in the flicker&gratings condition. Both individuals showed strong responses at the respective stimulation frequencies, with a maximum relative power change of 200% and 500% in subjects A and B, respectively. The identified IGFs (indicated by vertical dashed lines) were higher than the frequencies inducing the strongest flicker response in 20 out of 22 participants (exact Binomial Test against H0: p = 0.00012, probability of successes (IGF >flicker freq) 0.91, BayesFactorBF10 = 309.3). When averaged over all participants, the magnitude of the flicker response decreased systematically with frequency (Fig. 3C).

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

A, B, The response to the photic drive in the flicker condition and the corresponding topographies for two representative subjects. Spectra were estimated from the Time-Frequency representations (TFRs) of power averaged in the 0.25- to 1.75-s interval. Dashed vertical lines indicate the participants' Individual Gamma Frequencies (IGFs). The topographies (combined planar gradiometers) demonstrate a strong overlap with the ones in Figure 2. C, Grand average of the responses to the photic drive for each flicker frequency. On average, the magnitude of the flicker response decreases with increasing frequency, and is identifiable for stimulation below 80 Hz. D, Grand average flicker responses for frequencies from 52 to 90 Hz in steps of 4 Hz. The shaded areas, illustrating the standard deviation, indicate a substantial intersubject variability.

Figure 4A displays the power spectra in the flicker condition, estimated from the TFRs as explained above, averaged over all participants, as a function of stimulation frequency. These are equivalent to 3C. Diagonal values indicate the magnitude of the oscillatory responses (relative to baseline) at the stimulation frequencies, reaching values of up to 300% and decreasing monotonically with frequency. This confirms an upper limit for the stimulation of around 80 Hz. Off-diagonal values indicate oscillatory activity at frequencies different from the stimulation frequency. Figure 4B shows the same spectra after aligning to the IGFs, before averaging. Figure 4C,D display the spectra in the flicker&gratings condition (averaged in the 2.25- to 3.75-s interval), during which the photic drive was applied to the moving grating stimulus (see Fig. 1B). The induced gamma-band activity can be observed as the horizontal light red band at ∼60 Hz. When aligning the spectra to the IGF (Fig. 4D), we observe a decrease in the flicker response but no evidence for an amplification at or close to the IGF.

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

Average relative power change to the photic drive (y-axis) with respect to the driving frequencies (x-axis). A, The flicker condition. Note that the power changes mirror Figure 3C. Power decreases with increasing frequency, from a relative change of ∼3 at 52 Hz to ∼0.5 at 80 Hz. B, The flicker condition after the spectra were aligned to the Individual Gamma Frequency (IGF). C, The flicker&gratings condition. All spectra demonstrate both the flicker response and induced γ oscillation (observed as the light red horizontal band). Again, the amplitude of the rhythmic stimulation response appears to decrease with increasing frequency. D, The spectra for the flicker&gratings condition now aligned to the IGF. There is no indication that the rhythmic flicker captures the endogenous γ oscillations.

Magnitude of flicker response decreases as a function of frequency

The averaged TFRs of power in Figure 4 point to an approximately linear decrease in power of the flicker response with increasing frequency. Literature on neural resonance and entrainment, however, suggests the existence of a preferred rhythm at which oscillatory responses are amplified (Hutcheon and Yarom, 2000; Herrmann, 2001; Pikovsky et al., 2003; Notbohm et al., 2016; Gulbinaite et al., 2019). As argued in Pikovsky et al. (2003) phase-locking between the driving signal and the self-sustained oscillator is the most appropriate metric to investigate entrainment. Figure 5A,B depicts the PLV between the photodiode and the MEG signal at the SOI (planar gradiometers, not combined). This measure reveals a systematic decrease in phase-locking with increasing flicker frequency for both the flicker (orange) and flicker&gratings (blue) condition (A). The observed relationship is preserved when aligning the frequencies to the IGF (B, also see Table 1). Note the absence of increased phase-locking at the IGF. The magnitude of the flicker response, quantified by power change compared with baseline, as a function of frequency, is demonstrated in Figure 5C–F and depicts a similar relationship to the one observed for the PLV. The flicker condition (C, orange line) revealed a systematic decrease with frequency, whereas the flicker&gratings condition did show a peak at 56 Hz. However, this observed increase appeared to be caused by considerable variance between the power estimates of the individual participants (see Fig. 5E, each line graph depicts power estimates per individual participant). We again aligned the spectra to the IGF before computing the grand average (Fig. 5D). The absence of a peak at 0 Hz suggests no evidence for resonance at the IGF, confirming the peak at 56 Hz in C to be the result of intersubject variability. Indeed, simple linear regression models, fit individually to PLV and power as a function of frequency aligned to the IGF, separately for each condition, explain a considerable amount of the variance (see Table 1 and Fig. 5, dotted lines). We then identified the individual peak frequencies, eliciting the strongest response to the flicker in the flicker&gratings condition, and related those to the IGF, as seen in Figure 5F. As observed in the flicker condition, the frequency inducing the strongest response to the flicker was lower than the IGF in the majority of participants, i.e., 19 out of 22 (exact binomial test against H0: p = 0.0008, Bayes factor BF10 = 67.5).

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

Simple linear regression models: flicker response magnitude as a function of distance to IGF

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

Magnitude of the flicker response as a function of frequency in the flicker (orange) and flicker&gratings (blue) conditions. Shaded areas indicate the standard deviation. A, The phase-locking value (PLV) between the photodiode and the MEG signal over the SOIs as a function of driving frequency. B, The PLV between the photodiode and the MEG signals as a function of frequency after the spectra were aligned to the Individual Gamma Frequency (IGF). Again, the phase-locking decreases with increasing frequency (see Table 1 for a statistical quantification of the simple linear regression models). C, Relative power change with respect to baseline as a function of frequency. Generally, the power decreased with frequency, however, in the flicker&gratings condition, there is an apparent peak at ∼56 Hz. The shaded areas (standard deviation) indicate considerable variance between participants. D, Relative power change as a function of frequency after the individual spectra were aligned in frequency according to the IGF, demonstrating that responses to a photic drive at the IGF are not amplified. E, Relative power change as a function of frequency for each individual subject (N = 22), indicates that the peak at ∼56 Hz in C is driven by comparably high power in that frequency range in just a few individuals. F, Flicker frequency inducing highest power values versus IGF, demonstrating the IGF to be higher than the frequency inducing maximum power change in the majority of participants.

Gamma Oscillations and flicker response coexist

We initially hypothesized that entrainment of the gamma oscillations in the flicker&gratings condition would result in the photic drive capturing the oscillatory dynamics when the driving frequency was close to the IGF. Figure 6 depicts the TFRs of power relative to a 0.5-s baseline, for one representative subject (also shown in Figs. 2, 3A). The averaged trials for a photic drive at 52 Hz are shown in Figure 6A and separately for each flicker frequency in Figure 6B (figure created using function by Kumpulainen, 2020). The IGF (58 Hz for this subject) and the respective stimulation frequencies are indicated by dashed lines. The endogenous gamma oscillations, induced by the moving grating stimulus, are observed as the sustained power increase from 0 to 6 s, whereas the flicker response is demonstrated by a power increase at 2–4 s. The plots reveal that gamma oscillations persist at the IGF and coexist with the response to the photic drive, which is particularly apparent for stimulation at 52 Hz (Fig. 6A). Furthermore, the power increase at the flicker frequency does not appear to outlast termination of the drive at t = 4 s. In the subsequent step, we frequency-aligned the TFRs of power according to the IGF before averaging over participants. Again, the analyses were constrained to individuals with an IGF above 56 Hz (N = 22). The group averaged, aligned TFRs are shown in Figure 7 for frequencies ranging from IGF –6 Hz to IGF +16 Hz. The endogenous gamma oscillations are observed as the power increase extending from 0 to 6 s, and the flicker response as the power change in the 2- to 4-s interval marked by dashed lines, respectively. The photic stimulation induces a reliable response that decreases toward 12 Hz above the IGF. Despite the representation of the gamma oscillations being smoothed because of interindividual differences, the averaged aligned TFRs of power support the observations in the single subject data: both the gamma oscillations and flicker response coexist in the 2- to 4-s interval. Furthermore, there is no indication of the gamma power being reduced during the presentation of the flicker at frequencies close to, but different from, the IGF. In addition to the narrow-band gamma oscillations, the gratings elicited a rhythmic response at 4 Hz, i.e., the velocity of the concentric drift (data not shown). We did not find any evidence for an intermodulation between the frequency of the movement and the photic drive.

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

The Time-frequency representations (TFRs) of power for one representative subject, showing relative power change averaged over trials and SOIs in the flicker&gratings condition. A, Photic drive at 52 Hz. The moving grating stimuli were presented for 0–6 s, with the flicker superimposed from 2 to 4 s. Sustained γ-band activity is clearly observable throughout the presentation of the stimuli, with a power increase of 300% relative to baseline. Additionally, the rhythmic stimulation elicited a response at 52 Hz, which seems to coexist with the γ oscillations, indicating that the photic drive is unable to capture the dynamics of the γ oscillation. B, The plots for the frequencies from 52 to 90 Hz. Stimulation frequencies and Individual Gamma Frequency (IGF) (here 58 Hz) are indicated by horizontal dashed lines. The flicker induced responses up to 66 Hz in this participant. γ Oscillations persist in presence of flicker responses, suggesting that they coexist.

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

Grand average Time-frequency representations (TFRs) of power after aligning to the Individual Gamma Frequency (IGF) for each subject in the flicker&gratings condition. The stimulation frequencies (from −6 to 16 Hz relative to the IGF) are indicated by dashed horizontal lines. As suggested by the single subject TFRs in Figure 6, the endogenous γ oscillations and the flicker response seem to be coexistent. Thus, there is no obvious indication of the photic drive being able to capture the dynamics of the γ oscillations.

Frequency analyses with a longer time window confirm robustness of the reported results

To assess the robustness of our results, we repeated the frequency analyses in the flicker and flicker&gratings condition with a 2-s sliding time window. The longer window substantially increased the signal-to-noise ratio of the flicker response, to up to over 400% relative power change in the flicker condition and >600% in the flicker&gratings condition (data not shown). Besides that, the analyses replicated our reported main finding: a reduction in response magnitude (power) with increasing frequency, in both conditions, following the same trend as depicted in Figure 5C,D. The 2-s sliding time window did however not optimally capture the gamma power, which has a broader peak than the response to the photic drive. The 500-ms sliding window used in our reported analyses is therefore a good compromise, allowing both a reliable identification of a gamma peak frequency and a sufficiently high signal-to-noise ratio and frequency resolution of the flicker response (see Fig. 6A).

Oscillatory gamma dynamics cannot be captured by frequency entrainment

Synchronisation of neuronal oscillations by rhythmic stimulation could be conceptualized as the entrainment of a self-sustained oscillator by an external force (Notbohm et al., 2016; Helfrich et al., 2019). Frequency entrainment is reflected by a change in frequency of the ongoing oscillations toward the rhythm of the drive. Visual inspection of the TFRs of power in Figures 6, 7 do not indicate any modulation of the peak frequency of the gamma oscillations by the flicker response, suggesting that they do not synchronize. To quantify these observations, we investigated the power of the gamma oscillations before and during the photic drive (Fig. 8) in the flicker&gratings condition. A central assumption of oscillatory entrainment is the existence of a “synchronization region” in the frequency range around the endogenous frequency of the oscillator, the so-called Arnold tongue (Pikovsky et al., 2003). Driving frequencies falling inside this synchronization region, will be able to modulate the dynamics of the self-sustained oscillator (also see Hutt et al., 2018). With this in mind, the following analyses only included flicker frequencies in the vicinity of the IGF. For each participant, we considered the relative power change induced by the moving gratings in the 0.5- to 1.5-s interval (T1) before the flicker onset and in the 2.5- to 3.5-s interval (T2) in which both the moving gratings and the photic drive were present. We investigated this for stimulation frequencies below the IGF (averaged power for −6 and −4 Hz) and above (averaged power for +4 and +6 Hz). Assuming a symmetric Arnold tongue centered at the IGF, as shown for entrainment in the alpha-band (Notbohm et al., 2016), we expected a reduction in power at the IGF in interval T2 compared with interval T1 for both higher and lower driving frequencies, i.e., an effect of time, but not frequency. Figure 8 depicts power change at the IGF for the factors stimulation frequency (drive <IGF and drive >IGF) and time interval (T1 and T2), averaged over the SOIs for each subject. In accordance with the TFRs in Figure 7, there is no meaningful indication for gamma power being reduced during the T2 interval as compared with the T1 interval, affirming the coexistence of the two responses. A factorial repeated-measures ANOVA did not reveal any significant main effects of the factors time (T1 versus T2) and frequency (drive <IGF versus drive >IGF), but a significant interaction effect (F(1,21) = 5.09, p = 0.003, η2 = 0.003). These results were further investigated using a Bayesian repeated-measure ANOVA. The obtained Bayes factors (BF10) indicate that the variance in the data underlies the variability between participants, while the factor time (BF10 = 0.233) and both time and frequency (BF10 = 0.274) do not add any explanatory value. Evidence for the interaction effect time:frequency was found to be inconclusive (BF10 = 0.53), as was the main effect of frequency alone BF10 = 1.146). These results provide evidence against the expected reduction in gamma power during rhythmic photic stimulation at frequencies different from the IGF; suggesting that the flicker did not capture the oscillatory gamma dynamics.

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

Power change relative to baseline at Individual Gamma Frequency (IGF) in response to the moving grating stimuli before (T1; 0.5–1.5 s) and during application of the flicker (T2; 2.5–3.5 s), at frequencies below and above IGF (drive <IGF [−6, −4 Hz] and drive >IGF [+4, +6 Hz], respectively). Scatters demonstrate individual values, solid and dashed lines depict mean and standard deviation, respectively. The key finding is that power at T2 is not decreased compared with T1 for either of the frequency ranges, which is supported by a Bayesian repeated measures ANOVA (BF10 = 0.274).

Photic drive does not reliably modulate gamma phase

Synchronization of a self-sustained oscillator by an external force, cannot only be described by a change in frequency, but also “phase approximation” or “phase entrainment” (Pikovsky et al., 2003). This phenomenon is reflected by a constant phase angle between the two oscillators over extended intervals, so-called phase plateaus. These might occur when the frequency of the driver is close to the endogenous frequency of the oscillator, i.e., within its Arnold Tongue (Tass et al., 1998; Pikovsky et al., 2003; Notbohm et al., 2016). When approaching the edge of the synchronization region, episodes of constant phase angles are interrupted by so-called phase slips that emerge when the self-sustained oscillator briefly unlocks from the driving force and oscillates at its own frequency. These phase slips will be observed as steps between the phase plateaus. The phase plateau analysis was implemented to complement the PLV analysis shown in Figure 5. The PLV quantifies the average synchrony between photodiode and neuromagnetic signal over trials using a 500-ms sliding time window. We hypothesized that in the case of oscillatory entrainment, the gamma oscillator in the flicker&gratings condition would alternate between locking on to the photic drive for a few cycles and slipping back to its endogenous rhythm. Because of the short duration of the gamma cycle (∼17.2 ms for a 58-Hz IGF), this intermittency would be smeared out by the sliding window. As there was no endogenous gamma oscillator in the flicker condition, such an intermittency was not expected. To investigate phase entrainment of the gamma oscillations by the photic drive, we inspected the phase angle between the photodiode and one, individually selected, occipital gradiometer of interest per participant. The time series of the phase were estimated per trial, separately for the two sensors, using a sliding time-window Fourier transform approach (ΔT = 3 cycles = 3/fflickers; Hanning taper). Phase differences per trial were obtained by subtracting the unwrapped phase angle time series.

Phase angle between photodiode and MEG signal over time

Figure 9 illustrates the unwrapped phase angles between the MEG and photodiode signal during the photic drive at the IGF (here 58 Hz), in the flicker (A) and flicker&gratings condition (B), respectively, for the same representative participant shown in Figures 2A, 3A, 6. The colored line graphs depict individual trials. In both conditions, the MEG signal drifts apart from the photic drive, toward a maximum difference of 60 radians, i.e., a phase difference of ∼9.5 cycles, by the end of the trial (Fig. 9A and B, top panel). Interestingly, the direction of the phase angle appears to change during some of the trials, suggesting spectral instability of the gamma oscillations. Furthermore, the graphs demonstrate a substantial intertrial variability. This diffusion between trials, quantified for each participant as the standard deviation over trials at the end of the photic stimulation (t = 2 in flicker and t = 4 in flicker&gratings condition), converted from radian to ms, is juxtapositioned in Figure 9C for the two conditions. It can be readily seen that the phase angles between the stimulation and MEG signal fan out highly similarly in absence and presence of the endogenous gamma oscillations.

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

A, B, Phase angle between photodiode and the MEG signal (one gradiometer of interest) at the Individual Gamma Frequency (IGF), for one representative participant; colored lines depict individual trials. A, Phase angle θ in the flicker condition over duration of the flicker presentation (upper panel) and the first 250 ms (lower panel). The MEG signal drifts apart from the stimulation and can reach a maximum accumulated phase difference of 60 radians, i.e., 9.54 cycles, at the end of the stimulation and up to 15 radians, i.e., 2.39 cycles, in 250 ms. B, The increase in phase difference over the time of the stimulation for the flicker&gratings condition (upper panel) and in the first 250 ms (lower panel). The diffusion of the phase difference across trials is similar to the flicker condition. Moreover, there is no clear difference in the number and length of phase plateaus between conditions, implying that the presence of the γ oscillations does not facilitate entrainment at the IGF. C, Fanning out across trials as a function of frequency aligned to IGF. Trials diffuse to a highly similar extent in both conditions and across frequencies. D, Number of plateaus per trial as a function of frequency. While the flicker&gratings condition exhibits more plateaus for all flicker frequencies, there is no indication that stimulation at the IGF results in comparably strong synchronization.

Phase plateaus

Visual inspection of the first 0.25 s of the phase angle times series, depicted in Figure 9A,B, lower panel, does not suggest a relatively high number of phase plateaus in the flicker&gratings compared with the flicker condition, that would have been expected if the photic drive was able to entrain the endogenous gamma oscillator. Importantly, the graphs demonstrate the phase angles to reach values of over 2π, i.e., more than one cycle, within the duration of the first gamma cycle (17.2 ms), suggesting that even stimulation at the endogenous frequency of the oscillator cannot capture the gamma dynamics. To verify these observations for the entire sample, plateaus during stimulation at the IGF were identified based on the mean absolute gradient (≤0.01 radians/ms, see Eq. 3) over the duration of one cycle of stimulation, i.e., 18 consecutive samples for a flicker frequency of 58 Hz. Figure 9D shows the average number of plateaus per trial as a function of flicker frequency aligned to IGF, averaged over participants. The shaded areas indicate the standard deviation. While the flicker&gratings condition exhibits more phase plateaus than flicker for all stimulation frequencies, the number of plateaus decreases similarly in both conditions with increasing frequency. Importantly, stimulation at the IGF did not result in the highest number of plateaus in either condition. These results are in line with the reported frequency analyses: responses to the photic drive in flicker&gratings show strong similarity to the flicker condition despite the presence of the gamma oscillator. The results affirm the observations presented in Figure 5A,B.

The sources of the gamma oscillations and the flicker responses peak at different locations

The coexistence of the endogenous gamma oscillations and flicker response suggest that these two signals are generated by different neuronal populations; possibly in different regions. To test this assumption we localized the respective sources using LCMV spatial filters (Veen et al., 1992). The covariance matrix for the spatial filters was estimated based on the −0.75- to −0.25-s baseline in both conditions, the 0.75- to 1.25-s interval with the moving gratings in flicker&gratings and the invisible flicker in the flicker condition, as well as the 2.75–3.25 interval in the flicker&gratings condition in which the flicker was applied to the grating stimulus. Note that for each participant, one common filter was used for source estimation in both conditions. Power values at the IGF and flicker frequencies, averaged up to 78 Hz, respectively, for the flicker&gratings and flicker condition, were estimated based on the Fourier transform. To extract power at the IGF and flicker frequencies, power change was computed relative to the baseline interval at each of the 37,163 grid points using Equation 1. To isolate the flicker response on the flicker&gratings condition, the flicker&gratings interval was contrasted to the moving grating interval. Figure 10 illustrates the grand average of the source localization for the gamma oscillations (Fig. 10A), the invisible flicker response (Fig. 10B), and the response to the flickering gratings (Fig. 10C). Consistent with previous work, the responses originate from mid-occipital regions (Hoogenboom et al., 2006; Zhigalov et al., 2019). It is worth noting that the sources of the gamma oscillations and response to the invisible flicker are relatively focal, while the activity induced by the flickering gratings extends more broadly over visual cortex. Using the MNI to Talaraich mapping online tool by Biomag Suite Web (MNI2TAL Tool; see Lacadie et al., 2007, 2008), the peak of the gamma oscillations was located in the ventral part of the secondary visual cortex (V2; Brodmann area 18; MNI coordinates = [–6 mm –100 mm –8 mm], grand average). The peak sources of the flicker responses in both conditions were found in the calcarine fissure, at a 2 mm distance to the border of the primary visual cortex (V1) and V2 (in dorsal direction); suggesting that they are generated by neighboring, coherent sources in both hemispheres in and close to V1 (Belardinelli et al., 2012; MNI coordinates: flicker [6 mm –96 mm 12 mm]; flicker&gratings [6 mm –100 mm 0 mm]). To compare the peak locations between the sources in a lower dimensional space, the identified 3D coordinates were projected along their first principal component (Herrmann et al., 2011). Dependent sample t tests revealed a significant difference in location between the peak sources of the IGF and the invisible flicker responses, t(21) = −3.091, p = 0.017, Cohen's d = –0.845, 95% CI [–1.5–0.2], B10 = 8.2, as well as to the flickering gratings relative to gratings, t(21) = –2.633, p = 0.023, Cohen's d = –0.495, 95% CI [–0.89–0.09], B10 = 3.45; with the Bayes factors B10 revealing moderate evidence for the H1 (Quintana and Williams, 2018). There was no significant difference in location between the sources of the flicker responses in both conditions, t(21) = 0.732, p = 0.472, B10 = 0.28, with the Bayes factor providing moderate evidence for the H0. Note that all t values were Benjamini–Hochberg-corrected for multiple comparisons. In light of the coexistence of the two responses observed in Figures 6, 7, these results support the notion that gamma oscillations and flicker responses are generated by different neuronal populations.

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

Source estimates using the LCMV beamformer approach mapped on a standardized MNI brain. A, Source estimation of the visually induced γ oscillations (IGF; power change relative to baseline), with the peak of the source identified at MNI coordinates [–6 mm –100 mm –8 mm]. B, Source estimation of the flicker response (relative to baseline), with the average peak source at [6 mm –96 mm 12 mm] (in calcarine fissure). C, Source estimation of the flicker response in the flicker&gratings condition (relative to the gratings interval), with the average peak source at [6 mm –100 mm 0 mm] (in calcarine fissure). D, Coordinates of the identified peak sources for all participants (small scatters) and grand average (large scatters) for the IGF, and the flicker responses in the flicker and flicker&gratings condition (green, orange, and blue, respectively). The peak sources of the flicker responses are adjacent, while the γ sources tend to peak at inferior locations.

Discussion

In this MEG study, we explored resonance and entrainment in the human visual system in response to a rapid photic drive >50 Hz. Strong, sustained gamma oscillations were induced using moving grating stimuli (Hoogenboom et al., 2006, 2010; Muthukumaraswamy and Singh, 2013; Van Pelt and Fries, 2013) and used to identify each participant's gamma frequency. The superposition of the flicker and the gratings allowed us to investigate whether the flicker is able to entrain endogenous gamma oscillations. The photic drive induced responses for frequencies up to ∼80 Hz, both in presence and absence of grating-induced endogenous gamma oscillations. To our surprise, we did not find evidence for resonance, i.e., an amplification of an individually preferred frequency in the range of the rhythmic stimulation, in either condition, despite the IGF being above 50 Hz in all participants. Moreover, there was no indication that the endogenous gamma oscillations synchronized with the rhythmic stimulation, i.e., no evidence for entrainment. Despite their differences, the flicker responses in the two conditions show strong similarities in the phase and frequency measures, supporting the notion that the flicker response coexists with the grating-induced oscillations. In accordance with these results, source estimation using LCMV spatial filters (Veen et al., 1992), suggests that the neuronal sources of the flicker responses in both conditions and the endogenous gamma oscillations peak at different locations in visual cortex.

Flicker responses do not entrain the gamma oscillator

While the sources of the gamma oscillations and the response to the (nearly) invisible flicker did overlap in occipital cortex, their peak coordinates were found to be significantly different. Relative power change at the IGF peaked at sources inferior to the flicker responses in both conditions, and was located in the left V2 using the MNI2TAL online tool (see Lacadie et al., 2007, 2008). The flicker peak sources were located in the calcarine fissure, in close proximity to V1. These results are in line with the coexistence of the endogenous oscillations indicated by the time-frequency analyses and might be the result of the filter properties of synaptic transmission as the flicker response propagates in the visual system (see Kuffler, 1953; Hawken et al., 1996; Carandini et al., 1997; Ringach, 2004; Cormack, 2005; Shadlen and Movshon, 1999). Low-pass filtering at the transition from the thalamus to V1 (Connelly et al., 2015) might attenuate the photic drive at frequencies above 80 Hz, leading to an absence of measurable responses in this range. Low-pass filter properties in V1 in projections from granular layers (L4a, 4cα, and 4cβ) to supragranular (L2/3, 4b) and infragranular layers (L5,6; Hawken et al., 1996; Douglas and Martin, 2004; Fröhlich, 2016) might have prevented the flicker response to converge to the neuronal circuits generating the endogenous gamma rhythms. This idea is supported by intracranial recordings in macaques showing the strongest gamma synchronization in response to drifting grating stimuli in V1 in supragranular layers (L2/3 and 4B; Xing et al., 2012), whereas steady-state responses to a 60-Hz photic flicker have been localized in granular layer 4cα (Williams et al., 2004). While plausible, these interpretations are conjectural based on the present data. Recent findings by Drijvers et al. (2021), providing evidence for nonlinear integration of visual and auditory frequency tagging signals in frontal and temporal regions, challenge the notion that the flicker response might not propagate beyond V1. Pairing the current paradigm with intracranial recordings in non-human primates would allow to test the filtering properties without the limitations imposed by the inverse problem in the source localization of neuromagnetic signals (Baillet, 2013).

Flicker responses might not be wired to inhibitory interneurons orchestrating the endogenous gamma rhythm

Computational models, as the one demonstrated by Tiesinga (2012; also see Lee and Jones, 2013), would be suitable to investigate whether the grating-induced gamma oscillations and flicker response are likely to be generated by neuronal circuits whose wiring is not conducive to entrainment. As the properties of neuronal gamma oscillations have been repeatedly shown to depend on rhythmic inhibition imposed by inhibitory interneurons (Wilson and Cowan, 1972; Bartos et al., 2007; Buzsáki and Wang, 2012; Lozano-Soldevilla et al., 2014; Kujala et al., 2015), entrainment should only be achieved when the flicker response is able to modulate their activity. Indeed, Cardin et al. (2009) show resonance in the gamma range to optogenetic stimulation of fast-spiking interneurons, but not to stimulation of pyramidal cells (also see Tiesinga, 2012). We therefore suggest that the photic stimulation applied in our study drives the pyramidal cells in early visual cortex. As in the optogenetic study by Cardin et al. (2009), this drive is not sufficiently strong to entrain the GABAergic interneurons. This interpretation is contrasted to the findings of Adaikkan et al. (2019) who demonstrate that a non-invasive 40-Hz flicker evokes neuronal processes counteracting neuro-degeneration (Singer et al., 2018; Adaikkan et al., 2019). However, it should be noted that the authors understand entrainment as the neural response to rhythmic stimulation, rather than a synchronization of ongoing oscillations to an external drive (Adaikkan and Tsai, 2020). While our findings do not question the authors' compelling evidence that fast photic stimulation impacts neurocircuits and glia, the current study shows that it is not trivial to attribute these effects to entrainment of endogenous gamma oscillations.

Coexistence of flicker responses and oscillations versus oscillatory entrainment

The current study was inspired by studies reporting that a visual flicker in the alpha-band can capture the oscillatory dynamics of the visual system: resonance at distinct frequencies (Herrmann, 2001; Schwab et al., 2006; Gulbinaite et al., 2019; see Rager and Singer, 1998, for flicker responses in cat visual cortex), amplitude and phase effects outlasting the stimulation interval (Spaak et al., 2014; Otero et al., 2020) and an “Arnold Tongue” relationship between stimulation intensity, distance to the individual α frequency and flicker-response-synchrony (Notbohm et al., 2016). Unlike the works listed above, we did not find any indication for a synchronization or resonance of endogenous oscillations in the gamma-band to the visual stimulation. Recent studies applying photic stimulation in the alpha-band, have pointed to a coexistence of endogenous alpha oscillations and flicker responses, similar to the one we report here for the gamma-band. While retinotopic α modulation has been associated with suppression of unattended stimuli, allocating attention to a stimulus flickering in the alpha-band results in enhanced, phase-locked activity (Keitel et al., 2019; Gundlach et al., 2020; also see Antonov et al., 2020; Friedl and Keil, 2020 for stimulation at frequencies adjacent to the alpha-band). While the presented study does not allow nor aim to make generalized claims in favor or against neuronal entrainment, it is worth noting that the ability of rhythmic sensory stimulation to entrain endogenous oscillations is still a matter of debate.

Limitations and generalizability

Interpretation of the different locations of the peak sources

The results of the LCMV beamforming are in line with the notion that gamma oscillations and flicker response are generated by sources at different locations. Yet, because of the ill-posed inverse problem (Baillet, 2013) and the merging of coherent sources when using the LCMV approach (Belardinelli et al., 2012) these source estimates should be interpreted with caution. Figure 10 illustrates that the sources of the flicker response in the flicker&gratings condition extended more broadly over visual cortex than the sources of the gamma oscillations and invisible flicker response, which might be the result of the flickering rings stimulating different receptive fields (Gur and Nodderly, 1997). While our results suggest a coexistence of the gamma oscillations and flicker response, we do not exclude that they interact. These limitations do not seriously challenge our interpretation that the neuronal populations generating the flicker response do not entrain the activity of the neurons engaging in the endogenous gamma rhythm. First, it is reasonable to assume that the peak sources reflect the flicker response, which tends to be stronger than the endogenous gamma oscillations (see Figs. 6, 7). Second, the significant difference between the peak locations of the gamma oscillations and flicker response in the flicker&gratings condition provides circumstantial evidence for the notion that the two responses emerge from different neuronal populations, despite being elicited by the same stimulus; albeit there is also an overlap between the sources. Intracranial recordings in nonhuman primates or humans would be useful to substantiate this interpretation.

Strong flicker responses despite limited stimulation strength

The number of conditions that have been tested in this paradigm, i.e., 40 frequency × condition combinations, imposed limitations on the maximum number of trials per condition (N = 15) and the duration of the stimulation (2 s). Stimulation strength was limited to a contrast of 66% peak to trough, ensuring equal luminance across conditions. Because of these limitations, one might be concerned that the absence of oscillatory entrainment was caused by the limited magnitude of the photic drive. However, we found the flicker to induce strong responses of up to 400% in the flicker&gratings condition and over 200% in the flicker condition (see Figs. 4, 5). In light of these response magnitudes, we argue that the absence of evidence for entrainment cannot be explained by the photic drive being too weak.

Generalizability of the current findings to gamma oscillations associated with visual perception

The use of drifting gratings is a standard approach to induce strong narrow-band gamma oscillations in humans (Hoogenboom et al., 2006, 2010; Van Pelt et al., 2012; Muthukumaraswamy and Singh, 2013; Van Pelt and Fries, 2013; Michalareas et al., 2016) and nonhuman primates (Womelsdorf et al., 2006; Buffalo et al., 2011; Bosman et al., 2012). One might argue that the conclusions presented here only apply to these stimuli and that entrainment could have been achieved using more complex stimuli such as natural images or faces. We find this very unlikely for the following reasons: natural stimuli have been argued to induce gamma-band responses that are characterized by broadband activity (Ray and Maunsell, 2010; Hermes et al., 2015a,b; but also see Brunet et al., 2014; Bartoli et al., 2019; Brunet and Fries, 2019). This is likely explained by the fact that gamma power and frequency depend on stimulus properties such as contrast, size and orientation (Schadow et al., 2007; Ray and Maunsell, 2010; Jia et al., 2013; Muthukumaraswamy and Singh, 2013). As these factors vary greatly within a natural image, the net result of the oscillatory activity in the gamma-band is a broadband response. Moving gratings have been shown to induce stronger gamma oscillations than their stationary counterparts (Muthukumaraswamy and Singh, 2013; Perry et al., 2013) and were therefore chosen for the current paradigm. We expected the flicker responses to be substantially stronger than the grating-induced gamma oscillations, which is confirmed by Figures 6, 7. Had we relied on stationary gratings, the photic drive might have overshadowed weaker gamma-band activity. Moreover, the frequencies of the endogenous gamma rhythms have been found to be higher for moving than for stationary gratings (Muthukumaraswamy and Singh, 2013; Perry et al., 2013). As our study aimed to investigate entrainment by a flicker with minimal visibility, the IGFs had to be relatively high to be in the range of feasible stimulation frequencies. While the gratings' concentric drift in our study did induce a rhythmic response at 4 Hz, there was no evidence for an intermodulation with the flicker frequencies, nor an indication that the flicker&gratings condition was lacking spectral precision. In line with our findings, recent work by Lobo et al. (2021, biorxiv) using a superposition of a 60 Hz flicker and static gratings, has demonstrated the coexistence of the grating-induced broadband gamma activity and the flicker response. Furthermore, the authors report that behavior only correlated with the broadband gamma activity in absence of the flicker. When the photic drive was applied to the static grating, the flicker response correlated with reaction time, but behaviour did no longer relate to the broadband gamma activity. Another concern might be that grating stimuli do not engage downstream regions to the same extent as complex stimuli; as such they might be generated in specialized neuronal circuits. However, a number of studies in both human and non-human primates have demonstrated that attended as well as unattended gratings induce gamma oscillations that propagate to downstream areas along the ventral (V4 and inferotemporal cortex) and dorsal stream (area V5 and V7; Buffalo et al., 2011; Bosman et al., 2012; Bastos et al., 2015; Michalareas et al., 2016). For the reasons outlined above, we argue that moving grating stimuli created the optimal conditions to investigate gamma-band entrainment, as these induced strong, sustained, narrow-band gamma oscillations reflecting individual oscillatory dynamics (also see Hoogenboom et al., 2006; Van Pelt and Fries, 2013).

Conclusion

In conclusion, our results suggest that rapid photic stimulation does not entrain endogenous gamma oscillations and can therefore not be used as a tool to probe the causal role of gamma oscillations in cognition and perception. However, the approach can be applied in Rapid Frequency Tagging to track neuronal responses for instance, to investigate covert spatial attention (Zhigalov et al., 2019), multisensory integration (Drijvers et al., 2021), and parafoveal reading (Pan et al., 2020).

Footnotes

  • This work was supported by the James S. McDonnell Foundation Understanding Human Cognition Collaborative Award 220020448, the Wellcome Trust Investigator Award in Science 207550, the Biotechnology and Biological Sciences Research Council Grant BB/R018723/1, and a Royal Society Wolfson Research Merit Award (O.J.). We thank Prof. Veikko Jousmaki for providing the light-to-voltage converter and Jonathan L. Winter, Nina Salman, Ludwig Barbaro, and Roya Jalali for providing help with the MEG recordings and MRI scans. We also thank Dr. Simon Hanslmayr, Dr. Geoffrey Brookshire, and Dr. Florian Kasten for feedback on the project and manuscript.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Katharina Duecker at k.duecker{at}bham.ac.uk or Ole Jensen at o.jensen{at}bham.ac.uk

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    1. Adaikkan C,
    2. Tsai LH
    (2020) Gamma entrainment: impact on neurocircuits, glia, and therapeutic opportunities. Trends Neurosci 43:24–41. doi:10.1016/j.tins.2019.11.001 pmid:31836315
    OpenUrlCrossRefPubMed
  2. ↵
    1. Adaikkan C,
    2. Middleton SJ,
    3. Marco A,
    4. Pao PC,
    5. Mathys H,
    6. Kim DNW,
    7. Gao F,
    8. Young JZ,
    9. Suk HJ,
    10. Boyden ES,
    11. McHugh TJ,
    12. Tsai LH
    (2019) Gamma entrainment binds higher-order brain regions and offers neuroprotection. Neuron 102:929–943.e8. doi:10.1016/j.neuron.2019.04.011 pmid:31076275
    OpenUrlCrossRefPubMed
  3. ↵
    1. Antonov PA,
    2. Chakravarthi R,
    3. Andersen SK
    (2020) Too little, too late, and in the wrong place: alpha band activity does not reflect an active mechanism of selective attention. Neuroimage 219:117006. doi:10.1016/j.neuroimage.2020.117006 pmid:32485307
    OpenUrlCrossRefPubMed
  4. ↵
    1. Başar-Eroglu C,
    2. Strüber D,
    3. Schürmann M,
    4. Stadler M,
    5. Başar E
    (1996) Gamma-band responses in the brain: a short review of psychophysiological correlates and functional significance. Int J Psychophysiol 24:101–112. doi:10.1016/S0167-8760(96)00051-7
    OpenUrlCrossRefPubMed
  5. ↵
    1. Baillet S
    (2013) Forward and inverse problems of MEG/EEG. In: Encyclopedia of computational neuroscience, pp 1–8. New York: Springer.
  6. ↵
    1. Bartoli E,
    2. Bosking W,
    3. Chen Y,
    4. Li Y,
    5. Sheth SA,
    6. Beauchamp MS,
    7. Yoshor D,
    8. Foster BL
    (2019) Functionally distinct gamma range activity revealed by stimulus tuning in human visual cortex. Curr Biol 29:3345–3358. doi:10.1016/j.cub.2019.08.004 pmid:31588003
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bartos M,
    2. Vida I,
    3. Jonas P
    (2007) Synaptic mechanisms of synchronized gamma oscillations in inhibitory interneuron networks. Nat Rev Neurosci 8:45–56. doi:10.1038/nrn2044 pmid:17180162
    OpenUrlCrossRefPubMed
  8. ↵
    1. Bastos AM,
    2. Schoffelen JM
    (2016) A tutorial review of functional connectivity analysis methods and their interpretational pitfalls. Front Syst Neurosci 9:175. doi:10.3389/fnsys.2015.00175 pmid:26778976
    OpenUrlCrossRefPubMed
  9. ↵
    1. Bastos AM,
    2. Vezoli J,
    3. Bosman CA,
    4. Schoffelen JM,
    5. Oostenveld R,
    6. Dowdall JR,
    7. De Weerd P,
    8. Kennedy H,
    9. Fries P
    (2015) Visual areas exert feedforward and feedback influences through distinct frequency channels. Neuron 85:390–401. doi:10.1016/j.neuron.2014.12.018 pmid:25556836
    OpenUrlCrossRefPubMed
  10. ↵
    1. Belardinelli P,
    2. Ortiz E,
    3. Braun C
    (2012) Source activity correlation effects on lcmv beamformers in a realistic measurement environment. Comput Math Methods Med 2012:190513. doi:10.1155/2012/190513 pmid:22611439
    OpenUrlCrossRefPubMed
  11. ↵
    1. Besl PJ,
    2. McKay ND
    (1992) A method for registration of 3-D shapes. IEEE Trans Pattern Anal Mach Intell 14:239–256. doi:10.1109/34.121791
    OpenUrlCrossRef
  12. ↵
    1. Bosman CA,
    2. Schoffelen JM,
    3. Brunet N,
    4. Oostenveld R,
    5. Bastos AM,
    6. Womelsdorf T,
    7. Rubehn B,
    8. Stieglitz T,
    9. De Weerd P,
    10. Fries P
    (2012) Attentional stimulus selection through selective synchronization between monkey visual areas. Neuron 75:875–888. doi:10.1016/j.neuron.2012.06.037 pmid:22958827
    OpenUrlCrossRefPubMed
  13. ↵
    1. Brainard DH
    (1997) The psychophysics toolbox. Spat Vis 10:433–436. pmid:9176952
    OpenUrlCrossRefPubMed
  14. ↵
    1. Bressler SL
    (1990) The gamma wave: a cortical information carrier? Trends Neurosci 13:161–162. doi:10.1016/0166-2236(90)90039-d pmid:1693231
    OpenUrlCrossRefPubMed
  15. ↵
    1. Brosch M,
    2. Budinger E,
    3. Scheich H
    (2002) Stimulus-related gamma oscillations in primate auditory cortex. J Neurophysiol 87:2715–2725. doi:10.1152/jn.2002.87.6.2715 pmid:12037173
    OpenUrlCrossRefPubMed
  16. ↵
    1. Brunet NM,
    2. Fries P
    (2019) Human visual cortical gamma reflects natural image structure. Neuroimage 200:635–643. doi:10.1016/j.neuroimage.2019.06.051 pmid:31247299
    OpenUrlCrossRefPubMed
  17. ↵
    1. Brunet N,
    2. Vinck M,
    3. Bosman CA,
    4. Singer W,
    5. Fries P
    (2014) Gamma or no gamma, that is the question. Trends Cogn Sci 18:507–509. doi:10.1016/j.tics.2014.08.006 pmid:25199855
    OpenUrlCrossRefPubMed
  18. ↵
    1. Buffalo EA,
    2. Fries P,
    3. Landman R,
    4. Buschman TJ,
    5. Desimone R
    (2011) Laminar differences in gamma and alpha coherence in the ventral stream. Proc Natl Acad Sci USA 108:11262–11267. doi:10.1073/pnas.1011284108 pmid:21690410
    OpenUrlAbstract/FREE Full Text
  19. ↵
    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
  20. ↵
    1. Carandini M,
    2. Heeger DJ,
    3. Movshon JA
    (1997) Linearity and normalization in simple cells of the macaque primary visual cortex. J Neurosci 17:8621–8644. doi:10.1523/JNEUROSCI.17-21-08621.1997
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Cardin JA,
    2. Carlén M,
    3. Meletis K,
    4. Knoblich U,
    5. Zhang F,
    6. Deisseroth K,
    7. Tsai LH,
    8. Moore CI
    (2009) Driving fast-spiking cells induces gamma rhythm and controls sensory responses. Nature 459:663–667. doi:10.1038/nature08002 pmid:19396156
    OpenUrlCrossRefPubMed
  22. ↵
    1. Connelly WM,
    2. Laing M,
    3. Errington AC,
    4. Crunelli V
    (2015) The thalamus as a low pass filter: filtering at the cellular level does not equate with filtering at the network level. Front Neural Circuits 9:89. doi:10.3389/fncir.2015.00089 pmid:26834570
    OpenUrlCrossRefPubMed
  23. ↵
    1. Cormack LK
    (2005) Computational models of early human vision. In: Handbook of image and video processing, pp 325–345. New York: Elsevier.
  24. ↵
    1. Douglas RJ,
    2. Martin KA
    (2004) Neuronal circuits of the neocortex. Annu Rev Neurosci 27:419–451. doi:10.1146/annurev.neuro.27.070203.144152 pmid:15217339
    OpenUrlCrossRefPubMed
  25. ↵
    1. Drijvers L,
    2. Jensen O,
    3. Spaak E
    (2021) Rapid invisible frequency tagging reveals nonlinear integration of auditory and visual information. Hum Brain Mapp 42:1138–1152.
    OpenUrl
  26. ↵
    1. Eckhorn R,
    2. Bauer R,
    3. Jordan W,
    4. Brosch M,
    5. Kruse W,
    6. Munk M,
    7. Reitboeck HJ
    (1988) Coherent oscillations: a mechanism of feature linking in the visual cortex? Multiple electrode and correlation analyses in the cat. Biol Cybern 60:121–130. doi:10.1007/BF00202899 pmid:3228555
    OpenUrlCrossRefPubMed
  27. ↵
    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
  28. ↵
    1. Gray CM,
    2. Engel AK,
    3. König P,
    4. Singer W
    (1992) Synchronization of oscillatory neuronal responses in cat striate cortex: temporal properties. Vis Neurosci 8:337–347. doi:10.1017/s0952523800005071 pmid:1562569
    OpenUrlCrossRefPubMed
  29. ↵
    1. Friedl WM,
    2. Keil A
    (2020) Effects of experience on spatial frequency tuning in the visual system: behavioral, visuocortical, and alpha-band responses. J Cogn Neurosci 32:1153–1169. doi:10.1162/jocn_a_01524 pmid:31933434
    OpenUrlCrossRefPubMed
  30. ↵
    1. Fries P,
    2. Nikolić D,
    3. Singer W
    (2007) The gamma cycle. Trends Neurosci 30:309–316. doi:10.1016/j.tins.2007.05.005 pmid:17555828
    OpenUrlCrossRefPubMed
  31. ↵
    1. Fröhlich F
    (2016) Network neuroscience. San Diego: Academic Press.
  32. ↵
    1. Gray CM,
    2. Singer W
    (1989) Stimulus-specific neuronal oscillations in orientation columns of cat visual cortex. Proc Natl Acad Sci USA 86:1698–1702. doi:10.1073/pnas.86.5.1698 pmid:2922407
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Grützner C,
    2. Wibral M,
    3. Sun L,
    4. Rivolta D,
    5. Singer W,
    6. Maurer K,
    7. Uhlhaas P
    (2013) Deficits in high-(> 60 hz) gamma-band oscillations during visual processing in schizophrenia. Front Hum Neurosci 7:88. doi:10.3389/fnhum.2013.00088 pmid:23532620
    OpenUrlCrossRefPubMed
  34. ↵
    1. Gulbinaite R,
    2. Roozendaal DH,
    3. VanRullen R
    (2019) Attention differentially modulates the amplitude of resonance frequencies in the visual cortex. Neuroimage 203:116146. doi:10.1016/j.neuroimage.2019.116146 pmid:31493535
    OpenUrlCrossRefPubMed
  35. ↵
    1. Gundlach C,
    2. Moratti S,
    3. Forschack N,
    4. Müller M
    (2020) Spatial attentional selection modulates early visual stimulus processing independently of visual alpha modulations. Cereb Cortex 30:3686–3703. doi:10.1093/cercor/bhz335 pmid:31907512
    OpenUrlCrossRefPubMed
  36. ↵
    1. Gur M,
    2. Nodderly DM
    (1997) Visual receptive fields of neurons in primary visual cortex (v1) move in space with the eye movements of fixation. Vision Res 37:257–265. doi:10.1016/S0042-6989(96)00182-4
    OpenUrlCrossRefPubMed
  37. ↵
    1. Hawken MJ,
    2. Shapley RM,
    3. Grosof DH
    (1996) Temporal-frequency selectivity in monkey visual cortex. Vis Neurosci 13:477–492. doi:10.1017/s0952523800008154 pmid:8782375
    OpenUrlCrossRefPubMed
  38. ↵
    1. Helfrich RF,
    2. Breska A,
    3. Knight RT
    (2019) Neural entrainment and network resonance in support of top-down guided attention. Curr Opin Psychol 29:82–89. doi:10.1016/j.copsyc.2018.12.016 pmid:30690228
    OpenUrlCrossRefPubMed
  39. ↵
    1. Hermes D,
    2. Miller KJ,
    3. Wandell BA,
    4. Winawer J
    (2015a) Gamma oscillations in visual cortex: the stimulus matters. Trends Cogn Sci 19:57–58. doi:10.1016/j.tics.2014.12.009 pmid:25575448
    OpenUrlCrossRefPubMed
  40. ↵
    1. Hermes D,
    2. Miller K,
    3. Wandell B,
    4. Winawer J
    (2015b) Stimulus dependence of gamma oscillations in human visual cortex. Cereb Cortex 25:2951–2959. doi:10.1093/cercor/bhu091 pmid:24855114
    OpenUrlCrossRefPubMed
  41. ↵
    1. Herrmann B,
    2. Maess B,
    3. Hahne A,
    4. Schröger E,
    5. Friederici AD
    (2011) Syntactic and auditory spatial processing in the human temporal cortex: an MEG study. Neuroimage 57:624–633. doi:10.1016/j.neuroimage.2011.04.034 pmid:21554964
    OpenUrlCrossRefPubMed
  42. ↵
    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
  43. ↵
    1. Herrmann CS,
    2. Mecklinger A
    (2001) Gamma activity in human EEG is related to highspeed memory comparisons during object selective attention. Vis Cogn 8:593–608. doi:10.1080/13506280143000142
    OpenUrlCrossRef
  44. ↵
    1. Herrmann CS,
    2. Demiralp T
    (2005) Human EEG gamma oscillations in neuropsychiatric disorders. Clin Neurophysiol 116:2719–2733. doi:10.1016/j.clinph.2005.07.007 pmid:16253555
    OpenUrlCrossRefPubMed
  45. ↵
    1. Hoogenboom N,
    2. Schoffelen JM,
    3. Oostenveld R,
    4. Parkes LM,
    5. Fries P
    (2006) Localizing human visual gamma-band activity in frequency, time and space. Neuroimage 29:764–773. doi:10.1016/j.neuroimage.2005.08.043 pmid:16216533
    OpenUrlCrossRefPubMed
  46. ↵
    1. Hoogenboom N,
    2. Schoffelen JM,
    3. Oostenveld R,
    4. Fries P
    (2010) Visually induced gamma-band activity predicts speed of change detection in humans. Neuroimage 51:1162–1167. doi:10.1016/j.neuroimage.2010.03.041 pmid:20307670
    OpenUrlCrossRefPubMed
  47. ↵
    1. Hutcheon B,
    2. Yarom Y
    (2000) Resonance, oscillation and the intrinsic frequency preferences of neurons. Trends Neurosci 23:216–222. doi:10.1016/s0166-2236(00)01547-2 pmid:10782127
    OpenUrlCrossRefPubMed
  48. ↵
    1. Hutt A,
    2. Griffiths JD,
    3. Herrmann CS,
    4. Lefebvre J
    (2018) Effect of stimulation waveform on the non-linear entrainment of cortical alpha oscillations. Front Neurosci 12:376. doi:10.3389/fnins.2018.00376 pmid:29997467
    OpenUrlCrossRefPubMed
  49. ↵
    1. Iaccarino HF,
    2. Singer AC,
    3. Martorell AJ,
    4. Rudenko A,
    5. Gao F,
    6. Gillingham TZ,
    7. Mathys H,
    8. Seo J,
    9. Kritskiy O,
    10. Abdurrob F,
    11. Adaikkan C,
    12. Canter RG,
    13. Rueda R,
    14. Brown EN,
    15. Boyden ES,
    16. Tsai LH
    (2016) Gamma frequency entrainment attenuates amyloid load and modifies microglia. Nature 540:230–235. doi:10.1038/nature20587 pmid:27929004
    OpenUrlCrossRefPubMed
  50. ↵
    1. Jensen O,
    2. Kaiser J,
    3. Lachaux JP
    (2007) Human gamma-frequency oscillations associated with attention and memory. Trends Neurosci 30:317–324. doi:10.1016/j.tins.2007.05.001 pmid:17499860
    OpenUrlCrossRefPubMed
  51. ↵
    1. Jia X,
    2. Xing D,
    3. Kohn A
    (2013) No consistent relationship between gamma power and peak frequency in macaque primary visual cortex. J Neurosci 33:17–25. doi:10.1523/JNEUROSCI.1687-12.2013 pmid:23283318
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Keitel C,
    2. Keitel A,
    3. Benwell CS,
    4. Daube C,
    5. Thut G,
    6. Gross J
    (2019) Stimulus-driven brain rhythms within the alpha band: the attentional-modulation conundrum. J Neurosci 39:3119–3129. doi:10.1523/JNEUROSCI.1633-18.2019 pmid:30770401
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Kuffler SW
    (1953) Discharge patterns and functional organization of mammalian retina. J Neurophysiol 16:37–68. doi:10.1152/jn.1953.16.1.37 pmid:13035466
    OpenUrlCrossRefPubMed
  54. ↵
    1. Kujala J,
    2. Jung J,
    3. Bouvard S,
    4. Lecaignard F,
    5. Lothe A,
    6. Bouet R,
    7. Ciumas C,
    8. Ryvlin P,
    9. Jerbi K
    (2015) Gamma oscillations in V1 are correlated with GABA A receptor density: a multi-modal MEG and Flumazenil-PET study. Sci Rep 5:16347. doi:10.1038/srep16347 pmid:26572733
    OpenUrlCrossRefPubMed
  55. ↵
    1. Kumpulainen P
    (2020) tight_subplot(Nh, Nw, gap, marg_h, marg_w). Available at https://uk.mathworks.com/matlabcentral/fileexchange/27991-tight_subplot-nh-nw-gap-marg_h-marg_w.
  56. ↵
    1. Lacadie CM,
    2. Fulbright RK,
    3. Arora J,
    4. Constable RT,
    5. Papademetris X
    (2007) Brodmann areas defined in MNI space using new tracing tool in BioImage suite. In: Proceedings of the 14th Annual Meeting of the Organization for Human Brain Mapping, Chicago, IL, June, 10–14, Vol. 36, p 6494.
  57. ↵
    1. Lacadie CM,
    2. Fulbright RK,
    3. Rajeevan N,
    4. Constable RT,
    5. Papademetris X
    (2008) More accurate Talairach coordinates for neuroimaging using non-linear registration. Neuroimage 42:717–725. doi:10.1016/j.neuroimage.2008.04.240 pmid:18572418
    OpenUrlCrossRefPubMed
  58. ↵
    1. Lachaux JP,
    2. Rodriguez E,
    3. Martinerie J,
    4. Varela FJ
    (1999) Measuring phase synchrony in brain signals. Hum Brain Mapp 8:194–208. doi:10.1002/(SICI)1097-0193(1999)8:4<194::AID-HBM4>3.0.CO;2-C
    OpenUrlCrossRefPubMed
  59. ↵
    1. Lakatos P,
    2. Gross J,
    3. Thut G
    (2019) A new unifying account of the roles of neuronal entrainment. Curr Biol 29:R890–R905. doi:10.1016/j.cub.2019.07.075 pmid:31550478
    OpenUrlCrossRefPubMed
  60. ↵
    1. Lawrence MA
    (2016) ez: easy analysis and visualization of factorial experiments R package version 4.4-0.
  61. ↵
    1. Lee S,
    2. Jones SR
    (2013) Distinguishing mechanisms of gamma frequency oscillations in human current source signals using a computational model of a laminar neocortical network. Front Hum Neurosci 7:869. doi:10.3389/fnhum.2013.00869 pmid:24385958
    OpenUrlCrossRefPubMed
  62. ↵
    1. Lobo T,
    2. Brookes MJ,
    3. Bauer M
    (2021) Can the causal role of brain oscillations be studied through rhythmic brain stimulation? bioRxiv 2021.06.17.448493.
  63. ↵
    1. Lozano-Soldevilla D,
    2. ter Huurne N,
    3. Cools R,
    4. Jensen O
    (2014) GABAergic modulation of visual gamma and alpha oscillations and its consequences for working memory performance. Curr Biol 24:2878–2887. doi:10.1016/j.cub.2014.10.017 pmid:25454585
    OpenUrlCrossRefPubMed
  64. ↵
    1. Michalareas G,
    2. Vezoli J,
    3. Van Pelt S,
    4. Schoffelen JM,
    5. Kennedy H,
    6. Fries P
    (2016) Alpha-beta and gamma rhythms subserve feedback and feedforward influences among human visual cortical areas. Neuron 89:384–397. doi:10.1016/j.neuron.2015.12.018 pmid:26777277
    OpenUrlCrossRefPubMed
  65. ↵
    1. Morey RD,
    2. Rouder JN
    (2018) BayesFactor: computation of Bayes factors for common designs R package version 0.12-4.2 9.
  66. ↵
    1. Müller MM,
    2. Junghöfer M,
    3. Elbert T,
    4. Rochstroh B
    (1997) Visually induced gamma-band responses to coherent and incoherent motion: a replication study. Neuroreport 8:2575–2579. doi:10.1097/00001756-199707280-00031 pmid:9261830
    OpenUrlCrossRefPubMed
  67. ↵
    1. Muthukumaraswamy SD,
    2. Singh KD
    (2013) Visual gamma oscillations: the effects of stimulus type, visual field coverage and stimulus motion on MEG and EEG recordings. Neuroimage 69:223–230. doi:10.1016/j.neuroimage.2012.12.038 pmid:23274186
    OpenUrlCrossRefPubMed
  68. ↵
    1. Muthukumaraswamy SD,
    2. Singh KD,
    3. Swettenham JB,
    4. Jones DK
    (2010) Visual gamma oscillations and evoked responses: variability, repeatability and structural MRI correlates. Neuroimage 49:3349–3357. doi:10.1016/j.neuroimage.2009.11.045 pmid:19944770
    OpenUrlCrossRefPubMed
  69. ↵
    1. Nikolić D,
    2. Fries P,
    3. Singer W
    (2013) Gamma oscillations: precise temporal coordination without a metronome. Trends Cogn Sci 17:54–55. doi:10.1016/j.tics.2012.12.003 pmid:23287106
    OpenUrlCrossRefPubMed
  70. ↵
    1. Nolte G
    (2003) The magnetic lead field theorem in the quasi-static approximation and its use for magnetoencephalography forward calculation in realistic volume conductors. Phys Med Biol 48:3637–3652. doi:10.1088/0031-9155/48/22/002 pmid:14680264
    OpenUrlCrossRefPubMed
  71. ↵
    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
  72. ↵
    1. Oostenveld R,
    2. Fries P,
    3. Maris E,
    4. Schoffelen JM
    (2011) FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput Intell Neurosci 2011:156869. doi:10.1155/2011/156869 pmid:21253357
    OpenUrlCrossRefPubMed
  73. ↵
    1. Otero M,
    2. Prado-Gutiérrez P,
    3. Weinstein A,
    4. Escobar M-J,
    5. El-Deredy W
    (2020) Persistence of EEG alpha entrainment depends on stimulus phase at offset. Front Hum Neurosci 14:139. doi:10.3389/fnhum.2020.00139 pmid:32327989
    OpenUrlCrossRefPubMed
  74. ↵
    1. Pan Y,
    2. Frisson S,
    3. Jensen O
    (2020) Lexical parafoveal previewing predicts reading speed. bioRxiv 2020.10.05.326314.
  75. ↵
    1. Perry G,
    2. Hamandi K,
    3. Brindley LM,
    4. Muthukumaraswamy SD,
    5. Singh KD
    (2013) The properties of induced gamma oscillations in human visual cortex show individual variability in their dependence on stimulus size. Neuroimage 68:83–92. doi:10.1016/j.neuroimage.2012.11.043 pmid:23220427
    OpenUrlCrossRefPubMed
  76. ↵
    1. Pikovsky A,
    2. Kurths J,
    3. Rosenblum M,
    4. Kurths J
    (2003) Synchronization: a universal concept in nonlinear sciences. Vol. 12. Cambridge: Cambridge University Press.
  77. ↵
    1. Quintana DS,
    2. Williams DR
    (2018) Bayesian alternatives for common null-hypothesis significance tests in psychiatry: a non-technical guide using JASP. BMC Psychiatry 18:178. doi:10.1186/s12888-018-1761-4
    OpenUrlCrossRef
  78. ↵
    R Core Team (2020) R: a language and environment for statistical computing R foundation for statistical computing. Vienna: Austria.
  79. ↵
    1. Rager G,
    2. Singer W
    (1998) The response of cat visual cortex to flicker stimuli of variable frequency. Eur J Neurosci 10:1856–1877. doi:10.1046/j.1460-9568.1998.00197.x pmid:9751156
    OpenUrlCrossRefPubMed
  80. ↵
    1. Ray S,
    2. Maunsell JH
    (2010) Differences in gamma frequencies across visual cortex restrict their possible use in computation. Neuron 67:885–896. doi:10.1016/j.neuron.2010.08.004 pmid:20826318
    OpenUrlCrossRefPubMed
  81. ↵
    1. Ringach DL
    (2004) Mapping receptive fields in primary visual cortex. J Physiol 558:717–728. doi:10.1113/jphysiol.2004.065771 pmid:15155794
    OpenUrlCrossRefPubMed
  82. ↵
    1. Rodriguez E,
    2. George N,
    3. Lachaux JP,
    4. Martinerie J,
    5. Renault B,
    6. Varela FJ
    (1999) Perception's shadow: long-distance synchronization of human brain activity. Nature 397:430–433. doi:10.1038/17120 pmid:9989408
    OpenUrlCrossRefPubMed
  83. ↵
    1. Schadow J,
    2. Lenz D,
    3. Thaerig S,
    4. Busch NA,
    5. Fründ I,
    6. Rieger JW,
    7. Herrmann CS
    (2007) Stimulus intensity affects early sensory processing: visual contrast modulates evoked gamma-band activity in human EEG. Int J Psychophysiol 66:28–36. doi:10.1016/j.ijpsycho.2007.05.010 pmid:17599598
    OpenUrlCrossRefPubMed
  84. ↵
    1. Schwab K,
    2. Ligges C,
    3. Jungmann T,
    4. Hilgenfeld B,
    5. Haueisen J,
    6. Witte H
    (2006) Alpha entrainment in human electroencephalogram and magnetoencephalogram recordings. Neuroreport 17:1829–1833. doi:10.1097/01.wnr.0000246326.89308.ec pmid:17164673
    OpenUrlCrossRefPubMed
  85. ↵
    1. Self MW,
    2. Peters JC,
    3. Possel JK,
    4. Reithler J,
    5. Goebel R,
    6. Ris P,
    7. Jeurissen D,
    8. Reddy L,
    9. Claus S,
    10. Baayen JC,
    11. Roelfsema PR
    (2016) The effects of context and attention on spiking activity in human early visual cortex. PLoS Biol 14:e1002420. doi:10.1371/journal.pbio.1002420 pmid:27015604
    OpenUrlCrossRefPubMed
  86. ↵
    1. Shadlen MN,
    2. Movshon JA
    (1999) Synchrony unbound: a critical evaluation of the temporal binding hypothesis. Neuron 24:67–77. doi:10.1016/s0896-6273(00)80822-3 pmid:10677027
    OpenUrlCrossRefPubMed
  87. ↵
    1. Singer AC,
    2. Martorell AJ,
    3. Douglas JM,
    4. Abdurrob F,
    5. Attokaren MK,
    6. Tipton J,
    7. Mathys H,
    8. Adaikkan C,
    9. Tsai LH
    (2018) Noninvasive 40-hz light flicker to recruit microglia and reduce amyloid beta load. Nat Protoc 13:1850–1868. doi:10.1038/s41596-018-0021-x pmid:30072722
    OpenUrlCrossRefPubMed
  88. ↵
    1. Singer W
    (1999) Neuronal synchrony: a versatile code for the definition of relations? Neuron 24:49–65. doi:10.1016/s0896-6273(00)80821-1 pmid:10677026
    OpenUrlCrossRefPubMed
  89. ↵
    1. Singer W
    (2009) Distributed processing and temporal codes in neuronal networks. Cogn Neurodyn 3:189–196. doi:10.1007/s11571-009-9087-z
    OpenUrlCrossRefPubMed
  90. ↵
    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
  91. ↵
    1. Spaak E,
    2. de Lange FP,
    3. Jensen O
    (2014) Local entrainment of α oscillations by visual stimuli causes cyclic modulation of perception. J Neurosci 34:3536–3544. doi:10.1523/JNEUROSCI.4385-13.2014 pmid:24599454
    OpenUrlAbstract/FREE Full Text
  92. ↵
    1. Stenroos M,
    2. Hunold A,
    3. Eichardt R,
    4. Haueisen J
    (2012) Comparison of three- and single-shell volume conductor models in magnetoencephalography. Biomed Eng 57:311.
    OpenUrl
  93. ↵
    1. Tallon C,
    2. Bertrand O,
    3. Bouchet P,
    4. Pernier J
    (1995) Gamma-range activity evoked by coherent visual stimuli in humans. Eur J Neurosci 7:1285–1291. doi:10.1111/j.1460-9568.1995.tb01118.x pmid:7582101
    OpenUrlCrossRefPubMed
  94. ↵
    1. Tallon-Baudry C
    (2009) The roles of gamma-band oscillatory synchrony in human visual cognition. Front Biosci (Landmark Ed) 14:321–332. doi:10.2741/3246 pmid:19273069
    OpenUrlCrossRefPubMed
  95. ↵
    1. Tan HR,
    2. Gross J,
    3. Uhlhaas PJ
    (2016) MEG sensor and source measures of visually induced gamma-band oscillations are highly reliable. Neuroimage 137:34–44. doi:10.1016/j.neuroimage.2016.05.006 pmid:27153980
    OpenUrlCrossRefPubMed
  96. ↵
    1. Tass P,
    2. Rosenblum MG,
    3. Weule J,
    4. Kurths J,
    5. Pikovsky A,
    6. Volkmann J,
    7. Schnitzler A,
    8. Freund HJ
    (1998) Detection of n:m phase locking from noisy data: application to magnetoencephalography. Phys Rev Lett 81:3291–3294. doi:10.1103/PhysRevLett.81.3291
    OpenUrlCrossRef
  97. ↵
    1. Thut G,
    2. Schyns P,
    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:21811485
    OpenUrlCrossRefPubMed
  98. ↵
    1. Tiesinga PH
    (2012) Motifs in health and disease: the promise of circuit interrogation by optogenetics. Eur J Neurosci 36:2260–2272. doi:10.1111/j.1460-9568.2012.08186.x pmid:22805070
    OpenUrlCrossRefPubMed
  99. ↵
    1. Traub RD,
    2. Whittington M
    (2010) Cortical oscillations in health and disease. New York: Oxford University Press.
  100. ↵
    1. Traub RD,
    2. Whittington MA,
    3. Stanford IM,
    4. Jefferys JG
    (1996) A mechanism for generation of long-range synchronous fast oscillations in the cortex. Nature 383:621–624. doi:10.1038/383621a0 pmid:8857537
    OpenUrlCrossRefPubMed
  101. ↵
    1. Uhlhaas PJ,
    2. Singer W
    (2006) Neural synchrony in brain disorders: relevance for cognitive dysfunctions and pathophysiology. Neuron 52:155–168. doi:10.1016/j.neuron.2006.09.020 pmid:17015233
    OpenUrlCrossRefPubMed
  102. ↵
    1. Uhlhaas PJ,
    2. Pipa G,
    3. Lima B,
    4. Melloni L,
    5. Neuenschwander S,
    6. Nikolić D,
    7. Singer W
    (2009) Neural synchrony in cortical networks: history, concept and current status. Front Integr Neurosci 3:17. doi:10.3389/neuro.07.017.2009 pmid:19668703
    OpenUrlCrossRefPubMed
  103. ↵
    1. Van Pelt S,
    2. Fries P
    (2013) Visual stimulus eccentricity affects human gamma peak frequency. Neuroimage 78:439–447. doi:10.1016/j.neuroimage.2013.04.040 pmid:23611863
    OpenUrlCrossRefPubMed
  104. ↵
    1. Van Pelt S,
    2. Boomsma DI,
    3. Fries P
    (2012) Magnetoencephalography in twins reveals a strong genetic determination of the peak frequency of visually induced gamma-band synchronization. J Neurosci 32:3388–3392. doi:10.1523/JNEUROSCI.5592-11.2012 pmid:22399760
    OpenUrlAbstract/FREE Full Text
  105. ↵
    1. Varela F,
    2. Lachaux JP,
    3. Rodriguez E,
    4. Martinerie J
    (2001) The brainweb: phase synchronization and large-scale integration. Nat Rev Neurosci 2:229–239. doi:10.1038/35067550 pmid:11283746
    OpenUrlCrossRefPubMed
  106. ↵
    1. Veen B,
    2. Joseph J,
    3. Hecox K
    (1992) Localization of intra-cerebral sources of electrical activity via linearly constrained minimum variance spatial filtering. 1992 IEEE 6th SP Workshop on Statistical Signal and Array Processing, SSAP 1992, Conference Proceedings 44:526–529.
  107. ↵
    1. Von der Malsburg C
    (1999) The what and why of binding: the modeler's perspective. Neuron 24:95–104. doi:10.1016/S0896-6273(00)80825-9 pmid:10677030
    OpenUrlCrossRefPubMed
  108. ↵
    1. Wehr M,
    2. Laurent G
    (1996) Odour encoding by temporal sequences of firing in oscillating neural assemblies. Nature 384:162–166. doi:10.1038/384162a0 pmid:8906790
    OpenUrlCrossRefPubMed
  109. ↵
    1. Williams PE,
    2. Mechler F,
    3. Gordon J,
    4. Shapley R,
    5. Hawken MJ
    (2004) Erratum: entrainment to video displays in primary visual cortex of macaque and humans. J Neurosci 24:8278–8288. doi:10.1523/JNEUROSCI.2716-04.2004 pmid:15385611
    OpenUrlAbstract/FREE Full Text
  110. ↵
    1. Wilson HR,
    2. Cowan JD
    (1972) Excitatory and inhibitory interactions in localized populations of model neurons. Biophys J 12:1–24. doi:10.1016/S0006-3495(72)86068-5 pmid:4332108
    OpenUrlCrossRefPubMed
  111. ↵
    1. Womelsdorf T,
    2. Fries P,
    3. Mitra PP,
    4. Desimone R
    (2006) Gamma-band synchronization in visual cortex predicts speed of change detection. Nature 439:733–736. doi:10.1038/nature04258 pmid:16372022
    OpenUrlCrossRefPubMed
  112. ↵
    1. Xing D,
    2. Yeh CI,
    3. Burns S,
    4. Shapley RM
    (2012) Laminar analysis of visually evoked activity in the primary visual cortex. Proc Natl Acad Sci USA 109:13871–13876. doi:10.1073/pnas.1201478109 pmid:22872866
    OpenUrlAbstract/FREE Full Text
  113. ↵
    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. doi:10.1371/journal.pone.0013766 pmid:21072168
    OpenUrlCrossRefPubMed
  114. ↵
    1. Zhigalov A,
    2. Herring JD,
    3. Herpers J,
    4. Bergmann TO,
    5. Jensen O
    (2019) Probing cortical excitability using rapid frequency tagging. Neuroimage 195:59–66. doi:10.1016/j.neuroimage.2019.03.056
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 41 (31)
Journal of Neuroscience
Vol. 41, Issue 31
4 Aug 2021
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • 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.
No Evidence for Entrainment: Endogenous Gamma Oscillations and Rhythmic Flicker Responses Coexist in Visual Cortex
(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
No Evidence for Entrainment: Endogenous Gamma Oscillations and Rhythmic Flicker Responses Coexist in Visual Cortex
Katharina Duecker, Tjerk P. Gutteling, Christoph S. Herrmann, Ole Jensen
Journal of Neuroscience 4 August 2021, 41 (31) 6684-6698; DOI: 10.1523/JNEUROSCI.3134-20.2021

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
No Evidence for Entrainment: Endogenous Gamma Oscillations and Rhythmic Flicker Responses Coexist in Visual Cortex
Katharina Duecker, Tjerk P. Gutteling, Christoph S. Herrmann, Ole Jensen
Journal of Neuroscience 4 August 2021, 41 (31) 6684-6698; DOI: 10.1523/JNEUROSCI.3134-20.2021
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

  • entrainment
  • flicker
  • frequency tagging
  • gamma oscillations
  • magnetoencephalography
  • neuronal oscillations

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

Research Articles

  • Musical training facilitates exogenous temporal attention via delta phase entrainment within a sensorimotor network
  • Microglial Cytokines Mediate Plasticity Induced by 10 Hz Repetitive Magnetic Stimulation
  • Subgenual and hippocampal pathways in amygdala are set to balance affect and context processing
Show more Research Articles

Systems/Circuits

  • Subgenual and hippocampal pathways in amygdala are set to balance affect and context processing
  • Dorsolateral Striatum is a Bottleneck for Responding to Task-Relevant Stimuli in a Learned Whisker Detection Task in Mice
  • The Basolateral Amygdala Sends a Mixed (GABAergic and Glutamatergic) Projection to the Mediodorsal Thalamic Nucleus
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.