REVIEW article

Front. Integr. Neurosci., 27 November 2018

Volume 12 - 2018 | https://doi.org/10.3389/fnint.2018.00056

Population Coding and Correlated Variability in Electrosensory Pathways

  • Department of Physiology, McGill University, Montréal, QC, Canada

Abstract

The fact that perception and behavior depend on the simultaneous and coordinated activity of neural populations is well established. Understanding encoding through neuronal population activity is however complicated by the statistical dependencies between the activities of neurons, which can be present in terms of both their mean (signal correlations) and their response variability (noise correlations). Here, we review the state of knowledge regarding population coding and the influence of correlated variability in the electrosensory pathways of the weakly electric fish Apteronotus leptorhynchus. We summarize known population coding strategies at the peripheral level, which are largely unaffected by noise correlations. We then move on to the hindbrain, where existing data from the electrosensory lateral line lobe (ELL) shows the presence of noise correlations. We summarize the current knowledge regarding the mechanistic origins of noise correlations and known mechanisms of stimulus dependent correlation shaping in ELL. We finish by considering future directions for understanding population coding in the electrosensory pathways of weakly electric fish, highlighting the benefits of this model system for understanding the origins and impact of noise correlations on population coding.

Interpretation of Population Activity Requires Considering Correlations Between the Activities of Different Neurons

Understanding the concerted activity of neural populations remains a central problem in systems neuroscience. While simultaneous recordings of multiple neurons (i.e., population activity) has become increasingly feasible across animal models and brain areas, interpreting these data are often complicated. This is because neuronal activities are often not independent of one another, but rather show correlations. Such correlations have been found almost ubiquitously across species and brain areas (for review see ; ), and it has been shown that they can have substantial impact, of either detrimental or beneficial nature, on signal encoding performance at the population level (; Salinas and Sejnowski, 2011). Albeit their acknowledged importance (for review see Nirenberg and Latham, 2003; ; ; Salinas and Sejnowski, 2011; Rothschild and Mizrahi, 2015; ),the origins and the functional implications of correlations for coding remain poorly understood in many cases.

Here we review the current state of knowledge regarding how electrosensory neural populations encode behaviorally relevant stimuli in wave-type weakly electric fish. This model system is well described in terms of its anatomy and single cell physiology and readily accessible for in vivo recordings. As such, this system benefits from unique advantages when, e.g., considering realistic decoding of population activity and how this leads to behavior as well as for understanding population coding in more naturalistic experimental paradigms (e.g., when recording from unrestrained animals during active exploration and sensing).

Decomposing Correlations

When evaluating the statistical dependency between the spike trains of two neurons, their correlations (raw-correlation) can be decomposed into two types of correlations (Perkel et al., 1967). First, signal correlations, which are correlations between the mean activity of two neurons responding to stimuli. Second, noise correlations, which are correlations between the trial-to-trial variabilities of the neural responses to repeated presentations of a given stimulus (we will use the terms “noise correlations” and “correlated variability” as synonyms throughout this manuscript). It is important to note that simultaneous recordings are required to infer noise but not signal correlations. Theory predicts that the correlation structure (i.e., the relationship between signal and noise correlations) will determine their impact on information transmission (; ). To exemplify this, let us consider the simultaneous responses of a pair of neurons (Figure 1) to repeated presentations of two stimuli (black and light gray dots in Figures 1A,B). For each of the stimuli, the responses of the two neurons show variability and scatter around their mean (gray areas show the 95% probability distributions, white crosses indicate the means).

FIGURE 1

).

The neuron pair shown in Figure 1A has an opposing correlation structure: the average responses of the two neurons are positively correlated (i.e., both increase their mean response when stimulus 2 is presented vs. stimulus 1) and as such their signal correlations are positive (see insets, red curves). If the variabilities of the two neuron responses to repeated presentations of the same stimulus were independent, the probability distributions around the means would be circular in shape. Instead, they have an elliptical shape with the main axis being oriented from top left to bottom right. As such, the response variabilities are not independent but rather correlated. Indeed, whenever the response of neuron 1 is higher than its mean response, the response of neuron 2 tends to be lower than its mean response and vice versa. Thus the variabilities are negatively correlated, and noise correlations are negative (see insets, blue curves). This opposing correlation structure (positive signal and negative noise correlations, see orientation of red and blue arrows) is beneficial for stimulus encoding and by using a decision criterion (dotted line) it is possible to perfectly discriminate between population responses to the two different stimuli.

In Figure 1B, we show a different pair of neurons with a parallel correlation structure (i.e., positive signal and positive noise correlations). In this case, the distributions of the neural variabilities have an elliptical shape whose main axis is oriented from bottom left to upper right, which is parallel to the changes seen in mean responses to the different stimuli (white crosses, see orientation of red and blue arrows). As a result, the two response distributions show considerable overlap and the decision criterion (dotted line) results in impaired performance compared to the opposing correlation structure discussed above. This exemplifies that the presence of noise correlations can either be detrimental or beneficial for stimulus encoding and that their effect needs to be evaluated on a case by case basis.

The insets in Figures 1A,B show examples of how correlations are typically quantified. Cross correlograms (CCGs, inset top) quantify the number of coincident events per unit time relative to chance as a function of lag (i.e., the amount by which a spike train is shifted relative to the other). The integration and normalization of such CCGs produces a correlation coefficient that quantifies the correlation at infinite timescale (Shadlen and Newsome, 1998). Recently, however, the use of spike count correlations (inset bottom) has become more and more common. For this, the spike trains are binned into time windows of a defined width “t,” and the number of spikes falling within each bin is counted. The resultant spike count timeseries are compared by calculating the Pearson’s correlation coefficient. By reiterating the analysis with different width of the spike count window “t,” correlations can be analyzed at different timescales. Signal and noise correlations can be obtained from spike trains using standard computational methodology such as shuffle predictor and computing response residuals (Perkel et al., 1967).

In the absence of stimulation, raw-correlations between neural activities are often termed “baseline correlations.” These baseline correlations represent the limit that noise correlations will tend toward as stimulus amplitude goes to zero (). Therefore, it is expected that the presence of baseline correlations predicts the presence of noise correlation under stimulation.

Effects of Correlated Variability on Stimulus Encoding

One could argue that the detrimental effects of noise correlations toward stimulus discrimination in the example shown (Figure 1B) might seem minimal. This is because most of the responses will still be categorized correctly based on the decision criterion (Figure 1B) and only a minor part of the responses will be attributed to the wrong stimulus. It should, however, be noted that the shown example considers only two neurons and that perception is typically determined by integrating the activities of much larger neural populations. It was shown theoretically that small pairwise correlations can have strong effects on signal encoding when large neural populations are considered (Zohary et al., 1994; ; Schneidman et al., 2006) (Figure 1C). In the absence of noise correlations (i.e., when trial-to-trial variabilities are independent), the amount of information represented by the population activity grows as more neurons are considered for analysis (Figure 1C, solid line). The effects of trial-to-trial variability on the population information highly depend on the correlation structure. With an opposing structure (top dotted line; noise and signal correlations have an opposing sign), the growth of information quickly surpasses the independent case. In contrast, a parallel correlation structure (bottom dotted line; noise and signal correlations have the same sign) will lead to a reduction in information growth and quick saturation (Zohary et al., 1994). Thus, the presence of noise correlations alone does not suffice to assess their impact on signal encoding as they could be either detrimental (Zohary et al., 1994; Moreno-Bote et al., 2014) or beneficial (; Romo et al., 2003) and their actual effect highly depends on the correlation structure as well as on the subsequent decoding by downstream brain areas. In that regard, various different decoders can be used (Pouget et al., 2000) and many studies have assumed linear decoders (i.e., relevant quantities are estimated based on weighted linear sums of neuronal responses) (Seung and Sompolinsky, 1993; Zohary et al., 1994; Sanger, 1996; ; ; Pitkow et al., 2015). Such decoders are attractive because they are easy to implement and optimize on neural data to quantify the effects of correlations. However, knowledge gained from these require comparison to physiologically plausible decoding strategies which are, in general, nonlinear and thus can in theory extract much more information than linear decoding strategies (Shamir and Sompolinsky, 2004). The effects of various decoders on determining how correlations influence information transmission has been reviewed in detail elsewhere ().

Recent studies have shown that noise correlations are not static but can change in magnitude based on various factors such as the animal’s state (Poulet and Petersen, 2008; ; ; Vinck et al., 2015), the animal’s attention (Steinmetz et al., 2000; ), adaptation to stimuli (), or in a stimulus-dependent fashion (; Snyder et al., 2014; Tan et al., 2014; ; Zylberberg et al., 2016). The plasticity of noise correlations greatly complicates understanding their effect on information coding.

The Mechanistic Origins of Correlations in vivo Remain Poorly Understood

While correlations, both signal and noise, are found ubiquitously in the CNS, in many cases understanding their mechanistic origins remains elusive. On the one hand, it is clear that signal and noise correlations in pairs of neurons will highly depend on their input connectivity (Figure 1D). Signal correlations arise when two neurons receive inputs that encode the same signal and are also tuned to this signal (i.e., they both respond to this signal). Noise correlations are generally thought to arise because of shared neuronal input (Shadlen and Newsome, 1998; Renart et al., 2010; ). The activity of these common inputs (Figure 1D, solid lines) will influence the membrane potential of the receiving neurons similarly and thus introduce noise correlations between their spiking activities. Independent inputs (dotted lines) will usually dilute that stochasticity, thereby decreasing noise correlations potentially. Thus, noise correlations are likely to depend on the network architecture and activity (; ). For example, it was shown that the balance between excitatory and inhibitory inputs might be one key determinant of correlation magnitude (Renart et al., 2010; ).

There has been progress made toward understanding the mechanistic origins of neural correlations and their effects on coding, and it has become clear that such understanding will require detailed knowledge of the anatomical connections. Nonetheless, the connections of neural networks forming circuits in the brain and their inputs and outputs are often numerous, divers, and highly complicated. Toward this end, the electrosensory pathway benefits from a relatively simple and well-characterized anatomy, which should be advantageous for the investigation of the basic mechanisms underlying correlated variability.

Electrosensory Stimuli: Electrolocation and Social Interaction in Apteronotus

Wave-type weakly electric fish such as A. leptorhynchus emit a quasi-sinusoidal electric signal referred to as the “electric organ discharge” (EOD) thereby building up an oscillatory field surrounding their body (Figure 2A, top). It is important to realize that this ongoing EOD acts as a carrier signal during active sensing and that it is the perturbations of the EOD that carry information about the sensory environment. These are picked up by a distributed array of electroreceptors in the skin of the animal (). The animal can detect both amplitude modulations (AMs) as well as frequency modulations (FMs) of the EOD via separate pathways. In the following, we will focus on the AM coding pathway and will henceforth refer to the AM as the stimulus.

FIGURE 2

Electrosensory stimuli occur in different behavioral contexts. During prey capture (Nelson and MacIver, 1999), animals detect and localize (i.e., “electrolocation”) small prey items that cause weak and spatially localized AMs of the EOD (Figure 2A). Several studies have shown that the resulting pattern of stimulation carries important information about the distance, size, and conductivity of an object and the relative speed and angle of the motion between object and fish (Rasnow, 1996; Nelson and MacIver, 1999; Nelson et al., 2002; ; ; Pedraja et al., 2018).

Another behavioral context is that of interactions with conspecific fish (i.e., “electrocommunication”) (Ramcharitar et al., 2005; ; ). When two individuals are in close vicinity to one another, interaction between their EODs will create a sinusoidal stimulus (i.e., a beat) whose frequency is equal to the EOD frequency difference and ranges between a few Hz to several hundred Hz (Figure 2B). It is important to note that such stimuli are spatially diffuse and extend to most if not all the electroreceptors.

During social interaction, fish can emit short-term alterations of their EOD frequency with the purpose of active social communication (; Zupanc and Maler, 1993; Zakon et al., 2002; ; Zupanc et al., 2006; ). Such events are called chirps and always occur on top of the beat (Figure 2C). There are different types of chirps (Zakon et al., 2002) and in the following we will focus on so called “small chirps,” which are typically aggressive call signals. Chirps are produced through brief (ΔT: 13–16 ms) and small (ΔF: 30–50 Hz) excursions in the EOD frequency of one fish (Figure 2C, middle). As a result, the periodic signature of the AM is interrupted by a high-frequency transient that resets the phase of the AM (Figure 2C, bottom). Importantly, for a given chirp with a fixed ΔF and ΔT, the exact waveform of the chirp will look very different depending on the AM phase at which the chirp is emitted. Chirps will cause diverse responses in pyramidal cells (for review see ), but different chirp waveforms of the same chirp will give rise to similar behavioral responses (; Metzen et al., 2016a).

So far, we have only considered stimuli that consist of changes in the mean EOD amplitude. These are sometimes referred to as “first-order” stimuli. However, it is clear that rather than being stationary, fish move extensively during social interactions, thereby causing changes in the amplitude of the beat stimuli (). Such “second-order” stimuli have been termed “movement envelopes” (Yu et al., 2012; Stamper et al., 2013; Metzen and Chacron, 2014). As an example, let us consider one fish looming toward a conspecific (Figure 2D, top). As seen from the perspective of the moving fish, the EOD amplitude of its conspecific will grow during the looming motion, thereby causing an increase in the beat amplitude termed envelope (Figure 2D, middle and bottom). Such movement envelopes typically contain power at very low (<1 Hz) frequencies (; Metzen and Chacron, 2014). Movement envelopes will elicit behavioral responses in which the animal’s EOD frequency “tracks” the envelope stimulus (Metzen and Chacron, 2014).

Most studies investigating neuronal coding in weakly electric fish were performed in immobilized animals (note that the EOD persists after immobilization in species such as A. leptorhynchus). However, electrosensory behaviors consist of changes in the animal’s electric field and, as such, can be also be elicited in immobilized animals (). Investigators have taken advantage of this fact to gain better understanding as to the nature of the electrosensory neural circuits that give rise to behavior (). Recent studies have shown that the responses of electrosensory neurons to stimuli associated with different contexts strongly determine behavioral output (; , ; Metzen et al., 2016a, 2018; Metzen and Chacron, 2017). Some of these results pertaining to population coding are described below.

The Anatomy of the Electrosensory Pathway

At the peripheral level, electrosensory stimuli (AM of the EOD) are detected by about 16,000 tuberous electroreceptors or electrosensory afferents (EA) that are distributed across the animals’ body and embedded in its skin (). At baseline (i.e., in the absence of stimulation), each EA fires with a specific firing probability but phase locked to the EOD carrier wave and stimulation will cause changes in the firing probability (Scheich et al., 1973; ).

EAs project to the hindbrain (Figure 3, bottom left), the first processing station in the brain, where they trifurcate to make synaptic contact with pyramidal cells in the electrosensory lateral line lobe (ELL). The ELL is organized in three parallel somatotopic maps of the body surface: the lateral, the centro-lateral, and the centro-medial segment (LS, CLS, and CMS) (; ; ). All three segments are composed of columns as a repetitive motif (Figure 3, top left), with each column consisting of six different pyramidal neurons () (Figure 3, right). Three of these neurons (“on”) receive direct excitatory input from EAs and respond with increases in their firing rate to increases in the AM. The three other neurons (“off”) receive the EA input via an inhibitory granular interneuron (“gr”) and thus instead respond to increases in the AM with decreases in firing rate (; Saunders and Bastian, 1984; ). Pyramidal neurons are the sole output neurons of the hindbrain and project to the midbrain torus semicircularis (Figure 3, bottom left) where sensory information is further processed, and forwarded to various stages in the forebrain ultimately giving rise to behavior. There are three classes of pyramidal neurons: superficial, intermediate, and deep neurons, named after where their cell body is located within the pyramidal cell layer of the ELL (). The three classes differ in terms of their cell morphology, physiology, and connectivity. Superficial cells have low baseline firing rates, large apical dendritic arborization that extend widely through the molecular layer, and receive huge amounts of descending inputs (orange), deep cells have high baseline firing rates and receive little to no descending inputs (; ; ; ; ). The properties of intermediate cells are in between the others.

FIGURE 3

There are two major types of descending pathways to ELL that originate from higher brain areas (Figure 3, orange). Inputs from the nucleus praeminentialis (nP, commonly termed direct feedback pathway) as well as from the eminentia granularis posterior (EGP, commonly termed indirect feedback pathway) form parallel fibers in ELL that make contact with the apical dendrites of pyramidal neurons. Direct contacts of parallel fibers with pyramidal neurons are excitatory, indirect contacts through interneurons in the molecular layer are inhibitory. Both types of descending inputs can strongly affect the responses of single ELL pyramidal cells to stimulation (; ; ; ; ; Metzen et al., 2018). The detailed anatomy and function of these pathways have been reviewed elsewhere (; ).

Population Coding by Electrosensory Afferents

Correlation-Based Coding of Chirps by Electrosensory Afferents

As described above, chirp waveforms of a given chirp are very heterogeneous in nature depending on the phase of the AM at which they are emitted (“chirp phase”). Nonetheless, chirps occur with equal probability at any phase of the beat during electrocommunication (). The chirp waveforms can, depending on the chirp phase, consist of sharp increases, sharp decreases (Figure 4A, top traces, gray shading) or biphasic high-frequency transients in the AM. All these different chirp waveforms, however, were shown to elicit similar behavioral responses (Metzen et al., 2016a) suggesting that the social content of a given chirp is independent of the beat phase. At the level of EAs, the heterogeneous chirp waveforms will elicit very heterogeneous responses in the firing activities of single EAs (Figure 4A, middle curves) that have been well-characterized (, ; Walz et al., 2014). As a result, the responses of EAs, in terms of their firing rate, are highly different depending on the chirp phase (Figure 4B, dotted line) while the behavioral responses are ultimately not (Figure 4B, orange line). Chirp stimulation increased the similarity of the firing patterns of pairs of EAs, thereby causing an increase in their correlations (Figure 4A, bottom) (Metzen et al., 2016a). These increases in correlations were very similar across the different possible waveforms of a given chirp (i.e., correlation-based responses were invariant) (Figure 4B, solid black line). This correlation response closely resembled the behavioral invariance to different chirp waveforms (Figure 4B, compare black and orange solid lines).

FIGURE 4

Interestingly, these results were similar for both simultaneous and non-simultaneous recordings. As explained above, non-simultaneous recordings cannot be used to infer noise correlations but only signal correlations. This suggests that changes in EA’s correlations are primarily driven by changes in signal correlations and that noise correlations between EAs are either negligible or do not affect such coding. This is supported by other studies showing that correlations between EA baseline activities are negligible except at the EOD frequency and higher harmonics (). Based on recordings from several other processing stations (i.e., ELL and Torus semicircularis), Metzen et al. (2016a) show how EA correlations can be decoded in a physiologically plausible manner by combining and integrating parallel inputs along the ascending electrosensory pathway. Interestingly the correlation-based detectability of chirps depends on stimulus background: increasing beat frequency impairs detectability. Importantly the behavioral detection performance declines in parallel with correlation based coding performance (Metzen and Chacron, 2017). These results strongly suggest that correlations in EA activity are decoded by downstream brain areas to give rise to behavior.

Correlation-Based Coding of Envelopes by Electrosensory Afferents

Recent studies have focused on understanding how EA’s respond to envelopes through changes in firing rate (Metzen and Chacron, 2015). While envelopes do not elicit changes of the EA’s average firing rate (as compared to baseline), it was found that they caused changes in the similarity of firing patterns (Figure 4C, gray boxes) and thus in the correlation between EA’s. Therefore, they could be detected when analyzing the firing patterns of the EA population. In fact the correlation magnitude in the EA population nicely tracks the envelope waveform (Figure 4C, solid lines) (Metzen et al., 2015b), and there was a strong relationship between the envelope and correlation (Figure 4D). This indicates that also second-order stimulus features are encoded through neuronal correlations, which might give rise to previously observed behavioral responses by which the animal’s EOD frequency tracks the detailed time course of the envelope in an almost one-to-one fashion (Metzen and Chacron, 2014).

As found during EA correlation encoding of chirps, results were almost independent of whether simultaneous or non-simultaneous recordings were used for the analysis (Figure 4D, inset), suggesting that changes in correlation are primarily, if not exclusively, driven by changes in signal correlations between EA’s. Interestingly, theory predicts that such envelope coding by correlated activity was optimal for a given level of variability (Metzen et al., 2015a; ), a prediction that was verified experimentally (Metzen et al., 2015b). Moreover, such coding appears to be a general feature of sensory processing, with similar results found in the coding vestibular afferents (Metzen et al., 2015b) and acoustic processing () in non-human primates, as well as in LGN of the cats visual system ().

ELL Pyramidal Neurons Exhibit Correlated Variability: Mechanisms and Implications for Coding

ELL pyramidal neurons receive convergent input from EA projections: Anatomical studies have shown that up to about 65% of inputs are shared between pyramidal neurons in neighboring columns (). As such, it is expected that, unlike EA’s, pyramidal cells in the hindbrain will display both signal and noise correlations, i.e., exhibit correlated variability. In this section, we will start by reviewing how correlations between the baseline activities of ELL pyramidal cells arise and demonstrate that these serve as a good predictor of noise correlations under stimulation. We then move on and review the state of knowledge regarding the presence and plasticity of noise correlations in ELL.

Mechanisms Mediating Baseline Correlations

Under baseline conditions, ELL pyramidal neurons fire action potentials in an irregular pattern switching between bursts and single spikes (Figure 5A) (; Metzen et al., 2016b) with discharge rates in the range of a few to about 40 Hz (). The spike trains of simultaneously recorded neighboring pyramidal neurons are typically correlated in the absence of stimulation (; ). The magnitude of these baseline correlations (i.e., raw-correlations recorded in absence of stimulation; see section “Decomposing Correlations” ), will typically increase from small (<10 ms) to large (>1 s) time windows (Figure 5B). Similarly, when calculating a cross-correlogram (CCG) a prominent peak near lag zero is visible, but coincident events above chance level are also found at higher lags, e.g., at 50 ms or higher (Figure 5C). The magnitude of baseline correlations is independent of the difference in firing rate between the neurons in a pair and is stationary over time (). Baseline correlations are on average positive between pairs of the same type (i.e., on–on and off–off) and negative between opposite type pairs (on–off) (Figure 5D). By mapping the receptive fields (RF, the area on the skin or in the environment within which a stimulus causes a response in the neuron) of pyramidal neurons, were able to show a positive correlation between the amount of RF overlap of neuron pairs and their baseline correlation magnitude (Figure 5E). Such RF overlap is likely to be caused by shared EA input between pyramidal neurons ().

FIGURE 5

), data in (E) from ().

The average RF size and overlap between neighboring neurons was estimated to differ between ELL segments and decrease from the lateral segment (LS) over the centro lateral segment (CLS) to the centro medial segment (CMS) (). Based on the relationship between baseline correlation magnitude and RF overlap (Figure 5E) one would consequently predict baseline correlations to decrease from LS to CMS. This is not true, however, as the correlation magnitudes were found to be similar on average, in all three segments (Figure 6A) (). These similar magnitudes of baseline correlations most likely originate from different RF properties as described below.

FIGURE 6

). Statistical significance; n.s., not significant.

ELL RFs are organized in an antagonistic center-surround organization (Shumway, 1989; ). This means, for an on-type cell, stimulation within the RF center will cause an increase in firing rate, while stimulation within the RF surround will instead cause a decrease in firing rate. Thus, when considering the RFs of two neighboring ELL pyramidal neurons, one must consider up to eight different areas of RF overlap, depending on their spacing (Figure 6B). These areas will give either excitatory (“+”) or inhibitory (“-”) input to each of the two neurons, or will not project to a given neuron (“0”). Comparing the inputs, each neuron receives from each of the areas of overlap, one can expect these signals to be positively correlated (“+corr”), negatively correlated (“-corr”), or not correlated (“no corr”) (Figure 6B, labels). The balance between these inputs will determine the amount of input correlations. The correlation magnitude measured in a pair of pyramidal neurons will depend on their input correlations and the amount of correlation transfer in these neurons (Shea-Brown et al., 2008; ).

Using mathematical models and numerical simulations, it was found that changing the relative strength and size of the RF surround relative to that of the RF center impacts the magnitude of baseline correlations (Figure 6C). Interestingly, many and very different RF topographies and balances led to correlation magnitudes within the physiological range () (Figure 6C, red). As such, it is not only interactions between the RF centers but also interactions between the centers and the surrounds, as well as interactions between the RF surrounds themselves that contributes toward determining correlation magnitude. Therefore, the similar correlation magnitudes seen across the three ELL segments occur because decreases in RF center overlap when going from LS to CMS are effectively “compensated for” by the concomitant impact of signals from the RF surrounds (Figure 6D) (Shumway, 1989; ). As mentioned above, the presence of baseline correlations in the absence of stimulation strongly suggests that noise correlations will be present during stimulation.

Baseline Correlations Predict Noise Correlations Under Stimulation in ELL Pyramidal Cells

Under stimulation, pyramidal neurons will typically encode the stimulus waveform through changes in firing pattern while the overall firing rate changes only minimally on average (Figure 7A). This is generally attributed to gain control and the cancelation of redundant signals via descending pathways (i.e., indirect feedback) (, ; ; ; ) causing pyramidal neurons to adapt to both first- and second-order stimuli (; ; Zhang and Chacron, 2016). Simultaneously recorded spike trains of pyramidal neurons will display both signal and noise correlations (Figures 7B,C). Using previously published data (; ), we investigated the relationship between baseline correlations and noise correlations and found a strong relationship between them (Figure 7D). This confirms our earlier prediction that baseline correlations are a precursor for the presence of noise correlations under the assumption of weak stimulus amplitudes (). We note that the magnitude of noise correlations was systematically lower compared to the magnitude of baseline correlations (compare slope of fits to identity line) while the sign was preserved between same and opposite type pairs (Figure 7D, dots vs. triangles).

FIGURE 7

). (E) The magnitudes of signal and noise correlations were not systematically dependent on each other (red lines, fit to individual datasets, slopes were not significant). However, their sign seems to be preserved in general, i.e. correlation structure in ELL is on average parallel. (F,G) Raw correlations as a function of signal (F) and noise correlations (G). In both cases, strong and significant relationships were found indicating that both components contribute to the overall correlation coefficient. However, noise correlations vary over a larger range, and the relationship was stronger suggesting that the impact of noise correlations slightly outweighs that of signal correlations. Data in (A–C) re-analyzed from . Data in (D–G) reanalyzed from .

We further compared noise and signal correlations and found that the correlation structure in ELL is mostly parallel (i.e., signal and noise correlations typically have the same sign) (Figure 7E, compare also to Figures 1A,B). However, when considering either same or opposite-type pairs separately, there was no significant relationship between signal and noise correlation magnitudes (Figure 7E). As expected, both signal and noise correlations seem to contribute to the overall (raw-) correlation in ELL as, for both, significant relations were found (Figures 7F,G). The relationship between raw and noise correlations extends over a larger range (range of rnoise: -0.35 to 0.39) as compared to the relation between raw and signal correlations (range of rsignal: -0.28 to 0.37). Furthermore, the spread of the data is less (error areas of the fits are smaller). This could be seen as indication that noise correlations more strongly determine raw correlations between ELL spike trains than signal correlations are.

Based on the observed correlation structure (Figure 7E), one would predict that information transmission is compromised when assuming a decoder that relies on the linear sum of responses. It is, however, important to note that experimental data has shown that such a decoding scheme is most likely not completely accurate (Vonderschen and Chacron, 2011; ). Nevertheless, these results highlight the important fact that noise correlations between ELL pyramidal cell activities should not be assumed to be negligible and cannot be ignored when investigating population coding (; ; ; ).

Plasticity of Correlated Variability in ELL

Similar to what is reported in other brain areas, noise correlations in ELL are highly plastic (see also “Effects of Correlated Variability on Stimulus Encoding” section). Indeed, it was shown that their magnitude strongly depend on the stimulus’ spatial extent (; ; Simmonds and Chacron, 2015). Specifically, correlations in a given pair are low when using stimuli whose spatial extent mimics those caused by conspecifics (global stimulation; Figure 8A) compared to when using stimuli with the same temporal profile but whose spatial extent mimics those caused by prey (local stimulation; Figure 8B).

FIGURE 8

This effect is timescale specific: on short timescales (<10 ms) correlations slightly increased under global stimulation while on longer timescales (>10 ms), correlations strongly decreased () (Figures 8C,D, arrows). showed that this correlation shaping was in part due to changes in signal correlation, which reflects previously described changes in the response properties of single ELL pyramidal neurons (; ). Moreover, noise correlations were in general weaker under global stimulation than under local stimulation (Figures 8E,F, arrows).

It was shown, both by mathematical modeling and experimental manipulation, that the reduction of noise correlations under global stimulation is due to activation of the indirect feedback pathway (Simmonds and Chacron, 2015). This descending input pathway is diffuse and activated only under global but not local stimulation (). As explained above (see section “The Mechanistic Origins of Correlations in vivo Remain Poorly Understood”), noise correlations likely arise due to the shared noise in common inputs. For ELL, these are the afferent inputs from EAs (Figure 9A, bottom). As for the descending inputs, it is assumed that the granule cells within the EGP do not fire in the absence of stimulation but are active during stimulation. Further, it is assumed that the trial-to-trial variability in the granule cell firing activities contribute independent noise to the pyramidal cells due to the diffuse nature of the descending inputs (Figure 9A, left). Activation of the indirect feedback during global stimulation will therefore “dilute” noise correlations between pyramidal cells. Indeed, experimentally blocking this pathway during global stimulation led to an increase in noise correlations supporting this hypothesis (Figures 9B,C, arrows). As such, the descending inputs during global stimulation can be implicated as one functional component with which correlation plasticity is achieved in ELL. If and how this affects the processing of behavioral relevant signals is discussed below.

FIGURE 9

; ; Simmonds and Chacron, 2015). Statistical significance.

Future Directions

In the following, we highlight interesting future avenues of research on population coding in the electrosensory system.

Noise correlations have been shown within the CLS segment (; Simmonds and Chacron, 2015). Furthermore, the presence of baseline correlations was shown for all segments () which, with the re-analyzed data presented here (Figure 7D), strongly implies that noise correlations will be present during stimulation in all ELL segments. Systematic assessments of differences in noise correlations between pairs of different pyramidal neuron types (i.e., superficial, intermediate, and deep) have not been done so far. Baseline correlation magnitude was reported to weakly correlate with the average baseline firing rate in CLS pairs (), suggesting that pairs of deep cells tend to display slightly higher correlation magnitude. Whether this holds true for noise correlations under stimulation has not been investigated to date. Further studies are also needed to assess how pyramidal cell heterogeneities affect correlation plasticity. Specifically, as deep cells receive less descending inputs in comparison to intermediate and superficial neurons, we predict that the reduction in noise correlation during global stimulation will be less pronounced in these cell pairs.

The initial estimation of correlation structure (Figure 7E) suggests that noise correlations might influence signal encoding in a detrimental fashion with respect to many of the analytical tools that were used in the past. However, it is important to note that it remains to be shown directly if and how they influence signal encoding, which will require evaluation from the decoding perspective also. For this, recordings from areas downstream of ELL such as the midbrain torus semicircularis will be required. There, the diversity of cell classes and response properties (Vonderschen and Chacron, 2011; ; ; Sproule et al., 2015) could imply that different decoding strategies are used for different stimuli.

One important area of research concerns how natural electrocommunication stimuli (i.e., chirps) are encoded within the electrosensory pathway. As mentioned above, most studies have focused on how single EA’s or ELL pyramidal cells encode such stimuli. While extrapolations to the population level were attempted, the potential effects of noise correlations were generally neglected (; ). However, the fact that the baseline activities of LS pyramidal cells are correlated at timescales commensurate with those of chirps suggests that noise correlations need to be taken into account in future studies. Simultaneous recordings of ELL neurons during stimulation with chirps will be required to verify the above prediction and ascertain their effects on population coding. Further, the potential impact of noise correlations on population coding of chirps will require investigation of the decoding in downstream midbrain neurons. These integrate converging inputs from ELL pyramidal neurons and experimental studies have shown that some midbrain neurons, due to their non-linear integration of inputs from on- and off-type ELL pyramidal neurons, responded to chirps in a selective manner (Vonderschen and Chacron, 2011; ; Metzen et al., 2016a). Future studies should consider these more physiologically realistic decoding schemes, in particular we predict that pooling the activities of on- and off-type cells will reduce their overall responses to the beat, thereby making the response to the chirp more detectable.

With regard to electrosensory envelopes, previous studies have largely focused on understanding the encoding by single ELL pyramidal cells (; ; Zhang and Chacron, 2016; Metzen et al., 2018). Further studies are needed to understand the role and impact of noise correlations on population coding of envelopes. Here, responses of on- and off-type ELL pyramidal cells, while responding in andout of phase to first-order stimuli, respectively, actually respond largely in phase to second order stimuli such as envelopes (). From this, one would expect to find ELL signal correlations to be positive with regard to the envelope on average. However, the experimentally observed negative baseline correlations between opposite type pairs, would predict that these display negative noise correlations in response to envelopes. As such, noise correlations are predicted to be beneficial for envelope coding in such pairs. Further studies are needed to verify this prediction.

Finally, it should be noted that theoretical studies have suggested that noise correlations themselves could directly encode stimuli, therefore forming an independent channel of information transmission in the brain (). The documented correlation plasticity in ELL, together with the nature of electrosensory stimuli being intertwined with an active sensing carrier signal, could be an example in which such a correlation code is realized in the brain. If plasticity of noise correlations can be found due to stimulus attributes other than the stimulus spatial extend, is unclear, however, and remains to be investigated.

While many of the above discussed assessments of neuronal encoding are based on recordings from immobilized animals, it is important to note that weakly electric fish, on top of being able to display electrical behaviors when immobilized, show elaborate behaviors and astonishing cognitive abilities and are getting more and more attention for the study of various aspects of active sensing behaviors (Nelson and MacIver, 2006; ; von der Emde et al., 2010; , ; Pedraja et al., 2018). Recent technological advances such as electrophysiological recordings from freely moving aquatic animals are rapidly evolving (; Vinepinsky et al., 2017). Being able to perform such recordings in freely behaving electric fish will allow to combine the investigations of population coding aspects in active sensing contexts – two of the most prominent research streams in neuroscience. Specific questions will likely be: How are active sensing movements generated and controlled through neuronal populations, how do population codes contribute to decision making during active sensing movements, and how do neuron populations encode sensory signals discriminating re- and ex-afferent signals at the population level. Based on the vast body of knowledge regarding behavior, anatomy and physiology weakly electric fish promise to evolve into an exciting model system to study the neuronal control of active sensing behaviors.

Statements

Author contributions

VH and MC conceived the study, reviewed and edited the text and figures. VH curated and re-analyzed the data, prepared the figures, and wrote the initial draft.

Funding

This research was supported by the Deutsche Forschungsgemeins chaft DFG: HO 5912/1-1 (VH); and the Canadian Institutes of Health Research (MC).

Acknowledgments

We would like to thank Dr. Michael G. Metzen for providing digital artwork for Figures 4C,D.

Conflict 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.

References

Summary

Keywords

population coding, correlated variability, noise correlations, stimulus encoding, feedback, electric fish, correlation shaping, electrosensory lateral line lobe

Citation

Hofmann V and Chacron MJ (2018) Population Coding and Correlated Variability in Electrosensory Pathways. Front. Integr. Neurosci. 12:56. doi: 10.3389/fnint.2018.00056

Received

13 September 2018

Accepted

30 October 2018

Published

27 November 2018

Volume

12 - 2018

Edited by

Pedro E. Maldonado, Universidad de Chile, Chile

Reviewed by

Jose Luis Pena, Albert Einstein College of Medicine, United States; Luis Lemus, Universidad Nacional Autónoma de México, Mexico

Updates

Copyright

*Correspondence: Volker Hofmann,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics