2.1.

Animal Preparation

All experimental procedures were approved by the Weill Cornell Medical College Animal Care and Use Committee following NIH guidelines. Adult male Sprague-Dawley rats (220 to 380 g) were anesthetized with isoflurane in $70\%\mathrm{N}2:30\%\mathrm{O}2$, 4% induction, and 1.5 to 2% maintenance. Body temperature was maintained at 37°C with a regulated heating blanket (Harvard Apparatus, Holliston, Massachusetts). The heart rate, pO2, and end tail carbon dioxide (ETCO2) were carefully monitored with a small animal Capnograph (Surgivet, Waukesha, Wisconsin) and were maintained stable throughout the experiment (heart rate: 250 to $300\text{pulse}/\min$, pO2 $>90\%$, ETCO2 $\sim 25$ to 28 mm Hg). The head was fixed in a stereotaxic frame throughout experimentation.

2.2.

Dye Staining

Calcium dye, Oregon Green 488 BAPTA-1 AM (OGB-1, Life Technology, Grand Island, New York; excitation $470\pm 10\mathrm{nm}$, emission $>510\mathrm{nm}$), was chosen for this study following the work of Bouchard and colleagues.²⁰ The dye was first dissolved in 5 μl of 20% Pluronic F-127 in dimethylsulfoxide (DMSO) (w/v) and then in 50 μl of artificial cerebrospinal fluid (ACSF) (containing 125 mM NaCl, 5 mM KCl, 10 mM D-glucose, 10 mM HEPES, 2 mM MgSO4, and 2 mM CaCl2) to achieve a final concentration of $0.91\mathrm{mg}/\text{mL}$. CED was employed to stain a large area of the cortex with calcium dye. A glass pipette (outside tip diameter 50 to 100 μm) was filled with the calcium dye solution and inserted into the cortex to a depth of 1 mm from the surface through a small hole created in the skull. The dye was injected with a micropump (UMP-3, World Precision Instruments Inc., Sarasota, Florida) at a speed of $100\mathrm{nl}/\min$ for 80 min (8 μl in total). The glass pipette was removed 10 min following the injection allowing the animal to rest for 1 h to permit the dye to load into the cells. In order to evaluate loading efficacy, a subgroup of animals ($n=5$) was similarly labeled and used for histology. For in vivo imaging, after staining, two craniotomy windows were opened over both hemispheres, and the exposed brain was covered with silicon oil (12,500 centistoke) to preserve cortical moisture. For simultaneous VSD and calcium imaging, the VSD (RH-1691; excitation $625\pm 10\mathrm{nm}$, emission $>655\mathrm{nm}$) was dissolved in OGB-1 solution ($1\mathrm{mg}/\text{mL}$). Using the micropump and the same rate, 1 μl of the combined dye was injected into the cortex in order to stain a small region of the cortex.

To record the activity of astrocytic glia only (Fig. 7), we employed the staining protocol as described by Ghosh et al.²⁴ In brief, the calcium dye, rhod-2 (rhod-2, AM, Life Technology), was first dissolved in 5 μl of 20% Pluronic F-127 in DMSO (w/v) and then in 50 μl of ACSF ($0.91\mathrm{mg}/\text{mL}$). The dye solution was topically applied to the exposed cortical surface for 90 min and then the brain was washed with dye-free ACSF for 15 min. This technique allows staining of only the glia and real-time imaging of glial activity.

2.3.

Optical Imaging

For wide-field calcium imaging, a CCD camera (Dalsa camera in Imager 3001, Optical Imaging, Rehovot, Israel) using a tandem lens ($50\times 50\mathrm{mm}$) arrangement²⁵ was focused 300–400 μm below the cortical surface. A $470\pm 10$- or $570\pm 10\text{-}\mathrm{nm}$ light-emitting diode (LED) was employed as the illumination source for OGB-1 or rhod-2, respectively. The illumination was guided to the cortex with an optical fiber. A 510-nm long-pass or 590-nm long-pass filter was placed before the camera for OGB-1 or rhod-2 imaging, respectively.

The setup for simultaneous VSD and calcium imaging (OGB-1) was modified from our previous simultaneous VSD and IOS system²⁶ [Fig. 5(a)]. Two LEDs, $470\pm 10$ and $625\pm 10\mathrm{nm}$, coupled with a dichroic mirror were employed as illumination sources. The illumination was guided by optical fiber onto the brain surface. The fluorescence from the calcium dye (510 nm) and VSD ($>665\mathrm{nm}$) were wavelength-separated using a 570-nm dichroic mirror. The VSD images were captured using a MiCAM ULTIMA CMOS camera (SciMedia, Costa Mesa, California). The camera has a temporal resolution of 400 Hz and a spatial resolution of $100\times 100\text{pixels}$, with each pixel sampling an area of $50\mu \text{m}\times 50\mu \text{m}$. The calcium image was recorded with a CCD camera, with a temporal resolution of 110 Hz and a spatial resolution of $256\times 256\text{pixels}$, with each pixel sampling an area of $\sim 40\mu \text{m}\times 40\mu \text{m}$ (Dalsa camera in Imager 3001, Optical Imaging, Rehovot, Israel).

For simultaneous calcium (OGB-1) and IOS imaging, we upgraded Bouchard et al.’s setup [Fig. 6(a)]. Three LEDs at wavelengths of 470, 530, and 625 nm were employed as the illumination sources for calcium, total hemoglobin (Hbt), and deoxyhemoglobin (Hbr), respectively.^27,28 The LEDs were coupled with two dichroic mirrors (490 and 590 nm) and guided to the cortex using optical fibers. A 510-nm long-pass filter was placed before the camera. Since the hemodynamic signal is slower than the calcium signal, the illumination rate of the calcium signal was set at double that of the hemodynamic illumination, providing frame rates of 55 Hz for calcium and 27.5 Hz for IOS.

A modified Beer-Lambert law was used to calculate Hbr and Hbt changes from the 570- and 610-nm data using Eq. (1).^29,30

2.4.

Confocal Microscope and Histology

A confocal microscope was employed to examine the calcium dye loading at the single-cell level. Sulforhodamine 101 (SR101) $0.25\mathrm{mg}/\text{mL}$ was mixed with OGB-1 for injection in order to distinguish neurons and astroglia.³¹ Transcardiac perfusion fixation with 4% paraformaldehyde was performed following 0.01% sodium nitrite in phosphate buffered saline (pH 7.4) flush. The rat brain was further fixed in 4% paraformaldehyde at room temperature overnight. Coronal slices of 200 μm crossing the dye injection site was prepared. A confocal microscope (Leica TCS SP5, Buffalo Grove, Illinois) was used to image with single-cell resolution. Lasers of 488 and 543 nm were employed as the excitation wavelength for OGB-1 and SR101, respectively. The fluorescence emission was collected through 535-nm (green, for OGB-1) or 645-nm (red, for SR101) filters.

Standard hematoxyline-eosine (H&E) staining was employed to exam the potential damage created by the CED protocol. After perfusion, the rat brain was further fixed in 4% paraformaldehyde at room temperature for 48 h. The rat brain was then infiltrated and embedded with paraffin for histology processing. Standard H&E staining was performed on the 7-μm coronal slice crossing the dye injection site.

2.5.

Electrophysiology and Triggering of Epileptiform Events

To record local field activity, a glass microelectrode (2 to 4 MΩ) filled with saline was positioned $\sim 500\mu \text{m}$ deep from the surface. The local field potential (LFP) was amplified and band-pass filtered (1 to 500 Hz) using a DAB-S system (World Precision Instruments, Sarasota, Florida) and digitized at 1000 Hz by a CED Power 1401. The LFP data were recorded with a PC running the Spike2 software (Cambridge Electronic Design, Cambridge, United Kingdom). To trigger epileptiform events, seizure-like discharges were induced by injecting 4-Aminopyridine (4-AP, 15 mM, 0.5 μl) 300 to 500 μm into the cortex, through a second glass microelectrode positioned 500 μm away using a Nanoject II injector (Drummond Scientific, Broomall, Pennsylvania). Multifocal interictal-spike (IIS)-like events were triggered with either a topical application of 1-mM bicuculline (Fig. 2) or a focal injection (0.5 μl) 300 to 500 μm deep into the cortex (Fig. 5).

2.6.

Somatosensory Stimulation

Two stainless-steel electrodes (insulated to within 2 mm of the tip) were placed on the right hind-paw. A 2-s train of electrical pulses (5 Hz, 2 mA intensity, and 0.2 ms pulse width) was delivered in 10 stimulation trials with an intertrial interval of 30 s. The stimulation parameters and delivery were controlled by a stimulation box coupled to Master-9 software and hardware (A.M.P.I, Jerusalem, Israel). Each trial was of 3.5 s duration with stimulation onset at 0.5 s. Spatiotemporal calcium imaging data for all trials, in each animal, were averaged to create a mean trial and were subjected to statistical parametric mapping (SPM), in which each pixel’s time-series was regressed against a design matrix representing a direct current offset and a boxcar function of the same duration as the stimulation. This generated a z-score activation map where high values represented large evoked calcium signal increases due to hind-paw stimulation. Pixels possessing $>80\%$ of the maximum z-score were used to define a region of interest (ROI). The time-series from each associated pixel was then averaged.

The SPM is a standard approach to identify foci of significant changes in cerebral activity using SPM.³² The SPM approach is not parameterized, but the z-score output is sensitive to the choice of a design matrix. However, in cases such as ours, in which the stimulus onset and offset is well defined, changing design matrix parameters (e.g., using a Gaussian function instead of a boxcar) has qualitatively little effect on the spatial map of activation.

3.1.

Convection-Enhanced Staining of Calcium Dyes

Traditional multi-injection calcium staining creates inhomogeneous loading with islands of high concentration that decrease with distance from the injection site as the dye diffuses into the surrounding extracellular space [Fig. 1(b)] and multiple injections are employed to stain a large cortical area, which often results in an uneven staining [Fig. 1(b)]. CED, on the other hand, delivers a homogeneous concentration from the injection site all the way out to the advancing wavefront of the solvent as long as the pressure gradient is actively maintained [Fig. 1(a)]. Through a single injection, using our injection parameters, we were able to routinely stain a large cortical area, up to 6 mm from the injection [Fig. 1(c)].

Fig. 1

Convection-enhanced delivery and bulk-loading of calcium-sensitive dye. (a) Schematic diagram of convection-enhanced delivery of calcium-sensitive dye. The calcium-sensitive dye is injected into the cortex via the electrode. The dye is driven by the pressure and propagates in spherical space. A large volume of calcium dye is injected into the cortex and a large area of the cortex is stained (up to 6 mm). Examples of (b) the older multiple injection technique for bulk-loading and (c) convection-enhanced delivery (CED) technique stained brain. The images were taken $\sim 300\mu \text{m}$ below the surface. Note that the CED technique provides smooth loading to a broad area that slowly decreases with distance from the injection site.

Wide-field fluorescence imaging was employed to investigate the distribution of OGB-1 within the brain. In a coronal slice through the injection site [Fig. 2(b)], calcium fluorescence spread several millimeters past the injection site in all directions with higher concentrations at the injection site itself. At 1 to 2 mm from the injection site, the calcium fluorescence showed a much more even distribution [Figs 2(b), 2(c), and 2(d)]. Vertically, the dye clearly stained all cortical layers [Figs. 2(b), 2(c), and 2(e)]. The total volume of stained neocortical tissue is estimated to be 75 to $170{\mathrm{mm}}^3$.

We used curve fitting to analyze the fluorescent intensity distribution in Fig. 2(d) and found that the curve follows exponential decay with an equation of

Confocal microscopy was employed to verify the staining result at a cellular level. The confocal images were taken at different distances from the injection site, with green fluorescence (from OGB-1) overlaying red fluorescence (from SR101) from the same field of view [Fig. 2(f)]. The neurons were labeled with only OGB-1 (green cells), while astrocytes were co-labeled with both OGB-1 and SR101 [orange cells, Fig. 2(f)]. Close to the injection site (injection site and 1 mm), strong green fluorescence could be recorded from the background, indicating all compartments were heavily stained. The nucleus could be identified as dark areas since the calcium dye was not absorbed by the nucleus. A bright halo (slightly brighter than the background) surrounding the nucleus represents the soma. The glia could be clearly identified as orange cells. Far from the injection site (2, 3, and 4 mm), the background calcium fluorescence was weaker, and both the neurons and glia could be identified. The calcium fluorescence from the superficial neuropil could also be recorded. Clear labeling of neurons, glia, and neuropil indicates that the resulting calcium signal would be a mixture of activation from all three tissue compartments.

Standard H&E staining was employed to investigate the potential damage to cortical tissue by the CED protocol [Fig. 2(g)]. The coronal slice, $\sim 50\mu \text{m}$ away from the injection site, did not show obvious tissue damage.

We further investigated if neuronal activity could be recorded from the entire area of stained cortex. Bicuculline (BMI, 1mM) was topically applied to the whole brain and interictal spikes recorded from the entire extent of the exposed brain [Fig. 2(h), top]. The signal-to-noise ratio (s/n) was $>3$ in most areas [Fig. 2(h), bottom].

Fig. 2

Distribution of calcium dye in the cortex. (a) CED loading with Oregon Green 488 BAPTA-1 AM (OGB-1 AM) injected into one hemisphere (red dot, right panel) creates continuous loading up to the edge of the craniotomy. The entire left hemisphere is shown on the left (530-nm illumination) for spatial comparison. Two green lines indicate the locations of the coronal slice shown in (b) and (c). (b) Note that although the calcium fluorescence is higher immediately around the injection site, there is widespread staining up to $\sim 6\mathrm{mm}$ away in all directions. (c) In a coronal slice $\sim 3\mathrm{mm}$ away from the injection site, a 5- to 6-mm area of cortex remains stained in a homogeneous fashion. Also note that staining is continuous from layers 1 to 6 with abrupt cessation at the pial surface and ependyma of the ventricle where resistance is met. (d) The average horizontal distribution of calcium fluorescence in the cortex ($n=5$ rats). The red line shows an exponential decay fitting of the fluorescence curve. (e) The averaged vertical distribution of calcium fluorescence in the cortex ($n=5$ rats). (f) Confocal images taken from a coronal slice crossing the injection site. The images were taken at different distances from the injection site. Each image shows the fluorescence from the superficial 300 μm of cortex. The green fluorescence from OGB-1 is overlaid with the red fluorescence from sulforhodamine 101. The neurons are shown in green cell and the astrocytes are orange. The majority of superficial fluorescence comes from the neuropil. However, in the deeper layers, both neurons and glia are equally stained. (g) Standard H&E staining of a coronal slice $\sim 50\mu \text{m}$ from the injection site—top: the injection site and bottom: 3 to 6 mm away from the injection site. The black box in the left images indicates the area zoomed-in on the right. (h) An interictal spike induced with topical bath application of 1 mM BMI—top: the spatial distribution of calcium waveform during an interictal spike and bottom: the signal-to-noise ratio of each pixel during the ictal event. The red dot indicates the OGB-1 injection site.

3.2.

Calcium Imaging of Spontaneous Brain Activity

We also wished to evaluate whether the trauma from the CED method might alter normal cortical processing. First, we examined spontaneous brain activity, which also allowed us to measure the magnitude of the resulting signal and its s/n without any averaging-between events. Using LFP measurements of spontaneous synaptic activity, we measured simultaneous wide-field calcium fluctuations. LFP showed population activity of different amplitudes and frequencies during the spontaneous activity. The calcium (OGB-1) waveform from the same location showed similar fast-frequency activity riding on slower-frequency waves, and the fast activity appeared remarkably similar to the LFP signal (Fig. 3). A correlation coefficient between the LFP and calcium power traces from the brain ($0.46\pm 0.07$) was significantly higher than those on the bone [$0.16\pm 0.04$, $p<0.001$, Fig. 3(c)]. Indeed, the sensitivity of the technique is such that every spontaneous LFP fluctuation, each with its own morphology, also has different spatiotemporal calcium dynamics [Fig. 3(d)]. The sensitivity and reliability of the calcium signal to represent the synaptic activity of the LFP indicates that broad-field calcium imaging could be used as a high-resolution optical electrocorticography.

Fig. 3

Wide-field calcium imaging of spontaneous activity. (a) The field of view of the calcium imaging. A small green square adjacent to the tip of the local field potential (LFP) electrode on the calcium image indicates a region of interest (ROI) from which the calcium trace is derived. Note that the surrounding bone is not visible as the signal is too attenuated. (b) The simultaneously recorded LFP and calcium traces from two ROIs. The green ROI was adjacent to the LFP electrode and the blue ROI was taken from the adjacent bone. The left shows the raw LFP and calcium signals and the right shows the power of the signals. A black box over the LFP indicates the section zoomed in (d). Note that LFP clearly correlates with the faster activity in the calcium tracing. (c) The pseudo-color image of correlation coefficient between the power traces of LFP and calcium signal from each pixel. A high correlation is seen over the entire calcium labeled cortex whereas there is low correlation from the surrounding bone, indicating that the calcium signal effectively demonstrates synaptic activity throughout the cortex. (d) The spatial propagation of the calcium signal during three spontaneous spikes. Vertical bars over the LFP indicate the time points for the calcium images. Interframe interval: 36 ms. Note that the propagation patterns are different for each LFP event, each of which has its own morphology.

3.3.

Application of Wide-Field Calcium Imaging During Sensory Stimulation

In order to see if the CED technique produced calcium responses that were consistent with the multiple injection technique and to provide baseline normal calcium dynamics for comparison with pathological activity, we measured wide-field CED-loaded calcium changes (OGB-1) to hind-paw stimulation (2 s, 5 Hz, 0.2 ms) using SPM analysis [Fig. 4(a)]. A discrete area of robust activation was identified with an average maximal peak of a $3.4\pm 0.2\%$ change from baseline. Each stimulation pulse elicited a discernible transient change in the calcium signal amplitude, which underwent adaptation over time [$N=4$, Fig. 4(c)]. The spatial spread of the stimulation-induced brain activity can also be recorded with calcium imaging [Fig. 4(d)]. These results compare favorably to previous reports of calcium signal changes to somatosensory stimulation using multiple injection bulk-loading of calcium dye, where amplitude of stimulation has been reported at $\sim 8\%$,^19,33 providing further evidence that our single-injection technique is a reliable method by which to record wide-field calcium dynamics in the neocortex. The averaged IOS image (570 nm) 5 to 7 s post-stimulation showed hemodynamic changes in a similar area as the neuronal activity [Fig. 4(e)].

Fig. 4

Calcium imaging of somatosensory responses to hind-paw stimulation. (a) Statistical parametric map (SPM) of cortical calcium responses to 2 s hind-paw stimulation in a representative animal (hot colors denote large z-scores associated with robust increases in calcium signal). (b) Digital image of the cortical surface in the same animal as (a) coregistered with SPM-derived ROI (defined as pixels with $>80\%$ of maximum z-score, perimeter as white dashed line). Black triangle in (a) and (b) indicates location of single calcium dye injection in this animal. (c) Averaged calcium signal time-series ($N=4$) from SPM-derived ROIs, showing robust increases in response amplitude to stimulation and subsequent adaptation (as percentage change from baseline). Gray dashed line indicates stimulation onset. Error bars are SEM. Ten red bars indicate the time of the peak response to the stimulation and the images are shown in (d). (d) The spatial distribution of the stimulation-induced activity. Note that the first two stimuli induce activity in a broader area than later stimuli. (f) The spatial distribution of 530-nm intrinsic optical signal. The image is an average of the recording 5 to 7 s post the stimulation. The intrinsic signal showed similar spatial distribution as the SPM map.

3.4.

Wide-Field Calcium Dye and VSD Imaging

In order to investigate the major sources of the calcium signal from this technique, we investigated the degree of overlap between the calcium signal and the VSD signal. The VSD signal is known to represent presynaptic subthreshold activity, predominantly from dendrites.^13–15 For this experiment, we performed simultaneous calcium and VSD imaging of epileptiform IIS. As previously described in our laboratory,^34,35 IIS were induced with local intracortical bicuculline application. As is apparent (Fig. 5), the waveforms of VSD and calcium signals are qualitatively very similar. Both showed identical onset times and monophasic changes during IIS and were temporally locked to the LFP recorded at the same location [Fig. 5(b), ROI 1). The calcium waveforms were found to be significantly broader ($466\pm 110\mathrm{ms}$) than the VSD waveforms ($302\pm 97\mathrm{ms}$) ($p<0.001$, $n=30$ IIS from four rats), indicating longer decay kinetics. We fitted the decay of the fluorescence transient with a single-exponential function, yielding a decay time constant of $0.17\pm 0.03\mathrm{s}$ for calcium ($n=30$) and $0.09\pm 0.02\mathrm{s}$ for VSD ($n=30$). The normalized fast Fourier transform of LFP, VSD, and calcium signals were also highly comparable, with the exception of low-frequency ($<1\mathrm{Hz}$) power in calcium signals being greater than that of VSD, likely due to greater contributions from slower glial activity in the calcium data [Fig. 5(c)]. Very little low-frequency power was observed in LFP due to the fact that a 1-Hz pass filter was built into the amplifier. In the higher-frequency range ($>2\mathrm{Hz}$), the LFP spectrum is almost identical to VSD and calcium signals. Cross-correlation between the LFP spectrum ($>2\mathrm{Hz}$) and the VSD/calcium spectra showed strong correlation (LFP versus VSD: $R=0.99$, $p<0.001$, LFP versus Ca: $R=0.99$, $p<0.001$). In the low-frequency range ($<2\mathrm{Hz}$), there is no correlation between the LFP spectrum and the VSD/calcium spectra (LFP versus VSD: $R=0.05$, $p=0.88$, LFP versus Ca: $R=0.04$, $p=0.91$). The normalized calcium and VSD waveforms became approximately equivalent following the application of a 1-Hz high-pass filter [Fig. 5(d), 1 Hz, first-order Butterworth high-pass filter, using MATLAB®]. The amplitudes of VSD and calcium signal during 30 IISs in four rats are shown in Fig. 5(e). The data points followed a linear distribution, with $R^2=0.76$, indicating that the amplitude of VSD and calcium signals were highly correlated ($p<0.001$). Our data suggest that the majority of calcium signals closely reflects subthreshold synaptic activity recorded from the neuropil, which was very similar to the VSD signal but with a higher amplitude and higher s/n as well as a slower component, likely reflecting the slow dynamic of calcium clearance and slower glial waves.

We also recorded the VSD and OGB-1 signal from different locations far away from the BMI application site [Fig. 5(b), ROIs 2 and 3]. The amplitude of VSD and calcium signal decreased with the increasing distance.

Fig. 5

Simultaneous calcium and VSD imaging of an interictal event. (a) Experimental setup for simultaneous VSD and calcium imaging. Two light-emitting diodes (LEDs), $470\pm 10$ and $625\pm 10\mathrm{nm}$, coupled with a dichroic mirror are employed as illumination sources. The illumination is guided with fiber optics onto the cortex. The fluorescence of OGB-1 AM and RH-1691 are collected with a camera lens and separated with a 570-nm dichroic mirror. A $530\pm 10\text{-}\text{nm}$ filter is placed in front of the calcium imaging camera and a 665-nm long-pass filter (RG-665) is placed in front of the VSD camera. (b) The field of view imaged with VSD and calcium imaging. Three small white boxes on the image indicate the ROIs in which the VSD and calcium tracings were derived. The simultaneously recorded LFP, VSD, and calcium traces from the ROIs are shown at the bottom. (c) The frequency distributions of LFP, VSD, and calcium signal. The inner panel is a magnification of the low-frequency section. Note that the frequency distributions are largely overlapping and the VSD and calcium signals show highest power in the low-frequency range (0 to 1 Hz), while LFP correlates more strongly in the higher-frequency range. (d) The normalized 1-Hz high-pass filter (first-order Butterworth filter) of the VSD and calcium trace. Note that the VSD and calcium traces are almost identical. (e) Plot of the VSD fluorescence changes versus calcium fluorescence changes. Each dot indicates the peak amplitude of an interictal spike recorded by both imaging methods. A linear fit of the amplitudes was extrapolated.

3.5.

Simultaneous Calcium and Intrinsic Imaging of Hemodynamics and Oxygenation

In order for wide-field calcium imaging to be useful for the investigation of neurovascular coupling at a network level, it needs to be combined with similar wide-field hemodynamic monitoring. Using the propagation of focal neocortical seizures as an example of a pathologic brain state that involves a widespread network and that requires broad-field imaging for effective mapping, we employed focal intracortical injections of 4-AP as previously described by our laboratory.^26,36,37 The duration and morphology of the resulting ictal event was recorded with an LFP electrode placed in the 4-AP injection site, which showed the onset, evolution, and termination of the ictal-like events, which usually last $\sim 60$ to 80 s. The calcium signal (OGB-1) demonstrates fast activity riding on a wave of slower activity. The high-frequency component ($>1\mathrm{Hz}$) closely reflected the LFP waveform, while the low-frequency component ($>1\mathrm{Hz}$) appears to be non-neural in origin, likely reflecting glial activity [Fig. 6(c)]. In order to measure concomitant hemodynamic events, intrinsic optical signals at 530 and 625 nm were measured simultaneously. Using a modified Beer-Lambert law, we calculated the change in total hemoglobin (Hbt), deoxygenated hemoglobin (Hbr), and oxygenated hemoglobin (HbO) for correlation with the calcium signal [Fig. 6(c)]. During an ictal event, an increase in Hbt and HbO was combined with a prolonged decrease in Hbr, consistent with the absence of an epileptic dip and an oversupply of oxygenated hemoglobin, as previously described in a fraction of seizures using this model^36,37 [Fig. 6(c)]. Although high-frequency calcium signals (representing predominantly subthreshold neuronal activity) fluctuated with millisecond resolution, the low-frequency calcium signal and hemodynamic changes fluctuated with a temporal resolution of hundreds of milliseconds. With simultaneous low-frequency and high-frequency calcium imaging and multiwavelength intrinsic signal imaging, one can correlate wide-field spatiotemporal alterations in neurovascular coupling between neuronal and glial activity and blood volume and oxygenation changes. [Fig. 6(d)].

In order to better understand the signal source of the low-frequency component recorded by OGB-1, we selectively stained the astrocytic glia with a different calcium dye, rhod-2, AM, following the technique of Ghosh and colleagues.²⁴ The ictal waveform of the astrocytic calcium change appeared as a slowly propagating wavefront. These results indicate that the slowly progressive wave underlying the fast activity we recorded with OGB-1 likely arose from glial activation.

Fig. 6

Simultaneous calcium and multiwavelength intrinsic optical signal (IOS) of an ictal event. (a) Diagram of the setup for simultaneous calcium-sensitive dye and IOS imaging. 470-, 530, and 625-nm LEDs are coupled with dichroic mirror to serve as an illumination source. These LEDs are set to be alternately triggered by the camera frame indicator to capture calcium images at 55 Hz and IOS at 27.5 Hz. A 510-nm long-pass filter is placed before the camera to filter out the calcium illumination (470 nm). (b) The field of view of the imaged cortex. A small box indicates an ROI, in which the calcium and intrinsic optical imaging traces are derived. The left box is located adjacent to the 4-AP/LFP electrode, which is the source of the epileptic focus. (c) The waveforms of different recording modalities during an ictal event derived from the ROI shown in (b). A red box over the LFP indicates a section zoomed in (d). Note that the calcium trace has two components: the high-frequency component ($>1\mathrm{Hz}$) is closely related with LFP, but the low-frequency component ($<1\mathrm{Hz}$) does not reflect the dynamics of the LFP. (d) The spatiotemporal distribution of the high- and low-frequency calcium signals as well as the hemodynamic components. The vertical bars overlapping the LFP indicate the timing for the calcium imaging frames. Frame interval: 18 ms. The spatial distributions of low-frequency calcium signals and each of the hemodynamic signals are shown at the bottom.

4.1.

Convection-Enhanced Delivery

CED was originally developed as a technique to deliver macromolecules, such as chemotherapy, into the brain to treat infiltrative brain diseases, such as malignant tumors, that spread throughout the brain parenchyma.²³ In this method, the dissolved drug is driven through the extracellular space by maintaining a pressure gradient that distributes the drug spherically from the site of the injection into the parenchyma, creating a homogeneous drug concentration that does not dissipate as it would using only diffusion. Cerebrospinal fluid is replaced by the injected fluid, and the drug is delivered at a homogeneous dose. In adopting the technique to deliver calcium dyes, a critical component to ensuring reliable delivery has been the placement and depth of the injection electrode. We found delivery to be optimal when the electrode tip was placed $\sim 800\mathrm{um}$ below the surface (i.e., middle of the neocortex). This prevented upward dye leakage through the injection track and minimized the possibility of downward penetration into the underlying ventricle. After the electrode was placed in the cortex, we waited for $\sim 10\min$ before injecting, so that the tissue could bind to the glass electrode and form a close seal that could withstand the injection pressure. Another key determinant of dye-staining success was the injection speed. Higher injection speeds resulted in higher injection pressures and increased the likelihood of leakage. Conversely, low injection speeds resulted in a limited area of spread, the dye being subject to diffusion and unnecessary prolongation of the experiment. Given the thickness of the neocortex, the injection speed was found to be optimal at $100\mathrm{nl}/\min$.

4.2.

Convection-Enhanced Delivery and Targeted Bulk-Loading of Calcium Dye

Targeted bulk-loading is a staining approach in which cells are stained by a brief bolus injection of the acetoxymethyl ester of a calcium indicator dye into the extracellular space. The injected dye diffuses into the cells and is hydrolyzed by intracellular esterases.⁴⁰ This method can be used to stain an area of a few hundred micrometers in diameter and allow two-photon laser-scanning microscopy of calcium transients. This approach has been successfully applied in different species at various developmental stages (mice, rats, cats, and zebrafish^{16,31,41–44}). CED is a significant advancement to focal targeted delivery since a large area (several millimeters in diameter and 75 to $170{\mathrm{mm}}^3$ in volume) of cortex can be homogeneously stained with minimal trauma to the cortex as opposed to multiple small islands of staining with multiple cortical punctures (Fig. 1). Given minimal trauma, CED also increases the possibility for restaining, which enables long-term study.²¹ A clear advantage of CED over the targeted bulk-loading method is that the stained area is much larger. For targeted bulk-loading, the dye can only propagate a few hundred micrometers from the injection electrode, meaning that possible damage from the multiple cortical penetrations could affect subsequent recordings. With CED, this possibility is avoided since the dye can propagate several millimeters away from the injection electrode following a single injection.

4.3.

Difference between VSD and Calcium-Sensitive Dye

A VSD signal represents membrane voltage change in a mixture of all cells stained, including neurons, glia, and blood vessels. However, the contribution of glia and endothelia to the VSD signal is minimal. Although in vivo imaging of a single pixel samples multiple compartments of a population of neurons, including axons, dendrites, and cell bodies, studies have shown that the majority of the signal arises from dendrites, whose membrane area is largest.¹⁵ For this reason, the VSD signal waveform mainly reflects the subthreshold synaptic activity in superficial cortical layers and is a summation of activity in $\sim 230\text{-}\mu \text{m}$ area.^13,45 In vivo calcium imaging using the multisite injection technique also stains multiple cell compartments, such as cell bodies, dendrites, axons, and glia.⁴¹ Although with single-cell resolution or direct intracellular injection, the majority of the signal arises from suprathreshold action potentials triggered by voltage-gated calcium channels,^17,22 when imaged at low spatial resolution following bulk-loading, the calcium signal is a mixture of all cell compartments, and a large component arises from the neuropil, predominantly from axons.¹⁶ Thus, the signal we recorded mainly reflects strong synaptic activation in the neuropil of the superficial cortical layers.¹⁹ Hence, in our experimental preparation, VSD and calcium signals reflect very similar aspects of cortical activity (Fig. 5). However, calcium imaging has a much higher s/n ratio. Although the calcium dye can act as a buffer, the concentration of dye for the bulk-loading technique is lower, which is an advantage over intracellular injections. The other main difference between calcium and VSD imaging is that using calcium dye the slow glial component is much more apparent.¹⁸

The contribution of glial activation is quite different when comparing VSD with calcium-sensitive dyes. With VSDs, glial activity does not contribute significantly to the signal^14,15,17 because glial responses are weak (depolarized only 1 to 7 mV in response to visual stimulation).⁴⁶ However, the contribution of glia to the calcium signal is substantial, albeit on a much closer timescale than neuronal activity (Fig. 7). Schulz and colleagues demonstrated that the activation of glia introduces a prolonged low-frequency component to the calcium signal during electrical forepaw stimulation,¹⁹ which is similar to our recording. Therefore, under the experimental conditions we employ, calcium-dye-derived data appear to be a combination of high-frequency subthreshold synaptic activity derived from the neuropil and low-frequency glial activity. A high-pass filter (1 Hz) is sufficient to remove the glial component (Fig. 6), suggesting that calcium-sensitive dyes can provide the same two-dimensional activity maps as VSD imaging with significantly higher s/n ratio. Clearly, another main advantage of calcium imaging is the ability to image spiking in individual neurons, but this requires two-photon imaging of single cells, so the signal from the neuropil does not dominate the signal as it does with wide-field imaging.

Fig. 7

Astrocyte activity during an ictal event. (a) The field of view stained with rhod-2. Three small white dots indicate three ROIs in which the calcium tracings were derived. (b) The LFP and the astrocyte calcium activation recorded from the ROIs.

4.4.

Simultaneous Imaging of Neuronal Activity and Hemodynamic Changes

The study of neurovascular coupling requires, at a minimum, simultaneous recordings from both key components of the neurovascular unit, namely the neuronal as well as the hemodynamic element. Currently, a large number of methods are available to record hemodynamic changes on a wide variety of scales, from the movement of an individual blood cell or a single arteriole dilatation to the whole brain, such as fMRI, single photon emission computed tomography (SEPCT), and near-infrared spectroscopy. There are an equal number of methods available to record the voltage changes or electric currents on a wide variety of scales from single ion channel neuronal activity and intracellular recording to large-scale magnetoencephalography (MEG) and electroencephalography. However, the ability to record from a large area of cortex (several millimeters) with high temporal (1 ms) and spatial (50 micron) resolution can only be done with optical techniques. In this paper, we describe a reliable method to achieve this goal, which will be critical to understand neurovascular coupling for brain processes that involve wide-scale dynamic networks that cannot be easily averaged together to increase signat-to-noise ratio, such as epileptiform events or higher-order cognitive processing.