# Bursting

Bursting, or burst firing, is an extremely diverse general phenomenon of the activation patterns of neurons in the central nervous system and spinal cord where periods of rapid action potential spiking are followed by quiescent periods much longer than typical inter-spike intervals. Bursting is thought to be important in the operation of robust central pattern generators, the transmission of neural codes, and some neuropathologies such as epilepsy. The study of bursting both directly and in how it takes part in other neural phenomena has been very popular since the beginnings of cellular neuroscience and is closely tied to the fields of neural synchronization, neural coding, plasticity, and attention.
Observed bursts are named by the number of discrete action potentials they are composed of: a doublet is a two-spike burst, a triplet three and a quadruplet four. Neurons that are intrinsically prone to bursting behavior are referred to as bursters and this tendency to burst may be a product of the environment or the phenotype of the cell.

## Physiological context
### Overview

Neurons typically operate by firing single action potential spikes in relative isolation as discrete input postsynaptic potentials combine and drive the membrane potential across the threshold. Bursting can instead occur for many reasons, but neurons can be generally grouped as exhibiting input-driven or intrinsic bursting. Most cells will exhibit bursting if they are driven by a constant, subthreshold input and particular cells which are genotypically prone to bursting (called bursters) have complex feedback systems which will produce bursting patterns with less dependence on input and sometimes even in isolation.
In each case, the physiological system is often thought as being the action of two linked subsystems. The fast subsystem is responsible for each spike the neuron produces. The slow subsystem modulates the shape and intensity of these spikes before eventually triggering quiescence.
Input-driven bursting often encodes the intensity of input into the bursting frequency where a neuron then acts as an integrator. Intrinsic bursting is a more specialized phenomenon and is believed to play a much more diverse role in neural computation.

### Slow subsystem

Bursts differ from tonic firing, typically associated with Poisson distributed spike times for a given average firing rate, in that bursting involves a physiological "slow subsystem" that eventually depletes as the bursting continues and then must be replenished before the cell can burst again (compare refractory period). During the bursting event, this slow subsystem modulates the timing and intensity of the emitted spikes and is thought to be important in the computational aspects of the resulting burst pattern. There are many discovered mechanisms of slow subsystems including voltage- and Ca$^{2+}$-gated currents and spiking interplay between dendrites and the cell body.
The slow subsystem also is connected to endogenous bursting patterns in neurons, where the pattern can be maintained completely by internal mechanism without any synaptic input. This process also relies on calcium channels, which depolarize the neuron by allowing an influx of calcium ions. So long as internal calcium ion concentrations remain at an elevated level, the neuron will continue to undergo periods of rapid spiking. However, elevated calcium ion levels also trigger a second messenger cascade within the cell which lower calcium influx and promote calcium efflux and buffering. As calcium concentrations decline, the period of rapid bursting ceases, and the phase of quiescence begins. When calcium levels are low, the original calcium channels will reopen, restarting the process and creating a bursting pattern.

## Statistical detection

In isolation or in mathematical models bursting can be recognized since the environment and state of the neuron can be carefully observed and modulated. When observing neurons in the wild, however, bursting may be difficult to distinguish from normal firing patterns. In order to recognize bursting patterns in these contexts statistical methods are used to determine threshold parameters.
Bursting is characterized by a coefficient of variation (CV) of the interspike intervals (ISI) that is larger than one, or a Fano factor of the spike count that is larger than one, because bursting leads to spike patterns that are more irregular than a Poisson process (which has a CV and Fano factor equal to unity). Alternatively, the serial correlation coefficient of the ISI sequence is positive for bursting patterns, because in this case short ISIs tend to be followed by more short ISIs (at least if the bursts consist of more than two spikes).

## Mathematical models

Neuron behavior is often modeled as single-compartment, non-linear dynamical systems, where the neuron states represent physiological quantities such as membrane voltage, current flow, and the concentrations of various ions intra- and extracellularly. These models most generally take the singularly perturbed form
    fast subsystem: ${\dot {x}}=f(x,u)$
    slow subsystem: ${\dot {u}}=\mu g(x,u)$
where $f$ and $g$ are both Hodgkin–Huxley style relations, ${\dot {x}}$ is a vector representing the cell parameters relevant to the fast subsystem, ${\dot {u}}$ is a vector representing the parameters of the slow modulation subsystem, and $\mu \ll 1$ is the ratio of the time scales between the fast and slow subsystems.
Models of neuron dynamics generally exhibit a number of stable and unstable attractors in phase space which represent resting states. When the system is sufficiently perturbed by input stimuli it may follow a complex return path back to the stable attractor representing an action potential. In bursting neurons, these dynamic spaces bifurcate between quiescent and bursting modes according to the dynamics of the slow system. These two bifurcations may take many forms and the choice of bifurcation both from quiescent to bursting and bursting to quiescent can affect the behavioral aspects of the burster.
The complete classification of quiescent-to-bursting and bursting-to-quiescent bifurcations leads to 16 common forms and 120 possible forms if the dimensionality of the fast subsystem is not constrained. Of the most common 16, a few are well studied.
Common combinations of bifurcations :
| | saddle node on an invariant circle | saddle homoclinic orbit | supercritical Andronov-Hopf | fold limit cycle |
| saddle node (fold) | fold/ circle | fold/ homoclinic | fold/ Hopf | fold/ fold cycle |
| saddle node on an invariant circle | circle/ circle | circle/ homoclinic | circle/ Hopf | circle/ fold cycle |
| supercritical Andronov-Hopf | Hopf/ circle | Hopf/ homoclinic | Hopf/ Hopf | Hopf/ fold cycle |
| subcritical Andronov-Hopf | subHopf/ circle | subHopf/ homoclinic | subHopf/ Hopf | subHopf/ fold cycle |

### Square-wave burster

The fold/homoclinic, also called square-wave, burster is so named because the shape of the voltage trace during a burst looks similar to a square wave due to fast transitions between the resting state attractor and the spiking limit cycle.

## Purposes

Bursting is a very general phenomenon and is observed in many contexts in many neural systems. For this reason it is difficult to find a specific meaning or purpose for bursting and instead it plays many roles. In any given circuit observed bursts may play a part in any or all of the following mechanisms and may have a still more sophisticated impact on the network.

### Synaptic plasticity

Synaptic strengths between neurons follow changes that depend on spike timing and bursting. For excitatory synapses of the cortex, pairing an action potential in the pre-synaptic neuron with a burst in the post-synaptic neuron leads to long-term potentiation of the synaptic strength, while pairing an action potential in the pre-synaptic neuron with a single spike in the post-synaptic neuron leads to long-term depression of the synaptic strength. Such dependence of synaptic plasticity on the spike timing patterns is referred to as burst-dependent plasticity. Burst-dependent plasticity is observed with variations in multiple areas of the brain.

### Multiplexing and routing

Some neurons, sometimes called resonators, exhibit sensitivity for specific input frequencies and fire either more quickly or exclusively when stimulated at that frequency. Intrinsically bursting neurons can use this band-pass filtering effect in order to encode for specific destination neurons and multiplex signals along a single axon. More generally, due to short-term synaptic depression and facilitation specific synapses can be resonant for certain frequencies and thus become viable specific targets for bursting cells. When combined with burst-dependent long-term plasticity, such multiplexing can allow neurons to coordinate synaptic plasticity across hierarchical networks.

### Synchronization

Burst synchronization refers to the alignment of bursting and quiescent periods in interconnected neurons. In general, if a network of bursting neurons is linked it will eventually synchronize for most types of bursting. Synchronization can also appear in circuits containing no intrinsically bursting neurons, however its appearance and stability can often be improved by including intrinsically bursting cells in the network. Since synchronization is related to plasticity and memory via Hebbian plasticity and long-term potentiation the interplay with plasticity and intrinsic bursting is very important.

### Information content and channel robustness

Due to the all-or-nothing nature of action potentials, single spikes can only encode information in their interspike intervals (ISI). This is an inherently low fidelity method of transferring information as it depends on very accurate timing and is sensitive to noisy loss of signal: if just a single spike is mistimed or not properly received at the synapse it leads to a possibly unrecoverable loss in coding. Since intrinsic bursts are thought to be derived by a computational mechanism in the slow subsystem, each can represent a much larger amount of information in the specific shape of a single burst leading to far more robust transmission. Physiological models show that for a given input the interspike and interburst timings are much more variable than the timing of the burst shape itself which also implies that timing between events is a less robust way to encode information.
The expanded alphabet for communication enabled by considering burst patterns as discrete signals allows for a greater channel capacity in neuronal communications and provides a popular connection between neural coding and information theory.

## Example bursting neuron circuits
### Hippocampus

The subiculum, a component of the hippocampal formation, is thought to perform relaying of signals originating in the hippocampus to many other parts of the brain. In order to perform this function, it uses intrinsically bursting neurons to convert promising single stimuli into longer lasting burst patterns as a way to better focus attention on new stimuli and activate important processing circuits. Once these circuits have been activated, the subicular signal reverts to a single spiking mode.

### pre-Bötzinger complex

The pre-Bötzinger complex (preBötC) is located in ventrolateral medulla and is proposed to generate the rhythm underlying inspiratory efforts in mammals. Since the frequency that the lungs need to operate at can vary according to metabolic demand, preBötC activity is modulated over a wide range of frequencies and is able to entrain the respiratory system to meet metabolic demand. While pacemaker neurons do not necessarily require intrinsically bursting neurons the preBötC contains a heterogeneous population of both regular spiking and intrinsically bursting neurons. Intrinsically bursting neurons are thought to make the preBötC oscillations more robust to changing frequencies and the regularity of inspiratory efforts.

### Cerebellar cortex

Cerebellar Purkinje neurons have been proposed to have two distinct bursting modes: dendritically driven, by dendritic Ca$^{2+}$ spikes, and somatically driven, wherein the persistent Na$^+$ current is the burst initiator and the SK K$^+$ current is the burst terminator. Purkinje neurons may utilise these bursting forms in information coding to the deep cerebellar nuclei.