# Early movement restriction leads to maladaptive plasticity in the sensorimotor cortex and to movement disorders

Maxime S Delcour, Michaël Russier, Francis Castets, Nathalie Turle-Lorenzo, Marie-Hélène Canu, Florence Cayetanot, Mary F Barbe, Jacques-Olivier Coq

► To cite this version:

Maxime S Delcour, Michaël Russier, Francis Castets, Nathalie Turle-Lorenzo, Marie-Hélène Canu, et al.. Early movement restriction leads to maladaptive plasticity in the sensorimotor cortex and to movement disorders. Scientific Reports, 2018, 8 (1), 10.1038/s41598-018-34312-y. hal-01943923

HAL Id: hal-01943923

https://hal.sorbonne-universite.fr/hal-01943923

Submitted on 4 Dec 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

Received: 23 March 2018

Accepted: 16 October 2018

Published online: 05 November 2018

# Early movement restriction leads to maladaptive plasticity in the sensorimotor cortex and to movement disorders

Maxime Delcour<fnref n="1" />,⁷, Michaël Russier<fnref n="1" />,⁸, Francis Castets<fnref n="2" />, Nathalie Turle-Lorenzo³, Marie-Hélène Canu ID<fnref n="4" />, Florence Cayetanot⁵,⁹, Mary F Barbe⁶ & Jacques-Olivier Coq ID <fnref n="1" />,⁵

Motor control and body representations in the central nervous system are built, i.e., patterned, during development by sensorimotor experience and somatosensory feedback/reafference. Yet, early emergence of locomotor disorders remains a matter of debate, especially in the absence of brain damage. For instance, children with developmental coordination disorders (DCD) display deficits in planning, executing and controlling movements, concomitant with deficits in executive functions. Thus, are early sensorimotor atypicalities at the origin of long-lasting abnormal development of brain anatomy and functions? We hypothesize that degraded locomotor outcomes in adulthood originate as a consequence of early atypical sensorimotor experiences that induce developmental disorganization of sensorimotor circuitry. We showed recently that postnatal sensorimotor restriction (SMR), through hind limb immobilization from birth to one month, led to enduring digitigrade locomotion with ankle-knee overextension, degraded musculoskeletal tissues (e.g., gastrocnemius atrophy), and clear signs of spinal hyperreflexia in adult rats, suggestive of spasticity; each individual disorder likely interplaying in self-perpetuating cycles. In the present study, we investigated the impact of postnatal SMR on the anatomical and functional organization of hind limb representations in the sensorimotor cortex and processes representative of maladaptive neuroplasticity. We found that 28 days of daily SMR degraded the topographical organization of somatosensory hind limb maps, reduced both somatosensory and motor map areas devoted to the hind limb representation and altered neuronal response properties in the sensorimotor cortex several weeks after the cessation of SMR. We found no neuroanatomical histopathology in hind limb sensorimotor cortex, yet increased glutamatergic neurotransmission that matched clear signs of spasticity and hyperexcitability in the adult lumbar spinal network. Thus, even in the absence of a brain insult, movement disorders and brain dysfunction can emerge as a consequence of reduced and atypical patterns of motor outputs and somatosensory feedback that induce maladaptive neuroplasticity. Our results may contribute to understanding the inception and mechanisms underlying neurodevelopmental disorders, such as DCD.

<footnote n="1">¹Neurosciences Intégratives et Adaptatives, UMR 7260, CNRS, Aix-Marseille Université, 13331, Marseille, France.</footnote>

<footnote n="2">²Centre de Recherche en Neurobiologie et Neurophysiologie de Marseille UMR 7286, CNRS, Aix-Marseille Université, 13344, Marseille, France. ³FR 3512 Fédération 3C, Aix Marseille Université - CNRS, 13331, Marseille, France.</footnote>

<footnote n="4">⁴Université de Lille, EA 7369 « Activité Physique, Muscle et Santé » - URePSSS - Unité de Recherche Pluridisciplinaire Sport Santé Société, 59000, Lille, France. ⁵Institut de Neurosciences de la Timone, UMR 7289, CNRS, Aix-Marseille Université, 13385, Marseille, France. ⁶Department of Anatomy and Cell Biology, Lewis Katz School of Medicine, Temple University, Philadelphia, PA, 19140, USA. ⁷Present address: Equipe de Recherche en Réadaptation Sensorimotrice, Faculté de Médecine, Département de Physiologie, Université de Montréal, C.P. 6128, Montréal, H3C 3J7, Canada. ⁸Present address: Inserm UMR 1072, Unité de Neurobiologie des Canaux Ioniques et de la Synapse, Faculté de Médecine Secteur Nord, 13344, Marseille Cedex 15, France. ⁹Present address: UMR_S1158 Inserm-Sorbonne Université, Neurophysiologie Respiratoire Expérimentale et Clinique, Faculté de Médecine, 75636, Paris Cedex, France. Correspondence and requests for materials should be addressed to J.-O.C. (email: jacques-olivier.coq@univ-amu.fr)</footnote>
It is now understood that development of movement repertoires, motor control and body representations in sensorimotor circuitry are achieved through early spontaneous movements, sensorimotor experiences and reaf- ference in children¹⁻³ and rodents⁴,⁵. From this, it seems likely that early atypical sensorimotor experiences in children should lead to the emergence of atypical movements, motor control problems and disorganization of sen- sorimotor circuitry that persist into adulthood. The impact of disuse or immobilization (i.e., constraint-induced movement therapy) during child development has mainly been studied in the presence of brain lesions⁶. Only a very few studies have focused on the impact of early disuse occurring in the absence of brain damage. For exam- ple, swaddling, the ancient practice of wrapping infants in cloth so that limb movements are restricted, appears to delay the onset of several motor skills⁷,⁸. Swaddling has become popular again, such as in neonatal intensive care units using elastic cloth, as a means to reduce sudden infant death syndrome or crying, promote sleep, or improve muscle tone⁸,⁹. Disturbances in the planning, execution and control of body movements in the absence of brain damage are now termed developmental coordination disorders (DCD) and usually coexist with various deficits in executive functions in 5-6% of school-aged children¹⁰. Patients with DCD show reduced abilities to produce con- sistent movements¹¹, poor motor coordination and kinaesthetic acuity¹²,¹³, broad impairments in sensorimotor representations and perception¹³,¹⁴, each reflected by disrupted central networks¹⁵ Children with autism spectrum disorder (ASD) also exhibit gross or fine motor abnormalities, motor learning deficiencies and difficulties exe- cuting sequences of actions¹⁶. Most children with DCD or ASD show similar sensorimotor impairments, reduced physical activity and interactions with their environment, and atypical motor development, the latter detected as atypical spontaneous or general movements (GMs)¹⁷,¹⁸.

During typical development, the repertoire of GMs in limbs increases over time in variation, fluency, ampli-
tude and complexity into a continuous stream of small and elegant movements. Atypical GMs correspond to
rigid, cramped synchronized and stereotyped movements that exhibit limited fluency, variation and complex-
ity with increasing age³. Arising from spontaneous, self-generated and evoked movements during maturation,
early somatosensory feedback drives electrical activity patterning from the spinal cord to the cortex. This feed-
back guides the development and refinement of the anatomical and functional organization of sensorimotor
circuitry in rodents⁵,¹⁹,²⁰. Atypical, disturbed GMs reflect impaired connectivity and functional disorganization
in the brain²¹,²². Accordingly, we hypothesized that limited and abnormal patterns of somatosensory inputs
during development may lead to abnormal anatomical and functional organization of the sensorimotor circuitry
in adulthood as a result of maladaptive cortical plasticity. Such sensorimotor disorganization may in turn alter
somatosensory and bodily perceptions, motor outputs, and musculoskeletal structure and physiology.

We recently showed that transient postnatal sensorimotor restriction (SMR), experimentally produced using 28 days of transient (16 hours/day) hind limb immobilization of developing rats, lead to degraded locomotion on a treadmill that persisted for more than 30 days after cessation of the immobilization. This long-lasting deg- radation was characterized by reduced length, amplitude and velocity, overextended knees and ankles, and dig- itigrade locomotion that resembled true pes equinus or “toe walking”²³. SMR rats displayed not only increased variations in the kinematic parameters of treadmill locomotion, mentioned above, but also reduced variations of hind limb joint angles; these variations persisted over time²³. These results appear to recapitulate the “toe walk- ing” and increased spatiotemporal variations of movements observed in children with ASD²⁴,²⁵. SMR also led to musculoskeletal histopathology and increased stretch reflexes, suggestive of muscle hyperreflexia and spasticity. We postulated that reduced and atypical patterns of both motor outputs and somatosensory reafference during development likely contributed to the emergence of movement disorders and musculoskeletal pathologies that persisted into adulthood²³. To further explore this concept that adult movement disorders may have develop- mental origins even in the absence of brain damage, we investigated here the impact of postnatal SMR on the neuroanatomical and functional organization of hind limb somatosensory and motor representations in adult rats using microelectrode cortical mapping techniques, and brain histology and immunohis- tochemistry studies (SMR, n = 13; Cont, n = 18) from P90 to P120, i.e., about 60 to 90 days after cessation of the SMR. Since we did not provide anatomical or histological features of the region of cortex in which we recorded or stimulated, we will refer to somatosensory and motor maps originating from the cortex, as opposed to primary somatosensory and motor cortices²⁶.

Results

We used 17 rats that were exposed to sensorimotor restriction (SMR; Fig. 1) and 20 control rats of either sex. The rats that underwent gait testing at P30 and P65 were also used for hind limb musculoskeletal assessments (results of which have been previously published²³), as well as electrophysiological procedures (somatosensory maps: SMR, n = 9; Cont, n = 10; neural motor maps: SMR, n = 8; Cont, n = 9), in vivo microdialysis (SMR, n = 8; Cont, n = 6) and Western-blotting (SMR, n = 9; Cont, n = 6), and stereological brain histology and immunohistochemistry studies (SMR, n = 13; Cont, n = 18) from P90 to P120, i.e., about 60 to 90 days after cessation of the SMR. Since we did not provide anatomical or histological features of the region of cortex in which we recorded or stimulated, we will refer to somatosensory and motor maps originating from the cortex, as opposed to primary somatosensory and motor cortices²⁶.

**SMR degrades the somatosensory map organization and neuronal properties.** As shown in previous studies²⁷,²⁸, the somatosensory hind paw representation originating from the cortex displayed invariant organizational features despite inter-individual differences, called somatotopy. Briefly, the hind paw representation was usually located medial to the forepaw map and rostral to the tail and back/ventral representations. From rostral to caudal, the somatotopic organization corresponded to the progression of cortical sites from the toes, plantar pads of the sole to the heel, and leg (Fig. 2A,B). From lateral to medial in the rostral portion of the hind paw map, the toes were topographically represented from toe 1 to toe 5 (Fig. 2A–C). The hairy representation of
<figure><img src="image_2.png" /><figcaption>**Figure 1.** A young rat submitted to postnatal hind limb immobilization leading to sensorimotor restriction (SMR). To restrict movement, the feet of the pup were first tied together with medical tape and then attached to a cast made of epoxy stick. The proximal part of the hind limb was also taped to the cast. The hip joint was free to move perpendicular to the knee, ankle and toe joints that remained in an extended position during the casting period of 16 hours per day from P1 to P28.</figcaption></figure>

the toes was generally located medial (toe 1 to toe 3, innervated by the saphenous nerve) and lateral (toe 3 to toe 5, innervated by the sciatic nerve) to the glabrous representation of the toes (innervated by the sciatic nerve).

Compared to control rats, the overall somatotopy of the foot maps was degraded in adult SMR rats that had experienced daily hind limb immobilization for a month postnatally. Specifically, the somatosensory representation of contiguous skin surfaces of the foot was dramatically disrupted in SMR rats relative to control rats (Fig. 2B,D). The total area of the somatosensory hind paw representation was 1.6 times smaller in SMR rats than in controls ($t = 5.35; df = 17; p < 0.0001$; Fig. 2F). Analyzed together, the cortical neural areas serving the glabrous and hairy foot surfaces were smaller in SMR rats than in controls (Wilks: $F(2,33) = 3.81; p < 0.04$), as were the normalized neural areas of both glabrous and hairy foot surfaces (Wilks: $F(2,33) = 3.76; p < 0.05$; Fig. 2F). Interestingly, the normalized glabrous area was reduced by SMR ($t = 2.36; df = 17; p < 0.03$), whereas the normalized hairy representation of the foot did not differ between the two groups ($t = 1.28; df = 17; p > 0.1$; n.s.; Fig. 2F), suggestive of a specific effect of developmental SMR on the somatosensory representation of glabrous foot skin surfaces.

Illustrative of changes in neuronal properties after postnatal SMR, the proportion of cortical sites with multiple receptive fields (RFs) (i.e., located on three or more disjunctive foot territories) was 10 times greater in SMR rats (30.8%; 99/322) than in control rats (3.1%; 22/702). In contrast, the percent of double RFs did not differ between SMR (24.5%; 79/322) and control (34.5%; 242/702) rats. Concomitantly, the proportion of single RFs located on a single hind paw location was decreased in SMR rats (44.7%; 144/322), compared to controls (62.4%; 138/702; $\chi^2 = 194.57; p < 0.0001$). In addition to the increased number of multiple RFs after SMR, their size also increased (Fig. 2C,E). The mean size of all RFs was greater in SMR rats than in controls ($t = -2.22; df = 17; p < 0.04$; Fig. 2G). When the very large dorsal foot RFs were discarded to reduce variation, the RFs were approximately twice as large in SMR rats as in controls ($t = -2.30; df = 17; p < 0.03$; Fig. 2C,E,G). Also, the size of glabrous and hairy RFs differed significantly between both groups of rats (Wilks: $F(2,33) = 4.06; p < 0.03$). Specifically, RFs located on glabrous foot skin surfaces were about twice as large in SMR rats as in controls ($t = -2.22; df = 17; p < 0.04$), whereas the size of hairy foot RFs did not differ significantly between the two groups (Fig. 2G), confirming a specific impact of SMR on glabrous foot surfaces. These enlargements lead to overlapping glabrous RFs and less specific foot representations in SMR rats (Fig. 2D,E).

Figure 2D illustrates the patchy representation of glabrous and hairy foot surfaces after SMR. There were more RFs encompassing both glabrous and hairy foot skin surfaces in SMR rats (9.3%; 30/322) than in controls (2.4%; 17/702; $\chi^2 = 23.97; p < 0.0001$), more multiple RFs located on three or more foot subdivisions, and numerous RFs located on both sides of the hind paw. In fact, the delineations between toe and plantar pad representations were difficult in SMR rats, so that map sites with three or more RFs had to be pooled together in specific map sectors (colored in orange in Fig. 2B,D) of the somatosensory foot maps. The map sectors with three or more RFs were much larger in SMR than in control rats (Cont: $0.02 \pm 0.02 \text{ mm}^2$; SMR: $0.19 \pm 0.14 \text{ mm}^2; t = -3.89; df = 17; p < 0.001$) and the normalized area of such cortical sectors relative to the total map area (expressed as percentages) was about 16 times larger after SMR (Cont: $1.40 \pm 1.58$%; SMR: $22.22 \pm 15.42$%; $t = -4.26; df = 17; p < 0.0005$). Thus, SMR induced a drastic degradation of both the topographical organization of the somatosensory hind limb maps and cortical neuronal properties. Interestingly, there was an emergence of RFs covering both sides of the leg in SMR rats (11.8%; 13/110) yet no such double-sided RFs in controls (0%; 0/110; $\chi^2 = 13.94; p < 0.0002$). However, the area of the somatosensory map devoted to the leg representation did not differ significantly between SMR rats ($0.49 \pm 0.28 \text{ mm}^2$) and controls ($0.67 \pm 0.30 \text{ mm}^2; t = 1.19; df = 13; p > 0.1$; n.s. Fig. 2B,D).

Next, neural responsiveness originating from the cortex to light tactile stimulation was qualitatively classified on the basis of the magnitude of the signal-to-noise ratio, as in a previous study²⁷ using a three-level scale (weak, clear and excellent responses). The proportion of excellent responses were greater in SMR rats (32.9%; 106/322) than in controls (12.8%; 90/702), the percent of clear responses decreased in SMR rats (47.8%; 154/322) compared to controls (70.1%; 492/702; $\chi^2 = 61.24; p < 0.0001$), while weak responses did not differ between SMR rats (19.3%; 62/322) and controls (17.1%; 120/702). SMR seemed to increase the neural responsiveness to tactile stimulation, indicative of changes in neuronal properties and increased excitability in the somatosensory hind limb representation.
A

Lateral

Medial

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

B

Medial

Lateral

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

C

<figure><img src="image_7.png" /></figure>

CONTROL

<figure><img src="image_8.png" /></figure>

D

Medial

Lateral

<figure><img src="image_9.png" /></figure>

1 mm

SENSORIMOTOR
RESTRICTION

<figure><img src="image_10.png" /></figure>

E

<figure><img src="image_11.png" /></figure>

10 mm

<figure><img src="image_12.png" /></figure>

F

<figure><img src="image_13.png" /></figure>

G

<figure><img src="image_14.png" /></figure>

**Figure 2.** Impact of early sensorimotor restriction (SMR) on the topographical organization of somatosensory map devoted to the skin surface representation of the hind limb in control and SMR rats. (A) Illustration of the color code for the somatosensory maps depicted in (B), whose cortical sites display receptive fields (RFs) located on the different skin subdivisions of ventral (left) and dorsal (right) regions of the hind limb, e.g., toes 1 to 5, plantar pads, heel and leg (see Supplementary Methods for details). The hatched zones correspond to cortical sites whose RFs are located on hairy skin surfaces of the hind limb. Sectors in orange correspond to cortical sites whose RFs are located on three or more skin subdivisions of the foot. (B) Representative somatosensory maps serving the cortical representation of hind limb skin surfaces in two control rats. All penetration sites of the microelectrode are depicted by squares, circles or triangles. Note the somatotopic organization despite idiosyncratic differences. (C) Illustration of all RFs located on glabrous hind paw surfaces corresponding to somatosensory maps from the same control rats shown in (B). (D) Somatosensory hind limb maps in two representative SMR rats. Note the reduction of the extent of the somatosensory maps after SMR, concomitant to an increase in the size of glabrous
RFs illustrated for the same SMR rats shown in (E). Changes in map areal extent and RF size are plotted in (F and G), respectively. (E) Illustration of glabrous RFs corresponding to somatosensory maps from the same SMR rats depicted in (D). (F) Plots of the areal extent of the different subparts of somatosensory hind limb maps in absolute areas (mm²; left Y axis) or normalized values relative to the ventral or dorsal foot skin surfaces (right Y axis). (G) Plots of the size of all RFs and their subcategories that correspond to those illustrated in (F). Mean ± SD. \*p < 0.05; \*\*p < 0.01; \*\*\*p < 0.001, significantly different from controls.

Thus, SMR favored the emergence of numerous multiple, disjunctive and enlarged RFs, and a decrease in the hind limb map area, resulting in degradation of the somatotopy and topographical organization of the somatosensory cortex, as well as altered neuronal properties. SMR induced-map changes appeared to be specific to glabrous foot skin areas, and the appearance of RFs covering both sides of the hind limb or foot is indicative of degraded neuronal properties after SMR.

**SMR-related alterations of hind limb motor maps.** The representation of hind limb movements around joints did not appear topographically organized but scattered throughout the hind limb motor map (Fig. 3A,B), as previously reported<sup>28,29</sup>. In control rats, the motor hip representation occupied on average 70.9% of the total area of the hind limb motor map (3.88 ± 0.50 mm²; Fig. 3D), while 6.2% was knee, 13.9% was ankle, 1.5% was toes and 7.5% was devoted to complex movements involving several joints of the leg simultaneously. After SMR, we found a reduction in the total area of hind limb motor maps (t = 2.29; df = 15; *p* < 0.04; Fig. 3C,D). Analyzed together, the motor areas devoted to movements of the hip, knee, ankle, toes and multi-joints did not differ between either group (Wilks: F(5,11) = 1.11; *p*: n.s.; Fig. 3C,D). The relative areas devoted to the different joint representations relative to the total motor map area (expressed as percentages) were not significantly altered (χ² = 0.10; *p*: n.s.; Fig. 3D).

The overall threshold (i.e., minimal amount of current to evoke a movement) of neuronal motor responses was drastically reduced after SMR (Wilcoxon = 2.60; *p* < 0.007; Fig. 3E). Analyzed together, thresholds needed to evoke movements of joints of the leg did not differ between both groups of rats (Wilks: F(5,12) = 1.18; *p*: n.s.; Fig. 3D). Interestingly, the mean threshold required to elicit whisker, neck and forepaw movements outside of the boundaries of the motor hind paw map was significantly reduced (Wilcoxon = 2.19; *p* < 0.03; Fig. 3E), so that SMR appeared to exert a nonspecific role on the motor threshold changes in general. Although the relative representation of each joint within the motor map did not change in SMR rats, the total area for hind limb motor maps was drastically reduced. In addition, the reduction of the minimal stimulation thresholds after SMR revealed an increase in motor excitability.

**SMR changes the balance between glutamate and GABA in the sensorimotor cortex.** To assess the balance between excitation and inhibition in the hind paw sensorimotor cortical region, we analyzed dialysates obtained during *in vivo* microdialysis for levels of extracellular glutamate and GABA using gradient HPLC coupled to laser detection. The extracellular concentration of glutamate was greater in SMR rats than in controls (t = 2.20; df = 12; *p* < 0.05); whereas, the extracellular concentration of GABA did not differ between the two groups (Fig. 4A).

We then performed semi-quantitative Western blots from brains of additional animals, using specific antibodies against vesicular glutamate transporter (vGLUT1) and vesicular GABA transporter (vGAT), reliable markers of excitatory and inhibitory transmission, respectively. The amount of vGLUT1 increased approximately 90% in SMR rats relative to controls (Wilcoxon = 44; *p* < 0.0004; Fig. 4B,C). In contrast, the amount of vGAT did not differ between the two groups (Wilcoxon = 29.5; *p* = 0.34, n.s.; Fig. 4B,C). Thus, SMR seemed to induce an increase in glutamate release in the sensorimotor cortex of adult rats while GABAergic levels appear unaltered.

**SMR did not induce brain histopathology.** No microscope evidence of histopathology was observed in the somatosensory, motor and cingular cortices or underlying white matter, including in the corpus callosum of any SMR rat compared to control rats. Also, there were no differences between the two groups with respect to APP, caspase 3 or GFAP immunoreactivity, as shown in Supplementary Fig. 1. Nor were there differences in Fluoro-Jade staining that would be indicative of axonal injury in the brains or spinal cords of either group (data not shown).

The numerical density of NeuN positive cells did not differ significantly between the two groups of rats in either the hind paw somatosensory representation (Cont: 99.3 ± 22.4 × 1000/mm³; SMR: 86.8 ± 9.5 × 1000/mm³; Wilcoxon = 12; *p* = 0.71, n.s.) or hind paw motor representation (Cont: 96.7 ± 19.3 × 1000/mm³; SMR: 92.0 ± 16.9 × 1000/mm³; Wilcoxon = 12; *p* = 0.71, n.s.). In the same line, the numerical density of GAD positive cells did differ not between the two groups of rats in either the hind limb somatosensory representation (Cont: 7.0 ± 2.4 × 1000/mm³; SMR: 7.8 ± 2.4 × 1000/mm³; Wilcoxon = 7.5; *p* = 0.62, n.s.) or the hind limb motor representation (Cont: 7.0 ± 1.8 × 1000/mm³; SMR: 6.3 ± 2.2 × 1000/mm³; Wilcoxon = 13; *p* = 0.54, n.s.). Thus, SMR appeared to have no significant impact on brain neuroanatomy, compared to controls.

**Functional relationships between sets of variables using PCA and correlations.** In the present study, 105 variables obtained from the same rats (Cont: n = 18; SMR: n = 16) were used to describe 7 sets of variables: somatosensory (n = 23 variables) and motor (n = 17) hind limb map features, neurotransmission (n = 4) and neurohistopathology (n = 6), as well as locomotion at P30 (n = 24) and P65 (n = 24) and musculoskeletal histopathology (n = 7) as previously reported<sup>23</sup>. In that study, locomotion on a treadmill was assessed at P30 and P65, as was musculoskeletal histopathology in these same rats (see data from all variables in Supplementary Table 1).
<figure><img src="image_15.png" /><figcaption>A</figcaption></figure>

☐ Hindlimb joint movement

■ No response

▼ Vibrissae/head/forelimb/trunk

<caption>B</caption>

Rostral

<figure><img src="image_16.png" /></figure>

Caudal

<figure><img src="image_17.png" /></figure>

CONTROL

<figure><img src="image_18.png" /><figcaption>C</figcaption></figure>

1 mm

<figure><img src="image_19.png" /></figure>

SENSORIMOTOR
RESTRICTION

<figure><img src="image_20.png" /><figcaption>D</figcaption></figure>

Motor map area

<figure><img src="image_21.png" /><figcaption>E</figcaption></figure>

<figure><img src="image_22.png" /></figure>

**Figure 3.** Effects of postnatal sensorimotor restriction (SMR) on the musculotopic organization of the motor representation of hind limb movements in control and SMR rats. (A) Illustration of movements of different hind limb joints elicited by intracortical stimulation. The color code for movements of each joint corresponds to motor map sectors devoted to the corresponding movement. (B) Typical motor hind limb maps in two representative control rats. Note the patchy representation of hind limb movements and the lack of musculotopy of M1 maps. (C) Motor maps of hind limb movements in two representative SMR rats. Note the reduction of the extent of the M1 map after SMR, concomitant to a decrease in the overall threshold needed to evoke movements, changes plotted in (D and E), respectively. (D) Plots of the areal extent of the different subparts of motor hind limb maps in absolute areas (mm²; left Y axis), or relative areas of joint movements normalized to the total area of the hind limb motor map (%, right Y axis). (E) Plots of the motor thresholds (μA) that correspond to the minimal current needed to evoke a leg joint movement. Mean ± SD. \*p < 0.05; \*\*p < 0.01, significantly different from controls.

First, to define the variables that contributed the most to differentiating the two groups of rats, we computed a PCA to compress the dataset into a single score for each rat, based on coordinates along a PC (Supplementary Table 1). The overall score along PC1 summarized 24.4% of the variance of the dataset. The mean score for
A *In vivo* microdialysis

<figure><img src="image_23.png" /></figure>

B Western Blotting

<figure><img src="image_24.png" /></figure>

C

<figure><img src="image_25.png" /></figure>

**Figure 4.** Effects of postnatal sensorimotor restriction (SMR) on the excitation/inhibition balance in the hind limb area of the sensorimotor cortex. (A) Plots of glutamate and GABA extracellular levels assessed using *in vivo* microdialysis. (B) Typical immunoblots for vGLUT1 (vesicular transporter of glutamate) and vGAT (vesicular transporter of GABA) in the two groups of rats. Tubulin controls show that loading was equal for each lane. (C) Plots of glutamate and GABA intracellular levels assessed by Western blots of vGLUT1 and vGAT, respectively. Note that both measures indicate an increase in glutamate levels while GABA remained unchanged after SMR, compared to controls. Mean ± SEM. \*p < 0.05; \*\*p < 0.001.

SMR rats (−4.1 ± 4.0) differed significantly from the score for controls (3.7 ± 1.8; Wilcoxon = 323; *p* < 0.0001), indicating a clear separation of the two groups along PC1 (Supplementary Fig. 2A). The highest contributing variables included: 1) locomotion kinematics and joint angles at P30 (15 variables out of 24, including 9 variables with factor loadings along PC1 ≥ |0.76|), 2) locomotion at P65 (16/24 variables including 3 variables with factor loadings ≥ |0.76|), 3) somatosensory hind limb map features (14/23 variables including 7 variables |0.53| < factor loadings < |0.73|), and 4) musculoskeletal histopathological changes (6/7 variables including 2 variables with factor loadings ≥ |0.62|). In contrast, variables for motor hind limb map features (2/17 variables including the smallest factor loadings = |0.45|), neurotransmission (1/4 variables) and neurohistopathology (0/6 variables) appeared to contribute little to changes observed in PC1 (see Supplementary Table 1).

Second, to evaluate the relationships between the 7 sets of variables, individual animal scores for each set of variables (Supplementary Fig. 2B) were compared. We observed that somatosensory hind limb map features were highly correlated to locomotion kinematics at P30 (r = 0.84; *p* < 0.0002) and P65 (r = 0.84; *p* < 0.0007), to musculoskeletal histopathology (r = 0.75; *p* < 0.0001) and tended to correlate with neurotransmission findings (r = 0.54; *p* = 0.08; Table 1). Surprisingly, the neurotransmission variables correlated with musculoskeletal histopathology (r = 0.51; *p* < 0.04) and somatosensory map features, but not with the other sets of variables. The motor hind limb map features along the PC2 axis significantly correlated to the highest contributing variable, which was locomotion at P30 (r = 0.53; *p* < 0.03), and to musculoskeletal histopathology (r = 0.50; *p* < 0.04), but not to the other sets of variables (Table 1). Musculoskeletal histopathology was highly correlated with locomotion at both P30 (r = 0.75; *p* < 0.0001), as well as at P65, although to a lesser level (r = 0.56; *p* = 0.003). In contrast, locomotion
<table><caption>**Table 1.** Correlations between the 7 sets of variables, based on the individual scores along PC1 (x-axis) or PC2 (y-axis) of animals obtained by using PCA. Pairwise two-sided p-values: $^{a}p < 0.05$; $^{b}p < 0.001$; $^{†}$tendency, $p = 0.08$.</caption><thead><tr><th></th><th>Somatosensory map features PC1</th><th>Motor map features PC2</th><th>Neurotransmission PC1</th></tr></thead><tbody><tr><td>Locomotion P30 PC1</td><td><strong>0.84<sup>b</sup></strong></td><td><strong>0.53<sup>a</sup></strong></td><td>0.30</td></tr><tr><td>Locomotion P65 PC1</td><td><strong>0.84<sup>b</sup></strong></td><td>0.06</td><td>0.27</td></tr><tr><td>Musculoskeletal pathologies PC1</td><td><strong>0.75<sup>b</sup></strong></td><td><strong>0.50<sup>a</sup></strong></td><td><strong>0.51<sup>a</sup></strong></td></tr><tr><td><strong>Neurotransmission PC1</strong></td><td><strong>0.54<sup>t</sup></strong></td><td>0.28</td><td></td></tr><tr><td><strong>Neurohistopatho. PC1</strong></td><td>0.36</td><td>0.43</td><td>0.19</td></tr></tbody></table>

features at both ages were highly interrelated ($r = 0.82$; $p < 0.0001$), suggestive of a self-perpetuating interplay between these variables, as described previously$^{23}$. Neurohistopathological variables did not correlate significantly to any other variables, confirming the lack of impact of SMR on brain neurohistology. Unfortunately, we could not assess the relationships between somatosensory and motor map features as we could not map both cortices in the same animals due to the use of different anesthetics.

Thus, the most contributing variables differentiating SMR rats from controls were locomotion at P30 and P65, then somatosensory map features, as well as the circumference and belly length of the gastrocnemius (indicative of muscle atrophy). SMR-related somatosensory map reorganization was highly related to locomotion impairments, excitation/inhibition imbalance and musculoskeletal histopathology. Surprisingly, motor map features and neurotransmission contributed the least to group differentiation. Reorganization of the motor hind limb representation was also related to locomotion impairments, although only at P30, and to histopathology, although to a lesser extent than the degradation of the somatosensory hind limb representation. Peripheral histopathology and locomotion impairments were highly correlated, confirming the self-perpetuating interplay between these variables.

We also found that sex had no specific impact on the features explored here and previously$^{23}$.

## Discussion

To investigate a hypothesis of a developmental origin of adult disorders, we showed in the present study that reduced and abnormal patterns of somatosensory inputs after hind limb movement restriction during postnatal development had no significant impact on adult brain neuroanatomy yet did effect the functional organization of the sensorimotor cortex. Compared to controls, SMR led to: 1) a degradation of both somatosensory hind limb map organization and neuronal properties; 2) a decrease in the total motor map area devoted to hind limb movements concomitant with an overall increased excitability; 3) increased glutamate release, yet no changes in GABA levels in the sensorimotor cortex; yet 4) no neuroanatomical changes in the sensorimotor hind limb cortex.

To evaluate interactions underlying the observed neural maladaptive plasticity, we used PCA and correlational statistics. We found that SMR primarily effected locomotion at P30 and P65, as well as somatosensory map features and musculoskeletal features in decreasing order; yet had lesser effect on motor map features, neurotransmission or neurohistology. This meta-analysis of all data from the same rats also outlines that: 1) deterioration of somatosensory maps is tightly related to locomotion impairments, neurotransmission imbalance and musculoskeletal histopathology; 2) locomotor disorders and musculoskeletal pathology likely interplay tightly in a self-perpetuating cycle (see$^{23}$); whereas, 3) motor map alteration was weakly influenced by these variables and not by neurotransmission.

**Putative mechanisms of plasticity underlying the functional disorganization of the hind limb representations.** SMR induced a specific and drastic reduction of the somatosensory representation of the glabrous foot surfaces concomitant with a specific and marked enlargement of glabrous RFs, suggestive of a use-dependent (Hebbian) process of plasticity$^{30,31}$. The emergence of cortical sites with three or more disjunctive RFs and/or RFs located on both glabrous and hairy sides of the foot or leg in adult rats led to a greater overlap between glabrous RFs, a coarser-grained representation of the somatosensory foot representation. Most importantly, this emergence led to a dedifferentiation of the cortical representation of both single skin territories, such as toes and plantar pads, and glabrous/hairy surfaces that are usually represented segregated in control rats$^{27}$. In sum, SMR degraded the overall somatotopy and fine-detailed topography of the somatosensory hind limb representation and neuronal properties that correspond to degraded abilities in tactile discrimination$^{32}$. SMR induced a decrease in the motor leg map area while the relative cortical areas serving the different joint representations were unchanged, so that the overall “musculotopy” of hind limb movements was maintained despite inter-individual variations. Another key result after SMR was the decrease in the overall threshold of neurons in the motor cortex needed to evoke leg movements, a decrease that was not specific to the leg area. This threshold decrease corresponded to an overall increase in motor cortex excitability that fits with both increased cortical neuron responsiveness to tactile stimulation after SMR and increased glutamatergic transmission, yet unchanged GABA levels in the sensorimotor leg area. These results confirm and extend prior studies using similar hind limb immobilization during development$^{27,29,33}$, in which reorganization originating the sensorimotor cortex was extensively discussed, although the literature about chronic activity deprivation during development is limited.

We investigated the excitation-inhibition neurotransmission in the sensorimotor cortex because it appears to govern cortical map reorganization. Use-dependent regulation of glutamatergic excitation underlies expansion or retraction of somatosensory maps, while use-dependent GABAergic inhibition promotes cutaneous RF changes, corresponding to the modulation of neuronal tuning on the basis of synaptic potentiation/depression and
unmasking/masking of thalamocortical divergent and convergent connections, respectively³⁰,³¹,³⁴. In the same line, activity-dependent changes in excitation-inhibition balance leads to masking or unmasking of intracortical motor connections that acts upon motor map boundaries and thresholds, i.e., neuronal tuning³⁵,³⁶. Decreased hind limb sensorimotor maps and increased size of RFs did not match with the overall hyperexcitability found in the hind limb sensorimotor area, so that the hind limb map reorganization may be explained by alternative mechanisms of plasticity. The absence of neuroanatomical pathology after SMR allows us to rule out its impact on the neurotransmission changes and map reorganization. In addition, SMR induced a long-lasting increase in the stretch reflex evaluated with post-activation depression²³, that is supposed to reflect hyperexcitability within the lumbar spinal cord and spasticity³⁷,³⁸. The latter finding fits with the present observation of enduring cortical hyperexcitability.

In this long-term study, a slow homeostatic (non-Hebbian) increase in global synaptic strength (i.e., synaptic scaling) and/or intrinsic excitability, which in turn enhances responses to remaining inputs³¹,³⁹, may explain the increased cortical excitation observed in adult SMR rats. An alternative explanation for such an increase could reside in the persistence of dense dendritic arborisation that are usually pruned during the first few postnatal weeks⁴⁰. For instance, postnatal whisker removal reduces the pruning of supernumerary dendrites that are normally present in S1 layer V pyramidal cells⁴¹, thus preserving abnormally important levels of excitatory synapses involved in increased cortical excitability⁴². Thus, a co-existence of both Hebbian and homeostatic processes of plasticity may explain the observed differences in remodelling of somatosensory and motor maps in relation to hyperexcitability in the sensorimotor cortex.

**Putative interactions and processes underlying maladaptive plasticity.** Based on PCA and correlations within the same animals, our meta-analysis indicates that the SMR-related deterioration of the somatosensory map organization and neuronal properties was tightly related to locomotor disorders evaluated at P30 and P65. SMR led to severely degraded kinematics of gait on a treadmill in relation to knee and ankle overextension, and smaller footprints that revealed a digitigrade locomotion²³, which resembles "toe walking" or true pes equinus, observed in children with ASD²⁴. According to the literature³⁰,³¹,⁴³,⁴⁴, typical plantigrade locomotion in control rats²³,⁴⁵ produces asynchronous and distributed stimulation through the foot skin surfaces from proximal to distal parts during the stance up to the beginning of the swing. During typical motor development, spatiotemporally distributed cutaneous inputs of the foot in controls likely contribute to the differentiated representation of the different foot skin surfaces, such as individual toes and plantar pads. In contrast, we postulate that digitigrade locomotion during postnatal SMR provides synchronous and colocalized tactile inputs from several toe and/or pad surfaces, leading to enlarged RFs that cover abnormally large skin surfaces. This enlargement results in the dedifferentiation of the somatosensory representation of contiguous foot territories and finally to the disruption of the somatotopic and topographic organization of the somatosensory foot maps, also observed in animal models of focal hand dystonia⁴⁴,⁴⁶. In addition, unrepresented skin surfaces in control rats, such as the proximal phalanges of toes and skin surrounding pads, became represented in the somatosensory maps after SMR. The very large RFs covering the entire leg of SMR rats may result from synchronous activations of cutaneous mechanoreceptors during chronic co-contractions of antagonist leg muscles (see²⁷). Taken together, abnormal, specific tactile stimulation during digitigrade locomotion induced a maladaptive reorganization in the somatosensory cortex, in which only the glabrous representation was altered by SMR, not the hairy representation (see²⁷). Neural maladaptive plasticity of the somatosensory maps also appeared highly related to musculoskeletal histopathology and neurotransmission.

Those correlations do not indicate a directionality of the relation, so that it means that locomotor disorders and musculoskeletal histopathology may be attributed to the somatosensory map degradation. Yet, the high correlation between locomotion degradation and musculoskeletal histopathology after SMR confirms their interplay in a self-perpetuating cycle (Fig. 5), in which the spasticity related to increased stretch reflex also contributed to this cycle²³. Surprisingly, only locomotion disorders at P30 and musculoskeletal histopathology correlated weakly with the motor map reorganization. In addition, neurotransmission appeared correlated primarily with somatosensory map features and secondarily with musculoskeletal pathologies, but not with motor map features or locomotion. Taken together, our results emphasize the key role of the functional degradation of the somatosensory hind limb representation, which appears to mainly depend on abnormal patterns (quality) of somatosensory inputs resulting from both postnatal SMR and enduring degraded reafference resulting from disturbed locomotion rather than decreased levels (quantity) of inputs, which likely seems involved in the excitation/inhibition imbalance⁴⁷,⁴⁸. These results proceed from dynamic selection, competition and collaboration of primary neuronal groups on the basis of afferent inputs resulting from behaviors and experience that induce changes in the strength of synaptic connections, thus governing the topographic organization of cortical maps and the corresponding tuning of neurons, and leading to diverse repertoires of movements⁴⁹⁻⁵¹. Thus, the development and continuous refinement of neural connectivity is influenced interactively by feedback from neural targets⁵², supporting the interplay between neural connections, maps and functions intermingled in regulatory loops and self-perpetuating cycles sustained by experience.

Therefore, we postulate that early limited and abnormal patterns of movements through postnatal hind limb immobilization during the light active period provide atypical sensorimotor experience at the origin of adult locomotor disorders (e.g. digitigrade gait), increased stretch reflex (i.e., lumbar spinal cord hyperexcitability, spasticity, spasms, contractures...) and musculoskeletal tissue degradation, compared to typical sensorimotor development in controls (Fig. 5²³). Our meta-analysis suggests the existence of a self-perpetuating cycle, in which musculoskeletal pathologies likely contribute to muscle overactivity and locomotor disorders. Spasticity and related muscle disorders may contribute in turn to enhanced musculoskeletal and locomotion disturbances (Fig. 5²³). As a consequence, degraded movements due to spastic muscles and digitigrade locomotion likely produce atypical somatosensory feedback to the immature sensorimotor circuitry, including the lumbar spinal
<figure><img src="image_26.png" /></figure>

<figure><img src="image_27.png" /></figure>

**Figure 5.** Illustration of the interactions and processes underlying typical or atypical sensorimotor development in rats. (A) Schematic representation of the putative interactions between the different variables assessed (locomotion, musculoskeletal tissues, Hoffman-reflex, somatosensory and motor representations of the hind limb and excitation-inhibition neurotransmission in the sensorimotor cortex) that may occur during typical sensorimotor development in control rats. Daily activities provide early typical sensorimotor experience that allow muscle synergies and typical gait to interplay harmoniously. Normal movements during typical gait produce typical sensorimotor feedback to the immature brain, contributing to the development and refinement of connections, maps and functions, e.g., normal motor command to muscles that in turn contribute to muscle synergies and typical movements in a virtuous self-perpetuating cycle. Excitation-inhibition balance also participates in the refinement of neural maps and connections. (B) Schematic illustration of the deleterious impact of postnatal sensorimotor restriction (SMR) on neuromuscular interplay. Movement restriction during postnatal hind limb immobilization provides early atypical sensorimotor experience that leads to gait disorders, musculoskeletal pathologies and hyperreflexia, as signs of spasticity (see²³). In turn, gait disorders and spastic muscles provide atypical somatosensory feedback/reafference to the immature sensorimotor circuitry, mainly the spinal network and the primary somatosensory (S1) and motor (M1) cortices that intercommunicate. Such atypical somatosensory feedback likely degrades the functional organization of the somatosensory and motor maps and spinal network, thus producing abnormal motor commands that finally aggravate gait disorders, spasticity and musculoskeletal abnormalities into interrelated self-perpetuating cycles. Excitation-inhibition imbalances that tip the balance to hyperexcitability in sensorimotor hind limb cortices have more impact on the reorganization of somatosensory maps than on motor maps. Hyperexcitability in the lumbar spinal cord has been shown previously as an enduring hyperreflexia in adult SMR rats (see²³).

circuitry involved in automatic gait, and primary somatosensory and motor cortices associated with other sub-cortical structures in the sensorimotor loop¹⁹. In the absence of brain damage, this atypical somatosensory rea-ference probably triggers and maintains the functional disorganization of the sensorimotor maps and neuronal properties, and likely that of the lumbar spinal cord (i.e., hyperexcitability related to increased stretch reflex), involved in both motor commands during automatic locomotion and the emergence of muscle overactivity and paresis (Fig. 5³,⁵⁴). The disorganized sensorimotor circuitry may in turn produce atypical or aberrant motor com-mands, such as abnormal muscle synergies involved in gait disorders observed in rats after hind limb unloading⁵⁵ or in patients with stroke⁵⁶,⁵⁷. Thus, we propose the existence of three interplaying, self-perpetuating cycles: the first cycle is related to peripheral changes, including gait disorders, increased stretch reflex and musculoskeletal pathologies (Fig. 5). The second cycle involves the sensorimotor circuitry. The third one corresponds to inter-actions between the peripheral and central circuitry that are interconnected by the somatosensory reafference (altered by abnormal movements) and atypical motor commands that drive abnormal movements (Fig. 5).
**Functional implications.** The present study based on postnatal SMR confirms and emphasizes the pre-
ponderant role of individual, early experience in shaping the body and brain, and modulating their interplay
from maturation to adulthood, depending on physical and environmental constraints⁴,¹⁹,⁵². Limited and abnormal
patterns of somatosensory inputs arising from SMR in early development, in the absence of neuroanatomical
histopathology, induced long-lasting locomotor disorders, hind limb musculoskeletal histopathology, hyperre-
flexia and functional disorganization of the adult sensorimotor cortex. Taken together, this recapitulates several
symptoms and neurological features found in patients with ASD, and DCD in general¹⁰,¹⁸,²². Adults with DCD
that have no evident brain damage, show reduced abilities to produce consistent movements¹¹ and impairments
in kinesthetic acuity¹³. Children with ASD also exhibit gross and fine motor abnormalities, motor learning defi-
ciencies, inabilities in executing sequences of actions¹⁶ and deficits in tactile localization that may reflect broad
impairments in multisensory body representations¹³⁻¹⁵. Thus, these children share early aberrant sensorimo-
tor experience resulting from atypical spontaneous or general movements (GMs) observed in children with
DCD²,¹⁷,⁵⁸, along with its related brain disorganization, especially inabilities in processing and integrating sensory
information²,¹³,¹⁶. GMs are predictive of neural, motor and cognitive outcomes in children with DCD²,¹⁰,²¹,²²,⁵⁸.
Our animal studies on postnatal SMR supports the idea that developmental maladaptive plasticity may lead to
long-lasting motor disturbances, as observed in DCD and point out the pertinence of early rehabilitative pro-
grams. Further studies are needed for a better understanding of the possible impact of early SMR on behavioral
and cognitive abilities and on the postnatal development of the lumbar spinal circuitry.

**Materials and Methods**

All experiments and animal use were carried out in accordance with the guidelines laid down by NIH (NIH Publication #80-23) and EC Council Directive (2007/526/EEC). This research involving animals was approved by the Direction Départementale des Services Vétérinaires – Préfecture des Bouches du Rhône, France (permit # C13-055-18).

**Animals and hind limb movement restriction.** Sprague-Dawley rat pups of either sex from different litters were pseudo-randomly assigned to two groups: 1) a group subjected to transient hind limb immobilization from P1 to P28 for 16 hours per day, thus producing sensorimotor restriction (SMR); and 2) a control group (Cont). Feet of SMR pups were gently bound together with medical tape. Their hind limbs were then immobilized in an extended position and taped to a cast made of epoxy putty sticks (Fig. 1). The casts were well tolerated by the pups and mothers, and allowed the pups to move at the hip, urinate, defecate and receive maternal care. After casting, pups were returned to their mother and unrestrained littermates during the dark phase, a period of peak motor activity. During the day's light phase, the casts were removed so that pups could move freely for 8 h/day. Daily, while uncasted, the hind limb joints were passively moved by the investigators through their full range of motion. The size of the casts was adapted to the growth of the rats from P1 to P28 (see²³ for details).

**Electrophysiological mapping in cortical sensorimotor areas.** The topographical organization and neuronal properties of the primary somatosensory (S1) and motor (M1) cortices were examined in rats from 90 to 120 days-old, using either microelectrode electrical recordings or stimulation, respectively. Hind limb maps were derived from the left hemisphere, as these maps correspond to the right contralateral hind limb whose locomotion kinematics were previously recorded and reported²³.

*Somatosensory hind limb maps.* The somatosensory mapping procedure used was as previously described<sup>27,46</sup>
(for more details, see Supplementary Methods). Briefly, each animal was anaesthetized with sodium pentobarb-
bital (50 mg.Kg<sup>-1</sup>, i.p.); supplemental doses (5 mg.Kg<sup>-1</sup>, i.p.) were provided as needed to maintain a deep and
stable level of anaesthesia throughout the mapping session. A craniotomy was made over the S1 hind limb rep-
resentation located between −1 and +3 mm from bregma in the rostrocaudal axis and between 1 and 5 mm in
the mediolateral direction. Multiunit recordings were made using tungsten microelectrodes (1 MΩ at 1 KHz,
WPI, Sarasota, FL, USA) in the upper layer IV (650–700 µm) of the S1 cortex. The low-threshold receptive field
(RF) corresponded to the cutaneous area whose barely visible indentation or gentle hair movements elicited clear
bursts of S1 neuronal activity. High-threshold responses elicited by taps on the skin were classified as cutaneous
responses; stroking of hairs and manipulations of muscles and joints were classified as non-cutaneous responses.
Unresponsive cortical sites exhibited spontaneous discharges only. To elaborate maps of the hind limb representa-
tion we drew boundaries encompassing cortical sites whose RFs were restricted to a common hind paw subdivi-
sion (e.g. toe, plantar pad). When RFs were located on distinct and separate skin subdivisions of the foot, borders
were drawn midway between adjacent recording sites. When a single or multiple RF included different but adjoining
skin subdivisions of the hind limb, a boundary line was drawn that crossed the cortical sites. Map borders
were placed midway between responsive and unresponsive sites (Fig. 2B,D). When RFs were located on three
or more disjunctive RFs, we pooled together cortical sites in specific map sectors colored in orange (Fig. 2B,D).

*Motor hind limb maps.* Standard intracortical microstimulation techniques were used to derive detailed maps of hind limb movements<sup>28,29</sup> (see Supplementary Methods). Animals were anaesthetized with ketamine hydrochloride (70 mg.Kg<sup>-1</sup>, i.p.) and xylazine (5 mg.Kg<sup>-1</sup>, i.p.). Supplemental doses of ketamine (20 mg.Kg<sup>-1</sup>, i.p.) and acepromazine (0.02 mg.Kg<sup>-1</sup>, i.p.) were delivered as needed to maintain a deep and constant level of anaesthesia. A craniotomy was made over the M1 hind limb representation located between −1 and +3 mm from bregma in the rostrocaudal axis and between 1 and 5 mm in the mediolateral direction, indicative of large overlap of somatosensory and motor hind limb maps (See Supplementary Fig. 3). Microelectrodes (100 KΩ, FHC, Bowdoin, ME, USA) were advanced into layer V at a depth of about 1700 µm. At each electrode penetration, stimulation was initiated at the lowest intensity (5 µA) and was increased gradually until lower limb movement was evoked.
Movements of the toes, ankle, knee, hip and other joints were detected visually. If no movement was elicited at 60 µA, the site was defined as unresponsive. Maps were reconstructed offline. Detailed mapping procedures are provided in Supplementary Methods.

Somatosensory map areas were described by their absolute areas (mm²) for each rat and normalized relative to either the ventral or dorsal foot skin surfaces for somatosensory maps only. Relative areas of joint movements are normalized relative to the total motor hind limb map area for each rat. Average values were computed for both groups of rats. The absolute sizes of glabrous and hairy receptive fields (RFs), measured in mm², were normalized relative to the ventral and dorsal hind paw skin surfaces, respectively, and expressed as percentages. The relative RF areas measured in each rat were averaged, and mean RF sizes in each group were computed. Mean RF size per rat was computed after discarding the large RFs located on the dorsal surface of the hind paw²⁷.

**Excitatory and inhibitory neurotransmission.** To gain insight into the long-term impact of SMR on extra- and intracellular levels of glutamate and GABA, in subsets of rats, we performed *in vivo* microdialysis within the sensorimotor cortex located contralateral to the mapped side, at 2 hours before mapping. In other rats, Western blot analysis was used to quantify intracellular amounts of transporters for both glutamate (vGLUT1) and GABA (vGAT) in sensorimotor cortical tissues collected contralateral to the mapped side, at one hour after mapping.

*In vivo microdialysis.* Carnegie Medecin microdialysis probes (CMA/11, Phymep, France) were implanted within the hindpaw representation of the right S1-M1 area using stereotaxic coordinates (A: -1/0, L: + 2; H: 2). An implantation site that avoided large blood vessels; traces of blood were not found after probe withdrawal. The terminal ends of the probes were covered with a polycarbonate membrane. The membrane had a diameter of 0.24 mm, a length of 2 mm and a 6 KDa molecular mass cut-off. The membrane acted like a blood vessel, with extracellular molecules passing through it by diffusion gradient. A new probe was used for each rat. This procedure is detailed in Supplementary Methods.

*Western blotting and quantification.* The leg region from the right side of the sensorimotor cortex was collected and homogenized in TBS (Tris-HCl 50 mmol.L⁻¹ pH 7.4, NaCl 150 mmol.L⁻¹) containing protease inhibitors (complete EDTA free, Roche, Basel, Switzerland) using a Potter homogenizer and a ratio of 100 µg of brain pieces per mL of buffer. Homogenates were centrifuged at 500 g for 3 min at 4 °C. Supernatants were aliquoted and stored at -80 °C until use. Proteins present in each supernatant were resolved by SDS-PAGE and transferred onto nitrocellulose membranes. Proteins of interest were revealed with specific antibodies and then enhanced using chemiluminescence (ECL; Pierce Biotechnology, Rockford, IL, USA). For illustration purposes, image editing was performed using ImageJ software (http://rsb.info.nih.gov/ij/) and was limited to linear brightness/contrast adjustment (see original scans of nitrocellulose membranes in Supplementary Fig. 4). Quantification was performed using ImageJ. Integrative intensities minus background were plotted for each sample after normalization to 1 for the highest value per western blot (see Supplementary Methods).

**Brain immunohistology.** After mapping one side of the brain, and microdialysis or collection of tissues for western blot analysis of the contralateral side, the animals received terminal doses of anaesthesia (thiopentobarbital, 150 mg.Kg⁻¹), before undergoing transcardial perfusion with 4% paraformaldehyde in phosphate buffer (0.10 M, pH 7.4). The brains were then harvested, postfixed in the same fixative solution for 3 h, cryoprotected in sucrose 30%, and frozen for cryosectioning. Sections to be immunostained for glial fibrillary acidic protein (GFAP), amyloid-β precursor protein (APP), DAPI (a general cell nucleus marker), NeuN (a neuronal nucleus marker) and GAD (GABAergic interneurons expressing GAD67) consisted of 15 µm coronal sections (see²⁸; Supplementary Methods).

Counting of NeuN and GAD-67 cells was performed on the same side of the brain as used for electrophysiology, using a DMR microscope (Zeiss, Oberkochen, Germany) equipped with a motorized stage and Stereo Investigator software-controlled computer system (Microbright-Field Europe, Magdeburg, Germany). Cell counting was performed at x40 magnification using an optical fractionator method in semi serial sections (see Supplementary Methods). In both groups, the number of samples from the left and right hemispheres was roughly equal.

The heights of the genu of the corpus callosum, and S1, M1 and entorhinal cortices were measured in haematoxylin and eosin (H&E) stained sections using a 20x objective and the Auto-Width tool of Bioquant Image Analysis program in which the inner and outer boundaries of structures were traced, and then mean layer thicknesses of a length of genu or cortex were automatically generated. Alternate sections were stained with Fluoro-Jade using previously described methods²⁸.

**Data analysis.** Data normality and homogeneity of variance were determined with Shapiro test, Bartlett test and var.test before to apply parametric (Welch two sample t test according to interval confidence for independent variables and Wilks one-way MANOVA for dependent variables) or nonparametric (Wilcoxon rank test, Wilcoxon; and chi-square, χ²) tests with normal approximation with continuity correction or correlations (Pearson's product-moment with normal approximation with continuity correction - Yates correction to the chi-square test – and confidence interval or nonparametric Spearman rank-order with two-sided alternative hypothesis) using only R (The R Foundation for Statistical Computing, Institute for Statistics and Mathematics, Wien, Austria). The investigators were blinded to rearing conditions throughout the different experimental sessions until statistical comparisons were performed.
**Principal components analysis.** The high dimensionality of our dataset (105 variables) made it difficult to easily extract parameters that accounted for SMR-related changes. To identify variables that contributed the most to the overall variability of our data, we used principal components analysis (PCA) to reduce high dimensionality (FactomineR and Psych packages, The R Foundation for Statistical Computing). PCA is defined as an orthogonal linear transformation of a dataset into a reduced subspace with new coordinates, so that the variance is maximized on each new coordinate axis or principal component (PC), thus minimizing the loss of information. A normalization method, in which the data mean was adjusted to 0 and the standard deviation (SD) was adjusted to 1, was applied to allow comparisons of data with disparate values and variances⁵⁹. Then, in order to identify variables that most optimally separate the two groups, we performed a PCA with all dataset variables (n = 105). We usually kept the first two PCs with Kaiser's criterion over 1 that in combination explained at least 35 to 70% of the variance⁵⁹. In the PC analysis of this dataset, we kept the factor loadings over |0.35| of variables with significant contributions (p < 0.05) to reduce dataset dimensionality and to optimally separate animals of each group (Supplementary Table 1 and Supplementary Fig. 2A). Next, in order to assess possible relationships between sets of variables within the same animals, PCA was performed on each set (n = 7) of variables, i.e., somatosensory and motor foot map features (n = 23 and 17 variables, respectively), neurotransmission (n = 4) and neurohistopathology (n = 6), as well as 3 more variables: locomotion at P30 (n = 24) and P65 (n = 24) and musculoskeletal histopathologies (n = 7) previously reported²³. We thus obtained a score for each rat and for each set of variables along the PC that optimally differentiated the two groups of rats, depending on the set of variables (Supplementary Fig. 2B). The aim was to correlate the scores of each rat between the 7 sets of variables, i.e., the scores that optimally characterized each animal for a given set. We used these optimal scores to perform subsequent linear correlations (Pearson or Spearman correlations according to data normality) within the same animals by using R.

## References

1. Geertsen, S. S., Willerslev-Olsen, M., Lorentzen, J. & Nielsen, J. B. Development and aging of human spinal cord circuitries. *J. Neurophysiol.* **118**, 1133–1140 (2017).

2. Hadders-Algra, M. Neural substrate and clinical significance of general movements: an update. *Dev. Med. Child Neurol.* **60**, 39–46 (2018).

3. Hadders-Algra, M. Variation and variability: key words in human motor development. *Phys. Ther.* **90**, 1823–1837 (2010).

4. Brumley, M. R., Kauer, S. D. & Swann, H. E. Developmental plasticity of coordinated action patterns in the perinatal rat. *Dev. Psychobiol.* **57**, 409–420 (2015).

5. Tiriac, A. & Blumberg, M. S. Gating of reafference in the external cuneate nucleus during self-generated movements in wake but not sleep. *Elife* **5**, e18749 (2016).

6. Taub, E. Parallels between use of constraint-induced movement therapy to treat neurological motor disorders and amblyopia training. *Dev. Psychobiol.* **54**, 274–292 (2012).

7. Frenken, R. Psychology and history of swaddling. Part One: Antiquity until 15th century. *J. Psychohist.* **39**, 84–114 (2011).

8. Sleuwen, B. E. van et al. Swaddling: A Systematic Review. *Pediatrics* **120**, e1097–e1106 (2007).

9. Kitase, Y. et al. A new type of swaddling clothing improved development of preterm infants in neonatal intensive care units. *Early Hum. Dev.* **112**, 25–28 (2017).

10. Zwicker, J. G., Missiuna, C., Harris, S. R. & Boyd, L. A. Developmental coordination disorder: A review and update. *Eur. J. Paediatr. Neurol.* **16**, 573–581 (2012).

11. Du, W., Wilmut, K. & Barnett, A. L. Level walking in adults with and without Developmental Coordination Disorder: An analysis of movement variability. *Hum. Mov. Sci.* **43**, 9–14 (2015).

12. Coleman, R., Piek, J. P. & Livesey, D. J. A longitudinal study of motor ability and kinaesthetic acuity in young children at risk of developmental coordination disorder. *Hum. Mov. Sci.* **20**, 95–110 (2001).

13. Johnston, J. S., Begum Ali, J., Hill, E. L. & Bremner, A. J. Tactile localization performance in children with developmental coordination disorder (DCD) corresponds to their motor skill and not their cognitive ability. *Hum. Mov. Sci.* **53**, 72–83 (2017).

14. Gabbard, C. & Bobbio, T. The Inability To Mentally Represent Action May Be Associated With Performance Deficits in Children With Developmental Coordination Disorder. *Int. J. Neurosci.* **121**, 113–120 (2011).

15. Kashuk, S. R., Williams, J., Thorpe, G., Wilson, P. H. & Egan, G. F. Diminished motor imagery capability in adults with motor impairment: An fMRI mental rotation study. *Behav. Brain Res.* **334**, 86–96 (2017).

16. Chinello, A., Di Gangi, V. & Valenza, E. Persistent primary reflexes affect motor acts: Potential implications for autism spectrum disorder. *Res. Dev. Disabil*, https://doi.org/10.1016/j.ridd.2016.07.010 (2016).

17. Einspieler, C., Peharz, R. & Marschik, P. B. Fidgety movements – tiny in appearance, but huge in impact. *J. Pediatr. (Rio J.)* **92**, S64–S70 (2016).

18. Vaivre-Douret, L., Lalanne, C. & Golse, B. Developmental Coordination Disorder, An Umbrella Term for Motor Impairments in Children: Nature and Co-Morbid Disorders. *Front. Psychol.* **7**, e502 (2016).

19. Blumberg, M. S. & Dooley, J. C. Phantom Limbs, Neuroprosthetics, and the Developmental Origins of Embodiment. *Trends Neurosci.* **40**, 603–612 (2017).

20. Khazipov, R. & Buzsáki, G. Early patterns of electrical activity in the developing cortex. in *Oxford handbook of developmental behavioral neuroscience* 161–177 (Blumberg MS; Freeman JH; Robinson SR 2010).

21. Zuk, L. Fetal and Infant Spontaneous General Movements as Predictors of Developmental Disabilities. *Dev. Disabil. Res. Rev.* **17**, 93–101 (2011).

22. Einspieler, C. et al. Highlighting the first 5 months of life: General movements in infants later diagnosed with autism spectrum disorder or Rett Syndrome. *Res. Autism Spectr. Disord.* **8**, 286–291 (2014).

23. Delcour, M. et al. Early movement restriction leads to enduring disorders in muscle and locomotion. *Brain Pathol*, https://doi.org/10.1111/bpa.12594 (2018).

24. Vilensky, J. A., Damasio, A. R. & Maurer, R. G. Gait disturbances in patients with autistic behavior: a preliminary study. *Arch. Neurol.* **38**, 646–649 (1981).

25. Mostofsky, S. H. et al. Decreased connectivity and cerebellar activity in autism during motor task performance. *Brain J. Neurol.* **132**, 2413–2425 (2009).

26. Nielsen, J. B., Perez, M. A., Oudega, M., Enriquez-Denton, M. & Aimonetti, J.-M. Evaluation of transcranial magnetic stimulation for investigating transmission in descending motor tracts in the rat: TMS and transmission in rat descending motor tracts. *Eur. J. Neurosci.* **25**, 805–814 (2007).

27. Coj, J. O. et al. Impact of neonatal asphyxia and hind limb immobilization on musculoskeletal tissues and S1 map organization: implications for cerebral palsy. *Exp. Neurol.* **210**, 95–108 (2008).

28. Delcour, M. et al. Neuroanatomical, Sensorimotor and Cognitive Deficits in Adult Rats with White Matter Injury following Prenatal Ischemia. *Brain Pathol.* **22**, 1–16 (2012).
29. Strata, F., Coq, J. O., Byl, N. & Merzenich, M. M. Effects of sensorimotor restriction and anoxia on gait and motor cortex organization: implications for a rodent model of cerebral palsy. *Neuroscience* **129**, 141–156 (2004).

30. Xerri, C. Imprinting of idiosyncratic experience in cortical sensory maps: neural substrates of representational remodeling and correlative perceptual changes. *Behav. Brain Res.* **192**, 26–41 (2008).

31. Feldman, D. E. Synaptic mechanisms for plasticity in neocortex. *Annu. Rev. Neurosci.* **32**, 33–55 (2009).

32. Xerri, C., Bourgeon, S. & Coq, J. O. Perceptual context-dependent remodeling of the forepaw map in the SI cortex of rats trained on tactile discrimination. *Behav. Brain Res.* **162**, 207–221 (2005).

33. Marcuzzo, S. *et al*. Different effects of anoxia and hind-limb immobilization on sensorimotor development and cell numbers in the somatosensory cortex in rats. *Brain Dev.* **32**, 323–331 (2010).

34. Rosselet, C., Zennou-Azogui, Y. & Xerri, C. Nursing-induced somatosensory cortex plasticity: temporally decoupled changes in neuronal receptive field properties are accompanied by modifications in activity-dependent protein expression. *J. Neurosci. Off. J. Soc. Neurosci.* **26**, 10667–10676 (2006).

35. Sanes, J. N. & Donoghue, J. P. Plasticity and primary motor cortex. *Annu. Rev. Neurosci.* **23**, 393–415 (2000).

36. Young, N. A., Vuong, J. & Teskey, G. C. Development of motor maps in rats and their modulation by experience. *J. Neurophysiol.* **108**, 1309–1317 (2012).

37. Bos, R. *et al*. Activation of 5-HT2A receptors upregulates the function of the neuronal K-Cl cotransporter KCC2. *Proc. Natl. Acad. Sci. USA* **110**, 348–353 (2013).

38. Boulenguez, P. *et al*. Down-regulation of the potassium-chloride cotransporter KCC2 contributes to spasticity after spinal cord injury. *Nat. Med.* **16**, 302–307 (2010).

39. Maffei, A. & Turrigiano, G. The age of plasticity: developmental regulation of synaptic plasticity in neocortical microcircuits. *Prog. Brain Res.* **169**, 211–223 (2008).

40. Cline, H. & Haas, K. The regulation of dendritic arbor development and plasticity by glutamatergic synaptic input: a review of the synaptotrophic hypothesis. *J. Physiol.* **586**, 1509–1517 (2008).

41. Zuo, Y., Yang, G., Kwon, E. & Gan, W.-B. Long-term sensory deprivation prevents dendritic spine loss in primary somatosensory cortex. *Nature* **436**, 261–265 (2005).

42. Lee, L.-J., Chen, W.-J., Chuang, Y.-W. & Wang, Y.-C. Neonatal whisker trimming causes long-lasting changes in structure and function of the somatosensory system. *Exp. Neurol.* **219**, 524–532 (2009).

43. Blake, D. T. *et al*. Sensory representation abnormalities that parallel focal hand dystonia in a primate model. *Somatosens. Mot. Res.* **19**, 347–357 (2002).

44. Wang, X., Merzenich, M. M., Sameshima, K. & Jenkins, W. M. Remodelling of hand representation in adult cortex determined by timing of tactile stimulation. *Nature* **378**, 71–75 (1995).

45. Riesenfeld, A. Metatarsal robusticity in bipedal rats. *Am. J. Phys. Anthropol.* **36**, 229–233 (1972).

46. Coq, J. O. *et al*. Peripheral and central changes combine to induce motor behavioral deficits in a moderate repetition task. *Exp. Neurol.* **220**, 234–245 (2009).

47. Bland, S. T., Gonzales, R. A. & Schallert, T. Movement-related glutamate levels in rat hippocampus, striatum, and sensorimotor cortex. *Neurosci. Lett.* **277**, 119–122 (1999).

48. Orczyk, J. J. & Garraghty, P. E. Reconciling Homeostatic and Use-Dependent Plasticity in the Context of Somatosensory Deprivation. *Neural Plast.* **2015**, e290819 (2015).

49. Pearson, J. C., Finkel, L. H. & Edelman, G. M. Plasticity in the organization of adult cerebral cortical maps: a computer simulation based on neuronal group selection. *J. Neurosci. Off. J. Soc. Neurosci.* **7**, 4209–4223 (1987).

50. Edelman, G. M. Neural Darwinism: selection and reentrant signaling in higher brain function. *Neuron* **10**, 115–125 (1993).

51. Hadders-Algra, M. The neuronal group selection theory: promising principles for understanding and treating developmental motor disorders. *Dev. Med. Child Neurol.* **42**, 707–715 (2000).

52. Purves, D. *Body and Brain: A Trophic Theory of Neural Connections*. (Harvard University Press 1988).

53. Leonard, C., Hirschfeld, H. & Forssberg, H. Gait acquisition and reflex abnormalities in normal children and children with cerebral palsy. In *Posture and Gait: Development, Adaptation, and Modulation: Proceedings of the 9th International Symposium on Postural and Gait Research, Marseille, France*, 1988 33–45 (Excerpta Medica 1988).

54. Clarac, F., Brocard, F. & Vinay, L. The maturation of locomotor networks. *Prog. Brain Res.* **143**, 57–66 (2004).

55. Canu, M.-H., Garnier, C., Lepoutre, F.-X. & Falempin, M. A 3D analysis of hindlimb motion during treadmill locomotion in rats after a 14-day episode of simulated microgravity. *Behav. Brain Res.* **157**, 309–321 (2005).

56. Gracies, J.-M. Pathophysiology of spastic paresis. I: Paresis and soft tissue changes. *Muscle Nerve* **31**, 535–551 (2005).

57. Gracies, J.-M. Pathophysiology of spastic paresis. II: Emergence of muscle overactivity. *Muscle Nerve* **31**, 552–571 (2005).

58. Phagava, H. *et al*. General movements in infants with autism spectrum disorders. *Georgian Med. News* 100–105 (2008).

59. Field, A., Miles, J. & Field, Z. *Discovering Statistics Using R*. (SAGE Publications, 2012).

## Acknowledgements

This work was inspired by Prof Michael M Merzenich and Dr Fabrizio Strata. The authors are grateful to Prof J-M Gracies, L Bensoussan and V Triffeau, and Drs H Bras, M Kerzoncuf and F Boubred for helpful discussions and comments on the manuscript. The authors also wish to dedicate this work to Vicky Massicotte and Laurent Vinay, great neuroscientists recently deceased. This work was supported by the Centre National de la Recherche Scientifique (CNRS), the National Institute of Mental Health of the National Institutes of Health (Award Number P30 MH092177), la Fondation Paralysie Cérébrale (previously la Fondation Motrice), la Fondation NRJ-Institut de France, Région Provence-Alpes-Côte-d'Azur, the Cerebral Palsy Alliance, Aix-Marseille Université and Temple University.

## Author Contributions

J.O.C. planned and contributed to all the experiments, performed all statistics and P.C.A., and wrote the manuscript; M.D. and M.R. contributed to all the experiments; N.T.L. performed *in vivo* microdialysis; FrC performed western blotting in the sensorimotor cortex; M.F.B. performed brain histology, musculoskeletal histopathology and immunohistochemistry; J.O.C., M.F.B., M.D., M.H.C. and FloC contributed to write, edit and revise the manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-34312-y.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_28.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018