Abstract
A long-standing question in systems neuroscience is to what extent task-relevant features of neocortical processing are localized or distributed. Coordinated activity across the neocortex has been recently shown to drive complex behavior in the mouse, while activity in selected areas is canonically associated with specific functions (e.g., movements in the case of the motor cortex). Reach-to-grasp (RtG) movements are known to be dependent on motor circuits of the neocortex; however, the global activity of the neocortex during these movements has been largely unexplored in the mouse. Here, we characterized, using wide-field calcium imaging, these neocortex-wide dynamics in mice of either sex engaging in an RtG task. We demonstrate that, beyond motor regions, several areas, such as the visual and the retrosplenial cortices, also increase their activity levels during successful RtGs, and homologous regions across the ipsilateral hemisphere are also involved. Functional connectivity among neocortical areas increases transiently around movement onset and decreases during movement. Despite this global phenomenon, neural activity levels correlate with kinematics measures of successful RtGs in sensorimotor areas only. Our findings establish that distributed and localized neocortical dynamics co-orchestrate efficient control of complex movements.
SIGNIFICANCE STATEMENT Mammals rely on reaching and grasping movements for fine-scale interactions with the physical world. In the mouse, the motor cortex is critical for the execution of such behavior, yet little is known about the activity patterns across neocortical areas. Using the mesoscale-level networks as a model of cortical processing, we investigated the hypothesis that areas beyond the motor regions could participate in RtG planning and execution, and indeed a large network of areas is involved while performing RtGs. Movement kinematics correlates mostly with neural activity in sensorimotor areas. By demonstrating that distributed and localized neocortical dynamics for the execution of fine movements coexist in the mouse neocortex during RtG, we offer an unprecedented view on the neocortical correlates of mammalian motor control.
Introduction
Deciphering how the CNS produces goal-directed movements, such as the reach-to-grasp (RtG), is crucial for expanding our understanding of the neural computations underlying motor ability and animal cognition in general (Bayne et al., 2019). RtG is known to be dependent on neocortical circuits, and intensive research has provided detailed insight into how neural activity from motor areas of the neocortex contributes to the control of movements (Evarts, 1968; Georgopoulos et al., 1982; Graziano et al., 2002; Churchland et al., 2012). In rodents, medial prefrontal regions, such as the secondary motor cortex (also termed medial agranular cortex, precentral cortex, or frontal orienting field), have been shown to be critical for movement preparation (N. Li et al., 2016; Barthas and Kwan, 2017; Chen et al., 2017). Nevertheless, how neural systems give rise to behavior cannot be fully comprehended from segregated analysis of its components. Indeed, parallel to the patterns of neural activity at the local level, information processing occurring at the corticocortical level is necessary for complex behavior (Peters et al., 2014; Allen et al., 2017; Battaglia-Mayer and Caminiti, 2019).
However, despite the growing adoption of the mouse for elucidating the neural underpinnings of motor control (Ölveczky, 2011), the neocortex-wide dynamics during RtG remain largely unexplored. Advancements in optical methods, especially wide-field fluorescence microscopy, have allowed monitoring neuronal activity from almost the entire dorsal neocortex in awake, behaving mice (Allegra Mascaro et al., 2019; Montagni et al., 2019; Sancataldo et al., 2019; Ren and Komiyama, 2021). With this tool, it has been recently shown that neocortex-wide dynamics emerge during learning (Makino et al., 2017); nevertheless, the contribution of regions beyond the contralateral motor areas and the corticocortical interactions during skilled movements remain largely unexplored. Also, while movement kinematics are associated with neural activity in the contralateral motor areas, and information about such activity can also be used to predict several features of the executed movement (Prsa et al., 2017; C. Li et al., 2019), information on the relationship between RtG kinematics and neural activity across the neocortex is still scanty.
Our results indicate that global cortical activation emerges during successful RtGs, with functional connectivity (FC) across areas peaking around movement onset, whereas movement kinematics correlated only with the activity in a portion of the sensorimotor region.
Materials and Methods
Transgenic mice
All experimental procedures were authorized by the Italian Ministry of Health (authorization 127-2018-PR). A total of 9 mice (3-11 months of age) were used in this study. Six (3 male, 3 female) were GCaMP mice (C57BL/6J-Tg(Thy1-GCaMP6f)GP5.17Dkim/J, stock #025393, herein referred to as GCaMP6f) (Chen et al., 2013), whereas 3 (2 females, 1 male) GFP mice (Tg(Thy1-EGFP)MJrs, stock #007788) were used as controls to account for potential hemodynamic artifacts. Animals were housed on a 12 h light/dark cycle with ad libitum water. Mice undergoing behavioral training were food-restricted to 80%-90% of original body weight by limiting food intake to 2-3 g/d. Animals were monitored and weighed daily, and food ratios were increased if necessary. Veterinary staff monitored animals twice a week.
Surgery
The day before surgery, mice were given a subcutaneous injection of enrofloxacin (10 mg/kg) and each day thereafter for 3 d, to prevent infections. On the day of surgery, animals were anesthetized with isoflurane (3% induction, 1.5% maintenance) and placed on a stereotaxic apparatus (Kopf Instruments). The absence of tail reflex and toe-pinch reflex was tested to confirm that the mouse was adequately anesthetized. Throughout the surgery, the temperature was kept constant at 37°C by means of an electric heating pad controlled by a rectal temperature sensor (Stoelting). To prevent cornea dehydration, ophthalmic gel was applied over the eyes. The scalp was first scrubbed with ethanol and betadine, then a depilatory cream was applied to remove hairs, and then a few drops of a topic agent (lidocaine hydrochloride 2%) were applied as an analgesic measure. A circular piece of scalp was removed, and the anatomic points of reference (bregma and λ) were stained with a permanent marker. Two semicircular coverslip glasses, one per hemisphere, were glued on the skull above the cortex using a transparent dental cement. A custom-built headpost was placed ∼0.5 mm posterior to λ. At the end of the surgery, subcutaneous injections of analgesic and anti-inflammatory (carprofen, 5 mg/kg) drugs were given to facilitate recovery. To this end, 200 ml lactated Ringer's solution of the 0.9% saline was given to each mouse at the end of surgery. Mice were allowed to recover at least 1 week following surgery. During this period, animals were given ad libitum food and water.
Behavioral apparatus
The experimental apparatus (see Fig. 1A) comprised a stage (size 300 × 300 × 20 mm, Thorlabs) and walls to form an enclosure and secured above a lab jack (MLJ050/M, Thorlabs) and was inspired by the work of Guo et al. (2015). The enclosure is equipped with a base for the animal, a perch, a speaker to deliver the auditory cue, and a turntable attached to a servo-motor equipped with a contactless magnetic rotary encoder to ensure consistent pellet positioning (MX-28AT) and controlled by a controller board (ArbotiX-M. Robocontroller, Trousser Robotics). The servo is controlled via an ad hoc program. Chocolate pellets (10 mg, 5TUL Purified 10 mg pellets, catalog 1811529, TestDiet) were used to induce the RtG movements, thus serving as stimuli. The interstimulus interval was 6-20 s. The enclosure was placed within a custom-made sensory-isolated enclosed imaging recording chamber that is kept dark and has been insulated with acoustic foam to further reduce ambient sounds (RS components and woods).
The behavioral data were collected at 200 frames/s via a high-speed camera (CM3-U3-13YC-CS, Chameleon3, FLIR) with 2.8-8 mm varifocal objective lenses (LENS-30F2-V80CS, Fujinon), placed to the left of the head-fixed animal. A 470 nm LED illuminating the stage was used as light source to create a high-contrast image for processing and quantification. Movies were recorded using the FlyCapture software (FLIR).
Wide-field microscopy
Imaging was performed through the intact skull using a custom-made microscope. The microscope consisted of back-to-back 50 mm f/1.2 camera lenses (Nikon), separated by a FF495-Di03-50.8-D dichroic mirror (Semrock), mounted in a 60 mm cube (Thorlabs). To excite the GCaMP6f indicator, a 470 nm light source (LED, M470L3, Thorlabs) was deflected by a dichroic filter (DC FF 495-DI02 Semrock) on the objective (TL2X-SAP 2× Super Apochromatic Microscope Objective, 0.1 NA, 56.3 mm WD, Thorlabs). The fluorescence signal was selected by a bandpass filter (525/50 Semrock) and collected by a CMOS camera (ORCA-Flash 4.0 V2, Hamamatsu). Images were acquired at 25 Hz, with a resolution of 512 × 512 pixels with a FOV of ∼12 × 12 mm (depth 16-bit) via the HCImage Live software (Hamamatsu). The microscope thus allows a FOV embracing the entire dorsal prospect of the mouse neocortex. To reduce unwanted light scattering on the mouse eyes (which could serve as visual stimuli, inducing unwanted neural activity), an iris (ID15/M, Mounted Standard Iris, Ø15 mm Max Aperture, TR75/M. Post, Thorlabs) was placed 1-2 mm above the mouse skull. Imaging sessions were performed for each behavioral training session. The synchronization of the components of the experimental setup occurred via a common hardware trigger signal at the start of each block (for the definition, see Experimental design).
Experimental design
Following recovery from surgery for at least 1 week, mice were habituated to the experimental setup for ≥4 d (≥30 min a day/mouse). Then, mice were gently placed under the microscope objective and fixed (see Fig. 1B). Each session contained at least 5 blocks, and each block lasted between 65 s and 6 min, 50 s. Here, we analyzed the first 10 d starting from when the mouse performed the first successful RtG. The neural dynamics associated with both successful movements and with errors trials were analyzed. Successful RtGs attempts were defined as RtGs that ended with the pellet eaten, error RtGs attempts were defined as RtGs with reach that did not end with the pellet eaten because the animal missed the pellet or failed the grasp or dropped the pellet before eating it. We analyzed the neural dynamics for a total of 298 movements (n successful RtGs = 78, n error RtGs = 117, n successful RtGs performed by GFP mice = 103). To ensure consistency, only trajectories from successful RtG performed by GCaMP mice were analyzed.
Image processing and data analysis
For the movement data, image sequences were loaded onto FIJI (ImageJ) and the metacarpus position was tracked using the manual tracking plugin, to obtain the XY coordinates of the RtGs. The timing of the movement onset (MoveON) was defined as the frame for which the operator detected a variation in the metacarpus position. The metacarpi coordinates were then used to reconstruct the movement trajectories and compute the movement kinematics. We computed the following kinematics parameters that were correlated with the ones used to characterize the neural activity (see below): mean speed (mm/s), max speed (mm/s), duration (s), pathlength (mm), number of speed peaks, mean acceleration (mm/s2), and max acceleration (mm/s2). The number of movement speed peaks was used as a measure of movement smoothness (Lai et al., 2015) and was defined as the number of peaks above the mean movement speed.
For the neural data, image stacks for each animal collected from different sessions were registered using custom-made software, by taking into account the bregma and λ position.
Raw wide-field imaging data of GCaMP fluorescence entail a non-negligible source of contamination because of hemoglobin, whose absorption spectrum largely overlaps with that of GCaMP, making interpretation of neural activity from these wide-field data challenging (Ma et al., 2016a). In order to account for the hemodynamic contribution, we first acquired widefield fluorescence data from the neocortex of GFP mice trained to execute RtGs. We expected the fluorescence signal from these mice to be static (i.e., not responsive to neuronal calcium levels); however, as reported elsewhere (Allen et al., 2017), fluorescence fluctuations were present in GFP mice around the movement onset (see Movie 3). We used the non-negative matrix factorization (NMF) as a source separation technique to demix the neuronal calcium signal from the non–calcium-related signal like hemodynamics. Among other source separation algorithms, we used the NMF because it has been previously used to calcium signals (Pnevmatikakis et al., 2016; Saxena et al., 2020) and because the weights and components are non-negative as opposed to other techniques, such as independent component analysis or singular value decomposition. For each animal, we concatenated the ΔF/F0 fluorescence sequences of any valid RtG considering a time window ranging from −4 to 6 s around the movement onset (250 frames for each considered RtG). To avoid values of ΔF/F0 <0, which are required by the NMF, for those trials whose ΔF/F0 was <0, we found the minimum value for the entire stack in the time window and offset the trial by the same value so that the minimum value becomes zero. This image stack was then reshaped into a matrix F whose columns represent the number of frames in the stack (250 by the number of considered RtGs) and its rows represent the number of pixels of the each frame (40,000 = 200 × 200). We then used the NMF to factorize the fluorescent stack of each animal into two matrices W and H by setting the numbers of components to 40.
Where W is the component matrix whose columns equal the number of components and the same number of rows of F which represent the pixels of each frame, H, instead, represents the weight matrix and represents the coefficient to assign to each component to reconstruct the original fluorescence stack. We then analyzed each component of the W matrix by reshaping it into a 200 × 200 pixels image. We considered components to be non–calcium-related like hemodynamics if: (1) the component highest activation matches the superficial vasculature observed from the raw fluorescence image; (2) the components highest activation lays on the edge of the mask where the skull starts to form the temporal bones; and (3) the component was considered to be noisy (i.e., activation was low and spread over large portion of the cortical mantle). We then demixed the neuronal calcium signal from the non–calcium-related signal by setting to 0 the columns of W that met one of the three abovementioned conditions.
To remove the offset introduced by summing the minimum before applying the NMF, we subtracted the ΔF/F0 in the time window −2 to −0.5 s. The signal was then low-pass filtered (9 Hz) to reduce the unwanted contribution of respiratory and heartbeat components. The comparison between the raw and the resulting denoised signal for the GCaMP and GFP mice is in Movies 2 and 3, respectively.
To dissect the contribution of each cortical area for RtG, the processed stacks we also registered the cortex to the surface of the Allen Institute Mouse Brain Atlas (www.brain-map.org) projected to our plane of imaging. We further applied a mask to exclude the medial-most areas because they lay over the superior sagittal sinus; therefore, the neuronal signal in this region could be more easily affected by hemodynamic artifacts. Areas laying on the most lateral parts of the mouse cortex were also excluded. This parcellation of the neocortex created 22 areas for each hemisphere, for a total of 44 areas, plus the left neocortical hemisphere, the right neocortical hemisphere, and both hemispheres (CTX). The abbreviations and extended names for each areas are as follows: MOs, secondary motor area; RFA, rostral forelimb area; MOp, primary motor area; CFA, caudal forelimb area; SSp-bfd, primary somatosensory area, barrel field; SSp-ll, primary somatosensory area, lower limb; SSp-m, primary somatosensory area, mouth; SSp-n, primary somatosensory area, nose; SSp-tr, primary somatosensory area, trunk; SSp-ul, primary somatosensory area, upper limb; SSp-un, primary somatosensory area, unassigned; SSs, secondary somatosensory area; RSPagl, retrosplenial area, lateral agranular part; RSPd, retrosplenial area, dorsal part; VISrl, rostrolateral visual area; VISa, anterior visual area; VISal, anterolateral visual area; VISam, anteromedial visual area; VISl, lateral visual area; VISli, laterointermediate visual area; VISp, primary visual area; and VISpm, posteromedial visual area. The coordinates for CFA (AP: 0.2 mm; ML: 1.7 mm) and RFA (AP: 2.3 mm, ML: 0.9 mm) were determined herein as the average of the coordinates reported in recent papers (Tennant et al., 2011; Hira et al., 2013; Morandell and Huber, 2017; Wang et al., 2017). Throughout the text and figures, we added the suffix L and R to term cortical areas of the left (ipsilateral with respect to the paw executing the movement) or right (contralateral with respect to the paw executing the movement), respectively (e.g., MOs_L, MOs_R). For each block, the image stacks were processed to obtain the estimates of ΔF/F0. Briefly, the image stacks were processed considering the equation ΔF/F0 = (F – F0)/F0, where F defines the value of the fluorescence signal in a given moment and F0 defines mean fluorescence over each block. Neocortical areas were considered associated with movement if their activity during any moment of the time window crossed a threshold defined as the mean + 1.96 × SD of the neural signal occurring during the baseline (−2 to −1 s with respect to MoveON). The highest value reached at any moment in the time window −2 to 1.5 s around MoveON was defined as the PeakMax. Also, to statistically determine whether the transaction between the behavioral phases is associated with different neural activity levels, we used a two-way mixed-model ANOVA followed, if appropriate, by Tukey's HSD post hoc tests. To determine the temporal involvement of neocortical areas during the task, we computed the latency necessary to reach the 50% and the 100% of the PeakMax (latency to HalfPeakMax and latency to PeakMax, respectively) in each of the responsive areas. Finally, to quantify the rate of rise of the neural response in each area, we quantified the slope of the neural response as
To evaluate the FC between the areas during RtGs, cross-correlations between pairs of areas for each behavioral phase were computed at the single-trial level. The fluorescent time series were chunked according to the behavioral phase of interest and are presented as cross-correlation matrices (CCMs). We investigated the network membership of these regions by using a hierarchical clustering algorithm (HCA) that was based on Ward's linkage method (Ward, 1963), which minimizes the variance between the clusters. The proximity was then interpreted as an indirect measure of FC and represented graphically by means of dendrograms. All data are reported as mean ± SD if not stated otherwise. Sample size and appropriate statistical analyses are specified in each figure legend. Statistical significance was defined as α < 0.05 if not stated otherwise. Multiple comparisons were corrected by Bonferroni or Bonferroni–Hochberg procedures. No statistical methods were used to predetermine the sample size. All statistical analyses and related figures were created using either Python or R and were assembled in Inkscape.
Results
Wide-field calcium imaging during RtG
We conceived that, for a mouse to perform complex, goal-directed movements, such as RtG, a specific neural network within the neocortex could be recruited. To gain information about such dynamics, our setup integrates a large FOV microscope to simultaneously image the entirety of the neocortical surface and a behavioral apparatus for inducing, delivering, and videotaping the RtGs (Fig. 1A). Mice were trained to perform the task while head-fixed under the microscope (Fig. 1B). Once trained, animals executed RtG movements to eat chocolate pellets that were placed on a rotating table placed in front of them (Fig. 1C, left). The time-aligned image stacks obtained from the camera and from the microscope allowed the monitoring of neural activity across functional areas of the neocortex, throughout the behavioral task, with activity increasing in contralateral regions, and to a lesser extent for the ipsilateral sensorimotor regions (Fig. 1C,D; Movie 1).
Neural activity increases across several areas spanning both hemispheres during successful RtG behavior
Single-pixel analysis of cortex-wide activity indicates that, on average, for all successful RtG movements (Fig. 2) executed by all subjects, the neural activity increases bilaterally in the neocortical mantle. The contralateral hemisphere displayed higher activity compared with the ipsilateral one only in a small portion of the cortex; at the rostral portion of the secondary motor cortex, and in the sensorimotor cortex with peak activity centered around the upper limb somatosensory cortex (Fig. 2E). We averaged the different contributions of neocortical activation during four temporal windows (i.e., at 2-1.5 s and at 0.5-0.04 s before MoveON, at 0-0.5 s, and at 0.5 + 1.0 s after MoveON; Fig. 2A,B) showing that activity increases in the −0.52 to 0.04 s window and lasts up to 1 s. In particular, premovement activity irradiated from medial to lateral regions of both hemispheres to (Fig. 2D) and during movement, bilateral activation in more rostro-medial areas becomes apparent with more pronounced activity in the contralateral hemisphere ∼0.24-0.48 s after movement onset (Fig. 2A). Instead, later (>0.48 s) phases of movement activity are characterized by an increase in neural activity in the rostro-lateral regions of both hemispheres. When considering all windows together, it is clear that neural activity for successful RtGs spreads over a wide portion of the cortical mantle.
When considering unsuccessful RtGs, we observed a different activation pattern with respect to successful ones (Fig. 2F,G). In particular, neural activity during errors diverged from correct trials starting at 0.38 s and is characterized by low values of ΔF/F0. In addition, even if activity originated bilaterally in medial regions as in successful RtGs, the neural response was shorter and did not spread to the rostro-lateral regions as in successful RtGs.
To characterize the contribution of functional areas across the neocortex during successful RtG, we parcellated the acquired image stacks according to the Allen Institute Mouse Brain Atlas (Fig. 3A, left). The resulting parcellation allowed to extract the calcium signal for each neocortical area. The increase in the levels of neural activity encompassing the neocortical mantle is observed also on each hemisphere separately (Fig. 3B). Several areas displayed a vigorous increase in neural activity levels with movement execution, including the motor areas and the retrosplenial cortices of both hemispheres, whereas other areas, such as the SSp-bfd, displayed a more modest increase in activity (Fig. 3C,D). Overall, the activity levels across homologous areas of both hemispheres, including motor, somatosensory, retrosplenial, and visual areas, appeared comparable.
Calcium signals from neocortical areas increase before movement onset
We compared the temporal profile of the parcellated fluorescent calcium signal across the peri-event time window. Of the areas considered, all (n = 44) emerged as associated with movement (i.e., showed activity levels higher than the set threshold) with the primary and secondary motor areas displaying the highest PeakMax and other areas spread across both hemispheres, forming a distributed cortical network associated with RtG movements (Fig. 4A). Across all responsive areas, the HalfPeakMax was reached before movement onset, which indicates that the increases in the amplitude levels of the fluorescent signal anticipate the movement onset (Fig. 4B, left). In addition, across all responsive areas, we found that the first time point passing the threshold occurred well before movement onset (580 ± 113 ms). These results suggest that the temporal dynamics of the calcium signal detects variations in the levels of neural activity also in the RtG premovement phase. Moreover, the latency to reach the HalfPeakMax showed low variability across areas. In contrast, the latency to reach PeakMax occurred after the movement onset for all areas and differed across areas (Fig. 4B, right). In particular, somatosensory areas and the contralateral CFA reached the PeakMax later, compared with other areas. Next, we compared the neural activity levels across the three temporal windows (rest: −2 to 1.5 s; premove: −0.5 to 0.04 s; move: 0 to 0.5 s), across neocortical areas (Fig. 4C). We found that both the time window and the areas significantly affect the calcium responses [two-way mixed model ANOVA, effect of temporal window, F(1.89,6392.37) = 2824.784, p < 0.0001, temporal window × area interaction, F(81.13,6392.37) = 7.786, p ≤ 0.0001, Bonferroni post hoc tests]. Post hoc pairwise comparisons indicate statistically significant increases from rest to movement for several areas, including the MOs, the RFA, the MOp, the CFA, the RSPd, the SSp-ul, and the visual areas. These increases were found to be already present when comparing the rest versus the premove window (e.g., for the MOs, the RFA, the MOp, the CFA, the RSPd, and the SSp-ul) but were found not to be different for few areas, including the SSp-bfd and the several visual areas (for all pairwise comparisons, see the Extended Data Fig. 4-1). These findings suggest that a global network of areas is activated during movement execution, with the majority of them increasing their levels of activity before movement onset.
Figure 4-1
Pairwise correlations of neural activity levels across temporal windows for all neocortical areas. Download Figure 4-1, XLS file.
Movement onset transiently increases FC
Considering the increase in neural activity levels across the areas from rest to movement, we next asked whether such variation was associated with a modulation of FC across areas and whether FC was modulated throughout the behavioral phases. The hierarchical cluster analysis (HCA) revealed that near movement onset there is an increase in FC. Indeed, we found that the Ward distance was lowest around MoveON (−0.28 s, time window length = 0.48 s; Fig. 5, middle, left inset), indicating that the areas tend to form a single functional aggregate. This effect is illustrated by the CCM (Fig. 5, middle, right), indicating an overall increase in the correlation coefficients around MoveON. Conversely, during movement execution (“Grasp”), the FC decreases (Fig. 5, bottom, right), leading to the emergence of four distinct clusters (Fig. 5, bottom, left). The first cluster contained the caudal portion visual areas only and therefore could be regarded as a visual network. The second cluster included intermingled areas: that is, two subdivisions of the visual cortex (VISam and VISa), the retrosplenial cortices (RSPd and RSPagl), several somatosensory areas, as well as the CFA thus could be considered as a multimodal network. The third cluster was formed by two areas of the primary somatosensory cortex related to the mouth and nose representation (i.e., the SSp_m and the SSp_n). Finally, the last cluster is constituted by two somatosensory areas, including the one associated with the representation of the upper limbs (SS-p_ul) as well as an adjacent area with an unknown association (SS-p_un) as well as the primary and secondary motor areas and the RFA and could be regarded as the sensorimotor network.
Because the used parcellation may not account for subfield-specific FC, we next sought to determine whether such a scheme of communication is maintained in the context of single-pixel analysis. To this end, we used seeds of the contralateral hemisphere that affiliated with the four clusters emerging during the grasping phase to compute the single pixel seed correlations shown in Figure 5. Overall, the seed correlations are in good agreement with the results of the HCA, and both show that the transition from resting to movement execution was associated with a robust increase in the cross-correlations around movement onset (Fig. 5, right).
Our approach provides evidence for the existence of segregated functional networks emerging during motor execution. Overall, these FC results suggest that the corticocortical communication scheme during goal-directed forelimb movements can be decomposed into two parts: the earlier phase, with a transient increase in cross-correlation; and a latter phase, where separated networks emerge.
Movement kinematics correlates with neural activity levels in sensorimotor areas of both hemispheres
Considering the increase in activity observed across the neocortical areas, we sought to determine the extent of association between neural activity and movement parameters. To this end, we first reconstructed the spatial trajectories for the successful RtGs and computed the movement metrics (Fig. 6A–C). The movement distance was 54.54 ± 18.81 mm (Fig. 6D), the movements lasted 0.71 ± 0.26 s (Fig. 6E), while the number of speed peaks was 17.78 ± 7.99 (Fig. 6F).
Since during movement execution we observed the emergence of differentiated activation among areas, we asked to what extent each area contributed to specific kinematic parameters of RtG movement. To this end, we performed correlation analyses between neural activity measures in all responsive areas of the neocortex and the extracted parameters of movement kinematics (Fig. 7). Differently from the activity levels, which were broadly increased across the neocortical mantle, statistically significant correlations between neural activity and movement kinematics were found only for few areas across both hemispheres (Fig. 7A). Specifically, our analysis revealed a negative correlation between the amplitude of the signal (PeakMax) and movement smoothness in the MOp, the CFA, the SSp-ul, and the SSp-un of the ipsilateral hemisphere (number of speed peaks, ρ = −0.40, padj < 0.01; ρ = −0.49, padj < 0.01; ρ = −0.41, padj < 0.01; ρ = −0.39, padj = 0.02 for MOp, the CFA, the SSp-ul, and the SSp-un, respectively; Fig. 7B). Also, in the ipsilateral CFA, the PeakMax was inversely correlated with the movement pathlength (Fig. 7D). In the contralateral RFA, we found that the slope of the neural response was inversely correlated with movement pathlength (ρ = −0.54, padj < 0.01; Fig. 7C). In this area, the latency to reach the PeakMax correlated with movement smoothness (ρ = 0.46, padj < 0.01) and pathlength (ρ = 0.43, padj < 0.01; Fig. 7D), while the slope was inversely correlated with movement smoothness (ρ = −0.51, padj < 0.01; Fig. 7D). Similarly, the latency to reach the PeakMax in the contralateral MOp and in the CFA was anticorrelated with respect to movement speed (ρ = −0.47, padj < 0.01 and ρ = −0.43, padj < 0.01, for MOp and for the CFA, respectively; Fig. 7C). At variance with all the other areas, the degree of movement smoothness was positively correlated with the latency to reach the PeakMax in the contralateral SSp-ll (ρ = 0.41, padj <0.01). Therefore, while neural activity level increases neocortex-wide, a restricted group of rostral areas exhibit activity that correlates with kinematic parameters of RtGs.
Discussion
We provided an unprecedented view of the mesoscale level neocortical dynamics from mice performing an RtG task. To the best of our knowledge, this is the first reported step toward an optical characterization of neocortex-wide neural dynamics in rodents during voluntary RtG movement planning and execution. Our work provides novel evidence for the role played by the secondary motor areas during RtG and, by directly observing all areas of the neocortex, an integrative framework not reported previously.
Our results point the attention also to areas beyond the MOs, the MOp, and their subdivisions (RFA and CFA) toward other territories, such as the RSC, the SSp, and the visual areas, which were previously not recorded and studied simultaneously within the context of the RtG. For the RSC, while its participation to RtG in the mouse is demonstrated here for the first time, the contribution of RSC may be explained by its known involvement in the processing of visual information necessary to locate the target relevant for the behavioral tasks (Czajkowski et al., 2014). Indeed, RSC is involved in early phases of learning (Makino et al., 2017); however, no information was available about the corticocortical interactions (e.g., with the secondary motor areas) and about the relationship with kinematic parameters of movement.
The global activation found in the present work is consistent with recent reports on other types of behavioral tasks (Goard et al., 2016; Wekselblatt et al., 2016; Allen et al., 2017; Kyriakatos et al., 2017; Makino et al., 2017; Musall et al., 2019; Stringer et al., 2019; Salkoff et al., 2020); here we extended it to the RtG. In addition, we show that this activation is supported at the trial-by-trial level by an increase in FC occurring at movement onset. Indeed, our HCA revealed that such an increase in FC represents the transition between the premovement and postmovement onset FC clusters whose correlation maps have been observed previously in the mouse (Vanni et al., 2017; Hakon et al., 2018; Clancy et al., 2019) and closely related to the underlying structural connectivity (Hakon et al., 2018). Therefore, such widespread activation would raise the possibility that neocortical activation patterns for the production of successful RtG movements are distributed rather than localized.
In interpreting our results, we cannot exclude that whole-body movements may account for part of the widespread distribution of neocortical activity we observe during RtG. However, it is important to note that, because of the trial averaging around the movement onset, such movements have to be performed by each animal in a substantial number of trials to contribute to the spread of the response. Also, other parameters, such as arousal, could possibly contribute to the increased activity we observed in visual areas. Indeed, it was shown that the levels of activity in the visual cortex are related to arousal (Stringer et al., 2019).
Recent work emphasizes that it is critical to understand both functional specialization and functional integration schemes and how they relate to different aspects of complex behavior (Cohen and D'Esposito, 2016). Our work thus shows the first direct evidence indicating that, while neural activity increases globally during RtG, neuro-kinematic correlations are significant only for a subset of rostral areas. The present findings are in agreement with Salkoff et al. (2020), in which it was shown that these areas encode performance in a “simple” licking behavior. The present results suggest that forelimb movement planning and execution may involve large neural networks across the neocortex while the encoding of kinematic parameters occurs primarily in segregated areas. However, being correlational in nature, our work will be likely instrumental for guiding optogenetic interventions, to gain insight about the causal role of each neocortical area during voluntary movements. Nonetheless, such direct observation of global and local activation patterns for dexterous movements provides an integrated view of the nature of information processing in the mouse neocortex and encourages future investigations toward decoding of behavioral parameters from neural data, a computational approach that is increasingly used on fluorescent imaging data of neural activity (Prsa et al., 2017; C. Li et al., 2019; Salkoff et al., 2020).
Our data add experimental evidence to the relatively new idea of ipsilateral representation of movement parameters, supporting recent work in rodents, where recordings restricted to the posterior parietal cortex recently reported evidence for ipsilateral control of limb movement (Soma et al., 2019). Indeed, ipsilateral activity is also a computational feature of the human neocortex and is increasingly getting attention in the motor control community (Bundy et al., 2018; Bundy and Leuthardt, 2019).
While the correlation between movement kinematics and ipsilateral activity during RtG came as a surprise, there is a growing literature suggesting that neocortical dynamics concerning the generation of contralateral and ipsilateral limb movements are controlled by both hemispheres (Mutha et al., 2013; Schaffer and Sainburg, 2017; Heming et al., 2019). In this regard, our results on both contralateral and ipsilateral activation and correlations with movement kinematics for RtG control support this bihemispheric model of forelimb control.
Our findings extend these studies by providing evidence for a correlation between large scale, area-level activity, and specific metrics of RtG movements. Our findings also encourage to study, at the cellular level, the role of areas found to be active during movement (e.g., the visual areas) less intensively studied in the association to forelimb movements (Morandell and Huber, 2017; Ebina et al., 2018; Galinanes et al., 2018). In the future, it would also be relevant to assess the neocortical dynamics during RtG in a freely moving condition or during performance refinement (Makino et al., 2017; Whishaw et al., 2017; Bollu et al., 2019; Hwang et al., 2019).
Moreover, movement errors are associated with altered activity in the motor and parietal cortex (Diedrichsen et al., 2005); however, information on neocortex-wide dynamics of error movements is lacking in the mouse. Our wide-field imaging approach suggests that the widespread neural activity observed during successful RtGs is not generic to movement execution but rather the success of the completed goal-directed action.
There are some aspects that need to be considered when interpreting our results. The one-photon, wide-field calcium imaging approach is sensitive to light scattering potentially leading to artifacts in the signals extracted from the cortical areas, which could increase the risk of erroneously treating an area as movement-associated (Yang and Yuste, 2017). Such possibility is reduced by the evidence that areas even in close proximity (e.g., the barrel fields vs the upper limb portion of the somatosensory cortex) do not display the same levels of activity (Fig. 4A). We concede that the number of neocortical clusters emerging during movement could be influenced by the parcellation of the neocortex, which may introduce spurious segregations beyond those that are physiologically present. Importantly, the overlap between the results obtained from the seed-pixel analysis and the CCM supports the existence of at least four distinct neocortical networks during voluntary motor control. The analyzed neural activity is inferred from the fluctuations of fluorescence intensity throughout time, a signal composed also by non-neuronal sources, notoriously of hemodynamic origin, which represents another artifact that needs to be taken into account (Ma et al., 2016a,b). However, the hemodynamic component of the signal has a much slower rise and decay time with respect to the neuronal signal. This temporal feature of the hemodynamic component of the GCaMP6f signal reduces its contribution relative to the results reported here (Scott et al., 2018). In addition, our control experiments on GFP mice allowed us to obtain signal components that were exclusive for GCaMP mice thus further mitigating the contribution of non-neuronal signals in our results.
The transgenic line used here reports variation in the activity of excitatory neurons only, leaving the question open as to what is the role played by other neural populations, considering, for instance, that GABAergic neurons in the primary motor cortex signal the initiation of voluntary reaching movements (Estebanez et al., 2017).
The temporal kinetics (50 ms rise time and 150 ms decay time) intrinsic to the optical probe used here (GCaMP6f) (Chen et al., 2013; Dana et al., 2014) entails that the neural dynamics faster than 50 ms associated with RtG were not investigated in the present study. Moreover, because of the indicator decay time, it is hard to disentangle events that occur within a time window of 150 ms. This is also the reason why we do not extend our analysis past the peak in the calcium wave for the correlation analysis. Finally, the activity patterns observed here refer mostly to relatively slow (<13 Hz) cortical communication schemes; nonetheless, neural activity in the lower frequency range (including frequencies <4 Hz) encodes several parameters of movement kinematics (Bansal et al., 2011). The adaptation for such tasks for head-fixed mice is more recent (Guo et al., 2015); and recently, the work of the same group revealed how input from the thalamus drives activity in the primary motor cortex and RtG movement (Sauerbrei et al., 2020).
To overcome the temporal constraints associated with calcium imaging, a crucial step forward will be represented by the adoption of genetically encoded voltage indicators for the study of neocortex-wide dynamics during skilled movements, as very recent technological advances offer sufficiently stable and bright signal-to-noise ratio for in vivo applications (Abdelfattah et al., 2019; Adam et al., 2019; Villette et al., 2019). Even considering such limitations, our findings remain unexpected and provide experimental evidence to put into question some dogmas of the classical literature in motor control (Omrani et al., 2017), at least when investigating the modus operandi of the rodent neocortex. The present neocortex-wide characterization of the neural activity patterns associated with RtGs is a step toward all-optical interrogation of neural circuits associated with this behavior.
In conclusion, our investigation identified a global network of neocortical areas beyond those previously known to be associated with goal-directed arm movements and suggest that the operating mode of the mouse neocortex underlying unilateral goal-directed movements involves distant communication across both hemispheres, whereas the encoding of movement parameters occurs locally, and, in the context of our investigations, mostly in the sensorimotor areas.
Footnotes
This work was supported by European Union's Horizon 2020 Research and Innovation Program Grant Agreement 785907 (HBP SGA2); Italian Ministry for Education, University, and Research (Project FARE MIMIC); and European Research Council (European Union's Horizon 2020 Research and Innovation Program Grant Agreement 692943). We thank the staff of the mechanical and electronic workshops of European Laboratory for Non-Linear Spectroscopy for assistance; members of the F.S.P. laboratory for helpful discussions and suggestions; Costanza Campaioli and Chiara Caldini for help in developing earlier versions of the experimental setup; and Prof. Matteo Caleo for useful comments on the manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Francesco Saverio Pavone at pavone{at}lens.unifi.it