Detection of astrocyte regions of interest (ROI) containing Ca2+ fluctuations was performed in a semiautomated manner using the GECIquant program developed using the open source ImageJ analyses platform. The same procedure was followed for brain slice and in vivo data. The GECIquant program is implemented in Java based ImageJ script language and runs as a plugin on ImageJ. The input to GECIquant is a confocal 2D fluorescence image stack (8 or 16 or 32-bit) of arbitrary frame size, a user defined sampling rate and with time as the third dimension (t-stack). Data outputs of GECIquant include ROI intensity changes in time, ROI areas and centroid distances of each ROI from a reference ROI. Graphical outputs of GECIquant include ROI intensity kymographs and sub-stacks consisting of fluctuations. Supplementary Information 1 provides the script and Supplementary Information 2 provides a user manual for GECIquant.
Having analyzed all the data shown in this study, we observed three distinct types of spontaneous subcellular Ca2+ fluctuations within astrocytes, which we describe below and then clarify how they were detected within GECIquant. We classify Ca2+ fluctuations as: (1) somatic fluctuations that occur within the somata (these are restricted to somata and initial segments of processes arising from somata), (2) waves that occur exclusively within astrocyte processes and display repeated wave-like expansions and contractions of Ca2+, and (3) microdomain Ca2+ fluctuations that display highly restricted areas in astrocyte processes. These do not expand or contract as waves and remain restricted. The distinct areas covered by these three types of fluctuations are reported in the main text and inSupplementary Fig 5 . In this section, we describe a semi-automated method to accurately capture regions of interest (ROIs) for somatic, wave and microdomain Ca2+ fluctuations within astrocytes.
After an image series was acquired (e.g.Fig. 1 ), the x-y axis drift in the image stacks was stabilized using the Turboreg plugin in ImageJ. All ROIs were then detected using GECIquant. A scale was first assigned to image stacks, based on the confocal digital zoom setting. For most images, we used a 3x digital zoom, which corresponds to a scale of 0.23 μm per pixel. Briefly, a temporal projection of the movie stack was thresholded and the soma was detected with an area criterion of 30 μm2 to infinity within GECIquant. To do this, a temporal maximum intensity projection image was first generated by GECIquant from the image stack. The projection image was manually thresholded by the user with the default setting in ImageJ. Following thresholding, a polygon selection was manually drawn around the approximate astrocyte territory of interest, and the selection was added to the ImageJ ROI manager. Note that the assignment of territory was approximate and was not used for analysis except for the specific data set shown in Fig 2b . All ROIs falling within the range of 30 μm2 to infinity inside the polygon selection were detected by GECIquant and added to the ROI manager. An area range of 30 μm2 to infinity allowed detection of the astrocyte somata in all cases. The resulting detection was visually checked in every case.
To detect wave and microdomain ROIs, we first demarcated and deleted the soma from original image stacks using the clear selection feature in ImageJ. This was done because the increased basal fluorescence from the astrocyte soma relative to the processes prevented accurate thresholding of images for detection of ROIs within astrocyte processes. The ROI detection module in GECIquant was launched and the microdomain ROI option was selected. Microdomains and expanding wave ROIs were detected in separate analysis sessions. We used an area range of 0.5 to 4 μm2 to detect microdomains and an area range of 5 to 2000 μm2 for waves. These values were chosen after initial examination of the movie frames and by using several initial “best guess” test values as a guide. Other researchers who use GECIquant will also need to invest time initially to try several “best guess” values as a way to know what values will work best for the particular cell and fluctuation they are interested in measuring. The values we report here were appropriate for our experiments. For ROI detection, GECIquant generated a temporal maximum intensity projection image from the provided image stack with the deleted cell body. The projection image was manually thresholded by the user and a polygonal selection was manually drawn around the astrocyte of interest. GECIquant automatically detected microdomain and expanding wave ROIs based on the provided area criteria and the ROIs were added to the ImageJ ROI manager. Intensity values for each ROI were extracted in ImageJ and converted to dF/F values. For each ROI, basal F was determined during 50 s periods with no fluctuations. MiniAnalysis 6.0.07 (Synaptosoft) software was used to detect and measure amplitude, half width and frequency values for the somatic, wave and microdomain transients.
We comment on how we analyzed data for the experiments shown inFig 2 . First, for the analyses shown in Fig 2c , we made approximate ROIs that encompassed whole territories and then plotted the intensity of these regions over 300 s. From these traces, we measured the mean fluorescence intensity values over the 300 s period for each cell, and then averaged these values across all cells to generate the graphs in Fig 2c for WT and IP3R2−/− mice. In the case of the graphs shown in Fig 2d , we pooled the individual microdomain and wave Ca2+ fluctuations per cell, obtained the average value per cell of these pooled fluctuations and repeated this procedure for all cells. Then we averaged across all cells to generate the graphs that are shown in Fig 2c for WT and IP3R2−/− mice.
Graphs were made with Origin 8.1 and the figures assembled in CorelDraw 12 (Corel Corporation). No statistical methods were used to pre-determine sample sizes but our sample sizes are similar to those generally employed in the field. Randomization and blinding was not employed. Statistical comparisons were made using unpaired non parametric Mann-Whitney or unpaired parametric Student’s t tests as deemed appropriate after analyzing the raw data to ascertain whether they were normally-distributed using the Dallal and Wilkinson approximation to Lilliefors' method within Instat. When a statistical test was used, the precise P value and the test employed are reported in the text and/or figures legends. If the P was less than 0.00001, then it is reported as P < 0.00001. Otherwise, precise P values are provided in each case.
A methods checklist is available with thesupplementary materials .
Having analyzed all the data shown in this study, we observed three distinct types of spontaneous subcellular Ca2+ fluctuations within astrocytes, which we describe below and then clarify how they were detected within GECIquant. We classify Ca2+ fluctuations as: (1) somatic fluctuations that occur within the somata (these are restricted to somata and initial segments of processes arising from somata), (2) waves that occur exclusively within astrocyte processes and display repeated wave-like expansions and contractions of Ca2+, and (3) microdomain Ca2+ fluctuations that display highly restricted areas in astrocyte processes. These do not expand or contract as waves and remain restricted. The distinct areas covered by these three types of fluctuations are reported in the main text and in
After an image series was acquired (e.g.
To detect wave and microdomain ROIs, we first demarcated and deleted the soma from original image stacks using the clear selection feature in ImageJ. This was done because the increased basal fluorescence from the astrocyte soma relative to the processes prevented accurate thresholding of images for detection of ROIs within astrocyte processes. The ROI detection module in GECIquant was launched and the microdomain ROI option was selected. Microdomains and expanding wave ROIs were detected in separate analysis sessions. We used an area range of 0.5 to 4 μm2 to detect microdomains and an area range of 5 to 2000 μm2 for waves. These values were chosen after initial examination of the movie frames and by using several initial “best guess” test values as a guide. Other researchers who use GECIquant will also need to invest time initially to try several “best guess” values as a way to know what values will work best for the particular cell and fluctuation they are interested in measuring. The values we report here were appropriate for our experiments. For ROI detection, GECIquant generated a temporal maximum intensity projection image from the provided image stack with the deleted cell body. The projection image was manually thresholded by the user and a polygonal selection was manually drawn around the astrocyte of interest. GECIquant automatically detected microdomain and expanding wave ROIs based on the provided area criteria and the ROIs were added to the ImageJ ROI manager. Intensity values for each ROI were extracted in ImageJ and converted to dF/F values. For each ROI, basal F was determined during 50 s periods with no fluctuations. MiniAnalysis 6.0.07 (Synaptosoft) software was used to detect and measure amplitude, half width and frequency values for the somatic, wave and microdomain transients.
We comment on how we analyzed data for the experiments shown in
Graphs were made with Origin 8.1 and the figures assembled in CorelDraw 12 (Corel Corporation). No statistical methods were used to pre-determine sample sizes but our sample sizes are similar to those generally employed in the field. Randomization and blinding was not employed. Statistical comparisons were made using unpaired non parametric Mann-Whitney or unpaired parametric Student’s t tests as deemed appropriate after analyzing the raw data to ascertain whether they were normally-distributed using the Dallal and Wilkinson approximation to Lilliefors' method within Instat. When a statistical test was used, the precise P value and the test employed are reported in the text and/or figures legends. If the P was less than 0.00001, then it is reported as P < 0.00001. Otherwise, precise P values are provided in each case.
A methods checklist is available with the