Altered Motoneuron Properties Contribute to Motor Deficits in a Rabbit Hypoxia-Ischemia Model of Cerebral Palsy

Cerebral palsy (CP) is caused by a variety of factors attributed to early brain damage, resulting in permanently impaired motor control, marked by weakness and muscle stiffness. To find out if altered physiology of spinal motoneurons (MNs) could contribute to movement deficits, we performed whole-cell patch-clamp in neonatal rabbit spinal cord slices after developmental injury at 79% gestation. After preterm hypoxia-ischemia (HI), rabbits are born with motor deficits consistent with a spastic phenotype including hypertonia and hyperreflexia. There is a range in severity, thus kits are classified as severely affected, mildly affected, or unaffected based on modified Ashworth scores and other behavioral tests. At postnatal day (P)0–5, we recorded electrophysiological parameters of 40 MNs in transverse spinal cord slices using whole-cell patch-clamp. We found significant differences between groups (severe, mild, unaffected and sham control MNs). Severe HI MNs showed more sustained firing patterns, depolarized resting membrane potential, and fired action potentials at a higher frequency. These properties could contribute to muscle stiffness, a hallmark of spastic CP. Interestingly altered persistent inward currents (PICs) and morphology in severe HI MNs would dampen excitability (depolarized PIC onset and increased dendritic length). In summary, changes we observed in spinal MN physiology likely contribute to the severity of the phenotype, and therapeutic strategies for CP could target the excitability of spinal MNs.

Introduction

Cerebral palsy (CP) is not well understood, despite its prevalence and seriousness. There exist only a few evidence-based treatments for CP: the effectiveness of many currently used therapeutic strategies is unclear (Novak et al., 2013; Wimalasundera and Stevenson, 2016). Recent clinical advances include use of magnesium sulfate and hypothermia after hypoxic-ischemic encephalopathy to acutely reduce neural damage (Yager et al., 1993; Thoresen et al., 1996; Magee et al., 2011; Rouse and Gibbins, 2013), but little basic research is devoted to addressing symptoms after they arise. Part of the problem in treating CP may be the diversity of causes including neonatal stroke, placental insufficiency, preterm birth, inflammation, traumatic injury, difficulties during birth and many other contributing factors (MacLennan and International Cerebral Palsy Task Force, 1999; Graham et al., 2016). Another problem could be that modeling the condition in animals is complicated, and while rodent models are useful for the development of neuroprotective strategies, larger animal models are needed to study motor deficits (Clowry et al., 2014; Cavarsan et al., 2019).

Loss of corticospinal control of movement is considered causative of motor deficits in CP, but little investigation into the precise effect on spinal circuits has been conducted. A notable exception is the work of John H. Martin and colleagues, who have documented changes in corticospinal synaptic connectivity in specific spinal laminae and loss of cholinergic interneurons after either cortical silencing or lesioning the corticospinal tract (Martin et al., 1999; Li and Martin, 2000; Friel and Martin, 2005, 2007; Friel et al., 2012; Jiang et al., 2016, 2018). Another important study showed changes in parvalbumin-positive spinal interneurons after cortical silencing in development (Clowry et al., 2004; Clowry, 2007). Both of these interneuron classes (parvalbumin-positive and cholinergic) are synaptically connected to spinal MNs and could contribute to altered motor output after developmental injury. Based on these foundational studies, we hypothesized that altering development with hypoxia-ischemia (HI) injury would also alter the development of MNs, specifically the electrophysiological properties governing excitability in spinal MNs. We further hypothesized that changes in excitability would correspond/contribute to the severity of motor deficits. In short, that altered activity of spinal MNs could contribute to muscle stiffness and spasticity.

To assess changes in the intrinsic properties of spinal MNs, we used the rabbit HI model of CP (Derrick et al., 2004). It’s been shown in previous studies that HI injury during late gestation in rabbits can result in a variety of neurologic and muscular damage, including muscle stiffness (Derrick et al., 2004), loss of neurons in cortical layers 3 and 5, white matter injury, thinning of the corticospinal tract (Buser et al., 2010), cell death in the spinal cord and decreased numbers of spinal MNs (Drobyshevsky and Quinlan, 2017), increased sarcomere length, decreased muscle mass and hyperreflexia (Synowiec et al., 2019). There is also an increase in spinal monoamines which could increase the excitability of spinal neurons and thus promote spasticity (Bellot et al., 2014; Drobyshevsky et al., 2015). Thus, changes observed in spinal MNs in the rabbit model could be directly compared to motor deficits.

Changes in MN physiology are likely to contribute to motor impairment in CP, yet this has not been directly assessed in any animal models. Thus, we assessed electrophysiological parameters in spinal MNs in neonatal rabbits after sham surgery or hypoxic-ischemic insult during development.

Materials and Methods

All rabbits were used according to National Institutes of Health guide for the care and use of laboratory animals, and the University of Rhode Island’s, Northwestern University’s and Northshore University Health System’s Animal Care and Use Committee guidelines. Pregnant New Zealand White rabbits (Charles River Laboratories, Inc. Wilmington, MA, USA), underwent HI procedures as described previously (Derrick et al., 2004; Buser et al., 2010). Briefly, the procedure was performed at ~80% gestation (day 25 of gestation, or E25), a time when HI has been found to result in the greatest degree of white matter injury and corticospinal tract thinning. Dams were anesthetized, and the left femoral artery was isolated. A Fogarty balloon catheter was inserted into the femoral and advanced to the level of the descending aorta, above the uterine arteries and inflated for 40 min. Sham animals underwent the same procedures but without the inflation of the catheter. After the procedure, the dam recovered and later gave birth to kits with HI injuries. Categorization of the severity of the phenotype was performed by a blinded observer, using a modified Ashworth scale, observation/tests for activity, locomotion, posture, righting reflex, muscle tone (as described in Derrick et al., 2004). Kits could be given a maximum score of 6 (normal posture, righting and joint resistance). Since, there was a large variation in the severity of motor deficits, HI kits were divided into three groups: HI unaffected (scores were the same range as control kits, 5–6), HI mild (scores 3–4), HI severe (scores 1–2). One rabbit kit which was affected by HI but displayed a phenotype of hypotonia instead of hypertonia was removed from the data set. All other kits included in this study displayed a hypertonic phenotype if affected by HI.

Patch Clamp

Whole-cell patch-clamp was performed similarly to previously published work (Quinlan et al., 2011) from P0–5. Briefly, horizontal spinal cord slices 350 μm thick were obtained using a Leica 1,000 vibratome. Slices were incubated for 1 h at 30°C and recordings were performed at room temperature. During recording, slices were perfused with oxygenated (95% O₂ and 5% CO₂) modified Ringer’s solution containing (in mM): 111 NaCl, 3.09 KCl, 25.0 NaHCO₃, 1.10 KH₂PO₄, 1.26 MgSO₄, 2.52 CaCl₂, and 11.1 glucose at 2 ml/min. Whole-cell patch electrodes (1–3 MΩ) contained (in mM) 138 K-gluconate, 10 HEPES, 5 ATP-Mg, 0.3 GTP-Li and Texas Red dextran (150 μM, 3,000 MW). PICs were measured in voltage-clamp mode with holding-potential of −90 mV and depolarizing voltage ramps of both 36 mV/s and 11.25 mV/s bringing the cell to 0 mV in 2.5 s or 8 s, respectively and then back to the holding potential in the following 2.5 s or 8 s. Input resistance was measured from the slope of the leak current near the holding potential. Capacitance was measured with Multiclamp’s whole-cell capacitance compensation function. Resting membrane potential was measured in voltage-clamp as the voltage at which there is 0 pA of injected current in the descending ramp. In current clamp, frequency—current measurements were obtained from current ramps. The first spike on the current ramp was used to measure the properties of action potentials. The threshold voltage was defined as the voltage at which the action potential slope exceeds 10 V/s. Rate of rise and fall of the action potential were measured as peak and trough of the first derivative of the action potential. The duration of the action potential was measured at half-peak (defined as the midpoint between overshoot and threshold voltages). Depolarizing current steps of varying amplitude were used to find maximum firing rates (near depolarization block) and to measure after-spike afterhyperpolarization (in single spikes elicited near-threshold). Hyperpolarizing current steps (typically between −850 and −1,250 pA) were used to measure hyperpolarization-activated sag currents (I_H). Neuron selection: neurons were targeted in MN pools mainly from cervical and lumbar regions of the cord and were removed from the data set if their resting membrane potential was more depolarized than −45 mV in current clamp.

Imaging

After electrophysiological measurements were obtained, MNs were imaged to assess anatomical development, and photos were obtained of the electrode placement within the spinal cord slice, as shown in Figure 1. Images were acquired with a Nikon microscope fitted with a 40× water-dipping objective lens and two-photon excitation fluorescence microscopy performed with a galvanometer-based Coherent Chameleon Ultra II laser. To optimize the excitation of red/green fluorophores, the laser was tuned to 900 nm. 3D reconstructions of MNs were created using Neurolucida 360° software. It is likely that some processes extended past the surface of the slice and were excluded from reconstructions. However, since this was the case for all MNs in this study it is unlikely to have an impact on the findings.

FIGURE 1

Figure 1. Patch clamp of spinal motoneurons (MNs) with dye-filling via patch electrode as shown in panel (A). Scale bar = 100 μm. (B) Placement of patch electrode (at arrow) within the slice is captured with a photo. (C) Map of recorded MNs within medial and lateral motor pools. Red circles = HI Severe, blue triangles = HI Mild, yellow stars = HI unaffected, gray rectangles = Sham Controls.

Statistics

All variables were checked for normality and homogeneity (using Shapiro–Wilk and Levene’s tests). The variables that were parametric (normal and homogenous) were run with one-way Analysis of Variance (ANOVA) and then assessed post hoc with a Tukey test for between-group significance. The non-parametric variables were run with the Kruskal–Wallis test followed by Dunn’s test as a post hoc analysis to assess significance between groups, adjusting the p-value for multiple comparisons. The analysis was done using R software for determining the significance of parameters over groups, according to their injury classification (sham, HI unaffected, HI mildly affected, and HI severely affected). Injury classification, age of the kit (P0–P5), and spinal cord region (cervical, thoracic, lumbar or sacral) were all tested. Significance was determined by p-values ≤ 0.05.

Results

After HI surgery was performed in pregnant dams at 79% gestation, kits were born naturally about a week later. At ages P0, neonates were rated as unaffected, mildly affected or severely affected. Since, there was a large variation in the severity of motor deficits, HI kits were divided into three groups: HI unaffected, HI mild, HI severe. Experiments were all performed in the first 5 days of life. Over 40 spinal MNs were patched in transverse spinal cord slices, and over 40 parameters were measured from each. To determine the significance of the variables, a one-way ANOVA was performed to find differences among three factors: (1) injury classification (sham control, HI unaffected, HI mildly affected, and HI severely affected); (2) age (postnatal day 0–5); and (3) spinal region (cervical through sacral). All data, including mean, standard deviation, group size, and the p-value is included in table format (Tables 1–5 and Supplementary Tables S1, S2).

TABLE 1

Table 1. Frequency—current characteristics.

TABLE 2

Table 2. Persistent inward current characteristics (5 s ramp).

TABLE 3

Table 3. Persistent inward current characteristics (16 s ramp).

TABLE 4

Table 4. Soma morphology characteristics.

TABLE 5

Table 5. Dendrite characteristics.

HI MNs Show Sustained Firing and Higher Firing Frequency

In rabbit kits severely injured by HI, MNs had significantly increased sustained firing. The frequency current (F-I) relationship was measured using the current ramps, as shown in Figure 2. Depolarizing current ramps are used to evoke firing, and current at onset and offset of firing (I_ON and I_OFF) determine ΔI. In sham control MNs, ΔI was larger and always a positive value, indicating firing ceased at a higher current amplitude on the descending ramp than the current level that elicited firing on the ascending ramp (see Figure 2A). Severe HI MNs had a smaller, and usually negative ΔI, revealing increasingly sustained firing (see Figure 2B). Also, resting membrane potential was significantly more depolarized in HI MNs than sham controls (see Figure 2F). Another significant difference in severe HI MNs is the instantaneous firing rate, as shown in Figure 3. At the start of a depolarizing current step, the peak (instantaneous) firing rate is higher in severe HI MNs than sham controls (see Figure 3B). Sustained firing is also apparent in Figure 3C, in the second current step which evokes a brief burst of action potentials followed by depolarization block in both MNs. The severe HI MN recovers and resumes firing while the sham control MN remains in depolarization block. Mean values for both instantaneous and steady-state parameters are shown in bar graphs (Figures 3D,E). Significant changes in MN physiology were present in severely affected animals: post hoc analysis showed significant changes between sham and severe HI MNs in ΔI, RMP, and instantaneous firing rate. The instantaneous firing rate also reached significance in mild HI MNs vs. sham. No significant changes were present in HI unaffected MNs in any properties. While these parameters (ΔI, RMP, and instantaneous firing rate) suggest increased excitability in HI-injured MNs, there was no significant change in the threshold, I_ON, or I_OFF in HI MNs. Thus, MNs from severely affected kits should not be classified as hyperexcitable per se, since they begin firing with the same depolarizing input and at the same voltage threshold.

FIGURE 2

Figure 2. Severe HI MNs show more sustained firing than sham control MNs, and a more depolarized resting membrane potential. Control MNs (A) have larger values for ΔI compared to severe HI (B). Average ΔI and threshold are shown in panels (C) and (D) for all groups. The frequency current relationship (E) was not significantly different between groups. Resting membrane potential was significantly more depolarized in HI severe MNs than sham control MNs (F). Error bars = SEM. Scale bars in (A) = 20 mV (vertical) and 0.5 s (horizontal) and applies to panels (A,B).

FIGURE 3

Figure 3. Instantaneous firing frequency is increased in severe HI MNs. Depolarizing current steps (A) evoked action potentials (C) in a typical sham control (black trace) and severe HI (red trace) MNs. (B) Dotted lines show peak firing frequency of MNs in panel (C). (D) Average instantaneous and (E) steady-state firing frequency are shown for all groups. Error bars = SEM. Scale bar in (C) = 0.5 s and applies to panels (A–C).

Changes in Spike Properties and Subthreshold Responses Were Not Present

A complete analysis of I_H (sag and rebound currents), action potential parameters, and after-spike after-hyperpolarization (AHP) was performed and no significant differences in these parameters were found between groups. All data is included in Supplementary Tables S1, S2.

Persistent Inward Currents Suggest Excitability Is Dampened After HI

PICs were significantly affected by hypoxia-ischemia (HI), revealing that intrinsic excitability may be dampened. PICs were measured using both short (5 s) and long (16 s) protocols, which can preferentially activate and inactivate Na⁺ and Ca²⁺ mediated PICs. The different protocols yielded different results. For example, as shown in Figure 4, voltage dependence (PIC onset and PIC Max) was unchanged in the PICs evoked using a short 5-s voltage ramp HI severe MNs. Using longer voltage ramps (16 s), PIC onset was significantly depolarized in HI severe MNs compared to sham (see Tables 2, 3). Since the change in PIC onset was more pronounced in longer ramps this could suggest an altered balance of Na⁺ and Ca²⁺ channel activation or altered activation/inactivation of these channels (see “Discussion” section). Change in the magnitude of the PIC was not observed outright in either of the protocols: the magnitude of the currents was similar between groups (see Figure 4 and Tables 2, 3). Intrinsic properties including capacitance (which significantly increased in HI MNs compared to sham) and input resistance of all MNs are included in Table 2.

FIGURE 4

Figure 4. PICs are altered in HI injured MNs. (A) The typical current response to a 5-s voltage ramp in a sham (black trace) and HI severe (red) MN. (B) Leak-subtracted PICs from panel (A) are similar in amplitude. There is a trend for depolarized PIC onset (C) in short ramps after HI injury, which reaches significance in long ramps (E). PIC Max did not reach significance (D) and (F). Error bars = SEM.

Morphology Affected by HI Injury

Morphology of MNs was assessed in all patched neurons, as shown in Figure 5. As suggested by the significantly larger whole-cell capacitance, there were changes in MN morphology after HI injury. The soma size was unchanged: there were no significant differences between groups in soma’s largest cross-sectional area (Figure 5D) or other measurements of soma size (Table 4). There was, however, a significant increase in dendrite length in HI injured MNs compared to sham controls (Figure 5E), which could account for changes in electrical properties. Since, we recorded from motor pools throughout the spinal cord (cervical through sacral), there was a large amount of variability within our data set. Future studies in our lab focus on the analysis of specific motor pools. All data on dendritic morphology is included in Table 5.

FIGURE 5

Figure 5. Morphology is affected by HI. Typical sham control (A) and HI severe (B) MNs filled with dye during patch-clamp (electrodes visible on right). Average values of whole-cell capacitance (C), soma largest cross-sectional area (D), total dendrite length (E), and number of stem dendrites (F) are included for all neurons. Scale bar in (A) = 100 μm applies to (A,B). Error bars = SEM.

Age and Spinal Region

No significant changes in MN properties were found due to the spinal region. Postnatal age had a significant effect on the following properties: input resistance (decrease with age), action potential size (mV), rate of rise and rate of fall (all increase with age), 5 s PIC amplitude (increase with age), and normalized PIC (PIC/Cap; current density increased with age for both 5 and 16 s ramps). These results are in line with previous studies on embryonic and postnatal maturation of MNs.

Discussion

Summary

Electrophysiological properties of spinal MNs are altered by developmental HI injury, and the magnitude of changes is correlated to the severity of motor deficits. Specifically, these changes include increased sustained firing and a higher firing rate. These changes could indicate increased excitability and would contribute to muscle stiffness that is common in spastic CP. However, concomitant changes in PIC onset and longer dendritic length could serve to dampen excitability and could contribute to weakness. Since traditional views of CP largely view motor dysfunction as a result of the damaged motor cortex improperly signaling to spinal neurons, our new evidence suggests this is only part of the problem. Spinal MNs are not developing the same after HI and show an overall change in intrinsic properties. Whether these changes are directly due to the HI insult or indirectly due to downstream effects must be determined by future work.

Contribution of Spinal Motoneurons to Dysfunction in Cerebral Palsy

Here, we show that MNs show altered excitability after HI injury, including elevated resting potential, and more sustained firing. Previous work showed that after HI injury in rabbits there were also fewer spinal MNs, and spinal interneurons in lamina VII were undergoing apoptosis (Drobyshevsky and Quinlan, 2017). After the loss of corticospinal projections, it was recently found that the spinal cholinergic interneurons which give rise to C boutons on MNs are lost (Jiang et al., 2016, 2018, 2019). Taken together this data suggests spinal circuits are: (1) just as vulnerable to HI injury as the developing cortex; and (2) potentially functioning with fewer neurons and altered circuitry. In addition to fewer neurons, there is also atrophy of the muscles which could contribute to motor deficits in CP. In mice, rabbits, and humans, muscle atrophy appears along with losses in numbers of MNs (Han et al., 2013; Marciniak et al., 2015; Drobyshevsky and Quinlan, 2017; Brandenburg et al., 2018). Recent work has shown similar changes to muscle architecture in the rabbit HI model of CP to humans, including atrophy, muscle shortening, and longer sarcomere length. Increased muscle stiffness in rabbits affected by HI was found even after administration of anesthetic—indicating some muscle stiffness is derived from mechanical changes in the muscles, though a large component of the muscle stiffness was diminished with anesthetic thus was driven neurally (Synowiec et al., 2019). During development, both feedback and feed-forward signaling can regulate growth and maturation, processes that may be disrupted in CP in both MNs and muscle fibers. It is likely the loss of spinal interneurons, MNs and muscle fibers reduce coordination and strength in those with CP. Altered size and excitability of MNs is also a feature of other motor disorders including amyotrophic lateral sclerosis and spinal muscular atrophy (Quinlan et al., 2011, 2019; Gogliotti et al., 2012; Shoenfeld et al., 2014; Dukkipati et al., 2018). Since our data were collected from multiple MN pools throughout the spinal cord, from control and HI injured rabbit kits, relatively large variability in parameters was expected (Kanning et al., 2010), however consistent differences in electrical properties emerged in this data set. This suggests consistent changes are present across motor pools after injury and support further exploration of interventions for CP and other motor disorders that target spinal MNs. These treatments could include neuromodulators and therapies aimed at restoring balance between excitation and inhibition within spinal circuits for alleviation of spasticity.

Neuromodulation

The exact causes of the changes in MN physiology observed here are unclear, but they could result from the increase in spinal monoamines that occurs after developmental HI injury in both rodents and rabbits (Bellot et al., 2014; Drobyshevsky et al., 2015). Serotonin is generally thought of as a neurotransmitter and neuromodulator, but developmental disruption in 5HT is associated with neurological disorders including autism, Rett syndrome, Down’s syndrome and, more recently, CP (Bar-Peled et al., 1991; Whittle et al., 2007; Bellot et al., 2014; Yang et al., 2014; De Filippis et al., 2015; Drobyshevsky et al., 2015; Muller et al., 2016; Wirth et al., 2017). Serotonin increases MN excitability in neonatal and juvenile mice, rats and guinea pigs (Wang and Dun, 1990; Ziskind-Conhaim et al., 1993; Hsiao et al., 1997, 1998), and likely has the same effect on rabbit MNs. Depolarization of the resting membrane potential, increased action potential firing through hyperpolarization of the voltage threshold and enhanced PIC, increased action potential height and reduction of high-voltage activated Ca²⁺ entry are all associated with 5HT receptor activation in neonatal and adult MNs (Elliot and Wallis, 1992; Larkman and Kelly, 1992; Lindsay and Feldman, 1993; Bayliss et al., 1995; Hsiao et al., 1997, 1998; Inoue et al., 1999; Gilmore and Fedirchuk, 2004; Li et al., 2006). Therefore increased 5HT could have a direct impact on excitability, though in HI rabbits the increase in 5HT was accompanied by decreased mRNA for 5HT₂ receptors and increased mRNA for the SERT serotonin transporter (Drobyshevsky et al., 2015). In light of that finding, it is not clear that neurons remain responsive to 5HT. In the experiments here, all MNs were recorded in spinal cord slices incubated and perfused in standard oxygenated aCSF without any serotonergic drugs present. Therefore, HI MNs in vivo could show different levels of excitability since they would be in the presence of elevated 5HT, while our MNs were all recorded in the same standard physiological solution. Thus, the contribution of 5HT to the altered excitability observed here is restricted to its chronic effects on neuron development, namely morphological changes. Serotonin 5HT_1A and 5HT_2A receptor activation increases neurite outgrowth, dendritic branching, and spine formation (Bou-Flores et al., 2000; Fricker et al., 2005; Mogha et al., 2012), findings that align well with the present finding of increased dendritic length in the HI MNs. Future experiments will be needed to address the role of 5HT in enhancing MN excitability and its effects on synaptically-evoked action potentials. Synaptic events in dendrites would more strongly evoke PICs, though both altered dendritic morphology and elevated 5HT could dampen them.

Possible Mechanism of Altered MN Output

The mechanism for increased MN activity and thus muscle stiffness may be due to delayed Na⁺ channel inactivation. In neonatal MNs, Na⁺ channels generate the majority of the PIC and account for action potential initiation/repetitive firing. Specifically, Nav 1.1, 1.2 and 1.6 type Na⁺ channels (Boiko et al., 2001, 2003; Rush et al., 2005) inactivate faster than Ca²⁺ channels that contribute to PICs (Perrier and Hounsgaard, 2003; Li et al., 2006). Short voltage ramps preferentially measure the Na⁺ PIC for this reason: Na⁺ channels inactivate quickly enough that even on the descending ramp of the short protocol, there is no longer a region of negative slope (see Figure 4). Changes in Na⁺ channel inactivation in adult MNs along with postnatal development of the longer-lasting Ca²⁺ PIC makes typical adult MNs display more negative ΔI values and longer-lasting PICs (Harvey et al., 2006; Li et al., 2006; Quinlan et al., 2011). In the neonatal rabbit MNs, sham controls showed positive ΔI values that are quite typical for neonates, while HI MNs showed significantly more negative values. This could be due to slower Na⁺ channel inactivation or increased contribution of Ca²⁺ channels to the PICs after injury. Interestingly the more depolarized resting potential found in HI severe MN would serve to increase Na⁺ channel inactivation. To fully determine the altered mechanism of aberrant firing after HI injury future studies into the biophysical properties of Na⁺ channels and maturation of Ca²⁺ channel expression must be pursued.

Severity in Motor Deficits and Electrophysiology

Generally, unaffected and mildly affected MNs showed parameters that were intermediate between control MNs and severe MNs. There was only one category in which mildly affected MNs become statistically significantly different from sham controls (instantaneous firing frequency). Thus, electrophysiological changes were overwhelmingly in line with phenotype, suggesting aberrant MN properties contribute to the severity of the phenotype. It cannot be ruled out, however, that HI “unaffected” MNs may have a subtle phenotype that is not readily evident based on testing we performed here. Or perhaps abnormalities in these rabbits would develop in later in life: in humans CP patients, diagnosis of CP is not made until 18–24 months of life and the peak of spasticity occurs around 4 years of age (Hägglund and Wagner, 2008; Hadders-Algra, 2014; Novak et al., 2017). In the rabbit model, deficits have not been characterized past P18, and a detailed analysis of the progression of motor deficits from P0–18 is lacking. Future work is needed to assess the maturation of the MN properties in different groups, the potential contribution of delayed Na⁺ channel inactivation in CP, the progression of motor deficits with age, and the development of new therapeutic strategies that could target MNs.

Postnatal Maturation and Injury

The present study only extended from postnatal day 0–5, but even within this narrow window, significant changes were apparent in MN electrical parameters. As MNs undergo postnatal maturation, they grow larger, with more complex dendritic arborizations, and gain the ability to fire action potentials at higher rates, as reviewed in Carrascal et al. (2005). We found that age had a significant effect on 5 s PIC amplitude (increasing with age), normalized PIC amplitude (both 5 and 16 s PIC/Cap increased with age), input resistance (decreasing with age), and action potential size, rate of rise and rate of fall (all increase with age). While neuron size increases during this period, the amplitude of the PIC typically increases more than proportionally (Quinlan et al., 2011). In other words, voltage-gated ion channels are being inserted into the membrane at a faster rate than the cell is increasing in size, resulting in an increased normalized PIC amplitude with age. We suspect this parallels the ability of MNs to fire action potentials at higher rates during postnatal development, and the acquisition of coordinated motor control and weight-bearing in developing animals.

Conclusion

Changes in MN physiology after developmental injury are consistent with motor deficits in rabbits. This suggests not only brain injuries but also changes in the spinal cord contribute to impaired function in CP. Exploring both altered maturation of spinal neurons and loss of descending connectivity should be pursued to improve outcomes for individuals with CP.

This manuscript has been released as a Pre-Print at BioRxiv (Steele et al., 2019).

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Ethics Statement

The animal study was reviewed and approved by University of Rhode Island’s, Northwestern University’s and Northshore University Health System’s Animal Care and Use Committees.

Author Contributions

PS, CC, MG, and KQ wrote the manuscript. PS, CC, LD, MW, and KQ performed the experiments and analyzed the data. AD, MG, and KQ planned the experiments and interpreted the results. PS, CC, and NK performed statistical analysis. All authors approved the content and assisted in revision and review of manuscript.

Funding

This project was supported by National Institutes of Health, NINDS R01NS091278 to AD and NINDS R01NS104436 to MG and KQ.

Conflict of Interest

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

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2020.00069/full#supplementary-material.

References

Bar-Peled, O., Gross-Isseroff, R., Ben-Hur, H., Hoskins, I., Groner, Y., and Biegon, A. (1991). Fetal human brain exhibits a prenatal peak in the density of serotonin 5-HT1A receptors. Neurosci. Lett. 127, 173–176. doi: 10.1016/0304-3940(91)90787-t

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayliss, D. A., Umemiya, M., and Berger, A. J. (1995). Inhibition of N- and P-type calcium currents and the after-hyperpolarization in rat motoneurones by serotonin. J. Physiol. 485, 635–647. doi: 10.1113/jphysiol.1995.sp020758

PubMed Abstract | CrossRef Full Text | Google Scholar

Bellot, B., Peyronnet-Roux, J., Gire, C., Simeoni, U., Vinay, L., and Viemari, J. C. (2014). Deficits of brainstem and spinal cord functions after neonatal hypoxia-ischemia in mice. Pediatr. Res. 75, 723–730. doi: 10.1038/pr.2014.42

PubMed Abstract | CrossRef Full Text | Google Scholar

Boiko, T., Rasband, M. N., Rock, S., Caldwell, J. H., Mandel, G., Trimmer, J. S., et al. (2001). Compact myelin dictates the differential targeting of two sodium channel isoforms in the same axon. Neuron 30, 91–104. doi: 10.1016/s0896-6273(01)00265-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Boiko, T., Van Wart, A., Caldwell, J. H., Levinson, S. R., Trimmer, J. S., and Matthews, G. (2003). Functional specialization of the axon initial segment by isoform-specific sodium channel targeting. J. Neurosci. 23, 2306–2313. doi: 10.1523/jneurosci.23-06-02306.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bou-Flores, C., Lajard, A., Monteau, R., De Maeyer, E., Seif, I., Lanoir, J., et al. (2000). Abnormal phrenic motoneuron activity and morphology in neonatal monoamine oxidase a-deficient transgenic mice: possible role of a serotonin excess. J. Neurosci. 20, 4646–4656. doi: 10.1523/jneurosci.20-12-04646.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandenburg, J. E., Gransee, H. M., Fogarty, M. J., and Sieck, G. C. (2018). Differences in lumbar motor neuron pruning in an animal model of early onset spasticity. J. Neurophysiol. 120, 601–609. doi: 10.1152/jn.00186.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Buser, J. R., Segovia, K. N., Dean, J. M., Nelson, K., Beardsley, D., Gong, X., et al. (2010). Timing of appearance of late oligodendrocyte progenitors coincides with enhanced susceptibility of preterm rabbit cerebral white matter to hypoxia-ischemia. J. Cereb. Blood Flow Metab. 30, 1053–1065. doi: 10.1038/jcbfm.2009.286

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrascal, L., Nieto-gonzalez, L., Cameron, W. E., Torres, B., and Nunez-abades, P. A. (2005). Changes during the postnatal development in physiological and anatomical characteristics of rat motoneurons studied in vitro. Brain Res. Rev. 49, 377–387. doi: 10.1016/j.brainresrev.2005.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavarsan, C. F., Gorassini, M. A., and Quinlan, K. A. (2019). Animal models of developmental motor disorders: parallels to human motor dysfunction in cerebral palsy. J. Neurophysiol. 122, 1238–1253. doi: 10.1152/jn.00233.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Clowry, G. J. (2007). The dependence of spinal cord development on corticospinal input and its significance in understanding and treating spastic cerebral palsy. Neurosci. Biobehav. Rev. 31, 1114–1124. doi: 10.1016/j.neubiorev.2007.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Clowry, G. J., Basuodan, R., and Chan, F. (2014). What are the best animal models for testing early intervention in cerebral palsy? Front. Neurol. 5, 1–17. doi: 10.3389/fneur.2014.00258

PubMed Abstract | CrossRef Full Text | Google Scholar

Clowry, G. J., Davies, B. M., Upile, N. S., Gibson, C. L., and Bradley, P. M. (2004). Spinal cord plasticity in response to unilateral inhibition of the rat motor cortex during development: changes to gene expression, muscle afferents and the ipsilateral corticospinal projection. Eur. J. Neurosci. 20, 2555–2566. doi: 10.1111/j.1460-9568.2004.03713.x

PubMed Abstract | CrossRef Full Text | Google Scholar

De Filippis, B., Chiodi, V., Adriani, W., Lacivita, E., Mallozzi, C., Leopoldo, M., et al. (2015). Long-lasting beneficial effects of central serotonin receptor 7 stimulation in female mice modeling rett syndrome. Front. Behav. Neurosci. 9, 1–11. doi: 10.3389/fnbeh.2015.00086

PubMed Abstract | CrossRef Full Text | Google Scholar

Derrick, M., Luo, N. L., Bregman, J. C., Jilling, T., Ji, X., Fisher, K., et al. (2004). Preterm fetal hypoxia-ischemia causes hypertonia and motor deficits in the neonatal rabbit: a model for human cerebral palsy? J. Neurosci. 24, 24–34. doi: 10.1523/jneurosci.2816-03.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Drobyshevsky, A., and Quinlan, K. A. (2017). Spinal cord injury in hypertonic newborns after antenatal hypoxia-ischemia in a rabbit model of cerebral palsy. Exp. Neurol. 293, 13–26. doi: 10.1016/j.expneurol.2017.03.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Drobyshevsky, A., Takada, S. H., Luo, K., Derrick, M., Yu, L., Quinlan, K. A., et al. (2015). Elevated spinal monoamine neurotransmitters after antenatal hypoxia-ischemia in rabbit cerebral palsy model. J. Neurochem. 132, 394–402. doi: 10.1111/jnc.12997

PubMed Abstract | CrossRef Full Text | Google Scholar

Dukkipati, S. S., Garrett, T. L., and Elbasiouny, S. M. (2018). The vulnerability of spinal motoneurons and soma size plasticity in a mouse model of amyotrophic lateral sclerosis. J. Physiol. 596, 1723–1745. doi: 10.1113/jp275498

PubMed Abstract | CrossRef Full Text | Google Scholar

Elliot, P., and Wallis, D. I. (1992). Serotonin and l-norepinephrine as mediators of altered excitability in neonatal rat motoneurons studied in vitro. Neuroscience 47, 533–544. doi: 10.1016/0306-4522(92)90163-v

PubMed Abstract | CrossRef Full Text | Google Scholar

Fricker, A. D., Rios, C., Devi, L. A., and Gomes, I. (2005). Serotonin receptor activation leads to neurite outgrowth and neuronal survival. Mol. Brain Res. 138, 228–235. doi: 10.1016/j.molbrainres.2005.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Friel, K., Chakrabarty, S., Kuo, H. C., and Martin, J. (2012). Using motor behavior during an early critical period to restore skilled limb movement after damage to the corticospinal system during development. J. Neurosci. 32, 9265–9276. doi: 10.1523/jneurosci.1198-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Friel, K. M., and Martin, J. H. (2005). Role of sensory-motor cortex activity in postnatal development of corticospinal axon terminals in the cat. J. Comp. Neurol. 485, 43–56. doi: 10.1002/cne.20483

PubMed Abstract | CrossRef Full Text | Google Scholar

Friel, K. M., and Martin, J. H. (2007). Bilateral activity-dependent interactions in the developing corticospinal system. J. Neurosci. 27, 11083–11090. doi: 10.1523/jneurosci.2814-07.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilmore, J., and Fedirchuk, B. (2004). The excitability of lumbar motoneurones in the neonatal rat is increased by a hyperpolarization of their voltage threshold for activation by descending serotonergic fibres. J. Physiol. Lond. 558, 213–224. doi: 10.1113/jphysiol.2004.064717

PubMed Abstract | CrossRef Full Text | Google Scholar

Gogliotti, R. G., Quinlan, K. A., Barlow, C. B., Heier, C. R., Heckman, C. J., and DiDonato, C. J. (2012). Motor neuron rescue in spinal muscular atrophy mice demonstrates that sensory-motor defects are a consequence, not a cause, of motor neuron dysfunction. J. Neurosci. 32, 3818–3829. doi: 10.1523/jneurosci.5775-11.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, H. K., Rosenbaum, P., Paneth, N., Dan, B., Lin, J.-P., and Damiano, D. L. (2016). Cerebral palsy. Nat. Rev. Dis. Prim. 2:15082. doi: 10.1038/nrdp.2015.82

CrossRef Full Text | Google Scholar

Hadders-Algra, M. (2014). Early diagnosis and early intervention in cerebral palsy. Front. Neurol. 5:185. doi: 10.3389/fneur.2014.00185

PubMed Abstract | CrossRef Full Text | Google Scholar

Hägglund, G., and Wagner, P. (2008). Development of spasticity with age in a total population of children with cerebral palsy. BMC Musculoskelet. Disord. 9:150. doi: 10.1186/1471-2474-9-150

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Q., Feng, J., Qu, Y., Ding, Y., Wang, M., So, K.-F., et al. (2013). Spinal cord maturation and locomotion in mice with an isolated cortex. Neuroscience 253, 235–244. doi: 10.1016/j.neuroscience.2013.08.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Harvey, P. J., Li, X., Li, Y., and Bennett, D. J. (2006). 5-HT 2 receptor activation facilitates a persistent sodium current and repetitive firing in spinal motoneurons of rats with and without chronic spinal cord injury. J. Neurophysiol. 96, 1158–1170. doi: 10.1152/jn.01088.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsiao, C. F., Del Negro, C. A., Trueblood, P. R., and Chandler, S. H. (1998). Ionic basis for serotonin-induced bistable membrane properties in guinea pig trigeminal motoneurons. J. Neurophysiol. 79, 2847–2855. doi: 10.1152/jn.1998.79.6.2847

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsiao, C. F., Trueblood, P. R., Levine, M. S., and Chandler, S. H. (1997). Multiple effects of serotonin on membrane properties of trigeminal motoneurons in vitro. J. Neurophysiol. 77, 2910–2924. doi: 10.1152/jn.1997.77.6.2910

PubMed Abstract | CrossRef Full Text | Google Scholar

Inoue, T., Itoh, S., Kobayashi, M., Kang, Y., Matsuo, R., Wakisaka, S., et al. (1999). Serotonergic modulation of the hyperpolarizing spike afterpotential in rat jaw-closing motoneurons by PKA and PKC. J. Neurophysiol. 82, 626–637. doi: 10.1152/jn.1999.82.2.626

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y.-Q., Armada, K., and Martin, J. H. (2019). Neuronal activity and microglial activation support corticospinal tract and proprioceptive afferent sprouting in spinal circuits after a corticospinal system lesion. Exp. Neurol. 321:113015. doi: 10.1016/j.expneurol.2019.113015

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y. Q., Sarkar, A., Amer, A., and Martin, J. H. (2018). Transneuronal downregulation of the premotor cholinergic system after corticospinal tract loss. J. Neurosci. 38, 8329–8344. doi: 10.1523/jneurosci.3410-17.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y.-Q., Zaaimi, B., and Martin, J. H. (2016). Competition with primary sensory afferents drives remodeling of corticospinal axons in mature spinal motor circuits. J. Neurosci. 36, 193–203. doi: 10.1523/jneurosci.3441-15.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanning, K. C., Kaplan, A., and Henderson, C. E. (2010). Motor neuron diversity in development and disease. Annu. Rev. Neurosci. 33, 409–440. doi: 10.1146/annurev.neuro.051508.135722

PubMed Abstract | CrossRef Full Text | Google Scholar

Larkman, P. M., and Kelly, J. S. (1992). Ionic mechanisms mediating 5-hydroxytryptamine- and noradrenaline-evoked depolarization of adult rat facial motoneurones. J. Physiol. 456, 473–490. doi: 10.1113/jphysiol.1992.sp019347

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., and Martin, J. H. (2000). Postnatal development of differential projections from the caudal and rostral motor cortex subregions. Exp. Brain Res. 134, 187–198. doi: 10.1007/s002210000454

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Murray, K., Harvey, P. J., Ballou, E. W., and Bennett, D. J. (2006). Serotonin facilitates a persistent calcium current in motoneurons of rats with and without chronic spinal cord injury. J. Neurophysiol. 97, 1236–1246. doi: 10.1152/jn.00995.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindsay, A. D., and Feldman, J. L. (1993). Modulation of respiratory activity of neonatal rat phrenic motoneurones by serotonin. J. Physiol. Lond. 461, 213–233. doi: 10.1113/jphysiol.1993.sp019510

PubMed Abstract | CrossRef Full Text | Google Scholar

MacLennan, A., and International Cerebral Palsy Task Force (1999). A template for defining a causal relation between acute intrapartum events and cerebral palsy: international consensus statement. BMJ 319, 1054–1059. doi: 10.1136/bmj.319.7216.1054

PubMed Abstract | CrossRef Full Text | Google Scholar

Magee, L., Sawchuck, D., Synnes, A., van Dadelszen, P., Magnesium Sulphate For Fetal Neuroprotection Consensus Committee, and Maternal Fetal Medicine Committee. (2011). SOGC clinical practice guideline. Magnesium sulphate for fetal neuroprotection. J. Obstet. Gynaecol. Canada 33, 516–529. doi: 10.1016/S1701-2163(16)34886-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Marciniak, C., Li, X., and Zhou, P. (2015). An examination of motor unit number index in adults with cerebral palsy. J. Electromyogr. Kinesiol. 25, 444–450. doi: 10.1016/j.jelekin.2015.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, J. H., Kably, B., and Hacking, A. (1999). Activity-dependent development of cortical axon terminations in the spinal cord and brain stem. Exp. Brain Res. 125, 184–199. doi: 10.1007/s002210050673

PubMed Abstract | CrossRef Full Text | Google Scholar

Mogha, A., Guariglia, S. R., Debata, P. R., Wen, G. Y., and Banerjee, P. (2012). Serotonin 1A receptor-mediated signaling through ERK and PKCα a is essential for normal synaptogenesis in neonatal mouse hippocampus. Transl. Psychiatry 2:e66. doi: 10.1038/tp.2011.58

PubMed Abstract | CrossRef Full Text | Google Scholar

Muller, C. L., Anacker, A. M. S. J., and Veenstra-VanderWeele, J. (2016). The serotonin system in autism spectrum disorder: from biomarker to animal models. Neuroscience 321, 24–41. doi: 10.1016/j.neuroscience.2015.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Novak, I., Mcintyre, S., Morgan, C., Campbell, L., Dark, L., Morton, N., et al. (2013). A systematic review of interventions for children with cerebral palsy: state of the evidence. Dev. Med. Child Neurol. 55, 885–910. doi: 10.1111/dmcn.12246

PubMed Abstract | CrossRef Full Text | Google Scholar

Novak, I., Morgan, C., Adde, L., Blackman, J., Boyd, R. N., Brunstrom-Hernandez, J., et al. (2017). Early, accurate diagnosis and early intervention in cerebral palsy: advances in diagnosis and treatment. JAMA Pediatr. 171, 897–907. doi: 10.1001/jamapediatrics.2017.1689

PubMed Abstract | CrossRef Full Text | Google Scholar

Perrier, J. F., and Hounsgaard, J. (2003). 5-HT2 receptors promote plateau potentials in turtle spinal motoneurons by facilitating an L-type calcium current. J. Neurophysiol. 89, 954–959. doi: 10.1152/jn.00753.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinlan, K. A., Reedich, E., Arnold, W. D., Puritz, A., Cavarsan, C. F., Heckman, C., et al. (2019). Hyperexcitability precedes motoneuron loss in the Smn 2B/- mouse model of spinal muscular atrophy. J. Neurophysiol. 122, 1297–1311. doi: 10.1152/jn.00652.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinlan, K. A., Schuster, J. E., Fu, R., Siddique, T., and Heckman, C. J. (2011). Altered postnatal maturation of electrical properties in spinal motoneurons in a mouse model of amyotrophic lateral sclerosis. J. Physiol. 589, 2245–2260. doi: 10.1113/jphysiol.2010.200659

PubMed Abstract | CrossRef Full Text | Google Scholar

Rouse, D. J., and Gibbins, K. J. (2013). Magnesium sulfate for cerebral palsy prevention. Semin. Perinatol. 37, 414–416. doi: 10.1053/j.semperi.2013.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Rush, A. M., Dib-Hajj, S. D., and Waxman, S. G. (2005). Electrophysiological properties of two axonal sodium channels, Nav1.2 and Nav1.6, expressed in mouse spinal sensory neurones. J. Physiol. Lond. 564, 803–815. doi: 10.1113/jphysiol.2005.083089

PubMed Abstract | CrossRef Full Text | Google Scholar

Shoenfeld, L., Westenbroek, R. E., Fisher, E., Quinlan, K. A., Tysseling, V. M., Powers, R. K., et al. (2014). Soma size and Cav1.3 channel expression in vulnerable and resistant motoneuron populations of the SOD1G93A mouse model of ALS. Physiol. Rep. 2:e12113. doi: 10.14814/phy2.12113

PubMed Abstract | CrossRef Full Text | Google Scholar

Steele, P. R., Cavarsan, C. F., Dowaliby, L., Westefeld, M., Drobyshevsky, A., Gorassini, M. A., et al. (2019). Altered motoneuron properties contribute to motor deficits in a rabbit hypoxia ischemia model of cerebral palsy. bioRxiv [Preprint]. doi: 10.1101/817957

CrossRef Full Text | Google Scholar

Synowiec, S., Lu, J., Yu, L., Goussakov, I., Lieber, R., and Drobyshevsky, A. (2019). Spinal hyper-excitability and altered muscle structure contribute to muscle hypertonia in newborns after antenatal hypoxia-ischemia in a rabbit cerebral palsy model. Front. Neurol. 9, 1–18. doi: 10.3389/fneur.2018.01183

PubMed Abstract | CrossRef Full Text | Google Scholar

Thoresen, M., Bågenholm, R., Løberg, E. M., Apricena, F., and Kjellmer, I. (1996). Posthypoxic cooling of neonatal rats provides protection against brain injury. Arch. Dis. Child. 74, F3–F9. doi: 10.1136/fn.74.1.f3

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., and Dun, N. (1990). 5-Hyroxytryptamine responses in neonate rat motoneurones in vitro. J. Physiol. 430, 87–103. doi: 10.1113/jphysiol.1990.sp018283

PubMed Abstract | CrossRef Full Text | Google Scholar

Whittle, N., Sartori, S. B., Dierssen, M., Lubec, G., and Singewald, N. (2007). Fetal down syndrome brains exhibit aberrant levels of neurotransmitters critical for normal brain development. Pediatrics 120, e1465–e1471. doi: 10.1542/peds.2006-3448

PubMed Abstract | CrossRef Full Text | Google Scholar

Wimalasundera, N., and Stevenson, V. L. (2016). Cerebral palsy. Pract. Neurol. 16, 184–194. doi: 10.1136/practneurol-2015-001184

PubMed Abstract | CrossRef Full Text | Google Scholar

Wirth, A., Holst, K., and Ponimaskin, E. (2017). How serotonin receptors regulate morphogenic signalling in neurons. Prog. Neurobiol. 151, 35–56. doi: 10.1016/j.pneurobio.2016.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Yager, J., Towfighi, J., and Vannucci, R. C. (1993). Influence of mild hypothermia on hypoxic-ischemic brain damage in the immature rat. Pediatr. Res. 34, 525–529. doi: 10.1203/00006450-199310000-00029

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C. J., Tan, H. P., and Du, Y. J. (2014). The developmental disruptions of serotonin signaling may involved in autism during early brain development. Neuroscience 267, 1–10. doi: 10.1016/j.neuroscience.2014.02.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Ziskind-Conhaim, L., Seebach, B. S., and Gao, B. X. (1993). Changes in serotonin-induced potentials during spinal cord development. J. Neurophysiol. 69, 1338–1349. doi: 10.1152/jn.1993.69.4.1338

PubMed Abstract | CrossRef Full Text | Google Scholar