Abstract
θ-Band (4–12 Hz) activities in the frontal cortex have been thought to be a key mechanism of sustained attention and goal-related behaviors, forming a phase-coherent network with task-related sensory cortices for integrated neuronal ensembles. However, recent visual task studies found that selective attention attenuates stimulus-related θ power in the visual cortex, suggesting a functional dissociation of cortical θ oscillations. To investigate this contradictory behavior of cortical θ, a visual Go/No-Go task was performed with electroencephalogram (EEG) recording in C57BL/6J mice. During the No-Go period, transient θ oscillations were observed in both the frontal and visual cortices, but θ oscillations of the two areas were prominent in different trial epochs. By separating trial epochs based on subjects’ short-term performance, we found that frontal θ was prominent in good-performance epochs, while visual θ was prominent in bad-performance epochs, exhibiting a functional dissociation of cortical θ rhythms. Furthermore, the two θ rhythms also showed a heterogeneous pattern of phase-amplitude coupling with fast oscillations, reflecting their distinct architecture in underlying neuronal circuitry. Interestingly, in good-performance epochs, where visual θ was relatively weak, stronger fronto-visual long-range synchrony and shorter posterior-to-anterior temporal delay were found. These findings highlight a previously overlooked aspect of long-range synchrony between distinct oscillatory entities in the cerebral cortex and provide empirical evidence of a functional dissociation of cortical θ rhythms.
Significance Statement
Previous literature emphasized the pro-cognitive role of coherent oscillatory networks between distal brain regions, such as the fronto-visual θ synchrony. However, such a conceptual framework has been challenged as recent findings revealed distinct behavioral correlates of θ oscillations found in different cortical regions, especially in the frontal and visual cortices. Here, we show that frontal and visual θ represent distinct cortical processes and that the functional connectivity between them increases during sustained attention, especially when one of the two θ rhythms is relatively suppressed. The data presented here highlight a novel aspect of neural long-range synchrony between distinct cortical oscillators with distinct functional significance in task performance.
Introduction
In the mammalian brain, oscillatory networks in the lower frequency bands (<12 Hz) such as the θ band (4–12 Hz) are often thought to be a key mechanism of brain-wide interaction between distal regions (von Stein and Sarnthein, 2000). In particular, increased θ rhythmic activity of the frontal cortex is a hallmark of improved task performance and sustained attention, known to reflect the level of cognitive demand of the given task (Cavanagh and Frank, 2014; Clayton et al., 2015). In addition, θ oscillations of the frontal cortex have been known to form a phase-coherent network with task-related brain regions, including the sensory cortices (Benchenane et al., 2011). The emergence of a θ-oscillation network reflects active neuronal processing for cognitive functions such as sustained attention (Clayton et al., 2015), cognitive control (Cavanagh and Frank, 2014), and working memory (Sauseng et al., 2010). For instance, Liebe et al. (2012) reported local-field potentials and unit activities that were synchronous across the prefrontal and visual areas in the θ band during visual working memory tasks.
On the other hand, recent findings suggested that θ activities of the visual area are distinct from those of the frontal cortex and that visual θ is also associated with attention but in an opposite manner: stimulus-evoked θ power decreased with increased attention. For instance, Spyropoulos et al. (2018) recently reported that θ oscillations in the visual cortex elicited by visual stimuli decreased in amplitude when visual attention was directed. In addition, another recent study found arousal-reduced θ power in the high-order visual area (i.e., posterior parietal cortex) of the ferret brain (Stitt et al., 2018) despite its enhanced long-range θ synchrony with the prefrontal cortex during sustained attention reported by an earlier study (Sellers et al., 2016). In fact, a pioneering neuro-feedback study in the 1970s reported similar results of an association between occipital θ suppression and improved performance in a visual detection task (Beatty et al., 1974). There have also been several reports of enhanced θ in the posterior brain region being predictive of inattentiveness such as after sleep deprivation (Vyazovskiy and Tobler, 2005) or an attentional lapse (Mazaheri et al., 2009). Although the term “attention” of each study referred to different aspects of attention (e.g., internal or external attention; Chun et al., 2011) on different time scales (e.g., minute or second levels), evidence generally suggests that the employment of attention yields contrary patterns of θ oscillations in the frontal and visual cortices.
This attentional modulation of frontal and visual θ in opposite directions poses a question about the functionality of θ oscillations in different brain regions, especially in terms of their long-range functional connectivity. If cortical θ oscillations of different regions have different roles with opposing behavioral correlates, how could they form a coherent oscillatory network? A possible explanation of this counterintuitive interaction between cortical oscillators is that the oscillations of one region are suppressed to achieve synchronization with those of the other region. Here, we hypothesize that enhanced fronto-visual θ synchrony can occur during the suppression of frontal or visual θ.
To address this hypothesis, we designed a forced-paced (i.e., the opposite of self-paced) visual Go/No-Go task for head-fixed mice and analyzed the electroencephalogram (EEG) oscillations of each area. With this design, we aimed to capture the dynamic interaction between frontal and visual θ within a single protocol. In particular, in No-Go trials, we expected the recruitment of both frontal and visual θ oscillations, depending on the attentive state of the animal. This is because elaborate sustained attention (i.e., remaining alert while inhibiting a Go response) has been known to elicit control-related θ in the frontal cortex (Clayton et al., 2015; Sellers et al., 2016), while visual perception without attention has been reported to elicit θ oscillations in the visual area (Spyropoulos et al., 2018). By tracking short-term behavioral performance within a single experimental session, we investigated state-dependent changes of frontal and visual θ oscillations and their coupling with higher frequency oscillations (e.g., γ oscillations) to confirm their functional dissociation. Then, we characterized the property of functional connectivity between frontal and visual θ to show how distinct cortical oscillators are coupled to form a coherent long-range network.
Materials and Methods
All procedures were approved by the Institutional Animal Care and Use Committees at the corresponding author’s primary affiliation and complied with NIH guidelines. Detailed experimental procedures were as follows.
Animals and surgery
Five male C57BL/6 mice (four to eight weeks old at time of surgery, body weight M = 26.1 g, SD = 4.7 g) were used, maintained under a 12/12 h light/dark cycle (lights on at 8 A.M.) in a temperature-controlled and humidity-controlled environment. Food was available ad libitum, and water was mildly restricted (1 ml/d, additional water available as reward from task) from two weeks after surgery. Health status of the animals was assessed daily according to a systematic assessment protocol (Guo et al., 2014). All procedures were conducted in compliance with the Institutional Animal Care and Use Committee of corresponding author’s primary affiliation, conforming to the NIH Guide for Care and Use of Laboratory Animals (NIH Publication No. 86-23, revised 1985).
For stereotaxic surgery, mice were anesthetized with a ketamine and xylazine cocktail (120 and 6 mg/kg, respectively) by intraperitoneal injection and positioned on a stereotaxic apparatus (Model 957, Kopf Instruments). After shaving the head, an incision was made to expose the skull and implant sterilized screw (Asia Bolt) electrodes above the frontal (anterior-posterior, 1.9–2.1 mm; medial-lateral, ±2.1 mm from bregma) and visual cortices (anterior-posterior, –3.7 mm; medial-lateral, ±3.7 mm from bregma). For selection of recording sites, we referred to a previous study which targeted the mouse prelimbic area to record θ oscillations using surface electrodes in vivo (Sauer et al., 2015), and our pilot study using high-density surface EEG array (Choi et al., 2010) which showed that the selected mouse visual area exhibited most prominent motion-specific oscillatory responses to random-dot motion stimuli. Reference and ground electrodes were implanted on the occipital bone (Paxinos and Franklin, 2004). Dental cement (Vertex Self-Curing, Vertex Dental) was applied over the electrodes for fixation and two polycarbonate nuts (inner diameter 3 mm, Nippon Chemi-Con) were attached to the caudal edges of the cement, for use as anchors of the custom-built head-fixation apparatus. After surgery, mice were treated with antibiotics and analgesics.
Stimuli and apparatus
All experiments were conducted in a lightproof Faraday cage. Mice were restrained in a custom-built acrylic tube (3.5 cm in diameter) and head fixed by the polycarbonate nuts. For a reward, a water drop (3 μl) was released through a pressure equalizer tube (2 μm in diameter) with a programmable syringe pump (Fusion 200, Chemyx). Licking behavior was counted with a custom-built lickometer made of stainless-steel wire (3 μm in diameter) and a data acquisition system (cDAQ-9147, National Instruments). For punishment, a single puff of air was delivered to the animal’s cheek (air pressure = 15–20 psi, duration = 0.2 s), supplied with a compressed air tank and controlled by a solenoid valve. A schematic diagram of the experimental setup is shown in Figure 1A.
θ oscillations during Go/No-Go task in frontal and visual cortices. A, Schematic illustration of experimental setup. Mice performed Go/No-Go task in custom-built head restrainer. B, Experimental procedures. Task required subject to discriminate direction of random-dot motion (rightward: S+, downward: S–). Go response (i.e., licking) to S+ and S– motion resulted in reward and punishment, respectively. There was no particular reinforcing stimulus following No-Go response (i.e., no licking). The variation of the performance across days and within session was analyzed (Extended Data Fig. 1-1). C, Example EEG data of two trials. Upper, For S+ motion, licking response was counted as “hit,” and water reward was followed instantaneously. Motion stimulus was ceased at first lick (RT). Lower, For S– motion, No-Go response was counted as “CR.” Scale bars on bottom right indicate 500 ms and 300 µV. D, Summary of behavioral results. From all subjects (five mice), and all sessions (10 d), data of 8390 trials were collected. In this study, only the trials from CR (n = 2411) were analyzed. E, Time-frequency representation of grand-averaged event-related spectral perturbation (ERSP) of each cortical area. In both the frontal and visual cortices, prominent θ (4–12 Hz) oscillations were observed during 4 s of stimulus presentation (No-Go period). The stimulus-locked and response-locked spectrograms in all trials are summarized in Extended Data Figure 1-2. F, Temporal average (0.5–4 s) of θ oscillations in each cortical area. Numbers denote peak frequency of each channel. FA: false alarm, Fr: frontal cortex, Vis: visual cortex, Amptd: amplitude, Freq: frequency.
Extended Data Figure 1-1
Performances changing across- and within-session. A, Performance varies across sessions. Across the 10 d of sessions, overall behavioral performance (d’, M = 0.475, SE = 0.08) was significantly larger than chance level, Z = 5.845, p < 0.001. Learning effect, calculated by the difference of mean d’ between first 5 d (M = 0.399, SE = 0.07) and last 5 d (M = 0.551, SE = 0.21), was not significant, Z = 0.405, p = 0.686. Error bars represent ± 1 SEM. B, Performance varies within a session. Splitting trials of single session into three different groups, initial (first 1/3), middle (second 1/3), and final (last 1/3), change of behavioral performance within session was analyzed. In general, proportion of Go response trials (i.e., hit and FA) was higher at initial phase, and gradually decreased, suggesting the change of animals’ motivational level to water reward. Interestingly, despite the decrease of probability Go response, p(Go), overall probability of correct (i.e., hit and CR), p(correct), showed relatively small variation within session. C, Behavioral state changes across- and within-session. Changes of proportion of “attentive”-labeled trials (pc > 0.5) as a function of experimental session (left) and as a function of within-session progress (right) were analyzed. To calculate the effect size, one-way ANOVA (Kruskal–Wallis test) was performed. We found experimental session explained 17% of the variance, while the within-session progress explained 35% of the variance, suggesting the larger effect size of motivation-related factors relative to the learning-related factors. Error bars represent ± 1 SEM. D, Duration and transitions of behavioral state. Behavioral state showed transitions (i.e., transition from inattentive to attentive, and vice versa) and stationary phase with certain duration (inattentive: 6.60 trials in average, attentive: 7.81 trials in average; left). Using Markov-chain model with two-state process, this state transition was summarized in a state diagram (right). Download Figure 1-1, TIF file.
Extended Data Figure 1-2
Grand-averaged EEG responses during Go/No-Go task. A, Stimulus-locked grand-averaged amplitude spectrogram of four behavioral outcomes. Black dashed lines indicate the onset of visual stimulus. Hit and FA trials exhibited relatively large changes in oscillatory amplitudes, compared to those of Miss and CR trials. These large activities were considered as reinforcer-related activities (i.e., water reward, air puff punishment, motor movement, etc.) as they were locked to the moment of reaction (see C). B, Averaged (t = 0.5–4 s after motion onset) amplitude spectrum of spectrograms in A. C, Same as A, but response-locked. For the trials which do not contain behavioral response (CR, Miss), mean RTs of the session from hit and FA trials were used to align the epoch. D, Same as B, but response locked. Download Figure 1-2, TIF file.
For visual stimuli, a random-dot motion was generated and presented through MATLAB (MathWorks). The dot motion was composed of 100 small dots (full coherence, size = 0.5° in radius, field size = 60° in diameter, speed = 42°/s, dot life = 1.67 ± 0.33 s from a Gaussian distribution, orientation = 90° or 180°) and presented on a LCD screen (LG E1910PM-SN, 19 inches, 1280 × 1024 resolution), located 15 cm away from the animal’s eyes. As for the validity of using random-dot motions in the mouse model, previous studies have shown rodents are able to discriminate the direction of the random-dot motion (Douglas et al., 2006) and exhibit brain oscillatory responses similar to those found in humans (Han et al., 2017).
Training procedure
The behavioral training procedure was consisted of three stages: (1) habituation (7 d), (2) response shaping (2–3 d), and (3) conditioning (2–4 d). The habituation stage includes animal handling (10–15 min/d) and provides animals the opportunity to adapt to the head-fixation setup inside the experimental apparatus (30 min/d) without behavioral requirements. The response shaping stage was aimed to teach animals that water can be accessed by licking the tube. During this period, a water reward was given for every lick (lick port sensor time-out = 20 s). Once the animal started to robustly lick the tube, they advanced to the conditioning stage. The conditioning stage trained animals to associate a visual stimulus with a water reward. Each trial began with a monotone sound cue (10 kHz, 70 dB, 100 ms) and a white central fixation cross (size = 3°). After 2–3 s of the cue, a rightward moving-dots stimulus (S+, “Go” motion) was presented on the screen, which subsequently turned black for 30 s for the water-tube licking. Each stimulus did not last longer than 2 min, and 100 trials were given in 1 d. Our pilot study of reaction time (RT), validated through video recording, showed that intentional licking results in a typical RT distribution with a peak near 1 s, whereas spontaneous licking results in a scattered RT distribution over a 2 min time period. Also, we found that licking rates decreased over time, potentially due to thirst satiation. Thus, we considered the conditioning was successfully established when the RT of the first 80 trials reached a typical RT distribution. Overall, the percentage of mice licking in response to the visual stimulus on the last day of the conditioning stage was ∼70%.
Go/No-Go task
After mice learned to discriminate the direction of the random-dot motion on a screen, their task performance was tested by randomly presenting rightward (Go) or downward (“No-Go”) motion of random dots that stopped at licking or after 4 s when licking was not detected. Each mouse performed 10 experimental sessions on successive days. The task required the animals to discriminate the direction of the moving dots (Fig. 1B). Each trial began with a monotone sound cue (10 kHz, 70 dB, 100 ms) and a white central fixation cross (size = 3°). After 2–3 s of the cue, the stimulus was presented for 4 s with either a rightward (S+, Go motion) or downward (S–, No-Go motion) motion which was terminated when licking was detected. The direction of the motion of each trial was randomly determined (50%) with a restriction that the same stimulus was not presented more than four times in a row. A licking response during the S+, which was counted as a hit, resulted in instant reward. If the mouse did not take the reward within 20 s, the water in the tube was suctioned out. On the other hand, a licking response to the S– was counted as a FA, resulting in instant punishment and a 20-s timeout. Trials without licking responses were counted as a Miss or CR for the S+ or S–, respectively. To guide the animals to react more consciously, rather than reflexively, the response window started 0.5 s after the onset of motion (Fig. 1B). Intertrial intervals (it is; black screen) were 5–8 s long. Any spontaneous licking during the ITI or cue period before the response window reset the trial starting from the ITI.
The number of trials per session was determined by the animal’s performance that day (287 trials/d on average). Each session was terminated once the animal stopped licking for 30 successive trials, which we regarded as a lack of motivation for the task, and were excluded from the EEG analysis. We unexpectedly observed an obsessive licking behavior, during which the animal would indiscriminately lick the tube regardless of stimulus type. If animals licked for >30 trials in a row, those trials were excluded from EEG analysis. Overall, around 41% of trials were excluded.
Rating task performance
The ongoing task performance of mice was assessed by the ratio of correct answers [hit or correct rejection (CR)]. In human and animal studies, it has been known that the performance in cognitive tasks varies over time depending on the behavioral state, for instance due to an attentional lapse (Liebe et al., 2012; Esterman et al., 2013). To investigate the relationship between cortical θ activities and sustained attention, it was important to infer correctly the behavioral state, as either attentive or inattentive. This is particularly important in Go/No-Go tasks because mind wandering for a short period of time (i.e., attentional lapse) can be counted as a correct No-Go response. To rate each trial according to attentional state, the consistency of performance over time was analyzed. Because the span of attentional state in humans was reported to range between 2 and 7 min (Esterman et al., 2013) or exceed 10 min (Clayton et al., 2015), we assumed that short-term attentional state in mouse would be similar or shorter than that of human. Then, we assumed the level of attentional state is reflected through performance in the previous trials (in the order of a few minutes). Finally, we identified fluctuations in performance over 10 successive trials by calculating the probability of correct answers (pc; hit and CR), ranging from 0 (0 correct trials) to 1 (10 correct trials) with the chance level 0.5. pc of
k
-th trial was defined as follows,
(1)where the binary function,
, indicates the outcome of the trial as correct (1) or not (0). In our paradigm, 10 trials lasted ∼2–5 min. We subsequently classified the trials into attentive state (pc larger than the chance level) or inattentive state (pc equal to or smaller than the chance level).
EEG data acquisition and preprocessing
EEG responses were recorded at a 2-kHz sampling rate with microscrew electrodes implanted on the skull above the frontal and visual cortices using Grass 8-17C amplifiers (Grass Technologies, 60-Hz notch filters, gain = 50,000) and data acquisition system (cDAQ-9147). After bandpass filtering (1–200 Hz), the EEG data were processed on a trial-by-trial basis. To eliminate hemisphere-specific activities, EEG signals from bilateral electrodes were averaged in the time domain. To minimize the effect of early evoked activities (e.g., event-related potentials), only the EEG activities from 0.5 to 4 s after motion onset were used for statistical testing of oscillation amplitudes. We computed the mean and two-tailed 99% confidence intervals of the distribution of the overall amplitude spectrum, and the values lying outside of these confidence intervals were excluded. A time-frequency amplitude spectrogram was obtained by applying fast Fourier transform with a sliding Hanning window (bin size = 512 ms, step size = 100 ms) and by taking absolute real parts from the epochs, which resulted in a 1 Hz frequency resolution. For the time-frequency spectrogram visualization, the amplitudes of the pre-stimulus periods (2–0.5 s before stimulus onset) were averaged and subtracted in a frequency-wise manner.
Cross-frequency coupling (CFC) between θ and γ oscillations
To investigate CFC between θ and γ oscillations, we calculated the modulation index (MI) as described in Tort et al. (2010). Briefly: (1) raw EEG epochs were band passed (phase frequency: 4–12 Hz, amplitude frequencies: 6–200 Hz with 2.5-Hz step sliding window and 5-Hz width) with 5th order Butterworth filters. (2) The time series of the phase of θ oscillations denoted as
, and the time series of the amplitude of the γ oscillations denoted as
, were extracted via Hilbert transform. (3) The phase of the θ oscillations was divided into N bins from -π to π. In each bin, the average of
was calculated and then normalized by dividing the sum over the N bins and was denoted as Pj for the j-th bin (here, Pj equals the probability density function of the γ amplitude distribution over the θ phase ranging from 0 to 1). (4) An entropy index, H, was calculated by the definition of Shannon entropy, given by
(2)
if the γ oscillations are independent from θ phase, Pj is
; therefore
. (5) The statistical deviance of the measured distribution (P) from the uniform distribution was calculated by adopting Kullback–Leibler distance by the following formula,
(3)where
stands for the Kullback–Leibler distance. (6) Finally, MI was defined with a Z-score of observed entropy as the following formula indicates,
(4)
In the current study, we set N equal to 30 (i.e., bin size = 12°).
Generalized linear model (GLM) analysis of CFC
Single-trial GLM analyses were performed to investigate the increased level of θ-γ coupling during attentive state is modulated by either θ amplitude or attentive state, or both. The model was composed of the following equation,
(5)where
y
,
,
, and
ϵ
represent MI, θ amplitude, attentional state (inattentive = 0, attentive = 1), and error term, respectively. The term
represents a two-way interaction between θ amplitude and attentional state. GLM fitting was done by a MATLAB built-in function (fitlm.m).
Phase synchrony analysis of θ activities
To quantify the interaction between θ oscillations in the frontal and visual areas, phase synchrony was estimated by joint phase histograms and by calculating phase locking values (PLVs; Lachaux et al., 1999). To obtain joint phase histograms, a MATLAB built-in function for two-dimensional histogram (hist2.m) was used. To calculate PLVs, first, instantaneous phase of the θ oscillations was calculated via Hilbert transform and was denoted as
and
for the frontal and visual areas, respectively. Second, PLV was calculated in each stimulation period excluding the first 0.5 s (i.e., from
to
) as follows,
(6)where
. PLV ranges from 0 (i.e., absence of phase locking) and 1 (i.e., perfect phase locking).
A cross-correlation function (CCF)
of visual and frontal θ oscillations,
and
, respectively, was calculated in each trial excluding the first 0.5 s (i.e., from
to
), according to
(7)in which μ and σ indicate the mean and standard deviation. The time delay,
τ
, was obtained from;
(8)where
arg max
indicates the arguments of the maxima (i.e., the points that the CCFs are maximized). Positive values of
τ
indicate the temporal lead of visual θ over frontal θ and vice versa. The resolution of
τ
was 0.5 ms (2000-Hz sampling).
Statistical analysis
Non-parametric versions of t test (Wilcoxon signed-rank and rank-sum test), ANOVA (Kruskal–Wallis test), and circular-median test using CircStat Toolbox (Berens, 2009) were performed (
, two-tailed). For correlation analysis, Pearson’s method was used with the same α level. Sample size was not previously estimated for the experiments, and blinding or randomization was not needed in this work, as the design was decided by subjects’ performance.
Results
Changes of behavioral performance during Go/No-Go task
We trained mice (n = 5) to associate rightward motion (S+) with reward (water) and downward motion (S–) with punishment (air puff; Fig. 1A,B). The movement of visual stimulus was ceased at licking (i.e., black blank screen), as exemplified in Figure 1C. After training (for details, see Materials and Methods), we monitored the performance of the mice during the Go/No-Go test sessions for 10 successive days. During the test days, overall discrimination performance measured by D-prime (M ± SE, d’ = 0.475 ± 0.08) remained above the chance level (i.e., d’ = 0, Wilcoxon signed-rank test, Z = 5.845, p < 0.001; for details, see Extended Data Fig. 1-1). There was no statistically significant difference in the discrimination performance across the test days (Kruskal–Wallis test, χ2 = 1.179, p = 0.999). Moreover, a paired comparison between performance of the first half (1–5 d) versus the last half (6–10 d) did not yield any statistically significant difference (Wilcoxon signed-rank test, Z = 0.405, p = 0.686), indicating the absence of a learning effect (Extended Data Fig. 1-1A).
On the other hand, the animals showed clear changes in performance within a single session, probably due to the change of motivational factors (e.g., fatigue and quenched thirst; Extended Data Fig. 1-1B,C). The time-dependent patterns of performance were analyzed by dividing the session into early, middle, and late phases. When we classified the trials by hit (licking to S+), CR (no licking to S–), Miss (no licking to S+), and false alarm (FA; licking to S–; Fig. 1D, the ratios and numbers of classified trials), the rate of correct answers did not show any time dependency: sum of hit and CR ratios was 57.1%, 59.9%, and 57.4% in the early, middle, and late phases, respectively. Instead, the rates of licking decreased over time, shown through the sum of hit and FA ratios which was 64.7%, 51.6%, and 37.1% in the early, middle, and late phases, respectively, possibly due to the decreased motivation for the task (for details, see Extended Data Fig. 1-1).
Prominent θ rhythms in frontal and visual cortices during CR
All analyses of EEG data were focused on the time period from 0.5 s after stimulus onset until the cessation of stimulus (4 s after stimulus onset) during CR trials (n = 2411 from five mice). As expected, our results showed that during the stimulus period of CR trials, θ-band activities increased in both the frontal and visual cortices (Fig. 1E,F). Elicitation of γ activities in the low γ band (around 40 Hz) was also observed in the visual cortex (Fig. 1E,F). The peak frequency of baseline-corrected induced θ amplitudes was 7.62 Hz (FWHM = 6.26 Hz) and 6.64 Hz (FWHM = 4.96 Hz) in the frontal and visual cortices, respectively. The stimulus-locked and response-locked spectrograms in all trials are summarized in Extended Data Figure 1-2. Interestingly, although the behavior was the same (i.e., absence of licking response), Miss trials did not display such prominent θ-band activity (Extended Data Fig. 1-2).
Assessment of attentional state using short-term task performance
Because of the nature of the Go/No-Go task, some trial epochs without active engagement in the task (e.g., attentional lapse) could be counted as correct No-Go responses. To distinguish those by-accident CR trials from the intended, elaborate No-Go responses, we tracked the changes in attentional state within a single session. During a task, attentional states change and correspondingly, performance fluctuates over time (Clayton et al., 2015). In human and non-human primate studies, the short-term history of correctness (Liebe et al., 2012) or RT (Esterman et al., 2013) is often used to track the change of behavioral state. Similarly, we tracked the short-term proportion correct, pc, by defining it as the ratio of correct trials (i.e., hit and CR) across the 10 previous trials (see Materials and Methods for details). With a rough range of 2–5 min, pc monitors the momentarily fluctuating, varying trial-to-trial, attentional state. The distribution of pc for all CR trials exhibited a positively shifted Gaussian function (μ = 0.66, σ = 0.23,
= 0.98; Fig. 2A) whose mean was above the chance level (i.e., μ = 0.5; one-sample Z test, Z = 23.532, p < 0.001). Using the chance level as a threshold, we regarded the trial epochs as either attentive (i.e., pc above the chance level; pc
>
0.5) or inattentive states (i.e., pc
≤
0.5). For CR trials, the lowest value of pc was 0.2 and the highest value was 1, and attentive states were observed twice as frequently as inattentive states (1550 attentive states, 861 inattentive states). When we drew RT as a function of pc for licking cases (i.e., hit and FA), RT decreased in the case of hit trials (from 1.93 s for pc = 0.2–1.14 s for pc = 1, Z = 3.047, p = 0.002; Fig. 2B). On the other hand, RT in the case of FA trials did not show any dependence on pc. In this way each trial was labeled with pc, which successfully reflected the change of attentional state, for further analysis.
Task performance and θ amplitude. A, Distribution of pc (proportion correct) in all CR trials (n = 2411). pc was calculated across 10 previous trials. With a threshold of chance level (pc = 0.5), each trial was classified as either attentive (pc > 0.5) or inattentive (pc ≤ 0.5) state. B, Change of mean RT as function of pc. The RT of hit trials decreased as pc increased (Pearson’s r = –0.051, p < 0.05), suggesting the change of attentional level within the hit trial epochs. C, Scatterplot of frontal/visual θ amplitude with color coding of pc. Each gray dot represents the observed θ amplitude value of a single CR trial. Gray solid lines indicate the probability density function of each channel. Along with the changes of pc, the mean amplitudes of frontal and visual θ showed a negative relationship. D, Mean amplitude of frontal/visual θ as function of pc. The linear correlation between θ amplitude and pc was positive in the frontal (r = 0.076, p < 0.001) and negative in the visual (r = –0.099, p < 0.001), suggesting their distinct behavioral correlates. E, Time-frequency representation of oscillatory amplitudes highlighting state-dependent and region-dependent change of θ amplitudes. Black dashed lines indicate the start/end of the visual stimuli. Error bars represent ±1 SD (C) and SE (B, D) of the means. pc: probability correct. The changes in power spectral density and its peak frequency were calculated for different values of pc (Extended Data Fig. 2-1).
Extended Data Figure 2-1
Baseline-uncorrected amplitude- and cross-spectral densities. A, Mean amplitude- or cross-spectrum of CR trials (n = 2,411, t = 0.5–4 s) calculated from the frontal. Error bars represent ± 1 SEM. B, Same as A, from the visual. C, Same as A, from cross-spectrum of the frontal and visual. D, Mean frontal θ (4–12 Hz) amplitude. Error bars represent ± 1 SEM. E, Same as D, from the visual. F, Same as D, from cross-spectrum of the frontal and visual. G, Mean peak frequency of the frontal θ-amplitude spectrum. Peak values were calculated in trial-by-trial basis. Colored numbers illustrate the duration of one cycle of each frequency. Error bars represent ± 1 SEM. H, Same as G, from the visual. I, Same as G, from cross-spectrum of the frontal and visual. Download Figure 2-1, TIF file.
Opposite effects of attentional state on frontal and visual θ amplitudes
The primary aim of the experiment was to articulate the ongoing θ-activity signatures with respect to attention, which were reported to act in an opposite way in the frontal (Clayton et al., 2015) and visual (Spyropoulos et al., 2018) cortices. Thus, we were interested in the differences in ongoing θ activities between two attentional states. The scatterplot in Figure 2C shows a positive correlation between frontal and visual θ amplitudes. On the other hand, the mean of each pc presented as a colored error bar shows a negative correlation, with a movement in their distribution from the lower-right to the upper-left corner as pc is increased. This opposite behavior of frontal and visual θ is manifested in the plot of θ amplitudes as a function of pc (Fig. 2D). In the frontal cortex, the correlation between θ amplitude and pc was positive, reaching its maximum at pc = 1 (Pearson’s r = +0.08, p < 0.001). In contrast, the relationship was negative for visual θ, reaching its minimum at pc = 1 (Pearson’s r = –0.10, p < 0.001). The average θ amplitudes for attentive versus inattentive, analyzed in the time-frequency domain, demonstrates the opposite effect of attention on the frontal and visual θ amplitudes (Fig. 2E). In the frontal cortex, θ amplitude in attentive states was significantly larger compared to that in inattentive states (96.5 ± 0.8 μV in inattentive, 99.1 ± 0.7 μV in attentive; Wilcoxon rank-sum test, Z = 3.189, p < 0.001). On the other hand, in the visual cortex, the amplitude of θ in attentive states was significantly smaller compared to that in inattentive states (86.7 ± 1.0 μV in inattentive, 82.2 ± 0.9 μV in attentive; Wilcoxon rank-sum test, Z = –2.582, p < 0.001). These contrasting correlations suggest that frontal and visual θ activities might have different functional roles during the engagement of cognitive processes required for the No-Go response, raising a question about the neurologic signature of sustained attention.
Distinct pattern of phase-amplitude coupling with fast oscillations
To further investigate the role of θ in each region, we performed the phase-amplitude CFC analysis, whose characteristics are suggested to reflect the architecture and function of underlying neuronal circuits (Hyafil et al., 2015). The amplitude maps of fast oscillations plotted with respect to θ phase revealed strong but different CFC patterns across different pc values and regions (Fig. 3A). For example, in the high-performance epochs with higher pc, the amplitude of fast oscillations was strongly modulated by the phase of θ oscillation in the frontal cortex compared to in the low-performance epochs. To analyze the effect of behavioral state on CFC, we calculated the MI (Tort et al., 2010) as an indicator of CFC strength. As depicted in Figure 3B, the MI in the frontal cortex showed frequency-specific modulation of behavioral state, while the visual cortex did not show such modulation. Between θ and low γ (20–40 Hz), the attentional enhancement in the frontal CFC was statistically significant (inattentive: 6.1 × 10−4 ± 0.1 × 10−4, attentive: 6.6 × 10−4 ± 0.1 × 10−4; Wilcoxon rank-sum test, Z = 2.187, p = 0.029; Fig. 3C). For the faster frequency band (high γ, 80–160 Hz) in the frontal cortex, attentional reduction was observed (inattentive: 7.4 × 10−4 ± 0.2 × 10−4, attentive: 7.0 × 10−4 ± 0.1 × 10−4; Wilcoxon rank-sum test, Z = –2.546, p = 0.011). However, in the visual cortex, the attentional effect on CFC was not statistically significant for either low γ (Z = –1.443, p = 0.149) or high γ (Z = 1.121, p = 0.262). It should be taken into consideration that the CFC changed due to changing θ amplitudes, rather than behavioral performance, as previous studies have shown that CFC can be either dependent on the amplitude of ongoing θ oscillations (Canolty et al., 2006; Tort et al., 2008) or independent from it (Kendrick et al., 2011, their Fig. S9). To disentangle this possible confounding factor, an additional single-trial GLM analysis was performed. The results suggest that both factors, θ amplitude and attentional state, contribute to the change of CFC (Extended Data Fig. 3-1).
Task performance and θ-γ phase-amplitude coupling. A, Color-coded fast-frequency amplitudes as function of θ phase in frontal (top) and visual (bottom) and as function of pc (left to right). B, C, MI for local θ-γ coupling as function of pc (left: frontal, right: visual). In the frontal, the effect of attention on CFC was frequency-specific (increased CFC with low γ, decreased CFC with high γ). The attentional modulation of the visual coupling was not significant. The dependency of CFC on θ power or attentional state or their interaction was further analyzed using a general linear model, showing both θ power and its interaction with attentional state are significant independent parameters for CFC values (Extended Data Fig. 3-1). D, Grand-averaged amplitudes of fast oscillations as function of θ phase, showing noticeable regional difference of preferred θ phase (left: frontal, right: visual). E, Trial histogram of preferred θ phase of 20- to 40-Hz oscillations (left). Hexagons and dotted lines denote the peak and mean of each polar histogram, respectively. Mean preferred angle of each area (right). The distances between preferred angles over sites were both significant. Error bars represent ±1 SEM; *p < 0.05, ***p < 0.001. INA: inattentive, ATT: attentive.
Extended Data Figure 3-1
GLM analysis of CFC. To dissect the confounding effect of θ and attention on MI, MIs were further analyzed through GLM with three predictors: θ amplitude, attentive state, and two-way interaction of θ amplitude and attentive state. As a result, we found statistically significant main effect of θ amplitude on θ-low-γ coupling in both regions (frontal: t = 14.529, p < 0.001, visual: t = 9.902, p < 0.001) and on θ-high-γ coupling in the visual (t = 2.648, p < 0.01). The effect of θ amplitude on frontal θ-high-γ coupling did not reach statistical significance (t = 1.747, p = 0.081). Also, the effect of attentional state on MI was statistically significant only for visual θ-low-γ coupling (t = –2.864, p < 0.01) with non-significant effect on the other couplings (|t|s < 1.688, ps > .091). Most importantly, the two-way interaction between θ amplitude and attentional state was statistically significant for θ-low-γ coupling in both regions (ts > 2.294, ps < 0.05), but not for θ-high-γ coupling (|t|s < 1.301, ps > .193). Such two-way interaction with positive slope value (ß = 1.0 × 10–6 for frontal, ß = 1.3 × 10–6 for visual) suggests the effect of θ amplitude on θ-γ coupling becomes stronger with increased level of attention, especially for θ-low-γ coupling in both regions. In short, the change of CFC observed in this study can be best explained by a synergetic effect of θ amplitude and attentional state, rather than by single factor. Download Figure 3-1, PDF file.
These region-specific behaviors of CFC were further analyzed by investigating the θ phase at which γ amplitudes peaked (i.e., peak preference). This was done to check that the two θ oscillations are qualitatively distinct, as previous studies have shown that such peak phase difference implies differences in the underlying circuit motif (Voloh and Womelsdorf, 2016, their Fig. 4). When these phase angles were pooled separately for each region, their distributions differed (Fig. 3D, polar plots). In the frontal cortex, the preferred θ phase (i.e., a phase bin of θ where maximum low γ amplitude was detected) was observed in the middle of the rising phase (109.7 ± 2.0°), whereas in the visual cortex, the preferred θ phase was observed near the start of the falling phase (191.1 ± 1.7°). The distributions of the preferred θ phase of low γ were statistically different for the frontal cortex and the visual cortex (circular Kruskal–Wallis test, χ2 = 251.140, p < 0.001). The angle distance of the preferred θ between them was 81.4°, corresponding to 32.3 ms at 7 Hz. In addition to the results from distinct behavioral modulation (Fig. 2C–E), the result of our CFC analysis presented here provides more reliable evidence that the θ rhythms of each area have distinct functional roles in the cortical microcircuit.
Synchrony and lags of θ oscillations. A, Two-dimensional representation of joint phase histogram of frontal and visual θ oscillations as function of pc. The white color of each histogram indicates the higher probability of colocalization of the phase of two oscillations, suggesting a higher level of synchrony between two oscillations. The inset text indicates the maximum value of colocalization probability in each pc bin. B, Mean fronto-visual phase synchrony calculated by PLV as function of pc. The phase synchrony was higher for the attentive state than the inattentive state. Extended Data Figure 4-1A illustrates the phase-locking behaviors of two rhythms appeared as a transient locking. The phase synchrony is its maximal at the moment of frequency locking of two θ rhythms (Extended Data Fig. 4-1B), which is also dependent on the attentional state (Extended Data Fig. 4-1C). C, Cross-correlation analysis of example EEG epoch. Top, To calculate temporal lag (τ), raw time traces of EEG data from two channels were bandpass filtered in the θ band. Bottom, After filtering, the CCF was obtained and τ was defined as a temporal lag where the correlation coefficient reaches its maximum. Positive and negative values of τ indicate a lead of the visual and frontal θ, respectively. D, Trial histograms of τ between frontal and visual θ. The dots denote the peak of each histogram. Color-coded texts indicate the peak position from the lowest pc (green) and to the highest pc (red) bin. E, Mean τ as function of pc. The τ was smaller in the attentive state than in the inattentive state. Error bars represent ± 1 SEM; *p < 0.05, ***p < 0.001. INA: inattentive, ATT: attentive.
Extended Data Figure 4-1
Transient frequency matching during fronto-visual θ synchrony. A, Example single-trial data showing transient frequency matching and increased phase synchrony. To calculate frequency difference and phase synchrony, raw data were converted into Hilbert instantaneous angle after band-pass filtering. Frequency difference was obtained by taking derivate of unwrapped angle difference between two channels’ data. Fr: frontal, Vis: visual, Ang: angle, Diff: difference, Freq: frequency. B, Two-dimensional histogram of instantaneous frequency difference and phase synchrony during the No-Go period of all trials (irrespective of attentional state). High synchrony was observed when frequency difference is low. Note that synchrony value was calculated with 1-s width and 10-ms sliding temporal window, and instantaneous frequency difference was time averaged with same size sliding window.
: instantaneous frequency of frontal θ,
: instantaneous frequency of visual θ. C, Histogram of instantaneous frequency difference, drawn separately for each pc group. Probability of having transient frequency matching was higher during the attentive than during the inattentive state. Download Figure 4-1, TIF file.
Enhanced fronto-visual functional connectivity during attentive state
Phase synchrony analysis allows the examination of signals that are cooperative between brain areas. To investigate the dynamic interplay of the fronto-visual θ network as attention fluctuates, we first examined the phase relationship between two θ rhythms using a joint phase histogram (see Materials and Methods for details). Figure 4A highlights the systematic change of the instantaneous phase relationship as a function of pc. In general, the joint phase of frontal and visual θ showed the stronger colocalization along the diagonal axis regardless of behavioral state, which indicates their intimate phase relationship. Furthermore, the colocalization along the diagonal axis became more robust as pc increased, suggesting a higher level of phase synchrony between the two θ rhythms. The PLV (Lachaux et al., 1999), an amplitude-independent parameter, elucidates the brain state-dependent synchrony between frontal and visual θ oscillations (Fig. 4B). Compared to the inattentive state, the PLV was significantly higher in the attentive state (PLV = 0.614 ± .004 in inattentive, 0.624 ± .003 in attentive; Wilcoxon rank-sum test, Z = –2.313, p = 0.021). Our result of a fronto-visual θ synchrony associated with task performance is consistent with a previous finding in a non-human primate study (Liebe et al., 2012). While an opposite relationship between frontal and visual θ amplitudes was found in association with the behavioral state, a stronger fronto-visual long-range functional network was established, especially when frontal θ was relatively increased and when visual θ was relatively suppressed.
The asymmetrical joint phase distribution (Fig. 4A, bottom-right corner of each joint phase histogram) raised the possibility of a driver–responder relationship between frontal and visual θ. Therefore, we investigated the time lag, τ, of frontal θ with respect to visual θ by finding the peaks at the CCFs of the two θ rhythms. Figure 3C illustrates an example CCF of frontal and visual θ rhythms. τ was obtained at each trial and the distributions of τ were plotted for different pc values (Fig. 4D). Among all trials, a majority of the trials (84.61%) had positive τ, indicating that visual θ predominantly precedes frontal θ. In addition, as pc increased, the distribution of τ shifted to shorter phase lag. The peak τ in Figure 4D systematically decreased as the state changed from inattentive to attentive, with values of 14.83 and 4.33 ms at pc of 0.2 and 1, respectively. τ was 8.25 ± 0.24 and 6.97 ± 0.17 ms for inattentive and attentive states, respectively (Wilcoxon rank-sum test, Z = 4.832, p < 0.001). In addition, we found that the peak frequency of the cross-spectral density function was significantly higher for the attentive state (5.84 ± 0.03 Hz) compared to the inattentive state (5.76 ± 0.04 Hz; Wilcoxon rank-sum test, Z = 2.052, p = 0.040; Extended Data Fig. 2-1). The increase of synchrony and shift to shorter phase lag during the attentive state suggests the contribution of θ oscillations to task engagement and sustained attention via their long-range functional network.
It may seem odd that we calculated the phase synchrony after demonstrating that the two oscillators are distinct in peak frequency, since two oscillations must have transient periods of same frequency for phase locking to occur. Indeed, when we calculated synchrony for various instantaneous θ frequency, we found there were times during which the two θ rhythms oscillated in sync (Extended Data Fig. 4-1). Frequency difference was low when high phase synchrony was observed, and the probability of converging frequency difference (to 0) was found to be modulated by attention, which explains why PLV was higher during the attentive state.
Discussion
We found that the frontal and visual θ induced by the task had functionally distinct features. Whereas the frontal θ amplitude was positively correlated with behavioral performance, the visual θ amplitude was negatively correlated. In the frequency domain, the peaks of the two θ oscillations differed by almost 1 Hz, and their phases modulated fast cortical oscillations in different ways. Nonetheless, the frontal-visual θ synchrony increased as a function of performance despite the weakened amplitude of visual θ, implying a critical role of phase-locked oscillations in successful task engagement. It is unlikely that frontal and visual θ have common driving mechanisms; instead, it is more likely that the two θ have independent origins but are weakly coupled to each other. We should note that there is a distinction between phase locking (i.e., coupled oscillation with a constant phase difference) and phase trapping (i.e., un-locked oscillation with the same frequency) in nonlinear oscillator systems (Aronson et al., 1990). Indeed, the frontal and visual θ in our observation showed a phase transition from a stable phase-trapping system to a non-symmetric phase-locking system as the attention level increased. Although the differences in the values were small, we observed that there were abrupt changes in the interaction parameters (e.g., max in joint phase histogram, synchrony, and lag in Fig. 4) near pc = 0.8. Furthermore, the behavioral parameters (i.e., RT in Fig. 2B) and the neural dynamics parameters (i.e., θ amplitude in Fig. 2D; θ-γ coupling in Fig. 3A,C) presented their characteristic features near pc = 0.8. It is highly possible that the θ oscillations changed their dynamic trajectory at the highly task-engaged brain state in their phase portrait; however, determining whether this transition contributed to sustained attention or vice versa was beyond the scope of this study.
Our phase dynamics analysis showed that visual θ led frontal θ throughout the recording sessions. At the same time, frontal, not visual, θ increased its amplitude at higher pc accompanied by a reduction in its phase delay to visual θ. This confounding behavior at higher pc could be an influence of top-down directed frontal θ for cognitive control such as motor inhibition (Cavanagh and Frank, 2014) or phase resetting of θ oscillations by a stimulus in a highly task-engaged brain state (Rizzuto et al., 2003). In terms of the Communication through Coherence theory (Fries, 2015), the tight interplay of two oscillations renders neuronal communication more effective, precise, and selective. It is possible that a new common source arose around pc = 0.8 to enhance their synchrony and that the newly adjusted phase lag was simply differential conduction delay. However, the different patterns of θ-γ coupling suggest this was not the case. Recently, Zhang et al. (2018) reported that human cortical oscillations were regionally distinctive and had the properties of weakly coupled oscillations rather than those of oscillations driven by a common source or feedback loop. Although attentional modulation was not included in their study, it was consistent with our observations, in that the visual rhythm was faster than the frontal rhythm by ∼1 Hz and propagated to the anterior cortex. However, there was a mismatch in the propagation delay in the human visual to frontal cortex, which was tenfold slower compared to our observation in mice or in monkeys (Gregoriou et al., 2009). The leading of visual θ might reflect the bottom-up relay of sensory information, and the systematic reduction of the delay is functionally relevant to its purpose of faster relaying of sensory information. The delay is interpreted as a designated temporal shift for the effective relay of information, transmitting a presynaptic spike at the peak depolarization phases of postsynaptic neurons (Nowak and Bullier, 1997; Gregoriou et al., 2009). Metaphorically speaking, to selectively improve the traffic in one direction, it may be optimal for the green lights of successive crossways to be slightly delayed considering the amount of time it takes to drive through, rather than operating with zero-lag synchrony.
Then what is the generation mechanism of frontal θ observed here? Considering that our task requires an effort to withhold impulsive actions, frontal θ might originate from the interaction of the cortico-basal ganglia circuit (i.e., anterior-cingulate/medial-frontal and subthalamic nucleus) related to the conflict-related decision process (Cavanagh et al., 2011; Cavanagh and Frank, 2014; Cohen, 2014). Or the frontal θ may reflect communication between hippocampal structures and the frontal cortex (Jones and Wilson, 2005; Siapas et al., 2005), as the task required the subjects to maintain visuospatial information (i.e., direction of motion). Regarding the role of frontal θ oscillations in sustained attention, a large and consistent body of literature has discussed the top-down role of frontal θ in cognitive control (Cavanagh and Frank, 2014) and sustained attention (Clayton et al., 2015). All these works included task-engaged mental processes (e.g., keeping one’s goal in mind, overriding the prepotent response, monitoring the goal-related sensory cues). Our results provide further evidence for this top-down role of frontal θ oscillations, as we found that frontal θ increased during sustained attention in mice. Importantly, phase synchrony between frontal and visual θ increased in high-performance trials in a critical manner. Thus, frontal θ appears to be behaviorally relevant as its top-down control mode. It should be noted that the substantial work conducted in the frontal cortex has found γ oscillations to play an important role in pro-cognitive functions (Cho et al., 2015; Kim et al., 2016). On the neuronal spike levels, the spiking activities of the frontal cortex in the high-γ band were locked to θ oscillations in the frontal and posterior parietal cortices during sustained attention (Sellers et al., 2016). In our study, we found effective modulation of low γ, but not high γ, by frontal θ using CFC analysis. Continued work with neuronal recording or optogenetic interrogation will be needed to clarify the specific role of frontal θ on neuronal activities during cognitive processing.
The counteracting behavior of visual θ compared to frontal θ was counterintuitive, but it was consistently and significantly observed across animals and days. This phenomenon of “attentional reduction” in visual θ is in agreement with recent findings in monkeys (van Kerkoerle et al., 2014; Spyropoulos et al., 2018) and ferrets (Stitt et al., 2018). In humans, it has been reported that α activities (8–12 Hz) in the occipital cortex reflect functional inhibition of task-irrelevant cortical areas (Jensen and Mazaheri, 2010). To the best of our knowledge, wakeful relaxation related α activities have not been reported in mice unlike other cross-species frequency-preserved rhythms (Buzsáki et al., 2013). In addition, compared to baseline activity, visual θ in our study increased during stimulus presentation (Fig. 2E), whereas α usually decreases (Fries et al., 2008), which is similar to the finding of Spyropoulos et al. (2018). Therefore, it seems more plausible to regard the oscillatory activities observed in this study as visual θ rather than occipital α. However, the functional significance of such attentional reduction in the cortical microcircuit remains unknown.
What is the role of θ rhythmicity in the visual system? Recently, it has been suggested that the θ rhythmicity in the perceptual domain reflects the periodic sampling process of the sensory system (i.e., rhythmic sampling; VanRullen, 2016, 2018). In this framework, attention periodically fluctuates around the θ band to sample sensory stimuli, resulting in θ-rhythmic fluctuations of relevant neural activities and behavioral responses (Busch and VanRullen, 2010; Fiebelkorn et al., 2018; Kienitz et al., 2018). Still, the suppression of visual θ by attention is not explained, particularly when we consider that θ rhythms reflect the attentional sampling process. According to Spyropoulos et al., the attentional reduction of visual θ was suggested to reflect “more continuous processing” of attended stimuli, and a shift to the faster rhythm was predicted (Spyropoulos et al., 2018). Moreover, Harris et al., reported that the phase modulation of visual θ on visual detection performance was more effective for a non-attended than attended target in humans (Harris et al., 2018). These results suggest that visual θ might be in charge of processing outside the attentional spotlight rather than involved in neural processing in the focal attentional field.
While a growing number of studies have revealed the modulation effects of visual θ on neural excitability, the neural basis for the generation of visual θ remains unknown (VanRullen, 2018). Recently, a predominant feedforward directionality of visual θ was found in monkeys, from lower to higher levels of the visual cortex, and to the temporal cortex (Spyropoulos et al., 2018). In addition, a colocalization of stimulus-induced γ oscillations and visual θ oscillations was revealed (van Kerkoerle et al., 2014). Considering that the feedforward directionality of γ oscillations in the primary sensory cortex is relatively well known, it is believed that visual θ reflects bottom-up sensory processing rather than top-down influences from higher-order visual areas. Moreover, simultaneous recording of the visual thalamic area and the posterior-parietal cortex in the ferret brain also demonstrated a strong feedforward directionality in the θ frequency band (Stitt et al., 2018), but it is not clear whether it is generated from the thalamo-cortical network or instead is locally driven.
If visual θ propagates in a feedforward direction, how can it be modulated by top-down processes such as attention? One recent study using ferret thalamus recordings suggests an answer to this question. It found that spiking activities in the pulvinar are modulated by attention, and these neurons signal in the θ frequency band during sustained attention (Yu et al., 2018). This area exhibits bidirectional anatomic connections with higher-order visual cortices (Jones, 2001), and even with the prefrontal cortex (Romanski et al., 1997). As there has also been a report that unit activities in the primate visual area V4 are modulated by the external θ driver in the prefrontal cortex (Liebe et al., 2012), it seems highly plausible that the top-down modulation of θ is generated via a cortico-thalamic input, which in turn leads to a change of thalamo-cortical connectivity. Indeed, attentional modulation of neuronal activities in the visual cortices has attracted a wide range of interest among cognitive scientists. Considering the genetic techniques uniquely applicable to the mouse model, future studies may take advantage of mouse models to understand the functional architecture and underlying circuitry of visual θ rhythmicity in association with its attentional modulation.
Footnotes
The authors declare no competing financial interests.
This work was supported by the National Research Council of Science and Technology of Korea for the project Development of Solution for Diagnosis, Treatment, and Care System of Dementia Grant CRC-15-04-KIST, the K-DARPA project (2V0713) and by the National Research Foundation of Korea (NRF) Grant 2017R1A2B3012659.
↵* H.-B.H. and K.E.L. are first authors.
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.