Abstract
Summary
Brain microvessels possess the unique properties of a blood-brain barrier (BBB), tightly regulating the passage of molecules from the blood to the brain neuropil and vice versa. In models of brain injury, BBB dysfunction and the associated leakage of serum albumin to the neuropil have been shown to induce pathological plasticity, neuronal hyper-excitability, and seizures. The effect of neuronal activity on BBB function and whether it plays a role in plasticity in the healthy brain remain unclear. Here we show that neuronal activity induces modulation of microvascular permeability in the healthy brain and that it has a role in local network reorganization. Combining simultaneous electrophysiological recording and vascular imaging with transcriptomic analysis in rats, and functional and BBB-mapping MRI in human subjects we show that prolonged stimulation of the limb induces a focal increase in BBB permeability in the corresponding somatosensory cortex that is associated with long-term synaptic plasticity. We further show that the increased microvascular permeability depends on neuronal activity and involves caveolae-mediated transcytosis and transforming growth factor beta signaling. Our results reveal a role of BBB modulation in cortical plasticity in the healthy brain, highlighting the importance of neurovascular interactions for sensory experience and learning.
Introduction
The intricate components comprising the blood-brain barrier (BBB) regulate the bi-directional transport of molecules between the brain and the circulatory system and maintain the extracellular environment needed for normal function of the neurovascular unit (NVU) (Chow and Gu, 2015). Dysfunction of the BBB is common in many neurological disorders and is associated with impaired cross-BBB influx and efflux (Obermeier et al., 2013). For example, seizures were shown to increase cross-BBB influx of serum albumin – in a process partially mediated by excessive release of neuronal glutamate and subsequent activation of N-methyl-D-aspartate (NMDA) receptors (Vazana et al., 2016).
Emerging evidence suggests that cross-BBB influx/efflux may also change in response to physiological neuronal activity in the healthy brain: (1) whisker stimulation (1 hour, 2 Hz) was shown to increase BBB influx of insulin-like growth factor-I in rats (Nishijima et al., 2010); (2) physiological neuronal activity was shown to regulate BBB efflux transport in mice (Pulido et al., 2020); and (3) circadian rhythms were shown to alter BBB efflux in Drosophila (Zhang et al., 2018) and mice (Pulido et al., 2020). However, the mechanisms and the role of these physiological changes in BBB properties are not known.
Here, we sought to test whether physiological BBB modulation may involve mechanisms previously linked to pathological BBB changes, e.g., transcytosis of serum albumin (Ivens et al., 2007; Seiffert et al., 2004) and transforming growth factor β (TGF-β) signaling (David et al., 2009; Heinemann et al., 2012; Ivens et al., 2007; Weissberg et al., 2015). As TGF-β signaling has been linked to synaptic plasticity (Dahlmanns et al., 2022; Gradari et al., 2021; Patel and Weaver, 2021; Weissberg et al., 2015), we also explored the hypothesis that physiological BBB modulation may be associated with neuroplasticity. Lastly, we sought to provide the first evidence for physiological BBB modulation in humans.
Results
Stimulation increases BBB permeability
To test for changes in BBB permeability in response to stimulation we combined electrophysiological recordings with BBB imaging using intravital microscopy (Figure 1, see timeline of stimuli in Figure 1a). First, we localized vascular response to sensory stimulation using wide-field imaging of arterial diameter and monitoring change in total hemoglobin (HbT) signal during a 1 min stimulation of the limb (Figure 1b-d) (Bouchard et al., 2009). Fluorescent angiography before and immediately after 30 min of limb stimulation showed local extravasation of the tracer sodium fluorescein (NaFlu) around the responding blood vessels, indicative of BBB leakage (Figure 1f-g). Histological analysis confirmed the presence of NaFlu in the extravascular space around small vessels in the somatosensory cortex, contralateral (and not ipsilateral) to the stimulated limb (Figure 1m-n). Next, we intravenously injected a separate cohort of rats with either the albumin-binding dye Evans blue (EB) or albumin conjugated with Alexa Fluor 488 (Alexa488-Alb) without performing craniotomy. Brain extravasation of both dyes was found in the contralateral hemisphere, indicating the presence of serum albumin in the neuropil (Figure 1h-k). ELISA performed in cortical tissue at different time points after stimulation confirmed a higher concentration of albumin in the stimulated sensorimotor cortex early after stimulation (30 min) compared to non-stimulated sham-controls, that declined 4- and 24-hours post-stimulation (Figure 1l). As a positive control, brain tissues were also prepared from rats subjected to photothrombosis (PT)-induced ischemia (see methods). Albumin concentrations in the peri-ischemic cortex 24 hours following PT were significantly higher compared to sham and all time-points after stimulation (Figure S1d). These results suggest that the focal modulation of BBB permeability following stimulation is distinct from the extensive BBB dysfunction observed in injury models.
BBB permeability is associated with synaptic potentiation
To study changes in brain activity during repetitive stimulation of the limb, we measured somatosensory evoked potentials (SEP) by averaging neuronal response to a 1 min test stimulation before and after the 30 min stimulation train (6 Hz, 2 mA). SEP amplitudes and area under the curve were increased in stimulated, but not in control cortices (Figure 2a-e). Potentiation lasted until the end of each experiment (Figure 2f) and was therefore indicative of long-term potentiation (LTP).
To test whether stimulation-induced albumin extravasation has a role in LTP, we recorded local field potentials (LFP) before and after exposure of the cortex to 0.1 mM albumin (corresponding to 25% of serum concentration) (Ivens et al., 2007). Albumin application increased the amplitude and frequency of spontaneous activity (Figure 2g,k). Neuronal response to the stimulation showed synaptic potentiation in the presence of albumin compared to that recorded under aCSF (Figure 2h-j), indicating that cortical exposure to serum albumin results in synaptic potentiation.
BBB modulation is activity-dependent
To test the neurotransmission pathways involved in regulation of BBB permeability during stimulation we exposed the cortical window to selective antagonists of ionotropic glutamate receptors (Figure 3a). The AMPA/kainate receptor blocker CNQX prevented the stimulation-induced BBB modulation, whereas the NMDAR blocker AP5 did not (Figure 3b-c). CNQX also reversibly blocked synaptic activity and prevented the vascular response to neuronal activation, indicated by the absence of vasodilation (AKA “neurovascular coupling”; Figure 3d-e, Figure S2a, S3a,c). These results indicate that stimulation-induced BBB modulation was mediated by neuronal activation. In contrast, AP5 did not affect baseline SEP (Figure 3d-e), nor the vascular response (Figure S2a, S3b). Synaptic potentiation was blocked under both CNQX or AP5 (Figure 3d-e, g-h, Figure S2a), consistent with the crucial role of AMPA and NMDA receptors in synaptic plasticity. Next, we tested for molecular correlates of synaptic strength by comparing postsynaptic density 95 (PSD-95) expression levels in the cortex at different time points. PSD-95 is known to increase following LTP and AMPAR trafficking (El-Husseini et al., 2000; Holtmaat and Svoboda, 2009). Western blot analysis showed significantly higher levels of PSD-95 expression in stimulated rats as compared to sham, at both 30 min and 24 h after stimulation (Figure S2c-d).
Albumin extravasation and LTP require caveolae-mediated transcytosis
In endothelial cells, albumin is transported mainly via caveolae-mediated transcytosis (CMT) through binding to the gp60 receptor (Tiruppathi et al., 1997). To test whether CMT mediates stimulation-induced transcytosis of albumin (John et al., 2003) we perfused the exposed cortex with methyl-β-cyclodextrin (mβCD), which inhibits albumin uptake by disrupting caveolae (Schnitzer and Oh, 1994; Skotland et al., 2020). Cortical application of mβCD prevented stimulus-induced BBB opening (Figure 3c), without affecting vascular response to stimulation (Figure S3d). mβCD also prevented stimulus-induced LTP (Figure 3f-h), further supporting the role of CMT-mediated albumin extravasation in activity-dependent plasticity in vivo.
BBB modulation and LTP involve TGF-β signaling
Under pathological conditions, synaptogenesis and pathological plasticity that follow BBB dysfunction are mediated by the activation of TGF-βR1 signaling in astrocytes (Weissberg et al., 2015). TGF-β1 was also shown to directly induce BBB opening in endothelial cell cultures (McMillin et al., 2015), and to regulate hippocampal synaptic plasticity (Dahlmanns et al., 2022; Gradari et al., 2021). We therefore tested the effect of stimulation in the presence of the specific TGF-βR1 blocker SJN2511 (SJN, 0.3 mM (Gellibert et al., 2004; Weissberg et al., 2015)). Cortical exposure to SJN prevented stimulation-induced BBB opening (Figure 3b-c) with no effect on neurovascular coupling (Figure S3e). SJN also prevented stimulation-induced LTP (Figure 3f-h). These results demonstrate that TGF-β signaling is required for activity-dependent BBB modulation, and the consequent activity-dependent plasticity, but not for neurovascular coupling.
Regulation of BBB transport and plasticity genes
To explore gene-expression changes related to our stimulation protocol, we conducted bulk RNA-sequencing transcriptome analysis of tissue dissected from the sensorimotor area of the cortex contralateral and ipsilateral to stimulated limb. Tissue harvested 1 and 24 h post-stimulation (30 min) showed differentially expressed genes (DEGs) in the contralateral (“stimulated”) compared to the ipsilateral (“non-stimulated”) cortex (13.2% and 7.3% 24 and 1 h, respectively, Figure 4a, S4a). To assess the variability between paired samples of stimulated and non-stimulated cortex of each animal at two time points following stimulation (24 h vs 1 h) we used the Jensen–Shannon divergence metric (JSD) (Sudmant et al., 2015). JSD calculations using normalized counts indicated that the RNA expression 24 h after stimulation was significantly different from that of rats 1 h after stimulation (Figure S4b). Gene Ontology (GO) pathway enrichment analysis of DEGs from stimulated rats indicated DEGs were primarily involved in synaptic plasticity processes (Figure 4b). A significantly higher number of DEGs related to synaptic plasticity was found in the contralateral cortex of stimulated rats compared to ipsilateral (p<0.001, Figure 4d). No significant differences were found in BBB or inflammation related genes between the hemispheres, suggesting that the stimulation protocol is not associated with the transcriptional change observed in BBB dysfunction under pathological conditions (Cacheaux et al., 2009; Kim et al., 2017) (Figure 4d). In contrast, a larger amount of transporter genes such as solute carrier transporters (Slc) and ATP-binding cassette (ABCs) families were differentially expressed in the stimulated compared to the non-stimulated cortex (Figure S4c-d). Both DEG and enrichment analysis were consistent with activation of TGF-β signaling pathway (Figure 4d). Analysis of vascular cell-specific DEGs (Vanlandewijck et al., 2018) showed significantly more arterial smooth muscle cells specific genes in the stimulated cortex (Figure 4c, Figure S4e). These results support the roles of stimulation-induced modulation of BBB transcellular transport and TGF-β signaling in activity-dependent plasticity (Figure S5).
Physiological BBB modulation in the human brain
Finally, we tested for neuronal activity-induced BBB modulation in healthy human volunteers by measuring the blood-oxygenation-level dependent (BOLD) response to a 30 min stress-ball squeeze task and blood-to-brain transport using dynamic contrast-enhanced (DCE) MRI (Figure 5a-b). As expected, BOLD response was found in the pre- and post-central gyrus, corresponding with the primary motor and sensory cortices, respectively, of the hemisphere contralateral to the squeezing hand (Figure 5c-d). DCE-MRI showed higher BBB permeability to the gadolinium-based contrast agent gadoterate meglumine (Dotarem) in sensory/motor-related areas of subjects performing the task, as compared to controls (Figure 5e-g). Non activated cortical regions such as the medial frontal cortex and middle occipital gyrus did not differ in BBB permeability between task and controls (Figure 5g). Comparison of the spatial distribution of functional and permeability changes revealed a higher percentage of activated regions with modulated BBB in task performers than in control subjects (Figure 5e-g).
Discussion
This study demonstrates that prolonged physiological stimulation modulates BBB permeability in the rat and human cortex. In rats, we show that prolonged paw stimulation results in increased cross-BBB extravasation of both low- and high-molecular weight tracers. We also show that this extravasation is: (1) focal and localized to the cortical area corresponding to the stimulation; (2) mediated by AMPA neurotransmission, TGF-β signaling, and caveolae transcytosis; and (3) associated with synaptic potentiation and increased postsynaptic density. In humans, we present the first evidence for increased BBB permeability in response to prolonged limb activity. Together, our results suggest that stimulation-induced BBB modulation may play a role in synaptic potentiation.
Physiological changes in cross-BBB transport have been previously reported in response to whisker stimulation, circadian rhythms, and age (Nishijima et al., 2010; Pulido et al., 2020; Yang et al., 2020). Here, we examined whether the mechanisms involved in this process may include those identified in our studies of pathological BBB modulation. Previous studies from our group showed that brain injury is associated with albumin transcytosis, albumin binding to TGF-βR in astrocytes, and pathological TGF-β-mediated plasticity (Salar et al., 2016). Here, we show that prolonged physiological sensory-motor activation is also associated with albumin extravasation (although at a significantly lower concentration than under pathological conditions, Figure 1l and Figure S1d) and TGF-βR signaling. We further show that blocking caveolae-mediated albumin transcytosis or TGF-βR signaling inhibits stimulation-evoked increase in BBB permeability and synaptic potentiation (without affecting the physiological hemodynamic response). These results suggest that caveolae albumin transcytosis and TGF-β signaling are key mediators of activity-dependent BBB modulation and subsequent synaptic plasticity. These findings are supported by reports linking TGF-β signaling to: synaptic plasticity induced by environmental enrichment (Dahlmanns et al., 2022); enhanced spine and synapse formation (Patel and Weaver, 2021); and hippocampal neuroplasticity (Gradari et al., 2021). Moreover, our finding of caveolae-mediated albumin transcytosis in pial arterioles is in line with previous studies showing: caveolae-dependent transcytosis of albumin (Tecedor et al., 2013); caveolae abundance in arteriolar ECs and not capillary ECs (Chow et al., 2020); and the contribution of caveolae transcytosis to neurovascular coupling (Andreone et al., 2017; Chow et al., 2020). Notably, in our experiments, caveolae inhibition prevented BBB modulation yet did not significantly affect neurovascular coupling, likely due to the contribution of other vasodilatory pathways such as endothelial nitric oxide synthase.
As mentioned above, albumin is a known activator of TGF-β signaling, and TGF-β has a well-established role in neuroplasticity. Interestingly, emerging evidence suggests that TGF-β also increases cross-BBB transcytosis (Betterton et al., 2022; Kaplan et al., 2020; McMillin et al., 2015; Schumacher et al., 2023) Hence, we propose the following two-part hypothesis for the TGF-β/BBB-mediated synaptic potentiation observed in our experiments: (1) prolonged stimulation triggers TGF-β signaling and increased caveolae-mediated transcytosis of albumin; and (2) extravasated albumin induces further TGF-β signaling, leading to synaptogenesis and additional cross-BBB transport – in a self-reinforcing positive feedback loop. Future research is needed to examine the validity of this hypothesis.
Our study also shows that physiological BBB modulation requires AMPAR glutamate signaling and not NMDA. This finding is consistent with studies showing the dominant role of AMPAR in electrophysiological and hemodynamic responses to somatosensory stimulation (Graves et al., 2021; Gsell et al., 2006; Zhang et al., 2015). Moreover, the well-established role of AMPAR in neuroplasticity, long-term potentiation (LTP), and experience-dependent sensory cortical plasticity (Daw et al., 1993; Feldman, 2009; Watt et al., 2004) further supports the plasticity-related role of activity-dependent BBB modulation.
Our animal experiments show that a 30 min limb stimulation (at 6Hz and 2mA) increases cross-BBB influx, while a 1 min stimulation (of similar frequency and magnitude) does not. We believe that both types of stimulations fall within the physiological range because our continuous electrophysiological recordings showed no signs of epileptiform or otherwise pathological activity. Moreover, the recorded SEP levels were similar to those reported in previous physiological LTP studies in rats (Eckert and Abraham, 2010; Han et al., 2015; Mégevand et al., 2009) and humans (McGregor et al., 2016). Future studies are needed to explore the BBB modulating effects of additional stimulation protocols – with varying durations, frequencies, and magnitudes. Such studies may also elucidate the temporal and ultrastructural characteristics that differentiate between physiological and pathological BBB modulation.
In humans, skill acquisition often involves motor training sessions that last ≥30 minutes (Bengtsson et al., 2005; Classen et al., 1998) and result in physiological plasticity of sensory and motor systems (Classen et al., 1998; Draganski et al., 2004; Sagi et al., 2012). We, therefore, believe that the experimental task of our human study (30 minutes of repetitive squeezing of an elastic stress-ball) represents physiological activity, with neuronal activation in primarily motor and sensory areas (Halder et al., 2005). We further propose that the activity-related increase in BBB permeability we observed may be associated with learning.
A key limitation of our animal experiments is the fact they were performed under anesthesia, due to the complex nature of the experimental setup (i.e., simultaneous cortical imaging and electrophysiological recordings). Anesthetic agents can potentially alter neuronal activity, SEPs, CBF, and vascular responses (Aksenov et al., 2015; Lindauer et al., 1993; Masamoto and Kanno, 2012). To minimize these effects, we used ketamine-xylazine anesthesia, which unlike other anesthetics, was shown to maintain robust BOLD and SEP responses to neuronal activation (Franceschini et al., 2010; Shim et al., 2018). Another limitation of our animal study is the potentially non-specific effect of mβCD – an agent that disrupts caveola transport but may also lead to cholesterol depletion (Keller and Simons, 1998). To mitigate this issue, we used a very low mβCD concentration (10μM), orders of magnitude below the concentration reported to deplete cholesterol (Koudinov and Koudinova, 2001). Lastly, our animal study is limited by the inclusion of solely male rats. While our findings in humans did not point to sex-related differences in stimulation-evoked BBB modulation, larger animals and human studies are needed to examine this question.
To conclude, our study suggests that BBB modulation in response to physiological neuronal activity involves mechanisms previously identified in pathology-induced BBB changes (despite the significantly smaller increase in cross-BBB influx under physiological conditions). Our findings further suggest that physiological modulation of cross-BBB influx may play a role in synaptic plasticity, calling for future research into the dynamics of BBB function in health and disease.
Methods
Animal preparation and surgical procedures
All experimental procedures were approved by the Animal Care and Use Ethical Committee at the Ben-Gurion University of the Negev, Beer-Sheva, Israel.
Adult Sprague Dawley male rats (300-350 g; Harlan Laboratories) were kept under a 12:12 h light and dark regimen and supplied with drinking water and food ad libitum. Surgical procedures were performed as reported previously(Prager et al., 2010). Rats were deeply anesthetized by intraperitoneal administration of ketamine (100 mg/ml, 0.08 ml/100 g) and xylazine (20 mg/ml, 0.06 ml/100 g). The tail vein was cannulated, and animals were placed in a stereotactic frame under a SteREO Lumar V12 fluorescence microscope (Zeiss). Body temperature was continuously monitored and kept stable at 37 ± 0.5°C using a feedback-controlled heating pad (Physitemp). Heart rate, breath rate, and oxygen saturation levels were continuously monitored using MouseOx (STARR Life Sciences). A cranial section (4 mm caudal, 2 mm frontal, 5 mm lateral to bregma) was removed over the right sensorimotor cortex. The dura and arachnoid layers were removed, and the exposed cortex was continuously perfused with artificial CSF (aCSF (Prager et al., 2010), containing (in mM): 124 NaCl, 26 NaHCO3, 1.25 NaH2PO4, 2 MgSO4, 2 CaCl2, 3 KCl, and 10 glucose, pH 7.4).
Stimulation protocol
The left forelimb or hind limb of the rat was stimulated using Isolated Stimulator device (AD Instruments) attached with two subdermal needle electrodes (0.1 ms square pulses, 2-3 mA) at 6 Hz frequency. Test stimulation consisted of 360 pulses (60 s) and delivered before (as baseline) and after long-duration stimulation (30 min, referred throughout the text as ‘stimulation’). In control and albumin rats, only short-duration stimulations were performed. Under sham stimulation, electrodes were placed without delivering current.
Electrophysiology
In-vivo recordings
Stimulus-evoked potentials were recorded in the somatosensory cortex. Local field- and somatosensory evoked potentials (LFP, SEP) were recorded using a glass microelectrode (1.5 × 1.1 mm, 1-2 µm tip, 2-10 MΩ) filled with aCSF, which was inserted with an isolated, chlorinated silver wire connected to a headstage and amplifier (EXT-02B, npi electronic, Germany). The electrode was attached to a digital micromanipulator (MP-225, Sutter Instruments, CA, USA) and was carefully placed over the identified area of the sensorimotor response to the stimulation.
Analysis of electrophysiological recordings
Data acquisition was performed using PowerLab and LabChart (AD Instruments), for stimulation, and extracellular recordings. Signals were digitized (1 kHz) and filtered (high-pass 1 Hz, low pass 45 or 300 Hz). Analyses were performed using in house MATLAB scripts. Two measures were used to quantify synaptic plasticity: the maximal amplitude of the dominant negative peak and the absolute area under curve (AUC) of the SEP waveform (i.e., all positive and negative peaks between 10–160 ms post-stimulus period). For power analysis in albumin experiments, 10 min segments recorded before starting and after terminating albumin perfusion were used.
Drug Application
In some experiments the following blockers were added to the aCSF and continuously perfused to the cortical window for 50 min: D-(-)-2-Amino-5-phosphonopentanoic acid (AP5, 50 μM), 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX, 50 μM), methyl-β-cyclodextrin (MβCD, 10 μM), or SJN2511 (SJN, 0.3 mM (Weissberg et al., 2015)). To simulate BBB opening, bovine serum albumin (Alb, 0.1 mM (Seiffert et al., 2004)) was added to the aCSF and continuously perfused to the cortical window for 30 min.
In-vivo imaging
Dynamic imaging of regional cerebral blood flow, optical imaging of cortical oxygenation, neurovascular coupling (NVC), and vessel diameter measurements were carried out as previously described (Bouchard et al., 2009; Levi et al., 2012; Prager et al., 2010). Before, during, and after electrical limb stimulation, full-resolution (2560 × 2160 px) images were obtained for regional localization of the hemodynamic response, followed by acquisition of localized (512 × 512 px) images of cortical surface vessels (2 frames/s, EMCCD camera, Neo DC-152Q-COO-FI; Andor Technology) under 525/50 band pass filter (38 HE, Zeiss). Images were analyzed offline using in-house MATLAB (MathWorks) scripts.
Fluorescent angiography and BBB permeability measurements and analyses were performed as described previously (Prager et al., 2010; Vazana et al., 2016). The non-BBB permeable fluorescent dye sodium fluorescein (NaFlu) was injected intravenously (0.2 ml, 1 mg/ml). Images of cortical surface vessels were obtained before, during, and after tracer injection at 5 frames per second. Across all experiments, acquired images were the same size (512 × 512 pixel, ∼1×1 mm), centered above the responding arteriole. Images were analyzed offline using MATLAB as described (Vazana et al., 2016). Briefly, image registration and segmentation were performed to produce a binary image, separating blood vessels from extravascular regions. For each extravascular pixel, a time curve of signal intensity over time was constructed. To determine whether an extravascular pixel had tracer accumulation over time (due to BBB permeability), the pixel’s intensity curve was divided by that of the responding artery (i.e., the arterial input function, AIF, representing tracer input). This ratio was termed the BBB permeability index (PI), and extravascular pixels with PI > 1 were identified as pixels with tracer accumulation due to BBB permeability.
Histology and immunohistochemistry
Albumin extravasation was confirmed histologically in separate cohorts of rats that were anesthetized and stimulated without craniotomy surgery. Assessment of albumin extravasation was performed using a well-established approach that involves peripheral injection of either labeled-albumin (bovine serum albumin conjugated to Alexa Flour 488, Alexa488-Alb) or albumin-labeling dye (Evans blue, EB – a dye that binds to endogenous albumin and forms a fluorescent complex), followed by histological analysis of brain tissue (Ivens et al., 2007; Obermeier et al., 2013 ; Lapilover et al., 2012; Veksler et al., 2020, Ahishali et al., 2020). Since extravasated albumin is taken up by astrocytes (Ivens et al., 2007; Obermeier et al., 2013 and others), it can be visualized in the brain neuropil after brain removal and fixation (Ivens et al., 2007; Lapilover et al., 2012; Veksler et al., 2020, Ahishali et al., 2020). Five rats were injected with Alexa488-Alb (1.7 mg/ml) and five with EB (2%, 20 mg/ml, n=5). The injections were administered via the tail vein. Following injection, rats were transcardially perfused with cold phosphate-buffered saline (PBS), followed by PBS containing 4% paraformaldehyde (PFA). Brains were then removed, fixed overnight (4% PFA, 4°C), cryoprotected with sucrose gradient (10% followed by 20 and 30% sucrose in PBS) and frozen in optimal cutting temperature (OCT) compound. Coronal sections (40-μm thick) were obtained using a freezing microtome (Leica Biosystems) and imaged for dye extravasation using a fluorescence microscope (Axioskop 2; Zeiss) equipped with a CCD digital camera (AxioCam MRc 5; Zeiss).
For immunostaining, brain slices were incubated in blocking solution (5% donkey serum in 0.1% Triton X-100/tris-buffered saline (TBS)), then incubated in primary antibody, followed by secondary antibody. Finally, slices were mounted and stained for DAPI (DAPI-Fluoromount-G, Invitrogen) to label nuclei. Staining was performed against glial fibrillary acidic protein (mouse anti GFAP−Cy3™, Sigma-Aldrich, C9205), Secondary antibody was anti-mouse Alexa Fluor 488 (Thermo-Fisher, A21206).
Immunoassays
Enzyme-linked immunosorbent assay (ELISA) for albumin extravasation
To assess the degree of albumin extravasation following stimulation, forty adult Sprague Dawley male rats (200-300g) were randomly divided into five groups. Rats were anesthetized (ketamine (75 mg/kg) and xylazine (5 mg/kg)) and underwent either a 30 min stimulation, sham stimulation (electrodes placed without current delivery), or photothrombosis stroke (PT) as previously described (Lippmann et al., 2017; Schoknecht et al., 2014). Briefly, Rose Bengal was administered intravenously (20 mg/kg) and a halogen light beam was directed for 15 min onto the intact exposed skull over the right somatosensory cortex. Then, rats were deeply anaesthetized and transcardially perfused with cold PBS, and brains were extracted at three time points following stimulation (immediately (real and sham), 4 and 24 h post stimulation) and 24 h post stroke. The ipsi- and contra lateral somatosensory cortices were dissected, weighed, and snap-frozen on dry ice. Fresh frozen samples were then homogenized in PBS + 1% Triton X-100 (v/v). Tissue homogenates were centrifuged at 13000 g for 10 min and supernatants were stored at -80°C until further analysis. Albumin levels were determined in diluted samples, using the Rat Albumin ELISA kit (Bethyl Laboratories, TX, USA) according to the manufacturer’s instructions.
Western Blotting
Brain tissues were collected from stimulated and control rats as described above. Protein lysates were separated using 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS–PAGE, Bio-Rad) and then transferred to nitrocellulose membranes (Bio-Rad). Membranes were washed with TBS + 0.1% Tween20 (TBST, 3 × 5 min), then blocked with 5% BSA in TBS for 1 h. Membranes were then incubated overnight at 4°C with primary antibodies, followed by 2 h incubation with secondary antibodies. Protein bands were visualized with the enhanced chemiluminescent HRP substrate Crescendo Immobilon (Millipore). Analysis was performed in ImageJ, and band density levels were expressed as fold-changes to controls after normalization to β-actin and controls within the same blot. Primary antibodies used were mouse anti-PSD-95 (1:1500, Thermo Fisher Scientific, MA1-045), anti-β-actin (Loading control, 1:2000, Abcam, ab8226). Secondary antibodies were Goat anti-mouse (1:4000, Abcam, ab205719) and Goat anti-rabbit (1:4000, Abcam, ab97051).
RNA Extraction & Purification
Total RNA was extracted from 100-150 mg of flash-frozen rat somatosensory cortex tissue using TRIzol reagent (ThermoFisher Scientific #15596026) as previously described (Rio et al., 2010). Glycoblue (Invitrogen #AM9515) was added to assist with RNA pellet isolation. Pellets were washed with 75% ethanol and allowed to dry before being dissolved in 20 μL of RNase-DNase free water. 2 μL of 10x DNase I Buffer (New England Biolabs, #m0303s), 1 μL of rDNase I (ThermoFisher Scientific, AM2235), and DNase Inactivation Reagent (ThermoFisher Scientific, AM1907) were added to each RNA sample for purification.
RNA Sequencing & Bioinformatics
Library preparation was performed by the QB3-Berkeley Functional Genomics Laboratory at UC Berkeley. mRNA selection was done with the oligo-dT beads. Both total RNA and oligo dT bead selected mRNA quality was assessed on an Agilent 2100 Bioanalyzer. Libraries were prepared using the KAPA RNA Hyper Prep kit (Roche KK8540). Truncated universal stub adapters were ligated to cDNA fragments, which were then extended via PCR using unique dual indexing primers into full length Illumina adapters. Library quality was checked on an AATI (now Agilent) Fragment Analyzer. Libraries were then transferred to the QB3-Berkeley Vincent J. Coates Genomics Sequencing Laboratory, also at UC Berkeley. Library molarity was measured via quantitative PCR with the KAPA Library Quantification Kit (Roche KK4824) on a BioRad CFX Connect thermal cycler. Libraries were then pooled by molarity and sequenced on an Illumina NovaSeq 6000 S4 flowcell for 2 × 150 cycles, targeting at least 25M reads per sample. Fastq files were generated and demultiplexed using Illumina bcl2fastq2 v2.20. No other processing, trimming, or filtering was conducted. RNA-seq reads were aligned using STAR v2.9.1a against UCSC rat genome rn6. Downstream analysis was performed using FastQC v0.11.9, MultiQC v1.11, and featureCounts. R studio was used for data analysis, R v4.0.2. with the packages tidyverse v1.3.1, IHW, DESeq2. Gene ontology (GO) analysis of enriched genes in ipsi- and contra-lateral somatosensory cortices of stimulated rats was performed using Metascape (Zhou et al., 2019) (http://metascape.org/). Analysis of cell-specific and vascular zonation genes was performed as described (Vanlandewijck et al., 2018), using the database provided in (http://betsholtzlab.org/VascularSingleCells/database.html).
Magnetic Resonance Imaging
Participants
Male and female healthy individuals, aged 18-35, with no known neurological or psychiatric disorders were recruited to undergo MRI scanning while performing a motor task (n=6; 3 males and 3 females). MRI scans of 10 sex- and age-matched individuals (with no known neurological or psychiatric disorders) who did not perform the task were used as control data (n=10; 5 males and 5 females).
All procedures were approved by the Soroka University Medical Center Institutional Review Board. All participants gave their written informed consent before participation.
Data acquisition
MRI scans were performed in the imaging center at the Soroka Medical Center using a 3T Philips Ingenia MRI scanner. A high-resolution T1-weighted anatomical scan (3D gradient echo, TE/TR 3.7/8.2 ms, voxel size of 1 mm3, 256 × 256 acquisition matrix), followed by a T2-weighted scan (TE/TR = 90/3000 ms, voxel size 0.45 × 0.45 × 4 mm). Functional data were collected by using a gradient echo EPI, with voxel size of 3 × 3 × 3 mm, TE/TR = 35/2000 ms, 83 × 80 acquisition matrix, flip angle of 90°, and FOV 109 mm. For calculation of pre-contrast longitudinal relaxation times, variable flip angle method was used (3D T1w-FFE, TE/TR = 2/10 ms, acquisition matrix: 256 × 256, voxel size: 0.89 × 0.89 × 6 mm, flip angles: 5, 15, 20 and 25°). Dynamic contrast-enhanced (DCE) sequence was then acquired (Axial, 3D T1w-FFE, TE/TR = 2/4 ms, acquisition matrix: 192 × 187 (reconstructed to 256 × 256), voxel size: 0.9 × 0.9 × 6 mm, flip angle: 20°, ýt = 10 s, temporal repetitions: 100, total scan length: 16.7 min). An intravenous bolus injection (0.1 mmol/kg, 0.5 M) of the gadolinium-based contrast agent gadoterate meglumine (Dotarem, Guerbet, France) was administered using an automatic injector after 5 dynamic scans at a rate of 1.5 ml/s.
Preprocessing of functional data
fMRI data preprocessing was performed using the Statistical Parametric Mapping package (SPM12; Wellcome Trust Centre for Neuroimaging, London, UK) for MATLAB. The first ten functional images of each run series were discarded to allow stabilization of the magnet. Functional images were realigned and co-registered to the T1-weighted structural image, followed by segmentation and normalization to the Montreal Neurological Institute (MNI) space. Spatial smoothing using Gaussian kernel (FWHM of 6 mm) was performed.
fMRI Localizer Motor Task
fMRI was used to localize the brain area activated in response to the experimental motor task. During the fMRI scan, participants were asked to squeeze the stress-ball on cue (with green screen indicating Go and red screen indicating Stop) for a total of 5 experimental blocks. Task-related activation was evaluated for each cue, and voxels with statistically significant activation (t-statistic) were used to construct an ‘activation map’. The statistical threshold for clusters was set at p<0.05 and corrected for multiple comparisons using family-wise error (FWE) over the whole brain (uncorrected threshold p<0.001) to minimize detection of false positives (type I error).
BBB permeability quantification
Analysis was performed as reported (Serlin et al., 2019; Veksler et al., 2020). Briefly, preprocessing included image registration and normalization to MNI space using SPM12. BBB permeability was calculated using in-house MATLAB script. A linear fit was applied to the later stage of the scan and the slope of contrast agent concentration changes over time was calculated for each voxel. Positive slopes reflected contrast agent accumulation due to BBB modulation. To compensate for physiological (e.g., heart rate, blood flow) and technical (e.g., injection rate) variability, slopes were normalized to the slope of each subject’s transverse sinus. Values were considered to reflect modulated BBB when exceeding the 95th percentile of the corresponding mean cumulative distribution function from healthy young controls. Voxels that reflect modulated BBB values were color-coded and superimposed on the anatomical scans to illustrate the kinetics, localization, and extent of BBB pathology.
Statistical analysis
All statistical tests were performed using Prism9 (GraphPad Software) and are indicated in the figure legends. Tests were two-tailed and corrected for multiple comparisons with two-stage linear step-up false discovery rate (FDR) procedure of Benjamini, Krieger and Yekutieli whenever applicable. Paired tests were used for inter-hemispheric comparisons. Non-parametric tests were used for data which was not normally distributed. Differences in fluorescence in Figure 1j-k were tested using a nested t-test. Differences in BBB permeability index, area under the curve and maximum amplitude of the somatosensory evoked potential before and after stimulation, were tested using the Wilcoxon matched-pairs signed rank test. Differences between groups in ELISA, area under the curve and maximum amplitude of the SEP were analyzed using Kruskal-Wallis test with FDR correction. Differences in the amount of DEGs in Figures 4c-d, S4c-e were analyzed using chi-square contingency test. Differences in BBB permeability across regions in the human brain were analyzed using two-way ANOVA with FDR correction. P≤ 0.05 was defined as the statistical significance level.
Acknowledgements
This study was supported by the Israel Science Foundation (Awards No. 717/15 and 2254/20), Canadian Institutes of Health Research (CIHR Award No. PJT 148896), and the European Research Area Network Neuron Award (Award No. NDD 168164). We thank QB3 Genomics Sequencing Laboratory, UC Berkeley, Berkeley, CA, RRID:SCR_022170 for using their facilities.
Declaration of Interest
The authors declare no competing interests.
Data and Code Availability
Source Data for quantification described in the text or shown in Figs. 1–4 and Extended Data Figure 1–4 are available with the paper. Additional data are available upon reasonable requests. Code and custom MATLAB scripts for analyses used in this study is also available upon reasonable request from the corresponding authors.
References
- Effects of anesthesia on BOLD signal and neuronal activity in the somatosensory cortexJ. Cereb. Blood Flow Metab 35:1819–1826https://doi.org/10.1038/jcbfm.2015.130
- Blood-brain barrier permeability is regulated by lipid transport-dependent suppression of caveolae-mediated transcytosisNeuron 94https://doi.org/10.1016/J.NEURON.2017.03.043
- Extensive piano practicing has regionally specific effects on white matter developmentNat. Neurosci 8:1148–1150https://doi.org/10.1038/nn1516
- Regulation of Blood-Brain Barrier Transporters by Transforming Growth Factor-β/Activin Receptor-Like Kinase 1 Signaling: Relevance to the Brain Disposition of 3-Hydroxy-3-Methylglutaryl Coenzyme A Reductase Inhibitors (i.eStatins). Drug Metab. Dispos 50:942–956https://doi.org/10.1124/dmd.121.000781
- Ultra-fast multispectral optical imaging of cortical oxygenation, blood flow, and intracellular calcium dynamicsOpt. Express 17https://doi.org/10.1364/oe.17.015670
- Transcriptome profiling reveals TGF-β signaling involvement in epileptogenesisJ. Neurosci 29:8927–8935https://doi.org/10.1523/JNEUROSCI.0430-09.2009
- The Molecular Constituents of the Blood-Brain BarrierTrends Neurosci 38:598–608https://doi.org/10.1016/j.tins.2015.08.003
- Caveolae in CNS arterioles mediate neurovascular couplingNature 579:106–110https://doi.org/10.1038/s41586-020-2026-1
- Rapid plasticity of human cortical movement representation induced by practiceJ. Neurophysiol 79:1117–1123https://doi.org/10.1152/JN.1998.79.2.1117/ASSET/IMAGES/LARGE/JNP.JA47F4.JPEG
- Environmental enrichment recruits activin A to recalibrate neural activity in mouse hippocampusCereb. Cortex https://doi.org/10.1093/cercor/bhac092
- Astrocytic dysfunction in epileptogenesis: Consequence of altered potassium and glutamate homeostasis?J. Neurosci 29:10588–10599https://doi.org/10.1523/JNEUROSCI.2323-09.2009
- The role of NMDA receptors in information processingAnnu. Rev. Neurosci 16:207–222https://doi.org/10.1146/annurev.ne.16.030193.001231
- Changes in grey matter induced by trainingNature 427:311–312https://doi.org/10.1038/427311a
- Physiological effects of enriched environment exposure and LTP induction in the hippocampus in vivo do not transfer faithfully to in vitro slicesLearn. Mem 17:480–484https://doi.org/10.1101/lm.1822610
- PSD-95 involvement in maturation of excitatory synapsesScience (80-.) 290:1364–1368https://doi.org/10.1126/science.290.5495.1364
- Synaptic mechanisms for plasticity in neocortexAnnu. Rev. Neurosci 32:33–55https://doi.org/10.1146/annurev.neuro.051508.135516
- The effect of different anesthetics on neurovascular couplingNeuroimage 51:1367–1377https://doi.org/10.1016/j.neuroimage.2010.03.060
- Identification of 1,5-naphthyridine derivatives as a novel series of potent and selective TGF-β type I receptor inhibitorsJ. Med. Chem 47:4494–4506https://doi.org/10.1021/jm0400247
- The role of Smad2 in adult neuroplasticity as seen through hippocampal-dependent spatial learning/memory and neurogenesisJ. Neurosci 41:6836–6849https://doi.org/10.1523/JNEUROSCI.2619-20.2021
- Visualizing synaptic plasticity in vivo by large-scale imaging of endogenous ampa receptorsElife 10https://doi.org/10.7554/ELIFE.66809
- Differential effects of NMDA and AMPA glutamate receptors on functional magnetic resonance imaging signals and evoked neuronal activity during forepaw stimulation of the ratJ. Neurosci 26:8409–8416https://doi.org/10.1523/JNEUROSCI.4615-05.2006
- Electrophysiological evidence for cortical plasticity with movement repetitionEur. J. Neurosci 21:2271–2277https://doi.org/10.1111/J.1460-9568.2005.04045.X
- Long-term synaptic plasticity in rat barrel cortexCereb. Cortex 25:2741–2751https://doi.org/10.1093/cercor/bhu071
- Blood-brain barrier dysfunction, TGFβ signaling, and astrocyte dysfunction in epilepsyGlia 60:1251–1257https://doi.org/10.1002/glia.22311
- Experience-dependent structural synaptic plasticity in the mammalian brainNat. Rev. Neurosci 10:647–658https://doi.org/10.1038/nrn2699
- TGF-β receptor-mediated albumin uptake into astrocytes is involved in neocortical epileptogenesisBrain 130:535–547https://doi.org/10.1093/brain/awl317
- Quantitative analysis of albumin uptake and transport in the rat microvessel endothelial monolayerAm. J. Physiol. - Lung Cell. Mol. Physiol 284:L187–L196https://doi.org/10.1152/ajplung.00152.2002
- Neuronal regulation of the blood–brain barrier and neurovascular couplingNat. Rev. Neurosci 21:416–432https://doi.org/10.1038/s41583-020-0322-2
- Cholesterol is required for surface transport of influenza virus hemagglutininJ. Cell Biol 140:1357–1367https://doi.org/10.1083/jcb.140.6.1357
- TGFβ signaling is associated with changes in inflammatory gene expression and perineuronal net degradation around inhibitory neurons following various neurological insultsSci. Rep 7https://doi.org/10.1038/s41598-017-07394-3
- Essential role for cholesterol in synaptic plasticity and neuronal degenerationFASEB J 15:1858–1860https://doi.org/10.1096/fj.00-0815fje
- Stimulation of the Sphenopalatine Ganglion Induces Reperfusion and Blood-Brain Barrier Protection in the Photothrombotic Stroke ModelPLoS One 7https://doi.org/10.1371/journal.pone.0039636
- Characterization of CBF response to somatosensory stimulation: Model and influence of anestheticsAm. J. Physiol. - Hear. Circ. Physiol 264:223–1228https://doi.org/10.1152/ajpheart.1993.264.4.h1223
- Epileptiform activity and spreading depolarization in the blood-brain barrier-disrupted peri-infarct hippocampus are associated with impaired GABAergic inhibition and synaptic plasticityJ. Cereb. Blood Flow Metab 37:1803–1819https://doi.org/10.1177/0271678X16652631
- Anesthesia and the quantitative evaluation of neurovascular couplingJ. Cereb. Blood Flow Metab 32:1233–1247https://doi.org/10.1038/jcbfm.2012.50
- Functional Plasticity in Somatosensory Cortex Supports Motor Learning by ObservingCurr. Biol 26:921–927https://doi.org/10.1016/j.cub.2016.01.064
- TGFβ1 exacerbates blood-brain barrier permeability in a mouse model of hepatic encephalopathy via upregulation of MMP9 and downregulation of claudin-5Lab. Investig 95:903–913https://doi.org/10.1038/labinvest.2015.70
- Long-term plasticity in mouse sensorimotor circuits after rhythmic whisker stimulationJ. Neurosci 29:5326–5335https://doi.org/10.1523/JNEUROSCI.5965-08.2009
- Neuronal Activity Drives Localized Blood-Brain-Barrier Transport of Serum Insulin-like Growth Factor-I into the CNSNeuron 67:834–846https://doi.org/10.1016/j.neuron.2010.08.007
- Development, maintenance and disruption of the blood-brain barrierNat. Med 19:1584–1596https://doi.org/10.1038/nm.3407
- Astrocyte-derived small extracellular vesicles promote synapse formation via fibulin-2-mediated TGF-β signalingCell Rep 34https://doi.org/10.1016/j.celrep.2021.108829
- Dynamic in vivo imaging of cerebral blood flow and blood-brain barrier permeabilityNeuroimage 49:337–344https://doi.org/10.1016/j.neuroimage.2009.08.009
- Neuronal Activity Regulates Blood-Brain Barrier Efflux Transport through Endothelial Circadian GenesNeuron 108:937–952https://doi.org/10.1016/j.neuron.2020.09.002
- Purification of RNA using TRIzol (TRI Reagent)Cold Spring Harb. Protoc 5https://doi.org/10.1101/pdb.prot5439
- Learning in the Fast Lane: New Insights into NeuroplasticityNeuron 73:1195–1203https://doi.org/10.1016/j.neuron.2012.01.025
- Synaptic plasticity in area CA1 of rat hippocampal slices following intraventricular application of albuminNeurobiol. Dis 91:155–165https://doi.org/10.1016/j.nbd.2016.03.008
- Albondin-mediated capillary permeability to albumin. Differential role of receptors in endothelial transcytosis and endocytosis of native and modified albuminsJ. Biol. Chem 269:6072–6082https://doi.org/10.1016/s0021-9258(17)37571-3
- Monitoring stroke progression: In vivo imaging of cortical perfusion, blood-brain barrier permeability and cellular damage in the rat photothrombosis modelJ. Cereb. Blood Flow Metab 34:1791–1801https://doi.org/10.1038/jcbfm.2014.147
- TGF-Beta Modulates the Integrity of the Blood Brain Barrier In Vitro, and Is Associated with Metabolic Alterations in PericytesBiomedicines 11:1–19https://doi.org/10.3390/biomedicines11010214
- Lasting blood-brain barrier disruption induces epileptic focus in the rat somatosensory cortexJ. Neurosci 24:7829–7836https://doi.org/10.1523/JNEUROSCI.1751-04.2004
- Blood-Brain Barrier Leakage: A New Biomarker in Transient Ischemic AttacksStroke 50:1266–1269https://doi.org/10.1161/STROKEAHA.119.025247
- Mouse fMRI under ketamine and xylazine anesthesia: Robust contralateral somatosensory cortex activation in response to forepaw stimulationNeuroimage 177:30–44https://doi.org/10.1016/J.NEUROIMAGE.2018.04.062
- The role of lipid species in membranes and cancer-related changesCancer Metastasis Rev 39:343–360https://doi.org/10.1007/s10555-020-09872-z
- Meta-analysis of RNA-seq expression data across species, tissues and studiesGenome Biol 16https://doi.org/10.1186/s13059-015-0853-4
- CLN3 Loss Disturbs Membrane Microdomain Properties and Protein Transport in Brain Endothelial CellsJ. Neurosci 33:18065–18079https://doi.org/10.1523/JNEUROSCI.0498-13.2013
- Gp60 Activation Mediates Albumin Transcytosis in Endothelial Cells by Tyrosine Kinase-dependent PathwayJ. Biol. Chem 272:25968–25975https://doi.org/10.1074/JBC.272.41.25968
- A molecular atlas of cell types and zonation in the brain vasculatureNature 554:475–480https://doi.org/10.1038/nature25739
- Glutamate-mediated blood–brain barrier opening: Implications for neuroprotection and drug deliveryJ. Neurosci 36:7727–7739https://doi.org/10.1523/JNEUROSCI.0587-16.2016
- Slow blood-to-brain transport underlies enduring barrier dysfunction in American football playersBrain 143:1826–1842https://doi.org/10.1093/brain/awaa140
- A proportional but slower NMDA potentiation follows AMPA potentiation in LTPNat. Neurosci 7:518–524https://doi.org/10.1038/nn1220
- Albumin induces excitatory synaptogenesis through astrocytic TGF-β/ALK5 signaling in a model of acquired epilepsy following blood-brain barrier dysfunctionNeurobiol. Dis 78:115–125https://doi.org/10.1016/j.nbd.2015.02.029
- Physiological blood–brain transport is impaired with age by a shift in transcytosisNature 583:425–430https://doi.org/10.1038/s41586-020-2453-z
- A Circadian Clock in the Blood-Brain Barrier Regulates Xenobiotic EffluxCell 173:130–139https://doi.org/10.1016/j.cell.2018.02.017
- Visualization of NMDA receptor - Dependent AMPA receptor synaptic plasticity in vivoNat. Neurosci 18:402–409https://doi.org/10.1038/nn.3936
- Metascape provides a biologist-oriented resource for the analysis of systems-level datasetsNat. Commun 10:1–10https://doi.org/10.1038/s41467-019-09234-6
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Reviewed Preprint version 3:
- Version of Record published:
Copyright
© 2023, Swissa et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.