Learning alters theta amplitude, theta-gamma coupling and neuronal synchronization in inferotemporal cortex
© Kendrick et al; licensee BioMed Central Ltd. 2011
Received: 14 January 2011
Accepted: 9 June 2011
Published: 9 June 2011
How oscillatory brain rhythms alone, or in combination, influence cortical information processing to support learning has yet to be fully established. Local field potential and multi-unit neuronal activity recordings were made from 64-electrode arrays in the inferotemporal cortex of conscious sheep during and after visual discrimination learning of face or object pairs. A neural network model has been developed to simulate and aid functional interpretation of learning-evoked changes.
Following learning the amplitude of theta (4-8 Hz), but not gamma (30-70 Hz) oscillations was increased, as was the ratio of theta to gamma. Over 75% of electrodes showed significant coupling between theta phase and gamma amplitude (theta-nested gamma). The strength of this coupling was also increased following learning and this was not simply a consequence of increased theta amplitude. Actual discrimination performance was significantly correlated with theta and theta-gamma coupling changes. Neuronal activity was phase-locked with theta but learning had no effect on firing rates or the magnitude or latencies of visual evoked potentials during stimuli. The neural network model developed showed that a combination of fast and slow inhibitory interneurons could generate theta-nested gamma. By increasing N-methyl-D-aspartate receptor sensitivity in the model similar changes were produced as in inferotemporal cortex after learning. The model showed that these changes could potentiate the firing of downstream neurons by a temporal desynchronization of excitatory neuron output without increasing the firing frequencies of the latter. This desynchronization effect was confirmed in IT neuronal activity following learning and its magnitude was correlated with discrimination performance.
Face discrimination learning produces significant increases in both theta amplitude and the strength of theta-gamma coupling in the inferotemporal cortex which are correlated with behavioral performance. A network model which can reproduce these changes suggests that a key function of such learning-evoked alterations in theta and theta-nested gamma activity may be increased temporal desynchronization in neuronal firing leading to optimal timing of inputs to downstream neural networks potentiating their responses. In this way learning can produce potentiation in neural networks simply through altering the temporal pattern of their inputs.
The functions of both low and high frequency oscillations in the brain are the subject of considerable speculation . Low frequency theta oscillations (4-8 Hz) have been observed to increase in terms of power during working memory tasks [2, 3] and in power and phase-locked discharge of single neurons in a visual memory task . In hippocampus the phase of theta functions as the clock signal for timing of pyramidal neurons and long-term potentiation (theta peaks) and depotentiation (theta troughs) . These findings may reflect the patterns of synaptic plasticity and maintenance of the memory for a stimulus. Fast frequency gamma oscillations (30-70 Hz) can provide tighter control and coordination than lower frequency ones  and are hypothesised to be responsible for higher cognitive functions such as perceptual binding of visual features . Human electroencephalographic (EEG) recordings show event-related gamma activity indicating gamma as a signature of cortical networks underlying object representations. Modulation of oscillatory synchronization can also increase synaptic gain at postsynaptic target sites thereby potentiating responses to learned stimuli [9, 10].
Coupling between gamma amplitude and theta phase (theta-nested gamma) has been reported in both cortex and hippocampus [1, 11–13] and provides an effective combination for neuronal populations to communicate and integrate information during visual processing and learning. It may also provide a process of temporal segmentation that can maintain multiple working memory items . Altered coupling has been reported both in the context of human cognitive and perceptual tasks in the cortex  and in the rat hippocampus during item-context association learning , although how this might act to modulate neuronal activity has yet to be established.
There is still debate as to whether functionally important changes in theta or gamma involve amplitude or phase parameters, or both. Some studies report that theta phase rather than amplitude is correlated with cognitive processes, the so-called phase reset model [1, 15, 16], while others place more importance on coupling between theta amplitude and gamma frequency [11, 12]. The magnitude of both theta and gamma oscillations during encoding also appears to predict the efficacy of subsequent recall  and theta can both modulate gamma amplitude  and the firing of single neurons . The ratio of theta to gamma power has also recently been shown to be correlated with memory function in humans . It is clearly important therefore that changes in different theta and gamma parameters are investigated in a number of different learning contexts to help establish some general principles and also to aid development of neural network models which can further inform our interpretation of the outcome of these changes on neural encoding.
Face recognition learning is known to involve the inferotemporal cortex in humans, monkeys and sheep  and there is phase locking between neuronal activity and theta in this region in humans . We have therefore investigated the effects of face and object discrimination learning on theta and gamma oscillations and coupling and neuronal activity in sheep IT using 64-electrode recording arrays. Our results have identified learning-related changes in the amplitude of theta, the theta/gamma ratio and the coupling between theta phase and gamma amplitude. We have therefore also developed a neural network model which can effectively reproduce our electrophysiological findings. This model predicted that a consequence of these learning evoked changes in theta amplitude and theta-gamma coupling would be a potentiation of the firing of downstream neurons by desynchronizing the firing of excitatory neurons projecting to them. The presence of this predicted desynchronization effect following learning was then confirmed in multiunit activity recordings from the IT.
Visual discrimination performance during recordings
Overall local field potential and MUA data were collected from 51 separate blocks (Sheep A: 17, B: 24, C: 10) of visual discrimination trials (20-60 trials per block). During these trials the sheep were each presented with a total of 4 to 10 different face pairs and in addition two sheep were each presented with a non-face object pair (see Additional file 1, Figure S1). Successful learning was defined as the first block of 20 trials during which the animal achieved > 80% and then continued subsequently to perform at or above this criterion. To compare different electrophysiological parameters as a function of learning blocks of trials for each animal were sub-divided on the basis of whether the > 80% correct learning criterion for a particular face or object pair had been achieved or not. Respective mean ± sem discrimination performances on trial blocks during and after learning were: Sheep A: 57.8 ± 3.5% vs 89.4 ± 3.1%; Sheep B: 58.7 ± 8.5% vs 89.4 ± 1.5%; Sheep C: 66.8 ± 8.5% vs 90.5 ± 2.5%. There were no significant differences in response times made by the three animals during trials with errors as opposed to correct choices. After learning, while there was a slight tendency for response times to be shortened this was not significant in any animal. Respective mean ± sem response times for both correct and error choices on trial blocks during and after learning were: Sheep A: 1.89 ± 0.06s (correct) and 1.87 ± 0.09s (error) vs 1.77 ± 0.06s (correct) and 1.73 ± 0.04 (error); Sheep B: 2.42 ± 0.10 (correct) and 2.35 ± 0.10s (error) vs 2.31 ± 0.09s (correct) and 2.48 ± 0.12s (error); Sheep C: 2.99 ± 0.38s (correct) and 3.44 ± 0.38s (error) vs 2.37 ± 0.11s (correct) and 2.38 ± 0.26s (error) - p > 0.05 in all cases). This also indicates that in general the animals were equally motivated to perform the task during and after learning.
Following training, sheep generally learned to discriminate between new pairs of faces at > 80% correct in 20-80 trials although this was highly variable and learning could occur over time-periods of anything from 5-10 minutes to several days or more. In three cases (two in Sheep A and one in Sheep B) the > 80% criterion for novel face pairs was not reached even after 80-140 trials conducted across several days of recording sessions.
Theta and gamma oscillations in inferotemporal cortex
Effects of learning on theta and gamma oscillations
In addition, coupling between theta and gamma was also strengthened after learning in terms of a greater coherence (~7%) between the two frequencies (left IT: F 1,29 = 9.35, p = 0.005, PT, p < 0.001; right IT: F 1,41 = 14.4, p < 0.001, PT, p = 0.0011)(see Figure 2E). A significant effect of learning was also seen using the mean vector length method  (left IT: mean ± sem change during learning = 1.2 ± 10.1%, after learning = 11.1 ± 3.01% - F 1,29 = 4.37, p = 0.04, PT, p = 0.0189; right IT: mean ± sem change during learning = - 0.93 ± 5.3%, after learning = 22.1 ± 6.64% - F 1,41 = 11.19, p = 0.002, PT, p < 0.001 ) and modulation index method  (left IT: mean ± sem change during learning = 0.6 ± 7.07%, after learning = 26.0 ± 8.04% - F 1,29 = 4.39, p = 0.04, PT, p = 0.003; right IT: mean ± sem change during learning = 0.9 ± 4.67%, after learning = 28.8 ± 8.15% F 1,41 = 9.83, p = 0.003, PT, p < 0.001) for quantifying coupling strength. There was also significantly increased tightening of theta phase across electrodes after learning in the right IT (left IT: F 1,29 = 0.91, p = 0.348; right IT: F 1,41 = 8.09, p = 0.007, PT, p = 0.0014 - see Additional file 1, Figure S8a) where z-scores were ~3-fold higher in the right IT than in the left IT (see Figure 2F). A 3-way ANOVA adding hemisphere as a factor showed that while there was no overall effect of side (F 1,56 = 0.03, p = 0.871) or learning (F 1,56 = 0.24, p = 0.629) there was a small significant interaction between learning and hemisphere (F 1,56 = 3.89, p = 0.05) indicating that learning was affecting phase tightening differentially in the left and right IT.
We found no evidence for extensive theta-phase resetting in response to stimulus presentation with < 1.5% of recording electrodes showing a significant (p < 0.05) effect. There was also no significant increase in the associated phase reset z-score following learning in the left IT (F 1,29 = 1.30, p = 0.263), although in the right IT significance was just achieved with the ANOVA but not the permutation analysis (F 1,41 = 5.00, p = 0.031, PT, p = 0.062). However, these mean z-scores for theta phase reset were generally very low (from 1.5-2.1) (Figure 2G).
Overall, levels of theta synchronization across recording electrodes were higher in the right IT (> 95%) than in the left (~48%) (3-way ANOVA: F 1,56 = 189.7, p < 0.001, PT, p < 0.001) but with no effect of learning in either hemisphere (left IT: F 1,29 = 0.34, p = 0.657; right IT: F 1,41 = 0.55, p = 0.561) (see Figure 2H).
Correlations between altered theta and gamma oscillations and behavior
Correlations between theta and gamma parameters in IT and discrimination performance
Theta and gamma
n = 17
n = 24
n = 22
n = 10
n = 10
n = 34
n = 49
Theta phase tightening
Theta/gamma coupling parameters
Mean vector length
Finally, for the two learned face-pairs where we ran additional blocks of trials in Sheep A and B with the face images inverted, this resulted in a complete inability to discriminate between the two faces (55% and 45% correct responses respectively). As expected, the patterns of theta/gamma changes in both cases were similar to those seen in the face pairs prior to the learning criterion being achieved (data not shown).
Interdependence of theta amplitude and theta-gamma coupling changes
Visual evoked potentials and MUA responses
Following stimulus onset there were similar peak response latencies for the VEP, MUA and peak theta amplitude (overall mean ± sem across all recording sessions in the 3 animals: VEP: P100 = Right 133 ± 7 ms, Left = 118 ± 6 ms; N300 = Right 367 ± 13 ms Left 321 ± 15 ms; MUA: Right 266 ± 4 ms, Left 255 ± 4 ms; Theta = Right 265 ± 14 ms, Left 234 ± 11 ms) (see Figure 2A). However, as would be predicted from the lack of evidence for theta phase resetting in response to stimulus presentation, there was no significant correlation between the latency of peak theta amplitude and that of the MUA (left IT: r = 0.164; right IT: r = -0.02, P > 0.05 in both cases) or the N300 (left IT: r = 0.35; right IT = 0.16, P > 0.05 in both cases) although there was with the P100 in the left (r = 0.49, p = 0.006) but not the right IT (r = 0.27, p = 0.08). The patterns of the VEPs were visibly different between face and object pairs where the magnitude of the P100 component was consistently greater for face pairs than for non-face objects and face inversion also reduced its size compared to when upright faces were presented (see Additional file 1, Figure S4). Neither the response latencies nor the magnitudes of the P100 and N300 components of VEPs were influenced by learning (response latencies: P100 - excluding the two non-face pairs: left IT: during learning = 111.1 ± 7.4 ms vs after learning 123.4 ± 8.2 ms, F 1,27 = 0.383, p = 0.541; right IT: during learning = 131.8 ± 5.9 ms vs after learning 132.7 ± 11.9 ms, F 1,27 = 0.054, p = 0.817, N300: left IT: during learning = 310.2 ± 24.3 ms vs after learning 327.6 ± 19.9 ms, F 1,29 = 0.012, p = 0.914; right IT: during learning = 380.3 ± 16.6 ms vs after learning 355.8 ± 19.9 ms, F 1,41 = 0.178, p = 0.675, response magnitude: P100: left IT: during learning = 5.56 ± 0.09 μV vs after learning = 7.29 ± 1.5 μV, F 1,27 = 0.559, p = 0.461; right IT: during learning = 5.0 ± 0.9 μV vs after learning = 7.3 ± 1.4 μV, F 1,37 = 0.241, p = 0.627, N300: left IT: during learning = -12.52 ± 3.0 μV vs after learning = -15.04 ± 2.37 μV, F 1,29 = 0.328, p = 0.541; right IT: during learning = -11.5 ± 1.4 μV vs after learning = -14.01 ± 1.83 μV, F 1,41 = 0.67, p = 0.418).
A number of MUA recording channels showed significant (p < 0.05) phase locking with theta in each block of visual discrimination trials (during learning mean ± sem % across the theta range = left IT: 20.4 ± 2.4%; right IT 19.7 ± 1.7% and after learning = left IT: 30.7 ± 2.4%; right IT: 25 ± 2.4%). The slight increase after learning was significant in the left (F 1,26 = 4.91, p = 0.036, PT, p = 0.022) but not in the right IT (F 1,38 = 1.95, p = 0.171). Phase locking occurred in each 1 Hz bandwidth from 4-8 Hz frequencies both during and after learning (see Additional file 1, Figure S5).
Theta-nested gamma generated by a neuronal network
Effects of altering NMDA receptor sensitivity in the model
We first established that the model was able to reproduce patterns of theta and gamma activities observed in the IT. Simulations revealed similar changes in theta power during stimulus application and at the same latency (see Additional file 1, Figure S6 and Additional file 1, Figure S7). Theta activity was also strongly synchronised across the network and there was phase tightening during stimulus presentation (see Additional file 1, Figure S7 and Additional file 1, Figure S8).
Having validated the model's utility we next used it to investigate potential functional consequences of altering the ratio of theta to gamma to produce shallow nested gamma (as seen after learning) on communication between excitatory and downstream neurons in comparison with deeper nested gamma (similar to before learning) or where gamma activity was minimal. Figure 7B shows that the downstream neuron response during the stimulus is strongest when there is shallow nested gamma and there is increased theta amplitude and strong coupling between the two frequencies. With deeper nested gamma, excitatory neuron responses appear more highly synchronized and there is reduced theta/gamma coherence and a weaker downstream neuron response. When gamma is minimised to produce a very high ratio of theta to gamma there is reduced excitatory output and downstream neuron activity and theta/gamma coherence (Figure 7B). Thus for optimal coupling between gamma and theta, and to evoke maximal responses in the downstream neuron, gamma should be shallow nested on theta, producing a slightly increased theta-gamma ratio, as seen after learning.
Finally we also used the model to confirm findings in the IT that theta-gamma coupling changes were not necessarily dependent upon those of theta amplitude. To achieve this we systematically varied the strength of the connections between the slow and fast inhibitory neurons (GAsf) and within the slow inhibitory neurons (GAss). Simulations at different stimulus strengths showed that under these circumstances increased coupling strength could occur without any increase in theta amplitude (see Additional file 1, Figure S9).
Temporal desynchronization of neuronal firing
Since our model predicted that learning-induced changes in theta and its relationship to gamma should increase the impact of the firing of excitatory neurons on those downstream by desynchronizing their output we hypothesized that such desynchronization should also occur in MUA recordings from IT neurons. Despite the contribution of both inhibitory interneurons and output neurons to the MUA, Figure 6C shows that after learning there was indeed a significant overall desynchronization change across the 3 animals during the period of the first theta wave after stimulus onset (during learning synchronization index = 0.0871 ± 0.007 pre-stimulus and 0.0937 ± 0.01 during stimulus, t-test, t4 = -1.47, p = 0.22; after learning = 0.0861 ± 0.005 vs 0.0798 ± 0.006, t4 = 5.61, p = 0.005). The proportionate change in the synchronization index between during and after learning was significant (left IT: F 1,27 = 9.71, p = 0.004; right IT: F 1,39 = 15.17, p = 0.001 - Figure 6D). As in the model, levels of synchronization across the two hemispheres were also negatively correlated with the theta/gamma ratio (Pearson correlation, r = -0.32, p < 0.001, Figure 9C) and after learning there were significantly more 5 ms bins with spike activity across each electrode during stimulus period theta waves in the left IT (4.30 ± 0.11 per bin during learning vs 5.96 ± 0.42 per bin after learning, F 1,27 = 5.21, p = 0.031, PT, p = 0.02); using only electrodes showing stimulus-evoked increased firing rates in the recording arrays). This almost achieved significance in the right IT with the ANOVA test but did so with the permutation test (5.19 ± 0.30 per bin during learning vs 5.85 ± 0.27 per bin after learning, F 1,39 = 3.59, p = 0.066, PT, p = 0.04)(see Additional file 1, Figure S10). There was also a significant negative correlation between the number of bins with spike activity per theta wave and the magnitude of the synchronization index across both hemispheres (r = - 0.30, p < 0.001). Finally, actual visual discrimination performance across the three animals was significantly negatively correlated with the change in synchronization during the stimulus (left IT, r = - 0.42, p = 0.018; right IT, r = -0.35, p = 0.018 - see Additional file 1, Figure S10- individual animals A Right IT = -0.35; B left IT = -0.46*, right IT = -0.48*; C left = -0.37, right = -0.77* - *p < 0.05).
Overall, our results provide the first demonstration that both theta amplitude and theta-gamma coupling in IT are strongly and independently influenced by learning and may act to amplify and improve discriminability of inputs converging onto downstream neurons through a temporal desynchronization of neuronal firing. The magnitude of observed changes also correlates strongly with behavioral discrimination performance. The novel neural network model we have developed demonstrates that competitive and reciprocal coupling between fast and slow inhibitory interneurons and excitatory output neurons is important for generating theta-nested gamma and that learning-evoked changes in IT can be simulated by increasing NMDA receptor sensitivity, as in many other learning situations .
The time-course required for changes in theta and gamma correlates precisely with achievement of the > 80% learning criterion and since this time course was highly variable it effectively rules out any simple explanation of our findings in terms of elapsed time or stimulus repetition. Indeed, no changes were observed when animals failed to learn even after large numbers of trials. It is unlikely that there was differential attention to the visual stimuli presented following learning since animals were extensively trained to view the fixation stimulus before each trial was initiated, and latencies and magnitudes of VEPs as well as behavioral response times were unchanged. Sheep also have very limited eye movements and the receptive fields of IT neurons responding to faces and other visual stimuli are very large . While we cannot completely rule out some contribution of the P100 or N300 components of the VEP to increased theta amplitude during stimulus presentations the fact that neither of their latencies not their magnitudes were significantly influenced by learning make it unlikely that they were contributing to the large increase in theta amplitude found. Indeed, even the small number of recording electrodes showing increases in theta amplitude seen during stimulus presentation before the learning criterion had been reached did not achieve significance across the three animals.
To the best of our knowledge this is the first demonstration of extensive theta-nested gamma in IT and we confirmed both the presence of coupling and learning-evoked changes as well as correlations with behavioral performance using a coherence method as well as mean vector length  and modulation index  methods. Our model shows that it can be generated by a simple network of excitatory glutamatergic pyramidal neurons and slow and fast inhibitory GABAA receptor interneurons that are likely to be present in this as in other brain regions [23, 24]. In the model each module of one excitatory neuron and its two inhibitory neurons is capable of generating theta-nested gamma, so our findings are not dependent upon model size. That both slow and fast inhibitory interneurons are required for the generation of theta-nested gamma confirms a previous prediction , although differs from another study in the hippocampus suggesting that h-currents generated in oriens-lacunosum interneurons are important .
Although theta amplitude changes can enhance the strength of theta-gamma coupling by making theta phase more detectable [11, 12] we have shown that learning effects on theta amplitude and theta gamma coupling in IT can be dissociated in terms of which theta frequencies are involved. Thus, whereas behaviorally correlated changes in theta amplitude occurred across the 4-8 Hz range, those in coupling with gamma only occurred at 6 Hz (right IT) or 8 Hz (left IT). There were also no positive correlations between changes in theta amplitude and in theta-gamma coherence. Interestingly, in comparison with lower theta frequencies (4-6 Hz), higher ones (6-12 Hz) are resistant to cholinergic drugs . This might suggest a non-cholinergic mechanism for theta-nested gamma. In any event, learning evoked changes in theta amplitude and theta-gamma coupling in the IT clearly exhibit a degree of independence and may have separate functional significance. This is supported by a recent study reporting increased theta-gamma coupling in the hippocampus following item-context association learning in rats, but without changes in theta amplitude .
Learning evoked increases in the theta-gamma ratio in the present study were strongly correlated with both behavioural performance and also with desynchronization in both our network model and the IT. This appears contradictory with recent human findings that an increased theta-gamma ratio is found in individuals with significant cognitive impairment . However, our network model shows that when large increases in this ratio are generated (by reducing gamma power), similar to those reported in humans with cognitive impairment, then this results in reduced theta amplitude changes, weaker coupling between theta and gamma and also reduced firing in both excitatory output and downstream neurons in response to stimuli. It therefore seems probable that while small increases in the theta-gamma ratio are indicative of successful learning, very large increases are more likely to reflect learning impairment.
We have found no evidence to support a key role for resetting of theta phase in cortical information processing contrary to some previous studies [1, 15, 16]. The virtual absence of phase resetting supports the notion that event-related potentials are not generated by the phase resetting of oscillatory activity as has been proposed . A previous study has also failed to find evidence for phase-resetting by visual stimuli . An obvious role for phase resetting by an external stimulus is to synchronize oscillations across a wide network; however in our IT recordings we always found very high levels of synchronization in ongoing theta activity, so there would have been no real advantage in this respect if learning had resulted in stronger phase resetting.
The lack of any alteration in the amplitude of gamma oscillations following learning is in-line with our failure to find any change in stimulus-evoked neuronal firing rates in the IT following learning. Indeed, previous research in monkeys has also failed to find evidence for reward-associated learning changes in firing rates of individual IT neurons  and there is evidence in both sheep  and monkey  IT for population-based encoding. There was also no evidence for learning effects on latencies or magnitudes of the P100 and N300 components of the VEP in IT so VEP changes are unlikely to have contributed to any learning effects on phase locking . There was some indication that the magnitude of the P100 was increased in response to upright face pairs compared to inverted face pairs or upright non-face objects. This latter observation is similar to the face sensitive P100 in humans [32, 33] and the N170 component in monkeys  and humans . Some studies have reported increased magnitude of various VEP components as a function of face familiarity [34, 36, 37], and in response to visual tetanic stimulation , although not in the context of learning visual discrimination between pairs of objects used in the current study.
No differential learning effects were observed between the two brain hemispheres in the MUA and LFP parameters measured in current study. However, we have also analysed our LFP data using novel causality algorithm approaches and shown that learning reduces the strength of causal connections from the left to the right IT and increases the frequency of causal connections within the right IT . This increased connectivity may partly explain our observation here that only the right IT shows significant strengthening of the regularity of theta phase following learning. It is notable in this context too that in general the right IT showed significantly higher levels of theta synchronization than the left.
The presence of phase-locking between theta and MUA in the IT is in agreement with a recent study in humans  and phase locking has been reported in other neocortical areas in monkeys . There was also a tendency for this to increase in IT following learning and a recent study has reported increased phase locking of neuronal activity with theta in human temporal lobe neurons after learning . The presence of extensive phase locking between theta and neuronal spiking activity in the IT provides a potential mechanism whereby altered theta activity following learning could modulate neuronal firing to produce the desynchronization effects we have observed.
Although increased theta amplitude and theta gamma coupling might have been expected to lead to tighter control and greater synchronization in neuronal firing our model shows this is not the case. Indeed, a recent study on visual cortex neurons has also reported visual stimulus-evoked decorrelation . The large changes in theta activity observed would act to modulate the firing thresholds of the output neurons cells across the network. This would inevitably increase the variability in theta activity across the network leading to a wider range of firing thresholds and neurons being less likely to fire synchronously. This is supported by our model and IT findings showing that neural spike activity is more spread out in time during theta waves in the stimulus period following learning.
The prediction from our model that increased desynchronization in IT output neurons would enhance responses by downstream target neurons might also seem counter-intuitive. However, the more synchronized are the outputs from excitatory neurons converging onto a downstream neuron, the more information can potentially be lost as more excitatory post synaptic potentials (EPSPs) are generated than are necessary to cause the downstream neuron to fire. Where EPSPs generated are more separated in time they are less likely to be rendered impotent by refractory period limitations and contribute more efficiently towards eliciting responses by the downstream neuron. Spreading the temporal pattern of inputs reaching the downstream neuron would also enhance the information content of the inputs it is required to decode
It would clearly be difficult to test the above model prediction directly in vivo without being able to make simultaneous recordings from multiple connected neurons in say IT and the frontal cortex. However, our combined in vivo and model simulation findings do provide a mechanism for how learning induced changes in theta amplitude and theta-nested gamma could modulate temporal aspects of neuronal firing in neocortical networks such that downstream networks exhibit potentiated responses even in the absence of altered firing frequencies arriving from neuronal inputs.
Face discrimination learning produces significant increases in the magnitude of theta amplitude and the theta-gamma ratio and the strength of theta-gamma coupling in the inferotemporal cortex. Importantly these changes are all significantly correlated with actual behavioral performance and theta amplitude and theta-gamma coupling changes appear to occur independently. Interestingly, learning did not produce significant changes in IT neuronal firing frequencies although neuronal firing was often coupled to theta-phase. The network model which we have developed to reproduce these changes suggests that a key function of such learning-evoked changes in theta amplitude, the theta-gamma ratio and theta-gamma coupling may be to increase temporal desynchronization in neuronal firing leading to optimal timing of inputs to downstream neural networks and thereby potentiating their responses. The model's efficacy was supported by the fact that this temporal desynchronization was confirmed in our IT recordings following learning. An important functional consequence of the learning evoked changes in theta and gamma we have found in may therefore be to potentiate responses by neurons in IT projection regions to learned visual stimuli through a slight temporal desynchronization of firing by IT output neurons.
Animals and visual discrimination training
Three female sheep were used (Ovis aries, one Clun Forest and two Dorsets). The animals were trained initially over several months to perform operant-based face (sheep) or non-face (objects) discrimination tasks with a choice being made between two simultaneously presented pictures (side by side) only one of which was associated with a food reward. The position (left or right) of the rewarded picture was randomised in each trial. During stimulus presentations animals stood in a holding trolley and indicated their choice of picture by pressing one of two touch panels located in the front of the trolley with their nose. For correct responses the food reward was delivered automatically to a hopper between the two panels. The life-sized pictures were back projected onto a screen 0.5 m in front of the animal using a computer data projector. A white fixation spot on a black background was presented constantly in between trials to maintain attention and experimenters waited until the animals viewed this spot before triggering presentation of the image pairs (since sheep don't have extensive eye movements monitoring head position is generally sufficient to establish gaze direction). The stimulus images remained in view until the animal made an operant response (generally around 1-3 s). In each case successful learning of a face or object pair required that a performance criterion of > 80% correct choice over blocks of 20-40 presentation trials was achieved consistently (Chi-square, p = 0.05 for 20 trials and 0.01 for 40 trials). By the end of training animals were normally able to reach the > 80% correct criterion after 40-80 learning trials and to maintain this performance although in a few cases learning did not occur even after 80-140 trials. Some previously learned stimulus pairs (over periods ranging from 10 days to 9 months) were also presented during subsequent electrophysiological recording experiments although the animals were mainly presented with novel stimulus pairs and neurophysiological parameters recorded before and after the learning criterion was achieved.
For each sheep recordings were made in response to up to 10 different face or non face object pairs (Sheep A: 5 novel face and 1 novel object pair; B 7 novel face pairs, 3 previously learned face pairs and one previously learned object pair; Sheep C: 2 novel face pairs and 2 previously learned face pairs. In addition following learning effects of image inversion were recorded for one face pair in both Sheep A and Sheep B to assess whether any components of the VEP were sensitive to upright faces. Learning effects were monitored over between 80-189 trials and data was collected over blocks of 20-40 trials. For the face pairs Sheep A and B were discriminating between the faces of different socially familiar or unfamiliar individuals (face identity discrimination) whereas for Sheep C discrimination was between calm and stressed face expressions in the same animal (n = 3 pairs) or in different animals (n = 1 pair). With this latter animal the calm face was the rewarded stimulus. Where novel face or object pairs were being learned during recordings the > 80% performance criterion was normally achieved in 20-80 training trials. The face and object pairs used for each of the 3 sheep are shown in Additional file 1, Figure S1.
All animal experiments were performed in strict accordance with the UK 1986 Animals Scientific Procedures Act (including approval by the Babraham Institute Animal Welfare and Ethics Committee) and during them the animals were housed inside in individual pens and able to see and communicate with each other. Food and water were available ad libitum. Post-surgery all animals received both post-operative analgesia treatment to minimise discomfort and antibiotic treatment to prevent any possibility of infection.
Electrophysiological recordings and analyses of local field potentials (LFPs) and multiunit activity (MUA) in IT
Following initial behavioral training sheep were surgically implanted under general anesthesia (fluothane) and full aseptic conditions with either unilateral (one animal) or bilateral planar 64-electrode (for configuration see Additional file 1, Figure S6A) arrays (epoxylite coated, etched, tungsten wires with 250 μm spacing - total array area ~2 mm × 2 mm, electrode impedance ~0.2 MΩ, tip diameter ~ 1 μm, tip exposure length ~100 um) aimed at the IT. The electrode lengths varied by ~1 mm and so this combined with the tip exposure electrodes would have been recording activity across all cortical layers. Holes (0.7 cm diameter) were trephined in the skull and the dura beneath cut and reflected. Electrode arrays were placed 18-20 mm lateral to the midline, 35 mm posterior to the tip of the frontal pole and at a depth of 20-22 mm from the brain surface using a stereotaxic micromanipulator. Electrode depths and placements were calculated with reference to X-rays, as previously described . They were fixed in place with dental acrylic and stainless-steel screws attached to the skull. Two of these screws acted as reference electrodes, one for each array. Electrodes were connected to 34 pin female plugs (2 per array) also cemented in place on top of the skull.
LFP and MUA recording protocols
Starting 3 weeks after surgery the electrodes were connected via male plugs and ribbon cables to a 128 channel electrophysiological recording system (Cerebus 128 Data Acquisition System - Blackrock Microsystems, USA) and recordings made during performance of the different face and non-face pair operant discrimination tasks. This system allowed simultaneous recordings of both neuronal spike and local event-related (LFP) activity from each electrode. Typically, individual recording sessions lasted around 30 min and animals were presented with 2-6 blocks of 20-40 trials. There was at least a week between individual recording sessions in each animal.
The LFPs were sampled at 2 kHz and MUA spikes at 30 kHz (bandpass 0.3 Hz - 7.5KHz) and digitized for storage from ~3 seconds prior to the stimulus onset to ~3 seconds after the stimulus onset (stimulus durations were generally 1-3 s). Neural recordings from our data acquisition system consisted of two distinct large raw data files, one for the LFP and the other for the MUA. We used custom Spike 2 (Cambridge Electronic Design, Cambridge, UK) scripts to translate these into text files arranged either by trial or electrode prior to further analysis.
LFP data contaminated with noise artefacts, such as from animal chewing food, were excluded as were LFPs with unexpectedly high power. For LFPs, data were analyzed during a period of 1 second before and 1 second after stimulus onset. Trend was removed before spectral analysis. Any trials having more than 5 points outside the mean ± 5 standard deviation range were discarded before the analysis. The LFPs and MUA responses were all aligned to the onset of the visual stimuli. All analyses were carried out using custom written routines in Matlab (The Mathworks Inc, Natick, MA). Use of custom spike-sorting software revealed that 1-4 single neurons were contributing to the MUA at each electrode .
At the end of the experiments animals were euthanized with an intravenous injection of sodium pentobarbitone and the brains removed for subsequent histological confirmation of X-rays that array placements were within the IT cortex region. The general region where electrodes were located within the IT of these animals is also shown in .
Time dependent spectrum analysis
which has the form of a Gaussian function centred at f0, where f0 determines the wave numbers within the envelope. Here f0 = 0.849 (around 5.3 wave numbers) and this gives a real part where the peaks next to the central peak are half its amplitude.
The wavelet transform was applied to each individual LFP trial at each electrode (for a period of 1 s either side of stimulus onset) and a final time-dependent spectrum estimated as the trial-averaged scalograms (modulus square of the wavelet transform). When comparing pre- and during stimulus theta band activity we used the amplitude of the wavelet transform at 4-8 Hz and averaged it across this band. For the gamma band, amplitude in the 30-70 Hz frequency range was analysed. Theta and gamma power were also calculated, although in our freely behaving animals we found the amplitude measure to be less variable across trials and sessions whereas power was susceptible to abrupt changes making comparisons more difficult. To determine the significance of responsiveness at each electrode a post-hoc t-test (with Benjamini-Hochberg correction) was used to compare the amplitude of the wavelet transform in the 500 ms pre- with the 500 ms during-stimulus periods. Amplitude changes evoked by learning were normalized by subtracting the amplitude value of the pre-stimulus period and dividing by the maximal value for each electrode. The theta/gamma ratio was also calculated as the direct ratio between the theta amplitude and gamma amplitude (or theta power and gamma power).
Cross-frequency coupling between theta and gamma
where α is the significant level (e.g. α = 0.01) and K is the trial number which corresponds to the disjointed number of periodograms. The phase-locking index was then measured by the coherence in the range between 0 and 1, where values close to 1 indicate a strong cross-frequency modulation. A coherence calculation was obtained at all the pair-wise frequency combinations between the two bands and a Bonferroni correction applied to the multiple comparisons over all the frequency pairs.
We confirmed that this coherence method accurately determined the extent to which gamma amplitude changes were locked to theta phase using artificial data. We generated one theta wave and gamma wave and nested (added) them together using two sine waves of 5 Hz and 50 Hz which were linearly mixed. The gamma frequency sine wave (50 Hz) had an amplitude 1/5th of the theta frequency wave (5 Hz). Using a trial length of 500 ms, 30 trials were generated with a sampling frequency of 1 kHz. White noise was then added to the mixed sine waves with a signal to noise ratio equal to -5 dB. The coherence between the theta phase and gamma amplitude was maximal when gamma was nested directly on top of theta (Additional file 1, Figure S3) confirming that our coherence measure reliably measures the strength of coupling between theta phase and gamma amplitude. However in view of recent debates over the relative merits of different methods for measuring coupling between theta phase and gamma amplitude in electrophysiological data we additionally used the mean vector length  and the modulation index  approaches for our IT data and these produced similar results to those obtained with our coherence method.
For IT recordings, theta-gamma coupling analysis was performed for all the electrodes and at each electrode the theta/gamma values were calculated for all the pairs in the theta and gamma band. In the majority of cases 40 trials were analysed in each session. Where more than 40 trials were recorded only the first 40 were analysed and where less than 40 trials occurred the coherence was normalised to match that for 40 trials to avoid any bias due to differential numbers.
Theta phase reset
Therefore the hypothesis of uniformity could be rejected at a certain significance level if phase-locking was found for that electrode. The Z-score for the Rayleigh statistic is given as Z = NR2.
To determine if theta band waves exhibited phase-resetting with a locked phase over trials, we calculated the Z -scores as a function of time in the during-stimulus range across all trials. An electrode recording was considered to be exhibiting phase-locking if all the samples from the time of stimulus onset at a given frequency (4 - 8 Hz) passed the criterion of the Rayleigh test (p < 0.01) across two full oscillatory cycles. A comparison was made across all the electrodes in a recording array and a Bonferroni correction applied to compensate for type-I errors.
Theta phase synchronization
To assess whether there was synchronization of the LFP phases, the Rayleigh statistic was also used to calculate a Z-score across all electrodes. In each trial the Z-score for theta phase (4-8 Hz) was calculated for each time point in the 500 ms pre- and the 500 ms during-stimulus periods. If > 80% of all time points across the entire 1 s period showed significant phase-locking (p < 0.05), then LFPs were considered to be synchronized in that trial. For bilateral recordings, the left and right hemispheres were analysed separately.
We calculated the Z-scores for LFP phases in 500 ms pre-stimulus and 500 ms during stimulus periods across all the electrodes in the recording array. If the Z-score was significantly higher (t test, p < 0.05) in the during-stimulus period than in the pre-stimulus period then the phase was considered to be tightened. We used the percent change from the pre- to during-stimulus period to measure altered phase tightening.
Visual evoked potential (VEP)
The VEP was extracted from the LFPs by trial-averaging after aligning the data to stimulus onset. Two major peaks were identified from the VEP in the initial 500 ms of stimulus presentation: a positive peak at ~100 ms (P100) and a negative peak at ~300 ms (N300). We calculated the latency for these two peaks by finding the time corresponding to the maximum and minimum peak value respectively. The amplitudes of these two peaks were calculated as their peak values after subtracting the average baseline in the 100 ms before stimulus onset.
MUA phase-locking with theta and firing rates
For the analysis of MUA data a Gaussian kernel with width of 30 ms was convolved to the neuronal spike train. We used the maximum peak value in the initial 500 ms of stimulus presentation to characterise MUA response latency. To calculate phase-locking of MUA on each channel to theta, the same LFP phase data at theta band were used as in the analysis of theta phase-reset. At each given frequency from 4-8 Hz, (at 1 Hz intervals) a total length of 2 s (1 s before stimulus and 1 s after stimulus) recording data was used for both LFP and MUA. LFP phases corresponding to the MUA spike time data across all the trials were then tested for analysed for whether they passed the criterion of the Rayleigh test (p < 0.05). Where they did so, this showed the distribution of phase values was not uniform at that particular electrode and frequency and that there was significant phase-locking between MUA and the theta band wave.
For firing rate changes associated with a stimulus and overall calculation was made for each channel in each block of trials (1s before vs 1 s during a stimulus). Recording channels were also sub-divided into those showing an overall increase (excited) or decrease (inhibited) in firing rate during the stimulus.
MUA synchronization index
We measured synchronization in MUA data by counting the spikes within a brief time window (bin width, 2.5-10 ms). This is similar to a peri-stimulus time histogram (PSTH) over all MUA channels. In each trial we produced a PSTH over all the MUA channels and normalized it by the sum of the counts in all PSTH bins. If synchronization occurred in a certain time bin there should be a high spike count for that bin. Normalization was carried out to ensure that the influence of differential firing rates was removed. We then defined a MUA synchronization index as the sum of all the normalized spike counts which exceed half of the maximum value. We calculated the synchronization index choosing a bin width of 5 ms although we also used bin sizes of 2.5 ms and 10 ms and similar trends were observed.
These were performed at both the population and individual animal level. For each parameter measured 2-way ANOVAs (animals and learning as factors) were first carried out for the 3 animals where recordings were made in the right IT and the 2 animals in the left IT. To test for hemisphere differences 3-way ANOVAs (animals, learning and hemisphere as factors) were conducted. In all cases where data for a particular parameter failed the Shapiro-Wilks normality test an additional permutation analysis was carried out using 1000 replacements. Additionally, analyses within animals were carried out using t-tests corrected for multiple comparisons. Correlations between electrophysiological parameters and behavioral performance were carried out using Pearson tests.
For the fast GABAA channel, we chose τ I,f = 9 ms and α I,f = 1 ms-1. For the slow GABAA channel, τ I,s = 50 ms and α I,s = 0.2 ms-1. In the simulation, spikes in all presynaptic neurons are connected to a convergent neuron. For the background current of EX cells, we set where ξ i (t) is white noise with variance σ e = 0.01. For inhibitory cells, we set the background currents fixed and homogenous, for INf cells, and for INs cells.
Parameters and analytical methods used in application of the model
To generate post learning effects the following coefficient values were used for the different sites of AMPA (AM), NMDA (NM) and GABAA (GA) receptors (e = excitatory neuron, s = slow inhibitory neuron and f = fast inhibitory neuron): AMee = 0.02; AMef = 0.08; AMes = 0.0005; NMee = 0.0035; NMef = 0.001; NMes = 0.00055; GAff & GAss = 0.08; GAfe = 0.015; GAse = 0.06; GAsf = 0.03. For pre-learning only the values of two NMDA receptor coefficients were reduced: NMee to 0.002 and NMes to 0.0001. For parameters that generated different theta/gamma patterns (Figure 7), the shallow nested gamma was generated using the post learning parameters above. The parameters for deep nested gamma were the same as the shallow nested case except that GAfe = 0.045. The minimal gamma used the same parameters as the shallow nested gamma except that GAse and GAsf = 0.12. All the methods for calculating theta/gamma parameters were the same as for the data from IT recordings.
Research was supported by a Biotechnology and Biological Sciences (BBSRC) project grant (BBS/B/07691) and Search Foundation grant to KMK. JF was supported by Engineering and Physical Sciences research Council (EPSRC) and European Union grants. XZ was supported by The National Natural Science Foundation of China (Grant Nos. 10771155, 10971196). We thank Andrew Tate and Andrea Leigh for their help with carrying out behavioral and recording experiments and Dr Anne Segonds-Pichon for professional statistical advice on analysis of the IT data.
- Buzsaki G: Rhythms of the Brain. 2006, Oxford University Press, USAView Article
- Gevins A, Smith ME, McEvoy L, Yu D: High-resolution EEG mapping of cortical activation related to working memory: effects of task difficulty, type of processing and practice. Cereb Cortex. 1997, 7: 374-385. 10.1093/cercor/7.4.374.View ArticlePubMed
- Raghavachari S, Kahana MJ, Rizzuto DS, Caplan JB, Kirschen MP, Bourgeois B, Madsen R, Lisman JE: Gating of human theta oscillations by a working memory task. J Neurosci. 2001, 21: 3175-3183.PubMed
- Lee H, Simpson GV, Logothetis NK, Rainer G: Phase locking of single neuron activity to theta oscillations during working memory in monkey extrastriate visual cortex. Neuron. 2005, 45: 147-156. 10.1016/j.neuron.2004.12.025.View ArticlePubMed
- Holscher C, Anwyl R, Rowan MJ: Stimulation on the positive phase of hippocampal theta rhythm induces long-term potentiation that can be depotentiated by stimulation on the negative phase in area CA1 in vivo. J Neurosci. 1997, 17: 6470-6477.PubMed
- Jensen O, Kaiser J, Lachaux JP: Human gamma-frequency oscillations associated with attention and memory. Trends Neurosci. 2007, 30: 317-324. 10.1016/j.tins.2007.05.001.View ArticlePubMed
- Csibra G, Davis G, Spratling MW, Johnson MH: Gamma oscillations and object processing in the infant brain. Science. 2000, 290: 1582-1585.View ArticlePubMed
- Busch NA, Herrmann CS, Muller MM, Lenz D, Gruber T: A cross-laboratory study of event-related gamma activity in a standard object recognition paradigm. Neuroimage. 2006, 33: 1169-1177. 10.1016/j.neuroimage.2006.07.034.View ArticlePubMed
- Fries P, Reynolds JH, Rorie AE, Desimone R: Modulation of oscillatory neuronal synchronization by selective visual attention. Science. 2001, 291: 1560-1563. 10.1126/science.1055465.View ArticlePubMed
- Fries P: A mechanism for cognitive dynamics: neuronal communication through neuronal coherence. Trends Cog Sci. 2005, 9: 474-480. 10.1016/j.tics.2005.08.011.View Article
- Canolty RT, Edwards E, Dalal SS, Soltani M, Nagarajan SS, Kirsch HE, Berger MS, Barbaro NM, Knight RT: High gamma power is phase-locked to theta oscillations in human neocortex. Science. 2006, 313: 1626-1628. 10.1126/science.1128115.PubMed CentralView ArticlePubMed
- Tort ABL, Komorowski RW, Manns JR, Kopell NJ, Eichenbaum H: Theta-gamma coupling increases during the learning of item-context associations. Proc Natl Acad Sci USA. 2009, 106: 20942-20947. 10.1073/pnas.0911331106.PubMed CentralView ArticlePubMed
- Tort ABL, Komorowski R, Eichenbaum H, Kopell N: Measuring Phase-Amplitude Coupling Between Neuronal Oscillations of Different Frequencies. J Neurophysiol. 2010, 104: 1195-1210. 10.1152/jn.00106.2010.PubMed CentralView ArticlePubMed
- Jensen O: Maintenance of multiple working memory items by temporal segmentation. Neuroscience. 2006, 139: 237-249. 10.1016/j.neuroscience.2005.06.004.View ArticlePubMed
- Givens B: Stimulus-evoked resetting of the dentate theta rhythm: relation to working memory. NeuroReport. 1996, 8: 159-163. 10.1097/00001756-199612200-00032.View ArticlePubMed
- Rizzuto DS, Madsen JR, Bromfield EB, Schulze-Bonhage A, Seelig D, Aschenbrenner-Scheibe R, Kahana MJ: Reset of human neocortical oscillations during a working memory task. Proc Natl Acad Sci USA. 2003, 100: 7931-7936. 10.1073/pnas.0732061100.PubMed CentralView ArticlePubMed
- Sederberg PB, Kahana MJ, Howard MW, Donner EJ, Madsen JR: Theta and gamma oscillations during encoding predict subsequent recall. J Neurosci. 2003, 23: 10809-10814.PubMed
- Lakatos P, Shah AS, Knuth KH, Ulbert I, Karmos G, Schroeder CE: An oscillatory hierarchy controlling neuronal excitability and stimulus processing in the auditory cortex. J Neurophys. 2005, 94: 1904-1911. 10.1152/jn.00263.2005.View Article
- Moretti DV, Fracassi C, Pievani M, Geroldi C, Zanetti O, Sosta K, Rossini PM, Frisoni GB: Increase of theta/gamma ratio is associated with memory impairment. Clin Neurophysiol. 2009, 120: 295-303. 10.1016/j.clinph.2008.11.012.View ArticlePubMed
- Tate AJ, Fischer H, Leigh AE, Kendrick KM: Behavioural and neurophysiological evidence for face identity and face emotion processing in animals. Phil Trans Roy Soc B: Biol Sci. 2006, 361: 2155-2172. 10.1098/rstb.2006.1937.View Article
- Jacobs J, Kahana MJ, Ekstrom AD, Fried I: Brain oscillations control timing of single-neuron activity in humans. J Neurosci. 2007, 27: 3839-3844. 10.1523/JNEUROSCI.4636-06.2007.View ArticlePubMed
- Liu L, Wong TP, Pozza MF, Lingenhoehl K, Wang Y, Sheng M, Auberson YP, Wang YT: Role of NMDA receptor subtypes in governing the direction of hippocampal synaptic plasticity. Science. 2004, 304: 1021-1024. 10.1126/science.1096615.View ArticlePubMed
- Sceniak M, MacIver MB: Slow GABAA mediated synaptic transmission in rat visual cortex. BMC Neurosci. 2008, 9: 8-10.1186/1471-2202-9-8.PubMed CentralView ArticlePubMed
- Pearce RA: Physiological evidence for two distinct GABAA responses in rat hippocampus. Neuron. 1993, 10: 189-200. 10.1016/0896-6273(93)90310-N.View ArticlePubMed
- White JA, Banks MI, Pearce RA, Kopell NJ: Networks of interneurons with fast and slow γ-aminobutyric acid type A (GABAA) kinetics provide substrate for mixed gamma-theta rhythm. Proc Natl Acad Sci USA. 2000, 97: 8128-8133. 10.1073/pnas.100124097.PubMed CentralView ArticlePubMed
- Tort ABL, Rotstein HG, Dugladze T, Gloveli T, Kopell NJ: On the formation of gamma-coherent cell assemblies by oriens lacunosum-moleculare interneurons in the hippocampus. Proc Natl Acad Sci USA. 2007, 104: 3490-13495. 10.1073/pnas.0610163104.View Article
- Bland BH, Oddie SD: Theta band oscillation and synchrony in the hippocampal formation and associated structures: the case for its role in sensorimotor integration. Behav Brain Res. 2001, 127: 119-36. 10.1016/S0166-4328(01)00358-8.View ArticlePubMed
- Makeig S, Westerfield M, Jung TP, Enghoff S, Townsend J, Courchesne E, Sejnowski TJ: Dynamic brain sources of visual evoked responses. Science. 2002, 295: 690-694. 10.1126/science.1066168.View ArticlePubMed
- Mazaheri A, Jensen O: Posterior α activity is not phase-reset by visual stimuli. Proc Natl Acad Sci USA. 2006, 103: 2948-2952. 10.1073/pnas.0505785103.PubMed CentralView ArticlePubMed
- Rolls ET, Judge SJ, Sanghera MK: Activity of neurones in the inferotemporal cortex of the alert monkey. Brain Res. 1977, 130: 229-238. 10.1016/0006-8993(77)90272-4.View ArticlePubMed
- Tsao DY, Freiwald WA, Tootel RBH, Livingstone MS: A cortical region consisting entirely of face-selective cells. Science. 2006, 311: 670-674. 10.1126/science.1119983.PubMed CentralView ArticlePubMed
- Hermann MJ, Ehlis A-C, Ellgring H, Fallgatter AJ: Early stages (P100) of face perception in humans as measured with event-related potentials (ERPs). J Neural Transm. 2005, 112: 1073-1081. 10.1007/s00702-004-0250-8.View Article
- Liu J, Harris A, Kanwisher N: Stages of processing in face perception: an MEG study. Nat Neurosci. 2002, 5: 910-916. 10.1038/nn909.View ArticlePubMed
- Peissig JJ, Singer J, Kawasaki K, Sheinberg DL: Effects of long-term object familiarity on event-related potentials in the monkey. Cereb Cortex. 2007, 17: 1323-1334.View ArticlePubMed
- Bentin S, Allison T, Puce A, Perez E, McCarthy G: Electrophysiological studies of face perception in humans. J Cog Neurosci. 1996, 6: 551-565.View Article
- Scott LS, Tanaka JW, Sheinberg DL, Curran T: A re-evaluation of the electrophysiological correlates of expert object processing. J Cogn Neurosci. 2006, 18: 1453-1465. 10.1162/jocn.2006.18.9.1453.View ArticlePubMed
- Anderson B, Mruczek REB, Kawasaki K, Sheinberg D: Effects of familiarity of neural activity in monkey inferior temporal lobe. Cereb Cortex. 2008, 18: 2540-2552. 10.1093/cercor/bhn015.PubMed CentralView ArticlePubMed
- Teyler TJ, Hamm JP, Clapp WC, Johnson BW, Corballis MC, Kirk IJ: Long-term potentiation of human visual evoked potentials. Eur J Neurosci. 2005, 21: 2045-2050. 10.1111/j.1460-9568.2005.04007.x.PubMed CentralView ArticlePubMed
- Ge T, Kendrick KM, Feng J: A novel extended Granger causal model approach demonstrates brain hemispheric differences during face recognition learning. PLoS Comput Biol. 2009, 5 (11): e1000570-10.1371/journal.pcbi.1000570.PubMed CentralView ArticlePubMed
- Rutishauser U, Ross IB, Mamelak AN, Schuman EM: Human memory strength is predicted by theta-frequency phase-locking of single neurons. Nature. 2010, 464: 903-907. 10.1038/nature08860.View ArticlePubMed
- Smith MA, Kohn A: Spatial and temporal scales of neuronal correlation in primary visual cortex. J Neurosci. 2008, 28: 12591-12603. 10.1523/JNEUROSCI.2929-08.2008.PubMed CentralView ArticlePubMed
- Kendrick KM, Baldwin BA: Single unit recording in the conscious sheep. Methods in Neuroscience. Edited by: Conn PM. 1991, New York: Academic Press, 3-15.
- Horton PM, Nicol AU, Kendrick KM, Feng JF: Spike sorting based upon machine learning algorithms (SOMA). J Neurosci Meth. 2007, 160: 52-68. 10.1016/j.jneumeth.2006.08.013. 36View Article
- Torrence C, Compo GP: A Practical Guide to Wavelet Analysis. Bull Am Meteorological Soc. 1998, 79: 61-78. 10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2.View Article
- Fisher NI: Statistical Analysis of Circular Data. 1995, Cambridge University Press
- Halliday DM, Rosenberg JR, Amjad AM, Breeze P, Conway BA, Farmer SF: A framework for the analysis of mixed time series/point process data--Theory and application to the study of physiological tremor, single motor unit discharges and electromyograms. Prog Biophys Mol Biol. 1995, 64: 237-278. 10.1016/S0079-6107(96)00009-0.View ArticlePubMed
- Wang X-J: Synaptic basis of cortical persistent activity: the importance of NMDA receptors to working memory. J Neurosci. 1999, 19: 9587-9603.PubMed
- Wang X-J, Rinzel J: Alternating and synchronous rhythms in reciprocally inhibitory model neurons. Neural Comp. 1992, 4: 84-97. 10.1162/neco.19220.127.116.11.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.