Research Article  Open Access
Decoding Pigeon Behavior Outcomes Using Functional Connections among Local Field Potentials
Abstract
Recent studies indicate that the local field potential (LFP) carries information about an animal’s behavior, but issues regarding whether there are any relationships between the LFP functional networks and behavior tasks as well as whether it is possible to employ LFP network features to decode the behavioral outcome in a single trial remain unresolved. In this study, we developed a networkbased method to decode the behavioral outcomes in pigeons by using the functional connectivity strength values among LFPs recorded from the nidopallium caudolaterale (NCL). In our method, the functional connectivity strengths were first computed based on the synchronization likelihood. Second, the strength values were unwrapped into row vectors and their dimensions were then reduced by principal component analysis. Finally, the behavioral outcomes in single trials were decoded using leaveoneout combined with the nearest neighbor method. The results showed that the LFP functional network based on the gammaband was related to the goaldirected behavior of pigeons. Moreover, the accuracy of the network features (74 8%) was significantly higher than that of the power features (61 12%). The proposed method provides a powerful tool for decoding animal behavior outcomes using a neural functional network.
1. Introduction
The local field potential (LFP) is a lowpass filtered signal and it is considered to primarily reflect the sum of the slow local synaptic currents originating from an area around the electrode tip [1]. In recent years, the LFP has received increasing attention for several reasons. First, the LFP has been shown to reflect sensory and motorrelated signals that can be modulated by cognitive processes, and it provides additional information regarding single neuron activity [2]. Second, the LFP appears to be correlated more closely with the BOLD signal measured by fMRI than spike activity [3–5]. Third, the LFP has been shown to correspond to spike activity under certain behavioral or perceptual conditions [6]. Finally, the LFP is easy to record, even over long periods of time, so it may be an efficient candidate signal for the control of neural prostheses [7].
Recently, an increasing number of studies have employed the LFP to decode different behavioral outcomes and to control the hand’s trajectory via a brainmachine interface. For example, Scherberger et al. [2] showed that the LFP can predict an animal’s intended movements based on singletrial information. Hwang and Andersen [7] also showed that singlechannel LFP provided more accurate target information than singlechannel spikes, and decoders using LFPs outperformed decoders using only spikes. In addition, the LFP has been used in research and disease diagnostics for decades, such as Parkinsonian tremor identification [8] and the diagnosis of Alzheimer’s disease [9]. However, surprisingly little is known about the functional connectivity properties among LFPs and their relationships with behavioral outcomes.
In the last few years, functional connectivity has become a useful tool for exploring information related to directed functional interactions and the mechanism responsible for their integration in neuroscience. In practical applications, neural functional networks can be represented as a graph comprising statistical correlations between neural signals, where the nodes may be neurons or cortical areas, and the networks can be weighted or unweighted, directed or undirected. Xie et al. [10] demonstrated that the functional connectivity strengthened and the information transfer efficiency increased during a rat working memory task. Lu et al. [11] described a networkbased method for predicting the behavioral outcomes of single trials in rat memory tasks. In addition, many network features have also been introduced as neurophysiological biomarkers of disease in clinical application. For example, the smallworldness of functional networks can be used as a diagnostic criterion for Alzheimer’s disease [12–14]. Thus, neural functional network analysis may be a useful tool for decoding behavior in animals.
Evidence suggests that the prefrontal cortex (PFC) is critical for goaldirected behavior in mammals [15–17], where it organizes goaldirected actions. The avian nidopallium caudolaterale (NCL) resembles the PFC in terms of its anatomical connectivity [18], neurochemical organization [19], receptor architecture [20], and functional characteristics [21–23]. Studies have shown that [11, 24] neural activities can represent behavioral information in the mammalian PFC with different firing patterns. Several studies [10, 11, 25] have employed the features of spike functional networks to decode different behavioral outcomes. However, it is not clear whether the functional connectivity of LFPs from the avian NCL can be used to represent behavioral information. In addition, it is not known whether LFP network features can decode the behavioral outcome in a single trial.
To address these problems, we developed a networkbased feature extraction method to decode pigeon behavior outcomes in single trials during a goaldirected decisionmaking task. In the task, pigeons were trained to turn left, move forward, or turn right at the choice point in a maze with food at the ends of the reward arms. The LFP signals were recorded from the NCL during turning by pigeons. The LFP spectral properties were then analyzed and the feature band was extracted using a wavelet filter. A neural functional network was constructed using a synchronization likelihood (SL) index [26]. The strength values of the network connections were unwrapped into row vectors and their dimensions were reduced by principal component analysis (PCA) as the input for the decoder. Finally, the behavioral outcomes of single trials were decoded using leaveoneout combined with the nearest neighbor method. The proposed method is suitable for analyzing the characteristics of LFP functional networks and for decoding behavior in single trials.
2. Materials and Methods
2.1. Subjects
Six adult pigeons (Columba livia) of unknown sex (450550 g) were used in this study. All of the experiments were conducted in accordance with the Animals Act, 2006 (China), for the care and use of laboratory animals, and approved by the Life Science Ethical Review Committee of Zhengzhou University. Drug usage in the experiments complied with the Chinese Pharmacopoeia (2010 edition), approved by the Chinese Pharmacopoeia Commission. Details of behavior apparatus and experimental protocols of the pigeons are described elsewhere [27]; here we only do a brief introduction.
2.2. Behavioral Apparatus and Protocol
Figure 1 shows behavioral task and experiment procedures. Pigeons performed a goaldirected decisionmaking task in a plus maze, which included a waiting area, a choice arm, and three reward arms (Figure 1(a)). There was an automatic door between the waiting area and choice arm as well as four automatic food hampers at the waiting area and the ends of three reward arms (Figure 1(b)). Besides, four infrared sensors were located in the choice arm and in the three reward arms close to the intersection. The infrared sensor in the choice arm was defined as the turning begin sensor (TBS) and the infrared sensor in the reward arm was defined as the turning end sensor (TES). In the task, pigeons were trained to turn left, move forward, or turn right, which were random in each trial. But the times of each direction were approximately equal in a session. Pigeons ran from the waiting area to the reward arms and received a food reward, indicating the completion of a trial. The pigeons had to go back to the waiting area to start the next trial after the reward was consumed. If the trial was performed correctly, the pigeon received an additional reward in the waiting area. If a pigeon performed reliably in the trial, where the required correct attempt rate was more than 90% on two consecutive days, it was considered ready for electrode implantation. The experimental procedures are shown in Figure 1(c).
(a)
(b)
(c)
2.3. Surgery and Recordings
All surgical procedures were performed, while the animals were under general anesthesia. The head was placed in a stereotaxic holder, which was customized for pigeons, with the anterior fixation point (i.e., beak bar position) 45° below the horizontal axis of the instrument. Data were collected using 8 or 16channel nickelchromium microwire arrays (Hong Kong Plexon Inc., Hong Kong, China), which were chronically implanted in the left NCL area (a 5.5 mm; L 7.5 mm; D 2.54.5 mm), according to the atlas provided by Karten and Hodos [28]. Following 57 days of recovery, neural signals were recorded from the avian NCL with a Cerebus™ recording system (Blackrock Microsystem Inc., Salt Lake City, USA). The signals were amplified (4000x), filtered (0250 Hz), and sampled continuously at 2 kHz, before saving the files as LFPs. The LFPs were downsampled to 1 kHz after removing the baseline wander. To minimize potential contamination by power line noise, the LFPs were filtered using an adaptive common average reference method to remove the spatially correlated artifacts [29].
2.4. Data Analysis and Statistics
We recorded LFP signals from NCL in pigeons, while they performed the goaldirected decisionmaking task in a plus maze. In this study, we only analyzed the data with a completion time 60 s in a trial and correct trial data. The number of electrodes, electrode types, number of channels, and number of trials per pigeon are shown in Table 1. Figure 2 shows flowchart of neural functional network decoding. The decoding algorithm includes determination of dominant frequency band, functional network construction, feature extraction, and decoding, specific as follows. All data analyses were performed using MATLAB 7.14 (r2012a) by MathWorks (The MathWorks Inc., Natick, USA).

2.4.1. Frequency Spectral Analysis
Spectral analysis was used to assess the dominant frequency bands in the LFPs. To illustrate the temporal modulation of energy in different frequency bands, the LFP timefrequency spectrum was calculated using Complex Morlet’s wavelet with 1 Hz resolution. And the LFP power spectrum was performed using the multitaper method with the timebandwidth parameter . The dominant frequency band was extracted by a wavelet filter due to its advantage of wavelet transform in nonstationary signal analysis. A quadratic Bspline mother wavelet was used as the wavelet filter because it has compact support, while it is smooth, symmetric, and, most importantly, it can provide optimal timefrequency resolution [30, 31].
2.4.2. Network Construction
Using the characteristic frequency band, the functional connectivity of LFPs was constructed using the SL method. Although there is many methods for measure of functional connectivity, such as granger, crosscorrelation, multivariate regression model, and Stransform, the SL is arguably the most popular index for estimating functional connectivity in neurophysiological data [26, 32]. It provides a nonlinear estimate of the dynamical interdependencies and relies on the detection of simultaneously occurring patterns, which can be complex and very different among signals.
If we suppose that the microelectrode array records channels for LFP signals in a trial, , , then the corresponding dimensional delayed vectors at time are defined aswhere is the delay time. The probability that two embedded vectors from an LFP () are closer to each other than a given distance at time is given bywhere is the Euclidean distance operator; is the Theiler window, which is used to avoid an autocorrelation effect on the calculations; and is a window that sharpens the time resolution of the synchronization measure. The values of and should satisfy [33], and should be at least in the order of the autocorrelation time. is the Heaviside step function as follows:
For each of the LFP signals and each time , the critical distance is determined for : that is, , where denotes the percentage of reconstructed state vectors in that are sufficiently close to to be regarded as dynamically equivalent to them. For each discrete time pair , the number of channels, , where the embedded vectors and are closer together than this critical distance within the time window , iswhere denotes how many of the embedded signals “resemble” each other. Then, for each channel and discrete time pair , the SL, , is defined as
By averaging over all , we finally obtain the SL, :where describes how strongly channel at time is synchronized with all the other channels. An SL value for the th channel, , is considered by averaging for all times . For channels in LFP signals, [32], denotes that all LFPs are uncorrelated, and denotes the maximal synchronization of all LFPs.
Calculation the SL of all channels, the functional connectivity matrix is obtained bywhere and . In this study, the SL indexes for LFP signals are calculated by the HERMES toolbox [32].
2.4.3. Feature Extraction
To extract the network features of the LFP functional network, given that the functional connectivity matrix is symmetric, only the lower triangle of the matrix was considered. Then, the connection strength values, obtained from the lowertriangular matrix, are reshaped into a row vector: that is,where the length of the vector is .
For trials, a new functional connectivity matrix is obtained by combining with the equivalent connection strengths from all trials:where denotes the functional connectivity vector in the th trial, . The dimension of is , that is, trials edges. The dimension of is then reduced using the PCA algorithm. The first few principal components (PCs) where their sum accounts for more than 90% of the total variance in energy in the original data [34] are selected to generate the new functional connectivity matrix (), which is defined as a network feature. Here, the PCA algorithm was performed using the function princomp in MATLAB.
2.4.4. Decoding Algorithm
For the matrix , leaveoneoutput decoding was performed with the nearest neighbor method [35]. Specifically, for each trial, the event type was decoded based on the distribution of all the remaining trials. To illustrate the neural decoding of LFP, can be rewritten as follows:where denotes the th trial and is the number of the extracted PCs. We also assumed that PCs could encode common behavioral events in each trial, that is, turning left, moving forward, or turning right. Thus, for all trials, the behavior events can be expressed aswhere , denotes turning left, is moving forward, and is turning right.
In the decoding algorithm, we had pairs . If a trial was selected as the test dataset, then other trials were in the training dataset, and the training dataset was labeled as the class. Next, the class of the test dataset was decoded using the nearest neighbor (NN) method. We repeat this step until the classes of all trials are predicted. Finally, the predicted labels and their real labels are compared, and the accuracy is calculated as follows: , where is the number of correct predictions.
The NN algorithm is a nonparametric classification method that classifies points based on actual examples in the training set [35]. Each point in the test set is assigned the same class as the majority vote of its NN in the training set. Therefore, there is little or no prior knowledge about the distribution of data points. In this study, the dataset was split up into two partitions, where one partition was used as a training set and the other partition was used as a test set. The value was always less than the square root of the size of training samples [36]. The NN algorithm was performed using the function knnclassify in MATLAB.
2.4.5. Network Properties
To measure the properties of the functional networks during the goaldirected task, the global efficiency () and clustering coefficient () were calculated, which described the global properties of the network. Before calculating the network properties, traditional methods usually involve converting the network into a binary network, which is formed by a threshold. The value in the matrix larger than the threshold is set to be 1, while less than the threshold is 0. The threshold is difficult to choose [37]. In order to improve analysis, we directly analyzed the weighted functional connection matrix.
is a measure of the capacity of nodes to propagate information in parallel across a network, which is defined mathematically as [38]where is the number of nodes, is the shortest path length between node and , and . A network with a high suggests the existence of a strong association between nodes. Thus, information is communicated with a high level of efficiency throughout the network.
is a measure of the degree to which the nodes in a network tend to cluster together, which is computed by [38]where is the number of nodes, is the degree of node , and , where , . A larger value for denotes the more close the network connected and the higher the efficiency of information transfer between the nodes.
2.4.6. Statistical Analysis
For the network measures, the indexes used in the text and figures are expressed as the mean standard deviation (s.d.). Statistical differences were evaluated with a Wilcoxon ranksum test, where the significance level was set to 5%. values were considered to indicate significant differences at . Statistical analyses were conducted using the MATLAB Statistics and Machine Learning Toolbox.
3. Results
3.1. Spectral Properties of LFPs during the GoalDirected Task
First, we analyzed the LFPs and the spectral properties recorded from the pigeon NCL during the goaldirected decisionmaking task, where the results are shown in Figure 3. Illustration of both the event window and the implanting location of microelectrode array are showed by Figures 3(a) and 3(b), respectively. Figure 3(c) shows timefrequency spectrum of a single electrode and averaged timefrequency spectrum across all the electrodes during a trial. The signals below 8 Hz were filtered. Figure 3(d) shows the averaged timefrequency spectrum across all the trials for turning left, moving forward, and turning right. A 3 s window was used and time 0 indicated the moment when the pigeon was at 1 s before the TBS. Figure 3(e) shows examples of LFP waveforms (gammaband) during a trial. A dominant highfrequency oscillation was observed around the TES, which indicates that the highfrequency oscillations in the LFP may be related to the behavior of pigeons.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
The LFP power spectra were calculated for the baseline period (3 s before the gate opened) and the turning period (3 s after the TBS). Figure 3(f) shows the power spectra of a single electrode in all trials, which indicates that a broad band increased significantly in the turning area compared with the baseline in the waiting area located in a band centered around 55 Hz. At most of the recording sites, the LFP power during turning increased by 3.5fold in the gammaband (4060 Hz) of the frequency spectrum (Figure 3(g)). The increase in the power value was dependent on the pigeon (mean power increase: P010 = 33.9 dB, P600 = 12.1 dB, P601 = 15.0 dB, P605 = 7.6 dB, P609 = 3.5 dB, and P619 = 22.4 dB), but an increase in the gammaband power was prominent in all six pigeons.
3.2. Functional Connectivity Properties of LFP Networks
To measure the functional connectivity properties of LFP networks, we obtained the gammaband of the LFP using a wavelet filter. The mother wavelet is bior 2.2 and the decomposition scale was equal to 4. Therefore, the frequency range of the extracted gammaband was 3162 Hz. The LFP network was then constructed using the gammaband of the LFP. The functional connectivity features of the LFP network were measured by calculating the global efficiency and the clustering coefficient. We found that both and increased significantly in all six pigeons (Figure 4; Wilcoxon ranksum test, ). The results were similar as the power increased. Therefore, we propose that the LFP network may be related to the goaldirected behavior of pigeons.
(a)
(b)
To determine the decoding window for pigeon behavior outcomes, we calculated and for the LFP network before and after turning, that is, 1 s after TBS (TB) and 1 s before TES (TE), and the results are shown in Figure 5. Illustration of the different event windows was showed by Figure 5(a). Compared with TB, the connectivity of TE was stronger (Figure 5(b)). And and were significantly higher in TE than TB (Wilcoxon ranksum test, , Figure 5(c)). These results indicate that the functional connectivity was related more strongly to the goaldirected behavior of pigeons in TE than TB. Therefore, TE was used to decode the pigeon behavior outcome in a single trial.
(a)
(b)
(c)
3.3. Decoding Pigeon Behavior Outcomes
According to the results described previously, changes in the neuronal functional network during turning were related to the behavioral choices of pigeons. However, it was not clear whether these neuronal functional networks could be employed to decode the behavioral choices in a single trial. Thus, we decoded the pigeon behavior outcome in a single trial using the network features in the TE and we compared the results with those obtained using the power features, which are employed widely for behavior prediction. The power features were obtained as follows. The power of the gammaband from all channels was first calculated in the TE and for all trials; the power dimensions were reduced using the PCA algorithm. The power with decreased dimensions was defined as the power feature.
Figure 6 shows the decoding results obtained using network features and power features. As shown in Figure 6(a), the separability of the network features in the three directions was better than that of the power features. The accuracy of network features was significantly higher than that of the power features with different values of (Figure 6(b); ) except , which is denoted with arrow. Because the mean size of training samples is about 50 (see Table 1; ), the maximum value was set to 7. The largest accuracy was obtained with . Figure 6(c) shows the decoding accuracy using power feature and network features when for the six pigeons (power features ; network features ; Wilcoxon ranksum test, ). The results have shown that the LFP functional connections from the NCL could be used to decode the pigeon behavior outcomes of single trials, and the accuracy of the network features was better than that of the power features.
(a)
(b)
(c)
(d)
(e)
Next, we studied the time profile of the decoding performance based on the network features. Decoding was performed using a sliding window with a width of 1 s and a step size of 100 ms [10]. Time 0 indicated the moment when the pigeon was at the TES point. The accuracy peaked at the TES, except in P600, and the maximum of all pigeons was above the 80% (Figure 6(d)). The decoding accuracy differs significantly from the channel level (0.33) at the beginning, 0.7 s before the TES (Figure 6(e); Wilcoxon ranksum test, ), thereby confirming that the reward may also play a crucial role in the goaldirected behavior of pigeons.
4. Discussion
In the present study, we explored the functional connectivity properties of LFPs recorded from the pigeon NCL during a goaldirected decisionmaking task. We developed a networkbased feature extraction method based on the SL index to decode the behavior of animals. We found that the power of the gammaband increased significantly during turning by pigeons and the LFP functional network based on the gammaband correlated with the goaldirected behavior of the pigeons. Moreover, the functional connectivity strength values among LFPs could decode the behavioral choices of pigeons in an effective manner. The accuracy of network features was significantly higher than that of power features as well as the chance level. The results showed that the proposed method is effective for decoding the taskrelated behavior of animals based on the functional networks obtained from multielectrode recordings.
In this study, the experimental results were controlled by six parameters: embedding dimension (), embedding delay (), Theiler window (), and for the SL index, and the nearest neighbor in the NN method. The first five parameters directly affected the structures of the neuronal functional networks. In this study, the SL index was set according to previously described parameter settings in the HERMES toolbox [32], so here we do not provide full details of these parameters. The nearest neighbor, , controlled the decoding results, where we found that a smaller value for obtained better decoding performance in the experiment: that is, the performance was generally better when . However, comparing with the power features, the network features have better performance for all values.
Decoding animal behavior based on the activity patters of neurons is a key issue in neuroscience. Previous decoding approaches studied how relevant information about the stimuli or events is represented by the neuronal population activity [39]: for example, predicting the position of a rat in its environment from recordings in hippocampus [40], arm movements from motor cortex in human [41], and image presentations from spiking activity in human medial temporal lobe [39]. Here, we demonstrated the feasibility of decoding behavior outcomes of pigeon using LFP functional connections in the NCL. Functional connective network, which contains information about the correlations between neurons or channels, may be an efficient method to decode the behavior of the animals [37]. By applying the SL method to determine the correlation between all pairwise of LFP signals, we found that the functional connectivity strength values can distinguish the behavior outcomes of pigeons, and the accuracy was significantly higher than that of power features.
Similar to the proposed method, Rosenberg et al. [42] used the strength values of brain connections based on fMRI data to predict a subject’s attention capacity. By contrast, we used the functional connectivity strengths based on the LFP data obtained from the NCL to decode the behavioral outcomes of pigeon in this study. The dimensions of the connectivity strengths were also reduced by PCA, thereby improving the computational speed and avoiding the curse of dimensionality. Lu et al. [11] used the spike data from the PFC and hippocampal CA1 regions to build wholerecorded neuronal functional networks and then divided these networks into local neuronal circuit groups based on the maximization of modularity Q to predict rat behavior. Our study, which used the functional connectivity strength values among the LFP data to decode the pigeon behavior, expanded upon the work, and the results add to the growing body of evidence that neuronal functional network which exists in the brain of animals can decode the behavior of animals effectively.
Goaldirected decisionmaking is a complex behavior, requiring the subject to perceive its environment, learn about the significant of the environment, and then select where to go next [43]. Although rodents are classic model animals for the study of goaldirected behavior [15, 16], the wealth of available behavioral and neuroanatomical data renders the pigeons a highly suitable model system. Moreover, we have demonstrated experimentally that the NCL plays an important role during the goaldirected behavior of pigeons using the spike signals [27]. The results of this study are in line with the conclusion. The functional connections from LFP’s gammaband signals in the NCL decoded effectively the pigeon behavioral outcomes during the goaldirected task. But whether the functional connections from the spike signals could also decode the behavioral outcomes needs further study.
In addition, it should be noted that the number of recorded channels was limited in our experiments. The greater the number of channels is recorded, the clearer the relationships between the neural functional networks and behavior could be revealed [11]. But for the SL method, two important limitations for the practical use are its computational and memory costs of current implementations. As the number of channels increases, both the computational and the memory costs increase sharply. Moreover, the decoding approach used was relatively simple; the complex classification methods may get better results. And the decoding window might also not be optimal; we found that the accuracy of the network features in some pigeons, such as P601, P605, P609, and P619, still increased after the birds passed the TES. However, the method proposed in this study still provides meaningful information that may be suitable for future analysis of largescale neuronal networks. And the obtained results supported that the functional connections among LFPs recoded from the NCL can decode pigeon behavior outcomes.
In summary, we developed a networkbased method for decoding the behavior outcomes of pigeons. In the proposed method, the LFP functional connectivity strength values with reduced dimensions are set as the input for the decoder, which avoids setting thresholds and the curse of dimensionality. Moreover, the accuracy of the network features was significantly higher than that of the power features and the chance level. The results have shown that the NCL neurons contain the information about the goaldirected behavior of pigeons, and the functional connections among LFPs could decode pigeon behavior outcomes. The proposed method provides a powerful tool for decoding animal behavior outcomes using functional networks.
Conflicts of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Authors’ Contributions
Yan Chen and Xinyu Liu contributed equally to this work.
Acknowledgments
The work was supported in part by the National Natural Science Foundation of China (61673353, U1304602, and 61603344), Program for Science and Technology of Henan Province of China (182102210099 and 162102310167), and Open Foundation of Henan Key Laboratory of Brain Science and BrainComputer Interface Technology (HNBBL17005).
References
 S. Katzner, I. Nauhaus, A. Benucci, V. Bonin, D. L. Ringach, and M. Carandini, “Local origin of field potentials in visual cortex,” Neuron, vol. 61, no. 1, pp. 35–41, 2009. View at: Publisher Site  Google Scholar
 H. Scherberger, M. R. Jarvis, and R. A. Andersen, “Cortical local field potential encodes movement intentions in the posterior parietal cortex,” Neuron, vol. 46, no. 2, pp. 347–354, 2005. View at: Publisher Site  Google Scholar
 C. Magri, U. Schridde, Y. Murayama, S. Panzeri, and N. K. Logothetis, “The amplitude and timing of the BOLD signal reflects the relationship between local field potential power at different frequencies,” The Journal of Neuroscience, vol. 32, no. 4, pp. 1396–1407, 2012. View at: Publisher Site  Google Scholar
 G. Buzsáki, “Largescale recording of neuronal ensembles,” Nature Neuroscience, vol. 7, no. 5, pp. 446–451, 2004. View at: Publisher Site  Google Scholar
 N. K. Logothetis and B. A. Wandell, “Interpreting the BOLD signal,” Annual Review of Physiology, vol. 66, pp. 735–769, 2004. View at: Publisher Site  Google Scholar
 Q. Gaucher, J. Edeline, and B. Gourévitch, “How different are the local field potentials and spiking activities? Insights from multielectrodes arrays,” Journal of Physiology—Paris, vol. 106, no. 34, pp. 93–103, 2012. View at: Publisher Site  Google Scholar
 E. J. Hwang and R. A. Andersen, “The utility of multichannel local field potentials for brain–machine interfaces,” Journal of Neural Engineering, vol. 10, no. 4, p. 046005, 2013. View at: Publisher Site  Google Scholar
 E. Bakstein, J. Burgess, K. Warwick, V. Ruiz, T. Aziz, and J. Stein, “Parkinsonian tremor identification with multiple local field potential feature classification,” Journal of Neuroscience Methods, vol. 209, no. 2, pp. 320–330, 2012. View at: Publisher Site  Google Scholar
 D. S. Roy, A. Arons, T. I. Mitchell, M. Pignatelli, T. J. Ryan, and S. Tonegawa, “Memory retrieval by activating engram cells in mouse models of early Alzheimer's disease,” Nature, vol. 531, no. 7595, pp. 508–512, 2016. View at: Publisher Site  Google Scholar
 J. Xie, W. Bai, T. Liu, and X. Tian, “Functional connectivity among spike trains in neural assemblies during rat working memory task,” Behavioural Brain Research, vol. 274, pp. 248–257, 2014. View at: Publisher Site  Google Scholar
 H. Lu, S. Yang, L. Lin, B. Li, and H. Wei, “Prediction of rat behavior outcomes in memory tasks using functional connections among neurons,” PLoS ONE, vol. 8, no. 9, p. e74298, 2013. View at: Publisher Site  Google Scholar
 F. Bai, N. Shu, Y. Yuan et al., “Topologically convergent and divergent structural connectivity patterns between patients with remitted geriatric depression and amnestic mild cognitive impairment,” The Journal of Neuroscience, vol. 32, no. 12, pp. 4307–4318, 2012. View at: Publisher Site  Google Scholar
 B. M. Tijms, A. M. Wink, W. de Haan et al., “Alzheimer's disease: connecting findings from graph theoretical studies of brain networks,” Neurobiology of Aging, vol. 34, no. 8, pp. 2023–2036, 2013. View at: Publisher Site  Google Scholar
 J. Wang, X. Zuo, Z. Dai et al., “Disrupted functional brain connectome in individuals at risk for Alzheimer's disease,” Biological Psychiatry, vol. 73, no. 5, pp. 472–481, 2013. View at: Publisher Site  Google Scholar
 V. Hok, E. Save, P. P. LenckSantini, and B. Poucet, “Coding for spatial goals in the prelimbic/infralimbic area of the rat frontal cortex,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 102, no. 12, pp. 4602–4607, 2005. View at: Publisher Site  Google Scholar
 V. Hok, E. Chah, E. Save, and B. Poucet, “Prefrontal cortex focally modulates hippocampal place cell firing patterns,” The Journal of Neuroscience, vol. 33, no. 8, pp. 3443–3451, 2013. View at: Publisher Site  Google Scholar
 J. C. Churchwell, A. M. Morris, N. D. Musso, and R. P. Kesner, “Prefrontal and hippocampal contributions to encoding and retrieval of spatial memory,” Neurobiology of Learning and Memory, vol. 93, no. 3, pp. 415–421, 2010. View at: Publisher Site  Google Scholar
 S. Kröner and O. Güntürkün, “Afferent and efferent connections of the caudolateral neostriatum in the pigeon (columba livia): A retro and anterograde pathway tracing study,” Journal of Comparative Neurology, vol. 407, no. 2, pp. 228–260, 1999. View at: Publisher Site  Google Scholar
 T. Bast, B. Diekamp, C. Thiel, R. K. W. Schwarting, and O. Gntrkn, “Functional aspects of dopamine metabolism in the putative prefrontal cortex analogue and striatum of pigeons (Columba livia),” Journal of Comparative Neurology, vol. 446, no. 1, pp. 58–67, 2002. View at: Publisher Site  Google Scholar
 C. Herold, N. PalomeroGallagher, B. Hellmann et al., “The receptor architecture of the pigeons' nidopallium caudolaterale: An avian analogue to the mammalian prefrontal cortex,” Brain Structure & Function, vol. 216, no. 3, pp. 239–254, 2011. View at: Publisher Site  Google Scholar
 T. Kalt, B. Diekamp, and O. Güntürkün, “Single unit activity during a Go/NoGo task in the “prefrontal cortex” of pigeons,” Brain Research, vol. 839, no. 2, pp. 263–278, 1999. View at: Publisher Site  Google Scholar
 J. A. Kirsch, I. Vlachos, M. Hausmann et al., “Neuronal encoding of meaning: Establishing categoryselective response patterns in the avian 'prefrontal cortex',” Behavioural Brain Research, vol. 198, no. 1, pp. 214–223, 2009. View at: Publisher Site  Google Scholar
 S. Starosta, O. Güntürkün, and M. C. Stüttgen, “StimulusResponseOutcome Coding in the Pigeon Nidopallium Caudolaterale,” PLoS ONE, vol. 8, no. 2, Article ID e57407, 2013. View at: Publisher Site  Google Scholar
 W. J. Kargo, B. Szatmary, and D. A. Nitz, “Adaptation of prefrontal cortical firing patterns and their fidelity to changes in actionreward contingencies,” The Journal of Neuroscience, vol. 27, no. 13, pp. 3548–3559, 2007. View at: Publisher Site  Google Scholar
 A. S. Morcos and C. D. Harvey, “Historydependent variability in population dynamics during evidence accumulation in cortex,” Nature Neuroscience, vol. 19, no. 12, pp. 1672–1681, 2016. View at: Publisher Site  Google Scholar
 C. J. Stam and B. W. van Dijk, “Synchronization likelihood: an unbiased measure of generalized synchronization in multivariate data sets,” Physica D: Nonlinear Phenomena, vol. 163, no. 34, pp. 236–251, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 X. Liu, H. Wan, S. Li, Z. Shang, and L. Shi, “The role of nidopallium caudolaterale in the goaldirected behavior of pigeons,” Behavioural Brain Research, vol. 326, pp. 112–120, 2017. View at: Publisher Site  Google Scholar
 H. J. Karten and W. A. Hodos, A Stereotaxic atlas of the brain of the pigeon (Columba Livia), Johns Hopkins Press, Baltimore, 1967.
 X. Liu, H. Wan, S. Li, Y. Chen, and L. Shi, “Adaptive common average reference for in vivo multichannel local field potentials,” Biomedical Engineering Letters, vol. 7, no. 1, pp. 7–15, 2017. View at: Publisher Site  Google Scholar
 X. Liu, H. Wan, Z. Shang, and L. Shi, “Automatic extracellular spike denoising using wavelet neighbor coefficients and level dependency,” Neurocomputing, vol. 149, pp. 1407–1414, 2015. View at: Publisher Site  Google Scholar
 M. Ahmadi and R. Quian Quiroga, “Automatic denoising of singletrial evoked potentials,” NeuroImage, vol. 66, pp. 672–680, 2013. View at: Publisher Site  Google Scholar
 G. Niso, R. Bruña, E. Pereda et al., “HERMES: towards an integrated toolbox to characterize functional and effective brain connectivity,” Neuroinformatics, vol. 11, no. 4, pp. 405–434, 2013. View at: Publisher Site  Google Scholar
 T. Montez, K. LinkenkaerHansen, B. W. van Dijk, and C. J. Stam, “Synchronization likelihood with explicit timefrequency priors,” NeuroImage, vol. 33, no. 4, pp. 1117–1125, 2006. View at: Publisher Site  Google Scholar
 H. Abdi and L. J. Williams, “Principal component analysis,” Wiley Interdisciplinary Reviews: Computational Statistics, vol. 2, no. 4, pp. 433–459, 2010. View at: Publisher Site  Google Scholar
 L. Veit and A. Nieder, “Abstract rule neurons in the endbrain support intelligent behaviour in corvid songbirds,” Nature Communications, vol. 4, article no. 2878, 2013. View at: Publisher Site  Google Scholar
 T. Hastie and R. Tibshirani, “Discriminant adaptive nearest neighbor classification,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 18, no. 6, pp. 607–616, 1996. View at: Publisher Site  Google Scholar
 E. Bullmore and O. Sporns, “Complex brain networks: graph theoretical analysis of structural and functional systems,” Nature Reviews Neuroscience, vol. 10, no. 3, pp. 186–198, 2009. View at: Publisher Site  Google Scholar
 M. Rubinov and O. Sporns, “Complex network measures of brain connectivity: Uses and interpretations,” NeuroImage, vol. 52, no. 3, pp. 1059–1069, 2010. View at: Publisher Site  Google Scholar
 R. Q. Quiroga, L. Reddy, C. Koch, and I. Fried, “Decoding visual inputs from multiple neurons in the human temporal lobe,” Journal of Neurophysiology, vol. 98, no. 4, pp. 1997–2007, 2007. View at: Publisher Site  Google Scholar
 X. Liu, X. Hai, Z. Shang, and H. Wan, “Decoding movement trajectory of hippocampal place cells by particle filter,” Prog. Biochem. Biophys, vol. 43, no. 8, 2016. View at: Publisher Site  Google Scholar
 L. R. Hochberg, D. Bacher, B. Jarosiewicz et al., “Reach and grasp by people with tetraplegia using a neurally controlled robotic arm,” Nature, vol. 485, no. 7398, pp. 372–375, 2012. View at: Publisher Site  Google Scholar
 M. D. Rosenberg, E. S. Finn, D. Scheinost et al., “A neuromarker of sustained attention from wholebrain functional connectivity,” Nature Neuroscience, vol. 19, no. 1, pp. 165–171, 2015. View at: Publisher Site  Google Scholar
 M. R. Penner and S. J. Mizumori, “Neural systems analysis of decision making during goaldirected navigation,” Progress in Neurobiology, vol. 96, no. 1, pp. 96–135, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Yan Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.