Peripheral Somatosensory Stimulation Elicits Traveling Waves with State-Dependent, Opposite Trajectories.

Optical imaging was performed on n = 10 Thy1-GCaMP6f mice fitted with cranial windows, permitting observation of cortical calcium dynamics through the intact skull as previously described (19). Imaging data were collected during 5 min of right whisker stimulation (Fig. 1A) using computer-triggered air puffs in 10 s blocks (5 s of puffs at 2 Hz, followed by 5 s of rest). Mice were then anesthetized with ketamine/xylazine, and the same stimulation protocol was repeated. Evoked responses were block averaged for awake (n = 10 mice, 600 blocks) and anesthetized (n = 9 mice, 540 blocks) conditions (Movie S1).

Fig. 1.

Peripheral somatosensory stimulation elicits traveling waves with state-dependent, opposite trajectories. (A) A diagram of cortical field of view spanning frontal/motor (red), somatosensory (yellow), and visual (blue). Right whisker stimulation elicits activity in S1_W (solid yellow) and M1_W (solid red). (B) Still frames depicting group-averaged response to a single right whisker deflection during wake (n = 10 mice, 600 blocks) and ketamine/xylazine anesthesia (abbreviated “Ket,” n = 9 mice, 540 blocks). The frames come from corresponding Movie S1. (C) A group-averaged time series of global GCaMP activity. The vertical lines indicate individual puff stimuli.

Fig. 1B shows the group-averaged response to a single puff in each condition. In awake animals, right whisker stimulation elicits activity focally in the left S1 whisker barrel cortex (S1_W) followed by a traveling wave spreading rostrally from S1_W up to the whisker motor cortex, consistent with prior reports employing voltage-sensitive dye imaging in awake mice in a single hemisphere (17, 20). By expanding the field of view to include both hemispheres, we further demonstrate that this broad, rostrally propagating wave is mirrored in the contralateral sensorimotor cortex in awake animals (Fig. 1B). Under ketamine/xylazine anesthesia, whisker stimulation elicits a similar local response bilaterally in the sensorimotor cortex that is promptly subsumed by a global SW propagating caudally to the visual cortex. Block-averaged time series reveal that activity reverberates for several seconds after stimulation ceases in anesthetized but not awake animals (Fig. 1C and Movie S1). This suggests that spontaneous global SWs are briefly phase-reset by preceding stimulation, similar to prior reports (21, 22).

We characterized the spatiotemporal structure of cortical activity by examining dominant patterns of global coherence (SI Appendix, Fig. S1). We applied spatial singular value decomposition (SVD) to the frequency domain representation [i.e., space-frequency SVD (23)] of the group-averaged responses, yielding a set of orthogonal spatial modes containing global magnitude and phase information at each frequency (see Materials and Methods). Consistent with the patterns observed in Fig. 1B, opposing wave trajectories in wake (back-to-front) and anesthesia (front-to-back) are evident from the phase maps of the dominant spatial mode at 2 Hz in each condition (SI Appendix, Fig. S1A). Notably, the state-dependent changes in rostral/caudal directionality of stimulus-evoked waves mirrors previously reported findings in spontaneous activity during wake and anesthesia (24).

Rhythmic Stimulation Shifts the Dominant Frequency of SW Activity Globally Across the Dorsal Neocortex.

The emergence of broadly propagating, stimulation-evoked SWs suggests that rhythmic somatosensory input may induce resonant activity globally across the cortex under anesthesia. To investigate this further, we asked whether stimulation shifts the dominant frequency of cortical activity and to what extent stimulation-evoked activity is globally coherent.

We first examined spectral power in whole runs of alternating 5-s blocks with stimulation ON or OFF to capture both spontaneous and evoked activity. As in local field potential recordings (25), spontaneous resting-state GCaMP dynamics exhibit a 1/f power spectrum. Under ketamine anesthesia, a prominent spectral peak emerges at 1- to 1.5-Hz that corresponds to spontaneous front-to-back SWs (19, 24). Here, we likewise observe that spontaneous 1- to 1.5-Hz SWs occur under anesthesia with front-to-back topography during interstimulus OFF periods (SI Appendix, Fig. S1B).

Importantly, by stimulating at a frequency faster than the endogenous SW frequency, we are able to dissociate unique properties of spontaneous SWs emanating from the frontal cortex (1 to 1.5 Hz, SI Appendix, Fig. S1B) from evoked SWs initiated in S1 (2 Hz, SI Appendix, Fig. S1A), which may be coordinated by distinct pathways. In order to determine whether the dominant frequency of global SW activity changes during ON periods of rhythmic stimulation, we computed the global average brain signal within each mouse and generated power spectrograms using a short-time Fourier transform (STFT). We observed that, under ketamine/xylazine anesthesia, global dynamics alternate between primarily 1- to 1.5-Hz oscillations during 5-s stim-OFF periods (corresponding to spontaneous SWs) and 2-Hz oscillations during 5-s stim-ON periods (corresponding to the imposed stimulation rhythm) (Fig. 2A). Averaging STFT spectra across ON and OFF blocks in anesthetized animals (Fig. 2B), we observed that OFF block power spectra look identical to those we have previously reported for spontaneous activity under anesthesia (24), while stimulation during ON blocks elicits a relative increase in global power at 2 Hz (***P = 0.0001) in tandem with a significant decrease in 1- to 1.5-Hz power (****P = 6.7e-5). Thus, during rhythmic stimulation under anesthesia, spontaneous SW power is suppressed while evoked power at the imposed stimulation frequency is enhanced. In contrast, in the awake state when SWs are absent, 2-Hz global power increases (***P = 0.0001), but 1- to 1.5-Hz power is not significantly changed (^NSP = 0.079).

Fig. 2.

Rhythmic stimulation shifts the dominant frequency of SW activity globally across the dorsal neocortex. (A) Group-averaged spectrograms of global GCaMP power (0 to 8 Hz, y-axis) over the course of consecutive 10-s blocks consisting of 5-s stim ON and 5-s stim OFF (time in 5-s bins, x-axis). (B) Average power spectra for all stim-ON (magenta) and stim-OFF (black) periods, computed by averaging across all ON or OFF windows from spectrograms in A. The differences are statistically tested by paired t test. Note, spectral peaks at ∼4 Hz (wake) and ∼5.5 Hz (ketamine) correspond to heartrate. The spectra are not variance normalized to permit comparison of changes in absolute power. The broad peak widths reflect biological variance, as well as the brief period for computing the Fourier transform. (C) Cortical maps depicting spatial distribution of 2-Hz power, computed from the group-averaged block. (D) The global coherence of evoked (2-Hz) and spontaneous (1- to 1.5-Hz) activity, comparing wake versus ketamine. The global coherence is calculated within individual mice as the fractional power of the leading spatial mode obtained from space-frequency SVD (23, 26). The statistical significance is determined by paired t test.

Next, we examined whether stimulation-evoked resonance spreads more globally under anesthesia. First, we computed pixelwise power spectra over block-averaged time series and observed that, under anesthesia, evoked 2-Hz power is indeed globally distributed across the whole cortex. In contrast, in awake animals, 2-Hz power is largely restricted to the whisker barrel cortex and motor cortex (Fig. 2C). Similarly, the magnitude topography of the dominant spatial mode derived from space-frequency SVD shows more global effects under anesthesia (SI Appendix, Fig. S1A). To quantify this phenomenon further, we calculated global coherence of evoked (2-Hz) and spontaneous (1- to 1.5-Hz) SWs as the fractional power of the leading spatial mode obtained from space-frequency SVD [i.e., the proportion of global coherence explained by the first SVD component (23, 26)]. Indeed, global coherence of evoked 2-Hz waves is significantly greater under ketamine anesthesia than during wake (**P = 0.0028, Fig. 2D). As expected, with the emergence of spontaneous global SWs, the leading spatial mode of 1- to 1.5-Hz activity is more broadly distributed in magnitude maps (SI Appendix, Fig. S1B) and more globally coherent under ketamine anesthesia than in the awake state (Fig. 2D, ***P = 0.0005).

Unilateral Ablation of S1 Reduces the Global Coherence of Somatosensory-Evoked SWs.

The above-described results demonstrate that unilateral somatosensory stimulation under anesthesia elicits a focal response in S1 followed by a global, front-to-back SW. This sequence is consistent with several possibilities and does not necessarily indicate causality for S1 triggering the global SW (15, 27). If the global SW component is triggered in parallel with S1 activation via direct input from subcortical circuits onto frontal/prefrontal regions, one may expect the global SW component to be independent of S1 activation. Alternatively, the front-to-back SW may be initiated via direct or indirect connectivity between the frontal cortex and the S1 contralateral to stimulation. Lastly, given that unilateral stimulation elicits bilateral responses in S1 at baseline (albeit more weakly in ipsilateral S1), a third possibility is that paired activation of both somatosensory cortices—likely involving their recurrent interaction—is necessary for triggering global SWs.

To disambiguate these possibilities, we analyzed a previously collected dataset in which responses to 3-Hz forepaw stimulation under ketamine/xylazine anesthesia were examined before (“pre”) and 1 wk following (“post”) unilateral laser photothrombosis of forepaw sensory cortex (S1_FP) (28). These infarcts have well-circumscribed margins and are confined within the left hemisphere S1_FP cortex (Fig. 3A) with an estimated volume of 1.94 ± 0.33 mm³. Imaging was performed during 5 min of stimulation of the right or left forepaw, under anesthesia, using mild cutaneous electrical shock in 20-s blocks (5 s of shocks at 3 Hz, followed by 15 s of rest). Results represent an average of n = 12 mice, 180 blocks total per paw at each timepoint.

Fig. 3.

Unilateral ablation of S1 reduces global coherence of 3-Hz somatosensory-evoked SWs but not 1- to 1.5-Hz spontaneous SWs. (A) A diagram depicting photothrombosis of the left forepaw somatosensory cortex (S1_FP, plotted bilaterally in solid yellow) along with a Nissl-stained cross-section through the middle of the infarct used to calculate infarct volume. Cell death is isolated within S1_FP cortex. (B) The still frames depicting group-averaged response to right paw stimulation pre- (light green) and poststroke (dark green) as well as left paw stimulation pre- (pale purple) and poststroke (dark purple). The data represent an average of n = 12 mice, 180 blocks total per paw at each timepoint. The frames come from corresponding Movie S2 (right paw stimulation) and Movie S3 (left paw stimulation). (C and D) Dominant spatial mode phase maps of evoked (3-Hz) and spontaneous (1- to 1.5-Hz) SWs derived from space-frequency SVD of group-averaged blocks of left paw stimulation (C) and right paw stimulation (D). The early portions of global waves exhibit negative phase values and late portions exhibit positive values (units are in radians). (E and F) The global coherence of evoked and spontaneous SWs during runs of left (E) or right (F) paw stimulation, comparing pre- and poststroke within individual mice with a paired t test. The quantification of global coherence excluded the infarct (masked in black) to capture changes in surviving cortex.

At baseline, 3-Hz electrical stimulation of the right or left forepaw elicits contralateral S1_FP and motor responses followed by a global, front-to-back SW (Fig. 3B and Movies S2 and S3), similar to responses following whisker stimulation (Fig. 1). One week after left S1_FP stroke, right forepaw stimulation on average evokes a weaker left S1_FP response, as expected due to local tissue death (visible as a dark circular void in Fig. 3B). Furthermore, the subsequent global wave of activity evoked by right paw stimulations is also dramatically attenuated: secondary activation of the ipsilateral motor cortex and contralateral sensorimotor cortex are weaker, and the ensuing global front-to-back SW is reduced in amplitude.

Stimulating the left forepaw after left S1_FP stroke still elicits a similar amplitude response in the right S1_FP and motor cortex in the contralesional hemisphere on average; however, surprisingly, the subsequent global SW component is again weakened (Fig. 3B). While the dominant spatial modes exhibit similar propagation patterns pre- and poststroke (Fig. 3 C and D), global coherence of evoked SWs was significantly reduced after stroke for both contralesional (left paw *P = 0.02, Fig. 3E) and ipsilesional (right paw *P = 0.03, Fig. 3F) stimulation. Importantly, these changes are specific to somatosensory-evoked SWs. In contrast, global coherence of spontaneous SWs was strong and not significantly altered from prestroke baseline during OFF periods of left paw (^NSP =0.85, Fig. 3E) and right paw (^NSP = 0.15, Fig. 3F) stimulation runs. Thus, unilateral ablation of S1 reduces global coherence of somatosensory-evoked SWs but not spontaneous SWs.

Unilateral Ablation of S1 Disrupts Globally Resonant Evoked SWs in Both Hemispheres.

Next, we tested whether the spatial distribution of SW resonance is altered following unilateral stroke in S1_FP. As in Fig. 2C, we mapped group-averaged, 3-Hz evoked SW power during stimulation to the left (Fig. 4A) or right forepaw (Fig. 4B) pre- and poststroke. We then subtracted these maps to assess the relative change in 3-Hz evoked SW power and quantified total power change in the left (excluding the infarct) and right hemispheres (Fig. 4 C and D). Compared to baseline, right paw stimulation after stroke resulted in a significantly weaker increase in power at the stimulation frequency (3 Hz) in both hemispheres (Fig. 4 B and D, **P = 0.0094 in the left hemisphere, **P = 0.0075 in the right). Left paw stimulation provided more informative responses: 3-Hz evoked SW power was significantly reduced across the left (ipsilesional) hemisphere (Fig. 4 A and C **P = 0.0091). Notably, left paw stimulation generated 3-Hz evoked SW power within the right (contralesional) S1_FP at similar levels to baseline, but global power in the rest of the right hemisphere was reduced. On average, right hemisphere power was unchanged (Fig. 4D, ^NSP = 0.16).

Fig. 4.

Unilateral ablation of S1 disrupts global resonance of evoked 3-Hz SWs in both hemispheres. (A and B) Cortical maps of 3-Hz evoked SW power, computed from group-averaged blocks, pre- and poststroke (and ∆ between them) for left (A) and right (B) paw stimulation. (C and D) The quantification of 3-Hz power change averaged within left and right hemispheres for left (C) and right (D) paw stimulation pre- versus poststroke. A two-way ANOVA with Sidak’s multiple comparison test was used to confirm significant differences within individual mice. The quantification in C and D excludes the infarct to capture changes in surviving cortex.

In addition, we tested whether unilateral stroke in S1_FP would disrupt SWs evoked via other somatosensory modalities (SI Appendix, Fig. S2). Using a new cohort of mice, we recorded evoked responses to 3-Hz whisker stimulation before and after photothrombosis of neighboring forepaw S1 (n = 11 mice, 165 blocks total per timepoint). We found that whisker-evoked SWs were spared by injury in S1_FP. While there were focal areas of reduced evoked 3-Hz power in the left motor cortex (M1_W), global 3-Hz whisker evoked SW power was unchanged in both hemispheres overall (SI Appendix, Fig. S2 A and B, ^NSP = 0.70 in the left hemisphere, ^NSP = 0.16 in the right). Consistent with this finding, the dominant spatial mode evoked by 3-Hz whisker stimulation exhibits similar spatial topography and global coherence before and after photothrombosis in S1_FP (SI Appendix, Fig. S2 C and D, ^NSP = 0.34).

Lastly, we examined spontaneous SW power (1 to 1.5 Hz) during stim-OFF epochs, extracted from whole runs without block averaging. We found that spontaneous SWs exhibit less power perilesionally (i.e., in the ring of surviving cortex immediately contiguous to the lesion) but not globally (SI Appendix, Fig. S3). During right paw stim-OFF epochs, left hemisphere spontaneous SW power was significantly decreased after stroke (**P = 0.0036), primarily in the perilesional area, while right hemisphere SW power was unchanged overall (^NSP = 0.20). During left paw stim-OFF epochs, left hemisphere spontaneous SW power was significantly reduced (**P = 0.0044), again observed primarily in the perilesional area, while right hemisphere SW power was unchanged overall (^NSP = 0.45). Thus, unilateral ablation of S1 reduces global resonance of somatosensory-evoked SWs initiated in S1 in either hemisphere but not that of spontaneous SWs.

Mouse Model.

All procedures described below were approved by the Washington University Animal Studies Committee in compliance with the American Association for Accreditation of Laboratory Animal Care guidelines. Mice were raised in standard cages in a double-barrier mouse facility with a 12 h–12 h light/dark cycle and ad libitum access to food and water. All experiments used 10- to 12-wk-old males hemizygous for Thy1-GCaMP6f (JAX 024276) on a C57BL/6J background. Mice of this strain express GCaMP in pyramidal neurons, permitting widefield cortical imaging of excitatory calcium dynamics, which correlates well with multiunit activity (6869–70). Pups were genotyped by PCR prior to experiments to confirm the presence of the Thy1-GCaMP6f transgene using the forward primer 5′-CATCAGTGCAGCAGAGCTTC-3′ and reverse primer 5′-CAGCGTATCCACATAGCGTA-3′.

Cranial Windowing.

Mice received buprenorphine analgesia (0.03 mg/kg, subcutaneous) and were anesthetized with isoflurane (3% induction, 1.5% maintenance). Body temperature was maintained via a thermostatic heating pad. Mice were secured in a stereotactic frame. The scalp was shaved, sterilized with isopropyl alcohol and betadine scrub, locally washed with lidocaine, and then incised at midline and retracted. A custom Plexiglas window with pretapped screw holes for head fixation was attached to the skull using dental cement (C&B-Metabond, Parkell Inc.), completely containing the surgical opening.

Photothrombosis.

Photothrombosis was performed as previously described (71). This protocol induces focal, well-circumscribed, and highly reproducible lesions (72). Briefly, mice (n = 12) were injected with Rose Bengal (6.7 mg/kg, intraperitoneal), and a green diode-pumped solid-state laser (532 nm, 23 mW) collimated to a 1-mm spot was stereotactically centered on the Paxinos coordinates of the forepaw sensory cortex (−2.2 mm in X, +0.05 mm in Y versus bregma) for 10 min (73). Cessation of blood flow was confirmed via laser Doppler.

Infarct Volume Histology.

A separate cohort of three animals were killed 3 d after photothrombosis. Mice were deeply anesthetized with FatalPlus (Vortech Pharmaceuticals) and transcardially perfused with heparinized phosphate-buffered saline (PBS). The brains were removed and fixed in 4% paraformaldehyde for 24 h and transferred to 30% sucrose in PBS. After brains were saturated, they were snap frozen on dry ice and coronally sectioned (50 μm) on a sliding microtome. Sections were stored in 0.2 M PBS, 30% sucrose, and 30% ethylene glycol at −20 °C. For each brain, three sections (spaced 300 µm apart, spanning the infarct) were stained with cresyl violet. A blinded experimenter traced infarct margins in ImageJ to measure infarct area. These measures were multiplied by distance between sections to estimate infarct volume.

Optical Imaging System.

Widefield imaging of cortical calcium dynamics was performed as previously described (19). Briefly, sequential illumination was provided by four LEDs: the 454-nm (blue) LED was used for GCaMP excitation, while the 523-nm (green), 595-nm (yellow), and 640-nm (red) LEDs were used for hyperspectral oximetric imaging. The 523-nm LED was also used as an emission reference to remove hemoglobin confound in the fluorescence signal using a previously validated method (74). The acquisition framerate was set to 16 Hz per channel for whisker imaging studies and 16.81 Hz per channel for forepaw stroke experiments. For image acquisition, we used a cooled, frame-transfer Electron Multiplying Charge Coupled Device camera (iXon 897, Andor) with an 85-mm f/1.4 camera lens and 67-Hz framerate (Rokinon). The field of view was ∼1 cm², covering most of the dorsal convexity of the cerebral cortex (78 µm × 78µm pixels).

Optical Imaging Recordings.

Mice were acclimated to head fixation while secured in a comfortable black felt hammock until they resumed normal resting behavior (whisking, grooming, and relaxed posture). Mice were imaged while awake and/or while anesthetized via intraperitoneal injection with a mixture of ketamine (86.9 mg/kg) and xylazine (13.4 mg/kg). Animals were allowed 10 min for anesthetic induction with the plane of anesthesia measured by respiratory rate and loss of responsiveness to toe pinch. After induction, body temperature was maintained by a thermostatic heating pad with feedback via rectal probe (TCAT-2LV; HP-4M, Physitemp). Whisker stimulation in Figs. 1 and 2, S1 was performed in 10-min epochs using computer-controlled, 40-pounds-per-square-inch air puffs (Picospritzer, Parker Hannifin) in a block design (5 s of 2-Hz, 0.1-s puffs; 5 s rest; 60 blocks/10 min total per mouse, n = 10 mice for awake, n = 9 mice for anesthesia). Data from one ketamine whisker stimulation session was lost due to a file saving error. Stroke imaging sessions were conducted at baseline and 1 wk after photothrombosis in a separate set of experiments performed three years prior using 3-Hz stimulation. We note that evoked SW propagation speed and trajectory is similar regardless of the stimulation frequency or rhythmicity. At each imaging session, forepaw stimulation was delivered by transcutaneous electrical stimulation using microvascular clips (Roboz) in a block design (5 s rest; 5 s of 3-Hz, 1.0-mA, 0.3-ms electrical pulses; 10 s rest; 15 blocks/5 min per mouse × n = 12 mice). For SI Appendix, Fig. S2, whisker stimulation after forepaw stroke used the same block design as forepaw data in Figs. 3 and 4 (5 s rest; 5 s of 3-Hz, 0.1-s puffs; 10 s rest; 15 blocks/5 min per mouse × n = 11 mice).

Optical Imaging Signal Processing.

A binary brain mask was manually drawn in MATLAB for each recording session in each mouse. All subsequent analyses were performed on pixels labeled as brain. Image sequences from each mouse (as well as the brain mask for each mouse) were affine transformed to Paxinos atlas space using the positions of bregma and lambda (73) and then spatially and temporally detrended as previously described (19). The GCaMP6 fluorescence signal (%dF/F) was corrected for varying concentrations of absorptive hemoglobin using 523-nm LED reflectance.

Optical Imaging ROIs.

Regions of interest (ROIs) in the whisker barrel cortex (S1_W) and forepaw sensory cortex (S1_FP) were determined based on functional stimulus mapping from an independent cohort of 21 12-wk-old Thy1-GCaMP6f mice using the same stimulation paradigms as described above. Blocks were averaged within and between mice to generate an average peak response map. An evoked response ROI was defined by any pixels with amplitude >75% of the maximum pixel intensity of the peak response map. ROIs generated from this approach were used in subsequent analyses for averaging within functional territories. ROIs generated in this manner were affine registered to Paxinos atlas space and used to assess responses in subsequent analyses. Infarct ROIs were defined as pixels with >75% reduction of homotopic FC poststroke in awake resting-state data from the same animal. ROIs generated from this approach approximate post hoc measurements of infarct volume histology.

Space-Frequency SVD.

Frequency-specific propagation structure and global coherence were characterized via SVD of optical imaging data in the frequency domain. Performing SVD on multitaper spectral estimates provides a computationally efficient strategy to spectrally smooth and decompose time series data (23, 26). Thus, $$ K$$ Slepian tapers of bandwidth $$ W$$ were applied to pixelwise Fourier transforms. The resulting space-frequency matrix was decomposed into $$ K$$ spatially orthogonal complex modes at each frequency $$ f$$, each with coherence magnitude and phase information.

SVD was applied to full runs for spontaneous activity (time-bandwidth product $$ TW=12$$, $$ K=22$$ tapers; or block averages for evoked activity ($$ TW=3$$, $$ K=5$$)). Tapers results were robust to different parameter choices. Spontaneous and evoked were not directly compared. A single dominant spatial mode was observed in each condition; hence, analyses were restricted to the first SVD component. Note that complex singular vectors are unique up to a complex sign (i.e., a unit phase factor). Thus, dominant spatial modes were phase aligned by rotation in the complex plane according to a reference phase value (here, the global circular mean phase of the mode) prior to averaging across runs.

Following SVD, a measure of global coherence $$ C\left(f\right)$$ is provided by the fractional power of (i.e., proportion of total variance in the coherence spectrum that is explained by) the first mode:

Poststroke global coherence calculations excluded the infarct so as to capture changes in surviving cortex. Space-frequency SVD was implemented in MATLAB and incorporated software routines from the freely available Chronux package (75). Changes in global coherence between states were statistically tested within individual mice by paired two-tailed t test, with significance set at P < α = 0.05.

Power Spectral Analysis.

Global power spectrograms (Fig. 2A) were generated by computing the global average %dF/F signal across all brain pixels within each mouse/run, and then performing a STFT using a 5-s window with 0 overlap to capture the spectral signature of each individual 5-s ON or OFF period within each block. Spectrograms were computed within individual mice and then averaged across mice. Power spectra (Fig. 2B) were computed by averaging across all ON or OFF windows from spectrograms in Fig. 2A. Note, spectra were not variance normalized, allowing comparison of changes in absolute rather than relative band-limited power across conditions. Changes in power (at 2 Hz and summed between 1 to 1.5 Hz) between the ON and OFF spectra in Fig. 2B were analyzed by paired t test, with significance set at P < α = 0.05. Block-averaged power maps (Figs. 2C and 4 A and B and SI Appendix, Fig. S2A) were generated by averaging the %dF/F signal across all blocks within each mouse/run and then computing power at each pixel using Welch’s method. Power maps were averaged within each hemisphere for each mouse, and pre- versus poststroke changes were statistically tested by repeated measures one-way ANOVA with Sidak’s multiple comparison test (comparing within individual mice across time points) with significance set at an adjusted P < α = 0.05. Power was expressed in decibels/Hz using the following equation:

Experimental Design and Statistical Analysis.

Histograms and spectra are plotted as group mean ± 95% CI error bars. Individual data points represent individual mice. Sample size and specific time points were selected on the basis of previous studies examining network connectivity changes in mice (34, 41, 71). Prism 8 software was used to perform statistical testing. Tests used for each analysis are underlined above. Data were verified to be normally distributed prior to statistical testing using the D’Agostino & Pearson normality test. Throughout, significance of P values is depicted as *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001.