Astrocytes are a major cell type in the mammalian brain. They are not electrically excitable, but generate prominent Ca2+ signals related to a wide variety of critical functions. The mechanisms driving these Ca2+ events remain incompletely understood. In this study, we integrate Ca2+ imaging, quantitative data analysis, and mechanistic computational modeling to study the spatial and temporal heterogeneity of cortical astrocyte Ca2+ transients evoked by focal application of ATP in mouse brain slices. Based on experimental results, we tune a single-compartment mathematical model of IP3-dependent Ca2+ responses in astrocytes and use that model to study response heterogeneity. Using information from the experimental data and the underlying bifurcation structure of our mathematical model, we categorize all astrocyte Ca2+ responses into four general types based on their temporal characteristics: Single-Peak, Multi-Peak, Plateau, and Long-Lasting responses. We find that the distribution of experimentally-recorded response types depends on the location within an astrocyte, with somatic responses dominated by Single-Peak (SP) responses and large and small processes generating more Multi-Peak responses. On the other hand, response kinetics differ more between cells and trials than with location within a given astrocyte. We use the computational model to elucidate possible sources of Ca2+ response variability: (1) temporal dynamics of IP3, and (2) relative flux rates through Ca2+ channels and pumps. Our model also predicts the effects of blocking Ca2+ channels/pumps; for example, blocking store-operated Ca2+ (SOC) channels in the model eliminates Plateau and Long-Lasting responses (consistent with previous experimental observations). Finally, we propose that observed differences in response type distributions between astrocyte somas and processes can be attributed to systematic differences in IP3 rise durations and Ca2+ flux rates.

Complex interactions between networks of astrocytes and neurons are beginning to be appreciated, but remain poorly understood. Transgenic mice expressing fluorescent protein reporters of cellular activity, such as the GCaMP family of genetically encoded calcium indicators (GECIs), have been used to explore network behavior. However, in some cases, it may be desirable to use long-established rat models that closely mimic particular aspects of human conditions such as Parkinson's disease and the development of epilepsy following status epilepticus. Methods for expressing reporter proteins in the rat brain are relatively limited. Transgenic rat technologies exist but are fairly immature. Viral-mediated expression is robust but unstable, requires invasive injections, and only works well for fairly small genes (<5 kb). In utero electroporation (IUE) offers a valuable alternative. IUE is a proven method for transfecting populations of astrocytes and neurons in the rat brain without the strict limitations on transgene size. We built a toolset of IUE plasmids carrying GCaMP variants 3, 6s, or 6f driven by CAG and targeted to the cytosol or the plasma membrane. Because low baseline fluorescence of GCaMP can hinder identification of transfected cells, we included the option of co-expressing a cytosolic tdTomato protein. A binary system consisting of a plasmid carrying a piggyBac inverted terminal repeat (ITR)-flanked CAG-GCaMP-IRES-tdTomato cassette and a separate plasmid encoding for expression of piggyBac transposase was employed to stably express GCaMP and tdTomato. The plasmids were co-electroporated on embryonic days 13.5–14.5 and astrocytic and neuronal activity was subsequently imaged in acute or cultured brain slices prepared from the cortex or hippocampus. Large spontaneous transients were detected in slices obtained from rats of varying ages up to 127 days. In this report, we demonstrate the utility of this toolset for interrogating astrocytic and neuronal activity in the rat brain.

We study evoked calcium dynamics in astrocytes, a major cell type in the mammalian brain. Experimental evidence has shown that such dynamics are highly variable between different trials, cells, and cell subcompartments. Here we present a qualitative analysis of a recent mathematical model of astrocyte calcium responses. We show how the major response types are generated in the model as a result of the underlying bifurcation structure. By varying key channel parameters, mimicking blockers used by experimentalists, we manipulate this underlying bifurcation structure and predict how the distributions of responses can change. We find that store-operated calcium channels, plasma membrane bound channels with little activity during calcium transients, have a surprisingly strong effect, underscoring the importance of considering these channels in both experiments and mathematical settings. Variation in the maximum flow in different calcium channels is also shown to determine the range of stable oscillations, as well as set the range of frequencies of the oscillations. Further, by conducting a randomized search through the parameter space and recording the resulting calcium responses, we create a database that can be used by experimentalists to help estimate the underlying channel distribution of their cells.