Simultaneous calcium imaging and extracellular recording from the same neuron

Calcium imaging is a powerful method to record from many neurons simultaneously. But what do the recorded signals really mean?

This question can only be properly addressed by experiments which record both calcium signals and action potentials from the same neuron (ground truth recordings). These recordings are technically quite challenging. So we assembled several existing ground truth datasets, and in addition recorded ground truth datasets ourselves, totaling >200 neuronal recordings.

This blog blog posts contains raw movies together with recorded action potentials (black; also turn on your speakers for the spikes!) and the recorded ΔF/F of the calcium recording (blue). These ground truth data are a very direct way for everybody into calcium imaging to get an intuition about what is really going on. (Scroll down if you want to see recordings in zebrafish!)

Recording from a L2/3 neuron in visual cortex with GCaMP6f, tg(Emx1), from Huang et al., bioRxiv, 2019; a very beautiful recording. Replayed with 2x speed.

Recording from a L2/3 neuron in visual cortex with GCaMP6f, tg(Emx1), from Huang et al., bioRxiv, 2019. Stronger contamination from surrounding neuropil.

Recording from a L2/3 neuron in visual cortex with GCaMP6f, tg(Emx1), from Huang et al., bioRxiv, 2019. Note that single action potentials don’t seem to have any impact at all. – The negative transients in the calcium trace stem from center-surround neuropil decontamination (activity of the surround is subtracted).

Recording from a L2/3 neuron in visual cortex with GCaMP6s, tg(Emx1), from Huang et al., bioRxiv, 2019.

Recording from a L2/3 neuron in visual cortex with GCaMP6s, tg(Emx1), from Huang et al., bioRxiv, 2019.

Recording from a L2/3 neuron in visual cortex with GCaMP6f, virally induced, from Chen et al., Nature, 2013. From the left, you can see the shadow of the patch pipette used for recording of extracellular signals.

Something completely different: recording from a pyramidal neuron in Ca3 with R-CaMP1.07, virally induced, recorded by Stefano Carta, from Rupprecht et al., bioRxiv, 2020. What appears as single events are actually bursts of 5-15 action potentials with inter-spike-intervals of <6 ms.

A recording that I performed myself in adult zebrafish, in a subpart of the homolog of olfactory cortex (aDp) with GCaMP6f, tg(neuroD), in Rupprecht et al., bioRxiv, 2020. Around second 20, it is visible that even a single action potential can be seen in the calcium signal. However,this was not always the case in other neurons that I recorded from the same brain region.

Again a recording that I did in adult zebrafish, in the dorsal part of the dorsal telencephalon with GCaMP6f, tg(neuroD), in Rupprecht et al., bioRxiv, 2020.

What can you do if you want to detect single isolated action potentials with calcium imaging? GCaMP, due to its sigmoid non-linearity, is by often a bad choice and will be strongly biased towards bursts. Synthetic indicators, however, are very linear in the low-calcium regime. – This is a recording that I did myself in adult zebrafish, in a subpart of the homolog of olfactory cortex (pDp) with the injected synthetic indicator OGB-1 in Rupprecht et al., bioRxiv, 2020. Although the temporal resolution of the calcium recording is rather low, the indicator clearly responds to single action potentials. As another asset, the indicator not only fills the cytoplasm of the neuron in a ring-like shape, which makes neuropil-contamination much less of an issue compared to GCaMPs.

Another recording that I performed in adult zebrafish, in a subpart of the homolog of olfactory cortex (pDp) with the injected synthetic indicator Cal-520 in Rupprecht et al., bioRxiv, 2020. This indicator is much more sensitive compared to OGB-1, but also diffuses less well after bolus injection. – These two minutes of recording only contain 4 spikes (this brain region really is into low firing rates in general), but you can clearly see all of them. If this were a GCaMP recording, you would probably see only a flat line throughout the entire recording.

For more information, including all 20 datasets with >200 neurons (rather than these excerpts from 11 neurons), check out the following resources:

Posted in Calcium Imaging, Data analysis, electrophysiology, machine learning, Microscopy, Neuronal activity, zebrafish | Tagged , , , , | Leave a comment

Discrepancies between calcium imaging and extracellular ephys recordings

To record the activity from a population of neurons, calcium imaging and extracellular recordings with small electrodes are the two most widely used methods that are still able to disentangle the contributions from single units. Here, I would like to briefly mention two papers that try to connect these two approaches by comparing them more or less directly.

  1. Wei et al., A comparison of neuronal population dynamics measured with calcium imaging and electrophysiology, bioRxiv, 2019
    [Update, 2020-09-15: The paper just came out in PLoS Comp Biology just one day after the blog post!]
  2. Siegle, Ledochowitsch et al., Reconciling functional differences in populations of neurons recorded with two-photon imaging and electrophysiology, bioRxiv, 2020

Both papers compare calcium imaging datasets and extracellular ephys datasets, try to connect the results and point out the difficulties in reconciling the approaches.

Wei et al. use datasets recorded in mouse anterio-lateral motor cortex (ALM). They first focus on approaches to reconstruct spike rates from calcium imaging data (deconvolution) and find several limitations of this approach. On the other hand, they find that a forward model that transforms spiking activity to calcium fluorescence data can reconcile most of these differences.

The authors also provide a user-friendly website which can be used to explore the transformations between ephys and imaging data (also including datasets with simultaneous ephys-imaging): (Understanding the figures of the paper is however quite useful before exploring the website.)

From Wei et al., bioRxiv (2019), used under CC BY 4.0 license (excerpt from Fig. 7).

A large part of the paper focuses on high-level analyses (principal component analysis, decoding of behavior and decision). I would take it as an educational tale of caution which highlights wrong conclusions that could be made based on standard analyses. For example, the slow and variable decay times of calcium imaging data can lead to a dispersion of peak activity that was absent in the ephys data. This dispersion can smear out clearly timed activations of neuronal population into something more similar to a sequence (see Figure 7, of which I have pasted an excerpt above).

Siegle, Ledochowitsch et al. from the Allen Institute, rather than investigating the effects on higher-level population analyses, focus their attention on the effects seen in single neurons. When comparing a calcium imaging and an ephys dataset recorded in the same brain region (visual cortex V1) in mice that do the same standardized tasks, what differences can be seen in the firing properties of single neurons?

Due to the high standardization requirements at the Allen Institute, their datasets are probably uniquely qualified to be the basis for such a comparison. Interestingly, they find a couple of clear differences. For example, extracellular ephys data suggest typical firing rates of around 3 Hz (see Figure 2A), which is almost an order of magnitude higher than what has been recorded and estimated from calcium imaging data (see also Figure 7 in our preprint, which estimates spike rates of the same dataset).

The authors go to great lengths to use forward transformations (similar to Wei et al.) in order to reconcile differences seen for various response metrics (responsiveness, tuning preference, selectivity). However, their conclusion seemed to me quite a bit less optimistic compared to the Wei et al. paper. The authors go into more detail when discussing the potential reasons for the discrepancies, and focus on the recording methods themselves rather than on methods to transform between them. In particular their analysis of inter-spike-interval (ISI) violations in ephys recordings (which indicate that spikes from different neurons contaminate the recording of the neuron of interest) was, in my opinion, particularly interesting and convincing. I also really recommend to anybody the last paragraph of their discussion, from which I only want to cite their note of caution about extracellular ephys recordings:

From this study, we have learned that extracellular electrophysiology overestimates the fraction of neurons that elevate their activity in response to visual stimuli, in a manner that is consistent with the effects of selection bias and contamination. – Siegle, Ledochowitsch et al.

One of the reasons why I am writing about these two studies is that I have been working at the interface of calcium imaging and ephys myself, addressing the question, How much information about spike rates can we get out of calcium imaging data? Wei et al. and Siegle, Ledochowitsch et al. take a slightly broader perspective. And, in some way, they show how hard it is to reconcile two (methodological) perspectives on the same phenomenon. (I noticed this in my PhD lab as well, when it came to reconciling results from EM reconstructions of neuronal anatomy and calcium measurements of the same neurons.) Since almost any method in systems neuroscience is technically challenging, we often have in a single lab only a single perspective of a phenomenon, and I think it’s important to always be aware that the conclusions drawn from this perspective might be strongly biased.

In general, calcium imaging and extracellular ephys are extremely valuable tools to observe the living brain, and we better do anything to understand the properties and limitations of these tools. These studies sometimes might feel a bit like negative results and therefore not very attractive to those who want to advance neuroscience, and I therefore understand why not many are willing to undertake these projects. So I am glad to be able to highlight these two publications here.

Posted in Calcium Imaging, Data analysis, electrophysiology, Imaging, Network analysis, Neuronal activity, Reviews | Tagged , | 1 Comment

Alignment tools

This blog posts covers some tools and techniques that I’m typically using to align two-photon microscopes. If you’re an expert, you will probably find nothing new, but if you haven’t been doing this for years, this might offer you some pieces of inspiration.

Aligning a microscope is the process of optimizing the parts to produce better images than before. A useful overview of basic alignment procedures for two-photon microscopes has been put together in this #LabHacks blog post by Scientifica. It also includes safety advice (which I will not repeat here; keep in mind that lasers, especially pulsed IR-laser, are really dangerous, and all safety instructions of your institute and lab must always be obeyed!)

Microscope alignment still seems to be a secret art that most microscope users are afraid of. It is not by principle difficult to learn, but it is a practical skill and requires both patience and a mentor who shows where to buy stuff and how to touch the mechanical and optical elements. Most people working with microscopes are prevented from learning because they are afraid to touch anything. This blog post is intended to help lower the fear of optical alignment – by showing that the tools used for alignment can be very simple.

1. Adjust the beam to a given height above the optical breadboard

After the beam comes out of the laser, you often want to keep the beam in one single plane parallel to the table, in order to keep things simple. In other words, the distance of the beam from the optical table should be the same everywhere. One way to achieve this is to use mounted pinholes (e.g., these ones from Thorlabs). However, it is sometimes difficult to properly see where the beams hits the pinhole, which results in imprecise alignment and unnecessary uncertainty. When I worked with Robert Prevedel in 2013/2014, he showed me a simple trick which makes it very easy to adjust all beam positions to the same height. He used a small hex key and two washers to clamp it horizontally onto an inverted post with a screw, as shown below. The surfaces of the hex key lead to a very nice horizontal alignment, and the precise height indicated by the hex key can be used to (1) adjust the beam itself or (2) to consistently adjust the height of a set of pin holes. It is useful to have and very simple to make.

2. Printed-out resolution targets

Is the beam centered in a given aperture, for example the cage system of the microscope’s tube lens? Normally, I would use a threaded or cage-mounted iris (e.g., this one), but in other cases spatial constraints do not allow this, or the beam can only be viewed from an angle, and it can be difficult to judge whether the beam scattering from the iris is centered or not.

However, if the laser power is low enough not to burn paper, I simply paint resolution targets with Powerpoint or Inkscape and print them out on a sheet of paper resolution targets. When I’m in a hurry, I simply draw it with a pen.

For example, if I want to check whether a beam is collimated (that means, the beam does not change its diameter a lot over the distance), I use these alignment targets as a reference and as a guide for the eye.

Or I use some scotch to fix it into a given aperture, allowing me to check whether the beam is centered or not. Here illustrated for the aperture which is approximately located at the back focal plane of the objective. Does not require much work but is quite helpful.


3. Using mirrors to look around corners

One more practical problem in the above case is the viewing angle. Ideally, I would like to look at the alignment target from the top, but this would at the same time block the beam. To solve this (and many very similar problems), you can simply use a mirror shard. The photo below (left) shows a hand-held piece of mirror which allows to look at the alignment target in a relaxed way. It is difficult to see this from a single picture, but mirrors like this one (either hand-held or fixed in the setup) often make life much easier.

For more convenience, dental mirrors like this mirror (photo below, right) are designed to conveniently look around corners and are of great use to look at pinholes from angles that are difficult to get at without mirrors.

However, be very careful with hand-held and any other moving mirrors! By chance, you might reflect the laser beam into your own eyes! Always be careful and think three times whether there is any chance that any reflection might hit your eyes.

4. Retro-reflecting mirrors and lenses

Usually, the laser beam should hit a lens at its center and at a 90° angle. To ensure that the beam is centered, one can use a) a cage system with a pinhole, or b) a threaded pinhole that can be screwed onto the lens directly, or c) a printed-out removable resolution target (see above). To make sure that, in addition, the laser beam hits the lens with a 90° angle, it is helpful to use back-reflections of the beam. Since a small fraction of the beam is back-reflected by the lens surface, it should ideally coincide with the incoming beam. This can be checked with an IR detection card or, for a visible beam, a piece of paper held close to (but not hiding the) incoming beam path.

Sometimes it is necessary to align a caged element or something similar without lenses that provide back-reflecting surfaces. In this case, you can simply use a mirror that reflects the entire beam back. This mirror can be screwed into a thread of a cage system. More often, it suffices to put a small and flat mirror shard at the flat back of the caged system. This does not provide the highest precision alignment, but is usually good enough for most purposes.

Be aware that back-reflections that go directly back into the laser can make the laser unstable. If a pulsed laser stops pulsing, first check whether any back-reflections could be the reason for it.

5. Retro-reflecting gratings

The main problem with back-reflections from lenses and mirrors is that the beam is often small and coincides with the incoming beam, making it difficult to properly identify.

In my PhD lab, I inherited a really cool tool that was used for alignment of a Sutter MOM scope and which I was, unfortunately, unable to find afterwards on the internet. It is basically a mirror, but with a sort of grating scratched into the mirror surface. Due to interference, the back-reflections were not simply a single beam, but a symmetric diffraction pattern that extended over several centimeters and could be conveniently used for alignment – much more useful and easier to use than the back-reflections of a lens or mirror. I guess that any flat reflective gratings (maybe even compact discs? I have to tried that) could be used for this same purpose.

6. Wedge plates

At several points of the beam path (before entering the back of the objective; after exiting a beam expander) the laser beam should ideally be collimated. The standard method that I used for years was to print out a resolution target made of paper (described above). I used it to check whether the beam changes its diameter when propagating freely over several meters. To this end, it is often necessary to deflect the beam with a mirror that must be temporarily inserted into the beam path.

Fabian Voigt from my current lab in Zürich also showed me the more professional way to check for collimation, using a wedge plate. Wedge plates can be used for shearing interferometers (Thorlabs product, EO product). They generate an interference pattern which can be used to very precisely check the collimation of the beam.

Fabian also kindly pointed me to a paper (Tsai et al., 2015) which mentions that the pattern seen from a shearing interferometer can also be used to analyze less obvious properties of the beam, like coma and astigmatism aberrations (Okuda et al., 2000). It would be cool to have a look-up table of typical interferograms and the corresponding wavefront shapes and aberrations!

7. Alignment lasers

Another tool that Fabian showed to me was an alignment laser.

This alignment laser is basically a continuous wave-visible light laser that goes through a optical fiber and afterwards enters a beam coupler. I used a shearing interferometer (described above) to make sure that the outgoing beam was collimated, and then used the collimated beam for backward alignment.

Backward alignment means to insert the collimated beam of the alignment laser at the location of the microscope’s back focal plane and then aligning the microscope’s components in a backward manner (tube lens -> scan lens -> scanners -> etc.; instead of forward alignment, which starts at the pulsed laser and goes forward until it ends up at the objective). This is helpful for example when two or more separate incoming beams are combined. A second advantage of the alignment laser is that the laser beam is, unlike the near-IR pulsed laser, visible to the human eye.

8. Continuous wave (cw) mode for a pulsed Ti:Sa laser

Standard two-photon microscopes are based on pulsed lasers that have a center wavelength adjusted between 800 and 1000 nm. The light is therefore invisible to the human eye (except for some faint spectral components when the center wavelength is adjusted to 800 nm) and comes with high average power (>1 W, or often much more) and in ultra-short high-energy pulses, making it a very dangerous thing to deal with. IR viewer cards and IR viewers, or less expensive solutions based on simple cameras/webcams, can make the beam visible to our human eye, but I always found this very exhausting to work with. And being exhausted is not a good pre-condition for optical alignment, which requires the ability to stay focused and a sharp mind. Therefore, it would be great if it was possible to make the laser less dangerous: by preventing it from pulsing, by lowering the overall power, and by making it visible…

Fortunately, all of this is possible for the tunable wavelength Ti:Sa lasers (although not for the fixed 1040 nm lasers). For the Spectra Physics MaiTai lasers, the slightly old-fashioned user-interface (screenshot below) allows to manually lower the pump power (also called the “green power”). Change the wavelength to something visible (something between 700 and 750 nm), lower the green power until the pulsing (indicated by the green light in the main control window) stops and continue lowering the green power until the output IR power reaches something like 30-40 mW. This is still much more than the average laser pointer and still dangerous, but less so than a pulsed invisible beam.

It is however important to keep in mind that the control panel is not always 100% reliable. It switches between pulsing and non-pulsing as if these were binary states. But there is a sometimes continuous transition between the two, and sometimes the laser also becomes a bit unstable and oscillates between states – so, better watch the laser for a while before you start working with the beam.

I once talked to a laser engineer from Spectra Physics who told me that it would not be ideal to leave the laser in this non-pulsing state for too long (>> 1 hours). Probably this has something to do with inefficient energy conversion and the heat generated in this process. But I have to admit that I follow this recommendation simply because I don’t know it better. Lasers are complex things, and I mostly treat them in a pragmatic way, like a black box.

9. Improvised test samples

One of the most common mistakes that I keep seeing is to set up a microscope and then try to test it with a real sample, like a living brain of a transgenic mouse. That’s really not the best way to proceed. Instead, always test your microscope with a simple, bright and dead sample. (And, ideally in addtion also with sub-diffraction beads.)

For example, you can use colored plastic items (for example plastic slides), perfect to test how homogeneous the FOV is. Or pollen grains from any flowers outdoors during the pollen season. Or dead, small insects that you caught from your desk. Their surface will often be super bright under a two-photon microscope, especially if you remove the emission filters. Or simply a single strand of your hair (especially if you have still some pigments in them).

In general, most solid things are somehow autofluorescent, and it can be really fun to look at random things below the two-photon microscope. Pay attention not to kill the PMT detector (lower the initial PMT gain), since autofluorescence of natural things can be very bright.

10. More resources on (two-photon or other) optics alignment

Posted in Imaging, Microscopy | Tagged , , , , | 4 Comments

Matlab magic spells

Most neuroscientists who analyze their data themselves use either Matlab or Python or both – the use of R is much less common than in other fields of biology. I’ve been working with Matlab on a daily basis for >10 years, and I have started using Python regularly, although less frequently, about 5 years ago. I’d like to encourage everybody to learn Python, because the language is much more flexible, more beautiful and more powerful. But still I use Matlab often to get a first impression of data, because it is simply perfect for browsing intermediately sized (MB to GB) datasets and for doing simple things in general, also due to the plotting interfaces which I find are easier to use for data exploration than matplotlib/seaborn in Python.

Here, I would like to share a collection of useful Matlab commands which are not complicated or fancy but still nice to know. Don’t expect too much, and don’t read this if you have spent a couple of years with Matlab already. To be entirely honest, I’m writing this also for my own reference, because I keep forgetting the exact wording of some of the commands, and it’s nice to have everything in a single place. Without further ado, here is my list of favorite Matlab spells:

1. Leading zeros with sprintf()
2. Automatization with dir()
3. Exclude file patterns based on regular expressions
4. Callback functions to speed up image inspection
5. Regionprops for image segmentation
6. The curve fitting toolbox
7. Manually modifying colormap
8. Turning the background of a figure white
9. Make ticks point outwards instead of inwards
10. Rotate x-tick labels
11. Reverse y-axis direction
12. Save large figures to vectorized format (eps)
13. Get raw data from figures
14. Position the figure window at a given screen location
15. Profiling code with tic/toc

Continue reading

Posted in Data analysis | Tagged , , | Leave a comment

Annual report of my intuition about the brain (2019)

How does the brain work and how can we understand it? I want to make it a habit to report some of the thoughts about the brain that marked me most during the past twelve month at the end of each year – with the hope to advance and structure the progress in the part of my understanding of the brain that is not immediately reflected in journal publications. Enjoy the read! And check out previous year-end write-ups: 2018, 2019, 2020.

The purpose of neuronal spikes is not to represent …

It is a common goal of many neuroscientists to “decode” the activity patterns of the brain. But can the brain be simply understood by finding observables like visual objects or hand movements that are correlated with these activity patterns? Nobody would deny that neuronal activity is correlated with sensory inputs, motor outputs or internal states. But what does it help to look at these correlations? In other words, are these representations of sensory, motor or inner world key to understanding the brain? In the following, I will argue against this “representational” view of the brain.

First, what do I mean with “representational”? In the following, a neuron represents e.g. an external stimulus if its activity pattern is correlated with the occurrence of this external stimulus (for a stricter definition, see below). For example, a neuron in visual cortex fires when a drifting grating of a 55° orientation is presented to the animal. One could say that “the neuron codes for 55° drifting gratings” or “the activity of the neuron represents 55° drifting gratings”.

To give another example, a neuron in hippocampal CA1 fires reliably when the animal is at a certain location (“the neuron codes for this and that place” or “the neuron represents this and that place”). This representational view is probably the most intuitive approach when it comes to characterizing the activity of the brain. However, simply representing the internal or external world is not enough. What else does a brain do?

Before trying to answer this question, let’s speculate why the “representational” view is so deeply engrained in systems neuroscience. First, many thinkers and also neuroscientists have been biased towards a view of the brain as a passive receiving device because they have been sitting passively at a desk in front of an empty sheet of paper while shaping their thoughts. Second, the representational view is also very intuitive: finding something that can be observed, like a visual stimulus, to be reflected directly in neuronal activity seems to explain why this neuronal activity happened. Animals or humans subject to brain experiments are often told to hold still while waiting for a stimulus or cue. It is clear from this typical experimental design that it would be more likely to see “receptive fields”, which is a classic finding of representations in the brain, because the experimental design biased the experience of the subject towards passive receiving. Torsten Wiesel and David Hubel, who have pursued this line of thinking in the most compelling way in the early era of neuroscience, are probably the best examples for this approach:

Apart from that, artificial neuronal networks like the attractor networks put forward by John Hopfield and, more recently, deep convolutional networks, are typically based on classification tasks, e.g. to identify written numbers on an envelope, or to correctly categorize the object seen in an image.

For a device that is supposed to perform such a task, it is clearly sufficient to represent the inputs and the possible outputs. But this is not the main task of a brain.

… but to act upon something

The brain, unlike a pattern classifier, is in a closed loop with its environment, being moved by it and acting upon it at the same time. It is not a passive observer, but an involved device, an organism. The brain is not designed to represent its environment. Instead, it has evolved to act upon its environment.

How could a model of the brain account for this fact? An approach to circumvent the limitation of the “representational” view of the brain is to focus on its generative abilities. The brain is clearly capable of forming mental models of the external world, of other minds, of physical laws, and also of abstract concepts. “Predictive processing” is a term covering a variety of models also applied to the level of single neurons. In these models, predictions and prediction errors are not purely representational, but process information as a bi-directional interaction between sensory (bottom-up) and generative (top-down) content. However, these models consists of hierarchical networks that, again, perform classification tasks or other well-known toy problems from the machine learning community. Therefore, these models often appear to be representational models in disguise, rather than something completely different.

A researcher that inspired me to think into a slightly different direction is Romain Brette, who published an opinion article in 2018 called “Is coding a relevant metaphor for the brain?” (In this article, Brette also defines more precisely the meaning of “representations” and its relationship to “coding” and “causality”; but since I do not think there is a semantic agreement among neuroscientists, I decided not to use this more precise terminology here.) Brette’s article is first and foremost an article against the idea of “coding” and “representations” as useful concepts for understanding the brain. Naturally, it is a bit less clear what alternative he is actually suggesting and how this alternative could be implemented or understood. (In a comment to the opinion piece by Brette, Baltieri and Buckley convincingly argue that a set of predictive processing models that use active inference is an interesting alternative to the representative view and is close to the “subjective physics suggested by Brette.) Although some concepts in the opinion article remain a bit vague, Brette also touches upon the idea of understanding neurons not by analyzing their representational properties, but by analyzing their actions, i.e., their spikes:

“[C]oding variables are observables tied to the temporality of experiments, whereas spikes are timed actions that mediate coupling in a distributed dynamical system.” (Romain Brette)

From this paper and from posts on his blog, it seems to me that he is trying to circumvent the natural and historical representational bias of systems neuroscience by investigating systems which are simple and peripheric enough to be amenable to descriptions that are not based on representations, but actions/reactions. One area of his research, which I find particularly interesting, is the study of single-cellular organisms.

The basic unit of life is the single cell

During my PhD, I gradually went from large-scale calcium imaging to single-cell whole-cell recordings. I was fascinated by the complexity of active conductances and the variability between neurons, and as I mentioned in last year’s write-up, I also got interested in complex signal processing in a single neuron.

Therefore, investigating unicellular organisms in order to understand neurons and ultimately the brain was not so far from my intuition anyway. And although neurons do not move their processes a lot and seem to be rather immobile and static at first glance, there are many examples of single-cell organisms that exhibit quite intriguing behaviors and remind us of what single cells are capable of. Just watch this movie of this tear drop-shaped ciliate that uses a single long process to search for food:

And it’s just a single cell! Interestingly, these ciliates not only have an elaborate behavioral repertoire, but also use spikes to trigger these behaviors. (Check out this short review on “Cell Learning” to get an idea about what single-cell organisms are able to learn.)

At least in the beginning of life, the single cell was the basic unit of life, receiving information about the external world, using its own memory and acting upon the world. Of course, many cell types in humans are far more specialized and cannot really be considered a basic unit of life, since they serve a very specific task like forming a part of a muscle or sitting in the germline. However, information-processing cells like cells in the immune system or neurons in the brain, which live at the interface between inputs and outputs of an organism, are much more likely to be similar to these unicellular forms of life. They make most of their decisions on their own, based on local signals, without asking a higher entity what to do.

If we look at life with human eyes, the mind and the human body tell us to search for an understanding of neurons from the perspective of the whole brain (or the whole body). However, if we look at life from a greater distance, life shows its highest richness at the level of single cells, and it would therefore make sense to search for an understanding of neurons from the perspective of a single neuron.

The objective function of a single cell

What does it mean to understand a single cell? A very biologist way would be to make an as complete as possible list of all its behaviors and states when exposed to certain conditions. For example, a hypothetical cell senses a higher concentration of nutrients at one of its distal processes and will therefore move into this direction. Conversely, if it senses a lower concentration in this direction, it will move into the opposite direction. If the nutrient level is overall high, it will not move at all. If it is overall low, it will try to move away in a random direction. This is a biologist’s list of behaviors.

A more compressed and therefore better understanding would be in my opinion to write down the “goals” of this single cell. Or, to abuse an expression by mathematicians and machine learning people, to write down the objective function (or loss function) of the single cell. For the above-mentioned hypothetical cell, this would be simply: move to places with high nutrient concentrations.

If we shift back from ciliates in a pond to neurons in the brain, this leads to a very obvious, but also surprisingly difficult-to-answer question: what is the objective function of a single neuron?

Of course it is not clear whether it is reasonable to assume that the objective function or optimization goal is set on the single-cell level. But given that cognition is most likely a self-organized mess that is only coarsely orchestrated by reward and other feedback mechanisms, I would tend to say: no, it’s not unreasonable. (But check out below for more detailed criticism.)  Therefore, the question is, what are the goals of a single neuron? What does it want?

Let’s assume that this question is a good one. How can it be answered? I see three possible approaches:

  1. To investigate the evolution of self-sustaining cells into neurons of the central nervous system. Since objective functions of uni-cellular organisms seem to be rather accessible and objective functions of neurons seem to be difficult to guess, one could investigate intermediate stages and try to understand the evolutionary development.
  2. To understand the development of the brain and the rules that govern migration and initial wiring of a single neuron. This would possibly allow to make a list of development behaviors for a certain neuronal type and from that extract an objective function. This would be an objective function not of regular operation, but of development. That’s exactly what developmental neurobiology is actually doing.
  3. To observe the full local environment of a neuron (for example, its electrical inputs) and to try to extract rules that govern changes in its behavior (for example, synaptic strength or output spikes)

The third approach is probably the most difficult one because it is currently technically impossible to observe a significant fraction of the inputs of a neuron together with its internal dynamics. There have been many studies on the plasticity of very few synapses in artificial conditions in slices using paired recordings, but it is currently impossible to perform these paired recordings in a behaving animal. Even in slices, these experiments are daunting given the myriad of different neuronal cell types and the expected variability in these experiments. Effects like long-term potentiation, long-term depression, short-term facilitation, short-term depression, spike-timing dependent plasticity, immediate early genes, receptor trafficking and neuromodulatory effects are just some of many processes that are essential to understand what is happening to a single cell.

A different approach that, however, still resonates with the idea of single cells as actors has been put forward by Lansdell and Kording in 2018. Their idea, as I understand it, is that a neuron could try to estimate its own causal effect, that is, the effect of its actions. This is possible because its actions are a discontinuous function of their state, the membrane potential. If external conditions are almost identical in two situations but the neuron fires a spike only in one of the two situations, the neuron could extract from the received input the effect of this single spike. Therefore, the neuron would measure the causal impact of its actions on itself. This idea is very close to that of a unicellular organism that immediately feels the effects of its own actions.

But what could be the objective function of such a neuron? – For example, the goal of a neuron that receives recurrent feedback could be to find a regime of a certain level of recurrent excitation. Recurrence could be measured by estimating the effect of the spikes of the neuron on the received feedback, possibly also in a dendrite-specific manner. I could imagine that objective functions almost simple as that could make a neuronal net work, once it is embedded in a closed-loop environment with sensory input and motor output. Also, I think that this way of thinking could be close to ideas put forward by e.g. Sophie Denève about predictive processing in spike-based balanced recurrent networks (for example, check out this proposal), which heavily relies on the self-organizing properties of adaptive recurrent circuits with a couple of plasticity rules.

An additional strength of the idea to use spikes to measure a neuron’s causal impact is that it could explain why neurons fire stochastically: this way, they acquire information about their causal impact that is hidden in a regime of deterministic firing.


Clearly, these ideas need more refinement and restructuring, and it is obvious that some of the feasible experiments (spike-timing dependent plasticity, behavioral time scale plasticity) have been done already. But I still like the idea of reframing these experiments by analyzing a potential objective function of a single neuron.

Beyong that, it would be interesting to try to understand self-organized systems of artificial neurons that are not governed by global loss functions, but by neuron-intrinsic and possibly evolving objective functions. I’m really looking forward to seeing theoretical work and simulations that do not focus on the behavior and objectives of the whole network, but on the behavior and objectives of single units, of single cells.


Unicellular organisms and neurons are both single cells. Left: Unicellular ciliate lacrymaria olor, adapted from Wikicommons (CC 2.0). Right: Two-photon image of hippocampal pyramidal neurons in mice.

Appendix: the devil’s advocate

Criticism 1: Representations and coding seem to be still the best way to think about neuronal activity

Although I find the ideas above appealing conceptually, it is clear from past research that most of our understanding of the brain comes from correlating the observed brain activity with sensory or motor information, because it allows to conclude when certain brain areas or neuronal types are active and what happens if they fail. This functional anatomy of the brain might not reveal all the design principles of the brain, but it could be necessary to make an educated guess about those principles.

In addition, our human languages are based on concepts (“you”, “house”, “hope”, “cold”) rather than on dynamical processes. For example, the science of language, linguistics, has coined the term ‘signified’ which describes the content of such a concept as a central element of language. From this converging view from both naive usage and scientific study of languages, it seems to me likely that high- and low-level representations indeed exist and are accessible in the human mind – probably as an evolutionary byproduct of mental processes that were more directly connected to actions or sensation. – It is still interesting that some sensory modalities are more vaguely represented in thought and language (olfaction) than others (vision).

Criticism 2: It is so far technically impossible to record the perspective of a single neuron

It seems to be relatively easy, although not trivial, to monitor a unicellular ciliate and its environment. However, the environment of a neuron consists of all of its synaptic partners, which can be in the range of thousands, and any other source of electrical or chemical neuromodulation. A single neuron itself is widespread over such a large volume that is very far from being amenable to any electrophysiological or imaging technique that provides sufficient resolution in time and space.

The best candidate technique, voltage imaging of both soma and dendrites, cannot be used over periods longer than minutes due to bleaching, suffers from low signal-to-noise as well as difficult calibration and requires imaging speeds and photon yields that are impossible to achieve with any imaging technique and fluorescence sensor that I know of.

But researchers have proved before that imperfect methods can be used successfully. I’m curious about all technical developments and what they will bring.

Criticism 3: Regarding neurons as actors is a case of anthropomorphism and ignores the interactions in complex systems that can lead to emergence

If I scrutinize my own ideas, I realize that one reason why I like the idea of neurons as actors is that this idea is simple and connects to the human desire to have an intuitive, empathic understanding of things. For example, we tend to personify even the smallest animals like ants, bees or ciliates in movies, children books or youtube videos. We tend to personify items of our personal households like a cup, a table, or an old house; or our plants in the garden; or the ocean, the wind, a dead tree. On the other hand, I feel that we are largely unable to feel the same empathy with distributed systems like the mycelium of a mushroom that spreads over kilometers in a forest, or the distributed power supply grid, or a large fish swarm, or a deep artificial neuronal network. We can feel amazement, but no empathy. I can imagine to be single water molecule that vibrates, rotates and is tossed around by the Brownian motion of the environment and by the convection of the ocean, but I cannot imagine to be the emergent phenomenon of water waves breaking at the shoreline which are the complex environment of this single water particle.

I think that there is a tendency in human thinking to avoid complexity and to impose a simple human-like personality onto all actors. Therefore, we should observe ourselves and try to avoid our tendency to personify parts of complex systems when they actually cannot be disentangled from their environment. It is easy to tell a story about specific neurons as important actors and what they want to do and how they feel immediate feedback of their actions like we do when seizing an apple, but it is difficult to tell a story about complexity, which is over the top of our heads. But maybe this, which we can hardly understand, is the truth, and this is good to keep in mind.

Posted in Calcium Imaging, Data analysis, electrophysiology, machine learning, Network analysis, Neuronal activity, Review | Tagged , , , , , , , | 8 Comments

A practical guide for adaptive optics

There is no standard curriculum to learn practical procedures about microscopy: how to align a setup, how to identify misalignments, how to identify broken parts, where to buy components, how to check their performance, and much more. How to learn all of this?

Often, research papers contain close to zero information about these topics, and even if they do, the information is scattered in the methods sections or supplementary information. This is not a surprise since papers are mostly expected to be “research” and not engineering manuals or step-by-step protocols. This makes some papers quite impressive, but close to useless from a practical perspective.

Therefore, the most useful practical information about microscope construction and maintenance can be found in sources other than papers. Some examples:

  • On labrigger and other blogs (including my blog), useful but sometimes a bit anecdotal information is distributed over many pages, making it available via search machines.
  • A nice tutorial on alignment by Rainer Heintzmann can be found in the web – it’s an appendix to a book, but clearly more widely shared than the book itself.
  • Some projects, like the mesoSPIM project (by my labmate Fabian Voigt), provide a lot of useful practical information via dedicated websites or Github repositories.
  • Andrew York took advantage of his independence to provide less formal work-in-progress publications via Github, which allowed to focus more on practical and useful information, instead of trying to prove novelty and impact with beautiful figures.

These webpages are more flexible than PDF papers, can be easily updated and enhanced and allow for simple integration of animations or interactive elements. Since the webpage does not need to convince editors and reviewers with polished writing and impressive claims, it allows for more practical information to slip in.

Recently, I discovered another project that goes into a similar direction, initiated by the lab of Martin Booth, on the topic of adaptive optics. As the authors describe it themselves:

“The documents posted here will include tutorials, experimental protocols, and software. This will range from simple hints and tips through to extensive documentation of procedures. We intend to post material whenever it is ready, and documents will be updated with newer versions when we have them. We chose to set up the site in this way, to be free from the constraints of the traditional publishing process, which is ill-suited to the dissemination of this type of material, particularly when content could be frequently updated as our own approaches develop.” –

I really enjoyed the protocols uploaded so far. For example, despite quite some experience with microscopes, I have never aligned a confocal pinhole so far, and I found it interesting to read the instructions on how to do this, even though it’s a very simple procedure. I’m really looking forward to seeing this form of resource becoming more frequently used and also appreciated not only by other researchers, but also by funding sources.

Of course, all of this does not replace the best source for practical training in microscopy: doing it yourself, spending a lot of time solving problems – and having somebody in the lab who knows everything.

Posted in Imaging, Microscopy | Tagged , , | 1 Comment

Review: An artificial ground truth for calcium imaging

Selected paper: Charles, Song, Tank et al., Neural Anatomy and Optical Microscopy Simulation (NAOMi) for evaluating calcium imaging methods, bioRxiv (2019).

What is the paper about? Calcium imaging is a central method to observe neuronal activity in the brain of animal models. Many labs use rather complicated algorithms to extract meaningful information from these imaging data, but the results of those algorithms are often hard to judge: is this an artifact of the algorithm or the imaging method, or is it a biologically meaningful signal of a single neuron? The paper by Charles et al. addresses this question by simulating both neuronal activity in anatomically realistic neurons and the imaging process that makes these activity patterns visible. By knowing both the true simulated biological signals and the simulated imaging data, the procedure can benchmark analysis algorithms against a known ground truth. In addition, the method can be used to evaluate imaging modalities and their impact on resolution, movement artifacts and other factors for a specific use case.

More details: Neurons are simulated as realistic 3D structures with dendrites and axons, based on 3D EM and light microscopy data of neuronal tissue. Their activity is simulated as spikes. The spikes are transformed into slower calcium events, which in turn are transformed by the binding kinetics of the calcium indicator. Binding to the calcium indicator gives rise to a change in fluorescence, based on the Hill model of cooperative binding. Then, things are brought together by simulating the excitation laser beam and the light emitted by the fluorescent calcium sensors, assuming a typical signal to noise level and shot noise induced by the limited number of photons seen by the microscope. Altogether, this is a quite impressive simulation pipeline, covering anatomy, calcium sensors and binding kinetics as well as light simulation and relevant noise induced by instruments and other variables (for example, residual artifacts induced by the motion of the animal).

Evaluating demixing algorithms: In my opinion, the most interesting part of the paper is section 2.3 (“Evaluation of automated segmentation”). Here, the authors use their simulations to benchmark three commonly used algorithms for automated source extraction. These algorithms are used to automatically extract regions of interest (ROIs) from an imaging dataset, and are often used as a repeatable and hopefully more objective and more reliable replacement of the manual drawing of ROIs by a human being. Unfortunately, the authors do not benchmark those methods also against a human expert who manually selects ROIs.
Anyway, with respect to the three state-of-the-art algorithms, the evaluation using the calcium imaging simulations reveals a rather limited precision of all investigated methods. First, the different algorithms find rather different sets of neuronal activity units with relatively little overlap between the units found by each algorithm. In addition, also the absolute number of units found by each algorithm is quite variable (for a specific simulation of L2/3 imaging in visual cortex, the number of found components varies between 265 and 1091 for the different algorithms). More details can be found in the relevant section in the paper.
These are interesting findings. They do not disqualify demixing algorithms, but they should lower our trust that the large majority of the extracted components is based on activity signals of single neurons.

The devil’s advocate: The strength of the simulation in this paper (its realism) is also its weakness. The authors themselves state:

“Simulation-based approaches, however, often suffer from being either too simple or too complex.”

How can we be sure that they have found exactly the right level of detail? We can’t, or not very easily. There is an endless list of possible details which could have been omitted: For example, the calcium indicator concentration was set to 10 μM for the simulations, but these values might differ between neuronal types and depend on other factors as well. To give another example, the simulation assumes that the calcium concentration is constant across the dendritic tree, therefore omitting the known effect of localized calcium dynamics in dendritic spines.
And there are many more approximations (optics, physiology, anatomy) made by the authors that are difficult to judge. For example, I would have been glad to see a simulated excitation PSF at a certain tissue depth alongside with a measured PSF at the same depth (e.g., using beads injected into cortex), to make sure that the scattering simulation, which is too complex to be judged without investing a lot of time, is realistic or not. – However, these are details, and for most of the analyses performed by the paper (like the discussion of demixing algorithms), this level of detail of the simulation is probably not important – and I’m looking really forward to having this simulation tool publicly available. If it is user-friendly and easy to adapt, it could become a standard tool to check in silico what to expect for in vivo experiments.

Conclusion: Very interesting and possibly useful work, although it is difficult to understand the limitations and details of the simulation.

Further reading: A paper from he same lab on a related topic: Gauthier, J. L., Tank, D. W., Pillow, J. W. & Charles, A. S. Detecting and Correcting False Transients in Calcium Time-trace Inference.

Other paper reviews on this blog: on precise balance in the hippocampus [1]; on L5 apical dendrites [2].

Posted in Calcium Imaging, Data analysis, Imaging, Microscopy, Neuronal activity, Reviews | Tagged , , , , | 2 Comments

The cell-attached soundtrack of calcium imaging

Old-school electrophysiologists like to listen to the ephys signals during experiments. For example, this allows to precisely hear when the patch pipette approaches a target neuron. The technique is discussed in the Axon Guide: “Audio Monitor: Friend or Foe?”.

The following is something very similar: Calcium imaging of neuronal activity (GCaMP6f), and an audio track that is derived from a simultaneously performed cell-attached recording of the same neuron (detected action potentials are convolved with a particular sound event, like two metals hitting each other, the sound of an anvil, and the sound of gunshots – the “firing” neuron).

Three times the same calcium recording, but with different soundtracks:

CC-BY 3.0,

CC-BY 3.0,

CC-BY 3.0,

Posted in Calcium Imaging, electrophysiology, Imaging, Neuronal activity, zebrafish | Tagged , , , | Leave a comment

Post-publication review: Somato-dendritic coupling of L5 neurons in V1

It requires more than a quick look at the abstract and the figures to fully understand a research paper and its limitations. One way to get there is to write a summary or critical review of a paper. In a contribution to (informal) post-publication review, I will select neuroscience papers that are, in my opinion, worth the time needed to write up a review. other reviews: On precise balance in the hippocampus [1]; on simulations of calcium imaging datasets [3].

Selected papers:

1. Beaulieu-Laroche, Harnett et al., Widespread and Highly Correlated Somato-dendriticActivity in Cortical Layer 5 Neurons, Neuron (2019)

2. Francioni and Rochefort, High somato-dendritic coupling ofV1layer 5 neurons independent of behavioural state and visual stimulation, bioRxiv (2019)

L5 neuron

Simplified scheme of a L5 pyramidal neuron. From Wikipedia, used under CC BY 4.0 license (modified).

What are the papers about? The mammalian cortex is a layered structure that harbors different neuronal types in different layers. One of the most prevalent and most fascinating neuronal types is the layer 5 (L5) pyramidal neuron (schematic drawing above). While the cell body resides deep in layer 5, several hundred micrometers below the surface, it sends a so-called ‘apical trunk’ dendrite up to superficial layers, where it branches into a so-called ‘apical tuft’. The apical tuft receives input which is fundamentally different from input received at the soma. This, together with the long distance between apical tuft and soma, poses a simple, but difficult-to-answer question: how tightly are electrical activity at the soma, the apical trunk and the apical tuft coupled, and how do the electrical compartments interact? The two papers address exactly this question. – Their main finding is that calcium activity at the soma and the apical dendrite are very tightly coupled, no matter what happens to the mouse.

Why this question is important: Six months ago, I wrote a blog post which explained why it would be so interesting to see uncoupled activity of soma and the apical dendrite. Several studies reported – using indirect methods – that there was strong uncoupling happening during sleep or motor learning, but the few studies working on this topic that used direct methods have rather shown tight coupling between apical and somatic activity (Helmchen et al., 1999; Kerlin et al., 2018), with only few outlier events.

These two papers try to address this question again, in a different brain region, in various behavioral conditions, and probably more systematically and in a more targeted way than previous studies.

The method: Both teams performed simultaneous calcium imaging in two optical layers from L5 neurons in the primary visual cortex (V1) of transgenic GCaMP6 mice. The two layers covered either soma and apical trunk, or apical trunk and apical tuft.

Calcium does only partially reflect the electric events that happen in the respective location. Calcium events are an imperfect, non-linear and supra-threshold mirror of electrical activity, and, in addition, there are calcium events independent of electrical activity. Still, calcium imaging is the best available method to investigate somato-dendritic coupling, since theoretically better-suited methods (dendritic patching, voltage imaging) do not work reliably in awake animals for now.

I also would like to highlight the section “Limitations of the use of calcium imaging to assess dendritic activity in awake behaving mice” from the Francioni and Rochefort paper, which discusses in a bit more detail the challenges for dendritic calcium imaging. I think it is a very good idea to include this part in the main paper and not in the methods section, because this is an experiment that requires careful analysis due to motion artifacts, unwanted signal from spines or axonal boutons, imperfect calcium indicators and other problems.

Paper 1 (Harnett lab): The main result of the Beaulieu-Laroche et al. paper is that somatic activity and activity measured in the distal apical trunk (distance between dendritic and somatic regions: ca. 350 μm) are strongly correlated in V1. One would expect that this tight coupling becomes loose when the dendritic tree is bombarded by more input. Surprisingly, however, they find that this tight coupling is independent of visual stimulation or whether the mouse is running or sitting still.

I haven’t talked to the authors, but from this exploratory study design, which tested different conditions (visual input vs. no black screen, running vs. pausing), I would guess that they were hoping to see a behavioral regime where the tight somato-dendritic coupling breaks down, due to strong excitatory apical inputs, due to a modulation of the leakiness of the apical trunk, due to shunting inhibition or another weird phenomenon on the molecular level. The fact that nothing of this happens seems a bit disappointing from the perspective of the study director.

To better understand these calcium imaging measurements, the authors performed calibration experiments in vitro (slices), where they measured both calcium signals in somata and dendrites and patched the respective locations with micropipettes to record the electrical signals. Somehow, the whole-cell recordings did not wash out the GCaMP from the neurons (Fig. S3), which was a bit surprising to me. Using these slice experiments, they conclude that the high somato-dendritic correlations are either due to burst-firing somata that trigger (active) dendritic electrogenesis, or due to dendritic spikes that elicit somatic firing. The authors are unfortunately unable to discern these two possibilities.

To make things a bit more complex, a side-figure shows that single L5 spikes are not seen using the GCaMP6f reporter (but are seen using an OGB-1 calcium dye, see Fig. S3G,H). The authors find that single somatic L5 spikes also do not trigger dendritic electrogenesis, which seems to explain that the calcium signal is correlated also in this regime of single spikes. Indirectly, this suggests that there is no tight coupling of calcium activity between soma and apical dendrite if there is only one somatic spike.

Paper 2 (Rochefort lab): Francioni and Rochefort perform similar calcium imaging experiments. Due to technical reasons (using a piezo instead of a tunable lens), their two imaging planes are closer to each other (<200 μm). As a main result, the authors come to the conclusion that calcium activity is highly correlated not only between soma and (proximal) apical trunk, but also between the proximal and distal trunk and between distal trunk and the apical tuft; and also between processes within a single apical tuft. Including not only the apical trunk but also the apical tuft is an important extension compared to the paper from the Harnett lab. As another difference, they use for some of their experiments GCaMP6s instead of GCaMP6f to have a more sensitive readout of somatic activity.

Very similar to Beaulieu-Laroche et al., Francioni and Rochefort report that this tight somato-dendritic coupling is not affected by visual stimulation or by locomotion, also for the apical tuft. Actually, they report that “less than 3% of the total number of transients” were branch-specific, and that those few events were mostly low-amplitude signals.

The authors then analyze the data more carefully and find that quite some calcium events that can be seen in the soma do not propagate into the apical trunk, and many events in the apical trunk do not propagate into the apical tuft. They quantify the overall loss from somatic activity to the dendritic tuft to ca. 40%. In my opinion, this suggests that most of the calcium activity is generated first in the soma and then propagates into the apical compartment, where it is (according to Beaulieu-Laroche et al.) amplified by active conductances, but not always (as seen by Francioni and Rochefort). One remaining question is why the attenuation from soma to apical trunk has not also been observed by Beaulieu-Laroche et al., but this is probably due to technical reasons, for example the calcium indicator.

Consequences for understanding the role of apical dendrites: During the last 20 years, many people, in particular neuroscientists coming from slice physiology or from theory, have speculated about active roles of dendrites and the computational independence of dendritic segments. Although these two studies (hopefully) are only the beginning of more detailed research covering other brain regions and other behavioral conditions with similar methodological approaches, the findings seem to dampen the enthusiasm about dendritic processing a bit.

Francioni and Rochefort make it very clear, that “this result is consistent with two possible mechanisms: either (1) global apical tuft activation systematically triggers somatic action potentials, or (2) global tuft calcium transients are triggered by back-propagating action potentials alone or in conjunction with tuft synaptic inputs.”

Both options would be possible. In my opinion, however, their paper points rather towards the second option. First, they see that 40% of somatic events do not make it to the apical tuft; these are very likely back-propagating action potentials that get stuck in the middle of their journey in the apical dendrite. Second, given that they do not see the other phenomenon prominently (dendritic events that do not make it to the soma, <3% of events), it seems very likely that a very large fraction of all dendritic events are evoked by back-propagating action potentials.

This suggests a (to-be-confirmed) picture of L5 pyramidal neurons that are driven by basal instead of apical inputs. This is consistent with recent experiments that micro-dissected the dendritic tree, although of L2/3 neurons (Park, Papoutsi et al., bioRxiv, 2019). Somatic spikes would then propagate into the apical tree and, if elicited at high-enough frequency, recruit active dendritic sodium and calcium channels. The activation of the dendritic conductances could in addition depend on local apical input and play several roles, for example: 1) enhancement of the somatic spikes (bursting), 2) triggering of synaptic weight changes in the apical dendrite, 3) triggering of transmitter release to signal from apical locations to their respective pre-synapse. However, it would be unlikely that apical input by itself generates a calcium spike that activates the soma.

Conclusion: The studies by Francioni and Rochefort and Beaulieu-Laroche et al. contribute a lot to a better understanding of the function of apical dendrites of L5 pyramidal neurons. Even though the main result might seem like a partial disappointment at first glance, this kind of systematic study (instead of a study strongly driven by an hypothesis) seems to be the right thing to me. And more of the same, with different behavioral conditions, other cortical regions and possibly a more detailed analysis of the rare dendritic events seems to be the way to go. These two very interesting studies provide important information about somato-dendritic coupling, but the apical dendrite of L5 neurons still largely remains mysterious.

Posted in Calcium Imaging, electrophysiology, Microscopy, Reviews | Tagged , , , | 3 Comments

Post-publication review: “Precise excitation-inhibition balance controls gain and timing in the hippocampus”

It requires more than a quick look at the abstract and the figures to fully understand a research paper and its limitations. One way to go there is to write a summary or critical review of a paper. In a contribution to (informal) post-publication review, I will select neuroscience papers that are, in my opinion, worth the time needed to write up a review. On L5 apical dendrites [2]; on simulations of calcium imaging datasets [3].

Title: Bhatia et al., Precise excitation-inhibition balance controls gain and timing in the hippocampus, eLife (2019).

Target audience: Theoretical neuroscientists working on balanced networks and the role of inhibition; people working on the CA3 to CA1 feedforward circuit in hippocampus; people working on rules governing the connectivity of inhibitory interneurons; people interested in codes dependent on precise spike timing.

Why I chose this paper: I have worked on a similar topic at the interface of theoretical and experimental neuroscience, the so-called ‘precise balance of excitation and inhibition’ (Rupprecht and Friedrich, Neuron (2018)). Of course I’m interested what other experimental studies have found out about this important, although not very easily accessible concept.


From Bhatia et al., eLife (2019), used under CC BY 4.0 license (modified from Fig. 1).

Summary: The authors use a slice preparation to dissect the feedforward circuit from CA3 to CA1 in the mouse hippocampus (figure above; left part).

Two main findings.

The first finding (‘detailed balance’): They stimulate CA3 that expresses channelrhodopsin with an artificial set of patterned light stimuli (figure above; lower left and right parts). Then, they use whole-cell recordings to measure inhibitory and excitatory currents in single neurons in CA1. They find that even for very weak stimuli that presumably elicit few or even single action potentials only in CA3, an inhibitory current is evoked in CA1. Moreover, and even more interesting, the size of this inhibitory current matches the excitatory current not only on average, but also for individual stimuli. This feature has been termed ‘detailed balance’ by theoretical neuroscientists.


From Bhatia et al., eLife (2019), used under CC BY 4.0 license (modified from Fig. 2).

The second main finding: The balanced inputs as recorded with voltage-clamp end up, however, in a non-linear integration as observed with current-clamp recordings. In particular, the authors observe a supra-linear relationship which they describe as ‘subthreshold divisive normalization’. In addition, they notice that inhibition kicks in faster in the CA3 neuron when the excitatory (and inhibitory) inputs to the CA3 neuron are stronger. They use this finding and simulations to support the idea that the reduction of inhibitory delays is the underlying cause of the observed divisive normalization.

In the following, I will only discuss the first finding. The second finding is already broadly discussed in the paper and also by the reviews (which are partly accessible, thanks to eLife!). In general, I find the first finding (detailed balance) much more interesting, whereas the main focus of the paper is on the second finding.

Consequences of a detailed balance: A detailed balance in this feedforward path from CA3 to CA1 is something unexpected, because it would require very specific feedforward connectivity of interneurons matching the feedforward connectivity of excitatory connections. Importantly, the authors put this in the context of previous work on CA1, which found that almost any CA1 cell can be converted into a place cell for any spatial location:

“Precisely balanced networks, with all input subsets balanced, are well suited for input gating (Barron et al., 2017; Hennequin et al., 2017). The finding that most silent CA1 cells can be converted to place cells for arbitrary locations predicts the existence of an input gating mechanism (Lee et al., 2012), but the nature of this mechanism remains unknown. One prediction of precise balance is that inputs for multiple potential place fields may be balanced, and hence place field activity is gated ‘off’. Evoked depolarizations (Lee et al., 2012) or dendritic plateau potentials (Bittner et al., 2015; Bittner et al., 2017), which potentiate the subset of active synapses, that is, change the I/E ratio (Grienberger et al., 2017), can flip the gate ‘on’, thereby converting a silent cell to a place cell for that specific place field”, write Bhatia et al., eLife (2019).

The devil’s advocate: Detailed balanced means, as mentioned before, that the size of an inhibitory current matches the excitatory current not only on average, but for individual stimuli. It is important, as I have described in my own experimental work (Rupprecht et al., Neuron (2018)), that this statement of a detailed balance can only be made if the different individual stimuli are equally strong.

It is therefore not a good idea to assess detailed balance by pooling over all stimuli with 1, 3, 5, 7 and 9 squares of optogenetic activation. Instead, one should compare within a set of stimuli with e.g. 1 square. This is done, as far as I understand it, only in Fig. 2h (example shown for a single neuron) and, based on a current-clamp dataset, in Fig. 2/Supplement 2e-f. It would have been interesting to analyze this aspect in more depth.

More importantly, even when analyzing e.g. only stimuli that all of them correspond to the stimulation of one “square” with photoactivation, this does not mean that all those stimuli are equally strong – which seems to be asssumed in the modeling part (Fig. 2/Supplement 2i,k). As I understand it, the most parsimonious explanation of the data would not be a detailed balance of excitatory and inhibitory inputs in this feedforward circuit, but a distribution of activation strengths for each stimulus set.

For example, stimulation of grid position 3 may overall result in 17 action potentials in CA3, whereas stimulation of grid position 7 may result in 55 action potentials and grid position 9 in zero action potentials. Feedforward inhibition could then simply integrate all action potentials in CA3 and send out a global signal proportional to this number. Indeed, even a single interneuron could, in theory, do this job. Therefore, it seems that the results can be explained by simple and plausible assumptions that do not require detailed balance. Altogether, I think that the paper provides some evidence for a detailed balance in the CA3 to CA1 feedforward circuit, but given the unexpectedness of this result, the evidence is a bit too weak to make me believe the conclusion entirely.

To fully test for a detailed balance, one could try to further analyze the existing data. However, a more direct way to test it would be much more convincing. This could be done by stimulation of CA3 and recording currents in CA1 as before, but at the same time recording the precise number of action potentials elicited in CA3 neurons. The best and maybe only reliable way to do so would be to perform paired recordings, with one (or more) neuron(s) patched in CA3 that are triggered to fire a given number of action potentials, while recording from one (or more) neuron(s) in CA1 in voltage-clamp. This would require considerably more work and result in much lower throughput than the opto-activation experiments, but would allow to address the question of a detailed balance directly and with higher precision.

Conclusion: Overall, I liked the study because it uses experiments to address an important question from theoretical neuroscience. The finding of a balance of excitation and inhibition in CA1 even for (presumably) very few spikes elicited in CA3 is an interesting and unexpected finding which is clearly supported by the data. I’m a bit less convinced about the evidence that this balance is also detailed. However, this does not affect anything of the remaining paper (which constitutes its main part). This main part investigates the supra-linear summation of inputs and the information contained in the timing of inputs (which I do not want to discuss here in any detail).

More context: What is detailed and precise balance? Hennequin et al., Annual Review of Neuroscience (2017). – A broadly accessible explanation of the supra-linear summation presented in the paper: NCBS News.

Posted in electrophysiology, Neuronal activity, Reviews | Tagged , , , | 2 Comments