当前位置: 首页 > 医学版 > 期刊论文 > 基础医学 > 细胞生物学杂志 > 2005年 > 第4期 > 正文
编号:11167813
Transcriptional regulation of myotube fate specification and intrafusal muscle fiber morphogenesis
     1 Department of Pathology, Northwestern University, Chicago, IL 60611

    2 Neurology, Northwestern University, Chicago, IL 60611

    3 The Institute for Neuroscience, Northwestern University, Chicago, IL 60611

    Vertebrate muscle spindle stretch receptors are important for limb position sensation (proprioception) and stretch reflexes. The structurally complex stretch receptor arises from a single myotube, which is transformed into multiple intrafusal muscle fibers by sensory axon–dependent signal transduction that alters gene expression in the contacted myotubes. The sensory-derived signal transduction pathways that specify the fate of myotubes are very poorly understood. The zinc finger transcription factor, early growth response gene 3 (Egr3), is selectively expressed in sensory axon–contacted myotubes, and it is required for normal intrafusal muscle fiber differentiation and spindle development. Here, we show that overexpression of Egr3 in primary myotubes in vitro leads to the expression of a particular repertoire of genes, some of which we demonstrate are also regulated by Egr3 in developing intrafusal muscle fibers within spindles. Thus, our results identify a network of genes that are regulated by Egr3 and are involved in intrafusal muscle fiber differentiation. Moreover, we show that Egr3 mediates myotube fate specification that is induced by sensory innervation because skeletal myotubes that express Egr3 independent of other sensory axon regulation are transformed into muscle fibers with structural and molecular similarities to intrafusal muscle fibers. Hence, Egr3 is a target gene that is regulated by sensory innervation and that mediates gene expression involved in myotube fate specification and intrafusal muscle fiber morphogenesis.

    Abbreviations used in this paper: Arc, activity-regulated cytoskeletal-associated protein; DRG, dorsal root ganglion; E, embryonic gestational day; Egr3, early growth response gene 3; ErbB2, erythroblastic leukemia viral oncogene homologue 2; ERE, early growth response element; ERM, Ets-related protein; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; GDNF, glial-derived neurotrophic factor; GPR50, g-coupled protein receptor 50; Hey1, harry enhance of split 1; HSA, human skeletal actin; MOI, multiplicity of infection; NGFR (p75), low affinity nerve growth factor receptor; Nrg1, neuregulin-1; NT-3, neurotrophin-3; Prph1, peripherin 1; Pv, parvalbumin; Rhpn2, rhophilin 2; Sd-MyHC, slow developmental myosin heavy chain; SSTr2, somatostatin receptor type 2.

    Introduction

    Muscle spindles are stretch receptors that mediate myotactic stretch reflexes as well as limb and axial body position sensation (proprioception) in vertebrates. They consist of specialized, encapsulated muscle fibers (intrafusal muscle fibers) that are innervated by sensory (groups Ia and II) and motor (fusimotor) axons to form mechanoreceptors that relay skeletal muscle stretch sensation to the central nervous system. Late prenatal or early postnatal sensory denervation, but not motor denervation, of developing skeletal muscles leads to spindle degeneration, which indicates that normal spindle development is dependent on sensory, but not fusimotor, innervation (Zelena, 1957; Zelena and Soukup, 1973; Kucera and Walro, 1992). Spindle morphogenesis requires intact sensory innervation, and this has been convincingly demonstrated in mice that lack either neurotrophin-3 (NT-3) or its cognate, tyrosine kinase receptor/TrkC (Ernfors et al., 1994; Klein et al., 1994; Tessarollo et al., 1994). Proprioceptive sensory neuron survival during development depends on NT-3–tyrosine kinase receptor C signaling, and when it is abrogated in mutant mice Ia-afferents do not form properly to induce spindle morphogenesis in the periphery. Thus, Ia-afferent–derived signaling has been recognized as an essential aspect of muscle spindle induction for several decades, but not until recently have some specific signaling mediators been identified. For example, neuregulin-1 (Nrg1) appears to be an essential Ia-afferent–derived mediator of spindle morphogenesis, and although Ia-afferents establish contact with myotubes, muscle spindle morphogenesis and some intrafusal muscle fiber–specific genes are not appropriately induced in mice when sensory neurons are deficient in Nrg1 (Hippenmeyer et al., 2002). Similarly, muscle spindle morphogenesis is impaired in mice with a muscle-specific deficiency of erythroblastic leukemia viral oncogene homologue 2 (ErbB2), an essential component of the Nrg1 tyrosine kinase receptor signaling complex (Andrechek et al., 2002; Leu et al., 2003). Together, these results indicate that Ia-afferents induce spindle morphogenesis by releasing Nrg1 to engage ErbB2-dependent signal transduction in a subpopulation of myotubes that receive sensory myoneural contacts.

    Comparatively little is known about the gene regulatory networks that are engaged by ErbB2 signaling in Ia-afferent–contacted myotubes in order to induce their transformation to intrafusal muscle fibers. Presumably, de novo gene expression is required after sensory myoneural contact is established to specify the fate of myotubes that will become morphologically and biochemically distinct intrafusal muscle fibers; to mediate terminal Schwann cell differentiation that will generate the fusiform spindle capsule; and to regulate growth factors that are induced to establish and/or maintain specialized, spindle-related sensory and motor innervation. Accordingly, several transcriptional regulators are induced in myotubes shortly after they are contacted by Ia-afferents, including some Ets-related transcription factors, ER81 (Arber et al., 2000), Pea3 (Livet et al., 2002), Ets-related protein (ERM; Hippenmeyer et al., 2002), and the zinc finger transcription factor early growth response gene 3 (Egr3; Tourtellotte and Milbrandt, 1998). They are all regulated by Ia-afferent–derived Nrg1 after myotube contact (Hippenmeyer et al., 2002), indicating that they are candidate transcription factors that are involved in ErbB2-dependent and myotube-intrinsic gene expression related to spindle morphogenesis. However, whether Ets-related transcription factors regulate gene expression that is required for spindle morphogenesis has been difficult to discern. For example, ER81 is expressed in developing spindles, proprioceptive sensory neurons, and motor neurons. In ER81-deficient mice, the spindle number is altered in some muscles because ER81 has a role in the specification of some proprioceptive neurons (Kucera et al., 2002). In contrast, muscle spindles show no obvious defects in Pea3-deficient mice, despite the fact that it is also coordinately expressed by intrafusal muscle fibers, motor neurons, and sensory neurons (Livet et al., 2002). Finally, ERM expression is not restricted to spindles in the developing muscle, and its role in muscle spindle morphogenesis has not been studied because ERM-deficient mice die before muscle spindle induction (Hippenmeyer et al., 2002).

    Egr3 is a particularly interesting candidate effector molecule of Nrg1–ErbB2 signaling in myotubes. Egr3 is induced in Ia-afferent–contacted myotubes at a developmental time point that coincides with Ia-afferent innervation, and it is not expressed by sensory or motor neurons that innervate them. Moreover, Ia-afferent–contacted myotubes in Egr3-deficient mice fail to differentiate into intrafusal muscle fibers that express an intrafusal muscle fiber–specific, slow developmental myosin heavy chain (Sd-MyHC; Tourtellotte and Milbrandt, 1998; Tourtellotte et al., 2001), suggesting that their differentiation may be impaired. As one of the earliest recognized transcriptional regulators that are induced selectively in myotubes by Ia-afferent innervation, Egr3 may mediate some aspects of Nrg1–ErbB2 signaling that is relevant to myotube fate specification and intrafusal fiber differentiation. However, neither the target genes regulated by Egr3 in Ia-afferent–contacted myotubes nor its potential role in myotube fate specification have ever been directly examined. In this study, we identified many target genes that are regulated by Egr3 in primary myotubes in vitro and demonstrated that some of these genes are regulated by Egr3 in developing spindles in vivo. To examine whether Egr3 has any instructive role in specifying the fate of myotubes to become intrafusal muscle fibers, transgenic mice were generated to express Egr3 in all myotubes independent of Ia-afferent–Nrg1 signaling. These mice were not viable because their skeletal muscle fibers were entirely transformed to contain muscle fibers that were structurally and biochemically similar to intrafusal muscle fibers. Considered together, these results demonstrate that Ia-afferent–Nrg1 signaling in myotubes serves, at least in part, to regulate Egr3 as an important transcriptional regulator of myotube fate specification and intrafusal muscle fiber morphogenesis.

    Results

    Egr3-mediated gene regulation in primary murine myotubes

    Egr3 expression is regulated in myotubes by sensory innervation and Nrg1 signaling in vivo (Tourtellotte et al., 2001; Hippenmeyer et al., 2002; Jacobson et al., 2004). Ia-afferents contact myotubes in Egr3-deficient mice; in the absence of Egr3-mediated gene regulation, they do not properly differentiate into spindles, which leaves the mice with profound sensory ataxia and proprioceptive deficits (Tourtellotte and Milbrandt, 1998). To identify target genes that are potentially regulated by Egr3 in Ia-afferent–contacted myotubes, primary myoblasts were isolated from newborn, wild-type mice and were differentiated into multinucleated myotubes. One population of myotubes was infected with an adenovirus that coexpressed enhanced GFP (EGFP) and full-length, transcriptionally active wild-type Egr3 (Egr3WT), whereas another population of myotubes was infected with an adenovirus that coexpressed EGFP and transcriptionally inactive truncated Egr3 (Egr3Tr; Fig. 1 A). Primary myotubes that were infected with 100% efficiency expressed high levels of transgenic protein and showed no evidence of toxicity (Fig. 1 B). Endogenous Egr3 was not detected in cultured myotubes, similar to myotubes that neither receive Ia-afferent innervation nor become intrafusal muscle fibers in vivo (Fig. 1 B, bottom right; Tourtellotte and Milbrandt, 1998). To identify genes regulated by Egr3 in this primary myotube–cellular context, gene expression in the adenovirus-infected myotubes was compared using Affymetrix microarray analysis. The experiment was repeated twice, and a four-way analysis of the gene expression datasets identified 83 unique genes that were up-regulated and five unique genes that were down-regulated by Egr3WT in primary myotubes. Egr3 appears to be a transcriptional activator in the cellular contexts where it has been studied, which is consistent with the observation that many more genes were up-regulated by Egr3WT (Table I, 83 genes, 2–244-fold) than were down-regulated (Table II, five genes, two- to threefold). Although most of the 83 genes identified have not been evaluated further, up-regulation was confirmed in 15/15 genes examined by using a real-time PCR analysis on cDNA samples that were prepared from independently infected primary myotubes (Fig. 1 D, results from 12 genes shown; Table I).

    Egr3-mediated gene regulation in developing muscle spindles

    To examine whether Egr3 target genes identified in vitro may also be regulated by Egr3 in a similar cellular context in developing muscle spindles in vivo, we used in situ hybridization to characterize gene expression patterns in embryonic wild-type mouse muscles. Egr3 and Pea3 are transcription factors regulated by Ia-afferent–Nrg1 signaling in myotubes, and Ia-afferents contact myotubes in both Egr3-deficient and Pea3-deficient mice (Tourtellotte and Milbrandt, 1998; Livet et al., 2002). Pea3 is not necessary for spindle morphogenesis, but it is a useful marker for localizing Ia-afferent–contacted myotubes in the developing skeletal muscle (Hippenmeyer et al., 2002). We found that Pea3 was not regulated by Egr3 because Pea3 expression was identified in discretely labeled myotubes in Egr3-deficient mice (Hippenmeyer et al., 2002). However, Pea3-expressing myotubes were reduced by 66% in Egr3-deficient muscles at the embryonic gestational day (E) 16.5 (Fig. 2 A) and were reduced by 86% at E18.5 (Fig. 2 B). Because Pea3 is a marker of Ia-afferent–myotube contact, the decreased number of Pea3-expressing myotubes in Egr3-deficient mice is consistent with previous observations that demonstrate initial Ia-afferent contact of myotubes, progressive loss of afferents during prenatal development, and complete loss of afferents in adult Egr3-deficient mice (Tourtellotte and Milbrandt, 1998; Tourtellotte et al., 2001). Because some Ia-afferent–contacted, Pea3-expressing myotubes were present in prenatal Egr3-deficient muscles, it was possible to examine potential Egr3 target gene expression within them. Accordingly, Egr3-dependent gene expression in skeletal muscle should be restricted to Ia-afferent–contacted myotubes where Egr3 expression is restricted, and it should be abrogated in Pea3-expressing (Ia-afferent contacted), Egr3-deficient myotubes. We examined six target genes that were identified by Affymetrix microarray analysis (and were confirmed to be up-regulated in myotubes by real-time PCR analysis) using in situ hybridization on the muscle. Although three genes (rhophilin 2 [Rhpn2]; g-coupled protein receptor 50 [GPR50]; and activity-regulated cytoskeletal-associated protein [Arc]) showed no specific localization to muscle spindles (unpublished data), the three other genes, including two growth factor receptors (low affinity nerve growth factor receptor [NGFR (p75)]; and somatostatin receptor type 2 [SSTr2]) and the type III intermediate filament peripherin 1 (Prph1), appeared to be expressed selectively by spindles. In E18.5 wild-type skeletal muscle, NGFR (p75) (Fig. 3 A), SSTr2 (Fig. 3 B), and Prph1 (Fig. 3 C) expression was restricted to muscle spindles similar to the expression pattern of Egr3. In skeletal muscles from E18.5 Egr3-deficient mice, NGFR (p75), SSTr2, and Prph1 were not expressed in Pea3-expressing myotubes (Fig. 3, A–C, bottom). Hence, NGFR (p75), SSTr2, and Prph1 represent new intrafusal muscle fiber–restricted genes that are regulated by Egr3. These results indicate that the target genes identified in vitro may at least partially represent the Egr3-mediated gene regulatory network involved in spindle morphogenesis.

    A role for Egr3 in myotube fate specification and intrafusal muscle fiber morphogenesis

    During skeletal muscle development, Ia-afferents contact a relatively small number of myotubes to engage the Nrg1–ErbB2 signaling that specifies whether they will become intrafusal muscle fibers. Egr3 appears to have a role in regulating at least some target genes that are involved in spindle morphogenesis, but it is not known what role it may have in myotube fate specification that is mediated by Ia-afferent–Nrg1 signaling. To examine whether Egr3 has an instructive role in myotube fate specification, transgenic mice were generated using the human skeletal actin (HSA) promoter to enforce Egr3 expression in all skeletal myotubes independent of Ia-afferent–Nrg1 regulation (Fig. 4 A). It was not possible to generate transgenic founders that could be used to establish lines of transgenic mice because they were not viable. To facilitate the analysis of the transgenic founders, real-time PCR was used to examine the level of Egr3 expression from forelimb-derived RNA in every embryo that was generated. Animals were analyzed based upon their genotype and level of Egr3 expression by grouping them into those that did not carry the transgene (Tg–; wild type), those that carried the transgene but did not express transgenic Egr3 above endogenous levels (Tg+/), those that expressed low levels of transgenic Egr3 (Tg+/L; 3–10-fold increased expression relative to Tg– mice), and those that expressed high levels of transgenic Egr3 (Tg+/H; >10-fold increased expression relative to Tg– mice; Fig. 4 B). Late gestation embryos that expressed the transgene (Tg+/L and Tg+/H) were consistently smaller than their transgene negative (Tg–) or nontransgene-expressing (Tg+/) litter mates (Fig. 4 C). E18.5 transgene-expressing mice had beating hearts, flexion contractures of their extremities, no spontaneous or evoked motor activity, and were not viable at birth (Fig. 4 D). The level of Egr3 expression in particular founder mice correlated with the relative number of muscle fibers that expressed Egr3. Egr3 expression in E18.5 wild-type skeletal muscle was restricted to developing spindles, which were relatively few in number (Fig. 4 E, arrows), whereas Egr3 was expressed by most of the muscle fibers in transgene-expressing mice (Fig. 4 F).

    Transgene-expressing mice had poorly formed skeletal muscles that were generally more poorly formed in Tg+/H muscles than in Tg+/L muscles. As expected, muscles from mice that did not express the transgene (Tg+/) appeared normal. In E18.5 wild-type hindlimbs, discrete muscles were readily identified (Fig. 5 A), whereas distinct muscles were largely absent in transgene-expressing hindlimbs (Fig. 5 B). Wild-type muscle fibers contained skeletal myofilaments and subsarcolemmal nuclei that are typical of immature skeletal myofibers (Fig. 5, C and E). Occasional spindles were identified that consisted of intrafusal muscle fibers with scant myofibrillary architecture and internal nuclei. Sections through the equatorial plane of the spindles showed prominent annulospiral axons (Ia-afferents) circumscribing the intrafusal muscle fibers within rudimentary spindle capsules (Fig. 5 E, arrowhead). In contrast, all muscle fibers from transgene-expressing mice contained scant myofilaments with centrally located nuclei. Although most were single multinucleated fibers that had myofibrillary structure and nuclear features indistinguishable from intrafusal muscle fibers, many were aggregated into clusters that were similar to intrafusal muscle fibers within spindles (Fig. 5, D and F [arrow]). Most of the intrafusal-like muscle fibers were not encapsulated or innervated as determined from an analysis of serially sectioned, resin-embedded hindlimbs. Occasional normal-appearing spindles were identified in muscles from transgene-expressing mice that had the typical configuration of three to four intrafusal muscle fibers, which were circumscribed by annulospiral sensory axons and rudimentary capsules (Fig. 5 F, arrowhead). Thus, Egr3 expression in skeletal myotubes independent of other Ia-afferent signaling has the remarkable capacity to transform them into muscle fibers that are structurally similar to intrafusal muscle fibers, presumably orchestrated by the network of genes it regulates.

    Intrafusal-like muscle fibers in HSA–Egr3 transgenic muscle are not innervated

    Intrafusal-like muscle fibers in Egr3 transgene-expressing muscles do not have fusiform capsules, which may derive from terminal Schwann cells of innervating axons. To examine whether the transformed muscle fibers were innervated by either fusimotor or proprioceptive axons, parvalbumin (Pv) immunohistochemistry was used as a selective marker for proprioceptive neurons and Ia-afferents (Tourtellotte and Milbrandt, 1998), and ATP13 immunohistochemistry was used as a marker for both Ia-afferent and fusimotor axons in skeletal muscles (Dobretsov et al., 2003). There was no quantitative difference in the number of proprioceptive neurons in the Tg– or Tg+/H fifth lumbar dorsal root ganglion (DRG; Fig. 6 A), and there was no qualitative difference in the number of Pv+ neurons in any first through fourth cervical or first through fifth lumbar DRG that were examined, suggesting that the number of innervated and structurally normal spindles (Fig. 5 F, arrowhead) was normal in Tg+/H mice. However, whereas robust Pv+ Ia-afferent axons terminated on spindles in Tg– mice (and occasional normal-appearing spindles in Tg+/H mice), the intrafusal-like muscle fibers from Tg+/H mice showed no evidence of Ia-afferent innervation (Fig. 6 B). Similarly, the transformed muscle fibers were not able to sustain motor innervation as indicated by the near complete loss of motor neurons in the spinal cord and the absent ventral roots in Tg+/H mice (Fig. 6 C). In the periphery, ATP13 identified fusimotor and Ia-afferent axons that terminated selectively on normal spindles, whereas there was no evidence of such innervation in the transformed muscle fibers in Tg+/H muscles (Fig. 6 D).

    Pea3 is expressed in myotubes that have been contacted by Ia-afferents, but in Tg+/H muscle there was no change in the number of Pea3-expressing myotubes or the expression level of Pea3 in transgene-expressing muscle (Fig. 6 E). These results are consistent with our findings that Pea3 is not regulated by Egr3 (Fig. 2), that some morphologically intact spindles are innervated in transgene-expressing muscles (Fig. 5 F, arrowhead), and that the transformed muscle fibers are not innervated by Ia-afferents that would otherwise be expected to induce high levels of Pea3 expression. Thus, myotubes that express Egr3 in the absence of other Ia-afferent–derived signaling are morphologically similar to intrafusal muscle fibers, but they are devoid of intact fusiform capsules and are incapable of attracting and/or sustaining their innervation.

    Muscle fibers from HSA–Egr3 transgenic mice express genes similar to those of intrafusal muscle fibers

    Enforced expression of Egr3 in skeletal myotubes is sufficient to transform them into muscle fibers that structurally resemble intrafusal muscle fibers. To examine whether they also express genes that are characteristic of intrafusal muscle fibers, several well-established intrafusal muscle fiber markers and some newly identified markers were examined in HSA–Egr3 transgenic muscles. In wild-type muscle, Sd-MyHC is expressed by the Ia-afferent contact of myotubes as one of the earliest markers of intrafusal muscle fiber formation. Sd-MyHC is not expressed in Ia-afferent–contacted, Egr3-deficient myotubes, suggesting that they are impaired in their capacity to differentiate into intrafusal muscle fibers (Tourtellotte and Milbrandt, 1998; Tourtellotte et al., 2001). Sd-MyHC was localized to intrafusal muscle fibers in developing spindles in wild-type (Tg–) muscle as previously reported, whereas it was expressed by 15-fold more fibers in Tg+/H muscles, indicating that some myotubes had acquired at least one phenotype of intrafusal muscle fibers (Fig. 7 A). The neurotrophins NT-3 and glial-derived neurotrophic factor (GDNF) are preferentially expressed by intrafusal muscle fibers in E18.5 muscle and are required for their sensory and fusimotor innervation, respectively (Ernfors et al., 1994; Whitehead et al., 2005). Relative to Tg– muscles, neither GDNF nor NT-3 was up-regulated in Tg+/ muscles, but both were up-regulated 4- and 5.2-fold, respectively, in Tg+/H muscles (Fig. 7 B). The newly identified Egr3 target genes NGFR (p75), SSTr2, and Prph1 were not up-regulated in Tg+/ muscles relative to Tg– muscles, whereas they were all up-regulated in Tg+/H muscle 6.4-, 6.7-, and 4-fold, respectively (Fig. 7 C). Finally, several other genes identified as potential Egr3 target genes in primary myotubes (Fig. 1 D) were also up-regulated in Tg+/H muscles. ATP13, harry enhance of split 1 (Hey1), and 2810417MoRik were not up-regulated in Tg+/ muscles relative to Tg– muscles, whereas they were up-regulated 11.2-, 10.3-, and 7.3-fold, respectively. Together, these data demonstrate that enforced Egr3 expression in skeletal myotubes leads to their transformation into muscle fibers that are structurally similar to intrafusal muscle fibers and that express many genes characteristic of them.

    Discussion

    During skeletal muscle development, Ia-afferents specify the fate of myotubes by engaging novel gene transcription that transforms them into biochemically and functionally distinct intrafusal muscle fibers. Recent studies have refined our understanding of this classic nerve–muscle induction process by demonstrating that Nrg1 produced by Ia-afferents and signaling through the Nrg1 tyrosine kinase receptor ErbB2 in contacted myotubes is essential to engage the spindle morphogenesis program in vivo (Andrechek et al., 2002; Hippenmeyer et al., 2002; Leu et al., 2003) and in vitro (Jacobson et al., 2004). Moreover, Nrg1 signaling from Ia-afferents is necessary to induce several transcription factors (Egr3, Pea3, ER81, and ERM), which represent at least some proximal mediators of the gene regulatory networks that are engaged by Ia-afferent contact with myotubes. We focused on the role of Egr3 in mediating particular aspects of Ia-afferent (Nrg1–ErbB2) signal transduction because our earlier work indicated that Egr3 is required for normal spindle morphogenesis (Tourtellotte and Milbrandt, 1998). Its expression is restricted to Ia-afferent–contacted myotubes, and Egr3-deficient, Ia-afferent–contacted myotubes do not acquire characteristic structural features or phenotypic markers of intrafusal muscle fibers (Tourtellotte et al., 2001; Chen et al., 2002).

    Egr3 gene regulation in Ia-afferent–contacted myotubes

    Egr3 regulates a relatively large network of genes within myotubes, and although not all of the genes identified were characterized further in this study, some of them were regulated by Egr3 in developing intrafusal muscle fibers. We found that NGFR (p75), SSTr2, and Prph1 are regulated by Egr3 in intrafusal muscle fibers. Of the 83 genes identified by expression analysis in myotubes, we verified that 15/15 genes tested were up-regulated, suggesting that the list of Egr3 target genes has a high degree of specificity. Moreover, the Egr3-mediated gene regulatory network is evidently complex because the target genes have tremendous functional diversity, including processes such as transcriptional regulation, intracellular signal transduction, protein processing, and cytoskeletal organization. Interestingly, however, not all genes that were up-regulated by Egr3 in myotubes in vitro were Egr3 targets in the developing intrafusal muscle fibers in vivo because neither GPR50, Rphn2, nor Arc was expressed by them. Thus, enforced expression of Egr3 in myotubes leads to target gene expression that is not necessarily restricted to the myotube–cellular context. In support of this concept, studies in progress indicate that Arc is directly regulated by Egr3 in neurons, despite the fact that it is not expressed by intrafusal muscle fibers where Egr3 is expressed at high levels (unpublished data). These results indicate that the cellular context appears to be important in providing a permissive environment for expression of a particular repertoire of Egr3 target genes that are relevant to intrafusal muscle fiber morphogenesis. Whether the target genes identified are directly regulated by Egr3 or whether they are expressed by the activation of downstream mediators of Egr3 is currently not known. Moreover, it is not known whether any of these newly identified Egr3 target genes are essential for spindle morphogenesis. An analysis of additional Egr3 target genes by in situ hybridization will be necessary to further define the Nrg1–ErbB2–Egr3 regulatory axis that mediates intrafusal muscle fiber morphogenesis.

    Ia-afferents have a broad regulatory influence over spindle morphogenesis, and Egr3 regulation is apparently only part of the complex morphogenetic signaling pathway. Many of the genes that are up-regulated during the early aspects of muscle spindle morphogenesis, such as Egr3, ERM, ER81, and Pea3, are independently regulated (Hippenmeyer et al., 2002). Thus, multiple parallel regulatory networks likely exist, and they may have distinct roles during spindle morphogenesis. For example, the role of ERM during spindle morphogenesis has not been determined because ERM-deficient mice die before spindle morphogenesis is initiated, but its function in muscles is not likely to be restricted to spindle morphogenesis because it is expressed by intrafusal and nonintrafusal muscle fibers. Pea3, which neither regulates Egr3 nor is regulated by Egr3, is not necessary for spindle morphogenesis despite the fact that in muscles its expression is restricted to intrafusal muscle fibers (Livet et al., 2002). Pea3, which is also expressed in motor neurons and is regulated by GDNF, has a role in patterning motor innervation rather than in spindle morphogenesis itself (Haase et al., 2002; Livet et al., 2002). Similarly, ER81 has a role in sensory patterning, which may explain the partial disruption of spindle morphogenesis in ER81-deficient mice secondary to altered sensory innervation on which spindle morphogenesis depends (Kucera et al., 2002). In contrast, Egr3 appears to have a selective role in mediating some aspects of Ia-afferent signaling that are related to intrafusal muscle fiber differentiation because its restricted expression to spindles is required for their normal morphogenesis. Thus, Egr3-mediated gene regulation appears to be necessary for myotubes to acquire a normal intrafusal muscle fiber phenotype.

    Myotube fate specification mediated by Egr3 target gene regulation

    Although Egr3 expression is necessary for normal intrafusal muscle fiber morphogenesis, it also appears to be sufficient to drive some aspects of intrafusal muscle fiber differentiation from otherwise undifferentiated skeletal myotubes. Skeletal myotubes were transformed into muscle fibers that had structural and molecular similarities to intrafusal muscle fibers when Egr3 expression was enforced independently of Ia-afferent signaling. Moreover, several established intrafusal muscle fiber markers, including Sd-MyHC, NT-3, and GDNF were up-regulated in HSA–Egr3 transgene-expressing muscles. Similarly, three new Egr3 target genes that were selectively expressed by intrafusal muscle fibers (NGFR [p75], SSTr2, and Prph1) and three additional Egr3 target genes (ATP13, Hey1, and 2810417MoRik) were all up-regulated in transgene-expressing muscles. However, the transformed muscle fibers did not acquire the motor or sensory innervation that is characteristic of normal intrafusal muscle fibers, suggesting that Egr3 has a role that is distinct from other Ia-afferent–mediated pathways that may be required to establish and/or maintain these interactions. Moreover, the transformed myotubes did not acquire fusiform capsules that may normally be derived from terminal Schwann cells of the innervating sensory axons. Finally, somatic motor neurons were almost completely absent in HSA–Egr3 transgenic mice, indicating that the transformed muscle fibers were not capable of sustaining skeletomotor innervation that would have been present had they not expressed Egr3 and differentiated into intrafusal muscle fibers. It will be necessary to conditionally express Egr3 in skeletal myotubes in order to study the fate of motor and sensory innervation during development because the transgenic founders died and could be generated only in small numbers.

    Sensory control over the morphogenesis of many mechanoreceptors is recognized as a common theme in their ontogeny. A detailed understanding of the gene regulatory mechanisms that are governed by sensory axon innervation of muscle spindles may be applicable to mechanisms used during the ontogeny of other sensory mechanoreceptors. These studies define a set of Egr3-regulated target genes, some of which appear to be relevant to intrafusal muscle fiber morphogenesis. Egr3 is regulated by Nrg1–ErbB2 signaling, and it appears to mediate gene expression that is necessary and sufficient to specify the fate of Ia-afferent–contacted myotubes to become intrafusal muscle fibers and to mediate their morphogenesis.

    Materials and methods

    Animals and preparation of tissues

    Egr3-deficient mice were mated on a C57BL/6/129SvJ background and were generated and genotyped as previously described (Tourtellotte and Milbrandt, 1998; Whitehead et al., 2005).

    The HSA–Egr3 transgenic mice were generated by using the previously characterized HSA promoter to express Egr3 in skeletal myotubes (Muscat et al., 1992). The transgenic construct was derived from pBSX-HSAvpA (Crawford et al., 2000; provided by J. Chamberlain, University of Washington, Seattle, WA) by cloning a full-length, HA-tagged Egr3 cDNA into the NotI site of pBSX-HSAvpA. The construct was injected into B6SJL embryo pronuclei and transplanted into recipient females using standard procedures (provided by L. Doglio and R. Alvarez, Northwestern University, Chicago, IL). The transgene-expressing founder mice were not postnatally viable, and, therefore, recipient females were killed to analyze founder embryos at various gestational ages. Each embryo was genotyped using primers 5'-AAACTCGCCGAGAAGCTGCCGGTG-3' (sense) and 5'-TGTGCATCCATCGCTAGCTCGCTC-3' (antisense) to amplify a 230-bp fragment of the endogenous Egr3 genomic locus. A third primer, 5'-TCAGGGAGACGTGGAGGCCATGTA-3', was used to amplify a 413-bp fragment of the Egr3 cDNA transgene (nonintron containing) when paired with the sense probe in a multiplex PCR reaction.

    The day of vaginal plug for all of the mice was considered as E0.5. The embryos were isolated by cesarean section, were rapidly decapitated, were fixed at 4°C in 100 mM phosphate-buffered 4% PFA, and were either cryoprotected for 24 h in 30% phosphate-buffered sucrose and frozen or were directly embedded in paraffin. Serial paraffin sections (8 μm) or frozen sections (15 μm) were generated. Some tissues were embedded in resin, in which case the fixative was substituted for 2.4% glutaraldehyde and 1% PFA. The tissues were dehydrated in increasing concentrations of ethanol that was infiltrated with Embed 812 resin (Electron Microscopy Sciences) and heat cured at 60°C overnight. Resin sections were cut at 1 μm and were counterstained with toluidine blue (Sigma-Aldrich).

    Myoblast isolation and myotube culture

    Hindlimb skeletal muscle tissues were isolated from wild-type, postnatal (3–5) mice. Myoblasts were isolated by enzymatically dissociating (dispase, grade II, 2.4 U/ml; Boehringer), and the muscle tissues were triturated according to previously published procedures (Rando and Blau, 1994). Myoblasts were purified by multiple rounds of preplating on collagen-coated tissue culture plates (0.01% type I; Sigma-Aldrich) and expanded in growth medium (Ham's F10; Mediatech), 20% FBS, and 2.5 ng/ml basic FGF (Promega). Purified primary myoblasts were plated onto collagen-coated plates and were differentiated into multinucleated myotubes for 10 d using differentiation medium (Ham's F10 and 5% horse serum).

    Recombinant adenovirus preparation, characterization, and myotube infection

    Recombinant adenoviruses were generated using homologous recombination in Escherichia coli as previously described (He et al., 1998). To generate an adenovirus that expressed transcriptionally active Egr3, the full-length rat Egr3 cDNA (GenBank/EMBL/DDBJ accession no. NP058782, aa 1–387) was cloned into the BglII site of pAdTrackCMV. Similarly, a transcriptionally inactive, COOH-terminal truncation of Egr3 that lacked the entire three zinc finger DNA-binding and COOH-terminal domains (GenBank NP058782, aa 1–245) was cloned into the BglII site of pAdTrackCMV. Each recombinant pAdTrack shuttle vector, including pAdTrackCMV without a cDNA insert, was recombined by homologous recombination into the adenoviral genomic plasmid (pAdEasy) in BJ5183 E. coli to generate a recombinant, replication-deficient adenoviral genome plasmid. The three replication-deficient viruses (EGFP, Egr3WT, and Egr3Tr) were packaged and amplified in transfected HEK-293 cells, purified, concentrated on cesium chloride gradients, and titered using EGFP fluorescence and TCID50.

    NIH-3T3 fibroblasts were infected with the recombinant adenoviruses at a multiplicity of infection (MOI) of 10–20 to obtain 70–80% infection efficiency. Total cell lysates were obtained 24 h after infection, and a Western blot using an NH2-terminal Egr3, rabbit polyclonal antibody (O'Donovan et al., 1998) demonstrated a high expression of full-length and truncated Egr3 protein as expected. Myoblasts were cotransfected with an early growth response element (ERE) containing a firefly luciferase reporter plasmid (Swirnoff and Milbrandt, 1995) and a Renilla luciferase reporter plasmid (pRL-TK; Promega) for transfection efficiency control. 16 h after plasmid cotransfection, the cells were infected with either control EGFP, Egr3WT, or Egr3Tr adenoviruses at an MOI of 10 to obtain 80–90% infection efficiency. 8 h after infection, a dual-reporter luciferase assay was performed to quantify ERE activation by the virally produced proteins according to the manufacturer's specifications (Promega).

    Primary myotubes that differentiated for 10 d in vitro were infected with either Egr3WT or Egr3Tr adenovirus at an MOI of 100 to obtain 100% infection efficiency. Although myotubes were relatively resistant to infection and required a relatively high MOI for optimal infection, no evidence of toxicity was noted 24 h after infection at the time that total RNA was extracted from the infected myotubes.

    Gene expression profiling

    5 μg of total RNA was converted to cDNA using the superscript reverse transcriptase (Invitrogen) and the T7-Oligo (dT) promoter primer kit (Affymetrix, Inc.). The cDNA was purified using the GeneChip sample cleanup module (Affymetrix, Inc.) and was used for the in vitro synthesis of biotin-labeled cRNA using the GeneChip expression 3'-amplification reagents for IVT labeling (Affymetrix, Inc.) at 37°C for 16 h. cRNA was fragmented into 35–200–bp fragments using a magnesium acetate buffer (Affymetrix, Inc.). 10 μg of labeled cRNA were hybridized to GeneChip mouse expression arrays (430A and 430B; Affymetrix, Inc.) for 16 h at 45°C. The GeneChips were washed and stained according to the manufacturer's recommendations (Affymetrix, Inc.) using the GeneChips fluidics station (model 450; Affymetrix, Inc.). This procedure included washing the chips with phycoerythrin-streptavidin, performing signal amplification by a second staining with biotinalyted antistreptavidin, and performing a third staining with phycoerythrin-streptavidin. Each chip was scanned using the GeneChips scanner (model 3000; Affymetrix, Inc.). Signal intensity and detection calls were generated using the GeneChip operating software (Affymetrix, Inc.). The absolute intensity values of each chip were scaled to the same target intensity value of 150 in order to normalize the data for interarray comparisons. Four chip comparisons were generated using the samples from transcriptionally inactive Egr3 (Egr3Tr)-infected myotubes as the baseline. A list of up-regulated and down-regulated genes with a greater than twofold change was generated using an iterative comparison analysis (Chen et al., 2000). This method has been shown by the National Institutes of Health‘s (NIH) Microarray Consortium to increase specificity (>80% true positive rates) in the resultant gene list. The annotated MAGE-ML and image fields are available at http://arrayconsortium.tgen.org.

    Real-time PCR

    Total RNA was isolated from the primary myotubes or mouse muscle tissue using TRIzol (Invitrogen), and 0.5–1.0 μg was reverse transcribed using random octomer priming and powerscript reverse transcriptase according to the manufacturer's specifications (BD Biosciences). Real-time PCR was performed on a sequence detector (model SDS5700; Applied Biosystems) using Syber green fluorescence chemistry (Molecular Probes). Non–intron-spanning primers that amplified the coding sequence were designed for each target gene. For each primer pair, melt curves were generated to identify the optimal temperature to quantify the total fluorescence after each amplification cycle so that only the gene-specific amplicon fluorescence was measured during the PCR reaction. For each primer pair, serial dilutions of mouse genomic DNA were used to generate a standard curve, to characterize their linear dynamic range, and to provide relative DNA concentrations that correlated with gene expression at particular threshold cycle numbers in the linear range of amplification. For each cDNA sample, the relative level of the target gene and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) expression was interpolated from the gene-specific standard curves. The expression of each target gene was normalized to the level of GAPDH expression, and the normalized expression values were again normalized to a reference condition, such as gene expression in nontransgenic muscle or Egr3Tr-infected myotubes, to assess fold change in expression that is relative to that condition.

    Non–intron-spanning primer sequences were used for real-time PCR analysis. They amplified the following: Egr3 (GenBank NM_018781.1, nt 163–563), GAPDH (GenBank NM_001001303, nt 210–635), NGFR/p75 (GenBank BC038365, nt 1123–1369), SSTr2 (GenBank NM_009217, nt 603–859), Prph1 (GenBank NM_013639, nt 109–616), ATP13 (GenBank NM_144921, nt 3074–3293), Hey1 (GenBank NM_010423, nt 669–670), 2810417MoRik (GenBank NM_026516, nt 910–1209), Sbsn (GenBank NM_172205, nt 36–333), Arc (GenBank NM_018790, nt 1175–1483), GPR50 (GenBank AF065145, nt 256–567), Rhpn2 (GenBank NM_011872, nt 1821–2077), KLK7 (GenBank NM_011872, nt 94–396), Krtdap (GenBank AB011028, nt 8–297), Ntf3/NT-3 (GenBank NM_008742, nt 233–633), GDNF (GenBank D88351, nt 2552–2769), SNAP25 (GenBank NM_011428, nt 191–498), Ube3a (GenBank NM_173010, nt 964–1264), Clcn2 (GenBank NM_009900, nt 2710–3009), and Neud4 (GenBank NM_013874, nt 1327–1634).

    In situ hybridization and immunohistochemistry

    In situ hybridization was performed using digoxygenin-labeled riboprobes and previously published standard protocols (Darby, 2000). In situ hybridization probes were generated to label coding regions from the published mRNA sequences. The PCR-generated amplicons were polished (End-it; Epicentre), were subcloned into the EcoRV site of Bluescript (Stratagene), and were sequence verified. Sense and antisense riboprobes were synthesized using in vitro transcription and digoxygenin-labeled dUTP. In all cases, the sense probe was used on parallel tissue sections as a control for nonspecific labeling, and antisense probes were used to examine gene-specific expression.

    The probes used for in situ hybridization spanned the following coding sequences: Egr3 (GenBank NM018781, nt 890–1389), Pea3 (GenBank NM_008815, nt 533–1082), GPR50 (GenBank, AF065145, nt 157–656), Rhpn2 (GenBank NM_027897, nt 1381–1890), Arc (GenBank NM_018790, 577–1082), NGFR (p75) (GenBank BC038365, nt 253–767), SSTr2 (GenBank NM_009217, nt 685–1185), and Prph1 (GenBank NM_013639, nt 109–616). Wild-type muscle spindles were identified by their characteristic morphological appearance (spindle capsules surrounding muscle fibers with internal nuclei). In Egr3-deficient mice, prospective "spindles" were identified by Pea3 expression, which marks Ia-afferent–contacted myotubes and intrafusal muscle fibers.

    Immunoperoxidase histochemistry was performed using specific rabbit polyclonal antibodies that cross-reacted with the COOH-terminal region of Egr3 (sc-190; Santa Cruz Biotechnology, Inc.), the NH2-terminal region of Egr3 (O'Donovan et al., 1998; provided by J. Baraban, John Hopkins University, Baltimore, MD), Pv (R301; provided by K. Baimbridge, University of British Columbia, Canada), ATP13 (Upstate Biotechnology), and Sd-MyHC (S46; provided by D. Fischman, Cornell University, New York, NY). A biotinylated, antirabbit secondary antibody (Jackson ImmunoResearch Laboratories) and avidin–biotin complex histochemistry using diaminobenzidine as a chromagen was performed according to the manufacturer's specifications (Vector Laboratories). In some cases, an antirabbit, Cy3-conjugated secondary antibody was substituted (Jackson ImmunoResearch Laboratories). All photographs were acquired with a digital camera (model RT-Slider; Spot) attached to a microscope (model E600; Nikon), and the images were processed using Adobe Photoshop software.

    Neuron counts

    Stereological quantification of the total number of Pv+ neurons within the fifth lumbar DRG (n = 4 per genotype) was performed using the optical dissector method (StereoInvestigator, Microbrightfield). Every fourth section through the fifth lumbar DRG was analyzed. Contours containing a single DRG were optically sectioned using a 100x objective (NA 1.4) and an oil substage condenser. The sampling sites were 175 x 100 μm, and the counting frame boundaries were 80 x 60 x 6 μm in height. Only immunopositive neurons with sharply focused nuclei (containing at least one nucleolus) within the optical dissector counting frame boundaries were tallied. The total number of Pv+ neurons per DRG was estimated by the StereoInvestigator software.

    Statistical measures

    The Kolmogorov-Smirnov test was used to detect departures of the data from normality. F-tests were used to check for equal variances, and, in the case of unequal variances, the nonparametric Mann-Whitney U test was used. When variances were equal, the t test was used. The acceptance of significance was set to P < 0.05.

    Acknowledgments

    We thank the National Institute of Neurological Disorders and Stroke/National Institute of Mental Health Microarray Consortium (K. Ramsey) for performing the Affymetrix microarray analyses and the Northwestern University Transgenic and Targeted Mutagenesis Laboratory (L. Doglio and A. Alvarez) for pronuclear injections. The HSA promoter construct was provided by J. Chamberlain. The Egr3 NH2-terminal, Pv, and S46 antibodies were provided by J. Baraban, K. Baimbridge, and D. Fischman, respectively. K. Felio and M. McKnight provided technical assistance.

    This work was supported by the NIH (NS046468) and a Howard Hughes Faculty Scholar Award to W.G. Tourtellotte. L. Eldredge was supported by a predoctoral fellowship from the NIH (GM008061) and the NIH Medical Scientist Training Program (GM008152). J. Carter was supported by a predoctoral fellowship from the NIH (CA009560) and the NIH Medical Scientist Training Program (GM008152).

    References

    Andrechek, E.R., W.R. Hardy, A.A. Girgis-Gabardo, R.L. Perry, R. Butler, F.L. Graham, R.C. Kahn, M.A. Rudnicki, and W.J. Muller. 2002. ErbB2 is required for muscle spindle and myoblast cell survival. Mol. Cell. Biol. 22:4714–4722.

    Arber, S., D.R. Ladle, J.H. Lin, E. Frank, and T.M. Jessell. 2000. ETS gene Er81 controls the formation of functional connections between group Ia sensory afferents and motor neurons. Cell. 101:485–498.

    Chen, H.H., W.G. Tourtellotte, and E. Frank. 2002. Muscle spindle-derived neurotrophin 3 regulates synaptic connectivity between muscle sensory and motor neurons. J. Neurosci. 22:3512–3519.

    Chen, Y.W., P. Zhao, R. Borup, and E.P. Hoffman. 2000. Expression profiling in the muscular dystrophies: identification of novel aspects of molecular pathophysiology. J. Cell Biol. 151:1321–1336.

    Crawford, G.E., J.A. Faulkner, R.H. Crosbie, K.P. Campbell, S.C. Froehner, and J.S. Chamberlain. 2000. Assembly of the dystrophin-associated protein complex does not require the dystrophin COOH-terminal domain. J. Cell Biol. 150:1399–1410.

    Darby, I.A. 2000, ed. In situ hybridization protocols. Methods in Molecular Biology, vol. 123. Totowa, NJ: Humana Press.

    Dobretsov, M., S.L. Hastings, T.J. Sims, J.R. Stimers, and D. Romanovsky. 2003. Stretch receptor-associated expression of alpha 3 isoform of the Na+, K+-ATPase in rat peripheral nervous system. Neuroscience. 116:1069–1080.

    Ernfors, P., K.F. Lee, J. Kucera, and R. Jaenisch. 1994. Lack of neurotrophin-3 leads to deficiencies in the peripheral nervous system and loss of limb proprioceptive afferents. Cell. 77:503–512.

    Haase, G., E. Dessaud, A. Garces, B. de Bovis, M. Birling, P. Filippi, H. Schmalbruch, S. Arber, and O. deLapeyriere. 2002. GDNF acts through PEA3 to regulate cell body positioning and muscle innervation of specific motor neuron pools. Neuron. 35:893–905.

    He, T.C., S. Zhou, L.T. da Costa, J. Yu, K.W. Kinzler, and B. Vogelstein. 1998. A simplified system for generating recombinant adenoviruses. Proc. Natl. Acad. Sci. USA. 95:2509–2514.

    Hippenmeyer, S., N.A. Shneider, C. Birchmeier, S.J. Burden, T.M. Jessell, and S. Arber. 2002. A role for neuregulin1 signaling in muscle spindle differentiation. Neuron. 36:1035–1049.

    Jacobson, C., D. Duggan, and G. Fischbach. 2004. Neuregulin induces the expression of transcription factors and myosin heavy chains typical of muscle spindles in cultured human muscle. Proc. Natl. Acad. Sci. USA. 101:12218–12223.

    Klein, R., I. Silos-Santiago, R.J. Smeyne, S.A. Lira, R. Brambilla, S. Bryant, L. Zhang, W.D. Snider, and M. Barbacid. 1994. Disruption of the neurotrophin-3 receptor gene trkC eliminates la muscle afferents and results in abnormal movements. Nature. 368:249–251.

    Kucera, J., and J.M. Walro. 1992. Formation of muscle spindles in the absence of motor innervation. Neurosci. Lett. 145:47–50.

    Kucera, J., W. Cooney, A. Que, V. Szeder, H. Stancz-Szeder, and J. Walro. 2002. Formation of supernumerary muscle spindles at the expense of Golgi tendon organs in ER81-deficient mice. Dev. Dyn. 223:389–401.

    Leu, M., E. Bellmunt, M. Schwander, I. Farinas, H.R. Brenner, and U. Muller. 2003. Erbb2 regulates neuromuscular synapse formation and is essential for muscle spindle development. Development. 130:2291–2301.

    Livet, J., M. Sigrist, S. Stroebel, V. De Paola, S.R. Price, C.E. Henderson, T.M. Jessell, and S. Arber. 2002. ETS gene Pea3 controls the central position and terminal arborization of specific motor neuron pools. Neuron. 35:877–892.

    Muscat, G.E., S. Perry, H. Prentice, and L. Kedes. 1992. The human skeletal alpha-actin gene is regulated by a muscle-specific enhancer that binds three nuclear factors. Gene Expr. 2:111–126.

    O'Donovan, K.J., E.P. Wilkens, and J.M. Baraban. 1998. Sequential expression of Egr-1 and Egr-3 in hippocampal granule cells following electroconvulsive stimulation. J. Neurochem. 70:1241–1248.

    Rando, T.A., and H.M. Blau. 1994. Primary mouse myoblast purification, characterization, and transplantation for cell-mediated gene therapy. J. Cell Biol. 125:1275–1287.

    Swirnoff, A.H., and J. Milbrandt. 1995. DNA-binding specificity of NGFI-A and related zinc finger transcription factors. Mol. Cell. Biol. 15:2275–2287.

    Tessarollo, L., K.S. Vogel, M.E. Palko, S.W. Reid, and L.F. Parada. 1994. Targeted mutation in the neurotrophin-3 gene results in loss of muscle sensory neurons. Proc. Natl. Acad. Sci. USA. 91:11844–11848.

    Tourtellotte, W.G., and J. Milbrandt. 1998. Sensory ataxia and muscle spindle agenesis in mice lacking the transcription factor Egr3. Nat. Genet. 20:87–91.

    Tourtellotte, W.G., C. Keller-Peck, J. Milbrandt, and J. Kucera. 2001. The transcription factor Egr3 modulates sensory axon-myotube interactions during muscle spindle morphogenesis. Dev. Biol. 232:388–399.

    Whitehead, J., C. Keller-Peck, J. Kucera, and W.G. Tourtellotte. 2005. Glial cell-line derived neurotrophic factor-dependent fusimotor neuron survival during development. Mech. Dev. 122:27–41.

    Zelena, J. 1957. The morphogenic influences of innervation on the ontogenic development of muscle spindles. J. Embryol. Exp. Morphol. 5:283–292.

    Zelena, J., and T. Soukup. 1973. Development of muscle spindles deprived of fusimotor(Y'vonne Albert, Jennifer )