Mathematical relationships between spinal motoneuron properties
Figures

Detailed example for the process adopted to successively create the two final and datasets (right-side thick-solid contour rectangular boxes).
These final datasets were obtained from respectively three and three normalized global datasets of experimental data obtained from the literature (dashed-contour grey-filled boxes) and . The and datasets were first transformed (⊗ symbol) into two intermediary datasets (dotted-contour boxes) by converting the and values to equivalent values with two ‘inverse’ and power relationships (oval boxes with triple dots), which had been previously obtained from two unshown steps that had yielded the final and datasets. The two intermediary datasets were merged with the remaining global dataset to yield the final dataset, to which a power relationship of the form was fitted. If and p<0.01, an ‘inverse’ power relationship (oval box) was further fitted to this final dataset. In a similar approach, the three normalized global datasets , and were transformed with the three ‘inverse’ relationships into intermediary datasets, which were merged to yield the final dataset. An ‘inverse’ power relationship was further derived to be used in the creation of the final and datasets in the next taken steps.

Experimental (A) and unknown (B) relations between motoneuron (MN) and muscle unit (mU) properties.
(A) Bubble diagram representing the pairs of MN and/or mU properties that could be investigated in this study from the results provided by the 40 studies identified in our web search. MN and mU properties are represented by circle and square bubbles, respectively. Relationships between MN properties are represented by coloured connecting lines; the colours red, blue, green, yellow, and purple are consistent with the order in which the pairs were investigated (see Table 3 for mathematical relationships). Relationships between one MN and one mU property are represented by black dashed lines. (B) Bubble diagram representing the mathematical relationships proposed in this study between pairs of MN properties for which no concurrent experimental data has been measured to date.

Normalized global datasets.
These were obtained from the 19 studies reporting cat data that measured and investigated the 17 pairs of motoneuron (MN) properties reported in Figure 2A. For each pair, the property is read on the y-axis and on the x-axis. For information, power trendlines (red dotted curves) are fitted to the data of each dataset and reported in Table 3. The 95% confidence interval of the regression is also displayed for each dataset (green dotted lines). The studies are identified with the following symbols: • (Gustafsson, 1979; Gustafsson and Pinter, 1984a; Gustafsson and Pinter, 1984b), ○ (Munson et al., 1986), ▲ (Zengel et al., 1985), ∆ (Foehring et al., 1987), ■ (Cullheim, 1978), □ (Burke, 1968; Burke and ten Bruggencate, 1971; Burke et al., 1982), ◆(Krawitz et al., 2001), ◇ (Fleshman et al., 1981), + (Eccles et al., 1958b), ☓ (Kernell, 1966; Kernell and Zwaagstra, 1981; Kernell and Monster, 1981), - (Zwaagstra and Kernell, 1980), — (Sasaki, 1991), ✶ (Pinter and Vanden Noven, 1989). The axes are given in % of the maximum retrieved values in the studies consistently with ‘Methods’ section.
-
Figure 3—source data 1
All the cat datasets presented in Figure 3.
- https://cdn.elifesciences.org/articles/76489/elife-76489-fig3-data1-v2.xlsx

Normalized motoneuron (MN) size-related final datasets.
These were obtained from the 19 studies reporting cat data that concurrently measured at least two of the morphometric and electrophysiological properties listed in Table 1. For each pair, the property is read on the y-axis and on the x-axis. The power trendlines (red dotted curves) are fitted to each dataset and are reported in Table 3. The 95% confidence interval of the regression is also displayed for each dataset (green dotted lines). For each plot, the constitutive sub-datasets that were obtained from different global datasets are specified with the following symbols identifying the property : , , □, , and .

Normalized size-dependent behaviour of the motoneuron (MN) properties , and .
For displaying purposes, the MN properties are plotted in arbitrary units as power functions (intercept ) of : according to Table 3. The larger the MN size, the larger , , and in the order of increasing slopes, and the lower , and in the order of increasing slopes.

Fivefold cross-validation of the normalized mathematical relationships.
Here are reported for each dataset the average values across the five permutations of (A) the normalized maximum error (), (B) the normalized root mean square error (), and (C) coefficient of determination (, grey bars), which is compared with the coefficient of determination (, black bars) of the power trendline fitted to the log–log transformation of global experimental datasets directly.

Global datasets for rat and mouse species and predictions of the motoneuron (MN) properties with the cat mathematical relationships (Table 4).
They were obtained from the five studies reporting data on rats and the four studies presenting data on mice reported in Appendix 1—table 5 that measured the , , and pairs of MN electrophysiological properties. are given in nA, MΩ, ms, and nF, respectively. The experimental data (○ symbol) is fitted with a power trendline (red dotted curve) and compared to the predicted quantities obtained with the scaled cat relationships in Table 4 (■ symbol).
-
Figure 7—source data 1
Datasets used in the rat plot.
- https://cdn.elifesciences.org/articles/76489/elife-76489-fig7-data1-v2.xlsx
-
Figure 7—source data 2
Datasets used in the mouse plots.
- https://cdn.elifesciences.org/articles/76489/elife-76489-fig7-data2-v2.xlsx

Density histograms of the data distributions reported in the experimental studies included in the global datasets.
Depending on the typical range over which each property spans, the distributions are divided in steps of 10 or 20%. The frequency distribution is provided in percentage of the total number of reported data points in a study. Different studies are displayed with different colours in each graph. For each dataset , the frequency distributions of property and are provided in the left and right plot, respectively.

Assessment of data variability between the experimental studies that constitute the eight global datasets that include at least two experimental studies.
For each experimental study included in the global dataset , the range, mean, coefficient of variation , and the ratio of the experimental values measured in this study were computed. Then, the average (bars) and standard deviation (error bars) across experimental studies of these four metrics (range, mean, CoV, ) were calculated for each global dataset independently. For example, the first bar in the ‘ parameter’ plot represents the average across three studies (Kernell and Zwaagstra, 1981; Cullheim, 1978; Burke et al., 1982) of the range of values reported in these studies, for the global dataset . The computed standard deviations (error bars) express for each global dataset the inter-study variability of (1) the length of the identified bandwidth of the motoneuron (MN) pool (range metric), (2) the spread of values around the mean (CoV metric), (3) the skewness of the distributions, and (4) whether the distributions from different studies are centred (mean metric). A global dataset reporting narrow error bars for parameter for the four metrics infers that the experimental studies constituting this dataset measured similar distributions of property .

Distribution of the size of the experimental datasets constituting the global datasets for assessment of the variability in the input data.
The histogram is divided between global datasets (half vertical lines), grouped as final size-dependent datasets (full vertical lines). For each global dataset, the total number of data points is reported (label: ‘N’), while the size of the constitutive experimental studies is reported in decreasing order.

Assessment of data variability between the global datasets.
This figure compares the distributions of the , and properties between the normalized global datasets built in this study. For each property, the range, mean, coefficient of variation , and the ratio of its distribution in a global dataset is compared to the other global datasets it appears in. For each property, the standard deviation between global datasets of these four metrics is computed. A low value reflects the low variability of the property distribution between global datasets.
Tables
The motoneuron (MN) and muscle unit (mU) properties investigated in this study with their notations and SI base units.
is the size of the MN. As reproduced in Table 2, the MN size is adequately described by measures of the MN surface area and the soma diameter . and define the MN-specific electrical resistance properties of the MN and set the value of the MN-specific current threshold (Binder et al., 1996; Powers and Binder, 2001; Heckman and Enoka, 2012). C and (constant among the MN pool) define the capacitance properties of the MN and contribute to the definition of the MN membrane time constant τ (Gustafsson and Pinter, 1984a; Zengel et al., 1985). is the amplitude of the membrane voltage depolarization threshold relative to resting state required to elicit an action potential. is the corresponding electrical current causing a membrane depolarization of . AHP duration is defined in most studies as the duration between the action potential onset and the time at which the MN membrane potential meets the resting state after being hyperpolarized. ACV is the axonal conduction velocity of the elicited action potentials on the MN membrane. is the size of the mU. As indicated in Table 2, the mU size is adequately described by measures of (1) the sum of the cross-sectional areas (CSAs) of the fibres composing the mU , (2) the mean fibre CSA , (3) the innervation ratio IR, that is, the number of innervated fibres constituting the mU, and (4) the mU tetanic force . is the MU twitch force.
Properties | Notation | Unit | |
---|---|---|---|
MN properties | Size: Neuron surface area Soma diameter | ||
Resistance | |||
Specific resistance per unit area | |||
Capacitance | |||
Specific capacitance per unit area | |||
Time constant | |||
Rheobase (current recruitment threshold) | |||
Voltage threshold | |||
Afterhyperpolarization duration | |||
Axonal conduction velocity | |||
mU properties | Size: | ||
Total fibre cross-sectional area | |||
Mean fibre cross-sectional area | |||
Innervation ratio | |||
Tetanic force | |||
Twitch force |
Measurable indices of motoneuron (MN) and muscle unit (mU) sizes in mammals.
and are conceptual parameters which are adequately described by the measurable and linearly inter-related quantities reported in this table.
MN size (SMN) | mU size (SmU) |
---|---|
Fitted experimental data of pairs of motoneuron (MN) properties and subsequent normalized final size-related relationships.
For information, the , p-value, and the equation are reported for each fitted global dataset. The normalized MN-size dependent relationships are mathematically derived from the transformation of the global datasets and from the power trendline fitting of the final datasets (N data points) as described in ‘Methods’. The minimum and maximum values of , and defining the 95% confidence interval of the regression are also reported in parenthesis for each global and final dataset. The values reported in this table are consistent with the values obtained when directly fitting the normalized experimental datasets with power regressions (see Appendix 1—table 1).
MN property | (normalized global datasets) | (final MN-size dependent datasets) | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
Relationship | p-value | Reference studies | p-value | points | |||||||
4.0 (2.5; 6.4) | 0.7 (0.6; 0.8) | 0.58 | < 10-5 | Cullheim, 1978; Kernell and Zwaagstra, 1981; Burke et al., 1982 | 4.0 (2.5; 6.4) | 0.7 (0.6; 0.8) | 0.58 | < 10-5 | 109 | ||
6.1 · 103 (1.2 · 103; 3.2 · 104) | −1.2 (−1.6; −0.8) | 0.34 | < 10-5 | Zwaagstra and Kernell, 1980 | 2.5 · 104 (1.2 · 104; 5.0 · 104) | −1.5 (−1.7; −1.3) | 0.41 | < 10-5 | 492 | ||
1.5 · 104 (7.4 · 103; 2.9 · 104) | −1.4 (−1.5; −1.2) | 0.41 | < 10-5 | Eccles et al., 1958a; Zwaagstra and Kernell, 1980; Gustafsson and Pinter, 1984b; Foehring et al., 1987 | |||||||
1.5 · 105 (2.7 · 104; 7.9 · 105) | −2.1 (−2.5; −1.7) | 0.61 | < 10-5 | Kernell and Zwaagstra, 1981; Burke et al., 1982 | 9.6 · 105 (4.1 · 105; 2.3 · 106) | −2.4 (−2.6; −2.2) | 0.37 | < 10-5 | 745 | ||
6.3 · 105 (1.9 · 105; 2.1 · 106) | −2.3 (−2.6; −2.0) | 0.38 | < 10-5 | Kernell, 1966; Burke, 1968; Barrett and Crill, 1974; Kernell and Zwaagstra, 1981; Fleshman et al., 1981; Gustafsson and Pinter, 1984b; Sasaki, 1991 | |||||||
6.2 · 10−1 (4.1 · 10−1; 9.2 · 10−1) | 1.1 (0.9; 1.2) | 0.65 | < 10-5 | Gustafsson, 1979; Gustafsson and Pinter, 1984b; Foehring et al., 1987; Pinter and Vanden Noven, 1989; Sasaki, 1991 | |||||||
1.1 · 103 (0.8 · 103; 1.3 · 103) | −1.0 (−1.1; −0.9) | 0.37 | < 10-5 | Kernell, 1966; Fleshman et al., 1981; Gustafsson and Pinter, 1984a; Zengel et al., 1985; Munson et al., 1986; Foehring et al., 1987; Krawitz et al., 2001 | 9.0 · 10−4 (4.7 · 10−4; 1.7 · 10−3) | 2.5 (2.4; 2.7) | 0.37 | < 10-5 | 722 | ||
3.2 · 10−6 (1.3 · 10−7; 8.2 · 10−5) | 3.7 (3.0; 4.4) | 0.37 | < 10-5 | Kernell and Monster, 1981; Gustafsson and Pinter, 1984a | |||||||
2.5 · 104 (1.3 · 104; 4.8· 104) | −1.7 (−1.9; −1.6) | 0.60 | < 10-5 | Gustafsson and Pinter, 1984a | |||||||
2.4 · 102 (2.0 · 102; 3.9 · 102) | −0.4 (−0.4; −0.3) | 0.57 | < 10-5 | Gustafsson and Pinter, 1984b | 1.2 (0.7; 2.0) | 1.0 (0.9; 1.2) | 0.28 | < 10-5 | 444 | ||
2.9 · 101 (2.4 · 101; 3.5 · 101) | 0.3 (0.2; 0.3) | 0.51 | < 10-5 | Gustafsson and Pinter, 1984a | |||||||
2.8 · 102 (1.8 · 102; 4.4 · 102) | −0.4 (−0.5; −0.3) | 0.24 | < 10-5 | Gustafsson and Pinter, 1984b | |||||||
2.5 (0.7; 8.4) | 0.8 (0.5; 1.0) | 0.17 | < 10-5 | Gustafsson and Pinter, 1984b | |||||||
8.7 (7.2; 10.6) | 0.5 (0.4; 0.6) | 0.52 | < 10-5 | Burke and ten Bruggencate, 1971; Barrett and Crill, 1974; Gustafsson, 1979; Gustafsson and Pinter, 1984b; Zengel et al., 1985; Pinter and Vanden Noven, 1989; Sasaki, 1991 | 2.6 · 104 (1.5 · 104; 4.5 · 104) | −1.5 (−1.6; −1.4) | 0.46 | < 10-5 | 649 | ||
2.2 (1.3; 3.5) | 0.8 (0.7; 1.0) | 0.63 | < 10-5 | Gustafsson and Pinter, 1984b | |||||||
2.3 · 102 (1.9 · 102; 2.7 · 102) | −0.4 (−0.5; −0.3) | 0.72 | < 10-5 | Gustafsson and Pinter, 1984a | |||||||
1.2 · 104 (2.2 · 103; 6.6 · 104) | −1.3 (−1.7; −0.9) | 0.30 | < 10-5 | Gustafsson and Pinter, 1984b |
Mathematical empirical cat relationships between the motoneuron (MN) properties , and .
Each column provides the relationships between one and the eight other MN properties. If one property is known, the complete MN profile can be reconstructed by using the pertinent line in this table. All constants and properties are provided in SI base units (metres, seconds, ohms, farads, and amperes). The relationships involving were obtained from theoretical relationships involved in Rall’s cable theory (see ‘Discussion’).
Fitted experimental data of pairs of one muscle unit (mU) and one motoneuron (MN) property and subsequent relationships.
Species | (fitted relationships) | (final relationships) | ||||
---|---|---|---|---|---|---|
Relationship | p-value | Reference studies | c | |||
Rat | 3.4 | 0.45 | 6.10-3 | Kanda and Hashizume, 1992 | 2.4 | |
Cat | 9.4 | 0.43 | <10-5 | Knott et al., 1971 | 6.6 | |
7.2 | 0.37 | <10-5 | Mcphedran et al., 1965; Wuerker et al., 1965; Appelberg and Emonet-Dénand, 1967; Proske and Waite, 1974; Bagust, 1974; Jami and Petit, 1975; Stephens and Stuart, 1975; Burke et al., 1982; Emonet-Dénand et al., 1988 | 5.0 | ||
-1.3 | 0.27 | 6.10-5 | Dum and Kennedy, 1980 | 3.2 | ||
2.0 | 0.21 | 2.10-2 | Burke et al., 1982 | 2.0 | ||
Mouse | -2.1 | 0.42 | <10-5 | Manuel and Heckman, 2011; Martínez-Silva et al., 2018 | 5.1 | |
1.3 | 0.64 | 2.10-4 | Manuel and Heckman, 2011 | 3.3 | ||
1.0 | 0.80 | 6.10-2 | Manuel and Heckman, 2011 | 2.5 | ||
Mean ± sd | 3.8 ± 1.5 |
-
Table 5—source data 1
Numerical data used to derive the relationships presented in Table 5.
- https://cdn.elifesciences.org/articles/76489/elife-76489-table5-data1-v2.xlsx
values obtained for each experimental dataset when performing a linear regression analysis on directly (‘Linear’) and on the (‘Power’) and (‘Exponential’) transformations of .
The values returned by the three types of regression cannot be directly compared to estimate the best model. However, the power fit returned for relatively more experimental datasets than the linear and exponential fits.
Mathematical empirical normalized relationships between the motoneuron (MN) properties , and .
Each column provides the relationships between one and the six other MN properties. All constants and properties are normalized up to a theoretical 100% maximum value. To scale the normalized relationships for a specific mammalian species, the normalized intercept must be scaled with a pair of values for properties and obtained from the same motoneuron in that species.
Typical ranges of physiological values for and in cat rat and mouse species.
SMN is found to vary over an average -fold range, which sets the amplitude of the theoretical ranges. Absolute {min; max} reports the minimum and maximum values retrieved in the reference studies for and , while average {min; max} is obtained as the average across reference studies of minimum and maximum values retrieved per study.
Property | Unit | Absolute{min;max} | Average{min; max} | Reference studies | Theoretical range | ||
---|---|---|---|---|---|---|---|
Cat | [µm] | 2.2 | Kernell, 1966; Cullheim, 1978; Zwaagstra and Kernell, 1980; Kernell and Zwaagstra, 1981; Ulfhake and Kellerth, 1981; Zwaagstra and Kernell, 1981; Burke et al., 1982; Donselaar et al., 1986; Destombes et al., 1992 | [33; 79] | |||
[mm2] | 2.7 | Barrett and Crill, 1974; Ulfhake and Kellerth, 1981; Burke et al., 1982; Ulfhake and Kellerth, 1984; Ulfhake and Cullheim, 1988; Moschovakis et al., 1991 | [0.18; 0.44] | ||||
Rat | [µm] | 2.4 | Swett et al., 1986; Vult von Steyern et al., 1999; Copray and Kernell, 2000; Ishihara et al., 2001; Deardorff et al., 2013; Mierzejewska-Krzyżowska et al., 2014 | ||||
Mouse | [µm] | 2.5 | Vult von Steyern et al., 1999 | ||||
[mm2] | 4.4 | Amendola and Durand, 2008; Brandenburg et al., 2020 |
Typical ranges of physiological values for the motoneuron (MN) properties , and .
As described in ‘Methods’, is the average among reference studies of the ratios of minimum and maximum values; the properties experimentally vary over a -fold range. This ratio compares with the theoretical ratio (with taken from Table 3), which sets the amplitude of the theoretical ranges. Absolute and average {min; max} are obtained as described in the main sections.
Validation of the cat relationships (Table 4) against rat and mouse data.
nME, normalized maximum error; nRMSE, normalized root mean square error; , coefficient of determination between experimental and predicted quantities; , coefficient of determination of the power trendline directly fitted to the experimental data.
Animal | Dataset | Reference studies | ||||
---|---|---|---|---|---|---|
Rat | Gardiner, 1993; Bakels and Kernell, 1993; Lee and Heckman, 1998; Button et al., 2008; Turkin et al., 2010; Krutki et al., 2015 | 352 | 16 | 0.53 | 0.54 | |
Mouse | Delestrée et al., 2014; Martínez-Silva et al., 2018; Huh et al., 2021 | 385 | 16 | 0.46 | 0.46 | |
Manuel et al., 2009 | 1219 | 168 | 0.35 | 0.40 | ||
Manuel et al., 2009 | 1915 | 98 | 0.38 | 0.38 |
For one {A; B} global dataset | Training set | Test set |
---|---|---|
Permutation 1 | p1, p2, p3, p4 | p5 |
Permutation 2 | p2, p3, p4, p5 | p1 |
Permutation 3 | p1, p3, p4, p5 | p2 |
Permutation 4 | p1, p2, p4, p5 | p3 |
Permutation 5 | p1, p2, p3, p5 | p4 |
Datasets | Previous validation | New validation | ||||
---|---|---|---|---|---|---|
nME | nRMSE | R2 | nME | nRMSE | R2 | |
AHP = ka ∙ SMNa | 70 | 14 | 0,40 | 67 | 12 | 0,46 |
AHP = ka ∙ ACVa | 158 | 20 | 0,43 | 147 | 19 | 0,43 |
R = ka ∙ SMNa | 255 | 20 | 0,56 | 312 | 31 | 0,51 |
R = ka ∙ ACVa | 246 | 23 | 0,43 | 205 | 18 | 0,46 |
R = ka ∙ AHPa | 153 | 21 | 0,63 | 143 | 21 | 0,66 |
Ith = ka ∙ Ra | 410 | 19 | 0,37 | 300 | 19 | 0,36 |
Ith = ka ∙ ACVa | 171 | 26 | 0,34 | 155 | 20 | 0,42 |
Ith = ka ∙ AHPa | 84 | 15 | 0,59 | 78 | 17 | 0,61 |
C = ka ∙ Ra | 61 | 12 | 0,58 | 55 | 13 | 0,57 |
C = ka ∙ ITaℎ | 47 | 17 | 0,53 | 52 | 19 | 0,49 |
C = ka ∙ AHPa | 75 | 16 | 0,26 | 62 | 16 | 0,22 |
C = ka ∙ ACVa | 82 | 21 | 0,16 | 74 | 21 | 0,20 |
τ = ka ∙ Ra | 86 | 13 | 0,54 | 73 | 13 | 0,51 |
τ = ka ∙ AHPa | 88 | 14 | 0,62 | 77 | 13 | 0,64 |
τ = ka ∙ ITaℎ | 83 | 15 | 0,68 | 65 | 13 | 0,74 |
τ = ka ∙ ACVa | 81 | 17 | 0,38 | 84 | 16 | 0,44 |