Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE
PreviousNext
Articles, Systems/Circuits

Long-Term Recordings Improve the Detection of Weak Excitatory–Excitatory Connections in Rat Prefrontal Cortex

C. Daniela Schwindel, Karim Ali, Bruce L. McNaughton and Masami Tatsuno
Journal of Neuroscience 16 April 2014, 34 (16) 5454-5467; https://doi.org/10.1523/JNEUROSCI.4350-13.2014
C. Daniela Schwindel
Department of Neuroscience, Canadian Centre for Behavioral Neuroscience, The University of Lethbridge, Lethbridge, AB, T1K 3M4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Karim Ali
Department of Neuroscience, Canadian Centre for Behavioral Neuroscience, The University of Lethbridge, Lethbridge, AB, T1K 3M4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bruce L. McNaughton
Department of Neuroscience, Canadian Centre for Behavioral Neuroscience, The University of Lethbridge, Lethbridge, AB, T1K 3M4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Masami Tatsuno
Department of Neuroscience, Canadian Centre for Behavioral Neuroscience, The University of Lethbridge, Lethbridge, AB, T1K 3M4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Characterization of synaptic connectivity is essential to understanding neural circuit dynamics. For extracellularly recorded spike trains, indirect evidence for connectivity can be inferred from short-latency peaks in the correlogram between two neurons. Despite their predominance in cortex, however, significant interactions between excitatory neurons (E) have been hard to detect because of their intrinsic weakness. By taking advantage of long duration recordings, up to 25 h, from rat prefrontal cortex, we found that 7.6% of the recorded pyramidal neurons are connected. This corresponds to ∼70% of the local E–E connection probability that has been reported by paired intracellular recordings (11.6%). This value is significantly higher than previous reports from extracellular recordings, but still a substantial underestimate. Our analysis showed that long recording times and strict significance thresholds are necessary to detect weak connections while avoiding false-positive results, but will likely still leave many excitatory connections undetected. In addition, we found that hyper-reciprocity of connections in prefrontal cortex that was shown previously by paired intracellular recordings was only present in short-distance, but not in long distance (∼300 micrometers or more) interactions. As hyper-reciprocity is restricted to local clusters, it might be a minicolumnar effect. Given the current surge of interest in very high-density neural spike recording (e.g., NIH BRAIN Project) it is of paramount importance that we have statistically reliable methods for estimating connectivity from cross-correlation analysis available. We provide an important step in this direction.

  • cross correlations
  • extracellular recordings
  • local connectivity
  • medial prefrontal cortex

Introduction

Characterization of synaptic interactions on a large scale is essential to understand information processing in neural circuits. In an attempt to characterize local circuit dynamics, combined electrophysiological and imaging techniques have started to map out neocortical circuits (Bock et al., 2011; Ko et al., 2011). With recent developments in multielectrode recording technology, it has also become possible to monitor large numbers of neurons simultaneously with high temporal resolution (McNaughton et al., 1983; Csicsvari et al., 2003; Buzsaki, 2004). Even though this number is still small relative to the total number of neurons, one advantage of ensemble recordings is that it can greatly facilitate the study of spike-train interactions, because the number of neuron pairs increases as the square of the number of units recorded. The obtained spike trains are often analyzed by cross-correlations to infer synaptic interactions (Alonso and Martinez, 1998; Ostojic et al., 2009). The cross-correlogram, a well established technique to investigate temporal relationships between neural spikes, describes the covariance between the binned spike trains of two neurons at various time lags (Perkel et al., 1967b; Kirkwood, 1979; Aertsen and Gerstein, 1985; Brown et al., 2004). A plausible argument for synaptic connectivity can be made from the presence of short-latency peaks in the correlogram within the range of central glutamatergic EPSP and GABAergic IPSP rise times (Csicsvari et al., 1998; Barthó et al., 2004; Fujisawa et al., 2008).

Despite the potential usefulness of the cross-correlation method, the detection of excitatory interactions between pyramidal neurons (E–E interaction) has been difficult (Barthó et al., 2004; Fujisawa et al., 2008), because these interactions are generally weak (McNaughton, 1980; Mason et al., 1991; Deuchars et al., 1994; Markram et al., 1997; Thomson and Deuchars, 1997; Reyes and Sakmann, 1999; Thomson et al., 2002). Excitatory synaptic strength of intracellularly measured cortical connections has been shown to follow a lognormal distribution (Song et al., 2005), suggesting that a very small number of strong connections are embedded in a large number of weak connections. In addition, the duration of typical ensemble recordings is not long (up to a few hours). The consequence is a limited sampling of spike occurrences and also a poor sampling of the possible state-space occupancy distribution (Perkel et al., 1967a,b). Cells are tuned to respond to a given set of input vectors and some vectors in the set may never occur during the sampling epoch. Thus, for some connected cells in the sample, the input which brings them close enough to threshold to allow synaptic interaction may not occur during the sampling period. In other words, the effective contribution of cell A to the firing of cell B is not independent of the activity of other cells in the network. Hence, some connections may not be visible in spike cross-correlograms from a given sample epoch.

To explore these issues, we analyzed neural ensembles from rat prefrontal cortex and estimated the asymptotic detection probability of excitatory connections and investigated the nature of the distribution of the detected excitatory connections.

Materials and Methods

Recording experiment.

Twenty-five hour continuous, multineuron recording datasets were collected from three adult male Brown Norway/Fischer 344 hybrid rats that were trained to run a spatial sequence task (Euston et al., 2007) and were also subjected to a novel object exposure procedure (Tatsuno et al., 2006). The recording sessions were divided into two 12 h sleep/rest sessions interrupted by 1 h behavior, which was either the sequence task or the novel object exposure, to yield a total recording duration of 25 h per recording session. The animals were implanted with a hyperdrive with 12 independently movable tetrodes (Gothard et al., 1996; <1 MΩ impedance) in the left medial prefrontal cortex (3.2 mm AP, 1.3 mm ML relative to bregma) to record from dorsal anterior cingulate cortex and prelimbic cortex. The tetrodes were used to record extracellular single units. The anatomical positions of the tetrodes were confirmed histologically.

Two additional probes were placed 4–5 mm deep in the medial prefrontal cortex to record a differential reference signal. The animal's position was tracked at 60 frames/s using LEDs mounted on the headstage that were detected by a color camera that was mounted on the ceiling of the recording room (∼0.33 cm/pixel). The thresholded signals were recorded with a Cheetah Data Acquisition System (Neuralynx), digitized at 32 kHz and bandpass filtered between 600 Hz and 6 kHz. Local field potential data were filtered between 1 and 475 Hz and sampled at 2 kHz.

Animal care and surgeries were conducted in accordance with National Institutes of Health guidelines and approved Institutional Animal Care and Use Committee protocols. After surgery, rats were administered 26 mg of acetaminophen (children's Tylenol, McNeil) and also received 2.7 mg/ml acetaminophen in the drinking water for 1–2 d after surgery. In addition, they were given oral ampicillin on a 10 d on/off regimen for the duration of the experiment.

Spike sorting.

Spiking activity was analyzed offline using an automated spike-sorting algorithm (KlustaKwik by K.D. Harris, University College London, London, UK) to isolate units and separate them from noise. The resulting clusters were manually refined using cluster cutting software (MClust 3.0 by A.D. Redish, University of Minnesota, Minneapolis, MN, with customizations by P. Lipa, University of Arizona, Tucson, AZ, S. L. Cowen, University of Arizona, Tucson, AZ, and D. Euston, University of Lethbridge, Canada) in a multiparameter space including features, such as energy (area under the waveform), peak to trough distance, principal component, time (to control for stability of the recording of the unit over the entire recording duration), and cross correlograms. In addition, all units were verified by waveform visualization software (WaveformCutter 1.0 by S. Cowen). Only units with <0.3% interspike intervals falling within 2 ms refractory period were accepted.

During long recording times, the position of the electrode can shift and, consequently, the shape of the recorded waveform can change. This can lead to errors in spike sorting and can cause false temporal firing relationships between neuron pairs. For this reason, we only included neurons that show good isolation and little variance in their waveforms and Z-scored peak amplitudes over the entire length of recordings (Fig. 1a). As a measure for the stability of the neurons, we calculated the fractional change of firing rate between the first and last 4 h periods. We found that the fractional change for the majority (95%) of the neurons that were selected by the first criteria was <80%, suggesting that the neurons were stable for the 25 h period of recordings.

In addition, inspired by the recent finding by Mizuseki and Buzsaki (2013) that the distribution of firing rates follows a lognormal distribution in the hippocampal formation, we also calculated the distribution of firing rate of our data. Figure 2c shows the firing rate distribution during the first and last 4 h of recording. We found that not only does the firing rate follow a lognormal distribution, but also the shape of the distribution did not change between the two periods. We also confirmed that the firing rate calculated of the entire recording follows a lognormal distribution.

Neuron classification.

Neurons are classified into excitatory, inhibitory, or unclassified based on statistical dependency with other neurons. First, for each reference neuron, cross-correlations against all other neurons were calculated with a bin size of 1 ms and a window size of ±50 ms. Second, each cross-correlation was converted to a Z-score with its mean and SD. This alleviated the problem of spike rate differences. Next, we calculated the moving average (window size: 15 ms) of the Z-scored cross-correlations and subtracted it from the original signal to obtain a detrended Z-scored cross-correlation. This procedure removed modulations in intermediate temporal ranges that are slower than monosynaptic interactions. Correlations that showed putative common input (measured as a peak or trough encompassing lag [−1:1]) were excluded from further cell classification analysis. Cross-correlation pairs that did not contain relevant information were also removed [if none of the Z-scored correlation values within [−4, 4] exceeded 2 or −2, the pair was removed]. Finally, the average cross-correlation of the remaining pairs was calculated (Fig. 1b). Mean and 1 SD of the average cross-correlation for the bins outside of [−4, 4] ms were used as a threshold to determine whether a bin within [1, 4] ms has a relevant correlation signal for cell classification; if at least one bin exceeded the upper threshold (mean + 1 SD) and no bin undershot the lower threshold (mean − 1 SD), the reference neuron was classified as a putative excitatory neuron. If at least one bin undershot the lower threshold (mean − 1 SD) and no bin exceeded the upper threshold (mean + 1 SD), the reference neuron was classified as a putative inhibitory neuron. Otherwise, the reference neuron was left unclassified. We also inspected all individual cross-correlations visually to verify whether the decision of the algorithm we applied on the average detrended Z-scored cross-correlations was consistent with the judgment made by human observers. In summary, this procedure aims at classifying a reference neuron based on the average statistical influence of the reference neuron to all possible postsynaptic target neurons.

Cross-correlation analysis.

All simultaneously recorded neurons from each animal were cross correlated with each other with a bin size of 1 ms. The cross-correlation values were normalized by the firing rate of the reference neuron. This procedure resulted in n·(n − 1)/2 correlation pairs (n = number of recorded units; self-correlation was excluded). For each neuron pair, the spikes of one spike train were jittered randomly within an interval of [−5, 5] ms to break up the short-latency relationships between spikes. The bin size and jitter interval were chosen based on synaptic integration times in frontal cortex neurons in vivo (Léger et al., 2005). Then, the jittered spike trains were cross-correlated with a bin size of 1 ms. This procedure was repeated 1000 times to compute a surrogate dataset of cross-correlations. Various significance thresholds (α levels) were computed from the distribution of jittered data to detect peaks or troughs with short latency and duration (≤2 ms) in the cell pair correlations (Barthó et al., 2004; Fujisawa et al., 2008): α levels used are 0.5, 0.2, 0.1, 0.05, 0.02, 0.01, 5 × 10−3, 2 × 10−3, 1 × 10−3, 5 × 10−4, 2 × 10−4, 1 × 10−4, 5 × 10−5, 2 × 10−5, 1 × 10−3, and the global min/max value (referred to as absolute limit: abs) which corresponds to the overall maximum and minimum values of all bins in the surrogate dataset. If a single bin exceeded the significance level within the first 4 ms before or after lag 0, the correlation was scored as significant excitation and the pair was considered monosynaptically connected. Similarly, if two neighboring bins undershot the significance level within the first 4 ms before or after lag 0, the connection was scored as significant inhibition (Barthó et al., 2004). This way, significance was tested for both possible directions of interaction between the two neurons. To verify our method and temporal window within which we assess significance of the short-latency peaks, we repeated the same analysis as described above for 4 ms windows ∼±25 ms and ±50 ms.

The length of the recording period used for the cross correlogram calculation was varied from 1 to 25 h with 1 h increments. First, for each hour segment and for a given significance level α, we constructed a matrix of significant cross-correlations for all possible nonoverlapping blocks (i.e., 25 blocks for 1 h segments, 12 blocks for 2 h segments, etc.). This resulted in a matrix with an entry for each significant correlation type EE1, EE2, EI1, EI2, IE1, II1, EU1, and IU1 (EE1 represents unidirectional excitation between two excitatory neurons; EE2 represents reciprocal excitation between excitatory neurons; EI1 represents unidirectional excitation from an excitatory neuron onto an inhibitory neuron; EI2 represents reciprocal connection between an inhibitory and an excitatory neuron; IE1 represents unidirectional inhibition from an inhibitory neuron onto an excitatory neuron; II1 represents unidirectional inhibition between two inhibitory neurons; EU1 represents unidirectional excitation from an excitatory neuron onto an unclassified neuron; IU1 represents unidirectional inhibition from an inhibitory neuron onto an unclassified neuron). Because swapping reference- and target-neurons results in the mirror image of the original cross correlogram, only the upper triangle of the result matrix was considered. We then summed the number of significant correlations for each correlation type. To estimate the distribution of the number of detected correlations for each hour segment, we used bootstrapping with replacement (Mooney and Duval, 1993; Hoffman and McNaughton, 2002). The matrix was randomly resampled with replacement as many times as needed to obtain 1000 samples for each hour segment. For instance, for the 1 h segment, each upper triangle of the 25 result matrices was resampled 40 times. For the 25 h segment, only one result matrix was obtained and was therefore resampled 1000 times. The median, 84th and 16th percentile were calculated to characterize the distribution. For the tetrode-wise analysis, the connection matrix of each tetrode was resampled with replacement 1000 times. The result was summed across tetrodes and datasets. Again, the median, 84th and 16th percentile of the bootstrapped sample distribution were calculated.

The experimentally derived detection probability of excitatory connections pαn for the significance level α and the recording length of n hours was defined as follows: Embedded Image where #(EE1) is the number of significant unidirectionally connected pairs, #(EE2) is the number of significant reciprocally connected pairs, and N is the total number of neuron pairs.

Extrapolation by fitting the statistical power function using simulated annealing.

To investigate how the experimentally derived detection probability of weak E–E connections pαn may improve beyond the recording duration of 25 h, we applied curve fitting based on statistical power. To obtain a reasonable estimate, we used a framework for the case that both the data from the alternative hypothesis and the data from the null hypothesis follow the normal distribution with a known SD. The statistical power function, defined as the probability of rejecting the null hypothesis when the alternative is true, is written as follows: Embedded Image where the cutoff is the threshold value for rejecting the null hypothesis. When the alternative is true, μ1 is the mean of the alternative hypothesis (the mean peak cross-correlation value within −4 to +4 ms excluding 0), σ is the SD of the null and alternative hypotheses, n is the sample size (the length of recording in hours), and: Embedded Image Using the significance level α, the cutoff is written as follows: Embedded Image where μ0 is the mean of the null hypothesis (chance-level cross-correlation value estimated by spike-jittering) and: Embedded Image Plugging Equation 4 into Equation 2 yields the following: Embedded Image We conjectured that the probability of detecting a significant cross-correlation between two neurons, if a connection exists and the conditions for its making a contribution to postsynaptic spiking are realized, follows the same form of statistical power. The fitted detection probability of excitatory connections Cαn for the significance level α and the recording length of n hours can be obtained as follows: Embedded Image where CEEα is the asymptotic detection probability of excitatory connections pαn at the significance level α. Finally, to obtain a better fit, the significance level α and the SD σ in the right hand side in Equation 7 were treated as fitting parameters αparam and σparam, respectively. This gives the final fitting equation as follows: Embedded Image In summary, CEEα, αparam, and σparam in Equation 8 were optimized to fit the data of the experimentally derived detection probability of weak E–E connections pαn. This yields the fitted detection probability of excitatory connections Cαn. Note that in a standard testing of a normal mean with a known standard distribution, α and αparam are identical and they represent the significance level of the test. Similarly, σ and σparam are identical and they represent the standard distribution of the null and alternative hypotheses. In this study, however, αparam was separated from α and σparam was separated from σ, and they were treated as fitting parameters. This distinction was necessary to obtain a good fit to our experimental data. We speculate that it might be partly due to the fact that μ1 follows a skewed continuous distribution rather than a normal distribution. Given the fact that many different fitting functions were indeed able to fit our experimental data almost equally well, but would produce greatly different asymptotes, we chose a fitting function that is based on statistical power.

Actual curve fitting was conducted by simulated annealing (Metropolis et al., 1953; Kirkpatrick et al., 1983) using MATLAB Global Optimization Toolbox (MathWorks). Simulated annealing is a well established, general optimization method that has been used as a powerful optimizer for n-body problems including the classic traveling salesman problem (Metropolis et al., 1953; Kirkpatrick et al., 1983). The method can be pictured as the physical process of first heating a material and then lowering the temperature slowly, which corresponds to minimizing the system's energy. More precisely, at a sufficiently high temperature, the system is slightly perturbed and the change of energy, ΔE, is calculated. The new state is accepted if ΔE is negative. If the energy increases, the new state is accepted under a certain probability (p = exp [− ΔE⧸(kBT)], where kB = Boltzmann constant and T = temperature) to avoid that the system gets stuck in local minima. This process is repeated many times at the current temperature. The temperature is slowly decremented until the system is frozen (Kirkpatrick et al., 1983). The method may not find the optimum solution but is most likely to find one of the near-optimum solutions. We ran 1000 optimization trials to select the best parameter values that would give the least sum of squares. See Table 2 for optimization results for all α levels. The r-squares, R2AllEx and R2AllExTT, were computed to estimate the goodness of fit of the optimized function to the data.

Calculation of the effect size, h.

To compare the size of the short-latency peak across neuron pairs, we introduced the effect size h, which was defined as the normalized significant peak as follows: Embedded Image where hpeak is the height of the maximum peak in the cross correlogram in [−4,−1] ms and [1,4] ms (i.e., the maximum peak taken of this 8 ms window) for unidirectional excitatory connections (EE1) and the height of the maximum peaks in the cross correlogram in [−4,−1] ms and [1,4] ms for reciprocal excitatory connections (EE2; i.e., the maximum peak taken of each 4 ms window), and μjitter and σjitter are the mean and SD of the jittered data for the corresponding cross-correlation, respectively (Fujisawa et al., 2008).

Results

Long-term recordings of 25 h length were analyzed using cross-correlations to identify excitatory connections between prefrontal putative pyramidal neurons. Six datasets obtained from three rats that were subjected to both a repeated sequence running task and novel object experience, were analyzed (two novel-object datasets were described previously; Tatsuno et al., 2006). Table 1 shows the number of neuron pairs in each session. In total, we analyzed 237 neurons and 4787 correlation pairs. An example of units included in the analysis is shown in Figure 1a. We classified the neurons into excitatory (54%), inhibitory (8%), and unclassified neurons (38%) based on their overall short-latency effects on postsynaptic neurons (see Materials and Methods; Fig. 1b), as the different classes were not unambiguously differentiable based on their waveform characteristics. This may be due to our recording parameters, e.g., using a rather tight bandpass filter between 600 Hz and 6 kHz, which affects the shape of the recorded spikes.

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

The number of neuron pairs, significant excitatory connections, and the experimentally derived detection probability of excitatory connection pairs pαn in our data counted in three different ways for comparison with intracellular and extracellular recording studies

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

a, Stability of neurons and cleanness of clusters recorded on the same tetrode. Average waveforms (top rows), autocorrelograms (AC, middle left rows), cross correlogram (CC, middle right rows), and time plots (two bottom rows) of two excitatory neurons recorded on the same tetrode are depicted. The percentage of spikes within 2 ms of refractory period included in the cluster and the firing rate are indicated above the autocorrelograms. Both clusters show distinct waveforms (each waveform is depicted for a 1 ms window) and have clean refractory periods (0.02 and 0.03%, respectively, spikes within 2 ms of refractory period). Their cross correlogram is reciprocally significant (EE2) at a significance threshold of α = 10−5 (gray horizontal line). b, Classification of neurons. In the left column, the average Z-scored cross-correlations of three different neurons, without the detrending procedure, are depicted in blue. Right, The corresponding detrended average Z-scored cross-correlations are plotted. The black horizontal dashed line depicts the mean, the red horizontal dashed lines depict one SD above and below the mean, respectively, and the green line is the moving average of the Z-scored cross-correlation computed with a moving window of 15 ms. The vertical black dashed lines represent the limits of the window of interest between 1 and 4 ms, within which short-latency peaks or troughs are assessed for classification. The top row represents an unclassified neuron that, without detrending, would have been classified possibly falsely as excitatory. The middle row shows an inhibitory neuron that without detrending, would have been classified as excitatory. The last row represents an excitatory neuron that would have been left unclassified without detrending.

Changes of detection probability over various lengths of recordings

Significant excitatory connections between neuron pairs were detected by using a spike train jitter method (Barthó et al., 2004; Fujisawa et al., 2008). As the threshold for significant interactions is chosen somewhat arbitrarily by the experimenter and influences the size of the type I error (false-positives), we investigated a wide range of significance levels; 16 different α levels ranging from 0.5 to absolute limit (abs). The absolute limit is derived from the absolute maximum value of all bins in the jittered data for each neuron pair and represents the strictest α level. Note that it is not much smaller than 10−5 (the smallest used before abs).

Figure 2a,b shows four examples of cross correlograms between two putative pyramidal neurons with the absolute limit. Although no significant short-latency peak was detected with a 1 h recording (Fig. 2a,b, top row), it was detected with a 25 h recording (Fig. 2a,b, middle row). The reduced fluctuations in the correlogram reflect the improved signal-to-noise ratio and enable the detection of small peaks. As expected, longer recording times increase the chance of detecting weak excitatory connections. In addition, we also calculated the jitter-corrected cross-correlograms (Fig. 2a,b, bottom row; Hirabayashi et al., 2013; Mizuseki and Buzsaki, 2013). This is obtained by subtracting the mean jittered cross-correlations for each bin from the original cross-correlation. The results show that the significant peaks clearly stand out. The firing rate distribution of the analyzed dataset follows a lognormal distribution as previously shown (Mizuseki and Buzsaki, 2013) and does not change between the first and last 4 h of recordings (Fig. 2c). The distribution of the jitter corrected significant peaks follows a lognormal distribution as well (Fig. 2d), which is also consistent with previous findings (Mizuseki and Buzsaki, 2013).

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

Examples of cross correlograms between two putative pyramidal neurons that are unidirectionally connected (a, EE1 connections) or reciprocally connected (b, EE2 connections).Top, Cross-correlations for 1 h spike trains and in the second row, cross-correlations for the same neurons pairs for 25 h spike trains are shown with a window size of 20 ms and a bin size of 1 ms. The red horizontal lines depict the significance level (α = abs) determined by 1000 spike-jittered cross-correlations. Note the improved signal-to-noise ratio in the 25 h cross correlograms (middle row), which enables the detection of small, short latency peaks. The bottom row shows the jitter corrected correlograms for 25 h obtained by subtracting the mean of the jittered cross correlograms of each cell pair from the original cross correlogram (middle row). c, The distribution of firing rate probability among all neurons for the first 4 h (black) and last 4 h (pink) of recording. The firing rate distribution for both recording intervals is lognormal. Note the log scale on the x-axis. d, The distribution of the jitter corrected effect size for EE1 and EE2 connections. The distribution of the jitter corrected peaks of the significant cross-correlations is lognormal as previously shown by Mizuseki and Buzsaki (2013). The peaks of the neuron pairs shown in a and b are indicated in the plot.

The experimentally derived detection probability of excitatory connections pαn for the significance level α and the recording length of n hours was given by Equation 1. To make the results comparable to paired intracellular recordings between pyramidal neurons (Song et al., 2005), N was taken as the total number of possible excitatory neuron pairs (NAllEx = 1443; Table 1). In addition, because intracellular recordings are typically done in slices and are limited to a smaller volume of tissue, we also calculated N as the number of possible excitatory neuron pairs per tetrode (NAllExTT = 275; Table 1). Although the results for both NAllEx and NAllExTT are presented, we consider that NAllExTT is more directly comparable to previous intracellular recordings (Song et al., 2005). The results of detection probabilities normalized by all recorded neuron pairs (Barthó et al., 2004; Fujisawa et al., 2008; NAll = 4787), including connection types other than excitatory interactions are shown in Figure 9.

Figure 3a (top, data points) shows how the experimentally derived detection probabilities of excitatory connections pαn changes from 1 h recordings to 25 h recordings for four representative α levels: 0.05, 0.01, 10−5, and the absolute limit. The data points represent the median and the background band shows the range for the 84th and 16th percentile estimated by bootstrapping. The detection probabilities increased with longer recording times for both ways of calculating the connection probability. For α = 0.05, the detection probability increased from 6.7% (NAllEx) and 8% (NAllExTT) for 1 h recording segments to 16.4% (NAllEx) and 25% (NAllExTT) for 25 h recording segments, respectively. For α = 0.01, the detection probability increased from 1.6% (NAllEx) and 2.9% (NAllExTT) to 7.2% (NAllEx) and 17.6% (NAllExTT), respectively. For the stricter α values, 10−5 and the absolute limit, which show very similar trends for both ways of normalization, the detection probability increased from 0.1% (NAllEx) and 0.5% (NAllExTT), to 1.7% (NAllEx) and 7.6% (NAllExTT), respectively.

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

a, Top, The experimentally derived detection probability pαn of excitatory connections. The median (circles) and 16th and 84th percentile (shaded area) of the detection probability are plotted for recording lengths of 1–25 h. The data in the left plot was normalized by the total number of possible excitatory neuron pairs (NAllEx) and the data in the right plot by the total number of possible excitatory neuron pairs per tetrode (NAllExTT). The results for four different significance α levels (color coded) are shown: 0.05, 0.01, 10−5, abs (note: the curves for the latter two levels largely overlap). The detection probabilities (Cαn, solid lines) were extrapolated until asymptote (CEEα) as estimated by the fitting of the statistical power function. The gray dashed horizontal line at 11.6% denotes the local connectivity as detected by Song et al. (2005) between pyramidal neurons in layer 5 of rat visual cortex. The blue dashed line indicates the asymptote as approached by α = 10−5. All R2 values indicating the goodness of fit are between 0.98 and 0.99. Bottom, The fitted detection probabilities (Cαn) divided by CEEα are plotted for NAllEx normalization (left) and for NAllExTT normalization (right). b, The asymptotes (CEEα) are plotted against the α levels on a log–log scale. Open squares and asterisks correspond to CEEα normalized by NAllEx and by NAllExTT, respectively. Both ways of CEEα estimation seem to approach a lower limit. Triangles and circles correspond to CEEα normalized by NAllEx and by NAllExTT, respectively, for the alternative window at ±50 ms. CEEα approaches zero as α becomes stricter. This is expected because no monosynaptic interactions are expected to be detected in the time interval of ±50 ms.

Curve fitting and extrapolation of the relationship between statistical power and sample size

Next, we investigated how the detection probability of excitatory connections between putative pyramidal neurons would improve beyond the recording duration of 25 h. The fitting of Equation 8 was applied to the experimentally derived detection probabilities pαn with 16 different α levels (0.5 to the absolute limit) to obtain the fitted detection probabilities Cαn beyond 25 h. The optimization results are summarized in Table 2. Figure 3a (top) show the fitted results for four representative significance thresholds (α levels for 0.05, 0.01, 10−5, and abs). The corresponding power function, defined by Cαn/CEEα, where CEEα is the asymptotic value of pαn, was also obtained (Fig. 3a, bottom). We found that the fitted functions converge to different asymptotic values of the experimentally derived detection probability of excitatory connections pαn but the difference becomes very small for stricter α levels (10−5 and absolute limit). There indeed appears to exist a lower limit asymptote as the threshold becomes stringent (∼2%, NAllEx; ∼8.5%, NAllExTT), indicated by the blue dashed lines in the plots. The CEEα values for the stricter α levels (∼2% with NAllEx and ∼8.5% with NAllExTT) reach ∼20 and 70% of the local connectivity as reported by Song et al. (2005), respectively. As was discussed in the previous section, we consider that NAllExTT is more directly comparable to their study. Therefore, we conclude that ∼70% of local E–E connections could be detected by long extracellular recordings.

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

Optimization results for fitting the statistical power function (Eq. 8) to the experimental data pαn for all investigated α levels

For α levels of 0.05 and 0.01, the power function (Cαn/CEEα) approaches 1 more quickly than that of 10−5 and absolute limit (Fig. 3a, bottom); ∼90% of power could be achieved with 20 h of recording. However, their detection probability also increases quickly and exceeds Song et al.'s (2005) detection probability of 11.6%; α = 0.05 with NAllEx normalization and α = 0.05 and 0.01 with NAllExTT normalization (Fig. 3a, top). Because intracellular recordings have intrinsically a higher chance of detecting monosynaptic connections, as the presynaptic neurons are stimulated and even subthreshold responses in the postsynaptic neurons recorded, the overshoot for α levels 0.05 and 0.01 suggests that these acceptance criteria detect excessive false-positives.

To assess how the same asymptote could be possibly reached, CEEα are plotted against the α levels in percentage on a log–log scale (Fig. 3b). We found that as α levels get stricter, the slope of the consecutive data points for both ways of the CEEα estimation gets increasingly smaller (open square and asterisk correspond to NAllEx and NAllExTT normalization, respectively). This result indicates there is a lower limit asymptote. In contrast, when choosing a window of identical size at longer latency lags (±50 ms), the CEEα estimation approaches zero almost following a diagonal line (open circle and open triangle correspond to NAllEx and NAllExTT normalization, respectively). A similar result was also obtained for the latency lags of ±25 ms. These results demonstrate that the windows at the latency lags outside of monosynaptic interactions (e.g., ±25 ms or ±50 ms) do not contain monosynaptic interactions.

In the final form of the fitting function (Eq. 8), the significance level α was treated as a fitting parameter αparam to obtain a better fit. We confirmed that α and αparam were monotonically related, indicating that fitting αparam does not violate the original relationship of α = αparam significantly (Fig. 4). However, αparam deviates gradually when α gets smaller. This could be due to the fact that μ1 will follow a skewed distribution, which is not accounted for in Equation 8. Thus, Equation 8 would need to be rewritten accordingly but it is beyond the scope of the present paper. Together, these results show that our method using short-latency peaks of cross-correlations with strict α levels and using the approach of fitting a statistical power function is valid for detecting monosynaptic E–E interactions and estimating their asymptotic detection probability in local networks.

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

The relationship between the significance level α and the fitting parameter αparam, obtained from Equation 8. Asterisks depict the results obtained with NAllEx and open squares depict the results obtained with NAllExTT. αparam deviates gradually when α gets smaller. This is likely due to the fact that μ1 is not a single value, but a distribution, which is not accounted for in Equation 8. However, α and αparam maintained a monotonically decreasing relationship. This suggests that separating αparam from α does not violate the original relationship of α = αparam substantially.

False-positives and the distribution of effect size

We investigated whether the connections detected by less stringent α levels (e.g., 0.05 and 0.01) were contaminated by false-positives by plotting the distribution of the effect size h (Eq. 9). The distribution of the effect size within tetrodes (short range) and between tetrodes (long range) for four representative α levels (0.05, 0.01, 10−5, and abs) are plotted in Figure 5a and b, respectively. Both unidirectional (EE1) and reciprocal excitatory connections (EE2) are included in the analysis. The results show that the distributions for the stricter α levels (10−5 and abs) were almost identical and contained strong effects only (Fig. 5, right columns). On the contrary, less strict α levels (0.05 and 0.01) contained additional large counts of smaller effect sizes (Fig. 5, left columns). Subtraction of the most stringent distribution (α = abs) from less stringent groups (α = 10−5, 0.01, and 0.05) confirmed that the stricter α levels (10−5 and abs) are very similar to each other but less stringent α levels (0.05 and 0.01) have higher counts of weak effect sizes (Fig. 5, insets). This is likely due to peak counts being contaminated by a large proportion of false-positives. Because false-positives would occur randomly, their distribution is expected to be Gaussian-like, which is consistent with what we obtained (Fig. 5, insets). In addition, most of the stronger effects are found in short range interactions (within tetrodes; Fig. 5a), not in long-range interactions (between tetrodes; Fig. 5b), suggesting that the strongest connections are mostly local. We also investigated the distribution of the effect size for unidirectional (EE1) and bidirectional connections (EE2) separately; Figure 6a for EE1 (within tetrodes), Figure 6b for EE1 (between tetrodes), Figure 6c for EE2 (within tetrodes), and Figure 6d for EE2 (between tetrodes). We found that the majority of strong effects for both EE1 and EE2 are local and that the long-range excitatory connections are predominantly unidirectional (EE1).

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

The distribution of effect sizes (h) for short-range (a) and long-range (b) connections plotted for four different α values. The distribution of effect size for both unidirectional and reciprocal connections (EE1 and EE2) recorded within tetrodes (a) and between tetrodes (b) for four different α levels (0.05, 0.01, 10−5, and abs) are depicted. The number of peaks included, as well as the medians of the distributions are indicated above each histogram. The median is also plotted as a green vertical line. The results for α levels 10−5 and absolute limit are comparable. The distributions of effect size for α levels 0.05 and 0.01 are also similar. The insets show the difference between the most stringent group (α abs) and the less stringent groups (α levels 0.05 and 0.01). These counts correspond to weak effect sizes (<5) and most likely false-positives. The strong effect sizes (>10) are found in correlations between excitatory neurons recorded on the same tetrode, indicating that the strongest connections are mostly local.

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

The distribution of effect sizes (h) for unidirectional connections (EE1) within (a) and between (b) tetrodes and for reciprocal connections (EE2) within (c) and between (d) tetrodes for α levels 0.05, 0.01, 10−5, and abs. The significance levels as well as the median of the distributions are indicated above the histograms. The insets show the difference between higher α (0.05, 0.01, and 10−5) and α abs, indicating the amount of likely false-positives included in the counts. The strongest effect sizes of unidirectional connections seem to be distributed among short-range connections rather than long-range connections (a). The long-range connections appear to be predominantly single directed excitatory connections (b). As the significance level is determined by a jitter method over a distribution of 1000 jitter trials that can differ for different α levels, it is possible that not all pairs detected with α abs are included in the pairs detected with α 10−5 (even though they are mostly identical). In other words, occasionally a few more pairs can be detected with α abs than with α 10−5, which in the extreme case leads to no peaks detected at α 10−5, but still two peaks detected at α abs (d). Reciprocal excitatory connections seem to be mostly local (c) and according to α levels 10−5, and abs, there are hardly any reciprocal excitatory connections detectable between excitatory neurons recorded on different tetrodes (d). This also suggests that the peaks in α levels 0.01 and 0.05 are likely to be false-positives.

Hyper-reciprocity in mPFC

Further support of the conclusion that less stringent α levels (e.g., 0.01 and 0.05) are massively contaminated by false-positives came from comparing the counts of significant connection pairs at different α levels and the number of connection pairs predicted by a random connection assumption (Fig. 7). Given that a network is randomly connected by a connection probability pαn (Eq. 1), the expected number of unconnected pairs is N(1 − pαn)2, the expected number of unidirectionally connected pairs is 2Npαn (1 − pαn), and the expected number of reciprocally connected pairs is N (pαn)2. Song et al. (2005) reported that the count of reciprocally connected excitatory neuron pairs (EE2) in rat visual cortex is four times higher than the number predicted. Wang et al. (2006) reported 3.5 and 7.9 times more EE2 connections than predicted in visual cortex and prefrontal cortex, respectively, in young ferrets. The ratios of observed EE2 to predicted EE2 are comparable in visual cortex for both rodent studies (4 in rat and 3.5 in ferret). EE2 connections are twice as likely in ferret medial prefrontal cortex than in visual cortex (7.9/3.5 = 2.3) and if the relationship holds true for rat, we would expect the EE2 connection probability in rat medial prefrontal cortex to be up to eight times higher than predicted by a random connectivity assumption. Based on these considerations, we predicted that the EE2 connection in rat medial prefrontal cortex would be four to eight times higher than predicted. For the stringent α levels of the absolute limit and 10−5, we found that medians of observed EE2 connections are 4.6 and 5.6 times higher than the number predicted, respectively (Fig. 7, right columns). In addition, their 95% confidence intervals (Fig. 7, error bars) overlap with the predicted four to eight times range, suggesting that these results are consistent with what we predicted based on the numbers by Song et al. (2005) and Wang et al. (2006). However, for less stringent α levels, the difference between observed EE2 counts and predicted EE2 counts became lower than predicted; 3.2 times for α = 0.01, 2 times for α = 0.05 and 1.3 for α = 0.5 (Fig. 7, left three columns). Furthermore, their 95% confidence intervals do not overlap with the four to eight times range, except for α = 0.01, where there is a small overlap, indicating that the statistics of the observed EE2 connections were different for less stringent α levels. Because false-positives detected in cross correlograms are expected to occur randomly, the observed decrease of hyper-reciprocity suggests that more false-positives are included if less stringent α levels were used. We therefore conclude that stricter α levels together with long-term recordings are necessary for reliable detection of weak E–E interactions in extracellular recordings.

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

The counts of unconnected (NO), single directed excitatory (EE1), and reciprocal excitatory (EE2) neuron pairs for five representative α levels (left to right; 0.5, 0.05, 0.01, 10−5, and abs), relative to the prediction of the random-connection null hypothesis are shown. The connection probability pαn as well as the α levels are indicated above the individual bar plots. The absolute counts of connections are indicated above the bars. Error bars represent the 2.5th and 97.5th percentiles of the medians of 1000 bootstrapped distributions. The dashed horizontal lines at 4 and 8 indicate the expected counts of observed EE2 connections relative to the null hypothesis in the Song et al. (2005) and Wang et al. (2006) studies. For stricter α levels (abs and 10−5), reciprocally connected excitatory neuron pairs (EE2) with the NAllExTT normalization are 4.6–5.6 times more likely than expected. This is in the four to eight times range calculated from Song et al. (2005) and Wang et al. (2006). As α becomes less stringent (0.01, 0.05, and 0.5), the difference becomes less striking because more false-positive counts are likely to be included. The counts of single directed excitatory neuron pairs (EE1) and of unconnected neuron pairs (NO) are close to the expected numbers.

Excitatory connectivity in mPFC is predominantly local

The relationship between local excitatory connectivity (within tetrode connections, normalized by NAllExTT) and total excitatory connectivity (within and between tetrode connections, normalized by NAllEx) can be further investigated by plotting them against each other for different α levels (Fig. 8). For stricter α levels, observed connections were localized within the upper-left triangle, suggesting that predominantly local connections were detected. With increasingly relaxed α levels the difference between local and total connectivity disappeared (at α = 0.2, indicating that the difference was masked by an increasing number of false-positives. Interestingly, total connectivity became larger than local connectivity for α = 0.5. If there was a general tendency of underestimating the local connectivity and overestimating the total connectivity, then this could explain why the total connectivity is larger than the local connectivity for α = 0.5. One could speculate that over elimination of spikes during spike sorting resulting in decreased correlations (Cohen and Kohn, 2011) could lead to the observed effect as this would affect the detection of significant cross-correlations within tetrodes more than between tetrodes. Incapability of recording overlapping spikes within tetrodes may also enhance this tendency.

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

The number of excitatory connections per tetrode (local connectivity, normalized with NAllExTT) is plotted against the total number of excitatory connections (total connectivity, normalized with NAllEx) for all investigated α levels. The red diagonal dashed line indicates equal probability between local and total connectivity (including within and between tetrode connections). The detected connections are predominantly local.

Various types of excitatory connections are the most abundant in mPFC and their detection probability increases with sample size

In addition to E–E interactions (EE1 and EE2), our 25 h recordings also allowed us to investigate six additional connection categories (EI1, EI2, IE1, II1, EU1 and IU1, see Materials and Methods); Although, II2, reciprocal inhibition between putative inhibitory neurons, could be defined as well, we excluded this connection probability from our analysis because its cross-correlation cannot be easily distinguished from the case of common inhibitory input. The detection probabilities of connection categories involving excitatory neurons (EE1, EE2, EI1, EI2, IE1, and EU1) change with increased recording duration for two representative α levels (10−5 and abs; Fig. 9). Note that the normalization by all recorded neuron pairs (NAll) was used because all neuron types (excitatory, inhibitory, and unclassified) were included. We also found that EE1 (Fig. 9b) and EU1 (Fig. 9g) are the most abundant connection types followed by EE2 (Fig. 9c), EI1 (Fig. 9a), and EI2 (Fig. 9d), suggesting that excitatory interactions are the major portion of cortical connections in rat prefrontal cortex. Inhibitory interactions, IE1 (Fig. 9e), II1 (Fig. 9f), and IU1 (Fig. 9 h), are detected much less frequently.

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

Detection probabilities over different hour segments for different categories of connection types. Significance thresholds (α 10−5 and abs) are color coded and indicated at the bottom right of the figure panel. The median of the bootstrapped detection probabilities, normalized by the total number of neuron pairs (4787) including excitatory, inhibitory and unclassified neurons for 1 to 25 h are plotted for the two different significance thresholds for eight different types of connections: (a) unidirectional excitation from an excitatory neuron onto an inhibitory neuron, EI1; (b) unidirectional excitation between two excitatory neurons, EE1; (c) reciprocal excitation between excitatory neurons, EE2; (d) reciprocal connection between an inhibitory and an excitatory neuron, EI2; (e) unidirectional inhibition from an inhibitory neuron onto an excitatory neuron, IE1; (f) unidirectional inhibition between two inhibitory neurons, II1; (g) unidirectional excitation from an excitatory neuron onto an unclassified neuron, EU1; and (h) unidirectional inhibition from an inhibitory neuron onto an unclassified neuron, IU1. The most abundant connection types in mPFC are excitatory connections (a–d, g). Their detection probability increases with longer recording times except for EI1, which seem to be already reliably detected at shorter recording times (a). Inhibitory connections are rarely detected in our dataset and their detection probability does not improve significantly with longer recording times.

Discussion

Pyramidal neurons are the most abundant cell type in neocortex. Pyramid–pyramid connections (E–E interactions) provide the majority of intracortical and extracortical projections. However, most E–E synapses are extremely weak (McNaughton, 1980; Mason et al., 1991; Deuchars et al., 1994; Markram et al., 1997; Thomson and Deuchars, 1997; Reyes and Sakmann, 1999; Thomson et al., 2002). With long-term continuous recordings, we tried to capture weak connections between pyramidal neurons. As opposed to most electrophysiological studies that generally record for ∼0.5–2 h per session, we were able to detect excitatory interactions between putative pyramidal cells in 8.5% of cell pairs with recording lengths up to 25 h. This connection probability of 8.5%, estimated using the number of putative excitatory neuron pairs per tetrode, corresponds to ∼70% of the local E–E connections reported in Song et al. (2005). As was discussed, this normalization is considered to be most comparable to Song et al. (2005). The smaller connection probability (2%) found for the total excitatory connectivity indicates that on average the connection probability falls off with distance (Song et al., 2005; Fujisawa et al., 2008). In addition, 1 h fragments of our recordings yielded a consistent number of monosynaptically coupled pairs (0.2% for strictest α levels: abs and 10−5 and 1.04% for α 0.01; counts were normalized using all recorded neuron pairs, Fig. 9) to previously reported studies (0.2–0.8%; Csicsvari et al., 1998; Barthó et al., 2004; Maurer et al., 2006; Fujisawa et al., 2008).

The strength of a pyramidal–pyramidal connection is partly determined by the location of the synapse and the type of receptors activated (Deuchars et al., 1994). For instance, because lateral connections between pyramidal neurons in cortex are predominantly on distal parts of the dendrites, they are weakened as the potential propagates toward the soma. Those inputs are nearly only transferred if they coincide with other inputs (Deuchars et al., 1994). Therefore the E–E interactions detected by relatively short duration electrophysiological recordings are presumably the strongest ones. However, those are embedded in a much larger matrix of weak connections, which are difficult to capture with extracellular recordings that only detect action potential transmission. With regard to successful spike transmission, pyramidal–interneuron pairs have been shown to be more reliable (Mizumori et al., 1989; Marshall et al., 2002; Holmgren et al., 2003; Swadlow, 2003). Using a less stringent detection threshold, such as 1%, as used in many previous studies, is therefore less likely to be severely affected by false-positives. A higher convergence of excitatory inputs among pyramidal neurons is necessary to bring a postsynaptic pyramidal neuron to firing threshold (McNaughton et al., 1981; Markram et al., 1997). In this respect, long-term recordings dramatically increased the detection probability of the statistically rare E–E events.

In Figure 5, we reported that most connections are local, i.e., both presynaptic and postsynaptic neurons are within the recording radius of a single tetrode. In the hippocampus, monosynaptically coupled pairs of pyramidal neurons and interneurons were also more frequently observed on the same tetrode than on different tetrodes (Csicsvari et al., 1998; Maurer et al., 2006). Connection strength tends to cluster around a few neurons (Barthó et al., 2004; Song et al., 2005; Fujisawa et al., 2008), often referred to as hub neurons. We also observed that many connections converge onto a few neurons. The majority of detected connections are local (Figs. 10a, 11), recorded on the same tetrode (Fig. 10a, red dashed lines). A schematic diagram created from Figure 10a also confirmed that there exists a hub neuron (Fig. 10b, neuron #16). Its inhibitory property (a hub neuron is an inhibitory neuron) is consistent with previous studies (Barthó et al., 2004; Fujisawa et al., 2008).

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

a, The cross-correlogram matrix for all (n = 54) simultaneously recorded neurons in one session (rat 3, sequence task). Red dashed lines indicate correlations between neurons recorded on the same tetrode. Different colors denote different types of connections (see Fig. 9 legend) based on significant short-latency peaks or troughs in the cross correlogram at α abs. (b) Schematic mapping of the connections in the cross correlogram matrix in a between the neurons recorded on five different tetrodes. This animal had a “split bundle” hyperdrive with six recording tetrodes targeting the hippocampus (not included in this paper) and five tetrodes targeting the medial prefrontal cortex (b).

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

Summary of local and long range connections per dataset. The number of significant excitatory connections within (same TT) and between (different TT) tetrodes for the individual datasets for α abs. Numbers are given in absolute (y-axis) and percentage (written on bars) values.

In this study, we showed that the detection probability of weak E–E connections can be significantly improved by long-term recordings (25 h) due to the substantially improved signal-to-noise ratio reflected in smoother cross-correlation signals for 25 h compared with 1 h recordings (Fig. 2). We were also able to detect the hyper-reciprocity of excitatory connections in medial prefrontal cortex with our long-term extracellular recordings (Song et al., 2005; Wang et al., 2006; Fig. 7). This strongly suggests that connectivity is not random, but highly structured. Despite the fact that increased recording times improve the detection of connections among excitatory neurons, extracellular recordings still underestimate the local connectivity as is reported by paired intracellular recordings. At best, they identify ∼70% of local excitatory connections at an α level that is strict enough to exclude most type I errors (Fig. 3). This limitation can be partly attributed to the fact that extracellular recordings in behaving and resting animals are bound to the brain state spaces visited over the duration of recordings, which the experimenter has only limited control over. In contrast, through direct stimulation of the presynaptic neuron, any subthreshold response in the postsynaptic neuron can be monitored in paired intracellular recordings. In addition, our analysis indicates that, to estimate reliably an asymptotic value for the local connectivity, stable recordings of at least 20–25 h, ideally longer, have to be used. This can be a challenge under many experimental conditions. A further caveat of the correlation method is the fact that the nature of the measure is a correlation and therefore reflects synaptic interactions only indirectly. Consequently, the cross-correlation analysis of extracellularly recorded spike trains should not be considered a definitive tool for estimating connectivity and should be used with caution. Despite this caveat, the method can provide important preliminary information for within- and inter-regional circuit connectivity and changes in their connectivity due to experimental manipulations that may not be obtainable with other current methods.

The cross-correlogram, like all other statistical methods, has limitations in terms of its applicability to the interference of causal interactions from spike-train data. First, the method assumes that the spike trains are stationary, meaning that their stochastic properties do not change over time. This assumption is not always easy to justify; for example, with repeated stimulus presentations neurons tend to adapt their responses (Fairhall et al., 2001). Second, cross-correlations are affected by firing rate differences between the neurons (de la Rocha et al., 2007; Amari, 2009). Appropriate normalization, however, can alleviate this problem (Hirase et al., 2001). It has been proposed that the information-geometric measure could be an alternative correlation measure that is statistically orthogonal to the change of firing probability (Amari and Nagaoka, 2000; Amari, 2001). In addition, recent theoretical studies show that an information-geometric measure could be more directly related to synaptic interactions (Tatsuno and Okada, 2004; Tatsuno et al., 2009; Nie and Tatsuno, 2012) and that it can be applied to nonstationary data (Shimazaki et al., 2012). Other promising methods would include a Bayesian approach (Brown et al., 2004; Eldawlatly et al., 2010) and a regularized logistic approach (Zhao et al., 2012). These measures may perform well in identifying causal, nonlinear relationships between neurons. Further investigation of long-term electrophysiological data by such statistical methods would promote the detection of neural interaction and hence contribute to the understanding of the circuit dynamics underlying complex behavior.

Footnotes

  • This work was funded by an Alberta Innovates Health Solutions (AIHS) Polaris Award and National Institute of Mental Health R01MH46823-16 Grant awarded to B.L.M., an Innovates Centre of Research Excellence SCH001 Grant awarded to M.T., and an AIHS graduate Fellowship awarded to C.D.S. We thank Francesco Battaglia for helpful comments on a previous version of the paper.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Masami Tatsuno, Department of Neuroscience, Canadian Centre for Behavioural Neuroscience, The University of Lethbridge, 4401 University Drive W, Lethbridge, AB T1K 3M4, Canada. tatsuno{at}uleth.ca

This article is freely available online through the J Neurosci Author Open Choice option.

References

  1. ↵
    1. Aertsen AM,
    2. Gerstein GL
    (1985) Evaluation of neuronal connectivity: sensitivity of cross-correlation. Brain Res 340:341–354, doi:10.1016/0006-8993(85)90931-X, pmid:4027655.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Alonso JM,
    2. Martinez LM
    (1998) Functional connectivity between simple cells and complex cells in cat striate cortex. Nat Neurosci 1:395–403, doi:10.1038/1609, pmid:10196530.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Amari S
    (2009) Measure of correlation orthogonal to change in firing rate. Neural Comput 21:960–972, doi:10.1162/neco.2008.03-08-729, pmid:19018699.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Amari Si
    (2001) Information geometry on hierarchy of probability distributions. IEEE Trans Inform Theory 47:1701–1711, doi:10.1109/18.930911.
    OpenUrlCrossRef
  5. ↵
    1. Amari S,
    2. Nagaoka H
    (2000) Methods of information geometry (American Mathematical Society, Providence, RI).
  6. ↵
    1. Barthó P,
    2. Hirase H,
    3. Monconduit L,
    4. Zugaro M,
    5. Harris KD,
    6. Buzsáki G
    (2004) Characterization of neocortical principal cells and interneurons by network interactions and extracellular features. J Neurophysiol 92:600–608, doi:10.1152/jn.01170.2003, pmid:15056678.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Bock DD,
    2. Lee WC,
    3. Kerlin AM,
    4. Andermann ML,
    5. Hood G,
    6. Wetzel AW,
    7. Yurgenson S,
    8. Soucy ER,
    9. Kim HS,
    10. Reid RC
    (2011) Network anatomy and in vivo physiology of visual cortical neurons. Nature 471:177–182, doi:10.1038/nature09802, pmid:21390124.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Brown EN,
    2. Kass RE,
    3. Mitra PP
    (2004) Multiple neural spike train data analysis: state-of-the-art and future challenges. Nat Neurosci 7:456–461, doi:10.1038/nn1228, pmid:15114358.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Buzsáki G
    (2004) Large-scale recording of neuronal ensembles. Nat Neurosci 7:446–451, doi:10.1038/nn1233, pmid:15114356.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Cohen MR,
    2. Kohn A
    (2011) Measuring and interpreting neuronal correlations. Nat Neurosci 14:811–819, doi:10.1038/nn.2842, pmid:21709677.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Csicsvari J,
    2. Hirase H,
    3. Czurko A,
    4. Buzsáki G
    (1998) Reliability and state dependence of pyramidal cell-interneuron synapses in the hippocampus: an ensemble approach in the behaving rat. Neuron 21:179–189, doi:10.1016/S0896-6273(00)80525-5, pmid:9697862.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Csicsvari J,
    2. Henze DA,
    3. Jamieson B,
    4. Harris KD,
    5. Sirota A,
    6. Barthó P,
    7. Wise KD,
    8. Buzsáki G
    (2003) Massively parallel recording of unit and local field potentials with silicon-based electrodes. J Neurophysiol 90:1314–1323, doi:10.1152/jn.00116.2003, pmid:12904510.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. de la Rocha J,
    2. Doiron B,
    3. Shea-Brown E,
    4. Josić K,
    5. Reyes A
    (2007) Correlation between neural spike trains increases with firing rate. Nature 448:802–806, doi:10.1038/nature06028, pmid:17700699.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Deuchars J,
    2. West DC,
    3. Thomson AM
    (1994) Relationships between morphology and physiology of pyramid-pyramid single axon connections in rat neocortex in vitro. J Physiol 478:423–435, pmid:7965856.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Eldawlatly S,
    2. Zhou Y,
    3. Jin R,
    4. Oweiss KG
    (2010) On the use of dynamic Bayesian networks in reconstructing functional neuronal networks from spike train ensembles. Neural Comput 22:158–189, doi:10.1162/neco.2009.11-08-900, pmid:19852619.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Euston DR,
    2. Tatsuno M,
    3. McNaughton BL
    (2007) Fast-forward playback of recent memory sequences in prefrontal cortex during sleep. Science 318:1147–1150, doi:10.1126/science.1148979, pmid:18006749.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Fairhall AL,
    2. Lewen GD,
    3. Bialek W,
    4. de Ruyter Van Steveninck RR
    (2001) Efficiency and ambiguity in an adaptive neural code. Nature 412:787–792, doi:10.1038/35090500, pmid:11518957.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Fujisawa S,
    2. Amarasingham A,
    3. Harrison MT,
    4. Buzsáki G
    (2008) Behavior-dependent short-term assembly dynamics in the medial prefrontal cortex. Nat Neurosci 11:823–833, doi:10.1038/nn.2134, pmid:18516033.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Gothard KM,
    2. Skaggs WE,
    3. McNaughton BL
    (1996) Dynamics of mismatch correction in the hippocampal ensemble code for space: interaction between path integration and environmental cues. J Neurosci 16:8027–8040, pmid:8987829.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Hirabayashi T,
    2. Takeuchi D,
    3. Tamura K,
    4. Miyashita Y
    (2013) Microcircuits for hierarchical elaboration of object coding across primate temporal areas. Science 341:191–195, doi:10.1126/science.1236927, pmid:23846902.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Hirase H,
    2. Leinekugel X,
    3. Csicsvari J,
    4. Czurkó A,
    5. Buzsáki G
    (2001) Behavior-dependent states of the hippocampal network affect functional clustering of neurons. J Neurosci 21:RC145, pmid:11319243.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Hoffman KL,
    2. McNaughton BL
    (2002) Coordinated reactivation of distributed memory traces in primate neocortex. Science 297:2070–2073, doi:10.1126/science.1073538, pmid:12242447.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Holmgren C,
    2. Harkany T,
    3. Svennenfors B,
    4. Zilberter Y
    (2003) Pyramidal cell communication within local networks in layer 2/3 of rat neocortex. J Physiol 551:139–153, doi:10.1113/jphysiol.2003.044784, pmid:12813147.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Kirkpatrick S,
    2. Gelatt CD Jr.,
    3. Vecchi MP
    (1983) Optimization by simulated annealing. Science 220:671–680, doi:10.1126/science.220.4598.671, pmid:17813860.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Kirkwood PA
    (1979) On the use and interpretation of cross-correlations measurements in the mammalian central nervous system. J Neurosci Methods 1:107–132, doi:10.1016/0165-0270(79)90009-8, pmid:161793.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Ko H,
    2. Hofer SB,
    3. Pichler B,
    4. Buchanan KA,
    5. Sjöström PJ,
    6. Mrsic-Flogel TD
    (2011) Functional specificity of local synaptic connections in neocortical networks. Nature 473:87–91, doi:10.1038/nature09880, pmid:21478872.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Léger JF,
    2. Stern EA,
    3. Aertsen A,
    4. Heck D
    (2005) Synaptic integration in rat frontal cortex shaped by network activity. J Neurophysiol 93:281–293, doi:10.1152/jn.00067.2003, pmid:15306631.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Markram H,
    2. Lübke J,
    3. Frotscher M,
    4. Roth A,
    5. Sakmann B
    (1997) Physiology and anatomy of synaptic connections between thick tufted pyramidal neurones in the developing rat neocortex. J Physiol 500:409–440, pmid:9147328.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Marshall L,
    2. Henze DA,
    3. Hirase H,
    4. Leinekugel X,
    5. Dragoi G,
    6. Buzsáki G
    (2002) Hippocampal pyramidal cell-interneuron spike transmission is frequency dependent and responsible for place modulation of interneuron discharge. J Neurosci 22:RC197, pmid:11784809.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Mason A,
    2. Nicoll A,
    3. Stratford K
    (1991) Synaptic transmission between individual pyramidal neurons of the rat visual cortex in vitro. J Neurosci 11:72–84, pmid:1846012.
    OpenUrlAbstract
  31. ↵
    1. Maurer AP,
    2. Cowen SL,
    3. Burke SN,
    4. Barnes CA,
    5. McNaughton BL
    (2006) Phase precession in hippocampal interneurons showing strong functional coupling to individual pyramidal cells. J Neurosci 26:13485–13492, doi:10.1523/JNEUROSCI.2882-06.2006, pmid:17192431.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. McNaughton BL
    (1980) Evidence for two physiologically distinct perforant pathways to the fascia dentata. Brain Res 199:1–19, doi:10.1016/0006-8993(80)90226-7, pmid:7407615.
    OpenUrlCrossRefPubMed
  33. ↵
    1. McNaughton BL,
    2. Barnes CA,
    3. Andersen P
    (1981) Synaptic efficacy and EPSP summation in granule cells of rat fascia dentata studied in vitro. J Neurophysiol 46:952–966, pmid:7299453.
    OpenUrlFREE Full Text
  34. ↵
    1. McNaughton BL,
    2. O'Keefe J,
    3. Barnes CA
    (1983) The stereotrode: a new technique for simultaneous isolation of several single units in the central nervous system from multiple unit records. J Neurosci Methods 8:391–397, doi:10.1016/0165-0270(83)90097-3, pmid:6621101.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Metropolis N,
    2. Rosenbluth AW,
    3. Rosenbluth MN,
    4. Teller AH
    (1953) Equation of state calculations by fast computing machines. J Chem Phys 21:1087–1092, doi:10.1063/1.1699114.
    OpenUrlCrossRef
  36. ↵
    1. Mizumori SJ,
    2. McNaughton BL,
    3. Barnes CA
    (1989) A comparison of supramammillary and medial septal influences on hippocampal field potentials and single-unit activity. J Neurophysiol 61:15–31, pmid:2493075.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Mizuseki K,
    2. Buzsáki G
    (2013) Preconfigured, skewed distribution of firing rates in the hippocampus and entorhinal cortex. Cell Rep 4:1010–1021, doi:10.1016/j.celrep.2013.07.039, pmid:23994479.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Mooney CZ,
    2. Duval RD
    (1993) Bootstrapping: a non-parametrical approach to statistical inference (Sage Publications, Newbury Park, CA).
  39. ↵
    1. Nie Y,
    2. Tatsuno M
    (2012) Information-geometric measures for estimation of connection weight under correlated inputs. Neural Comput 24:3213–3245, doi:10.1162/NECO_a_00367, pmid:22970877.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Ostojic S,
    2. Brunel N,
    3. Hakim V
    (2009) How connectivity, background activity, and synaptic properties shape the cross-correlation between spike trains. J Neurosci 29:10234–10253, doi:10.1523/JNEUROSCI.1275-09.2009, pmid:19692598.
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Perkel DH,
    2. Gerstein GL,
    3. Moore GP
    (1967a) Neuronal spike trains and stochastic point processes: I. The single spike train. Biophys J 7:391–418, doi:10.1016/S0006-3495(67)86596-2, pmid:4292791.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Perkel DH,
    2. Gerstein GL,
    3. Moore GP
    (1967b) Neuronal spike trains and stochastic point processes: II. Simultaneous spike trains. Biophys J 7:419–440, doi:10.1016/S0006-3495(67)86597-4, pmid:4292792.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Reyes A,
    2. Sakmann B
    (1999) Developmental switch in the short-term modification of unitary EPSPs evoked in layer 2/3 and layer 5 pyramidal neurons of rat neocortex. J Neurosci 19:3827–3835, pmid:10234015.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Shimazaki H,
    2. Amari S,
    3. Brown EN,
    4. Grün S
    (2012) State-space analysis of time-varying higher-order spike correlation for multiple neural spike train data. PLoS Comput Biol 8:e1002385, doi:10.1371/journal.pcbi.1002385, pmid:22412358.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Song S,
    2. Sjöström PJ,
    3. Reigl M,
    4. Nelson S,
    5. Chklovskii DB
    (2005) Highly nonrandom features of synaptic connectivity in local cortical circuits. PLoS Biol 3:e68, doi:10.1371/journal.pbio.0030068, pmid:15737062.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Swadlow HA
    (2003) Fast-spike interneurons and feedforward inhibition in awake sensory neocortex. Cereb Cortex 13:25–32, doi:10.1093/cercor/13.1.25, pmid:12466212.
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Tatsuno M,
    2. Okada M
    (2004) Investigation of possible neural architectures underlying information-geometric measures. Neural Comput 16:737–765, doi:10.1162/089976604322860686, pmid:15025828.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Tatsuno M,
    2. Lipa P,
    3. McNaughton BL
    (2006) Methodological considerations on the use of template matching to study long-lasting memory trace replay. J Neurosci 26:10727–10742, doi:10.1523/JNEUROSCI.3317-06.2006, pmid:17050712.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Tatsuno M,
    2. Fellous JM,
    3. Amari S
    (2009) Information-geometric measures as robust estimators of connection strengths and external inputs. Neural Comput 21:2309–2335, doi:10.1162/neco.2009.04-08-748, pmid:19538092.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Thomson AM,
    2. Deuchars J
    (1997) Synaptic interactions in neocortical local circuits: dual intracellular recordings in vitro. Cereb Cortex 7:510–522, doi:10.1093/cercor/7.6.510, pmid:9276176.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Thomson AM,
    2. West DC,
    3. Wang Y,
    4. Bannister AP
    (2002) Synaptic connections and small circuits involving excitatory and inhibitory neurons in layers 2–5 of adult rat and cat neocortex: triple intracellular recordings and biocytin labelling in vitro. Cereb Cortex 12:936–953, doi:10.1093/cercor/12.9.936, pmid:12183393.
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Wang Y,
    2. Markram H,
    3. Goodman PH,
    4. Berger TK,
    5. Ma J,
    6. Goldman-Rakic PS
    (2006) Heterogeneity in the pyramidal network of the medial prefrontal cortex. Nat Neurosci 9:534–542, doi:10.1038/nn1670, pmid:16547512.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Zhao M,
    2. Batista A,
    3. Cunningham JP,
    4. Chestek C,
    5. Rivera-Alvidrez Z,
    6. Kalmar R,
    7. Ryu S,
    8. Shenoy K,
    9. Iyengar S
    (2012) An L(1)-regularized logistic model for detecting short-term neuronal interactions. J Comput Neurosci 32:479–497, doi:10.1007/s10827-011-0365-5, pmid:22038503.
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 34 (16)
Journal of Neuroscience
Vol. 34, Issue 16
16 Apr 2014
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Long-Term Recordings Improve the Detection of Weak Excitatory–Excitatory Connections in Rat Prefrontal Cortex
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Long-Term Recordings Improve the Detection of Weak Excitatory–Excitatory Connections in Rat Prefrontal Cortex
C. Daniela Schwindel, Karim Ali, Bruce L. McNaughton, Masami Tatsuno
Journal of Neuroscience 16 April 2014, 34 (16) 5454-5467; DOI: 10.1523/JNEUROSCI.4350-13.2014

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Long-Term Recordings Improve the Detection of Weak Excitatory–Excitatory Connections in Rat Prefrontal Cortex
C. Daniela Schwindel, Karim Ali, Bruce L. McNaughton, Masami Tatsuno
Journal of Neuroscience 16 April 2014, 34 (16) 5454-5467; DOI: 10.1523/JNEUROSCI.4350-13.2014
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • cross correlations
  • extracellular recordings
  • local connectivity
  • medial prefrontal cortex

Responses to this article

Respond to this article

Jump to comment:

  • Re:Shadowing effect?
    Masami Tatsuno
    Published on: 29 June 2014
  • Shadowing effect?
    Jianing Yu
    Published on: 17 April 2014
  • Published on: (29 June 2014)
    Page navigation anchor for Re:Shadowing effect?
    Re:Shadowing effect?
    • Masami Tatsuno, Faculty
    • Other Contributors:
      • C. Daniela Schwindel, Karim Ali, Bruce L. McNaughton

    I thank Dr. Jianing Yu for the comment on our article. We have carefully read the papers about the shadowing effect. Unlike the very high firing rate neurons in the globus pallidus that Bar-Gad and his colleagues recorded, our cortical neurons have relatively low firing rate. This indicates that the shadowing effect would be less severe in our study.

    We have identified eight EE2 neuron pairs. Among them, seven n...

    Show More

    I thank Dr. Jianing Yu for the comment on our article. We have carefully read the papers about the shadowing effect. Unlike the very high firing rate neurons in the globus pallidus that Bar-Gad and his colleagues recorded, our cortical neurons have relatively low firing rate. This indicates that the shadowing effect would be less severe in our study.

    We have identified eight EE2 neuron pairs. Among them, seven neuron pairs were recorded from the same tetrodes. To assess the influence of the shadowing effect on those seven EE2 interactions, we have recalculated the cross-correlogram with the shadowing effect taken into account. For each of the shadowing windows of 0.5 ms, 1.0 ms and 1.5 ms, 1000 shuffled cross -correlograms were calculated from spike-jittered neuron pairs. Although the size of the shadowing window of the Cheetah system was 1.0 ms, we performed the calculation up to 1.5 ms.

    We did observe a mild shadowing effect in all of the seven data sets. The amplitude of the shadowing effect increased as the window size was increased from 0.5 ms to 1.5 ms. However, we also confirmed that the peaks of the original EE2 connection within plus/minus [1-4] ms bins were always higher than the maximum values of the shadowing effect in the corresponding bins. There was one case in which a rather large shadowing effect arose but the artifacts were outside of the plus/minus [1-4] ms range.

    Taken together, we concluded that the shadowing effect did not affect the detection of EE2 connectivity in our study.

    We appreciate the important and informative feedback. We will take the shadowing effect into account when we apply the cross-correlation analysis in the future.

    Masami Tatsuno (Corresponding author)

    Conflict of Interest:

    None declared

    Show Less
    Competing Interests: None declared.
  • Published on: (17 April 2014)
    Page navigation anchor for Shadowing effect?
    Shadowing effect?
    • Jianing Yu

    The authors might want to study the two papers listed below, which discuss the shadowing effect: the distortion of cross-correlation if units were from the same electrode. In short, overlapping spikes cannot be resolved for both cells, thereby producing an artificial decrease of correlation near the zero-lag of the cross-correlograms. The peaks surrounding the decrease reflect a gradual decrease of such inability. In additi...

    Show More

    The authors might want to study the two papers listed below, which discuss the shadowing effect: the distortion of cross-correlation if units were from the same electrode. In short, overlapping spikes cannot be resolved for both cells, thereby producing an artificial decrease of correlation near the zero-lag of the cross-correlograms. The peaks surrounding the decrease reflect a gradual decrease of such inability. In addition, synchronous firing of cortical cells could enhance the peaks surrounding zero-lag.

    In the current paper, example data that were supposed to show "connections" look extremely similar to what Bar-Gad and colleagues believed to be an artifact (compare Fig. 1, 2 in Schwindel et al., to Fig. 2 in Bar-Gad 2003). As Schwindel et al. noticed, such "connections" were most likely to occur when units were from the same electrode/tetrode (Fig. 10, 11 ). I am wondering what the authors' thoughts are on this matter.

    References

    1. Functional correlations between neighboring neurons in the primate globus pallidus are weak or nonexistent. Bar-Gad I, Heimer G, Ritov Y, Bergman H. J Neurosci. 2003 May 15;23(10):4012-6.

    2. Failure in identification of overlapping spikes from multiple neuron activity causes artificial correlations. Bar-Gad I, Ritov Y, Vaadia E, Bergman H. J Neurosci Methods. 2001 May 30;107(1-2):1-13.

    Conflict of Interest:

    None declared

    Show Less
    Competing Interests: None declared.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Memory Retrieval Has a Dynamic Influence on the Maintenance Mechanisms That Are Sensitive to ζ-Inhibitory Peptide (ZIP)
  • Neurophysiological Evidence for a Cortical Contribution to the Wakefulness-Related Drive to Breathe Explaining Hypocapnia-Resistant Ventilation in Humans
  • Monomeric Alpha-Synuclein Exerts a Physiological Role on Brain ATP Synthase
Show more Articles

Systems/Circuits

  • Exploring Neural Dynamics in the Auditory Telencephalon of Crows Using Functional Ultrasound Imaging
  • Oddball Evoked Deviant Responses Reflect Complex Context-Dependent Expectations in Mouse V1
  • V2b Neurons Act via Multiple Targets to Produce in Phase Inhibition during Locomotion
Show more Systems/Circuits
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Accessibility
(JNeurosci logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.