Neuronal Population Activity in Spinal Motor Circuits: Greater Than the Sum of Its Parts
Rune W. Berg- 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. 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 ((Vm − Vthres)/σ) 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 ((Vm − Vthres)/σ) 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. 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. 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. 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
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.
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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.
Goulding, M. (2009). Circuits controlling vertebrate locomotion: moving in a new direction. Nat. Rev. Neurosci. 10, 507–518. doi: 10.1038/nrn2608
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
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
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
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.
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
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
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
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.
Humphries, M. D. (2016). The Goldilocks zone in neural circuits. eLife 5:e22735. doi: 10.7554/eLife.22735
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
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.
Kiehn, O. (2016). Decoding the organization of spinal circuits that control locomotion. Nat. Rev. Neurosci. 17, 224–238. doi: 10.1038/nrn.2016.9
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
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
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
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
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
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
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
Liu, Y.-Y., Slotine, J.-J., and Barabási, A.-L. (2011). Controllability of complex networks. Nature 473, 167–173. doi: 10.1038/nature10011
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
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
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.
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
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
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
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
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
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
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
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
Pozo, K., and Goda, Y. (2010). Unraveling mechanisms of homeostatic synaptic plasticity. Neuron 66, 337–351. doi: 10.1016/j.neuron.2010.04.028
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
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
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
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
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
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
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.
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
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
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
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
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.
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
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
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.
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
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
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
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
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
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
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
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), UruguayReviewed by:
Stephen M. Johnson, University of Wisconsin-Madison, United StatesLorenzo 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