Spinal lumbar dI2 interneurons contribute to stability of bipedal stepping
Abstract
Peripheral and intraspinal feedback is required to shape and update the output of spinal networks that execute motor behavior. We report that lumbar dI2 spinal interneurons in chicks receive synaptic input from afferents and premotor neurons. These interneurons innervate contralateral premotor networks in the lumbar and brachial spinal cord, and their ascending projections innervate the cerebellum. These findings suggest that dI2 neurons function as interneurons in local lumbar circuits, are involved in lumbo-brachial coupling, and that part of them deliver peripheral and intraspinal feedback to the cerebellum. Silencing of dI2 neurons leads to destabilized stepping in posthatching day 8 hatchlings, with occasional collapses, variable step profiles, and a wide-base walking gait, suggesting that dI2 neurons may contribute to the stabilization of the bipedal gait.
Introduction
The spinal cord integrates and relays the somatosensory inputs required for further execution of complex motor behaviors. Interneurons (INs) that differentiate at the ventral progenitor domain, V3-V0, are involved in the control of rhythmic motor activity, alternating between the left and right limbs, as well as between the flexor and extensor muscles (Lai et al., 2016; Osseward and Pfaff, 2019; Alaynick et al., 2011). Some of the dorsally born INs, dI1, dI3, and dI6, migrate ventrally and are also assembled within circuitries that control motor activity (Yuengert et al., 2015; Bui et al., 2013; Andersson et al., 2012), while other dorsal progenitor neurons, dI4 and dI5, give rise to INs that mediate somatosensation (Lai et al., 2016).
dI2 neurons originate in the dorsal spinal cord. The progenitor pdI2 cells, topographically positioned between the adjacent dorsally located dI1 and ventrally located dI3 neurons, express Ngn1, Ngn2, Olig3, and Pax3 transcription factors (TFs). Early postmitotic dI2 neurons undergo ventral migration and are defined by the combinatorial expression of Foxd3+/Lhx1+/Pou4f1+ TFs (Alaynick et al., 2011; Morikawa et al., 2009; Francius et al., 2013). Importantly, none of these TFs is specific to dI2 neurons; rather, their combinatorial expression defines dI2. The lack of dI2-specific cell fate markers and the dynamic expression of the above TFs in other INs causes ambiguity regarding the molecular profile and outcome of late postmitotic dI2 neurons. Using intersectional genetics, we have shown previously that dI2 neurons are commissurally projecting neurons (Avraham et al., 2009). The lack of dI2-specific TFs impeded their genetic targeting. Hence, little is known about the wiring and physiological function of dI2 neurons.
The maintenance of stability and the coordination, precision, and timing of movements are regulated and modulated by the cerebellum. Anatomical and electrophysiological studies of cats and rodents revealed two major pathways ascending from neurons in the lumbar spinal cord to the cerebellum: the dorsal spinocerebellar tract (DSCT) and ventral spinocerebellar tract (VSCT). DSCT neurons are considered to relay mainly proprioceptive information, while VSCT neurons are thought to relay internal spinal network information to the cerebellum along with proprioceptive data (Jankowska and Hammar, 2013; Spanne and Jörntell, 2013; Stecina et al., 2013; Jiang et al., 2015). While subpopulations of DSCT neurons are genetically accessible (Hantman and Jessell, 2010), the genetic inaccessibility of VSCT neurons hinders efforts to reveal their actual contribution to the regulatory functions of the cerebellum in locomotion and other motor behaviors.
In the present study, we investigated the possible functions of dI2 neurons in chick motor behavior. Employing intersectional genetics in the chick spinal cord, we targeted dI2 neurons and found evidence implicating them in the control of stability during locomotion. There are several advantages to performing these studies in chicks. The patterning of neurons within the spinal cord (Jessell, 2000) and the spinocerebellar tracts (Furue et al., 2010; Furue et al., 2011; Uehara et al., 2012) is conserved between mammals and birds. In addition, chicks use bipedal locomotion (evolved in humans and birds) that can be examined soon after hatching. To decode the circuitry and function of spinal INs, we developed a unique circuit-deciphering toolbox that enables neuron-specific targeting and tracing of circuits in the chick embryo (Hadas et al., 2014), and we utilized kinematic analysis of overground bipedal stepping by the hatched chicks following silencing of dI2 neurons.
Our studies revealed that lumbar dl2 neurons receive synaptic inputs from inhibitory and excitatory premotor neurons (pre-MNs) and relay output to the cerebellar granular layer, pre-MNs in the contralateral spinal cord and the contralateral dI2 cells. In a kinematic analysis of overground stepping by posthatching day (P) 8 hatchlings after inhibition of dI2 neuronal activity by expression of the tetanus toxin light chain gene, the genetically manipulated hatchlings showed an unstable gait, demonstrating that dI2 neurons play a role in shaping and stabilizing the bipedal gait.
Results
To define potential VSCT neurons within spinal INs, we defined the following criteria: (1) soma location in accordance with precerebellar neurons at the lumbar level, which were previously revealed by retrograde labeling experiments of the chick cerebellar lobes (Furue et al., 2010; Furue et al., 2011; Uehara et al., 2012), (2) commissural neurons, (3) excitatory neurons, and (4) non-pre-MNs (Lai et al., 2016; Osseward and Pfaff, 2019; Alaynick et al., 2011). Based on these criteria, dI1c and dI2 neurons are likely candidates (Bermingham et al., 2001; Yuengert et al., 2015; Figure 1—figure supplement 1). This is further supported by Sakai et al., 2012, who demonstrated that in the embryonic day (E) 12 chick, the dI1 and dI2 axons project to the hindbrain and toward the cerebellum. In the current study, we focused on deciphering the circuitry and function of dI2 neurons and their possible association with VSCT.
dI2 INs are mainly excitatory neurons with commissural axonal projections
To label dI2 neurons, axons, and terminals in the chick spinal cord, we used intersection between enhancers of two TFs expressed by dI2 – Ngn1 and Foxd3 – via the expression of two recombinases (Cre and FLPo) and double conditional reporters (Figure 1—figure supplement 1A). We have shown previously that the combination of these enhancers reliably labels dI2 neurons (Avraham et al., 2009; Hadas et al., 2014). The recombinases and the double conditional reporter plasmids were delivered via spatially restricted electroporation to the lumbar spinal cord at HH18. At E5, early postmitotic dI2 neurons migrate ventrally from the dorsolateral to the midlateral spinal cord (Figure 1A). As they migrate ventrally, at E6, dI2 neurons assume a midlateral position along the dorsoventral axis (Figure 1B). Subsequently, dI2 neurons migrate medially, and at E17, comparable to postnatal day 4 (P4) of mice, most of them (70.7 and 71.5% at the sciatic and the crural levels, respectively) occupy lamina VII (Figure 1C and F). At all rostrocaudal levels and embryonic stages, dI2 axons cross the floor plate (Figure 1A–D). After crossing, dI2 axons extend rostrally for a few segments in the ventral funiculus (VF) and subsequently turn into the lateral funiculus (LF) (Avraham et al., 2009; Figure 1C and D). Collaterals originating from the crossed VF and LF tracts invade the contralateral spinal cord (Figure 1C and D; white arrows).
A recent study suggested that dI2 neurons at early stages of development in mice (E9.5–E13.5, comparable to chick E4–8) can be divided into several subclasses based on their genetic signature and degree of maturation (Delile et al., 2019). To assess the diversity of dI2 neurons in the chick, the expression of dI2 TFs in dI2::GFP cells was analyzed at E5 before and during ventral migration, at E6 and E14. The early postmitotic dI2::GFP cells at E5 were a homogenous population defined by Foxd3+/Lhx1+/Pou4f1+/Pax2- (Figure 1—figure supplement 2A, B, E). dI2 neurons that underwent ventral migration at E5, as well as at E6 and E14, express variable combinations of Lhx1, Pou4f1, and FoxP1/2/4 (Figure 1—figure supplement 2C–E). At E14, approximately 50% of dI2::GFP cells did not express any of the tested TFs (Figure 1—figure supplement 2E), suggesting that the early expression of TFs is required for cell fate acquisition, axon guidance, and target recognition, while their expression is not required after the establishment of the circuitry, as shown for other spinal INs (Bikoff et al., 2016). Interestingly, approximately 12% of ventrally migrating dI2 neurons (from E5 to E14) expressed Pax2 (Figure 1—figure supplement 2D and E). Pax2 is associated with a GABAergic inhibitory phenotype (Cheng et al., 2004), suggesting that a subpopulation of dI2 are inhibitory neurons. The distribution of excitatory and inhibitory dI2 neurons is also apparent at E17. In situ hybridization on cross-sections of the E17 dI2::GFP-labeled lumbar spinal cord using the vGlut2 probe revealed that 73% were vGlut2+, while the VIAAT probe, which labels GABAergic and glycinergic inhibitory neurons, measured 27% VIAAT+ dI2 neurons (Figure 1G–I; N = 308 neurons from two embryos). Similar percentages of Gad2- and Slc6a5-dI2-expressing cells were also found in mice (Delile et al., 2019). At E13–17 at the caudal lumbar level and at the level of the sciatic plexus, most dI2 neurons are located in the medial aspect of lamina VII. Approximately 91% of dI2 neurons are small-diameter neurons located in the lateral dorsal aspect of lamina VII, and 9% are large-diameter neurons. At the lumbar sciatic plexus level, large-diameter dI2 neurons are located mostly at the ventral aspect of lamina VII (Figure 1F) and at the level of the crural plexus in the ventral and dorsal aspect of lamina VII (Figure 1—figure supplement 3A). The location of the lumbar dI2 neurons, mainly lamina VII, is also apparent in posthatching chicks (P8, Figure 1—figure supplement 3C and D). Importantly, large-diameter dI2 neurons were apparent only at the lumbar level (Figure 1F, Figure 1—figure supplement 3A). The division of large- and small-diameter lumbar dI2 neurons was not reflected in the expression of the tested TFs or in a specific neurotransmitter phenotype; the inhibitory/excitatory ratios were 0.297 ± 0.13 and 0.343 ± 0.02, for large- and small-diameter dI2 neurons, respectively. Hence, dI2 neurons consist of several subpopulations, as has been shown in other spinal INs (Bikoff et al., 2016; Delile et al., 2019; Sweeney et al., 2018).
Subpopulation of dI2 neurons project to the cerebellum
To study the supraspinal targets of dI2 neurons, axonal and synaptic reporters were expressed in lumbar dI2 neurons (Figure 2A). At stage HH18 (E3), dI2 enhancers were co-electroporated with the double conditional axonal reporter membrane-tethered Cherry and the synaptic reporter SV2-GFP (Figure 1—figure supplement 1A). Expression in the lumbar spinal cord was attained by using thin electrodes positioned near the lumbar segments. At E17, the stage in which the internal granule layer is formed in the chick cerebellum, the axons and synapses of dI2 neurons were studied. dI2 axons cross the spinal cord at the floor plate at the segmental level, ascend to the cerebellum, enter through the superior cerebellar peduncle, and cross back to the ipsilateral side of the cerebellum (Figure 2B). Synaptic boutons were noticeable in the granule layer at the ipsilateral and contralateral sides of the anterior cerebellar lobules (Figure 2C). Synaptic boutons were also present in the central cerebellar nuclei (Figure 2—figure supplement 1A).
The difference in soma size between dorsally and ventrally located dI2 neurons prompted us to test which dI2 neurons project to the cerebellum. dI2 neurons and precerebellar neurons were colabeled by genetic targeting of dI2 at early stages of embryogenesis (HH18), coupled with intracerebellar injection of replication-defective HSV-LacZ at E15 or PRV-Cherry (Figure 2A). Cholera toxin subunit B (CTB) was coinjected with PRV-Cherry to verify primary infection of precerebellar neurons. Spinal neurons retrogradely labeled from the cerebellum consist of double-crossed VSCT neurons and ipsilaterally projecting DSCT neurons. However, dI2 neurons, double labeled by genetic targeting and retrograde labeling from the cerebellum, are all VSCT neurons since dI2 neurons are commissural neurons. The soma distributions of precerebellar neurons, dI2 neurons, and dI2 synapses overlapped at the sciatic level (Figure 2E, G and H) and to a lesser extent at the crural level (Figure 2—figure supplement 1B–D). We found that large-diameter dI2 neurons were mostly colabeled, most of them in the ventral aspect of lamina VII (Figure 2D). Interestingly, many of the cherry+ and LacZ+ neurons were contacted by dI2 axons (Figure 2F), suggesting that dI2 neurons innervate precerebellar neurons.
Segmental crossing at the lumbar level and recrossing back to the ipsilateral side at the cerebellum are characteristic of the VSCT projection pattern. To measure the proportion of precerebellar neurons in dI2, we used sparse labeling of lumbar dI2 neurons coupled with whole-mount imaging of the E13 spinal cord from the sacral level to the cerebellum utilizing light-sheet microscopy and iDISCO imaging (Belle et al., 2014; Renier et al., 2014; Figure 3; Video 1, Video 2, Video 3, Video 4). 85 neurons were labeled at the lumbar level, 17 of which were large-diameter neurons. The axons of all dI2 neurons cross the midline. Longitudinally projecting axons were apparent at the contralateral VF and LF (Figure 3A–F; Video 1, Video 2, Video 3). In the gray matter, axonal ramifications originating from collateral branching were apparent along the entire extent of the spinal cord (Figure 3A–C; Video 1, Video 2, Video 3). Bifurcation of the longitudinal axons was also apparent. The bifurcating branch elongated rostrally for a few segments and subsequently turned transversely into the spinal cord (Figure 3—figure supplement 1A). We counted the numbers of longitudinally projecting axons at different levels (Figure 3A and D-H): 22 axons at LS1 (Figure 3D), 30 at T1 (Figure 3E), 21 at C9 (Figure 3F), 8 at the rostral brain stem (Figure 3G), and 7 entering the cerebellum via the superior cerebellar peduncle (Figure 3H). The elevated number of axons in T1 levels likely reflects longitudinal bifurcations. Considering that the spinal cord was observed at E13, a relatively early stage of development, it is likely that additional axons enter the cerebellum at later stages of development.
To evaluate the complexity of the branching pattern, we reconstructed the axonal projection and branching pattern of two large-diameter dI2 neurons at the lumbar and brachial levels (Figure 3B and C). Numerous collaterals that penetrate the spinal cord along its entire length were evident. Importantly, VSCT dI2 neurons projected to spinal targets at the lumbar, thoracic, and brachial spinal levels and to the brain stem and cerebellum (Figure 3).
To measure the proportion of dI2 in VSCT neurons, we labeled VSCT axons with GFP and dI2 axons with Cherry (for experimental design, see Figure 3—figure supplement 1B and C). The number of axons expressing the reporters at the contralateral superior cerebellar peduncle was scored. 10% of the VSCT axons belonged to dI2 neurons (Figure 3—figure supplement 1D). Thus, the large-diameter dI2 neurons constitute 10% of the VSCT neurons, consistent with the anatomical observation that the VSCT comprises a heterogeneous population of INs (Jankowska and Hammar, 2013; Stecina et al., 2013).
Mapping the synaptic input and output of dI2 neurons
To obtain the connectome of dI2 neurons, we employed enhancer-mediated synaptic labeling of presynaptic neurons coupled with soma labeling of postsynaptic neurons. We used three criteria for assessing synaptic contact: (1) the likelihood of connectivity was examined by spatial overlap of axonal terminals from the presumed presynaptic neurons and the somata of the postsynaptic neurons; (2) synaptic boutons were detected on the somatodendritic membrane of postsynaptic neuron; and (3) colabeling was observed between the presynaptic reporter and synaptotagmin (syn). We used confocal imaging and 3D reconstitution to score overlap (Figure 4—figure supplement 1A, B).
As a proof of concept, we tested the colabeling of dI2::SV2-GFP and syn in dI2 to contralateral pre-MN synapses. Of 144 genetically labeled boutons, 121 (84%) were syn+. The syn- boutons were significantly smaller (Figure 4—figure supplement 1C). We set a volume threshold (0.07 µm3, Figure 4—figure supplement 1C), and small-volume SV2-GFP boutons were not considered synapses in the study. Using these criteria, we mapped the putative pre-dI2 and post-dI2 neurons.
dI2 neurons receive synaptic input from pre-MNs and sensory neurons
To assess the synaptic input to dI2 neurons, we investigated their synaptic connectivity with the following: (1) dorsal root ganglion (DRG) neurons (Figure 4A). (2) Ipsilateral pre-MNs. General ipsilateral pre-MNs were labeled by injecting a PRV-Cherry virus into the ipsilateral hindlimb musculature (Hadas et al., 2014; Figure 4B). Two genetically defined classes of pre-MNs were examined: dI1i excitatory INs (Figure 4C, Figure 4—figure supplement 2) and the V1 inhibitory pre-MN population (Bikoff et al., 2016; Gosgnach et al., 2006; Figure 4D). (3) Reticulospinal tract neurons (Figure 4E). dI2, DRG, V1, and dI1 neurons were labeled using specific enhancers (Figure 1—figure supplement 1A).
A density profile of the axons of DRG neurons (Figure 4A, Figure 4—figure supplement 3A–D) was aligned with the density plots of the dl2 somata. Overlap between the axonal terminals of DRG neurons was evident (Figure 4A, Figure 4—figure supplement 3C and D). Contact between DRG axons and dI2 neurons was mainly apparent in the dorsal dI2 neurons, while the ventral dI2 neurons received little to no input from DRG neurons (2.8 ± 2.4 vs. 16.9 ± 11.3 contacts per neuron for ventral and dorsal dI2 neurons, respectively; p<1e-5; Figure 4A, Figure 4—figure supplement 3E). In contrast, large- and small-diameter dI2 neurons did not exhibit a significant difference in DRG axon contacts (10.4 ± 14.9 vs. 7.8 ± 5.7 contacts per neuron for large and small dI2 neurons, respectively; p=0.4) (Figure 4—figure supplement 3F).
The density plot of PRV-labeled pre-MNs overlapped with the density plots of the dl2 somata (Figure 4B), and the axonal terminals of pre-MNs were visible on dI2 somata and dendrites (Figure 4B′,B′′), suggesting that pre-MNs contacted dI2 neurons. To solidify the evidence for pre-MN/dI2 connectivity, we used synaptic reporters expressed in genetically identified dI1 and V1 pre-MNs. Excitatory dI1 synapses and inhibitory V1 synapses overlapped with the density plots of the dl2 somata (Figure 4C and D). Synaptic connections, evaluated by boutons found on dI2 dendrites and somata, were apparent from V1 and dI1i (Figure 4B–D, Figure 4—figure supplement 3G).
Serotonergic neurons are the main reticulospinal input to VSCT in cat (Hammar et al., 2004; Hammar and Maxwell, 2002). Serotonergic synapses were concentrated on motor neurons and were not observed on dI2 neurons (Figure 4E, Figure 4—figure supplement 3H). Double labeling of 5-HT and dI2 neurons did not reveal any synaptic input. The lack of synaptic serotonergic input may be related to the difference in species or may suggest that other, non-dI2 VSCT neurons located adjacent to motoneurons are contacted by the reticulospinal neurons. The analysis of synaptic inputs supports the concept that dI2 neurons constitute part of the VSCT. These cells receive input from sensory afferents and inhibitory and excitatory pre-MNs and project to the cerebellum.
dI2 neurons innervate contralateral lumbar and brachial pre-MNs and dI2 neurons
Axon collaterals of dI2 invade the gray matter along the entire length of the spinal cord, as revealed by whole-mount staining of spinal cords electroporated with an alkaline phosphatase reporter (dI2::AP) (Figure 5A), cross-sections of dI2 neurons expressing membrane-tethered EGFP (Figure 5B), and light-sheet microscopy analysis (Figure 3—figure supplement 1; Video 1, Video 2, Video 3). The region innervated by dI2 collaterals (arrow in Figure 5B) overlaps with that of the V0 and V1 pre-MNs (Lai et al., 2016; Griener et al., 2015) as well as with that of the contralateral dI2 neurons (Figures 1 and 5B). To assess the potential spinal targets of dI2 neurons, we inspected the degree of overlap between dI2 synapses and dI2 somata (Figure 5C), ipsilateral pre-MNs (Figure 5D), and contralateral pre-MNs (Figure 5E). The alignment revealed an overlap of dI2 synapses with ipsilateral/contralateral pre-MNs and dI2 neurons (Figure 5C–E), supporting their potential connectivity. Labeling of dI2 synapses coupled with labeling of the above neuronal population showed dI2 synaptic boutons on pre-MNs and dI2 neurons at the lumbar level (Figure 5C–E, Figure 5—figure supplement 1A–C).
The pattern of dI2 collaterals along the entire rostrocaudal axis (Figure 3A–C, Figure 5A) suggests that dI2 neurons innervate contralateral pre-MNs and dI2 neurons at multiple levels. To test this hypothesis, labeling of lumbar dI2 neurons was coupled with labeling of brachial pre-MNs and dI2 somata by injecting PRV into the wing musculature or electroporating a reporter into brachial dI2 neurons, respectively (Figure 5F, Figure 5—figure supplement 1D). dI2 synapses overlapped with the putative targets, and synaptic boutons originating from lumbar-level dI2 neurons were apparent on dI2 neurons and on the contralateral and ipsilateral pre-MNs of the wings (Figure 5F, Figure 5—figure supplement 1D and E).
Neuronal and synaptic labeling experiments showed that lumbar dI2 neurons innervate the cerebellum, lumbar and brachial pre-MNs, and contralateral dI2 neurons. Hence, dI2 neurons may relay peripheral and intraspinal information to the cerebellum and to the contralateral lumbar and brachial motor control centers.
Silencing of dI2 neurons impairs the stability of bipedal stepping
The synaptic input to dI2 neurons and their putative targets implicates them as relaying information about motor activity to the contralateral spinal cord and the cerebellum. Thus, we hypothesized that manipulation of their neuronal activity may affect the dynamic profile of stepping. To study the physiological role of dI2 neurons, we silenced their activity by expressing the tetanus toxin (TeTX) light chain gene, which blocks synaptic transmission (Yamamoto et al., 2003), in the bilateral lumbar dI2 neurons. EGFP was cotargeted in a 2/1 TeTX/EGFP ratio for post hoc analysis of the efficacy of electroporation (Supplementary file 1). Chicks expressing EGFP in dI2 neurons and chicks that did not undergo electroporation were used as controls. To maximize the number of targeted dI2 neurons, we combined genetic targeting with the Foxd3 enhancer and spatial placement of the electrodes at the dorsal lumbar spinal cord (Figure 1—figure supplement 1A). Embryos were electroporated at HH18. Upon hatching, chicks were trained for targeted overground locomotion. The gait parameters of four controls and 5 TeTX-treated chicks were measured while chicks were walking toward their imprinting trainer along a horizontal track (6–20 walking sessions, 5–8 strides each, per chick).
To test whether silencing of dI2 neurons impairs posthatching development and muscle strength, the chicks were weighed at P8, and their foot grip strength was evaluated on the same day. All chicks were of comparable weight (average – 144.7 ± 12.1 g; Supplementary file 1). As a functional measure of foot grip, we tested the ability of the chicks to maintain balance on a tilted mesh surface. TeTX-manipulated chicks and control chicks maintained balance on the tilted surface up to 63–70°, with no apparent statistically significant differences (Supplementary file 1). Thus, manipulation of dI2 neuronal activity did not impair the development or balance or muscle strength.
Analysis of overground locomotion in the control and TeTX-treated chicks revealed no significant differences in swing velocity or striding pattern. A 180° out-of-phase pattern was found during stepping in all the manipulated and control chicks (Figure 6—figure supplement 1A, Table 1). However, substantial differences were scored in stability parameters: TeTX chicks exhibited whole-body collapses during stepping (Figures 6B and C and 7A), a wide-base gait (Table 2), and variable limb movements (Figure 6A, D and E; Figure 7B and C; Figure 7—figure supplement 1; Table 3).
-
Figure 6—source data 1
- https://cdn.elifesciences.org/articles/62001/elife-62001-fig6-data1-v2.xlsx
-
Figure 6—source data 2
- https://cdn.elifesciences.org/articles/62001/elife-62001-fig6-data2-v2.xlsx
-
Figure 6—source data 3
- https://cdn.elifesciences.org/articles/62001/elife-62001-fig6-data3-v2.docx
Whole-body collapses
A collapse was scored as a decline in knee height below 85% of the average knee height at the stance phase of the step (arrow in Figure 6C). We measured the number of collapses in 50–190 steps. In control chicks, collapses occurred in 0.53% ± 0.92% of the steps. In TeTX-manipulated chicks, we observed collapses in 19.46% ± 8.3% of the steps, which was significantly different from the rate in controls (Figure 7A). As some collapses were followed by an overextension (‘overshoot’ in leg elevation), also manifested in the profile of the knee height trajectory during the swing phase (Figure 6D), we further studied the relationship between the two phenomena. In general, there was high variability between chicks in this aspect.
Most collapses (64.9% ± 19.7%) were not preceded or followed by overextension. About 22% of the collapses (22.47% ± 21.5% of collapses, e.g., arrowhead in Figure 6C) were followed by overextension, suggesting a postcollapse compensation in the extensor drive. The rest of the collapses (12.63% ± 7.56%) were preceded by overextension.
Wide-base stepping
A wide-base stance is typical of an unbalanced ataxic gait. The stride width was measured between the two feet during the double stance phase of stepping. The mean stride in TeTX-manipulated chicks 1, 2, 4, and 5 was significantly wider than that in the control chicks, while the width in TeTX3 was similar to that in the controls (Table 2).
Variable limb movements
In stable gait, limb trajectories are consistent from stride to stride. Therefore, we compared the trajectories of knee height and angle of the TMP joint during the swing phase of stepping between control and TeTX-manipulated chicks. Plots of the knee height and TMP angle trajectories during the normalized swing in all the analyzed steps of each chick are shown superimposed in Figure 6D and E, respectively. These data demonstrate that the range of changes in TeTX-manipulated chicks was higher than that in control chicks.
Further analyses revealed that, overall, the control group showed lower knee height and TMP angle ranges than the TeTX-treated group, even though there were differences within groups (Figure 7—figure supplement 1). The average knee height range of the combined control chicks (1.981 ± 0.33) was significantly lower than the range of the combined TeTX-treated chicks (3.109 ± 0.74) (Figure 7B). A similar comparison of the combined ranges of angular excursions of the TMP joint during the normalized swing revealed that the average angle of the control group (49.34 ± 16.03) was significantly lower than the average of the TeTX-treated chicks (77 ± 22.15; Figure 7C).
Since the increased range of changes could be due to the effects of the substantial increase in body collapses during stepping (Figure 7A, see also Figure 6), we excluded steps featuring whole-body collapses and reanalyzed the data. The data summarized in Table 3 show that the significant difference between controls and the TeTX-treated chicks in the range of the knee height and the TMP angle excursions was maintained. Thus, the increase in irregularity in the TeTX-treated chicks is not caused exclusively by the body collapses of the TeTX-treated chicks.
Loss of balance can also arise from slipping, which can originate from a shallow landing angle (Clark and Higham, 2011). The landing angle was characterized as the angle between the ground and the imaginary line connecting the knee joint (which is located near the chicken’s center of mass; Daley and Biewener, 2006) to the TMP joint at the end of the swing phase (Figure 7—figure supplement 2A). Thus, we analyzed the landing angles of the manipulated and control chickens. No significant differences in landing angles were detected (Figure 7—figure supplement 2B). Additionally, we compared the landing angle before a collapse to landing angles not preceding a collapse. We found that the angles preceding a collapse were not smaller and even tended to be slightly larger than the landing angles that did not precede a collapse (Figure 7—figure supplement 2C, p=0.02). These results argue against the possibility that the increased frequency of collapses in the manipulated chickens stemmed from slipping and sliding events.
Overall, the kinematic parameters of the dI2-TeTX-treated chicks demonstrate a reduction in stability during locomotion, indicating a possible role of dI2 in the stabilization of bipedal stepping.
Discussion
The VSCT is thought to provide peripheral and intrinsic spinal information to the cerebellum to shape and update the output of spinal networks that execute motor behavior. The lack of genetic access to VSCT neurons hampers efforts to elucidate their role in locomotion. Using a genetic toolbox to dissect the circuitry and manipulate neuronal activity in the chick spinal cord, we studied spinal INs with VSCT characteristics. The main finding in our study is that dI2 neurons in the chick lumbar spinal cord are commissural neurons that innervate pre-MNs at the contralateral lumbar and brachial spinal levels and granule neurons in the ipsilateral cerebellum. Hence, a subpopulation of dI2 neurons form part of the avian VSCT. Targeted silencing of dI2 neurons leads to impaired stepping in P8 hatchlings. We described the spatial distribution of subpopulations of dI2 neurons, deciphered their connectomes, mapped the trajectory of their projections to the cerebellum, and suggested possible mechanisms for the gait perturbation resulting from their genetic silencing, as discussed below.
The connectome of dI2 neurons
Using the intersection between genetic drivers and spatially restricted delivery of reporters to define lumbar and brachial neurons, we identified several targets of dI2 lumbar neurons. Lumbar dI2 neurons innervate contralateral lumbar dI2 neurons as well as commissural and noncommissural lumbar pre-MNs. This connectivity may influence the bilateral spinal output circuitry at the lumbar cord (e.g., Bras et al., 1988; Jankowska and Hammar, 2013). Moreover, the ascending axons of lumbar dI2 neurons give off gray matter collaterals innervating contralateral dI2 neurons and commissural and noncommissural pre-MNs throughout the brachial spinal cord (Figure 7D and E). Therefore, lumbar dI2 neurons may also contribute to the inter-enlargement coupling described between the segments of the spinal cord that move the legs and the wings (e.g., Valenzuela et al., 1990; Ruder et al., 2016 for forelimb to hindlimb coupling connectivity in mice).
We demonstrated that lumbar dI2 cells receive sensory innervation, premotor inhibitory and excitatory innervation, and innervation from contralateral lumbar dI2 cells (Figure 7D and E). Thus, lumbar dI2 neurons can provide the cerebellum and the contralateral pre-MNs with proprioceptive information, copies of motor commands delivered from the ipsilateral pre-MNs, and integrated information from contralateral dI2 neurons (Figure 7D and E).
Our wiring-decoding studies are based on the availability of enhancer elements that direct expression in specific spinal INs. The lack of enhancer elements for known pre-MNs, such as V2, precluded their analysis. Future experiments using identified regulatory elements that direct expression in two other pre-MN populations, the dI3 and V0 neurons (Avraham et al., 2010b; Gard et al., 2017), will reveal the extent of premotor information relayed by dI2 neurons.
dI2 subpopulations
Our study reveals two anatomically distinct subpopulations of dI2 neurons: precerebellar projection dI2 cells, which also innervate spinal targets along the entire extent of the spinal cord, and propriospinal dI2 cells, which innervate targets within the lumbar level. The laminar and medial/lateral positions of the two dI2 populations are similar, and we did not find subtype-specific expression of the known dI2 TFs. Contrary to our findings, a recent study (Osseward et al., 2021) has shown that tract and propriospinal INs in the mouse spinal cord differ in the localization of their somata along the mediolateral axis and their transcription of TFs. Tract neurons reside in the lateral spinal cord and express group-N TFs, while propriospinal neurons settle at the medial spinal cord and express group-Z TFs. Several reasons may explain the discrepancy between our results and those of Osseward et al. The analysis in mice (Osseward et al., 2021) was applied to numerous cardinal populations of neurons but excluded dI2 neurons. Hence, dI2 may represent an exception. In addition, our study revealed that precerebellar dI2 neurons share two wiring patterns: tract and propriospinal neurons. Applying single-cell RNA sequencing to dI2 neurons will reveal whether short- and long-range targeting in dI2 neurons are characterized by distinct transcriptomes or by a shared N + Z transcriptome.
Physiological role of dI2 neurons
The unclear genetic origin of physiologically equivalent lumbar VSCT neurons has prevented a better understanding of their role in hindlimb locomotion. Our wiring and neuronal-silencing studies implicated dI2 as a significant contributor to the regularity and stability of locomotion in P8 hatchlings. The kinematic analysis of TeTX-treated hatchlings revealed imbalanced locomotion with occasional collapses, increased stride variability, a wide-base gait, and variable limb movements during stepping.
The mechanisms accounting for impaired stepping following dI2 neuron silencing are still unknown. One of the possible mechanisms is that silencing dI2 neurons perturbs the delivery of peripheral and intrinsic feedback to the cerebellum, leading to unreliable updating of the motor output produced by the locomotor networks, thereby impairing bipedal stepping. Another possible mechanism is based on the similarity of the gait instabilities of TeTX-treated hatchlings to ataxic motor disorders. Mammalian VSCT neurons receive descending input from reticulospinal, rubrospinal, and vestibulospinal pathways (Bras et al., 1988; Jankowska and Hammar, 2013). Neurons from the lateral vestibular nucleus have been reported to innervate extensor motoneurons at the lumbar level, as well as INs residing at medial lamina VII (Murray et al., 2018), the location where dI2 neurons were found to reside in our study. Thus, the vestibulospinal tract may convey input directly to the ipsilateral motor neurons and indirectly to contralateral motor neurons through dI2 neurons that innervate contralateral pre-MNs.
The local spinal connections between dI2 neurons and the contralateral pre-MNs and contralateral dI2 neurons may serve an important component of coordinated limb movements. dI2 synapses were found on both ipsilaterally and contralaterally projecting pre-MNs, both within their segmental level and at the brachial level, which regulates the movement of the wings. Thus, dI2 neurons may affect the motor output of the contralateral and ipsilateral sides of the cord by contacting commissural pre-MNs. Specific targeting of dI2 subpopulations – the precerebellar versus the propriospinal dI2 cells – is necessary to determine the relative contribution of dI2 subpopulations to the impaired stepping phenotype. However, there is no available genetic technique for differentially targeting the two subpopulations. In addition, the fact that the precerebellar dI2 neurons also innervate the lumbar spinal cord precludes the use of retrogradely target-derived neuronal activity modifiers.
In summary, our mapping studies of dI2 neurons and their connectomes, followed by characterization of the effects of their silencing on bipedal stepping, offer new insights into the function of dI2 neurons in vertebrates. We suggest that lumbar dI2 neurons not only relay sensory and intrinsic spinal network information to the cerebellum but also act as active mediators of motor functions at the lumbar segments and at the wing-controlling brachial segments of the spinal cord. Further circuit-deciphering studies of the constituents of subpopulations of dI2 cells, their targets, and their descending inputs are required to extend our understanding of the function of dI2 subpopulations in motor control.
Materials and methods
Animals
Fertilized white leghorn chicken eggs (Gil-Guy Farm, Israel) were incubated under standard conditions at 38°C. All experiments involving animals followed the designated policies of the Experiments in Animals Ethics Committee and were performed with its approval.
3D reconstruction and density plot analysis
Request a detailed protocolThe programs for both 3D reconstruction and the density plot analysis were written in MATLAB. The density plots were generated based on cross-sectional images converted to a standard form. The background was subtracted, and the cells were filtered automatically based on their soma size or using a manual approach. Subsequently, two-dimensional kernel density estimation was obtained using the MATLAB function ‘kde2d.’ Finally, unless indicated otherwise, a contour plot was drawn for density values between 20% and 80% of the estimated density range, with six contour lines.
In ovo electroporation
Request a detailed protocolA DNA solution of 5 mg/mL was injected into the lumen of the neural tube at HH stage 17–18 (E2.75–E3). Electroporation was performed using 3 × 50 ms pulses at 25–30 V, applied across the embryo using a 0.5 mm tungsten wire and a BTX electroporator (ECM 830). Following electroporation, 150–300 μL of antibiotic solution containing 100 unit/mL penicillin in Hanks’ balanced salt solution (Biological Industry, Beit-Haemek) was added on top of the embryos. Embryos were incubated for 3–19 days prior to further treatment or analysis.
Immunohistochemistry and in situ hybridization
Request a detailed protocolEmbryos were fixed overnight at 4°C in 4% paraformaldehyde/0.1 M phosphate buffer, washed twice with phosphate buffered saline (PBS), incubated in 30% sucrose/PBS for 24 hr, and embedded in optimal cutting temperature (OCT) compound (Scigen, Grandad, USA). Sections with a thickness of 20 μm were cut on a cryostat. These sections were collected on Superfrost Plus slides and kept at −20°C. For 100 μm sections, spinal cords were isolated from the fixed embryos and subsequently embedded in warm 5% agar (in PBS), and 100 μm sections (E12–E17) were cut with a vibratome. Sections were collected in wells (free-floating technique) and processed for immunolabeling.
The following primary antibodies were used: rabbit polyclonal GFP antibody 1:1000 (Molecular Probes, Eugene, OR, USA), mouse anti-GFP 1:100, goat anti-GFP 1:300 (Abcam), rabbit anti-RFP 1:1000 (Acris), goat anti-ChAT antibody 1:300 (Cemicon, Temecula, CA, USA), mouse anti-synaptotagmin antibody 1:100 (ASV30), mouse anti-Lhx1/5 1:100 (4F2), mouse anti-FoxP4 1:50 (hybridoma bank, University of Iowa, Iowa City, USA), mouse anti-Brn3a 1:50 (Mercury), rabbit anti-Pax2 antibody 1:50 (Abcam), chicken anti-lacZ antibody 1:300 (Abcam), rabbit anti-Calbindin 1:200 (Swant), rabbit anti-VGLUT2 antibody (Synaptic Systems, Göttingen, Germany), goat anti-FoxP2:1000 (Abcam), anti-FoxP1:100 (ABR Synaptic), and rabbit anti-5-HT (Abcam). The following secondary antibodies were used: Alexa Fluor 488/647-conjugated AffiniPure donkey anti-mouse, anti-rabbit, and anti-goat (Jackson) and Rhodamine Red-X-conjugated donkey anti-mouse and anti-rabbit (Jackson). Images were taken under a microscope (Eclipse Ni; Nikon) with a digital camera (Zyla sCMOS; Andor) or captured using the integrated camera of a confocal microscope (FV1000; Olympus).
In situ hybridization was performed as previously described (Avraham et al., 2010a). The following probes were employed: Foxd3, vGlut2, and VIAAT probes were amplified from the cDNA of E6 chick embryos using the following primers. Foxd3: forward TCATCACCATGGCCATCCTG, reverse GCTGGGCTCGGATTTCACGAT. vGlut2: forward GGAAGATGGGAAGCCCATGG, and reverse GAAGTCGGCAATTTGTCCCC. VIAAT: forward CTGAACGTCACCAACGCCATCC, reverse GGGTAGGAGAGCAAGGCTTTG. The T7 RNA polymerase cis-binding sequence was added to the reverse primers.
Laminar division
Request a detailed protocolThe standard forms of the spinal cord (for the crural, sciatic, and brachial plexus levels) were computationally divided into polygons for the different laminae (Martin, 1979). The number of neurons or synapses inside each lamina border was scored using their coordinates.
Light-sheet microscopy dI2::mCherry was electroporated into the embryos at HH stage 17–18. Embryos were removed at E13, and the spinal cord and cerebellum were isolated. The tissue was cleared using the iDISCO technique as described (Renier et al., 2014). The mCherry-expressing neurons were stained by application of an anti-RFP antibody followed by Rhodamine Red-X-conjugated donkey secondary antibody. Each staining step included 3 days of incubation with the antibody and subsequent washing for 2 days. Then, the cleared tissue was divided into three segments: a lumbar spinal segment, a brachial spinal segment, and a segment including the brainstem and cerebellum. Each sample was placed in a quartz imaging chamber (LaVision BioTec) and scanned by a LaVision Ultramicroscope II light-sheet microscope operated by ImspectorPro software (LaVision BioTec). An Andor Neo sCMOS camera was used for 16-bit image acquisition. The imaging was performed at 2× magnification with a 0.5–1 µm step size and a green excitation filter (peak – 525 nm/width − 50 nm). Then, for 3D reconstruction and analysis of the samples, the resulting image z-stacks were converted to IMS format using Imaris File Converter (version 9.5). Because of the sample size, several z-stacks were required for full acquisition of each sample; they were stitched together into one z-stack by Imaris Stitcher. Then, the files were uploaded to Imaris (9.6 version) for advanced visualization and analysis. dI2 axons were tracked using the filament tracer feature in semiautomatic mode. After tracking, Imaris was used to generate videos and snapshots describing different features of the analyzed samples. Finally, text, arrows, and other symbols were added using Adobe AfterEffects software.
Synaptic marker validation
Request a detailed protocolValidation of SV2-GFP reporter specificity was performed by using Imaris software. High-resolution confocal images of spinal cord sections, with clear SV2-GFP reporter expression and synaptotagmin (syn) immunolabeling, were used to quantify the degree of overlap of GFP+ terminals and syn+-labeled boutons. Both signals were three-dimensionally reconstructed, and we used the automatic quantification abilities of Imaris, further validated by additional manual counting, to quantify the number of GFP+ presynaptic terminals containing at least one syn+ bouton. In addition, the volume of GFP+ presynaptic terminals was documented to explore a possible dependence between terminal volume and syn+ bouton containment.
AP staining
Request a detailed protocolThe treated embryos were fixed with 4% paraformaldehyde–PBS for 24 hr at 4°C and washed twice with PBS for 30 min at 4°C. The fixed embryos were incubated at 65°C in PBS for 8–16 hr to inactivate endogenous AP activity. The treated embryos were washed with 100 mM Tris–Cl (pH 9.5) containing 100 mM NaCl and 50 mM MgCl2, and the residual placental alkaline phosphatase activity was visualized by incubating the embryos with NBT/BCIP (Roche) in the same buffer at 4°C for 24 hr. After extensively washing the embryos with PBS–5 mM EDTA, the spinal cord was imaged.
PRV infection and CTB retrograde labeling
Request a detailed protocolFrom the attenuated PRV Bartha strain, we used two isogenic recombinants that express enhanced GFP (PRV152) and monomeric red fluorescent protein (PRV614). The viruses were harvested from Vero cell cultures at titers of 4 × 108, 7 × 108, and 1 × 109 plaque-forming units (PFU/mL). Viral stocks were stored at −80°C. Injections of 3 μL of PRV152 or PRV614 were made into the thigh, pectoralis, or distal wing musculature of E13 or E14 chick embryos using a Hamilton syringe (Hamilton; Reno, NV, USA) equipped with a 33-gauge needle. The embryos were incubated for 36–40 hr and sacrificed for analysis. For spinocerebellar projecting neuron labeling, we used a replication-defective HSV (TK-) that contains a lacZ reporter. The virus was injected into the cerebellum of E12–15 embryos in ovo, and the embryos were incubated for another 40–48 hr. Alternatively, CTB conjugated to Alexa Fluor 647 (Thermo Fisher) was injected into the cerebellum of E12–15 embryos together with the virus for visualization of both precerebellular neurons and the upstream neurons.
Force test
Request a detailed protocolMuscle strength was evaluated using the measurement of the slope at which the chicks fell from a mesh surface as it was gradually tilted up from the horizontal. This test was repeated for each chick at least three times, and the average falling angle was calculated.
Analysis of left-right phase
Request a detailed protocolStride duration was measured as the time from right toe-off/foot-off to the next right toe-off (as a complete stride cycle for the right leg), and the ‘half-cycle’ duration was measured as the time from right toe-off to the time of left toe-off. The following formula was used to calculate the phase: ((LeftToeOff_1 - RightToeOff_1)/(RightToeOff_2 – RightToeOff_1))*360.
Behavioral tests and analysis
Request a detailed protocolThe embryos were bilaterally electroporated and then allowed to develop and hatch in a properly humidified and heated incubator. Afterwards, within 32 hr after hatching, the hatchling chicks were imprinted on the trainer. The P8 chicks were filmed in slow motion (240 fps) while freely walking (side and top views). The following parameters were scored: (1) weight; (2) foot grip strength; and (3) kinematic parameters during overground locomotion: (a) swing velocity, (b) swing and stance duration, (c) phase of footfalls, (d) heights of the knee and TMP joints, (e) angles of the TMP and ankle joints, (f) stride width (distance between feet during the double stance phase), and (g) landing angle.
Using semiautomated MATLAB-based tracking software (Hedrick, 2008), several points of interest were encoded. The leg joints as well as the eye and the tail were tracked. The position of these reference points was used for computational analysis using in-house MATLAB code for calculating different basic locomotion parameters (e.g., stick diagrams, velocity, joint trajectory, angles, range, and elevation), step patterns, and degrees of similarity. The landing angle was calculated as the angle between the imaginary line connecting the knee and the TMP joints and the ground, at the end of the swing phase. Dunnett’s test (Dunnett, 1955) was used to perform multiple comparisons of group means following one-way ANOVA. Circular statistics were used for analyses of angular data utilizing Oriana (KCS, version 4).
Data availability
All data generated or analysed during this study are included in the manuscript and the supporting files.
References
-
Motor and dorsal root ganglion axons serve as choice points for the ipsilateral turning of di3 axonsThe Journal of Neuroscience 30:15546–15557.https://doi.org/10.1523/JNEUROSCI.2380-10.2010
-
Deciphering axonal pathways of genetically defined groups of neurons in the chick neural tube utilizing in ovo electroporationJournal of Visualized Experiments 1:1792.https://doi.org/10.3791/1792
-
Demonstration of initial axon collaterals of cells of origin of the ventral spinocerebellar tract in the catThe Journal of Comparative Neurology 273:584–592.https://doi.org/10.1002/cne.902730412
-
Slipping, sliding and stability: Locomotor strategies for overcoming low-friction surfacesThe Journal of Experimental Biology 214:1369–1378.https://doi.org/10.1242/jeb.051136
-
A multiple comparison procedure for comparing several treatments with a controlJournal of the American Statistical Association 50:1096–1121.https://doi.org/10.1080/01621459.1955.10501294
-
Recent advances in the use of neurotropic viruses for circuit analysisCurrent Opinion in Neurobiology 13:603–606.https://doi.org/10.1016/j.conb.2003.08.001
-
Spinocerebellar projections from the cervical and lumbosacral enlargements in the chicken spinal cordBrain, Behavior and Evolution 76:271–278.https://doi.org/10.1159/000321910
-
Probing diversity within subpopulations of locomotor-related v0 interneuronsDevelopmental Neurobiology 75:1189–1203.https://doi.org/10.1002/dneu.22277
-
A “tool box” for deciphering neuronal circuits in the developing chick spinal cordNucleic Acids Research 42:e148.https://doi.org/10.1093/nar/gku750
-
Serotoninergic and noradrenergic axons make contacts with neurons of the ventral spinocerebellar tract in the catThe Journal of Comparative Neurology 443:310–319.https://doi.org/10.1002/cne.10134
-
Clarke’s column neurons as the focus of a corticospinal corollary circuitNature Neuroscience 13:1233–1239.https://doi.org/10.1038/nn.2637
-
Interactions between spinal interneurons and ventral spinocerebellar tract neuronsThe Journal of Physiology 591:5445–5451.https://doi.org/10.1113/jphysiol.2012.248740
-
Neuronal specification in the spinal cord: Inductive signals and transcriptional codesNature Reviews. Genetics 1:20–29.https://doi.org/10.1038/35049541
-
Direct and indirect spino-cerebellar pathways: Shared ideas but different functions in motor controlFrontiers in Computational Neuroscience 9:75.https://doi.org/10.3389/fncom.2015.00075
-
Making sense out of spinal cord somatosensory developmentDevelopment 143:3434–3448.https://doi.org/10.1242/dev.139592
-
A cytoarchitectonic scheme for the spinal cord of the domestic fowl, Gallus gallus domesticus: Lumbar regionActa Morphologica Neerlando-Scandinavica 17:105–117.
-
Cell type and circuit modules in the spinal cordCurrent Opinion in Neurobiology 56:175–184.https://doi.org/10.1016/j.conb.2019.03.003
-
Information to cerebellum on spinal motor networks mediated by the dorsal spinocerebellar tractThe Journal of Physiology 591:5433–5443.https://doi.org/10.1113/jphysiol.2012.249110
-
Laterality of spinocerebellar neurons in the chicken spinal cordThe Journal of Veterinary Medical Science 74:495–498.https://doi.org/10.1292/jvms.11-0398
-
Stimulation of the brainstem reticular formation evokes locomotor activity in embryonic chicken (in ovoBrain Research. Developmental Brain Research 56:13–18.https://doi.org/10.1016/0165-3806(90)90158-u
-
Reversible suppression of glutamatergic neurotransmission of cerebellar granule cells in vivo by genetically manipulated expression of tetanus neurotoxin light chainThe Journal of Neuroscience 23:6759–6767.
Article and author information
Author details
Funding
Israel Science Foundation (1400/16)
- Avihu Klar
United States - Israel Binational Science Foundation (2017/172)
- Avihu Klar
The Avraham and Ida Baruch Endowment Fund
- Avihu Klar
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank Haya Falk for PRV purification; Alona Katzir, Cole Bendor, Mevaseret Avital, Sapir Shevah, Eitan Yisraeli, Ruth Segal, Fedaa Bazan, and Eden Kimchi for technical assistance; Nadav Yayon for assistance with the light sheet microscopy; and Michael O’Donovan for comments on the manuscript. This work was supported by grants to AK from the Israel Science Foundation (grant no. 1400/16), the US–Israel Binational Science Foundation (grant no. 2017/172), and the Avraham and Ida Baruch endowment fund.
Ethics
All experiments involved with animals were conducted in accordance with the designated Experiments in Animals Ethic Committee policies and under its approval.
Copyright
© 2021, Haimson et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,048
- views
-
- 199
- downloads
-
- 3
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Premature infants with bronchopulmonary dysplasia (BPD) have impaired alveolar gas exchange due to alveolar simplification and dysmorphic pulmonary vasculature. Advances in clinical care have improved survival for infants with BPD, but the overall incidence of BPD remains unchanged because we lack specific therapies to prevent this disease. Recent work has suggested a role for increased transforming growth factor-beta (TGFβ) signaling and myofibroblast populations in BPD pathogenesis, but the functional significance of each remains unclear. Here, we utilize multiple murine models of alveolar simplification and comparative single-cell RNA sequencing to identify shared mechanisms that could contribute to BPD pathogenesis. Single-cell RNA sequencing reveals a profound loss of myofibroblasts in two models of BPD and identifies gene expression signatures of increased TGFβ signaling, cell cycle arrest, and impaired proliferation in myofibroblasts. Using pharmacologic and genetic approaches, we find no evidence that increased TGFβ signaling in the lung mesenchyme contributes to alveolar simplification. In contrast, this is likely a failed compensatory response, since none of our approaches to inhibit TGFβ signaling protect mice from alveolar simplification due to hyperoxia while several make simplification worse. In contrast, we find that impaired myofibroblast proliferation is a central feature in several murine models of BPD, and we show that inhibiting myofibroblast proliferation is sufficient to cause pathologic alveolar simplification. Our results underscore the importance of impaired myofibroblast proliferation as a central feature of alveolar simplification and suggest that efforts to reverse this process could have therapeutic value in BPD.
-
- Developmental Biology
Cells called alveolar myofibroblasts, which have a central role in the development of the lung after birth, receive an orchestrated input from a range of different signaling pathways.