Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro

eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleResearch Article: Methods/New Tools, Novel Tools and Methods

Information Theoretic Approaches to Deciphering the Neural Code with Functional Fluorescence Imaging

Jason R. Climer and Daniel A. Dombeck
eNeuro 25 August 2021, 8 (5) ENEURO.0266-21.2021; DOI: https://doi.org/10.1523/ENEURO.0266-21.2021
Jason R. Climer
Department of Neurobiology, Northwestern University, Evanston, 60208 IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel A. Dombeck
Department of Neurobiology, Northwestern University, Evanston, 60208 IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daniel A. Dombeck
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Information theoretic metrics have proven useful in quantifying the relationship between behaviorally relevant parameters and neuronal activity with relatively few assumptions. However, these metrics are typically applied to action potential (AP) recordings and were not designed for the slow timescales and variable amplitudes typical of functional fluorescence recordings (e.g., calcium imaging). The lack of research guidelines on how to apply and interpret these metrics with fluorescence traces means the neuroscience community has yet to realize the power of information theoretic metrics. Here, we used computational methods to create mock AP traces with known amounts of information. From these, we generated fluorescence traces and examined the ability of different information metrics to recover the known information values. We provide guidelines for how to use information metrics when applying them to functional fluorescence and demonstrate their appropriate application to GCaMP6f population recordings from mouse hippocampal neurons imaged during virtual navigation.

  • behavior
  • calcium imaging
  • in vivo
  • information theory
  • place cells
  • two-photon microscopy

Significance Statement

Functional fluorescence imaging and information theoretic quantification could provide a powerful new combination of tools to study neural correlates of behavior, but functional fluorescence signals represent altered versions of the underlying physiological events. Therefore, it is unclear whether or how information metrics can be applied to functional fluorescence imaging data. Here, we performed an in-depth simulation study to examine the application of the widely used bits per second and bits per action potential (AP) metrics of mutual information (MI) to functional florescence recordings. We provide guidelines for how to use information metrics when applying them to functional fluorescence and demonstrate their appropriate application to GCaMP6f population recordings from mouse hippocampal neurons imaged during virtual navigation.

Introduction

Neurons encode parameters important for animal behavior, at least in part, through the rate of production of action potentials (APs). Evidence for this can be found from electrophysiological AP recordings of orientation tuning in the visual system (Hubel and Wiesel, 2009), chemical sensing in the olfactory system (Leveteau and MacLeod, 1966; Wachowiak and Shipley, 2006), and spatial encoding in the hippocampus (O’Keefe, 1976). Key to deciphering the neural code, therefore, is defining metrics to quantify the relationship between behavioral parameter spaces and a neuron’s spiking rate. There are many metrics used for quantification, and are often used to compare neural responses across conditions or in neurons with complex responses. The underlying assumptions of the different metrics then become important factors to consider when determining which one to use.

Information theory is growing in popularity in the neuroscience community, largely because it provides a means to quantify rate coding with relatively few assumptions. One useful information theoretic measure is mutual information (MI), which is typically measured in bits per unit time, and describes the increase in predictability of the neural response when behavioral parameters are known. Formally, MI is the information about one variable that can be extracted from another, such as the information about behavior that can be derived from observing neural activity. MI can be applied to neurons with widely varying response properties because it (1) is a nonlinear metric, not requiring the linearity assumptions of correlation metrics (Grubb and Thompson, 2003; Kropff et al., 2015; Hinman et al., 2016); (2) does not assume a response shape, as is typical with Gaussian field mapping metrics (Soo et al., 2011; Kraus et al., 2015; Tang, 2016) or metrics using exponential or polynomial curve fitting (Hinman et al., 2016); and (3) uses the full time trace or shape of the mean response profile, rather than defining receptive fields with thresholding (Niell and Stryker, 2008; Pastalkova et al., 2008; Harvey et al., 2009).

However, MI can be nontrivial to estimate from neural and behavioral recordings and its estimation is an ongoing area of research (Kraskov et al., 2004; Gao et al., 2017; Belghazi et al., 2018; Timme and Lapish, 2018).

Here, we focus on the most widely used estimator of MI in neuroscience, the SMGM estimator developed by Skaggs, McNaughton, Gothard, and Markus (Skaggs et al., 1993), although as a point of comparison, we also consider the binned estimator (Timme and Lapish, 2018) and a separate technique developed by Kraskov, Stogbauer, and Grassberger (KSG; Kraskov et al., 2004). The binned estimator estimates the joint probability distribution using a 2D histogram of neural response versus behavioral variable; this transforms continuous variables into discrete values (Timme and Lapish, 2018). KSG estimates MI by examining the distance between data-points in the neural activity-behavioral parameter space. The SMGM estimator, on the other hand, relies on the assumption that AP firing follows an inhomogeneous Poisson process. The SMGM estimator therefore requires binning of only the behavioral variable(s), in contrast to the binned estimator. The profile of firing rates versus behavioral variable is then used to estimate the MI.

The relative simplicity of the SMGM estimator has added to its popularity and widespread use in neuroscience applications for estimating behavioral information contained in single unit AP recordings. This metric has proven useful in quantifying rate coding in place cells (Knierim et al., 1995; Markus et al., 1995; Lee et al., 2006; Poucet and Sargolini, 2013), complex spatial responses of hippocampal interneurons (Frank et al., 2001; Wilent and Nitz, 2007), odor sequence cells (Allen et al., 2016), time cells (MacDonald et al., 2013), head direction cells (Stackman and Taube, 1998), speed cells (Fyhn et al., 2002), and face differential neurons (Nguyen et al., 2013, 2014), and has been used across multiple different species (Yartsev and Ulanovsky, 2013; Hazama and Tamura, 2019; Mankin et al., 2019). Furthermore, as a single neuron metric, it provides statistical power for comparisons. Thus, it has been used to quantify differences in rate coding across different brain regions (Simonnet and Brecht, 2019) and across experimental interventions such as lesions (Calton et al., 2003; Liu et al., 2004), inactivations (Huang et al., 2009; Koenig et al., 2011; Brandon et al., 2011; Hok et al., 2013), and applications of drugs (Robbe and Buzsáki, 2009; Newman et al., 2014). Further, it has been used to examine differences in encoding across different behaviors (Zinyuk et al., 2000; Park et al., 2011; Aronov and Tank, 2014) and disease states (Zhou et al., 2007; Gerrard et al., 2008; Fu et al., 2017). SMGM information is often normalized from measuring bits per unit time to instead measure bits per AP. This creates a measure sensitive only to the selectivity of a neuron, and not its average firing rate. Thus, SMGM is a powerful tool for measuring the neural code in electrophysiological recordings of APs.

The power of MI estimators has yet to be fully exploited by the neuroscience community. For example, the estimators have not yet been widely used to compare encoding properties of large numbers of genetically identified neurons, or to quantify information content of other discrete signaling events such as synaptic inputs; both of which are difficult to study using electrophysiological methods. In vivo imaging of functional indicators has emerged as an important tool, largely because it possesses these capabilities. For example, using fluorescent calcium indicators, the functional properties of large populations of neurons can be simultaneously recorded in rodents (Dombeck et al., 2007; Ziv et al., 2013; Stirman et al., 2016; Sheffield et al., 2017; Radvansky and Dombeck, 2018; Stringer et al., 2019), zebrafish (Ahrens et al., 2013), or invertebrates such as Caenorhabditis elegans (Nguyen et al., 2016) and Drosophila (Keller and Ahrens, 2015; Mann et al., 2017). Furthermore, in vivo imaging can assure the genetic identity of the recorded neurons (Khoshkhoo et al., 2017; Sheffield et al., 2017; Jing et al., 2018a,b,c) and can access subcellular structures, allowing for functional recordings from synapses and dendrites using different functional fluorescent indicators (Sheffield and Dombeck, 2015; Scholl et al., 2017; Sheffield et al., 2017; Jing et al., 2018d; Marvin et al., 2018, 2019; Adoff et al., 2021).

However, these indicators generate signals that are different from the underlying quantal events. For example, somatic calcium indicators reveal intensity variations that are correlated with somatic AP firing rates but are a smoothed and varying amplitude version of the AP train. This transformation from AP train to fluorescence trace is an active area of research (Dana et al., 2018; Greenberg et al., 2018; Éltes et al., 2019), but it is often approximated by convolving the AP train with a kernel, which defines the indicator’s response to a single AP. The shape of the kernel is a function of the indicator expression level, intracellular calcium buffering, amount of calcium influx, efflux rates, background fluorescence, resting calcium concentration, and other factors. When measured in pyramidal neurons, average kernels typically take the shape of a sharp increase in fluorescence followed by an exponential decay to baseline (Yaksi and Friedrich, 2006; Chen et al., 2013; Park et al., 2013; Dana et al., 2018; Pachitariu et al., 2018). Therefore, while functional fluorescence imaging and information theoretic quantification may prove to be a powerful new combination of tools to study neural correlates of behavior, it is critical to remember that functional fluorescence signals represent altered versions of the underlying physiological events.

Caution is then needed when applying information metrics to continuous functional fluorescence traces, yet the imaging community is already beginning to use information metrics, particularly SMGM. This metric has been applied to somatic calcium responses to compare the information content of the same neurons across different behavioral epochs (Heys and Dombeck, 2018), across different populations of neurons in different brain regions (Hainmueller and Bartos, 2018), across different genetically identified neural populations (Khoshkhoo et al., 2017), or to examine encoding by subcellular structures (Rashid et al., 2020), or to classify the significance of encoding particular parameters by individual neurons (Kinsky et al., 2018; Mau et al., 2018; Rashid et al., 2020).

However, it is essential to recognize some of the assumptions underlying these information metrics are violated by functional florescence recordings. All three metrics (SMGM, KSG, and binned estimation) assume stationarity in the neural response, which is violated by the elongated time responses and relatively slow fluctuations of the fluorescence intensity of the reporters. When applied to spiking data, there is also a change in units: rather than AP counts, functional fluorescence traces are typically plotted in units of florescence change with respect to baseline (ΔF/F). One possible solution to these issues would be to deconvolve calcium traces to recover APs; however, deconvolution is an active area of research, and the accuracy of these methods has recently been questioned (Evans et al., 2019). Ideally, the calcium traces could be used directly to measure spiking information, without the need for such an in between, potentially error inducing, step.

Quantifying the effects of the above violations on measurements of information using functional fluorescence recordings with an analytical solution is particularly challenging with behaviorally modulated neural recording data. However, a more tractable means of quantifying the effects would be to use a simulation study to measure the induced biases and changes in measurement quality (Morris et al., 2019). This strategy makes use of pseudo-randomly generated AP traces and has the advantage that the ground truth parameters of the simulations are known, while variability because of behavior and other features can be incorporated (Cohen and Kohn, 2011; Climer et al., 2013, 2015; Østergaard et al., 2018).

To provide the field with guidelines for the use of information metrics applied to functional fluorescence recording data, we used computational simulation methods to create a library of ten thousand mock neurons whose spiking output carry an exact, known (ground-truth) amount of information about the animal’s spatial location in its environment. We used real behavioral data (available at https://doi.org/10.7910/DVN/SCQYKR) of spatial position over time from mice navigating in virtual linear tracks and then simulated the spatial firing patterns of the mock neurons using an inhomogeneous Poisson process framework (Brown et al., 2003; Paninski, 2004; Climer et al., 2013). We then simulated fluorescent calcium responses for each neuron in each session by convolving the AP trains with calcium kernels for different indicators, primarily GCamp6f (Chen et al., 2013), and then we added noise. MI metrics (between spatial location and the neural signals) were then applied to the spiking or fluorescence traces to quantify the performance of the metrics for estimating information. We provide a user toolbox (found at https://github.com/DombeckLab/infoTheory), which consists of MATLAB functions to generate libraries of model neurons with known amounts of information, to generate spiking or fluorescence time-series from those model neurons, and to estimate neuron information from real or model spiking or fluorescence time-series datasets using the three metrics considered here (SMGM, binned estimator, KSG). We focused on testing the performance of the SMGM method, and then compared its performance to the binned estimation and KSG methods, which do not have the underlying Poisson assumption required for the SMGM approach. We also applied a deconvolution algorithm to test its performance. We then applied this analysis to real datasets of hippocampal neuron populations from mice navigating in virtual linear tracks. We quantified the spatial information content of the populations and then performed Bayesian decoding of mouse position from different information containing subsets of this population. Interestingly, we found that the population quantile with the lowest information values were still able to decode mouse position to the closest quarter of the track. Thus, we provide new findings about the neural code for space that were made possible by the information metrics and guidelines that we introduce here.

The SMGM method applied directly to the mean ΔF/F intensity map appeared to best recover the ground truth information. We provide guidelines for the use of the SMGM metric when applied to functional fluorescence recordings and demonstrate the appropriate application of these guidelines to GCaMP6f population recordings from hippocampal neurons in mice navigating virtual linear tracks.

Materials and Methods

Toolbox and data availability

We provide a user toolbox (freely available at https://github.com/DombeckLab/infoTheory), which consists of MATLAB functions to generate libraries of model neurons with known amounts of information, to generate spiking or fluorescence time-series from those model neurons, and to estimate neuron information from real or model spiking or fluorescence time-series datasets using the three metrics considered here (SMGM, binned estimator, KSG). This toolbox also contains tools to generate mock neurons using a binned distribution, avoiding the Poisson assumption of SMGM. Behavioral data used to generate the random traces is freely available at https://doi.org/10.7910/DVN/SCQYKR.

Construction of AP trains with known ground truth information

To construct mock neurons with ground truth information, we adapted the differential form of the AP information, in bits per AP (Eq. 6). To create a rate map, we first selected an average firing rate and target ground truth information. The mean rate (λ¯ ) was always between 0.1 and 30 Hz, the information in bits per AP (IAPE ) between 0 and 6 bits/AP, and the information in bits per second (IsE ) between 0 and 24. To more evenly sample each of these, we first randomly selected the bits per second (IsE ) or bits per AP (IAPE ) to target. If the information target was in bits per AP, both the information (IAPE ) and mean firing rate (λ¯ ) were chosen uniformly. Because the information in bits per second IsE=λ¯IAPE , the bits per second information was not uniformly sampled in this case. If the target was to be in bits per second, both the bits per AP (IAPE ) and SMGM bits per second (IsE ) measures were first chosen uniformly. Because the rate λ¯=ISE/IAPE , this was not chosen uniformly. This procedure was repeated to maintain the bounds on λ¯ , resulting in a non-uniform sampling of information. The final distribution (Fig. 1C) was spread acceptably for further analysis.

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

Procedures for generating a library of 10,000 neurons with known amounts of information. A, Five splines with a gradient of ground truth information (IAPE ) representing the steps in generating a continuous rate map (λ(x) ) matching the desired target information, in this case, 2 bits/AP. Red Xs indicate control nodes that were moved to change the shape of the spline and minimize the squared error to the target information. B, Cross-section of the error surface around the solution point as a function of the position of node 3, and the trajectory taken by the solver to minimize the error and arrive at the target. C, Histograms of ground truth information resulting from repeating the procedure in A, B 10,000 times to target a range of ground truth information values in bits per second (ISE ). D, Splines representing λ(x) ) and bits per AP (IAPE at the solution point for a low (IAPE=0.04 bits/AP, left) and high (IAPE , 2 bits/AP, right) information neuron. E, Steps to generate mock AP and functional fluorescence data. (1) An example real behavior trace from a mouse running on a linear track that was used to generate the simulated spiking. (2) The behavior in combination with the rate maps generated in A–D were used to generate an instantaneous firing rate trace. (3) The instantaneous rate was used to pseudorandomly generate APs, as shown in this mock raster. (4) The AP raster was convolved with the GCaMP6f kernel (red, inset), and noise was added to generate a mock ΔFF trace. (5) Large numbers of these traces were generated and used to assess the effects of many simulation parameters on the estimators. F–L, Spiking and fluorescence activity patterns generated from the example simulated neurons shown in D and using a mean firing rate of 1 Hz. F, Behavioral trace in blue with AP raster shown in red. G, Lap-by-lap raster of the neurons’ firing versus mouse track position. H, Lap by lap binned, firing rates versus mouse track position for the neurons. I, AP raster (red) and mock calcium traces for the same behavioral period shown in F. J, Lap by lap mean binned fluorescence versus mouse position for the neurons. K, Binned average firing rate (λi , black) and fluorescence intensity (fi , green) maps for the two neurons. These maps were used for information analyses.

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

Quantification of the precision of the SMGM bits per second metric using APs or functional fluorescence recordings. A, Three representative mock neurons spanning the range of ground truth information values in bits per second (ISE ). From top to bottom for each, Mouse track position versus time, AP raster, fluorescence calcium trace (green), and firing rate map (λi , black) and change in fluorescence map (fi , green). B–D, The ground truth bits per second values are well recovered when measured from AP traces. B, Information measured from AP data using the SMGM bits per second metric (IsÊ ) versus ground truth information (IsE ). Each dot is a single mock neuron, the gray dashed line is the unity line (perfect measurement), the pink line is the line of best fit. Red circles show the examples in A. C, Percentage error for the information measurements shown in B. D, Heat map of percentage error measurements shown in C. Black lines are 2 SDs, the white line is the mean. E–G, Effects of applying the SMGM bits per second metric to fluorescence traces. E, Information measured from mock GCaMP6f traces using the SMGM bits per second metric (IsF̂ ) versus ground truth information (IsE ). F, Percentage error for the information measurements shown in E. G, Heat map of percentage error measurements shown in F. H, Representative mock kernels mimicking responses from different indicators. I–K, The effect of kernel height on estimating ground truth information (IsE ) using the SMGM bits per second metric (IsF̂ ). Kernel height for the kernels shown in H are indicated by colored triangles. I, Percentage error as a function of kernel height. J, Heat map of percentage error measurements shown in I with mean (white) and 2 SDs (black). K, The average percentage error as a function of kernel height and ground truth information in SMGM bits per second (IsE ). L–N, The effect of kernel width on estimating ground truth information (IsE ) using the SMGM bits per second metric (IsF̂ ). Kernel widths for the kernels shown in H are indicated by colored triangles. L, Percentage error as a function of kernel width. M, Heat map of percentage error measurements shown in L with mean (white) and 2 SDs (black). N, The average percentage error as a function of kernel width. Recording density affected the metrics (Extended Data Figure 2-1). Changing the kernel to common indicators yielded qualitatively similar, but quantitatively different results (Extended Data Figure 2-2).

Extended Data Figure 2-1

Effect of recording density on information metrics. A, The mean error (top) and absolute error (bottom) between the ground truth information (IsE ) and measured information in SMGM bits per second from mock APs (IsF̂ ) as a function of the mean firing rate (λ¯ ). B, As A, but as a function of number of laps. C, D, As A, B, but for measurements from mock fluorescence traces. E, F, As A, B, but for the SMGM bits per AP metric. G, H, As C, D, but for the SMGM bits per AP metric. Download Figure 2-1, EPS file.

Extended Data Figure 2-2

Effects of applying the SMGM bits per second metric to fluorescence traces from different common functional indicators. Top, Information measured from mock traces using the SMGM bits per second metric (IsF̂ ) versus ground truth information (IsE ). Each dot is a single mock neuron, the gray dashed line is the unity line (perfect measurement), the pink line is the line of best fit. Middle, Percentage error density plots. The white line is the mean, the black are ± 1 SD. Bottom, Summary statistics and estimates for the scaling factor c. Download Figure 2-2, EPS file.

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

Quantification of the precision of the SMGM bits per AP metric using APs or functional fluorescence recordings. A, Three representative mock neurons spanning the range of ground truth information values in bits per AP (IAPE ). From top to bottom for each, Mouse track position versus time, AP raster, fluorescence calcium trace (green), and firing rate map (λi , black) and change in fluorescence map (fi , green). B–D, The ground truth bits per AP values are well recovered when measured from AP traces. B, Information measured from AP data using the SMGM bits per AP metric (IAPÊ ) versus ground truth information (IAPE ). Each dot is a single mock neuron, the gray dashed line is the unity line (perfect measurement). Red circles show the examples in A. C, Percentage error for the information measurements shown in B. D, Heat map of percentage error measurements shown in C. Black lines are 2 SDs, the white line is the mean. E–G, Effects of applying the SMGM bits per AP metric to fluorescence traces. E, Information measured from mock GCaMP6f traces using the SMGM bits per AP metric (IAPF̂ ) versus ground truth information (IAPE ). F, Percentage error for the information measurements shown in E. G, Heat map of percentage error measurements shown in F. H–J, The effect of kernel height on estimating ground truth information (IAPE ) using the SMGM bits per second metric (IAPF̂ ). Kernel height for the kernels shown in Figure 2H are indicated by colored triangles. H, Percentage error as a function of kernel height. I, Heat map of percentage error measurements shown in H with mean (white) and 2 SDs (black). J, The average percentage error as a function of kernel height and ground truth information in bits per AP (IAPE ). K–M, The effect of kernel width on estimating ground truth information (IAPE ) using the SMGM bits per AP metric (IAPF̂ ). Kernel widths for the kernels shown in Figure 2H are indicated by colored triangles. L, Percentage error as a function of kernel width. M, Heat map of percentage error measurements shown in L with mean (white) and 2 SDs (black). N, The average percentage error as a function of kernel width. Changing the kernel to common indicators yielded qualitatively similar, but quantitatively different results (Extended Data Figure 3-1). These errors could not be resolved by changing the bin width (Extended Data Figure 3-2). Addition of a nonlinearity further distorted the measured information (Extended Data Figure 3-3).

Extended Data Figure 3-1

Effects of applying the SMGM bits per AP metric to fluorescence traces from different common functional indicators. Top, Information measured from mock traces using the SMGM bits per AP metric (IAPF̂ ) versus ground truth information (IAPE ). Each dot is a single mock neuron, the gray dashed line is the unity line (perfect measurement), the pink line is the exponential fit. Middle, Percentage error density plots. The white line is the mean, the black are ±1 SD. Bottom, Summary statistics and percent error cutoffs. Download Figure 3-1, EPS file.

Extended Data Figure 3-2

Effect of number of bins on the SMGM metrics. A, The mean percentage error (top) and the SD (bottom) for the bits per second measure (ISF̂ ) applied to 10,000 mock gCamp6f traces, with ground truth information on the x-axis and 2–60 bins (3-m track) on the y-axis. B, As A, for the bits per AP measure (IAPF̂ ). C, The mean percentage error (top) and the SD (bottom) for the bits per second measure (ISF̂ ) applied to 20,000 mock fluorescence traces with differing kernel width. Kernel width is on the x-axis, number of bins is on the y-axis. D, As C, but for the bits per AP measure (IAPF̂ ). Download Figure 3-2, EPS file.

Extended Data Figure 3-3

Effects of a sigmoid nonlinearity between ΔF/F and firing rate. Here, we applied a log-sigmoid nonlinearity to the 10,000 mock GCaMP6f time-series traces and then measured information using the fluorescence SMGM metrics. A, Nonlinearity applied to AP-to-florescence trace transformation. B, Information measured from AP data using the SMGM bits per second (ISF̂ ) versus ground truth information (ISE ). Each dot is a single mock neuron, the gray dashed line is the unity line (perfect measurement). C, Percentage error for the information measurements shown in B. D, Heat map of percentage error measurements shown in C. Black lines are 2 SDs, the white line is the mean. E–G, As B–D, but for the bits per AP measured (IAPF̂ ) versus the ground truth information in bits per AP (IAPE ). Download Figure 3-3, EPS file.

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

Alternative techniques for measuring MI from functional fluorescence traces. A–E, top, Information measured from mock GCaMP6f traces versus ground truth information. The gray line is the unity line, the pink line is the best fit saturating exponential. Middle, Percentage error for the information measurements shown on top. Bottom, Heat map of percentage error measurements shown in middle. A, FOOPSI deconvolved traces using the SMGM bits per second metric (Isd̂ ). B, FOOPSI deconvolved traces using the SMGM bits per AP metric (IAPd̂ ). The regularization coefficient had little effect on these results (Extended Data Figure 4-1). C, The KSG measure applied to GCaMP6f traces. D, The binned estimator applied to GCaMP6f traces using uniform bins. E, The binned estimator applied to GCaMP6f traces using equal occupancy bins. The binned estimators were less distorted on the raw AP traces (Extended Data Figure 4-2). F, Table of summary statistics for each measure. P exponential is the p value from the χ2 test used to determine whether a saturating exponential fit is better than a linear fit for the measured versus ground truth information plots. An analytic solution yielded qualitatively similar, but quantitatively disparate results (Extended Data Figure 4-3).

Extended Data Figure 4-1

Effect of changing the regularization coefficient in deconvolution (Vogelstein et al., 2010; Friedrich et al., 2017) on the measured information. A, The measured information using the fluorescence bits per second metric applied to 10,000 mock GCaMP6f traces using regularization coefficients between 0 and 3. B, As A, but for the fluorescence bits per AP metric. Download Figure 4-1, EPS file.

Extended Data Figure 4-2

The binned estimator applied to AP traces and then compared to ground truth information. (Top) Information measured from mock AP traces vs ground truth information. The gray line is the unity line, the pink line is the best fit saturating exponential. (Middle) Percentage error for the information measurements shown on top (same scale as shown in Figure 4). (Bottom) Heat map of percentage error measurements shown in middle. (A) The Binned Estimator applied to AP traces using uniform bins. (B) The Binned Estimator applied to AP traces using equal occupancy bins. Download Figure 4-2, EPS file.

Extended Data Figure 4-3

Information for neurons with Gaussian rate maps. A, blue, Information measured from mock GCaMP6f traces using the SMGM bits per second metric (IsF̂ ) versus ground truth information (IsE ). Red, Information approximated using the analytic approximation in Equation 9. Each dot is a single mock neuron with a Gaussian rate map. The gray dashed line is the unity line (perfect measurement). B, Percentage error for the information measurements shown in A. C, blue, Information measured from mock GCaMP6f traces using the SMGM bits per AP metric (IAPF̂ ) versus ground truth information (IAPE ). Red, Information approximated using the analytic approximation in Equation 10. Each dot is a single mock neuron with a Gaussian rate map. The gray dashed line is the unity line (perfect measurement). D, Percentage error for the information measurements shown in C. Download Figure 4-3, EPS file.

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

Application of SMGM information metrics to functional fluorescence imaging data from hippocampus during spatial behavior. A, Example field of hippocampal pyramidal neurons expressing GCaMP6f and imaged during linear track navigation. Active cell ROIs shown in yellow; traces for green cells shown in B. B, Fluorescence DF/F traces (green) from two neurons in the field shown in A and the track position during the recording (blue). C, Distribution of information values using the fluorescence SMGM bits per second metric (IsF̂ , top) and the fluorescence SMGM bits per AP metric (IAPF̂ , bottom). The gray line indicates the recommended cutoff for reliability using GCaMP6f. D, Plot of IsF̂ versus IAPF̂ for each neuron. Place cells indicated in red and nonplace cells in blue. E, Example non-place cells spanning the information ranges shown in C. Spatial fluorescence map (fi ) shown on left, and average change in fluorescence per track traversal on right. F, Same as E, but for place cells. G, Bayesian decoding of mouse’s track position using different subpopulations of neurons for one example session. From top to bottom, All active neurons, all nonplace cells, and place cells, the first through third quantiles of the SMGM bits per second formulation (IsF̂ ), and the first through third quantiles of the SMGM bits per AP formation (IAPF̂ ). The white dashed line indicates the ground truth position of the animal, the color map indicates the decoded position probability (peak-normalized posterior distribution). H, Decoding accuracy (absolute position decoding error in units of % of track) pooled over all sessions for each neuron group indicated in G. Black bars indicate significant differences by Holm–Bonferroni corrected rank-sum tests (α=0.05 ). Consistent results were obtained when measuring information from real spiking data and simulated florescence traces (Extended Data Fig. 5-1).

Extended Data Figure 5-1

The SMGM estimators as applied to real AP data from a real spiking dataset from hippocampal neurons in mice running on a behavioral track (Chen et al., 2016; Grosmark and Buzsaki, 2016; Grossmark et al., 2016). A, Example real place cell. From top to bottom, Rat track position versus time, real AP raster, mock fluorescence calcium trace generated from real AP trace by convolving APs with GCaMP6f kernel and adding noise (green), and firing rate map (λi , black) and change in mock fluorescence map (fi , green). B, Plot of IsF̂ versus IAPF̂ for each neuron. Place cells indicated in red and nonplace cells in blue. C, The SMGM bits per second metric applied to the real AP traces (ISÊ ) versus the mock fluorescence traces (ISF̂ ; generated from the real AP traces). D, The percentage difference between the SMGM bits per second metric applied to the real AP traces (ISÊ ) and the mock fluorescence traces (ISF̂ ). The mean and SD are indicated in black. In magenta, the mean and SD of 5000 of the mock neuron traces seen in Figures 2, 3, sampled to have the same firing rates as the real neurons and shortened to the same mean session duration. E, Density plot for the data shown in D. F–H, As C–E, but for the bits per AP metric (IAPÊ vs IAPF̂ ). Download Figure 5-1, EPS file.

The rate maps were constructed by spline interpolating across five control points with two anchored at each end of the track, and taking the exponential for each point, and then normalizing by the numerically calculated integral (Fig. 1A,D). To create a map matching the target information, we began with a random spline. The “y” (relative rate) initial position of each node was chosen from a standard normal distribution and the initial “x” (track position) of the three center nodes was chosen uniformly. The nodes were then systematically moved using the MATLAB built in optimizer ‘fmincon’ with constraints preventing the crossing of the center nodes and keeping them on the track, and the ‘OptimalityTolerance’ option set to 0 (Fig. 1A). This was accomplished using the ‘genExpSpline’ function in the toolbox.

We then randomly selected behavioral traces (see below, Behavior) and concatenated sessions until a total time randomly chosen between 3 and 60 min was reached (Figs. 1E,F, 2A, 3A). This was accomplished using the ‘loadBehaviorT’ function in the toolbox. The track positions were normalized and used to build a conditional intensity function (CIF) from the rate function above. The CIF was normalized to match an expected mean rate over the entire session, and the MATLAB built-in ‘poissrnd’ function was used to generate AP times, sampled at 1 kHz. The was accomplished using the ‘genSpikeTrain’ function in the toolbox. Finally, the AP times were binned according to the counts within mock imaging frames sampled at 30 Hz.

Simulated ΔFF traces

To construct the ΔFF traces (Figs. 1E,J,K, 2A, 3A), we first created a single AP response kernel from the peak-normalized sum of two exponentials: g(t)=e−at−e−bt(ab)ab−a−(ab)bb−a,

where t is the time since the AP and a and b are chosen to minimize (1−g(τrise))2+(0.5−g(τrise+τfall))2 where τrise is the rise time in seconds and τfall is the half-fall time in seconds. Deviations in τrise and τfall from baseline were also measured. The kernel g(t) was then multiplied by the indicator height. The kernel parameters were generated using the ‘fluorescenceKernel’ function, and evaluated using the ‘doubleExp’ function in the toolbox.

The GCaMP6f, GCaMP6s, and jRGECO1a heights, rise and fall times were measured as responses to single APs in vivo (Kalko et al., 2011; Chen et al., 2013; Dana et al., 2019): other kernels (Fig. 2H; Extended Data Figs. 2-2, 3-1) were approximated from other experiments presented in the references (seen in Table 1).

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

Properties of indicator kernels used

To define the width of the kernel (Figs. 2L–N, 3K–M), we considered the kernel as a low pass filtered version of the APs. If we normalize the filter to mean 1, it has the Fourier transform (1a+2πfĵ−1b+2πfĵ)(abb−a) . The kernel width was defined as the −3-dB (50%) cutoff period of this filter: f−1=−a2−b2+a4+14a2b2+b42π2 . For the simulations with different width kernels, a kernel width was chosen between 0.01 and 10 s, a rise time between 0.001 and 1 s, and a fall time between the rise time and 2 s. Then, a and b were chosen to minimize the squared error between these three targets using the built in MATLAB optimizer ‘fminsearch.’

White noise with a SD of 0.15 ΔFF was then added to the mock fluorescence traces.

Nonlinearity

In our linear simulations used throughout this work, the fluorescence kernels associated with a fast sequence of APs were approximated to sum linearly. In real cultured neurons, a summation nonlinearity has been observed such that sequences of APs do not generate a linear summation in ΔF/F (Dana et al., 2019). To simulate this nonlinearity, the ΔFF trace was then further transformed as: ΔFF′=sign(ΔFF)*6.2641+e−3.251Re(log10(ΔFF)).

This equation was arrived at by fitting the measured responses in Dana et al. (2019; their Fig. 2C), which can be compared with the nonlinearity used here (Extended Data Fig. 3-3A).

Deconvolution

Deconvolution was performed using the previously described FOOPSI algorithm (Vogelstein et al., 2010; Friedrich et al., 2017). The regularization coefficient was set at 0.02154, which maximized the correlation between the deconvolved trace and the true spike train in a random sample of 500 simulated traces: all other parameters were optimized for each trace. Because the example regularization coefficient provided by Friedrich et al., 2017 was 2.4, we also measured information values at 100 different values for the regularization coefficient between 0 and 3; this had little effect on the measured information (Extended Data Fig. 4-1).

KSG estimator

The previously described second KSG estimator (Kraskov et al., 2004) was used using the fifth nearest neighbor distance.

Binned estimators

The binned MI estimators were used (Timme and Lapish, 2018). The activity trace was divided into 10 bins, either evenly across the span of the activity (uniform binned) or variably so the bins contained the same number of samples (occupancy binned). Position was similarly divided into 60 bins.

Gaussian simulations

To compare the analytic approximation to our numerical method, the numerical techniques had to be applied to place cells with Gaussian rate maps. The same target information, firing rates, and behavior were used as for our original 10,000 simulations with spline rate maps. However, instead the rate map was chosen as a Gaussian with width σ=e12(−1−2IAPElog(2)−log(2π)) . For the numeric simulations, the true amount of information was calculated using a numeric integrator. The instantaneous rate was calculated using the normal distribution PDF. This was normalized and used to generate a spike train and florescence trace as above.

Bayesian decoding

The Bayesian decoder used here (Fig. 5G,H) was adapted from a previously described method (Zhang et al., 1998). Decoding was performed on the likelihood that a significant transient occurred in a time frame, trained on the first 80% of the session and tested on the last 20%. The session was divided into Δt = 0.1 s bins. The conditional likelihood that an animal is in position xi given the number of active frames during a time window ( n) is p(xi|n)=pX(xi)(∏j=1Mfi,jnj)e−Δt∑j=1Mfi,j, Where pX(xi) is the (marginal) probability that the animal is in the ith spatial bin during a time sample, fi,j is the average rate of significant frames by the jth neuron in the ith spatial bin, nj is the number of significant frames observed during the time window in neuron j, and M is the total number of neurons. The decoded position was selected as the one with maximum conditional likelihood.

Animals

Ten- to 12-week-old male C57BL/6 mice (20–30 g) were individually housed under a reverse 12/12 h light/dark cycle, all experiments were conducted during the dark phase. All experiments were approved by the Northwestern University Animal Care and Use committee.

Behavior

We used a previously described virtual reality set-up and task (Heys et al., 2014; Sheffield and Dombeck, 2015; Sheffield et al., 2017), some of the behavior sessions used here has previously appeared in these studies. Briefly, water scheduled, head fixed mice were trained to run on a cylindrical treadmill down a 3-m virtual track to receive a water (4 μl) reward at the end of the track, and were subsequently teleported to the beginning of the track after a 1.5-s delay. Behavioral sessions were included if the animal ran at least 20 laps containing a continuous 40-cm run for which the velocity was over 7 cm/s during a 5- to 30-min session.

Mouse surgery and virus injected

We performed population calcium imaging of CA1 neurons as described previously (Sheffield and Dombeck, 2015; Sheffield et al., 2017). Briefly, 30 nl of AAV1-SynFCaMP6f (University of Pennsylvania Vector Core, 1.5 × 1013 GC/ml) was injected through a small craniotomy over the right hippocampus (1.8 mm lateral, 2.3 mm caudal of bregma; 1.25 mm below the surface of the brain) under isoflurane (1–2%) anesthesia. 7 d later, a hippocampal window and head plate was implanted as described previously (Dombeck et al., 2010).

Two-photon imaging

Imaging was performed as previously described (Sheffield and Dombeck, 2015; Sheffield et al., 2017). Scanimage four was used for microscope control and acquisition (Pologruto et al., 2003). Time series movies 1024 or 512 × 256 pixels) were acquired at 50 Hz. A Digidata1440A (Molecular Devices) with Clampex 10.3 synchronized position on the linear track, reward timing, and the timing of image frames.

Image processing, region of interest (ROI) selection, and calcium transient analysis

Images were processed as previously described (Sheffield and Dombeck, 2015; Sheffield et al., 2017), with minor modifications. Briefly, rigid motion correction was performed using cross-correlation as in (Dombeck et al., 2010; Miri et al., 2011; Sheffield and Dombeck, 2015), but here using a fast Fourier transform approximation on the full video. ROIs were defined as previously described (Mukamel et al., 2009; μ = 0.6, 150 principal/independent components, SD threshold = 2.5, SD smoothing width = 1, area limits = 100–1200 pixels). ΔFF traces were generated by normalizing around the eighth percentile of a 3-s sliding window. Significant transients from both experimental and mock fluorescence traces were selected by comparing the ratio of amplitudes and durations of positive to negative going transients with a false positive rate <0.01% (Dombeck et al., 2010). Mock traces used the histograms generated from the mock gCaMP6f traces (Extended Data Fig. 2-3) or from the specific matching indicator traces (Extended Data Figs. 2-2, 3-1): experimental data histograms were built separately. All subsequent analyses were run using these significant transients.

Behavior analysis

The mean virtual track velocity was defined as the total virtual track distance covered during the session divided by the total duration of the session; slow and stop periods were included in this metric. All other analyses were restricted to long running periods, where the animal exceeded a virtual track velocity of 4 cm/s and ran continuously for at least 40 cm.

Defining place fields

Place fields were defined by first creating the spatial fluorescence intensity map (fi ) with the 300-cm track divided into 60 5-cm bins. This map was smoothed via a 3-bin boxcar. Transients identified during run periods were shuffled in order and to random intervals to create 1000 bootstrapped intensity maps. Candidate fields were defined as regions of the original fluorescence map with values >99% of the bootstrapped maps. Fields were then retained if they were between 20 and 120 cm wide: significant place cells retained at least one field that satisfied these criteria.

Results

The SMGM information metrics

Here, we review the derivation of the SMGM information metrics and the underlying assumptions. For illustrative purposes throughout this manuscript, we use the example of spatial encoding in which the firing pattern of neurons carry information about the animal’s location along a linear track; however, the derivations, equations and conclusions generalize to encoded variables over other domains and dimensionalities.

Consider a random variable X representing the positions an animal might take, with x being its value measured at one time sample. The positions are subdivided into N spatial bins, such that x can take on the values {1,2,...,N} . For our analyses, N=60 . Consider a random variable Y representing the number of APs a neuron might fire, where y is the count measured within a time sample. y can take on the values of {0,1,...,+∞} . X and Y are both discrete. If X and Y both obey the assumption that each time sample is independent (i.e., they are stationary), then the MI (I, in bits per sample) between X and Y is expressed as follows: I(X;Y)=∑i=1N∑y=0+∞pX,Y(xi,y)log2pX,Y(xi,y)pX(xi)pY(y), (1)

where pX(xi) is the (marginal) probability that the animal is in the ith spatial bin during a time sample, pY(y) is the probability that the neuron fires y APs in the time sample, and pX,Y(xi,y) is the joint probability that the neuron fires y APs and is in the ith bin. Recall that pX,Y(xi,y)=pY|X(y|xi)pX(xi) , where pY|X(y|xi) is the conditional probability that the neuron fires y APs given that the animal is in the ith spatial bin. We can thus rewrite Equation 1 as follows: I(X;Y)=∑i=1N∑y=0+∞pY|X(y|xi)pX(xi)log2pY|X(y|xi)pY(y). (2)With the further assumption that the firing of the neuron follows Poisson statistics, we can then estimate the MI as follows: let the AP rate (AP/s or Hz) in a single bin be λi , and the average across the session be λ¯ . For an arbitrarily small time window Δt , the probability that an AP occurs in that window is Pr(Y=1|x=i)=λiΔt , with the probability that an AP occurs regardless of position as Pr(Y=1)=λ¯Δt . We can thus rewrite Equation 2 as: I(X;Y)=∑i=1NλiΔtpX(xi)log2λiλ¯. (3)By integrating over 1 s (∫01I(X;Y)dΔt ), we obtain the first key SMGM metric for spatial information as measured by AP firing, which is in units of bits per second: IsÊ=∑i=1NλipX(xi)log2λiλ¯. (4)For notation, we will use a carrot (̂) to indicate an information value that is measured from experiment, the superscript (E in this case) to show the source of the data, and a subscript to show the units/formula used (bits per second in this case). Thus, IsÊ is the information measured via electrophysiology in bits per second. This metric is linearly dependent on the average firing rate of the neuron, and this dependence is often removed through normalization by the average firing rate to obtain the second key metric of spatial information as measured by AP firing, which is in units of bits/(s/Hz), or more commonly, bits per AP: IAPÊ=1λ¯∑i=1NλipX(xi)log2λiλ¯. (5)Therefore, these two key metrics of spatial information are defined completely by quantities that can be experimentally measured: the mean firing rate (λ¯) from the AP counts over the duration of the recording, the AP firing rate in the ith bin from the average rate map (λi) , and the probability that the animal is in the ith spatial bin from the normalized occupancy map (pX(xi) ). The quantity of and noise in these measurements affects the quality of the metric: in particular, undersampling because of low firing rates or low trial counts induces a substantial positive bias (Treves and Panzeri, 1995).

In the derivation of these metrics, there are two key assumptions that are violated by functional fluorescence recordings. First, the recordings do not follow Poisson statistics: instead of discrete counts of APs ( y), the functional fluorescence traces consists of a continuous relative change in fluorescence (ΔF/F ), and instead of a firing rate map (λi) measured in Hz, average intensity maps in units of ΔF/F are generated. The stationarity assumption is also violated: because of the slow decay, a time sample of the fluorescence traces depend on the previous samples. The violation of these assumptions by functional fluorescence recording will affect the precision and induce biases in the SMGM information metrics. Since these effects have not previously been addressed or quantified, we measured these biases here using a simulation study.

Building a ground truth library of 10,000 neurons with known values of information

To create a neuron with a known, ground truth information value, it was necessary to generate a continuous (i.e., infinitesimally small bins) rate map (λ(x)) matching the desired information. To do this, we first normalized the track length to 1 and assumed the animal’s occupancy map to be spatially uniform (pX(xi)=1N) . We then created an exponentiated cubic spline with five randomly positioned nodes (Fig. 1A) to build a starting continuous map of the normalized instantaneous firing rate, λλ¯(x) , with the integral normalized to 1. We calculated the ground truth amount of information in bits per AP as follows: IAPE=∫01λλ¯(x)log2(λλ¯(x))dx. (6)The locations of the five nodes were then systematically varied (see Materials and Methods) to minimize the squared error between the value calculated in Equation 6 and a target amount of information (Fig. 1A,B), in the end resulting in a mean error of 5.1 × 10−9 bits/AP and a mean absolute error of 1.5 × 10−7 bits/AP. The rate map at this convergence point was used for further analysis. This procedure was repeated to generate 10,000 mock neurons with a range of (known and ground truth) information values. Note that the value in Equation 5 cannot be higher than when all the APs arrive in one spatial bin; the rate in that bin is Nλ¯ . If we assume uniform occupancy (pX(xi)=1N) , then the maximum measurable information is log2N , in our case, 5.9 bits/AP with N=60 bins. Thus, the information values considered here range between 0 and 6 bits/AP (Fig. 1C). We chose a mean firing rate (λ¯) for the neurons between 0.1 and 30 Hz, a range observed for a variety of different cortical and hippocampal neurons during behavior (Shafi et al., 2007; DeWeese et al., 2008; O’Connor et al., 2010; Roxin et al., 2011; Buzsáki and Mizuseki, 2014). From Equations 4, 5, the ground truth information in bits per second is IsE=λ¯IAPE . IsE for these choices resulted in ground truth information values between IsE=0 and IsE=24 bits per second (Fig. 1C). Example low (IAPE=0.04 bits/AP) and mid (IAPE=2 bits/AP) rate maps are shown in Figure 1D.

These rate maps provided a basis for generating mock AP firing data (and functional fluorescence data, see below). Under real experimental conditions, recording duration and bin sizes are finite and animal occupancy maps (pX(xi) ) are not spatially uniform. These experimental limitations add error to the estimate of a neuron’s ground truth information value. Therefore, to accurately re-create these limitations in our simulation study, we used real behavior datasets from head-restrained mice running along a 3-m virtual linear track for water rewards (acquired as in Sheffield and Dombeck, 2015; Sheffield et al., 2017). Unless otherwise indicated, all values reported will be the mean±standarddeviation . We selected at random from a library of 574 behavior sessions from mice navigating along familiar tracks and concatenated and truncated these sessions to create behavior sessions uniformly sampled up to 60 min in duration (average 30.2±17.1 min), resulting in an average 132±71.2 laps per session and an average running speed of 19.3 + 3.87 cm/s (Fig. 1E1). This behavior, the average firing rate (λ¯) , and the normalized rate map (λλ¯(x) ) from the mock neurons were used to create an instantaneous firing rate trace (Fig. 1E2), sampled at 1 kHz, from which AP times were generated assuming Poisson firing statistics (Fig. 1E3). An example mock of spiking in response to behavior for low (0.04 bits/AP) and mid (2 bits/AP) information neurons can be seen in Figure 1F–H. From these spiking responses, we then generated mock fluorescence traces by convolving the raster with a double-exponential kernel matching the rise and fall times for GCaMP6f (Chen et al., 2013; Fig. 1E4) and adding random Gaussian noise to model shot noise. Mock fluorescence traces for the two example neurons in Figure 1F–H can be seen in Figure 1I,J. The mock AP and fluorescence traces were used to create session mean spatial maps, of binned firing rate (λi in Hz) and change in fluorescence (fi in ΔF/F ), for information analyses (Fig. 1K). By repeating this process, we built a large dataset of spiking and fluorescence traces, generated from our library of mock neurons with known amounts of information and using real animal spatial behavior. With tens of thousands of these mock neuron recordings, we could then assess the effects of many simulation parameters on the information values determined from the metrics including firing rate, session duration, fluorescence kernel shape, and ground truth information value.

Quantification of the accuracy and precision of the SMGM bits per second metric using functional fluorescence recordings

We first applied the SMGM bits per second metric (IsÊ ) to our mock AP recording traces to verify that they can recover our ground-truth information values given finite recording durations and bin sizes, and non-uniform animal occupancy maps (pX(xi) ). Figure 2A shows three mock neurons with ground truth information values of IsE=3 , 15, and 23 bits per second. When the SMGM bits per second metric (IsÊ ) was applied to the AP traces from these example neurons, the information was well recovered, with IsÊ=2.8 , 15, and 24 bits per second, respectively. The results from these examples also held across the full 10,000 mock neuron library (Figs. 2B–D), as a linear fit (y-intercept=0.093±0.040 , intercept p = 4.6 × 10−6 bits per second and slope=0.97±0.0030 , slope p ≪ 0.01) explained nearly all the variance (R2 = 0.97), the average error was 0.22±1.25 bits per second (1.0±0.69% error) and the absolute error was 0.64±1.05 bits per second (8.4±0.69% error). There is a substantial positive bias for the lowest firing rates and smallest number of trials (Extended Data Fig. 2-1A,B) which has been previously well characterized (Treves and Panzeri, 1995), with average errors exceeding +10% for <6 min of recording, mean rate under 0.6 Hz, and under 11 trials. Thus, the SMGM bits per second metric (IsÊ ) recovers the ground-truth information well using AP recordings, with the only error coming from finite recording time and variable animal behavior.

We next discuss the changes to the SMGM bits per second metric (IsÊ ) commonly used for application to functional fluorescence traces (Hainmueller and Bartos, 2018; Heys and Dombeck, 2018), and explore the implications of these changes. Most simply, the mean firing rate (λ¯) and the mean firing rate in a spatial bin (λi ) are replaced by the mean change in fluorescence (f¯) and the mean change in fluorescence in a bin (fi) . Making these substitutions in Equation 4 results in the information as measured by functional fluorescence: IsF̂=∑i=1NfipX(xi)log2fif¯. (7)The fluorescence map fi differs from the firing rate map λi in two ways. First, the fluorescence map is approximated by the firing rate map scaled by a factor c, dependent on the height and width of the kernel and measured in units of ΔFFHz ; and second, it is smoothed by the kernel (Fig. 1E4). If we discount the latter for a moment and focus on the scaling, f≈cλ , we can see that substituting λ with cλ in Equation 4 results in IsF̂=cIsÊ . The units for IsF̂ are no longer in bits per second, as it has previously been reported (Hainmueller and Bartos, 2018), but are instead in units of bitsΔFFsecHz or bitsΔFFAP , which are difficult to interpret (see below, Guidelines for application of information metrics to functional fluorescence imaging data). The effect of smoothing is difficult to analytically since it both alters c by changing the average intensity and distorts the firing rate map. Therefore, to fully quantify the impact of convolving an AP recording with a functional fluorescence kernel on recovering ground truth information, we used our mock fluorescence traces.

We applied a GCaMP6f modeled kernel to the 10,000 mock AP traces to generate 10,000 mock fluorescence calcium traces. Figure 2A shows the fluorescence traces generated from three mock neurons with ground truth information values of IsE=3 , 15, and 23 bits per second. The effects of the convolution can be seen in the differences in scaling and shape between the fluorescence maps fi and the firing rate maps λi . When the fluorescence metric (IsF̂ ) was applied to the fluorescence traces from these example neurons, the information recovered was IsF̂=0.13 , 0.47, and 1.1 bitsΔFFAP , respectively, indicating significant deviation from the ground truth information values assuming the units are comparable. The results from these examples also held across the full 10,000 mock neuron library (Fig. 2E–G), as there was a clear scaling of the ground truth information and a consistent underestimation with a mean error of −11.1±6.7 AU (−96.0±1.3% error). The best-fit line of the measured information (IsF̂ ) versus the ground truth information (IsE) had an intercept near 0 (0.0029±0.0016 bitsΔFFAP , p = 0.07). The slope of this fit was 0.039±12e-4 ΔFFHz (p ≪ 0.01), which provides a measure of the scaling factor ( c). This error was not corrected for with denser sampling: it remained consistent even at high firing rates and many trials (Extended Data Fig. 2-1C,D). In addition to this scaling effect caused by c, smoothing of the rate map could induce nonlinearity in the relationship between and IsE . To test for such an effect, we fit the measured information in Figure 2E with a saturating exponential and compared the fits using a likelihood ratio test: the exponential did not significantly improve the fit (χ12=0.093 , p = 0.76), which indicates that smoothing by the kernel does not induce significant nonlinearities. c is dependent on the height and width (the integral) of the kernel and was measured here as 0.039 ± 12e-4 ΔFFHz . The consistent, negative bias observed in estimating information with IsF̂ (Fig. 2E) would be easy to correct for assuming the c factor, and therefore the kernel, were similar across all measured neurons. This point is considered further below in the Guidelines for application of information metrics to functional fluorescence imaging data. We conclude that ground truth information, as measured by the fluorescence SMGM bits per second metric (IsF̂ ), is transformed into different units and is linearly scaled by a factor c dependent on the height and width of the kernel.

The amplitude (height) of the change in fluorescence can vary across indicators and conditions. The height of the kernel, given a constant kernel width, should linearly scale c and the error in estimating information with IsF̂ . To explicitly test this prediction, we simulated an additional 5000 fluorescence traces with kernels of varying height (0–3 ΔF/F ; Fig. 2I–K), but that maintain the same shape and width (from the GCaMP6f kernel), and then measured the percent error in estimating information with IsF̂ . As observed above for the GCaMP6f example (Fig. 2E–G), the percent error in estimating information with IsF̂ shows little dependence on ground truth information (Fig. 2I–K). However, as a function of the height of the kernel, the percent error (averaged over all ground truth information values) in estimating information with IsF̂ is fit well with an increasing linear function (intercept = −99.8 ± 0.42%, intercept p ≪ 0.01, slope 20.7 ± 0.14%/ΔFF , slope p ≪ 0.01, R2=0.80 ; Fig. 2I,J). Over the wide array of available functional fluorescent indicators in use today (Fig. 2H), this leads to differences in error because of differences in transient height of the indicator used alone. For the indicators shown in Figure 2H, there is an average height of 0.603 ± 0.10 SD ΔF/F : the error spans from −95.8% for the kernel height reported for gCamp6f (0.19 ΔF/F ) to −88.2% for gCamp7f (0.56 ΔF/F ). It should be noted that fluorescence (ΔF/F ) is always reported here as a fractional change, not as a percentage (% ΔF/F ); if a kernel height of 19% ΔF/F is used, the units would again change. Thus, as expected, the percent error in estimating information with the SMGM bits per second estimator (IsF̂ ) scales linearly with the height of the kernel.

The width of the kernel can vary widely across fluorescent indicators (Fig. 2H), with “faster” indicators boasting shorter rise and fall times. The combined effect of a longer rise and fall time is to smooth and delay the AP train; in other words, it acts as a causal low-pass filter. The cutoff period of this low pass filter provides a measurement of the effective width of the kernel (see Materials and Methods). The effect of such differences in kernel shape on the error in estimating information with IsF̂ is difficult to measure analytically. We therefore simulated an additional 5000 fluorescence traces with kernels of different kernel widths (but constant height of the GCaMP6f kernel), resulting in a range of kernel durations (rise times: 1 ms to 1 s, fall times longer than the rise time up to 2 s), and then we measured the percent error in estimating information with IsF̂ . Similarly, as observed above for the GCaMP6f and varying kernel height examples (Fig. 2E–K), the percent error in estimating information with IsF̂ shows little dependence on ground truth information (Fig. 2N). Interestingly, the percent error (averaged over all ground truth information values) in estimating information with IsF̂ shows a complex nonlinear response as a function of the width of the kernel (Fig. 2L,M). The error increases up to a kernel width of ∼3 s, at which point it saturates at approximately −85% error. This arises from an interaction between changing the average value of the original AP trace and flattening the average fluorescence map (fi ). Over the wide array of available functional fluorescent indicators in use today, this leads to differences in error because of differences in width of the indicator used alone. For example, an average error of −97.1 ± 0.63% were observed for iGluSnfR, the shortest indicator considered here at 0.52 s. For gCamp6s, the slowest indicator examined (2.54 s), the average was −89.6 ± 4.6%. To estimate the percent errors for these five indicators considering differences in both height and duration, we used these five kernels to generate mock fluorescence traces from the 10,000 neurons in Figure 2B–G. The resulting distributions, estimated c values, and mean and absolute errors can be seen in Extended Data Figure 2-2. In summary, we conclude that information, as measured by the fluorescence SMGM bits per second metric (IsF̂ ), is transformed into different units and is linearly scaled by a factor ( c) dependent on the height and width of the kernel, with c linearly dependent on height and nonlinearly dependent on width. The error induced by these transformations changes substantially over the range of kernel values of the different functional indicators widely used today, and therefore these are important factors to consider when designing and interpreting functional imaging experiments (for further discussion, see below, Guidelines for application of information metrics to functional fluorescence imaging data).

Quantification of the accuracy and precision of the SMGM bits per AP metric using functional fluorescence recordings

The SMGM metric is commonly normalized by the mean rate to obtain a measurement in units of bits per AP. We thus applied the SMGM bits per AP metric (IAPÊ ) to our mock AP recording traces to verify that they can recover our ground-truth information values. Figure 3A shows three mock neurons with ground truth information values IAPE=0.05 , 1.8, and 4.2 bits/AP. When the SMGM bits per AP metric (IAPÊ ) was applied to the AP traces from these example neurons, the information was well recovered, with IAPÊ=0.06 , 1.8, and 4.2 bits/AP, respectively. The results from these examples also held across the full 10,000 mock neuron library (Fig. 3B–D), as a linear fit (y-intercept=0.087±0.029 , intercept p = 2.8e-184 bits per second and slope = 0.93 ± 0.0010, slope p ≪ 0.01) explained nearly all the variance (R2 = 0.99), the average error was −0.071 ± 0.23 bits/AP (3.2 ± 5.9% error) and the absolute error was 0.13 ± 0.21 bits per second (8.1 ± 9% error). However, the data were better fit with a saturating exponential (χ12=1.6e3 , p ≪ 0.01) converging to 5.8 bits/AP as it approached the limit because of the finite bin count. There is a substantial positive bias for the lowest firing rates and smallest number of trials (Extended Data Fig. 2-1E,F) which has been previously well characterized (Treves and Panzeri, 1995). Thus, the SMGM bits per-AP metric (IAPÊ ) recovers the ground-truth information well using AP recordings (except at the largest ground truth information values), with the primary error coming from finite recording time and variable animal behavior.

We next discuss the changes needed to apply the SMGM bits per AP metric (IAPÊ ) to functional fluorescence traces and explore the implications of these changes. Most simply, the mean firing rate (λ¯) and the mean firing rate in a spatial bin (λi ) are replaced by the mean change in fluorescence (f¯) and the mean change in fluorescence in a bin (fi) . Making these substitutions in Equation 5 results in the information as measured by functional fluorescence: IAPF̂=1f¯∑i=1NfipX(xi)log2fif¯. (8)

As discussed above, the fluorescence map (fi ) can be approximated as a scaled version of the rate, that is, f=cλ and f¯=cλ¯ . Thus, under this approximation, the c factors in Equation 8 cancel, leading to IAPF̂ equivalent to IAPÊ , with the same units of bits/AP. This, of course, ignores the fact that the kernel smooths the rate map, leading to a bias in the metric that is difficult to quantify analytically.

We then applied the fluorescence SMGM bits per AP metric (IAPF̂ ) to our 10,000 mock GCaMP6f traces. Figure 3A shows the fluorescence traces generated from three mock neurons with ground truth information of IAPE=0.05 , 1.8, and 4.2 bits/AP. When the fluorescence metric (IAPF̂ ) was applied to the fluorescence traces in these examples, the information recovered was 0.04, 1.8, and 3.5 bits/AP, indicating some deviations, especially for the highest information neuron. These results held for the 10,000 mock neuron library (Fig. 3E–G). At low information values, there was little bias, but at higher information values the information recovered was substantially lower than the ground truth information. The mean resulting error was −0.38 ± 0.58 bits/AP (−9.7 ± 27.8%) and absolute error of 0.39 (12.9 ± 26.4%). This error was better fit with a saturating exponential than a linear fit (χ12=1.6e3 , p ≪ 0.01), with the average error <5% up to ground truth information of 1.8 bits/AP and <10% up to 3.0 bits/AP. At ground truth information values higher than 3 bits/AP, the average error was −1.06 ± 0.595 (−22.5 ± 9.44%) and absolute error was 1.07 ± 0.589 bits/AP (22.6 ± 9.21%). This error persisted even with denser sampling: it remained consistent even at high firing rates and many trials (Extended Data Fig. 2-1E,F). Thus, the indicator induces relatively little error at lower information values (<3 bits/AP), but the smoothing effect of the kernel induces a nonlinear, negative bias to the estimator, particularly at ground truth information values over 3 bits/AP.

Although the height of the kernel can vary between different functional fluorescence indicators (Fig. 2H), these height variations linearly scale the fluorescence map. Thus, since IAPF̂ involves normalization by the mean change in fluorescence (f¯) , IAPF̂ should not depend on kernel height. To explicitly test this prediction, we used the 5000 fluorescence traces described in the previous section (Quantification of the accuracy and precision of the SMGM bits per second metric using functional fluorescence recordings), with kernels of varying height (0–3 ΔF/F), but that maintain the same shape and width (from the GCaMP6f kernel). Then, we measured the percent error in estimating information with IAPF̂ (Fig. 3H–J). Unlike for the SMGM bits per second metric, the percent error (averaged over all ground truth information values) in estimating information with IAPF̂ shows little or no dependence on the height of the kernel (p = 0.43), but a nonlinear dependence on ground truth information as in Figure 3E–G, with no significant difference in the parameters of the saturating exponential fit (χ22=1.67 , p = 0.43). Thus, as expected, the percent error in estimating information with the SMGM bits per AP metric (IAPF̂ ) does not vary with the height of the kernel.

With little effect of kernel height on IAPF̂ , the width of the kernel likely drives biases in the metric. We thus used the 5000 fluorescence traces generated from a range of different kernel durations (rise times: 1 ms to 1 s, fall times longer than the rise time up to 2 s, but constant height of the GCaMP6f kernel) from the previous section (Quantification of the accuracy and precision of the SMGM bits per second metric using functional fluorescence recordings), and then we measured the percent error in estimating information with IAPF̂ (Fig. 3K–M). Similarly, as observed above for GCaMP6f and the varying kernel height examples (Fig. 3E–J), the percent error in estimating information with IsF̂ shows a nonlinear dependence on ground truth information (Fig. 3E–G). The percent error (averaged over all ground truth information values) showed a nonlinear response as a function of the width of the kernel (Fig. 3K,L), with a steep increase in error for kernel widths >∼1 s. Even for kernel widths <∼1 s, the percent error was strongly dependent on the ground truth information value, with steep increases in error for values more than ∼2.5–3 bits/AP (Fig. 3M). Thus, as the kernels gets wider, there is more negative bias at lower and lower information measured. The resulting errors are thus larger for wider kernel indicators, for example, with a kernel width the same as gCaMP6s (2.54s), the error exceeds −17% even at low (<0.25 bits/AP) information, with average errors of −0.86 ± 1.0 bits/AP (−31 ± 19% error) and absolute errors of 0.87 ± 1.0 bits/AP (32.6 ± 16% error). In contrast, with a kernel width the same as iGluSnfR (0.52 s), the average error exceeded 5% at 3 bits/AP and 10% at 3.7 bits/AP with a mean error of −0.57±1.00 (−8.0 ± 14%) bits/AP and absolute error of 0.41 ± 0.56 bits/AP (11 ± 11%). To estimate the percent errors for the five indicators shown in Figure 2H, taking into account differences in both height and duration, we used the five kernels to generate mock fluorescence traces from the 10k neurons in Figure 3B–G. The resulting distributions, mean and absolute errors, and error thresholds can be seen in Extended Data Figure 3-1.

Since the known information values in our library of 10,000 mock neurons were determined using the SMGM metric, which includes the assumption that neuron firing follows an inhomogeneous Poisson process, we next investigated whether the biases observed between AP and fluorescence metrics (ISÊ vs ISF̂ and IAPÊ vs IAPF̂ ; Figs. 2, 3) in our mock neuron datasets were also observed in real neuron recordings (i.e., real spiking that could deviate from Poisson firing).We therefore measured information in a real spiking dataset from hippocampal neurons in rats running on a behavioral track (Chen et al., 2016; Grosmark and Buzsáki, 2016; Grossmark et al., 2016). We generated mock fluorescence traces as we did with simulated AP trains from our mock neurons, compared the information measured from APs versus fluorescence (ISÊ vs ISF̂ and IAPÊ vs IAPF̂ ) in the real neuron recordings and found that the biases were largely consistent with the simulated mock neuron datasets (Extended Data Fig. 5-1).

In summary, we conclude that ground truth information, as measured by the fluorescence SMGM bits per AP metric (IAPF̂ ), retains the units and insensitivity to height scaling of the electrophysiological metric (IAPÊ ), but is nonlinearly biased by the smoothing of the fluorescence map dictated by the width of the kernel. The estimation errors strongly depended on both the width of the kernel and the information value being measured. Since these parameters change substantially over the different functional indicators and different neuron types and behaviors that are commonly used today, they are important factors to consider when designing and interpreting functional imaging experiments (see below for further discussion).

Nonlinearity introduces further biases

The results presented in the previous two sections rely on the approximation that ΔF/F scales linearly with the firing rate, which is not strictly true in practice (Dana et al., 2018; Greenberg et al., 2018; Éltes et al., 2019). Calcium imaging can be more responsive to bursts of APs rather than isolated spikes, and saturates at high firing rates. As an example for how to examine how nonlinearities between ΔF/F and firing rate could affect the fluorescence SMGM metrics (ISF̂ and IAPF̂ ), we applied a log-sigmoid nonlinearity (Extended Data Fig. 3-3A) to the 10,000 mock GCaMP6f time-series traces described above, based on the real behavior of GCaMP6f in cultured neurons (Dana et al., 2019; see Materials and Methods). While the resulting measurements (Extended Data Fig. 3-3) of ground truth information, as measured by the fluorescence SMGM metrics, are largely consistent with the results observed when using the linear assumption (Figs. 2, 3), some quantitative difference can be seen. Thus, even a relatively simple nonlinearity between ΔF/F and firing rate can add distortions to the amount of information measured using the fluorescence SMGM approach.

Deconvolution may not be sufficient to eliminate biases

The framework presented here for comparing ground truth information with information measured with the SMGM metrics can be extended to test the efficacy of other strategies for extracting MI. In particular, a perfect AP inference method would alleviate the problems associated with applying the SMGM metrics to functional fluorescence recordings. To test the utility of such a strategy in measuring information, we applied a popular deconvolution algorithm, FOOPSI (Vogelstein et al., 2010; see Materials and Methods), to the same 10,000 mock GCaMP6f time-series traces described above. Importantly, this deconvolution algorithm (and other available algorithms) does not recover traces of relative spike probability or exact spikes times, but instead produces sparse traces with arbitrary units, that have non-zero values estimating the relative “intensity” of spike production over time ( d). This signal can be thought of as a scaled estimate of the number of spikes per time bin, and thus the average intensity map will have some similar properties to the florescence intensity maps, that is, we would expect the intensity maps from deconvolution to approximate the relative firing rate scaled by some factor c , which has arbitrary units.

We then measured information in these deconvolved d-traces using the SMGM metrics (Isd̂ and IAPd̂ ), which are identical to (IsF̂ and IAPF̂ ), except SMGM is applied to d-traces instead of functional fluorescence traces. When using the SMGM bits per second measure (Isd̂ ; Fig. 4A), we found a clear scaling of the ground truth information. The scaling factor was very small (c=1.15e-3±1.8e-5 AU), resulting in low predicted information (mean error −11.4 ± 6.92 AU, mean % error −99.8 ± 0.38%). This error was larger than when we measured information directly from the florescence traces using IsF̂ (Fig. 2E; 96.0% absolute error, rank-sum p ≪ 0.01; c = 0.0390). It is worth noting that the deconvolved trace d can be arbitrarily scaled, so in a sense this error is arbitrary. However, these are the results from the scaling chosen by a widely used deconvolution algorithm and the large error emphasize that the scale of d can have a large effect on the bits per second measure (Isd̂ ).

Assuming that the intensity map of the deconvolved d-traces are a scaled version of the true rate maps, we could measure information using the SMGM bits per AP metric IAPd̂ without changing units (Fig. 4B). Compared with the SMGM bits per AP metric applied to florescence (IAPF ), on average there was some reduction in the nonlinearity at higher ground truth information (IAPE ) values when using IAPd̂ , resulting in linear fits closer to the unitary line (IAPd̂slope=1.02±0.0027,R2=0.76 vs IAPF̂ slope = 0.78 ± 0.0011, R2 = 0.93). However, information measured with IAPd̂ was still better fit with a saturating exponential (χ12 = 2.3e3, p ∼ 0) converging to a saturation value of 5.51 bits/AP (compared with 5.78 for IAPF̂ ), as expected since the algorithm is not expected to resolve spikes at orders of magnitude shorter timescales than the kernel. This resulted in a positive bias at lower levels of ground truth information. For ground truth information values below 3 bits/AP, the average error for IAPd̂ was 0.556 ± 0.50 bits/AP (68.9 ± 104%) compared with −0.060 ± 0.13 bits/AP (2.75 ± 0.36%) for IAPF̂ . For ground truth information values above 3 bits/AP, the average error for IAPd̂ was −0.127 ± 0.76 bits/AP (−1.0 ± 16.8%) as compared with −1.04 ± 0.59 bits/AP (−0.22 ± 9.4%) for IAPF̂ . Overall, there was more error when the SMGM bits/AP metric was applied to deconvolved data compared with when applied directly to fluorescence traces [IAPd̂ mean absolute error 0.60 ± 0.47 bits/AP (52.7 ± 89.1%) versus IAPF̂ was 0.13 ± 0.21 bits per second (8.1 ± 9% error), rank-sum p ≪ 0.01]. Thus, when comparing the recovery of ground truth information from functional fluorescence traces using either direct application of the SMGM metrics (IsF̂ and IAPF̂ ) or the application of the SMGM metrics to deconvolved z-traces (Isd̂ and IAPd̂ ), we found better recovery using the direct application approach (IsF̂ and IAPF̂ ).

The KSG and binned estimators are poor estimators of MI in functional florescence data

In addition to SMGM, the KSG and binned estimation metrics have been developed for estimating MI between variables. These other two metrics produce information measured in bits per second, so they are only comparable to the SMGM bits per second estimator (IsF̂ ). The KSG metric uses the kth nearest neighbor distances between points in the neural response and behavioral variable space to estimate information (Kraskov et al., 2004). The binned estimation metric uses discrete bins to estimate the full multidimensional joint probability distribution (p(X,Y) in Eq. 1) to estimate MI (Timme and Lapish, 2018). These two metrics estimate information across time samples and therefore are dependent on firing rate like the SMGM bits per second metric considered above. Further, the binned estimator is sensitive to the precise method used for data binning, and thus, we have used two commonly applied binning methods: uniform and occupancy-based bins.

We applied the KSG, binned estimator (uniform bins), and binned estimator (occupancy binned; Materials and Methods) to the same 10,000 mock GCaMP6f time-series traces and behavioral data used to assess the SMGM approach. These methods all behaved similarly when applied to our simulations (Fig. 4C–E), so they will be discussed together here. The information values measured by these techniques correlated with ground truth information in bits per second (ISE ). Interestingly, unlike the SMGM bits per second metric (IsF̂ ; Fig. 2), the KSG and binned estimator results were better fit with a saturating exponential than with a linear fit (χ12=1.70e4 , 3.12e4, and 3.14e4, respectively, p ∼ 0). The KSG and binned estimator methods overestimated the information at lower ground truth (IsE ) values and saturated quickly at higher values. For ground truth information (IsE ) values below 10 bits per second, the mean absolute errors were 10.3 ± 7.82 bits per second (465 ± 4390%), 22.2 ± 12.6 bits per second (821.7 ± 3353%), and 24.2 ± 14.0 bits per second (911 ± 3926%) for the KSG, binned estimator (uniform bins), and binned estimator (occupancy binned), respectively (Fig. 4F). This is in comparison to the 0.35 ± 0.59 bitsΔFFAP (14 ± 103%) mean absolute error found using the SMGM bits per second metric (IsF̂ ) for ground truth information values <10 bits per second. For ground truth information (IsE ) values >10 bits per second, the mean absolute errors were 13.2 ± 8.66 bits per second (85.0 ± 64.4%), 25.4 ± 16.2 bits per second (165.2 ± 118.3%), and 28.9 ± 17.34 bits per second (186.3 ± 124.9%) for the KSG, binned estimator (uniform bins), and binned estimator (occupancy binned), respectively. This is in comparison to the 0.97 ± 1.15 bits per second (5.85 ± 6.73%) mean absolute errors found using the SMGM bits per second metric (IsF̂ ) for ground truth information values >10 bits per second. Over the full range of ground truth values, mean absolute errors of 11.9 ± 8.42 bits per second (254.8%), 23.5 ± 14.8 bits per second (456.7%) and 26.7 ± 16.1 bits per second (509.0%) were found for the KSG, binned estimator (uniform bins), and binned estimator (occupancy binned), respectively, an order magnitude larger than the 2.4 ± 2.97 AU (22 ± 27% error) error seen using the SMGM bits per second metric (IsF̂ ). As a control, we applied the binned estimator s to AP traces and compared the estimated information to the ground truth information to verify that the large errors observed (Fig. 4D,E) were caused when the estimators were applied to fluorescence data (rather than simply a difference between the binned estimator, which do not rely on a Poisson firing assumption, and the ground truth information established using SMGM, which does rely on a Poisson firing assumption). We found the errors when applying the binned estimators to AP traces were relatively small [mean absolute error 2.72 ± 3.38 bits per second (41.1%) and 2.70 ± 3.33 bits per second (41.0%) for the uniform and occupancy-based binning; respectively; Extended Data Fig. 4-2]. Therefore, when comparing the recovery of ground truth information from functional fluorescence traces using either the SMGM metric IsF̂ or the KSG and binned estimator metrics, we found better recovery using the SMGM approach (IsF̂ ).

An analytic approximation can reproduce some qualitative, but not quantitative, results of the numeric solutions

Some of the general features of the relationship between ground truth information and fluorescence SMGM metrics can also be seen using an analytic approximation. For example, if we approximate the rate map as a Gaussian firing field with mean rate λ¯ , simplify the kernel approximation to a single exponential with falloff τ, assume constant, normalized movement speed v, and assume that the convolution between the kernel and the mean rate map is nearly Gaussian, we can approximate the relationship between the bits per second ground truth information ISE and measured fluorescence information ISF as ISF≈−Avτλ¯log(4−ISEλ¯+2πev2τ2)log(4), (9)and between the bits per AP ground truth information IAPE and measured fluorescence information IAPF as IAPF≈−log(4−IAPE+2eπv2τ2)log(4). (10)

Similar to the numerical solution, the analytic approximation provided by these equations (Eqs. 9, 10) predict that the fluorescence bits per second metric is dominated by a prefactor (Aτvλ− in the analytical case), and that the fluorescence bits per AP metric saturates at larger information values. Our numerical solutions provide more accurate measures for the magnitude of these effects, and for the magnitude of information values themselves, given that they include the more accurate double exponential kernel, signal noise, and the realistic nonstationary speed, position and fluorescence signals. These quantitative differences can be seen in Extended Data Figure 4-3, where we directly compared this analytic approximation to our numerical approach by simulating 10,000 neurons with Gaussian rate maps (λ(x) ) with known (ground truth) information. We found a significant difference in the slope of the bits per second estimator (c, 0.041 vs 0.0036 ΔF/FHz for the numeric vs analytical, respectively), likely because of the non-stationarities present in behavior and florescence signals. For the fluorescence bits per AP measure, the analytic approximation predicts a large positive bias for ground truth values up to ∼4 bits/AP. This is in contrast to the numeric solution, which has <10% error for ground truth information values below 3.12 bits/AP.

Guidelines for application of information metrics to functional fluorescence imaging data

Taken together, the above results suggest that across the information metrics applied directly to functional fluorescence traces, the SMGM metrics provide the most reliable and interpretable information measurements. We thus suggest the following guidelines for use and interpretation of the SMGM metrics as applied to fluorescence MI metrics (IsF̂ and IAPF̂ ) defined in Equations 7, 8.

The SMGM bits per second metric (IsF̂ ) is likely attractive to imaging researchers because the units suggest that precise knowledge of AP numbers and times are not required for its use. However, there are several challenges when applying the SMGM bits per second metric to functional fluorescence imaging data. First, the substitution of the change in fluorescence map ( f) for the AP firing rate map ( λ) introduces a change in units, from bits per second to ΔFFbitsAP , which is difficult to interpret and relate back to bits per second. Second, the transformation of AP firing rate to change in fluorescence can be approximated by a c scaling factor (f=cλ ), which is measured in ΔFFHz , a quantity that is unknown a priori. If c is not consistent between the neurons of a population of interest, then the information values will be scaled differently and cannot be directly compared (Fig. 2). Since c is dependent on the width and height of the indicator response to a single AP (the kernel), it can vary from neuron to neuron based on difference in indicator expression level, intracellular calcium buffering, and many other factors (Park and Dunlap, 1998; Aponte et al., 2008; Helmchen and Tank, 2015; Greenberg et al., 2018). More research will be needed to measure these parameters (Chen et al., 2012, 2013), and thus c, across neurons. Some results suggest that there may be non-trivial amounts of variability within a populations of neurons (Greenberg et al., 2018; Éltes et al., 2019). With the impact of c on the SMGM bits per second metric, and the possible variability of c across a population of neurons, how can researchers properly extract useful measurements of information using IsF̂ ?

Guideline 1

First, we note that if experimental measurements reveal small and acceptable variations in c across the neurons of interest, then the information values derived from IsF̂ can be normalized by this factor to recover information values in units of SMGM bits per second (independent of ΔFF ) that can be compared across neurons.

Under the assumption of a consistent kernel, approximations for c for common indicators can be found in Extended Data Figure 2-2.

Guideline 2

Further, given small variations in c across the neurons of interest, the ratio of IsF̂ between neurons in the population provides a meaningful metric for comparisons. For example, such ratios could be used to divide a population of neurons accurately into groups based on their information values (e.g., three quantiles of information) or compare the information values between different functional subtypes of neurons (e.g., between place and non-place cells).

Guideline 3

The metric can still be useful even if experimental measurements reveal large and unacceptable variations in c across the neurons of interest, or if experimental measurements of c do not exist. In such cases, since it is reasonable to assume that c is consistent in the same neuron over time, comparisons across the same neuron can provide meaningful insights by using a ratio of IsF̂ measured (from the same neuron) across different conditions. For example, quantifying the neuron-by-neuron ratio of IsF̂ between different behavioral states or conditions of an animal or task, such as between goal-directed versus non-goal-directed running down a linear track, would allow researchers to make conclusions such as the following: “The population of neurons in region Z carries X+/−Y times more information during goal-directed than non-goal-directed running.”

Therefore, we conclude that with careful consideration of the (known or unknown) variability of the fluorescence response kernel ( c), IsF̂ can be used to extract useful measurements of information, either direct measurements of information across a population of neurons (with known and similar c), ratios of information between different neurons of a population (with known and similar c) or differences across different conditions within the same neuron (with c unknown or different across neurons).

The SMGM fluorescence bits per AP metric (IAPF̂ ) results in the same units as the AP-based metric IAPÊ , and therefore may provide imaging researchers with information values that are relatively easy to interpret. However, this similarity in units is somewhat misleading since the number and timing of APs are not directly measured with functional fluorescence traces and the asymmetric and relatively slow dynamics of fluorescence indicators leads to shifting and smoothing of the AP rate map ( λ). This issue can have the effect of inducing a significant negative bias in information measurements, especially at high information values and with functional indicators with wider kernels (Fig. 3). This is the most important factor to consider when determining how researchers can properly extract useful measurements of information using IAPF̂ . The shifting and smoothing of the AP rate map by fluorescence effectively leads to crosstalk between adjacent spatial bins. Therefore, it is critical to consider the size of the spatial bins in relation to the spatial shift and smoothing induced by the indicator (effectively the kernel plotted in space, rather than time, using the animal’s average running velocity to transform from time to space). It is reasonable then to counteract the spatial shift and smoothing effect by using larger bin sizes, but this only works up to a point since larger bins limit the maximum amount of information possible to measure and may negatively bias the information values near this upper limit, even for AP-based recordings (Fig. 3B–D). Researchers could potentially optimize the recovery of ground truth information by appropriately selecting bin size for a particular indicator (see Extended Data Fig. 3-2A,B).

In practice, using gCaMP6f and the rodent spatial behavior and spatial bin sizes (5 cm) used here, our analysis suggests that IAPF̂ provides reasonable measurements of information for neurons with values up to 3 bits/AP (Fig. 3E–G), since this is the point where the absolute error exceeds 10% [comparable to the mean absolute error when measuring information from AP data (8.4%)]. Equivalent thresholds for other common indicators are shown in Extended Data Fig. 3-1. The error is exacerbated by slower indicators and thus more accurate measurements of information will result from using the fastest, narrowest kernel indicators available, assuming signal-to-noise and detection efficiency are comparable across the different width indicators.

Guideline 4

We conclude that with careful consideration of the size of the spatial bins in relation to the spatial shift and smoothing induced by the indicator, IAPF̂ can be used to extract useful measurements of information, most accurately for neurons with <3 bits/AP under recording conditions similar to those considered here.

Previous research quantifying information in bits per AP using IAPÊ have found that the majority of neurons carry information in this range (<3 bits/AP; Knierim et al., 1995; Markus et al., 1995; Lee et al., 2006; Poucet and Sargolini, 2013; Bourboulou et al., 2019), with a few exceptions (Ji and Wilson, 2007). Although these levels are dependent on the number of bins and bin dwell time, IAPF̂ should be widely applicable to quantifying information throughout the brain during behavior.

Example: application of information metrics to functional fluorescence imaging data from hippocampus during spatial behavior

In this section, we demonstrate use of the above guidelines for proper application and interpretation of the SMGM fluorescence MI metrics (IsF̂ and IAPF̂ ) defined in Equations 7, 8. We applied these metrics to functional fluorescence recordings (gCaMP6f) from pyramidal neurons in CA1 of the hippocampus acquired during mouse spatial behavior.

CA1 neurons expressing gCaMP6f (viral transfection, Camk2a promoter) were imaged with two-photon microscopy through a chronic imaging window during mouse navigation along a familiar 1D virtual linear track, as described previously (Fig. 5A,B; Dombeck et al., 2010; Sheffield et al., 2017; Radvansky and Dombeck, 2018). Eight fields of view from four mice were recorded in eight total sessions (recording duration 8.8 ± 1.3 min, number of traversals/session: 29 ± 2.5, 3.6 ± 0.3 laps/min, 3-m-long track). From these eight sessions, 1500 neurons were identified from our segmentation algorithm (see Materials and Methods), and analysis was restricted to the 964 neurons that displayed at least one calcium transient on at least 1/3 of the traversals during the session. Among these 964 neurons, 304 (31.5%) had significant place fields and were thus identified as place cells (see Materials and Methods), while the remaining 660 (68.5%) did not pass a place field test and were thus identified as non-place cells.

By applying Equation 7 (using 5-cm sized spatial bins), we found a continuum of spatial information values measured by the fluorescence SMGM bits per second metric (IsF̂ ) across the 964 CA1 neurons, with an average value of IsF̂=0.14±0.0040 bitsΔFFHz (Fig. 5C). The units for IsF̂ make direct use and interpretation of these values difficult, however, because these recordings were all from pyramidal neurons in a single area, here for illustrative purposes, we will presume that variations in c (discussed above) across the 964 neurons of interest are small and acceptable, with the absolute value of c unknown. This allows for comparisons of information ratios across different subsets of the population. For example, place cells had 2.8±0.20 times more information than non-place cells using the SMGM bits per second metric (IsF̂=0.18±0.0047 for place vs 0.063 ± 0.0029 bitsΔFFHz for non-place, rank-sum p = 1.7e-63; Fig. 5D–F), although there was substantial overlap in information between the populations (see distributions in Fig. 5D and individual examples in Fig. 5E,F). This also allows for accurate division of the 964 neurons into three quantiles based on information values, which we use below for spatial location decoding.

By applying Equation 8 (using 5-cm sized spatial bins), we found a continuum of spatial information values measured by the fluorescence SMGM bits per AP metric (IAPF̂ ) across the 964 CA1 neurons, with an average value of IAPF̂=1.65±0.023 bits/AP (Fig. 5C). The units for IAPF̂ allow for direct use and interpretation of these values, and notably, because most (97%) of the neurons had values <3 bits/AP, a mean absolute error of <10% can be assumed across the distribution of SMGM bits per AP values. When applied to the place and non-place cell populations, we found that place cells had higher information than nonplace cells using the fluorescence SMGM bits per AP metric (IAPF̂=1.8±0.026 and 1.35 ± 0.042 bits/AP, rank-sum p = 4.6e-21; Fig. 5D–F). This is consistent with mock fluorescence traces generated from real neuron AP datasets (Extended Data Fig. 5-1B).

As a demonstration of the usefulness of using information metrics to analyze large functional fluorescence population recordings, we explored the accuracy of decoding the animal’s track position using different subsets of neurons. We divided the 964 neurons into nine groups: all neurons, place cells, non-place cells, three quantiles based on the fluorescence SMGM bits per second metric, and three quantiles based on the fluorescence SMGM bits per AP metric. We then used a Bayesian decoder of the animals’ position (see Materials and Methods) separately for each of the nine neuron groups in each of the eight sessions (Fig. 5G,H). An individual session decoding example can be seen in Figure 5G. We quantified decoding accuracy using the absolute position decoding error (% of track), and pooled this measure across sessions for each neuron group (Fig. 5H). The means and standard errors for each group are: all neurons (7.33 ± 2.5%), place cells (6.97 ± 1.9%), non-place cells (20.9 ± 1.8%), SMGM bits per second Q1 (21.9 ± 1.5%), SMGM bits per second Q2 (13.2 ± 2.4%), SMGM bits per second Q3 (8.97 ± 2.4%), SMGM bits per AP Q1 (17.6 ± 2.7%), SMGM bits per AP Q2 (17.8 ± 3.1%), SMGM bits per AP Q3 (10.4 ± 3.0%). Interestingly, even the lowest quantile information groups still could be used to determine animal track location to within ∼1/5 of the track. This supports the idea that the hippocampal code for space is carried by a large population of active neurons (Meshulam et al., 2016), and not just by a select subpopulation with the highest information or most well-defined tuning curves. As could be expected, place cells encoded the position of the animal better than nonplace cells and better than the lowest quantile information groups (Holm–Bonferroni corrected rank-sum, α=0.05 ), and neurons in the higher quantiles provided more accurate decoding. Thus, the fluorescence information metrics provide a means to compare the relative contribution of hippocampal neurons with different information values to decoding animal position.

Discussion

Here, we performed an in-depth simulation study to examine the application of the SMGM bits per second and SMGM bits per AP metrics of MI to functional fluorescence recordings. Since these metrics were designed for AP recordings and since functional fluorescence recordings violate some of the assumptions that these metrics are based on, it was unclear whether and how the metrics could be used for functional fluorescence recordings. We created a library of ten thousand mock neurons whose AP output carried ground-truth amounts of information about the animal’s spatial location, and by using real behavioral recording data from mice navigating in virtual linear tracks, we simulated the spatial firing patterns of the mock neurons. We then simulated fluorescent calcium responses for each neuron in each session by convolving the AP trains with calcium kernels for different indicators, primarily GCamp6f (although see Extended Data Figs. 2-2, 3-1 for results from other indicators), and then added noise.

We then derived fluorescence versions of the SMGM bits per second (IsF̂ ) and SMGM bits per AP metrics (IAPF̂ ; Eqs. 7, 8) and applied them to the fluorescence traces to quantify the performance of the metrics for estimating information. We found that ground truth information, as measured by the fluorescence SMGM bits per second metric (IsF̂ ), was transformed into different units and was linearly scaled by a factor ( c) dependent on the height and width of the kernel, with c linearly dependent on height and nonlinearly dependent on width. The error induced by these transformations changed substantially over the range of kernel values of the different functional indicators widely used today, and therefore are important factors to consider when designing and interpreting functional imaging experiments. We then found that ground truth information, as measured by the fluorescence SMGM bits per AP metric (IAPF̂ ), retains the units and insensitivity to height scaling of the electrophysiological metric (IAPÊ ), but is nonlinearly biased by the smoothing of the fluorescence map dictated by the width of the kernel. The estimation errors strongly depended on both the width of the kernel and the information value being measured. Importantly, since these parameters change substantially over the different functional indicators and different neuron types and behaviors that are commonly used today, they are important factors to consider when designing and interpreting functional imaging experiments. For example, even for the same indicator, the shape of the kernel is a function of intracellular calcium buffering, indicator concentration, the amount of calcium influx, the efflux rates, background fluorescence and resting calcium concentration, which can all vary across different cells. Additionally, the results presented here rely on the approximation that ΔF/F scales linearly with the firing rate, which is not strictly true in practice. We show in Extended Data Figure 3-3 that even a relatively simple nonlinearity between ΔF/F and firing rate can add distortions to the amount of information measured using the fluorescence SMGM approach. This relationship between ΔF/F and firing rate can vary across different indicators and, since the Toolbox can be used to vary this relationship, users can further explore this source of bias.

In our approach, the known information values in our library of 10,000 mock neurons were determined using the SMGM metric, which includes the assumption that neuron firing follows an inhomogeneous Poisson process. It is important to remember that the SMGM metric, which has been applied to spiking data extensively over the past few decades, requires the use of a Poisson estimate of spiking probability, i.e., the Poisson assumption is built into the original metric. In practice, even spiking data violates this and other assumptions of the SMGM metric since real neurons do not strictly follow Poisson statistics (for example, they can display neural hysteresis) and animal behavior is non-stationary. Here, we are building from this existing framework and adding and testing whether it is possible to apply the metric to functional fluorescence datasets. Even still, the Poisson assumption could have contributed to some of the biases found when evaluating the fluorescence SMGM metrics with respect to ground truth information. We explored this potential source of bias further using two different analyses. First, in Extended Data Figure 4-2, we applied the binned estimators (which do not rely on a Poisson firing assumption) to AP traces and compared the estimated information to the ground truth information (which was established using the SMGM metric that does rely on a Poisson firing assumption). We found the errors to be relatively small, particularly in comparison to the errors induced by the binned estimators when applied to fluorescence traces (Fig. 4D,E). Second, in Extended Data Figure 5-1, we used a real spiking dataset from hippocampal neurons in mice running on a behavioral track (i.e., real spiking neurons that can deviate from Poisson firing) and generated mock fluorescence traces from the AP traces. When we compared the information measured from the AP traces to the fluorescence traces, we found biases that were largely consistent with those observed in Figures 2, 3 from our simulated mock neuron datasets. Taken together, these analyses indicate that any biases resulting from the Poisson assumption in the simulation procedure appear to be small, particularly with respect to the biases introduced when AP traces are transformed into functional fluorescence traces. Finally, in the Toolbox, we also include code to generate mock neurons using a binned distribution, avoiding the Poisson assumption of SMGM. Thus, users can further explore sources of bias using a different ground truth dataset.

Using our mock fluorescence traces, we also asked whether an AP estimation method could relieve the biases in the SMGM metrics. Applying the SMGM bits per second metric (Isd̂ ) to AP estimation traces from a deconvolution algorithm (FOOPSI) resulted in a low c value for recovered versus ground truth information. When the SMGM bits per AP measure was applied (IAPd̂ ), the resulting measurements of information were still nonlinear (compared with IsF̂ ), with a positive bias at lower values of ground truth information. Overall, applying FOOPSI to fluorescence traces led to a poorer recovery of ground truth information using SMGM compared with direct application of SMGM to the florescence traces (IAPF̂ ). Importantly, this result from deconvolution is only specific to GCaMP6f, and conclusions should not be drawn about other indicators or situations; users will be able to use the Toolbox to explore this area further. We also tested other metrics to measure MI directly from the fluorescence time traces [KSG, binned estimator (uniform bins), and binned estimator (occupancy binned)] and found these alternatives produced highly variable, saturating measurements of recovered versus ground truth information. This was in contrast to the SMGM bits per second measure (IsF̂ ) which produced a linearly scaled bias with lower error.

Taken together, we find that the SMGM bits per AP metric can well recover the MI between spiking and behavior. The SMGM bits per second metric is scaled such that comparisons should be limited to within populations of well characterized neurons or for within neuron comparisons, e.g., ratios of information across conditions. In general, researchers should use caution when applying measures developed for AP data in fluorescence recordings: there’s no guarantee that the assumptions that support the measures hold for fluorescence data, and this can lead to difficult to interpret and biased results.

Acknowledgments

Acknowledgements: We thank B. Kath and members of the Dombeck Lab for comments on this manuscript and V. Jayaraman, R. Kerr, D. Kim, L. Looger, and K. Svoboda from the GENIE Project (Janelia, Howard Hughes Medical Institute) for GCaMP6.

Footnotes

  • The authors declare no competing financial interests.

  • This work was supported by The McKnight Foundation, Northwestern University, The Chicago Biomedical Consortium with support from the Searle Funds at The Chicago Community Trust, National Institutes of Health Grants R01MH101297 and T32AG020506, and the National Science Foundation Grant CRCNS1516235.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    Adoff MD, Climer JR, Davoudi H, Marvin JS, Looger LL, Dombeck DA (2021) The functional organization of excitatory synaptic input to place cells. Nat Commun 12:3558.
    OpenUrl
  2. ↵
    Ahrens MB, Orger MB, Robson DN, Li JM, Keller PJ (2013) Whole-brain functional imaging at cellular resolution using light-sheet microscopy. Nat Methods 10:413–420. doi:10.1038/nmeth.2434 pmid:23524393
    OpenUrlCrossRefPubMed
  3. ↵
    Allen TA, Salz DM, McKenzie S, Fortin NJ (2016) Nonspatial sequence coding in CA1 neurons. J Neurosci 36:1547–1563. doi:10.1523/JNEUROSCI.2874-15.2016 pmid:26843637
    OpenUrlAbstract/FREE Full Text
  4. ↵
    Aponte Y, Bischofberger J, Jonas P (2008) Efficient Ca2+ buffering in fast-spiking basket cells of rat hippocampus. J Physiol 586:2061–2075. doi:10.1113/jphysiol.2007.147298 pmid:18276734
    OpenUrlCrossRefPubMed
  5. ↵
    Aronov D, Tank DW (2014) Engagement of neural circuits underlying 2D spatial navigation in a rodent virtual reality system. Neuron 84:442–456. doi:10.1016/j.neuron.2014.08.042 pmid:25374363
    OpenUrlCrossRefPubMed
  6. ↵
    Belghazi MI, Baratin A, Rajeswar S, Ozair S, Bengio Y, Courville A, Hjelm RD (2018) MINE: mutual information neural estimation. 35th Int Conf Mach Learn ICML 2018 2:864–873.
    OpenUrl
  7. ↵
    Bourboulou R, Marti G, Michon FX, El Feghaly E, Nouguier M, Robbe D, Koenig J, Epsztein J (2019) Dynamic control of hippocampal spatial coding resolution by local visual cues. Elife 8:e44487. doi:10.7554/eLife.44487
    OpenUrlCrossRef
  8. ↵
    Brandon MP, Bogaard AR, Libby CP, Connerney MA, Gupta K, Hasselmo ME (2011) Reduction of theta rhythm dissociates grid cell spatial periodicity from directional tuning. Science 332:595–599. doi:10.1126/science.1201652 pmid:21527714
    OpenUrlAbstract/FREE Full Text
  9. ↵
    Brown EN, Barbieri R, Eden UT, Frank L (2003) Likelihood methods for neural spike train data analysis. In: Computational neuroscience: a comprehensive approach (Feng J, ed), pp 252–281. Boca Raton: CRC Press.
  10. ↵
    Buzsáki G, Mizuseki K (2014) The log-dynamic brain: how skewed distributions affect network operations. Nat Rev Neurosci 15:264–278. doi:10.1038/nrn3687 pmid:24569488
    OpenUrlCrossRefPubMed
  11. ↵
    Calton JL, Stackman RW, Goodridge JP, Archey WB, Dudchenko PA, Taube JS (2003) Hippocampal place cell instability after lesions of the head direction cell network. J Neurosci 23:9719–9731. pmid:14585999
    OpenUrlAbstract/FREE Full Text
  12. ↵
    Chen Q, Cichon J, Wang W, Qiu L, Lee SJR, Campbell NR, DeStefino N, Goard MJ, Fu Z, Yasuda R, Looger LL, Arenkiel BR, Gan WB, Feng G (2012) Imaging neural activity using Thy1-GCaMP transgenic mice. Neuron 76:297–308. doi:10.1016/j.neuron.2012.07.011 pmid:23083733
    OpenUrlCrossRefPubMed
  13. ↵
    Chen TW, Wardill TJ, Sun Y, Pulver SR, Renninger SL, Baohan A, Schreiter ER, Kerr RA, Orger MB, Jayaraman V, Looger LL, Svoboda K, Kim DS (2013) Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499:295–300. doi:10.1038/nature12354 pmid:23868258
    OpenUrlCrossRefPubMed
  14. ↵
    Chen Z, Grosmark AD, Penagos H, Wilson MA (2016) Uncovering representations of sleep-associated hippocampal ensemble spike activity. Sci Rep 6:32193. doi:10.1038/srep32193 pmid:27573200
    OpenUrlCrossRefPubMed
  15. ↵
    Climer JR, Newman EL, Hasselmo ME (2013) Phase coding by grid cells in unconstrained environments: two-dimensional phase precession. Eur J Neurosci 38:2526–2541. doi:10.1111/ejn.12256 pmid:23718553
    OpenUrlCrossRefPubMed
  16. ↵
    Climer JR, DiTullio R, Newman EL, Hasselmo ME, Eden UT (2015) Examination of rhythmicity of extracellularly recorded neurons in the entorhinal cortex. Hippocampus 25:460–473. doi:10.1002/hipo.22383 pmid:25331248
    OpenUrlCrossRefPubMed
  17. ↵
    Cohen MR, Kohn A (2011) Measuring and interpreting neuronal correlations. Nat Neurosci 14:811–819. doi:10.1038/nn.2842 pmid:21709677
    OpenUrlCrossRefPubMed
  18. ↵
    Dana H, Sun Y, Mohar B, Hulse B, Hasseman JP, Tsegaye G, Tsang A, Wong A, Patel R, Macklin JJ, Chen Y, Konnerth A, Jayaraman V, Looger LL, Schreiter ER, Svoboda K, Kim DS, H, Dana Y, Sun B, Mohar B, et al. (2018) High-performance GFP-based calcium indicators for imaging activity in neuronal populations and microcompartments. bioRxiv 16:434589.
    OpenUrlCrossRef
  19. ↵
    Dana H, Sun Y, Mohar B, Hulse BK, Kerlin AM, Hasseman JP, Tsegaye G, Tsang A, Wong A, Patel R, Macklin JJ, Chen Y, Konnerth A, Jayaraman V, Looger LL, Schreiter ER, Svoboda K, Kim DS (2019) High-performance calcium sensors for imaging activity in neuronal populations and microcompartments. Nat Methods 16:649–657. doi:10.1038/s41592-019-0435-6 pmid:31209382
    OpenUrlCrossRefPubMed
  20. ↵
    DeWeese MR, Zador AM, Hromádka T (2008) Sparse representation of sounds in the unanesthetized auditory cortex. PLoS Biol 6:e16. doi:10.1371/journal.pbio.0060016
    OpenUrlCrossRefPubMed
  21. ↵
    Dombeck DA, Khabbaz AN, Collman F, Adelman TL, Tank DW (2007) Imaging large-scale neural activity with cellular resolution in awake, mobile mice. Neuron 56:43–57. doi:10.1016/j.neuron.2007.08.003
    OpenUrlCrossRefPubMed
  22. ↵
    Dombeck DA, Harvey CD, Tian L, Looger LL, Tank DW (2010) Functional imaging of hippocampal place cells at cellular resolution during virtual navigation. Nat Neurosci 13:1433–1440. doi:10.1038/nn.2648 pmid:20890294
    OpenUrlCrossRefPubMed
  23. ↵
    Éltes T, Szoboszlay M, Kerti‐Szigeti K, Nusser Z (2019) Improved spike inference accuracy by estimating the peak amplitude of unitary [Ca2+] transients in weakly GCaMP6f‐expressing hippocampal pyramidal cells. J Physiol 597:2925–2947. doi:10.1113/JP277681 pmid:31006863
    OpenUrlCrossRefPubMed
  24. ↵
    Evans MH, Petersen RS, Humphries MD (2019) On the use of calcium deconvolution algorithms in practical contexts. bioRxiv 871137.
  25. ↵
    Frank LM, Brown EN, Wilson MA (2001) A comparison of the firing properties of putative excitatory and inhibitory neurons from CA1 and the entorhinal cortex. J Neurophysiol 86:2029–2040. doi:10.1152/jn.2001.86.4.2029 pmid:11600659
    OpenUrlCrossRefPubMed
  26. ↵
    Friedrich J, Zhou P, Paninski L (2017) Fast online deconvolution of calcium imaging data. PLoS Comput Biol 13:e1005423. doi:10.1371/journal.pcbi.1005423 pmid:28291787
    OpenUrlCrossRefPubMed
  27. ↵
    Fu H, Rodriguez GA, Herman M, Emrani S, Nahmani E, Barrett G, Figueroa HY, Goldberg E, Hussaini SA, Duff KE (2017) Tau pathology induces excitatory neuron loss, grid cell dysfunction, and spatial memory deficits reminiscent of early Alzheimer’s disease. Neuron 93:533–541.e5. doi:10.1016/j.neuron.2016.12.023 pmid:28111080
    OpenUrlCrossRefPubMed
  28. ↵
    Fyhn M, Molden S, Hollup S, Moser M-B, Moser EI (2002) Hippocampal neurons responding to first-time dislocation of a target object. Neuron 35:555–566. doi:10.1016/S0896-6273(02)00784-5 pmid:12165476
    OpenUrlCrossRefPubMed
  29. ↵
    Gao W, Kannan S, Oh S, Viswanath P (2017) Estimating mutual information for discrete-continuous mixtures. Advances in neural information processing systems, pp 5987–5998. Curren Associates, Inc: Long Beach California.
  30. ↵
    Gerrard JL, Burke SN, McNaughton BL, Barnes CA (2008) Sequence reactivation in the hippocampus is impaired in aged rats. J Neurosci 28:7883–7890. doi:10.1523/JNEUROSCI.1265-08.2008 pmid:18667620
    OpenUrlAbstract/FREE Full Text
  31. ↵
    Greenberg DS, Wallace DJ, Voit K-M, Wuertenberger S, Czubayko U, Monsees A, Handa T, Vogelstein JT, Seifert R, Groemping Y, Kerr JN (2018) Accurate action potential inference from a calcium sensor protein through biophysical modeling. bioRxiv 479055.
  32. ↵
    Grosmark AD, Buzsáki G (2016) Diversity in neural firing dynamics supports both rigid and learned hippocampal sequences. Science 351:1440–1443. doi:10.1126/science.aad1935 pmid:27013730
    OpenUrlAbstract/FREE Full Text
  33. ↵
    Grossmark AD, Long J, Buzsáki G (2016) Recordings from hippocampal area CA1, PRE, during and POST novel spatial learning. doi: 10.6080/K0862DC5.doi:10.6080/K0862DC5
    OpenUrlCrossRef
  34. ↵
    Grubb MS, Thompson ID (2003) Quantitative characterization of visual response properties in the mouse dorsal lateral geniculate nucleus. J Neurophysiol 90:3594–3607. doi:10.1152/jn.00699.2003 pmid:12944530
    OpenUrlCrossRefPubMed
  35. ↵
    Hainmueller T, Bartos M (2018) Parallel emergence of stable and dynamic memory engrams in the hippocampus. Nature 558:292–296. doi:10.1038/s41586-018-0191-2 pmid:29875406
    OpenUrlCrossRefPubMed
  36. ↵
    Harvey CD, Collman F, Dombeck DA, Tank DW (2009) Intracellular dynamics of hippocampal place cells during virtual navigation. Nature 461:941–946. doi:10.1038/nature08499 pmid:19829374
    OpenUrlCrossRefPubMed
  37. ↵
    Hazama Y, Tamura R (2019) Effects of self-locomotion on the activity of place cells in the hippocampus of a freely behaving monkey. Neurosci Lett 701:32–37. doi:10.1016/j.neulet.2019.02.009 pmid:30738872
    OpenUrlCrossRefPubMed
  38. ↵
    Helmchen F, Tank DW (2015) A single-compartment model of calcium dynamics in nerve terminals and dendrites. Cold Spring Harb Protoc 2015:155–167. doi:10.1101/pdb.top085910 pmid:25646507
    OpenUrlCrossRefPubMed
  39. ↵
    Heys JG, Dombeck DA (2018) Evidence for a subcircuit in medial entorhinal cortex representing elapsed time during immobility. Nat Neurosci 21:1574–1582. doi:10.1038/s41593-018-0252-8 pmid:30349104
    OpenUrlCrossRefPubMed
  40. ↵
    Heys JG, Rangarajan KV, Dombeck DA (2014) The functional micro-organization of grid cells revealed by cellular-resolution imaging. Neuron 84:1079–1090. doi:10.1016/j.neuron.2014.10.048 pmid:25467986
    OpenUrlCrossRefPubMed
  41. ↵
    Hinman JR, Brandon MP, Climer JR, Chapman GW, Hasselmo ME (2016) Multiple running speed signals in medial entorhinal cortex. Neuron 91:666–679. doi:10.1016/j.neuron.2016.06.027 pmid:27427460
    OpenUrlCrossRefPubMed
  42. ↵
    Hok V, Chah E, Save E, Poucet B (2013) Prefrontal cortex focally modulates hippocampal place cell firing patterns. J Neurosci 33:3443–3451. doi:10.1523/JNEUROSCI.3427-12.2013 pmid:23426672
    OpenUrlAbstract/FREE Full Text
  43. ↵
    Huang Y, Brandon MP, Griffin AL, Hasselmo ME, Eden UT (2009) Decoding movement trajectories through a T-maze using point process filters applied to place field data from rat hippocampal region CA1. Neural Comput 21:3305–3334. doi:10.1162/neco.2009.10-08-893 pmid:19764871
    OpenUrlCrossRefPubMed
  44. ↵
    Hubel DH, Wiesel TN (2009) Republication of The J of Physiology (1959) 148, 574-591: receptive fields of single neurones in the cat’s striate cortex. 1959. J Physiol 587:2721–2732. doi:10.1113/jphysiol.2009.174151 pmid:19525558
    OpenUrlCrossRefPubMed
  45. ↵
    Ji D, Wilson MA (2007) Coordinated memory replay in the visual cortex and hippocampus during sleep. Nat Neurosci 10:100–107. doi:10.1038/nn1825 pmid:17173043
    OpenUrlCrossRefPubMed
  46. ↵
    Jing M, Zhang P, Wang G, Feng J, Mesik L, Zeng J, Jiang H, Wang S, Looby JC, Guagliardo NA, Langma LW, Lu J, Zuo Y, Talmage DA, Role LW, Barrett PQ, Zhang LI, Luo M, Song Y, Zhu JJ, et al. (2018a) Dendritic inhibition in the hippocampus supports fear learning. Nat Biotechnol 36:726–737. doi:10.1038/nbt.4184
    OpenUrlCrossRefPubMed
  47. ↵
    Jing M, Zhang P, Wang G, Feng J, Mesik L, Zeng J, Jiang H, Wang S, Looby JC, Guagliardo NA, Langma LW, Lu J, Zuo Y, Talmage DA, Role LW, Barrett PQ, Zhang LI, Luo M, Song Y, Zhu JJ, et al. (2018b) State-dependent subnetworks of parvalbumin-expressing interneurons in neocortex. Nat Biotechnol 36:726–737. doi:10.1038/nbt.4184
    OpenUrlCrossRefPubMed
  48. ↵
    Jing M, Zhang P, Wang G, Feng J, Mesik L, Zeng J, Jiang H, Wang S, Looby JC, Guagliardo NA, Langma LW, Lu J, Zuo Y, Talmage DA, Role LW, Barrett PQ, Zhang LI, Luo M, Song Y, Zhu JJ, et al. (2018c) Calcium dynamics in dendrites of hippocampal CA1 interneurons in awake mice. Nat Biotechnol 36:726–737. doi:10.1038/nbt.4184
    OpenUrlCrossRefPubMed
  49. ↵
    Jing M, Zhang P, Wang G, Feng J, Mesik L, Zeng J, Jiang H, Wang S, Looby JC, Guagliardo NA, Langma LW, Lu J, Zuo Y, Talmage DA, Role LW, Barrett PQ, Zhang LI, Luo M, Song Y, Zhu JJ, et al. (2018d) A genetically encoded fluorescent acetylcholine indicator for in vitro and in vivo studies. Nat Biotechnol 36:726–737. doi:10.1038/nbt.4184 pmid:29985477
    OpenUrlCrossRefPubMed
  50. ↵
    Kalko EV, Dukas R, Ratcliffe JM, Teeling EC, Haven N, Fattu JM, Bates ME, Simmons JA, Riquimaroux H, Surlykke A, Bouffard FH, Lee DN, Dear SP, Horiuchi TK, Krishnaprasad PS, Moss CF, Schuller G, Brudzynski SM, Syme DA, Hollingworth S, et al. (2011) An expanded palette of genetically encoded Ca2+ indicators. Science 333:1888–1891.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    Keller PJ, Ahrens MB (2015) Visualizing whole-brain activity and development at the single-cell level using light-sheet microscopy. Neuron 85:462–483. doi:10.1016/j.neuron.2014.12.039 pmid:25654253
    OpenUrlCrossRefPubMed
  52. ↵
    Khoshkhoo S, Vogt D, Sohal VS (2017) Dynamic, cell-type-specific roles for GABAergic interneurons in a mouse model of optogenetically inducible seizures. Neuron 93:291–298. doi:10.1016/j.neuron.2016.11.043 pmid:28041880
    OpenUrlCrossRefPubMed
  53. ↵
    Kinsky NR, Sullivan DW, Mau W, Hasselmo ME, Eichenbaum HB (2018) Hippocampal place fields maintain a coherent and flexible map across long timescales. Curr Biol 28:3578–3588.e6. doi:10.1016/j.cub.2018.09.037 pmid:30393037
    OpenUrlCrossRefPubMed
  54. ↵
    Knierim JJ, Kudrimoti HS, McNaughton BL (1995) Place cells, head direction cells, and the learning of landmark stability. J Neurosci 15:1648–1659. pmid:7891125
    OpenUrlAbstract/FREE Full Text
  55. ↵
    Koenig J, Linder AN, Leutgeb JK, Leutgeb S (2011) The spatial periodicity of grid cells is not sustained during reduced theta oscillations. Science 332:592–595. doi:10.1126/science.1201685 pmid:21527713
    OpenUrlAbstract/FREE Full Text
  56. ↵
    Kraskov A, Stögbauer H, Grassberger P (2004) Estimating mutual information. Phys Rev E Stat Physics Plasmas Fluids Relat Interdiscip Top 69:16.
    OpenUrl
  57. ↵
    Kraus BJ, Brandon MP, Robinson RJ, Connerney MA, Hasselmo ME, Eichenbaum HB (2015) During running in place, grid cells integrate elapsed time and distance run. Neuron 88:578–589. doi:10.1016/j.neuron.2015.09.031 pmid:26539893
    OpenUrlCrossRefPubMed
  58. ↵
    Kropff E, Carmichael JE, Moser M-B, Moser EI (2015) Speed cells in the medial entorhinal cortex. Nature 523:419–424. doi:10.1038/nature14622 pmid:26176924
    OpenUrlCrossRefPubMed
  59. ↵
    Lee I, Griffin AL, Zilli EA, Eichenbaum HB, Hasselmo ME (2006) Gradual translocation of spatial correlates of neuronal firing in the hippocampus toward prospective reward locations. Neuron 51:639–650. doi:10.1016/j.neuron.2006.06.033 pmid:16950161
    OpenUrlCrossRefPubMed
  60. ↵
    Leveteau J, MacLeod P (1966) Olfactory discrimination in the rabbit olfactory glomerulus. Science 153:175–176. doi:10.1126/science.153.3732.175 pmid:5940356
    OpenUrlAbstract/FREE Full Text
  61. ↵
    Liu P, Jerrard LE, Bilkey DK, Jarrard LE, Bilkey DK (2004) Excitotoxic lesions of the pre- and parasubiculum disrupt the place fields of hippocampal pyramidal cells. Hippocampus 14:107–116. doi:10.1002/hipo.10161 pmid:15058488
    OpenUrlCrossRefPubMed
  62. ↵
    MacDonald CJ, Carrow S, Place R, Eichenbaum HB (2013) Distinct hippocampal time cell sequences represent odor memories in immobilized rats. J Neurosci 33:14607–14616. doi:10.1523/JNEUROSCI.1537-13.2013 pmid:24005311
    OpenUrlAbstract/FREE Full Text
  63. ↵
    Mankin EA, Thurley K, Chenani A, Haas OV, Debs L, Henke J, Galinato M, Leutgeb JK, Leutgeb S, Leibold C (2019) The hippocampal code for space in Mongolian gerbils. Hippocampus 29:787–715. doi:10.1002/hipo.23075 pmid:30746805
    OpenUrlCrossRefPubMed
  64. ↵
    Mann K, Gallen CL, Clandinin TR (2017) Whole-brain calcium imaging reveals an intrinsic functional network in Drosophila. Curr Biol 27:2389–2396.e4. doi:10.1016/j.cub.2017.06.076 pmid:28756955
    OpenUrlCrossRefPubMed
  65. ↵
    Markus EJ, Qin YL, Leonard B, Skaggs WEW, McNaughton BBL, Barnes CA (1995) Interactions between location and task affect the spatial and directional firing of hippocampal neurons. J Neurosci 15:7079–7094. doi:10.1523/JNEUROSCI.15-11-07079.1995 pmid:7472463
    OpenUrlAbstract/FREE Full Text
  66. ↵
    Marvin JS, Scholl B, Wilson DE, Podgorski K, Kazemipour A, Müller JA, Schoch S, Quiroz FJU, Rebola N, Bao H, Little JP, Tkachuk AN, Cai E, Hantman AW, Wang SSH, DePiero VJ, Borghuis BG, Chapman ER, Dietrich D, DiGregorio DA, et al. (2018) Stability, affinity, and chromatic variants of the glutamate sensor iGluSnFR. Nat Methods 15:936–939. doi:10.1038/s41592-018-0171-3 pmid:30377363
    OpenUrlCrossRefPubMed
  67. ↵
    Marvin JS, Shimoda Y, Magloire V, Leite M, Kawashima T, Jensen TP, Kolb I, Knott EL, Novak O, Podgorski K, Leidenheimer NJ, Rusakov DA, Ahrens MB, Kullmann DM, Looger LL (2019) A genetically encoded fluorescent sensor for in vivo imaging of GABA. Nat Methods 16:763–770. doi:10.1038/s41592-019-0471-2 pmid:31308547
    OpenUrlCrossRefPubMed
  68. ↵
    Mau W, Sullivan DW, Kinsky NR, Hasselmo ME, Howard MW, Eichenbaum H (2018) The same hippocampal CA1 population simultaneously codes temporal information over multiple timescales. Curr Biol 28:1499–1508.e4. doi:10.1016/j.cub.2018.03.051 pmid:29706516
    OpenUrlCrossRefPubMed
  69. ↵
    Meshulam L, Gauthier JL, Brody CD, Tank DW, Bialek W (2016) Collective behavior of place and non-place neurons in the hippocampal network. Neuron 96:1178–1192.e4.
    OpenUrl
  70. ↵
    Miri A, Daie K, Burdine RD, Aksay E, Tank DW (2011) Regression-based identification of behavior-encoding neurons during large-scale optical imaging of neural activity at cellular resolution. J Neurophysiol 105:964–980. doi:10.1152/jn.00702.2010 pmid:21084686
    OpenUrlCrossRefPubMed
  71. ↵
    Morris TP, White IR, Crowther MJ (2019) Using simulation studies to evaluate statistical methods. Stat Med 38:2074–2102. doi:10.1002/sim.8086 pmid:30652356
    OpenUrlCrossRefPubMed
  72. ↵
    Mukamel EA, Nimmerjahn A, Schnitzer MJ (2009) Automated analysis of cellular signals from large-scale calcium imaging data. Neuron 63:747–760. doi:10.1016/j.neuron.2009.08.009 pmid:19778505
    OpenUrlCrossRefPubMed
  73. ↵
    Newman EL, Climer JR, Hasselmo ME (2014) Grid cell spatial tuning reduced following systemic muscarinic receptor blockade. Hippocampus 24:643–655. doi:10.1002/hipo.22253 pmid:24493379
    OpenUrlCrossRefPubMed
  74. ↵
    Nguyen MN, Hori E, Matsumoto J, Tran AH, Ono T, Nishijo H (2013) Neuronal responses to face-like stimuli in the monkey pulvinar. Eur J Neurosci 37:35–51. doi:10.1111/ejn.12020 pmid:23121157
    OpenUrlCrossRefPubMed
  75. ↵
    Nguyen MN, Matsumoto J, Hori E, Maior RS, Tomaz C, Tran AH, Ono T, Nishijo H (2014) Neuronal responses to face-like and facial stimuli in the monkey superior colliculus. Front Behav Neurosci 8:85. doi:10.3389/fnbeh.2014.00085 pmid:24672448
    OpenUrlCrossRefPubMed
  76. ↵
    Nguyen JP, Shipley FB, Linder AN, Plummer GS, Liu M, Setru SU, Shaevitz JW, Leifer AM (2016) Whole-brain calcium imaging with cellular resolution in freely behaving Caenorhabditis elegans. Proc Natl Acad Sci USA 113:E1074–E1081. doi:10.1073/pnas.1507110112 pmid:26712014
    OpenUrlAbstract/FREE Full Text
  77. ↵
    Niell CM, Stryker MP (2008) Highly selective receptive fields in mouse visual cortex. J Neurosci 28:7520–7536. doi:10.1523/JNEUROSCI.0623-08.2008 pmid:18650330
    OpenUrlAbstract/FREE Full Text
  78. ↵
    O’Connor DH, Peron SP, Huber D, Svoboda K (2010) Neural activity in barrel cortex underlying vibrissa-based object localization in mice. Neuron 67:1048–1061. doi:10.1016/j.neuron.2010.08.026 pmid:20869600
    OpenUrlCrossRefPubMed
  79. ↵
    O’Keefe J (1976) Place units in the hippocampus of the freely moving rat. Exp Neurol 51:78–109.
    OpenUrlCrossRefPubMed
  80. ↵
    Østergaard J, Kramer MA, Eden UT (2018) Capturing spike variability in noisy Izhikevich neurons using point process generalized linear models. Neural Comput 30:125–148. doi:10.1162/neco_a_01030 pmid:29064782
    OpenUrlCrossRefPubMed
  81. ↵
    Pachitariu M, Stringer C, Harris KD (2018) Robustness of spike deconvolution for neuronal calcium imaging. J Neurosci 38:7976–7985. doi:10.1523/JNEUROSCI.3339-17.2018 pmid:30082416
    OpenUrlAbstract/FREE Full Text
  82. ↵
    Paninski L (2004) Maximum likelihood estimation of cascade point-process neural encoding models. Netw Comput Neural Syst 15:243–262. doi:10.1088/0954-898X_15_4_002
    OpenUrlCrossRef
  83. ↵
    Park D, Dunlap K (1998) Dynamic regulation of calcium influx by G-proteins, action potential waveform, and neuronal firing frequency. J Neurosci 18:6757–6766. pmid:9712647
    OpenUrlAbstract/FREE Full Text
  84. ↵
    Park E, Dvorak D, Fenton AA (2011) Ensemble place codes in hippocampus: CA1, CA3, and dentate gyrus place cells have multiple place fields in large environments. PLoS One 6:e22349. doi:10.1371/journal.pone.0022349 pmid:21789250
    OpenUrlCrossRefPubMed
  85. ↵
    Park IJ, Bobkov YV, Ache BW, Principe JC (2013) Quantifying bursting neuron activity from calcium signals using blind deconvolution. J Neurosci Methods 218:196–205. doi:10.1016/j.jneumeth.2013.05.007 pmid:23711821
    OpenUrlCrossRefPubMed
  86. ↵
    Pastalkova E, Itskov V, Amarasingham A, Buzsáki G (2008) Internally generated cell assembly sequences in the rat hippocampus. Science 321:1322–1327. doi:10.1126/science.1159775 pmid:18772431
    OpenUrlAbstract/FREE Full Text
  87. ↵
    Pologruto TA, Sabatini BL, Svoboda K (2003) ScanImage: flexible software for operating laser scanning microscopes. Biomed Eng Online 2:13. doi:10.1186/1475-925X-2-13 pmid:12801419
    OpenUrlCrossRefPubMed
  88. ↵
    Poucet B, Sargolini F (2013) A trace of your place. Science 340:35–36. doi:10.1126/science.1237567 pmid:23559237
    OpenUrlAbstract/FREE Full Text
  89. ↵
    Radvansky BA, Dombeck DA (2018) An olfactory virtual reality system for mice. Nat Commun 9:839. doi:10.1038/s41467-018-03262-4 pmid:29483530
    OpenUrlCrossRefPubMed
  90. ↵
    Rashid SK, Pedrosa V, Dufour MA, Moore JJ, Chavlis S, Delatorre RG, Poirazi P, Clopath C, Basu J (2020) The dendritic spatial code: branch-specific place tuning and its experience-dependent decoupling. bioRxiv 2020.01.24.916643.
  91. ↵
    Robbe D, Buzsáki G (2009) Alteration of theta timescale dynamics of hippocampal place cells by a cannabinoid is associated with memory impairment. J Neurosci 29:12597–12605. doi:10.1523/JNEUROSCI.2407-09.2009 pmid:19812334
    OpenUrlAbstract/FREE Full Text
  92. ↵
    Roxin A, Brunel N, Hansel D, Mongillo G, van Vreeswijk C (2011) On the distribution of firing rates in networks of cortical neurons. J Neurosci 31:16217–16226. doi:10.1523/JNEUROSCI.1677-11.2011 pmid:22072673
    OpenUrlAbstract/FREE Full Text
  93. ↵
    Scholl B, Wilson DE, Fitzpatrick D (2017) Local order within global disorder: synaptic architecture of visual space. Neuron 96:1127–1138.e4. doi:10.1016/j.neuron.2017.10.017 pmid:29103806
    OpenUrlCrossRefPubMed
  94. ↵
    Shafi M, Zhou Y, Quintana J, Chow C, Fuster J, Bodner M (2007) Variability in neuronal activity in primate cortex during working memory tasks. Neuroscience 146:1082–1108. doi:10.1016/j.neuroscience.2006.12.072 pmid:17418956
    OpenUrlCrossRefPubMed
  95. ↵
    Sheffield MEJ, Dombeck DA (2015) Calcium transient prevalence across the dendritic arbour predicts place field properties. Nature 517:200–204. doi:10.1038/nature13871 pmid:25363782
    OpenUrlCrossRefPubMed
  96. ↵
    Sheffield MEJ, Adoff MD, Dombeck DA (2017) Increased prevalence of calcium transients across the dendritic arbor during place field formation. Neuron 96:490–504.e5. doi:10.1016/j.neuron.2017.09.029 pmid:29024668
    OpenUrlCrossRefPubMed
  97. ↵
    Simonnet J, Brecht M (2019) Burst firing and spatial coding in subicular principal cells. J Neurosci 39:3651–3662. doi:10.1523/JNEUROSCI.1656-18.2019 pmid:30819796
    OpenUrlAbstract/FREE Full Text
  98. ↵
    Skaggs WE, McNaughton BL, Gothard KM (1993) An information-theoretic approach to deciphering the hippocampal code 1030–1037. In: Advances in neural information processing systems 5 (Hanson SJ, Cowan JD, Giles CL, eds). Burlington: Morgan-Kaufmann.
  99. ↵
    Soo FS, Schwartz GW, Sadeghi K, Berry MJ (2011) Fine spatial information represented in a population of retinal ganglion cells. J Neurosci 31:2145–2155. doi:10.1523/JNEUROSCI.5129-10.2011 pmid:21307251
    OpenUrlAbstract/FREE Full Text
  100. ↵
    Stackman RW, Taube JS (1998) Firing properties of rat lateral mammillary single units: head direction, head pitch, and angular head velocity. J Neurosci 18:9020–9037. pmid:9787007
    OpenUrlAbstract/FREE Full Text
  101. ↵
    Stirman JN, Smith IT, Kudenov MW, Smith SL (2016) Wide field-of-view, multi-region, two-photon imaging of neuronal activity in the mammalian brain. Nat Biotechnol 34:857–862. doi:10.1038/nbt.3594 pmid:27347754
    OpenUrlCrossRefPubMed
  102. ↵
    Stringer C, Pachitariu M, Steinmetz N, Carandini M, Harris KD (2019) High-dimensional geometry of population responses in visual cortex. Nature 571:361–365. doi:10.1038/s41586-019-1346-5 pmid:31243367
    OpenUrlCrossRefPubMed
  103. ↵
    Tang J, Ardila Jimenez SC, Chakraborty S, Schultz SR (2016) Visual receptive field properties of neurons in the mouse LGN. PLoS One 11:e0146017.
    OpenUrlCrossRefPubMed
  104. ↵
    Timme NM, Lapish C (2018) A tutorial for information theory in neuroscience. eNeuro 5:ENEURO.0052-18.2018. doi:10.1523/ENEURO.0052-18.2018
    OpenUrlAbstract/FREE Full Text
  105. ↵
    Treves A, Panzeri S (1995) The upward bias in measures of information derived from limited data samples. Neural Comput 7:399–407. doi:10.1162/neco.1995.7.2.399
    OpenUrlCrossRef
  106. ↵
    Vogelstein JT, Packer AM, Machado TA, Sippy T, Babadi B, Yuste R, Paninski L (2010) Fast nonnegative deconvolution for spike train inference from population calcium imaging. J Neurophysiol 104:3691–3704. doi:10.1152/jn.01073.2009 pmid:20554834
    OpenUrlCrossRefPubMed
  107. ↵
    Wachowiak M, Shipley MT (2006) Coding and synaptic processing of sensory information in the glomerular layer of the olfactory bulb. Semin Cell Dev Biol 17:411–423. doi:10.1016/j.semcdb.2006.04.007 pmid:16765614
    OpenUrlCrossRefPubMed
  108. ↵
    Wilent WB, Nitz DA (2007) Discrete place fields of hippocampal formation interneurons. J Neurophysiol 97:4152–4161. doi:10.1152/jn.01200.2006 pmid:17392415
    OpenUrlCrossRefPubMed
  109. ↵
    Yaksi E, Friedrich RW (2006) Reconstruction of firing rate changes across neuronal populations by temporally deconvolved Ca2+ imaging. Nat Methods 3:377–383. doi:10.1038/nmeth874 pmid:16628208
    OpenUrlCrossRefPubMed
  110. ↵
    Yartsev MM, Ulanovsky N (2013) Representation of three-dimensional space in the hippocampus of flying bats. Science 340:367–372. doi:10.1126/science.1235338 pmid:23599496
    OpenUrlAbstract/FREE Full Text
  111. ↵
    Zhang K, Ginzburg I, McNaughton BL, Sejnowski TJ (1998) Interpreting neuronal population activity by reconstruction: unified framework with application to hippocampal place cells. J Neurophysiol 79:1017–1044. doi:10.1152/jn.1998.79.2.1017 pmid:9463459
    OpenUrlCrossRefPubMed
  112. ↵
    Zhou JL, Shatskikh TN, Liu X, Holmes GL (2007) Impaired single cell firing and long-term potentiation parallels memory impairment following recurrent seizures. Eur J Neurosci 25:3667–3677. doi:10.1111/j.1460-9568.2007.05598.x pmid:17610586
    OpenUrlCrossRefPubMed
  113. ↵
    Zinyuk L, Kubik S, Kaminsky Y, Fenton AA, Bures J (2000) Understanding hippocampal activity by using purposeful behavior: place navigation induces place cell discharge in both task-relevant and task-irrelevant spatial reference frames. Proc Natl Acad Sci USA 97:3771–3776. doi:10.1073/pnas.050576397 pmid:10716713
    OpenUrlAbstract/FREE Full Text
  114. ↵
    Ziv Y, Burns LD, Cocker ED, Hamel EO, Ghosh KK, Kitch LJ, El Gamal A, Schnitzer MJ (2013) Long-term dynamics of CA1 hippocampal place codes. Nat Neurosci 16:264–266. doi:10.1038/nn.3329 pmid:23396101
    OpenUrlCrossRefPubMed

Synthesis

Reviewing Editor: Michaël Zugaro, CNRS, Collège de France, Inserm

Decisions are customarily a result of the Reviewing Editor and the peer reviewers coming together and discussing their recommendations until a consensus is reached. When revisions are invited, a fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision will be listed below. The following reviewer(s) agreed to reveal their identity: Cathie Ventalon.

This paper examines whether information metrics can be used directly on calcium imaging data, without prior spike deconvolution. The authors use extensive numerical simulations to show the conditions at which information metrics can be properly used on calcium traces. They also provide a toolbox for readers that would like to explore these questions on their own data. This represent a necessary advance in the field, since several papers have already started to apply information metrics directly on functional imaging data.

The authors have addressed all the questions raised during the previous review process.

The only remaining issue is that the toolbox does not seem to run properly (error while running the function genExpSpline). The authors should simply conduct a review of the code and fix the (possible) remaining issues before the paper can be accepted.

Author Response

Author Response:

We appreciate the reviewer(s) finding this error in our toolbox code. There was a bug in the code that has now been fixed.

Back to top

In this issue

eneuro: 8 (5)
eNeuro
Vol. 8, Issue 5
September/October 2021
  • Table of Contents
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this eNeuro 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.
Information Theoretic Approaches to Deciphering the Neural Code with Functional Fluorescence Imaging
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
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
Information Theoretic Approaches to Deciphering the Neural Code with Functional Fluorescence Imaging
Jason R. Climer, Daniel A. Dombeck
eNeuro 25 August 2021, 8 (5) ENEURO.0266-21.2021; DOI: 10.1523/ENEURO.0266-21.2021

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Information Theoretic Approaches to Deciphering the Neural Code with Functional Fluorescence Imaging
Jason R. Climer, Daniel A. Dombeck
eNeuro 25 August 2021, 8 (5) ENEURO.0266-21.2021; DOI: 10.1523/ENEURO.0266-21.2021
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
    • Synthesis
    • Author Response
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • behavior
  • calcium imaging
  • in vivo
  • information theory
  • place cells
  • two-photon microscopy

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Article: Methods/New Tools

  • A New Tool for Quantifying Mouse Facial Expressions
  • Validation of a New Coil Array Tailored for Dog Functional Magnetic Resonance Imaging Studies
  • Photothrombotic Middle Cerebral Artery Occlusion in Mice: A Novel Model of Ischemic Stroke
Show more Research Article: Methods/New Tools

Novel Tools and Methods

  • Behavioral and Functional Brain Activity Alterations Induced by TMS Coils with Different Spatial Distributions
  • Bicistronic expression of a high-performance calcium indicator and opsin for all-optical stimulation and imaging at cellular resolution
  • Synthetic Data Resource and Benchmarks for Time Cell Analysis and Detection Algorithms
Show more Novel Tools and Methods

Subjects

  • Novel Tools and Methods

  • Home
  • Alerts
  • 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
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro 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 eNeuro.