Limb stimulation modulates blood-brain barrier permeability.

a. The experimental paradigm in control animals (Ctrl, top panel) and animals that underwent a 30 min stimulation (Stim, bottom panel). Gray lines indicate 1 min test stimulations, lightning indicates the timing of BBB imaging (using injection of sodium fluorescein, NaFlu). b. Top: local field potential (LFP) trace from the corresponding region in L2/3 sensorimotor cortex before during (greyed area) and after 30 min limb stimulation. Bottom: 5 s excerpts of the above trace (left-to-right: before during and after stimulation). c. Image of the cortical window over the rat sensorimotor cortex, followed by total hemoglobin (HbT) concentration maps showing the evolution of the hemodynamic response during and after stimulation. Red rectangle marks the responding region magnified in d. d. Images of the responding arteriole (marked with red arrows) in the rat’s cortex before (left) during (middle) and after stimulation (right). e. The diameter of the responding arteriole during 1 min stimulation. Gray area corresponds to time of stimulation. f. NaFlu permeability maps before (left) and after 30 min stimulation (right) showing tracer accumulation around a responding arteriole. g. Permeability index is higher after stimulation compared to baseline (n=13 rats, mean ± SEM, Wilcoxon, p<0.001). h-i. Fluorescence images of cortical sections of stimulated rats injected with the albumin-binding dye Evans blue (h) and Alexa-488-Alb (i). j-k. Total fluorescence of EB (j) and Alexa488-Alb (k) after stimulation in the contralateral hemisphere compared to the ipsilateral hemisphere. (EB n=4 rats, 32 sections, nested t-test, p=0.0296; Alexa488-Alb n=5 rats, 20 sections, nested t-test, p=0.0229; mean ± SEM). l. Albumin concentration in the contralateral hemisphere relative to the ipsilateral in three different timepoints after stimulation, compared to sham stimulation. (0.5 and 4 h post-stimulation n=8, sham n=7, mean ± SEM, Kruskal-Wallis with FDR correction, p=0.0242, q=0.0406). m. Cortical sections of the area of limb representation from both hemispheres of a stimulated rat (left), and higher magnification images (right) n. Fluorescence image of a cortical small vessel in the stimulated hemisphere showing extravascular accumulation of EB and NaFlu.

Stimulation and cortical perfusion of serum albumin induce LTP.

a. Top: LFP trace from the rat L2/3 sensorimotor cortex before (left) during (middle) and after (right) test stimulation (1 min, 6 Hz, 2 mA). Bottom: LFP trace from the same rat before (left) during (middle) and after (right) test stimulation (1 min, 6 Hz, 2 mA), administered following a 30 min stimulation (6 Hz, 2 mA). b. The somatosensory evoked potential (SEP) of test stimulation at baseline (0 min, BL) and after 30 and 60 minutes (blue, red, black, respectively, each averaged over 360 stimuli). c. SEP in response to test stimulation at baseline (blue) and following a 30 min stimulation (red). d. SEP in response to test stimulation at baseline (blue) and three time points following a 30 min stimulation (60 min, red; 90 min, green; 120 min, purple). e. Maximum amplitude of the SEP (absolute values) following a 30 min stimulation compared to baseline (n=6 rats, mean ± SEM, Wilcoxon, p=0.0312). f. Area under curve (AUC) of the SEP following a 30 min stimulation compared to baseline. (n=6 rats, mean ± SEM, Wilcoxon, p=0.0312). g. Top: 1-hour LFP trace from a representative rat. Bottom: 5 s magnifications of the above trace at selected time points (noted by asterisks). Left to right: baseline activity; during test stimulation; following cortical application of 0.1 mM albumin (Alb); during test stimulation post-Alb. h. SEP amplitude during test stimulation at baseline (normal aCSF, blue) and following 0.1 mM Alb (red). i. Maximum amplitude of the SEP during test stimulation post-Alb compared to baseline (n=5 rats, mean ± SEM, Wilcoxon, p=0.0312). j. AUC of the SEP post-Alb compared to baseline. (n=5 rats, mean ± SEM, Wilcoxon, p=0.0312). k. Power spectrum density of 10 min spontaneous activity before (blue) and post-Alb (red).

Stimulation-induced plasticity is associated with BBB modulation.

a. Timeline of the experimental protocol for stimulations and imaging with blockers application (CNQX/AP5 (50 μM)/mβCD (10 μM)/SJN (0.3 mM). b. NaFlu permeability maps of the cortical window before (control, left), and after CNQX/AP5/SJN + 30 min stimulation (right). c. PI following blockers application before and after stimulation compared to a 30 min stimulation. (% change from baseline, 30 min stim n=6, AP5 n=4, CNQX n=5, mβCD n=6, SJN n=5; mean ± SEM, Kruskal-Wallis with FDR correction, *p<0.05). d. SEP amplitude in response to test stimulation. baseline (blue); following application of CNQX (left, red) or AP5 (right, red); following CNQX + 30 min stimulation (left, black); and following AP5 + 30 min stimulation (right, black); e. Left: Max SEP amplitude to test stimulation following CNQX compared to baseline. (mean ± SEM, n=5, paired t-test, p=0.0176). Right: Max SEP amplitude following AP5 compared to baseline. f. SEP amplitude in response to test stimulation. baseline (blue); following SJN (black); following SJN + 30 min stimulation (red); and following mβCD (right). g. Max SEP amplitude following 30 min stimulation or albumin application compared to blockers, and blockers + 30 min stimulation. h. AUC of the SEP for test stimulation following 30 min stimulation or albumin application compared to blockers and blockers + 30 min stimulation. (g-h: % change from baseline, mean ± SEM, 30 min stim n=5, Alb n=4, AP5 n=5, CNQX n=5, SJN n=4, mβCD n=5; Kruskal-Wallis with FDR correction, *p<0.05).

Neuronal activity regulates BBB transport and synaptic plasticity genes.

a. Scatter plot of gene expression from RNA-seq in the contralateral cortex 24 vs. 1 h after 30 min stimulation. The x axis represents the log fold change, and the y axis represents the mean expression levels. Blue dots indicate statistically significant differentially expressed genes (DEGs) by Wald Test (n=8 rats per group). b. Top Gene Ontologies (GO) enriched terms in the contralateral cortices of rats 24 vs. 1 h after stimulation. c. Vascular cell specific DEGs to pericytes (PC), arterial smooth muscle cells (aSMC), and venous smooth muscle cells (vSMC). (PC n=42, aSMC n=70, p=0.0017, Chi-square). d. Scatter plot of DEGs divided by groups of interest: BBB properties, NVU properties, Synaptic plasticity, and inflammatory related genes in the contralateral (red) vs. ipsilateral (black) cortices of stimulated rats. Circles represent genes with no significant differences between 1 and 24 h post-stimulation. Upward and downward triangles indicate significantly up- and down-regulated genes, respectively (synaptic plasticity n=53, ***p<0.001, Chi-square).

Cortical activation in fMRI co-localizes with BBB modulation.

a. Subjects were given an elastic stress-ball to squeeze continuously for the length of the session (30 min). b. Timeline of the experimental protocol for task performance and MRI sequences. c. Activation map for the localizer task displayed over the inflated brain of an exemplary subject. (p<0.05, Family wise error (FWE) corrected, neurologic convention). d. Activation maps for the localizer task displayed over anatomical axial slices of an exemplary subject. (p<0.05, FWE corrected, radiologic convention). e. Left: Superimposed masks of BBB modulated voxels (red) and fMRI activation (yellow). Right: Co-localized voxels map on the same slice (cyan). f. BBB modulation in the primary motor cortex (M1, precentral gyrus, PrG) and primary somatosensory cortex (S1, postcentral gyrus, PoG) for subjects performing the task compared to controls (% of voxels in each region with BBB leakage). (Task n=6, controls n=10, i–ipsi, c-contra, mean ± SEM, *p<0.05, two-way ANOVA with FDR correction). g. Heatmap of BBB modulated voxels percentage in motor/sensory related areas and some non-activated cortical regions of task vs. controls (+Task n=6, -Task n=10, i–ipsi, c-contra, mean ± SEM, *p<0.05, **p<0.001, ***p<0.0001, ns – non significant, two-way ANOVA with FDR correction).