Skip to main content

REVIEW article

Front. Neural Circuits, 19 December 2017
This article is part of the Research Topic Spinal Cord Functions and Plasticity: Lessons from non-mammalian Vertebrates View all 18 articles

Neuronal Population Activity in Spinal Motor Circuits: Greater Than the Sum of Its Parts

  • Department of Neuroscience, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark

The core elements of stereotypical movements such as locomotion, scratching and breathing are generated by networks in the lower brainstem and the spinal cord. Ensemble activities in spinal motor networks had until recently been merely a black box, but with the emergence of ultra-thin Silicon multi-electrode technology it was possible to reveal the spiking activity of larger parts of the network. A series of experiments revealed unexpected features of spinal networks, such as multiple spiking regimes and lognormal firing rate distributions. The lognormality renders the widespread idea of a typical firing rate ± standard deviation an ill-suited description, and therefore these findings define a new arithmetic of motor networks. Focusing on the population activity behind motor pattern generation this review summarizes this advance and discusses its implications.

1. Introduction

We often assume that neurons can be categorized in homogenous and genetically well-defined groups, where each member behaves in more or less the same manner. In spinal motor research, this notion is especially appealing since genetic tools have advanced the field to the forefront of neuroscience (Machado et al., 2015; Bikoff et al., 2016; Gabitto et al., 2016; Sternfeld et al., 2017) and the cellular identity is helpful in the search for potential specialization. However, such genetic reductionism carries weaknesses as well as strengths. The sole focus on cellular subtypes has the risk of failing to see the collective properties of the network. The strategy of isolating an identified population to study its impact on behavior, for instance with genetic knock-out or optogenetics, rests on the radical assumption that the impact of this population only has a feedforward influence. Nevertheless, circuits operate by a delicate interaction between neurons of different genetic origin most likely with pervasive recurrent connectivity, where it may be fruitless to assign a role to any one member (Yuste, 2015). The intricacy of control of one population by another has recently become evident from control theory, i.e., the study of manipulation of dynamics on complex networks. The behavior of such a network is difficult to control by manipulation of selected nodes and this strategy often has counter-intuitive effects (Liu et al., 2011). A complementary investigation of the collective population dynamics is therefore recommendable in concert with genetic analysis. However, population spiking activity is a challenge not only to analyze, but also to acquire and have only rarely been done in spinal motor systems (Berg et al., 2009; Auyong et al., 2011). To achieve recordings of the population spiking activity requires multiple electrode arrays in preparations, which are mechanically stable. These challenges are part of the reason for the scarce literature and slow progress on collective properties of spinal motor networks.

Arithmetic of Population Activity

In a recent set of experiments, however, the mechanical stability of the turtle spinal cord (Stein, 2005) was used to investigate the neuronal population spiking activity associated with motor pattern generation within spinal circuitry (Vestergaard and Berg, 2015; Petersen and Berg, 2016). Custom-design silicon electrode arrays were used to study populations of interneurons and motoneurons in the medial-ventral portions of the lumbar region (Petersen and Berg, 2017), in order to probe the concerted activity. This yielded insight to the complexity of population spiking during motor behavior (rhythmic hindlimb scratching Petersen and Berg, 2016). Neurons were not easily sorted into discrete classes, but was rather statistically distributed as a whole, without any identifiable clustering of simple spiking behavior. Rather, the heterogeneous population revealed firing rates, which were scattered in a smooth and continuous fashion with no indication of a multimodal distribution. Further, the distribution was not the expected normal distribution with a symmetric spread around some mean, but was instead strongly skewed with a fat-tail and approximately a lognormal distribution, i.e., a normal distribution on a logarithmic x-axis (Figure 1F). Where does this lognormality come from? As it turns out, there is a simple explanation for the skewed ‘participation’ among spinal neurons. Lognormality can arise by transforming normally distributed variables with a non-linearity (Roxin et al., 2011). Here, the non-linearity is the transformation of synaptic input to a firing rate output, i.e., the gain-function, and the normally-distributed variable is the synaptic input. The net synaptic input represents the sum of excitation and inhibition, whose relative amount (or ratio) varies from neuron to neuron. Some neurons receive more excitatory than inhibitory connections (excitation-dominated) and others receive more inhibition (inhibition-dominated). Although the distribution of synaptic connections is unknown (Koulakov et al., 2009; Roxin, 2011) the resulting distribution of membrane potentials over the population is Gaussian (normal) (Figure 1). When passing such normal distribution through a non-linear response function the firing rate distribution becomes skewed. Whereas normal-distributions are typically associated with the addition of random variation (synaptic currents and potentials are roughly additive), lognormal-distributions are results of multiplicative effects (Limpert et al., 2001). The widely held belief that firing rates across the population are scattered with a mean value ± standard deviation has thus proven to be unsuitable. The network has an embedded supra-linear element. As a consequence, the addition of multiple inputs amounts to more than the linear sum of these inputs. These observations define a novel arithmetic of spinal population activity, which have previously escaped attention. As we will see below, the results also tie together population dynamics, regimes of spiking with stability and flexibility and reconciles neuronal activity in spinal cord with that found in cortex and elsewhere (Mizuseki and Buzsáki, 2013; Wohrer et al., 2013; Buzsáki and Mizuseki, 2014).

FIGURE 1
www.frontiersin.org

Figure 1. Lognormality in firing rates across neuronal population stems from a nonlinearity of the gain-curve. (A) Some neurons receive more inhibition (red) and others receive more excitation (green, B). This creates various mean Vm (cf. bottom traces) with a normal distribution across the population. (C) A normal-distributed input will be transformed to a normal-distributed output, if the gain-curve is linear. If the gain curve is exponential the output distribution will be skewed (D). (E) Population spiking activity (~300 neurons, raster at bottom) was recorded during rhythmic motor behavior (4 motor nerves shown at top). Cycle indicated by green shaded region. Neurons sorted according to phase. (F) Spike count rate distribution across the population is approximately lognormal. (G) There is a nonlinear relationship between sub-threshold membrane potential and firing rate. (H) Subthreshold membrane potential strongly fluctuates with a normal distribution. (I) Distribution of the mean membrane potential subtracted threshold and normalized to standard deviation ((VmVthres)/σ) is normal across population. The typical distance of the mean to the threshold is 3σ (arrow). Adapted with permission (Roxin et al., 2011; Vestergaard and Berg, 2015; Petersen and Berg, 2016).

Input Is “Normal” and Output Is “Lognormal”

Neurons receive a mixture of excitatory and inhibitory synaptic connections. The precise number of connections and strength of the contacts are subject to statistical variation. Some neurons have more inhibition than excitation (Figure 1A) whereas other neurons are biased toward more excitation (Figure 1B). As a result, the mean synaptic current and membrane potential are not the same from neuron to neuron, but are rather distributed with a mean and a spread. What is the shape of this distribution? Since synaptic inputs are additive we expect a normal distribution according to the central limit theorem, both regarding net synaptic current as well as membrane potential in between spikes (Roxin et al., 2011).

Let us assume for now that the population distribution of means is normally distributed and consider two possible neuronal response functions (gain-curves). The most generic gain-curve is linear and will transform the normal distribution into a normal firing-rate distribution across the population (Figure 1C). However, if the gain-curve is non-linear the transformation will become skewed (Figure 1D). An exponential gain-curve will give an exact lognormal distribution. Other non-linearities, e.g., a power-law, will also result in skewed and “lognormal-like” distributions, although they will not be strictly lognormal. Therefore firing rate distributions are intricately linked to input distributions as well as the gain-curve.

How can a gain-curve become supra-linear? The traditional gain-curve, i.e., the frequency response to an injected current (F-I-curve) is zero below rheobase and linear above (Gerstner et al., 2014). This description is based on early intracellular measurements performed on neurons in the absence of synaptic input. Nevertheless, real input consists of rapid fluctuations and increase in synaptic conductance, which will not only change the rheobase (Grigonis et al., 2016), but also bend the gain-curve (Silver, 2010). A combination of a threshold and fluctuations in the membrane potential (or current) will thus result in an ‘expansive’ non-linearity (Hansel and van Vreeswijk, 2002; Miller and Troyer, 2002; Priebe and Ferster, 2008). The strength of such non-linearity has an inverse relationship with the degree of fluctuations (Vestergaard and Berg, 2015). Until recently nevertheless, it had not been addressed how these different elements affect the population activity in rhythm generating circuits of the spinal cord.

Lognormality of Spinal Population Activity

The neuronal ensemble activity of the rhythm generating circuitry of the lumbar spinal cord was therefore investigated in order to address the question (Petersen and Berg, 2016). The rhythmic hindlimb scratching of the turtle preparation was used as a model for stereotypical movement (Stein et al., 1998; Hao et al., 2014). This movement can be generated purely by the lumbar spinal cord network without the confounding factors of supraspinal input. The lack of corticospinal and bulbuspinal input also deprive the spinal circuits for certain excitatory and inhibitory input as well as neuromodulatory input such as serotonergic, adrenergic and peptidergic input. Nevertheless, several distinct behaviors can be evoked by somatic touch with no overt difference from natural behaviors (Keifer and Stein, 1983; Stein, 2005). This preparation (Petersen and Berg, 2017) offers better mechanical stability compared to in-vivo and the anoxia-resistance of this adult reptilian nervous system permits more physiologically intact activity than similar experiments in the mammalian nervous system. Thus, it was possible to monitor the spiking activity of hundreds of neurons (~300) while simultaneously recording the neuronal intracellular activity together with the motor nerve output to various hindlimb muscles (Figure 1E). The distribution of firing rates across the population closely resembled a lognormal distribution-note the log-scale (Figure 1F). Such lognormal-like distributions are also present in various other parts of the nervous system (O'Connor et al., 2010; Mizuseki and Buzsáki, 2013; Buzsáki and Mizuseki, 2014) and could represent a ubiquitous feature of neuronal networks. In order to verify the cause of this lognormality, as discussed above, the gain-curve was estimated in the subthreshold spiking (Figure 1G). The gain-curve had an expansive non-linearity similar to an exponential as well as a power-law. This non-linearity is likely due to the presence of synaptic fluctuations (Figure 1H) combined with a threshold mechanism.

Are the fluctuations, and therefore the gain-curves, different for neuron to neuron? Yes, but it was demonstrated that the non-linearity had a rather weak inverse dependence on the size of the synaptic fluctuations (Vestergaard and Berg, 2015; Petersen and Berg, 2016). This weak effect suggests that although different neurons have different levels of fluctuations in input, the non-linearity of the gain-curve is largely conserved across the population of neurons. As a consequence, the normally distributed mean membrane potential with the non-linear gain-curve together offer an explanation of the lognormality in firing rate distribution. To further verify the constancy of gain-curve across members, the membrane potentials were analyzed across a population of neurons. Since neurons have different electrotonic morphologies and thus fluctuations and thresholds, the distance of the mean membrane potential from the threshold was obtained and normalized by the size of their fluctuations (σ). This distribution ((VmVthres)/σ) turned out to also be normal across population (Figure 1I). Interestingly, this suggests a universal normalization of membrane potential distribution for the individual neuron to have a preferred distance to threshold of 3σ (arrow, Figure 1I). This normalization also predicts a preserved gain-curve across the population. In toto, these two elements, the normal input and the expansive non-linearity, lend support to a mechanistic explanation of the lognormal firing rate distribution (Figure 1D).

Regimes of Spiking: Regular and Irregular

So far, the analysis has focused on the sub-threshold spiking, where the gain-curve is non-linear and the spiking is driven by subthreshold fluctuations. How typical is this type of spiking and is there another regime of spiking? Indeed, there is another type of spiking referred to as mean-driven spiking, where the mean input current is larger than the rheobase (Renart et al., 2007; Gerstner et al., 2014). Here, the gain-curve is linear or even sub-linear (Figure 2A). The firing-rate distribution is symmetric rather than skewed (Figure 1C) although some spinal neurons have shown an increase in the gain-curve slope for much higher injected currents, i.e., the ‘secondary firing range’ due to persistent inward currents (Heckman et al., 2008; Meehan et al., 2010). The inter-spike intervals here are more affected by the after-hyperpolarizations (Matthews, 1996) and therefore the spiking is more regular and at higher firing rate (Figure 2B). This is in contrast with the fluctuation-driven spiking, which is irregular and at lower rates (Figure 2C). Some neurons spent most of the time below threshold (96% in the shown sample neuron), even without counting the inter-burst intervals, whereas other neurons spent much less time (35% in the shown sample neuron, Figure 2D). The majority of the neurons spent most of the time below threshold, i.e., in the fluctuation-driven regime (right, Figure 2D), which may be relevant for muscular control since this regime represents 90% of the force modulation in mice (Manuel and Heckman, 2011). Since the threshold has a dependence on the firing rate (Grigonis and Alaburda, 2017) the initial threshold, i.e., the threshold of the first spike in a trial, was used in this analysis. It should also be noted that the impact of intrinsic properties, such as the spike frequency adaptation, on the spiking dynamics would be difficult to assess during the these intense synaptic input. Nevertheless, intrinsic properties work on slower timescales and therefore they mainly affect spiking in a regular manner. Hence, irregularity is an indicator that can be utilized to quantify the fraction of an ensemble found in either of the regimes using spike times from extracellular recordings. Multi-electrode arrays record the extracellular potential of many neurons, and can thus more easily capture data from a larger population. Such experiments show that the neuronal population is almost equally divided between the two regimes (Petersen and Berg, 2016; see below, Figure 4).

FIGURE 2
www.frontiersin.org

Figure 2. Two regimes of spiking during motor behavior. (A) The gain curve of a single neuron has two regimes. Fluctuations in input allow sub-threshold spiking and cause a smooth transition across threshold (vertical broken line). Larger noise causes more smooth transition. (B) A spinal neuron during rhythmic motor activity has mean-driven spiking, i.e., the mean membrane potential (blue line) is above threshold (broken line, left). Trace selected from a locomotion cycle (right) where most spikes have low irregularity (CV2, green trace). (C) Another neuron has fluctuation-driven spiking, i.e., the mean membrane potential (blue) is below threshold (broken line, left). Whole cycle indicates that majority of spikes are fluctuation-driven (right) and high irregularity (green curve). (D) Amount of time spent in the two regimes is quantified by the cumulative time spent below the initial threshold (left). Neuron in (B) has only 35% of time below threshold whereas neuron in (C) has 96%. Most of the neurons had a majority of time below threshold, i.e., in fluctuation-driven regime (total n = 68 neurons, right). Adapted with permission (Petersen and Berg, 2016).

Absent Clustering and Cellular Identity

The identity of interneurons can be delineated by their genetic origin (Goulding, 2009; Hinckley et al., 2015; Kiehn, 2016). The spinal cord has the fortunate architecture that cells are developmentally segregated primarily in the dorsoventral (DV) axis of its gray matter. Therefore it is possible to probe the physiological diversity of spinal neurons by recording along the DV-axis. Custom-designed multi-channel electrodes (~256) were thus implanted in the lumbar spinal cord (Petersen and Berg, 2016) and the ensemble activity recorded (Figures 3A,B). Using trilateration and source separation combined with the shank depth it was possible to tease apart their location in the DV-axis (Figures 3C,D). The irregularity of individual members, as quantified using the CV2-measure (Holt et al., 1996), demonstrated no difference in the distribution for different location in the DV-axis (Figures 3E,F). Their distributions had simple Gaussian shapes with variance much larger than the difference in mean. This suggests that the spiking activity was equally irregular for neurons despite diversity in genetic identity. Similar observation was obtained for the firing rate distributions. These distributions were all lognormal and independent of location, indicating an absence of location-specific clustering of these simple features. Hence, there was no obvious link between genetic identity and these simple neurophysiological aspects that were characterized. This is not in conflict with another investigation that demonstrated less rhythmicity in more dorsal units (Berkowitz, 2001).

FIGURE 3
www.frontiersin.org

Figure 3. Spike patterns lack topological clustering in the dorso-ventral axis. (A) The multi-electrode arrays are inserted in the spinal cord gray matter along the axis of high cellular segregation (D-V direction). (B) Spiking pattern of sorted units (~300 neurons) concurrent with motor nerves during rhythmic movement. (C) Layout of electrodes on a shank. (D) The recording sites (black dots) are used to locate the strongest sources along the shank using trilateration. Three sample units indicated in colors, all units in gray. (E) The irregularity of spiking vs. depths in DV-axis. Each circle is the mean CV2 of a neuron. (F) The population distribution of mean spiking irregularity (CV2) as a function of depth in DV-axis. A difference between the distributions could not be detected. Adapted with permission (Petersen and Berg, 2016).

Purpose of Lognormality in Spinal Population Activity?

A long-standing question in theoretic neuroscience has been how neuronal networks maintain self-sustained activity while avoiding run-away excitation (Roxin et al., 2004; Kumar et al., 2008; Renart et al., 2010; Vogels et al., 2011). This is also an open question in spinal research. How do spinal circuits generate activity so reliably with high sensitivity and yet stay clear of seizures or saturation? Neural circuits most likely reach this ‘Goldilocks zone (Humphries, 2016) by stabilizing excitation via recurrent inhibition. The purpose of such an arrangement could be to ensure sensitivity to smaller input and curbing the response to strong input. The gain-curve has a sigmoidal shape, where the left part has a supra-linear summation that can amplify weak input (Rubin et al., 2015). The input can both represent sensory input or internally generated activity within the network itself. Stronger input will move up the curve to sub-linear summation, which will attenuate and stabilize activity. Such enhancement of weak input while curbing strong activity is beneficial for extending the dynamic range of the network operation. The network will always be able to increase output, although it is strenuous for stronger input. The extension of dynamical range can thus be accomplished by the sigmoidal shape of the gain curve, where each part represents the two spiking regimes. How many neurons are in the either the supra- or sub-linear part of the gain curve? The supra-linear part is exactly the fluctuation-driven regime and the sub-linear part represents the mean-driven spiking. It was therefore possible to answer this question by quantifying the fraction of the population in either of these regimes via the irregularity of spiking. Using the CV2 metric for local irregularity the time of each neuron spiking irregularly (CV2 > 0.5) the inverse cumulative distribution for the whole population measures the fraction of neurons that spike at least a certain fraction of time (x-axis) in the irregular regime (Figure 4A). The time that half of the population (broken line) spend in the fluctuation-driven regime was remarkably close to 50%. This number was conserved across animals (TIF50, inset). Similar analysis was performed for the number of fluctuation-driven spikes (spikes in fluctuation regime, SIF50) rather than time, since these might be different when the firing rate changes. Nevertheless, the analysis gave qualitatively similar results (Figure 4B). This could indicate a kind of homeostasis for network spiking activity in redundancy with cellular mechanisms, such as homeostatic plasticity (Kline et al., 2007; Pozo and Goda, 2010). Irregularity of spiking of motor neurons is not expected to be a great disadvantage since it introduces only a mild decrease in steadiness of muscle force, due to temporal summation within the muscle fibers and spatial summation of over neighboring muscle fibers, which both work as low-pass filters (Dideriksen et al., 2012).

FIGURE 4
www.frontiersin.org

Figure 4. Fraction of population vs. time and spikes in fluctuation-driven regime. (A) Fraction of minimal time spent in the fluctuation regime (x-axis) vs. fraction of neuronal population. About half of the population spends at least 50% in with CV2 > 0.5 (broken line), i.e., Time In Fluctuation driven-regime, TIF50. Similar numbers for TIF50 was found for all animals tested (inset, n = 5, sample animal indicated ⋆). (B) Fraction of neurons vs. spikes in fluctuation-driven regime, and the similar metric Spikes in Fluctuation-driven regime SIF50. Similar numbers slightly above 50% were found for all 5 animals (inset). Adapted with permission (Petersen and Berg, 2016).

Future Directions

An outstanding question in spinal motor research is how rhythm generation takes place in the spinal circuits (Whelan, 2010; McLean and Dougherty, 2015; Sternfeld et al., 2017). The recent findings of the ensemble activity, which I discussed in this review, do not explain rhythm generation, only certain concerted properties associated with it. Pacemaker neurons are widely believed to be responsible for rhythm generation or at least to be involved as a supporting mechanism for e.g., increasing the robustness of the oscillation (Purvis et al., 2007). Nevertheless, new findings in the respiratory field has cast doubt on pacemaker-hypothesis (Feldman et al., 2013; Feldman and Kam, 2015) although it is an on-going debate. The finding presented here that a larger portion of the population is in the fluctuation-driven regime is interesting in this context, since irregular spiking most likely originates from balanced excitation and inhibition (E/I) (Shadlen and Newsome, 1998; Petersen and Berg, 2016). Balanced E/I leaves the membrane potential close to threshold while increasing the membrane conductance sometimes several fold (Destexhe et al., 2003; Alaburda et al., 2005; Berg et al., 2007). Large increases in conductance are consequential since it can shunt the intrinsic properties (Kolind et al., 2012) and thus in principle “confiscate” the specialization of neurons (Berg and Hounsgaard, 2009). Nevertheless, conductance increase is more likely to only distort the kinetics of the pacemaker dynamics, yet still present a challenge to the pacemaker hypothesis. There are alternatives to the pacemaker-hypothesis such as network mechanisms (Grillner, 2006; Brocard et al., 2010; Li, 2011; Ramirez et al., 2011), i.e., an emerging feature on a level above single neurons (Yuste, 2015), but this hypothesis needs refinement. If the rhythm is generated by the population rather than individual pacemaker cells, ensemble recordings—as presented in this review—are crucial for understanding rhythm generation (Shenoy et al., 2013; Bruno et al., 2015, 2017). Therefore, the fact that experiments on population spiking activity are so rare suggests a research area worth while of pursuing. Optical imaging of the calcium activity of neuronal populations in the spinal cord (Ritter et al., 2001; Wilson et al., 2007; Kwan et al., 2010; Johannssen and Helmchen, 2013; Renninger and Orger, 2013; Weinger et al., 2015) are an important tool to include, especially with DNA-encoded calcium indicators, which can specifically label subpopulations of neurons (Muto et al., 2011; Machado et al., 2015). Obstacles in optical imaging are the lower temporal resolution as well as the relativity in signal (ΔF/F), which each prevents an analysis, similar to the one presented here, to be performed. Functional connectivity in association with multi-electrode recordings can assist in mapping the network in spinal cord behind motor pattern generation, which is largely unknown. Although multi-electrode recordings in mammalian spinal cord remain a challenge (Berg et al., 2009; Auyong et al., 2011) a natural next step is to verify these findings in the mammalian spinal cord.

Author Contributions

The author confirms being the sole contributor of this work and approved it for publication.

Funding

Funded by Danish Council for Independent Research Medical Sciences and the Dynamical Systems Interdisciplinary Network, University of Copenhagen.

Conflict of Interest Statement

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Thanks to Peter Petersen and the reviewers for comments on and earlier version of the manuscript.

References

Alaburda, A., Russo, R., MacAulay, N., and Hounsgaard, J. (2005). Periodic high-conductance states in spinal neurons during scratch-like network activity in adult turtles. J. Neurosci. 25, 6316–6321. doi: 10.1523/JNEUROSCI.0843-05.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Auyong, N., Ollivier-lanvin, K., and Lemay, M. A. (2011). Population spatiotemporal dynamics of spinal intermediate zone interneurons during air-stepping in adult spinal cats. J Neurophysiol 106, 1943–1953. doi: 10.1152/jn.00258.2011.

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg, R. W., Alaburda, A., and Hounsgaard, J. (2007). Balanced inhibition and excitation drive spike activity in spinal half-centers. Science 315, 390–393. doi: 10.1126/science.1134960

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg, R. W., Chen, M.-T., Huang, H.-C., Hsiao, M.-C., and Cheng, H. (2009). A method for unit recording in the lumbar spinal cord during locomotion of the conscious adult rat. J. Neurosci. Method. 182, 49–54. doi: 10.1016/j.jneumeth.2009.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg, R. W., and Hounsgaard, J. (2009). Signaling in large-scale neural networks. Cogn. Proces. 10(Suppl. 1), S9–S15. doi: 10.1007/s10339-008-0238-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Berkowitz, A. (2001). Rhythmicity of spinal neurons activated during each form of fictive scratching in spinal turtles. J. Neurophysiol. 86, 1026–1036. doi: 10.1152/jn.2001.86.2.1026

PubMed Abstract | CrossRef Full Text | Google Scholar

Bikoff, J. B., Gabitto, M. I., Rivard, A. F., Drobac, E., Machado, T. A., Miri, A., et al. (2016). Spinal inhibitory interneuron diversity delineates variant motor microcircuits. Cell 165, 207–219. doi: 10.1016/j.cell.2016.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Brocard, F., Tazerart, S., and Vinay, L. (2010). Do pacemakers drive the central pattern generator for locomotion in mammals? Neuroscientist 16, 139–155. doi: 10.1177/1073858409346339

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruno, A. M., Frost, W. N., and Humphries, M. D. (2015). Modular deconstruction reveals the dynamical and physical building blocks of a locomotion motor program. Neuron 86, 304–318. doi: 10.1016/j.neuron.2015.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruno, A. M., Frost, W. N., and Humphries, M. D. (2017). A spiral attractor network drives locomotion in Aplysia. eLIFE 6:e27342. doi: 10.7554/eLife.27342

PubMed Abstract | CrossRef Full Text | Google Scholar

Buzsáki, G., and Mizuseki, K. (2014). The log-dynamic brain: how skewed distributions affect network operations. Nat. Rev. Neurosci. 15, 264–278. doi: 10.1038/nrn3687

PubMed Abstract | CrossRef Full Text | Google Scholar

Destexhe, A., Rudolph, M., and Paré, D. (2003). The high-conductance state of neocortical neurons in vivo. Nat. Rev. Neurosci. 4, 739–751. doi: 10.1038/nrn1198

PubMed Abstract | CrossRef Full Text | Google Scholar

Dideriksen, J. L., Negro, F., Enoka, R. M., and Farina, D. (2012). Motor unit recruitment strategies and muscle properties determine the influence of synaptic noise on force steadiness. J. neurophysiol. 107, 3357–3369. doi: 10.1152/jn.00938.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Feldman, J. L., Del Negro, C. A., and Gray, P. A. (2013). Understanding the rhythm of breathing: so near, yet so far. Annu. Rev. Physiol. 75, 423–452. doi: 10.1146/annurev-physiol-040510-130049

PubMed Abstract | CrossRef Full Text | Google Scholar

Feldman, J. L., and Kam, K. (2015). Facing the challenge of mammalian neural microcircuits: taking a few breaths may help. J. Physiol. 593, 3–23. doi: 10.1113/jphysiol.2014.277632

PubMed Abstract | CrossRef Full Text | Google Scholar

Gabitto, M. I., Pakman, A., Bikoff, J. B., Abbott, L. F., Jessell, T. M., and Paninski, L. (2016). Bayesian sparse regression analysis documents the diversity of spinal inhibitory interneurons. Cell 165, 220–233. doi: 10.1016/j.cell.2016.01.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerstner, W., Kistler, W. M., Naud, R., and Paninski, L. (2014). Neuronal Dynamics: From Single Neurons to Networks and Models of Cognition. New York, NY: Cambridge University Press.

Google Scholar

Goulding, M. (2009). Circuits controlling vertebrate locomotion: moving in a new direction. Nat. Rev. Neurosci. 10, 507–518. doi: 10.1038/nrn2608

PubMed Abstract | CrossRef Full Text | Google Scholar

Grigonis, R., and Alaburda, A. (2017). Spike threshold dynamics in spinal motoneurons during scratching and swimming. J. Physiol. 595, 5843–5855. doi: 10.1113/JP274434

PubMed Abstract | CrossRef Full Text | Google Scholar

Grigonis, R., Guzulaitis, R., Buisas, R., and Alaburda, A. (2016). The influence of increased membrane conductance on response properties of spinal motoneurons. Brain Res. 1648, 110–118. doi: 10.1016/j.brainres.2016.07.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Grillner, S. (2006). Biological pattern generation: the cellular and computational logic of networks in motion. Neuron 52, 751–766. doi: 10.1016/j.neuron.2006.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansel, D., and van Vreeswijk, C. (2002). How noise contributes to contrast invariance of orientation tuning in cat visual cortex. J. Neurosci. 22, 5118–5128.

PubMed Abstract | Google Scholar

Hao, Z., Meier, M., and Berkowitz, A. (2014). Rostral spinal cord segments are sufficient to generate a rhythm for both locomotion and scratching but affect their hip extensor phases differently. J. Neurophysiol. 112, 147–155. doi: 10.1152/jn.00119.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Heckman, C. J., Johnson, M., Mottram, C., and Schuster, J. (2008). Persistent inward currents in spinal motoneurons and their influence on human motoneuron firing patterns. Neuroscientist 14, 264–275. doi: 10.1177/1073858408314986

PubMed Abstract | CrossRef Full Text | Google Scholar

Hinckley, C. A., Alaynick, W. A., Gallarda, B. W., Hayashi, M., Hilde, K. L., Driscoll, S. P., et al. (2015). Spinal locomotor circuits develop using hierarchical rules based on motorneuron position and identity. Neuron 87, 1008–1021. doi: 10.1016/j.neuron.2015.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Holt, G. R., Softky, W. R., Koch, C., and Douglas, R. J. (1996). Comparison of discharge variability in vitro and in vivo in cat visual cortex neurons. J. Neurophysiol. 75, 1806–1814.

PubMed Abstract | Google Scholar

Humphries, M. D. (2016). The Goldilocks zone in neural circuits. eLife 5:e22735. doi: 10.7554/eLife.22735

PubMed Abstract | CrossRef Full Text | Google Scholar

Johannssen, H. C., and Helmchen, F. (2013). Two-photon imaging of spinal cord cellular networks. Exp. Neurol. 242, 18–26. doi: 10.1016/j.expneurol.2012.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Keifer, J., and Stein, P. S. (1983). In vitro motor program for the rostral scratch reflex generated by the turtle spinal cord. Brain Res. 266, 148–151.

PubMed Abstract | Google Scholar

Kiehn, O. (2016). Decoding the organization of spinal circuits that control locomotion. Nat. Rev. Neurosci. 17, 224–238. doi: 10.1038/nrn.2016.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kline, D. D., Ramirez-Navarro, A., and Kunze, D. L. (2007). Adaptive depression in synaptic transmission in the nucleus of the solitary tract after in vivo chronic intermittent hypoxia: evidence for homeostatic plasticity. J. Neurosci. 27, 4663–4673. doi: 10.1523/JNEUROSCI.4946-06.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolind, J., Hounsgaard, J., and Berg, R. W. (2012). Opposing effects of intrinsic conductance and correlated synaptic input on Vm-fluctuations during network activity. Front. Comput. Neurosci. 6:40. doi: 10.3389/fncom.2012.00040

PubMed Abstract | CrossRef Full Text | Google Scholar

Koulakov, A. A., Hromádka, T., and Zador, A. M. (2009). Correlated connectivity and the distribution of firing rates in the neocortex. J. Neurosci. 29, 3685–3694. doi: 10.1523/JNEUROSCI.4500-08.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, A., Schrader, S., Aertsen, A., and Rotter, S. (2008). The high-conductance state of cortical networks. Neural. Comput. 20, 1–43. doi: 10.1162/neco.2008.20.1.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwan, A. C., Dietz, S. B., Zhong, G., Harris-Warrick, R. M., and Webb, W. W. (2010). Spatiotemporal dynamics of rhythmic spinal interneurons measured with two-photon calcium imaging and coherence analysis. J. Neurophysiol. 104, 3323–3333. doi: 10.1152/jn.00679.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W.-C. (2011). Generation of locomotion rhythms without inhibition in vertebrates: the search for pacemaker neurons. Int. Compar. Biol. 51, 879–889. doi: 10.1093/icb/icr021

PubMed Abstract | CrossRef Full Text | Google Scholar

Limpert, E., Stahel, W. A., and Abbt, M. (2001). Log-normal distributions across the Sciences: keys and clues. BioScience 51:341. doi: 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2

CrossRef Full Text | Google Scholar

Liu, Y.-Y., Slotine, J.-J., and Barabási, A.-L. (2011). Controllability of complex networks. Nature 473, 167–173. doi: 10.1038/nature10011

PubMed Abstract | CrossRef Full Text | Google Scholar

Machado, T. A., Pnevmatikakis, E., Paninski, L., Jessell, T. M., and Miri, A. (2015). Primacy of flexor locomotor pattern revealed by ancestral reversion of motor neuron identity. Cell 162, 338–350. doi: 10.1016/j.cell.2015.06.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Manuel, M., and Heckman, C. J. (2011). Adult mouse motor units develop almost all of their force in the subprimary range: a new all-or-none strategy for force recruitment? J Neurosci 31, 15188–15194. doi: 10.1523/JNEUROSCI.2893-11.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Matthews, P. B. (1996). Relationship of firing intervals of human motor units to the trajectory of post-spike after-hyperpolarization and synaptic noise. J. Physiol. 492(Pt 2), 597–628.

PubMed Abstract | Google Scholar

McLean, D. L., and Dougherty, K. J. (2015). Peeling back the layers of locomotor control in the spinal cord. Curr. Opin. Neurobiol. 33, 63–70. doi: 10.1016/j.conb.2015.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Meehan, C. F., Sukiasyan, N., Zhang, M., Nielsen, J. B., and Hultborn, H. (2010). Intrinsic properties of mouse lumbar motoneurons revealed by intracellular recording in vivo. J. Neurophysiol. 103, 2599–2610. doi: 10.1152/jn.00668.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, K. D., and Troyer, T. W. (2002). Neural noise can explain expansive, power-law nonlinearities in neural response functions. J. neurophysiol. 87, 653–659. doi: 10.1152/jn.00425.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Mizuseki, K., and Buzsáki, G. (2013). Preconfigured, skewed distribution of firing rates in the hippocampus and entorhinal cortex. Cell Reports 4, 1010–1021. doi: 10.1016/j.celrep.2013.07.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Muto, A., Ohkura, M., Kotani, T., Higashijima, S.-I., Nakai, J., and Kawakami, K. (2011). Genetic visualization with an improved GCaMP calcium indicator reveals spatiotemporal activation of the spinal motor neurons in zebrafish. Proc. Natl. Acad. Sci. U.S.A. 108, 5425–5430. doi: 10.1073/pnas.1000887108

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Connor, D. H., Peron, S. P., Huber, D., and 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersen, P. C., and Berg, R. W. (2016). Lognormal firing rate distribution reveals prominent fluctuation-driven regime in spinal motor networks. eLife 5:e18805. doi: 10.7554/eLife.18805

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersen, P. C., and Berg, R. W. (2017). Spinal cord preparation from adult red-eared turtles for electrophysiological recordings during motor activity. Bio. Protocol. 7:e2381. doi: 10.21769/BioProtoc.2381

CrossRef Full Text | Google Scholar

Pozo, K., and Goda, Y. (2010). Unraveling mechanisms of homeostatic synaptic plasticity. Neuron 66, 337–351. doi: 10.1016/j.neuron.2010.04.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Priebe, N. J., and Ferster, D. (2008). Inhibition, spike threshold, and stimulus selectivity in primary visual cortex. Neuron 57, 482–497. doi: 10.1016/j.neuron.2008.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Purvis, L. K., Smith, J. C., Koizumi, H., and Butera, R. J. (2007). Intrinsic bursters increase the robustness of rhythm generation in an excitatory network. J. Neurophysiol. 97, 1515–1526. doi: 10.1152/jn.00908.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramirez, J. M., Koch, H., Garcia, A. J., Doi, A., and Zanella, S. (2011). The role of spiking and bursting pacemakers in the neuronal control of breathing. J. Biol. Phys. 37, 241–261. doi: 10.1007/s10867-011-9214-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Renart, A., de la Rocha, J., Bartho, P., Hollender, L., Parga, N., Reyes, A. D., et al. (2010). The asynchronous state in cortical circuits. Science 327, 587–590. doi: 10.1126/science.1179850

PubMed Abstract | CrossRef Full Text | Google Scholar

Renart, A., Moreno-Bote, R., Wang, X.-J., and Parga, N. (2007). Mean-driven and fluctuation-driven persistent activity in recurrent networks. Neural. Comput. 19, 1–46. doi: 10.1162/neco.2007.19.1.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Renninger, S. L., and Orger, M. B. (2013). Two-photon imaging of neural population activity in zebrafish. Methods 62, 255–267. doi: 10.1016/j.ymeth.2013.05.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Ritter, D. A., Bhatt, D. H., and Fetcho, J. R. (2001). In vivo imaging of zebrafish reveals differences in the spinal networks for escape and swimming movements. J. Neurosci. 21, 8956–8965.

PubMed Abstract | Google Scholar

Roxin, A. (2011). The role of degree distribution in shaping the dynamics in networks of sparsely connected spiking neurons. Front. Comput. Neurosci. 5:8. doi: 10.3389/fncom.2011.00008

PubMed Abstract | CrossRef Full Text | Google Scholar

Roxin, A., Brunel, N., Hansel, D., Mongillo, G., and 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Roxin, A., Riecke, H., and Solla, S. A. (2004). Self-sustained activity in a small-world network of excitable neurons. Phys. Rev. Lett. 92, 198101–198101. doi: 10.1103/PhysRevLett.92.198101

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubin, D. B., Hooser, S. D. V., and Miller, K. D. (2015). The stabilized supralinear network : a unifying circuit motif underlying multi-input integration in sensory cortex. Neuron 85, 402–417. doi: 10.1016/j.neuron.2014.12.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Shadlen, M. N., and Newsome, W. T. (1998). The variable discharge of cortical neurons: implications for connectivity, computation, and information coding. J. Neurosci. 18, 3870–3896.

PubMed Abstract | Google Scholar

Shenoy, K. V., Sahani, M., and Churchland, M. M. (2013). Cortical control of arm movements: a dynamical systems perspective. Annu. Rev. Neurosci. 36, 337–359. doi: 10.1146/annurev-neuro-062111-150509

PubMed Abstract | CrossRef Full Text | Google Scholar

Silver, R. A. (2010). Neuronal arithmetic. Nat. Rev. Neurosci. 11, 474–489. doi: 10.1038/nrn2864

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, P. (2005). Neuronal control of turtle hindlimb motor rhythms. J. Comp. Physiol. A. Neuroethol. Sens. Neural. Behav. Physiol. 191, 213–229. doi: 10.1007/s00359-004-0568-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, P. S., McCullough, M. L., and Currie, S. N. (1998). Spinal motor patterns in the turtle. Ann. N.Y. Acad. Sci. 860, 142–154.

PubMed Abstract | Google Scholar

Sternfeld, M. J., Hinckley, C. A., Moore, N. J., Pankratz, M. T., Hilde, K. L., Driscoll, S. P., et al. (2017). Speed and segmentation control mechanisms characterized in rhythmically-active circuits created from spinal neurons produced from genetically-tagged embryonic stem cells. eLife 6:e21540. doi: 10.7554/eLife.21540

PubMed Abstract | CrossRef Full Text | Google Scholar

Vestergaard, M., and Berg, R. W. (2015). Divisive gain modulation of motoneurons by inhibition optimizes muscular control. J. Neurosci. 35, 3711–3723. doi: 10.1523/JNEUROSCI.3899-14.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogels, T. P., Sprekeler, H., Zenke, F., Clopath, C., and Gerstner, W. (2011). Inhibitory plasticity balances excitation and inhibition in sensory pathways and memory networks. Science 334, 1569–1573. doi: 10.1126/science.1211095

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinger, J. G., Greenberg, M. L., Matheu, M. P., Parker, I., Walsh, C. M., Lane, T. E., et al. (2015). Two-photon imaging of cellular dynamics in the mouse spinal cord. J. Visual. Exp. 96:52580. doi: 10.3791/52580

CrossRef Full Text | Google Scholar

Whelan, P. J. (2010). Shining light into the black box of spinal locomotor networks. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 365, 2383–2395. doi: 10.1098/rstb.2009.0322

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, J. M., Dombeck, D. A., Diaz-Rios, M., Harris-Warrick, R. M., and Brownstone, R. M. (2007). Two-photon calcium imaging of network activity in XFP-Expressing eurons in the mouse. J. Neurophysiol. 97, 3118–3125. doi: 10.1152/jn.01207.2006

PubMed Abstract | CrossRef Full Text

Wohrer, A., Humphries, M. D., and Machens, C. K. (2013). Population-wide distributions of neural activity during perceptual decision-making. Prog. Neurobiol. 103, 156–193. doi: 10.1016/j.pneurobio.2012.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuste, R. (2015). From the neuron doctrine to neural networks. Nat. Rev. Neurosci. 16, 487–497. doi: 10.1038/nrn3962

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: ensemble, spinal cord, population, CPG, motor, lognormal, fluctuation-driven, homeostasis

Citation: Berg RW (2017) Neuronal Population Activity in Spinal Motor Circuits: Greater Than the Sum of Its Parts. Front. Neural Circuits 11:103. doi: 10.3389/fncir.2017.00103

Received: 01 October 2017; Accepted: 29 November 2017;
Published: 19 December 2017.

Edited by:

Raúl E. Russo, Instituto de Investigaciones Biológicas Clemente Estable (IIBCE), Uruguay

Reviewed by:

Stephen M. Johnson, University of Wisconsin-Madison, United States
Lorenzo Cangiano, University of Pisa, Italy

Copyright © 2017 Berg. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Rune W. Berg, runeb@sund.ku.dk

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