Generation of the primary hair follicle pattern.

Hair follicles are spaced apart from one another at regular intervals through the skin. Although follicles are predominantly epidermal structures, classical tissue recombination experiments indicated that the underlying dermis defines their location during development. Although many molecules involved in hair follicle formation have been identified, the molecular interactions that determine the emergent property of pattern formation have remained elusive. We have used embryonic skin cultures to dissect signaling responses and patterning outcomes as the skin spatially organizes itself. We find that ectodysplasin receptor (Edar)-bone morphogenetic protein (BMP) signaling and transcriptional interactions are central to generation of the primary hair follicle pattern, with restriction of responsiveness, rather than localization of an inducing ligand, being the key driver in this process. The crux of this patterning mechanism is rapid Edar-positive feedback in the epidermis coupled with induction of dermal BMP4/7. The BMPs in turn repress epidermal Edar and hence follicle fate. Edar activation also induces connective tissue growth factor, an inhibitor of BMP signaling, allowing BMP action only at a distance from their site of synthesis. Consistent with this model, transgenic hyperactivation of Edar signaling leads to widespread overproduction of hair follicles. This Edar-BMP activation-inhibition mechanism appears to operate alongside a labile prepattern, suggesting that Edar-mediated stabilization of beta-catenin active foci is a key event in determining definitive follicle locations.

A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion

An epidemic model is formulated by a reactionâeuro"diffusion system where the spatial pattern formation is driven by cross-diffusion. The reaction terms describe the local dynamics of susceptible and infected species, whereas the diffusion terms account for the spatial distribution dynamics. For both self-diffusion and cross-diffusion, nonlinear constitutive assumptions are suggested. To simulate the pattern formation two finite volume formulations are proposed, which employ a conservative and a non-conservative discretization, respectively. An efficient simulation is obtained by a fully adaptive multiresolution strategy. Numerical examples illustrate the impact of the cross-diffusion on the pattern formation.

Robustness of Drosophila early embryo and wing imaginal disc development

Within a developing organism, cells require information on where they are in order to differentiate into the correct cell-type. Pattern formation is the process by which cells acquire and process positional cues and thus determine their fate. This can be achieved by the production and release of a diffusible signaling molecule, called a morphogen, which forms a concentration gradient: exposure to different morphogen levels leads to the activation of specific signaling pathways. Thus, in response to the morphogen gradient, cells start to express different sets of genes, forming domains characterized by a unique combination of differentially expressed genes. As a result, a pattern of cell fates and specification emerges.Though morphogens have been known for decades, it is not yet clear how these gradients form and are interpreted in order to yield highly robust patterns of gene expression. During my PhD thesis, I investigated the properties of Bicoid (Bcd) and Decapentaplegic (Dpp), two morphogens involved in the patterning of the anterior-posterior axis of Drosophila embryo and wing primordium, respectively. In particular, I have been interested in understanding how the pattern proportions are maintained across embryos of different sizes or within a growing tissue. This property is commonly referred to as scaling and is essential for yielding functional organs or organisms. In order to tackle these questions, I analysed fluorescence images showing the pattern of gene expression domains in the early embryo and wing imaginal disc. After characterizing the extent of these domains in a quantitative and systematic manner, I introduced and applied a new scaling measure in order to assess how well proportions are maintained. I found that scaling emerged as a universal property both in early embryos (at least far away from the Bcd source) and in wing imaginal discs (across different developmental stages). Since we were also interested in understanding the mechanisms underlying scaling and how it is transmitted from the morphogen to the target genes down in the signaling cascade, I also quantified scaling in mutant flies where this property could be disrupted. While scaling is largely conserved in embryos with altered bcd dosage, my modeling suggests that Bcd trapping by the nuclei as well as pre-steady state decoding of the morphogen gradient are essential to ensure precise and scaled patterning of the Bcd signaling cascade. In the wing imaginal disc, it appears that as the disc grows, the Dpp response expands and scales with the tissue size. Interestingly, scaling is not perfect at all positions in the field. The scaling of the target gene domains is best where they have a function; Spalt, for example, scales best at the position in the anterior compartment where it helps to form one of the anterior veins of the wing. Analysis of mutants for pentagone, a transcriptional target of Dpp that encodes a secreted feedback regulator of the pathway, indicates that Pentagone plays a key role in scaling the Dpp gradient activity.

Analysis of a finite volume method for a cross-diffusion model in population dynamics

The main goal of this paper is to propose a convergent finite volume method for a reactionâeuro"diffusion system with cross-diffusion. First, we sketch an existence proof for a class of cross-diffusion systems. Then the standard two-point finite volume fluxes are used in combination with a nonlinear positivity-preserving approximation of the cross-diffusion coefficients. Existence and uniqueness of the approximate solution are addressed, and it is also shown that the scheme converges to the corresponding weak solution for the studied model. Furthermore, we provide a stability analysis to study pattern-formation phenomena, and we perform two-dimensional numerical examples which exhibit formation of nonuniform spatial patterns. From the simulations it is also found that experimental rates of convergence are slightly below second order. The convergence proof uses two ingredients of interest for various applications, namely the discrete Sobolev embedding inequalities with general boundary conditions and a space-time $L^1$ compactness argument that mimics the compactness lemma due to Kruzhkov. The proofs of these results are given in the Appendix.

Secondary effects of catalytic diesel particulate filters: copper-induced formation of PCDD/Fs.

Potential risks of a secondary formation of polychlorinated dibenzodioxins/furans (PCDD/Fs) were assessed for two cordierite-based, wall-through diesel particulate filters (DPFs) for which soot combustion was either catalyzed with an iron- or a copper-based fuel additive. A heavy duty diesel engine was used as test platform, applying the eight-stage ISO 8178/4 C1 cycle. DPF applications neither affected the engine performance, nor did they increase NO, NO2, CO, and CO2 emissions. The latter is a metric for fuel consumption. THC emissions decreased by about 40% when deploying DPFs. PCDD/F emissions, with a focus on tetra- to octachlorinated congeners, were compared under standard and worst case conditions (enhanced chlorine uptake). The iron-catalyzed DPF neither increased PCDD/F emissions, nor did it change the congener pattern, even when traces of chlorine became available. In case of copper, PCDD/F emissions increased by up to 3 orders of magnitude from 22 to 200 to 12 700 pg I-TEQ/L with fuels of < 2, 14, and 110 microg/g chlorine, respectively. Mainly lower chlorinated DD/Fs were formed. Based on these substantial effects on PCDD/F emissions, the copper-catalyzed DPF system was not approved for workplace applications, whereas the iron system fulfilled all the specifications of the Swiss procedures for DPF approval (VERT).

Mechanotransduction, PROX1, and FOXC2 cooperate to control connexin37 and calcineurin during lymphatic-valve formation.

Lymphatic valves are essential for efficient lymphatic transport, but the mechanisms of early lymphatic-valve morphogenesis and the role of biomechanical forces are not well understood. We found that the transcription factors PROX1 and FOXC2, highly expressed from the onset of valve formation, mediate segregation of lymphatic-valve-forming cells and cell mechanosensory responses to shear stress in vitro. Mechanistically, PROX1, FOXC2, and flow coordinately control expression of the gap junction protein connexin37 and activation of calcineurin/NFAT signaling. Connexin37 and calcineurin are required for the assembly and delimitation of lymphatic valve territory during development and for its postnatal maintenance. We propose a model in which regionally increased levels/activation states of transcription factors cooperate with mechanotransduction to induce a discrete cell-signaling pattern and morphogenetic event, such as formation of lymphatic valves. Our results also provide molecular insights into the role of endothelial cell identity in the regulation of vascular mechanotransduction.

Innate signaling promotes formation of regulatory nitric oxide-producing dendritic cells limiting T-cell expansion in experimental autoimmune myocarditis.

BACKGROUND: Activation of innate pattern-recognition receptors promotes CD4+ T-cell-mediated autoimmune myocarditis and subsequent inflammatory cardiomyopathy. Mechanisms that counterregulate exaggerated heart-specific autoimmunity are poorly understood. METHODS AND RESULTS: Experimental autoimmune myocarditis was induced in BALB/c mice by immunization with α-myosin heavy chain peptide and complete Freund's adjuvant. Together with interferon-γ, heat-killed Mycobacterium tuberculosis, an essential component of complete Freund's adjuvant, converted CD11b(hi)CD11c(-) monocytes into tumor necrosis factor-α- and nitric oxide synthase 2-producing dendritic cells (TipDCs). Heat-killed M. tuberculosis stimulated production of nitric oxide synthase 2 via Toll-like receptor 2-mediated nuclear factor-κB activation. TipDCs limited antigen-specific T-cell expansion through nitric oxide synthase 2-dependent nitric oxide production. Moreover, they promoted nitric oxide synthase 2 production in hematopoietic and stromal cells in a paracrine manner. Consequently, nitric oxide synthase 2 production by both radiosensitive hematopoietic and radioresistant stromal cells prevented exacerbation of autoimmune myocarditis in vivo. CONCLUSIONS: Innate Toll-like receptor 2 stimulation promotes formation of regulatory TipDCs, which confine autoreactive T-cell responses in experimental autoimmune myocarditis via nitric oxide. Therefore, activation of innate pattern-recognition receptors is critical not only for disease induction but also for counterregulatory mechanisms, protecting the heart from exaggerated autoimmunity.

Formation of a spatial memory trace : a behavioral, cellular and molecular study

THESIS ABSTRACTThis thesis project was aimed at studying the molecular mechanisms underlying learning and memory formation, in particular as they relate to the metabolic coupling between astrocytes and neurons. For that, changes in the metabolic activity of different mice brain regions after 1 or 9 days of training in an eight-arm radial maze were assessed by (14C) 2-deoxyglucose (2DG) autoradiography. Significant differences in the areas engaged during the behavioral task at day 1 (when animals are confronted for the first time to the learning task) and at day 9 (when animals are highly performing) have been identified. These areas include the hippocampus, the fornix, the parietal cortex, the laterodorsal thalamic nucleus and the mammillary bodies at day 1 ; and the anterior cingulate, the retrosplenial cortex and the dorsal striatum at day 9. Two of these cerebral regions (those presenting the greatest changes at day 1 and day 9: the hippocampus and the retrosplenial cortex, respectively) were microdissected by laser capture microscopy and selected genes related to neuron-glia metabolic coupling, glucose metabolism and synaptic plasticity were analyzed by RT-PCR. 2DG and gene expression analysis were performed at three different times: 1) immediately after the end of the behavioral paradigm, 2) 45 minutes and 3) 6 hours after training. The main goal of this study was the identification of the metabolic adaptations following the learning task. Gene expression results demonstrate that the learning task profoundly modulates the pattern of gene expression in time, meaning that these two cerebral regions with high 2DG signal (hippocampus and retrosplenial cortex) have adapted their metabolic molecular machinery in consequence. Almost all studied genes show a higher expression in the hippocampus at day 1 compared to day 9, while an increased expression was found in the retrosplenial cortex at day 9. We can observe these molecular adaptations with a short delay of 45 minutes after the end of the task. However, 6 hours after training a high gene expression was found at day 9 (compared to day 1) in both regions, suggesting that only one day of training is not sufficient to detect transcriptional modifications several hours after the task. Thus, gene expression data match 2DG results indicating a transfer of information in time (from day 1 to day 9) and in space (from the hippocampus to the retrosplenial cortex), and this at a cellular and a molecular level. Moreover, learning seems to modify the neuron-glia metabolic coupling, since several genes involved in this coupling are induced. These results also suggest a role of glia in neuronal plasticity.RESUME DU TRAVAIL DE THESECe projet de thèse a eu pour but l'étude des mécanismes moléculaires qui sont impliqués dans l'apprentissage et la mémoire et, en particulier, à les mettre en rapport avec le couplage métabolique existant entre les astrocytes et les neurones. Pour cela, des changements de l'activité métabolique dans différentes régions du cerveau des souris après 1 ou 9 jours d'entraînement dans un labyrinthe radial à huit-bras ont été évalués par autoradiographie au 2-désoxyglucose (2DG). Des différences significatives dans les régions engagées pendant la tâche comportementale au jour 1 (quand les animaux sont confrontés pour la première fois à la tâche) et au jour 9 (quand les animaux ont déjà appris) ont été identifiés. Ces régions incluent, au jour 1, l'hippocampe, le fornix, le cortex pariétal, le noyau thalamic laterodorsal et les corps mamillaires; et, au jour 9, le cingulaire antérieur, le cortex retrosplenial et le striatum dorsal. Deux de ces régions cérébrales (celles présentant les plus grands changements à jour 1 et à jour 9: l'hippocampe et le cortex retrosplenial, respectivement) ont été découpées par microdissection au laser et quelques gènes liés au couplage métabolique neurone-glie, au métabolisme du glucose et à la plasticité synaptique ont été analysées par RT-PCR. L'étude 2DG et l'analyse de l'expression de gènes ont été exécutés à trois temps différents: 1) juste après entraînement, 2) 45 minutes et 3) 6 heures après la fin de la tâche. L'objectif principal de cette étude était l'identification des adaptations métaboliques suivant la tâche d'apprentissage. Les résultats de l'expression de gènes démontrent que la tâche d'apprentissage module profondément le profile d'expression des gènes dans le temps, signifiant que ces deux régions cérébrales avec un signal 2DG élevé (l'hippocampe et le cortex retrosplenial) ont adapté leurs « machines moléculaires » en conséquence. Presque tous les gènes étudiés montrent une expression plus élevée dans l'hippocampe au jour 1 comparé au jour 9, alors qu'une expression accrue a été trouvée dans le cortex retrosplenial au jour 9. Nous pouvons observer ces adaptations moléculaires avec un retard court de 45 minutes après la fin de la tâche. Cependant, 6 heures après l'entraînement, une expression de gènes élevée a été trouvée au jour 9 (comparé à jour 1) dans les deux régions, suggérant que seulement un jour d'entraînement ne suffit pas pour détecter des modifications transcriptionelles plusieurs heures après la tâche. Ainsi, les données d'expression de gènes corroborent les résultats 2DG indiquant un transfert d'information dans le temps (de jour 1 à jour 9) et dans l'espace (de l'hippocampe au cortex retrosplenial), et ceci à un niveau cellulaire et moléculaire. D'ailleurs, la tâche d'apprentissage semble modifier le couplage métabolique neurone-glie, puisque de nombreux gènes impliqués dans ce couplage sont induits. Ces observations suggèrent un rôle important de la glie dans les mécanismes de plasticité du système nerveux.

Edar/Eda interactions regulate enamel knot formation in tooth morphogenesis.

tabby and downless mutant mice have apparently identical defects in teeth, hair and sweat glands. Recently, genes responsible for these spontaneous mutations have been identified. downless (Dl) encodes Edar, a novel member of the tumour necrosis factor (TNF) receptor family, containing the characteristic extracellular cysteine rich fold, a single transmembrane region and a death homology domain close to the C terminus. tabby (Ta) encodes ectodysplasin-A (Eda) a type II membrane protein of the TNF ligand family containing an internal collagen-like domain. As predicted by the similarity in adult mutant phenotype and the structure of the proteins, we demonstrate that Eda and Edar specifically interact in vitro. We have compared the expression pattern of Dl and Ta in mouse development, taking the tooth as our model system, and find that they are not expressed in adjacent cells as would have been expected. Teeth develop by a well recorded series of epithelial-mesenchymal interactions, similar to those in hair follicle and sweat gland development, the structures found to be defective in tabby and downless mice. We have analysed the downless mutant teeth in detail, and have traced the defect in cusp morphology back to initial defects in the structure of the tooth enamel knot at E13. Significantly, the defect is distinct from that of the tabby mutant. In the tabby mutant, there is a recognisable but small enamel knot, whereas in the downless mutant the knot is absent, but enamel knot cells are organised into a different shape, the enamel rope, showing altered expression of signalling factors (Shh, Fgf4, Bmp4 and Wnt10b). By adding a soluble form of Edar to tooth germs, we were able to mimic the tabby enamel knot phenotype, demonstrating the involvement of endogenous Eda in tooth development. We could not, however, reproduce the downless phenotype, suggesting the existence of yet another ligand or receptor, or of ligand-independent activation mechanisms for Edar. Changes in the structure of the enamel knot signalling centre in downless tooth germs provide functional data directly linking the enamel knot with tooth cusp morphogenesis. We also show that the Lef1 pathway, thought to be involved in these mutants, functions independently in a parallel pathway.

Spatial pattern and process at the plant neighbourhood scale: insights from communities dominated by the clonal grass Elymus repens (L.) Gould.

Question: How do clonal traits of a locally dominant grass (Elymus repens (L.) Gould.) respond to soil heterogeneity and shape spatial patterns of its tillers? How do tiller spatial patterns constrain seedling recruitment within the community?Locations: Artificial banks of the River Rhone, France.Material and Methods: We examined 45 vegetation patches dominated by Elymus repens. During a first phase we tested relationships between soil variables and three clonal traits (spacer length, number of clumping tillers and branching rate), and between the same clonal traits and spatial patterns (i.e. density and degree of spatial aggregation) of tillers at a very fine scale. During a second phase, we performed a sowing experiment to investigate effects of density and spatial patterns of E. repens on recruitment of eight species selected from the regional species pool.Results: Clonal traits had clear effects - especially spacer length - on densification and aggregation of E. repens tillers and, at the same time, a clear response of these same clonal traits as soil granulometry changed. The density and degree of aggregation of E. repens tillers was positively correlated to total seedling cover and diversity at the finest spatial scales.Conclusions: Spatial patterning of a dominant perennial grass responds to soil heterogeneity through modifications of its clonal morphology as a trade-off between phalanx and guerrilla forms. In turn, spatial patterns have strong effects on abundance and diversity of seedlings. Spatial patterns of tillers most probably led to formation of endogenous gaps in which the recruitment of new plant individuals was enhanced. Interestingly, we also observed more idiosyncratic effects of tiller spatial patterns on seedling cover and diversity when focusing on different growth forms of the sown species.

Citrobacter rodentium Requires Intestinal Neural Wiskott-Aldrich Syndrome Protein (N-WASp) for Actin Pedestal Formation, Colonization and Epithelial Barrier Disruption In Vivo

Background: Citrobacter rodentium is a natural mouse pathogen that is genetically closelyrelated to the human enteric pathogens enteropathogenic and enterohemorrhagic E. coli.Among the repertoire of conserved virulence factors that these pathogens deliver via typeIII secretion, Tir and EspF are responsible for the formation of characteristic actin-richpedestals and disruption of tight junction integrity, respectively. There is evidence In Vitrothese effectors accomplish this, at least in part, by subverting the normal host cellularfunctions of N-WASP, a critical regulator of branched chain actin assembly. Although NWASPhas been shown to be involved in pedestal formation In Vitro, the requirements ofN-WASP-mediated actin pedestals for intestinal colonization by attaching/effacing (A/E)pathogens In Vivo is not known. Furthermore, it is not known whether N-WASP is requiredfor EspF-mediated tight junction disruption. Methods: To investigate the role of N-WASPin the gut epithelium, we generated mice with intestine-specific deletion of N-WASP(iNWKO), by mating mice homozygous for a floxed N-WASP allele (N-WASPL2L/L2L) tomice expressing Cre recombinase under the villin promoter. Separately housed groups ofWT and iNWKO mice were inoculated with 5x108 GFP-expressing C. rodentium by intragastriclavage. Stool was collected 2, 4, 7, and 12 days after infection, and recoverablecolony forming units (CFUs) of C. rodentium were quantified by plating serial dilutions ofhomogenized stool on MacConkey's agar. GFP+ colonies were counted after 24 hoursincubation at 37°C. The presence of actin pedestals was investigated by electron microscopy(EM), and tight junction morphology was assessed by immunofluorescence staining ofoccludin, ZO-1 and claudin-2. Results: C. rodentium infection did not result in mortalityin WT or iNWKO mice. Compared to controls, iNWKO mice exhibited higher levels ofbacterial shedding during the first 4 days of infection (day 4 average: WT 5.2x104 CFU/gvs. iNWKO 4.7x105 CFU/g, p=0.08), followed by a more rapid clearance of C. rodentium, (day7-12 average: WT 2x106 CFU/g vs. iNWKO 2.7x105, p=0.01). EM and immunofluorescencerevealed the complete lack of actin pedestals in iNWKO mice and no mucosa-associatedGFP+ C. rodentium by day 7. WT controls exhibited tight junction disruption, reflected byaltered distribution of ZO-1, whereas iNWKO mice had no change in the pattern of ZO-1.Conclusion: Intestinal N-WASP is required for actin pedestal formation by C. rodentium InVivo, and ablation of N-WASP is associated with more rapid bacterial clearance and decreasedability of C. rodentium to disrupt intercellular junctions.

The primary in vivo immune response to Mls-1 (Mtv-7 sag). Route of injection determines the immune response pattern.

Injection of cells expressing the retroviral superantigen Mls-1 (Mtv-7 sag) into adult Mls-1- mice induces a strong immune response including both T- and B-cell activation. This model was used for studying qualitative aspects of the immune response in normal mice with a defined antigen-presenting cell (the B cell) and without the use of adjuvant. BALB/c mice were injected locally or systemically with Mls-1-expressing spleen cells from Mls-1-congenic BALB.D2 mice. Intravenous injection led to an initially strong expansion of Mls-1-reactive V beta 6+ CD4+ cells mainly in the spleen, to a large degree explained by the trapping of reactive cells, and a rapid down-regulation of interleukin-2 (IL-2) and interferon-gamma (IFN-gamma) production, consistent with the proposed tolerogenic property of B cells as antigen-presenting cells. However, these mice developed a slowly appearing but persistent B-cell response dominated by IgG1-producing cells, suggesting a shift in lymphokines produced rather than complete unresponsiveness. Subcutaneous injection into the hind footpad with the same number of cells led to a strong local response in the draining lymph node, characterized by a dramatic increase of V beta 6+ CD4+ T cells, local production of IL-2 and IFN-gamma and a strong but short-lived antibody response dominated by IgG2a-producing cells, characteristic of a T-helper type 1 (Th1) type of response. Both routes of injection led ultimately to deletion of reactive T cells and anergy, as defined by the inability to produce IL-2 upon in vitro stimulation with Mls-1. It is concluded that Mls-1 presented by B cells induces qualitatively different responses in vivo dependent on the route of injection. We propose that the different responses result from the migration of the injected cells to different micro-anatomical sites in the lymphoid tissue. Furthermore, these results suggest that B cells may function as professional antigen-presenting cells in vivo present in an appropriate environment.

Learning-induced gene expression in the heads of two Nasonia species that differ in long-term memory formation.

BACKGROUND: Cellular processes underlying memory formation are evolutionary conserved, but natural variation in memory dynamics between animal species or populations is common. The genetic basis of this fascinating phenomenon is poorly understood. Closely related species of Nasonia parasitic wasps differ in long-term memory (LTM) formation: N. vitripennis will form transcription-dependent LTM after a single conditioning trial, whereas the closely-related species N. giraulti will not. Genes that were differentially expressed (DE) after conditioning in N. vitripennis, but not in N. giraulti, were identified as candidate genes that may regulate LTM formation. RESULTS: RNA was collected from heads of both species before and immediately, 4 or 24 hours after conditioning, with 3 replicates per time point. It was sequenced strand-specifically, which allows distinguishing sense from antisense transcripts and improves the quality of expression analyses. We determined conditioning-induced DE compared to naïve controls for both species. These expression patterns were then analysed with GO enrichment analyses for each species and time point, which demonstrated an enrichment of signalling-related genes immediately after conditioning in N. vitripennis only. Analyses of known LTM genes and genes with an opposing expression pattern between the two species revealed additional candidate genes for the difference in LTM formation. These include genes from various signalling cascades, including several members of the Ras and PI3 kinase signalling pathways, and glutamate receptors. Interestingly, several other known LTM genes were exclusively differentially expressed in N. giraulti, which may indicate an LTM-inhibitory mechanism. Among the DE transcripts were also antisense transcripts. Furthermore, antisense transcripts aligning to a number of known memory genes were detected, which may have a role in regulating these genes. CONCLUSION: This study is the first to describe and compare expression patterns of both protein-coding and antisense transcripts, at different time points after conditioning, of two closely related animal species that differ in LTM formation. Several candidate genes that may regulate differences in LTM have been identified. This transcriptome analysis is a valuable resource for future in-depth studies to elucidate the role of candidate genes and antisense transcription in natural variation in LTM formation.

Spatiotemporal neuromodulation therapies engaging muscle synergies improve motor control after spinal cord injury.

Electrical neuromodulation of lumbar segments improves motor control after spinal cord injury in animal models and humans. However, the physiological principles underlying the effect of this intervention remain poorly understood, which has limited the therapeutic approach to continuous stimulation applied to restricted spinal cord locations. Here we developed stimulation protocols that reproduce the natural dynamics of motoneuron activation during locomotion. For this, we computed the spatiotemporal activation pattern of muscle synergies during locomotion in healthy rats. Computer simulations identified optimal electrode locations to target each synergy through the recruitment of proprioceptive feedback circuits. This framework steered the design of spatially selective spinal implants and real-time control software that modulate extensor and flexor synergies with precise temporal resolution. Spatiotemporal neuromodulation therapies improved gait quality, weight-bearing capacity, endurance and skilled locomotion in several rodent models of spinal cord injury. These new concepts are directly translatable to strategies to improve motor control in humans.